Study on Effect of Different Porosity on Thermally Driven Heat Transfer in a Centrifugal Force Field

Similar documents
The Analysis of Convection Experiment

Large Scale Electrostatic Precipitator

Principles of Food and Bioprocess Engineering (FS 231) Solutions to Example Problems on Heat Transfer

Energy configuration optimization of submerged propeller in oxidation ditch based on CFD

Lattice Boltzmann simulation of nucleate boiling in micro-pillar structured surface

PERFORMANCE OF HEAVY-DUTY PLANETARY GEARS

Numerical simulation of water flow in an axial flow pump with adjustable guide vanes

Simulation and experiment of the effect of clearance of impeller wear-rings on the performance of centrifugal pump

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

High resolution entropy stable scheme for shallow water equations

829. An adaptive method for inertia force identification in cantilever under moving mass

Air Age Equation Parameterized by Ventilation Grouped Time WU Wen-zhong

Week 9 Chapter 10 Section 1-5

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

Turbulence classification of load data by the frequency and severity of wind gusts. Oscar Moñux, DEWI GmbH Kevin Bleibler, DEWI GmbH

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate

A Comparative Investigation into Aerodynamic Performances of Two Set Finned Bodies with Circular and Non Circular Cross Sections

Supplementary Information for Observation of Parity-Time Symmetry in. Optically Induced Atomic Lattices

FORCED CONVECTION HEAT TRANSFER FROM A RECTANGULAR CYLINDER: EFFECT OF ASPECT RATIO

modeling of equilibrium and dynamic multi-component adsorption in a two-layered fixed bed for purification of hydrogen from methane reforming products

STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS

Temperature. Chapter Heat Engine

5th International Conference on Measurement, Instrumentation and Automation (ICMIA 2016) Star image identification uninfluenced by rotation

MEASUREMENT OF MOMENT OF INERTIA

1. Int In e t rnal Losses: tak ta e k place plac in the inner passages adding heat to the flow medium 2. External losses:

Assessment of Site Amplification Effect from Input Energy Spectra of Strong Ground Motion

Electrical double layer: revisit based on boundary conditions

Numerical Transient Heat Conduction Experiment

Fluid structure interaction analysis of a high-pressure regulating valve of a 600-MW ultra-supercritical steam turbine

The unsteady flow characteristic research on the initial period flow of micro channel

Study on Active Micro-vibration Isolation System with Linear Motor Actuator. Gong-yu PAN, Wen-yan GU and Dong LI

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

Outline. Unit Eight Calculations with Entropy. The Second Law. Second Law Notes. Uses of Entropy. Entropy is a Property.

2010 Black Engineering Building, Department of Mechanical Engineering. Iowa State University, Ames, IA, 50011

STUDY OF A THREE-AXIS PIEZORESISTIVE ACCELEROMETER WITH UNIFORM AXIAL SENSITIVITIES

Simulation of Gas Exothermic Chemical Reaction in Porous Media Reactor with Lattice Boltzmann Method

The Order Relation and Trace Inequalities for. Hermitian Operators

NON LINEAR ANALYSIS OF STRUCTURES ACCORDING TO NEW EUROPEAN DESIGN CODE

Numerical analysis on the agglomeration behavior of fine particles in plane jet

Operating conditions of a mine fan under conditions of variable resistance

ERROR RESEARCH ON A HEPA FILTER MEDIA TESTING SYSTEM OF MPPS(MOST PENETRATION PARTICLE SIZE) EFFICIENCY

REAL-TIME DETERMINATION OF INDOOR CONTAMINANT SOURCE LOCATION AND STRENGTH, PART II: WITH TWO SENSORS. Beijing , China,

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

An Improved multiple fractal algorithm

Orientation Model of Elite Education and Mass Education

A Numerical Study of Heat Transfer and Fluid Flow past Single Tube

International Power, Electronics and Materials Engineering Conference (IPEMEC 2015)

Simulation of Flow Pattern in Open Channels with Sudden Expansions

The Study of Teaching-learning-based Optimization Algorithm

Experiment 1 Mass, volume and density

Airflow and Contaminant Simulation with CONTAM

Relationship between Refractive Index and Molar Concentration of Multi-Component Solutions Zhu Xingyu 1, a, Mai Tiancheng 2, b and Zhao Zilong 2, c

EFFECT OF THE SIZE OF GRADED BAFFLES ON THE PERFORMANCE OF CHANNEL HEAT EXCHANGERS

First Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force.

Investigation of a New Monte Carlo Method for the Transitional Gas Flow

Inductance Calculation for Conductors of Arbitrary Shape

Homogeneous model: Horizontal pipe and horizontal well. Flow loops can't duplicate field conditions. Daniel D. Joseph. April 2001

NUMERICAL MODEL FOR NON-DARCY FLOW THROUGH COARSE POROUS MEDIA USING THE MOVING PARTICLE SIMULATION METHOD

Back Analysis of Self-Compacting Concrete Rheological Parameters Based on H-B Model

Constitutive Modelling of Superplastic AA-5083

Analysis of the Magnetomotive Force of a Three-Phase Winding with Concentrated Coils and Different Symmetry Features

INTERROGATING THE FLOW BEHAVIOUR IN A NOVEL MAGNETIC DESICCANT VENTILATION SYSTEM USING COMPUTATIONAL FLUID DYNAMICS (CFD)

Turbulent Heat Transfer Modelling in a Vessel Stirred by a Pitched Blade Turbine Impeller

NUMERICAL SIMULATION OF FLOW OVER STEPPED SPILLWAYS

Modal Strain Energy Decomposition Method for Damage Detection of an Offshore Structure Using Modal Testing Information

Physics 2102 Spring 2007 Lecture 10 Current and Resistance

A large scale tsunami run-up simulation and numerical evaluation of fluid force during tsunami by using a particle method

Natural Convection in a Rectangular Enclosure with Colliding Boundary Layers

Thermal State Investigations of Two-stage Thruster with Anode Layer

Module 3: The Whole-Process Perspective for Thermochemical Hydrogen

A New Scrambling Evaluation Scheme based on Spatial Distribution Entropy and Centroid Difference of Bit-plane

Fast Simulation of Pyroshock Responses of a Conical Structure Using Rotation-Superposition Method

THE NEAR-WALL INFLUENCE ON THE FLOW AROUND A SINGLE SQUARE CYLINDER.

A PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.

CHEMICAL ENGINEERING

Optimization of the thermodynamic model of a solar driven Aqua - ammonia absorption refrigeration system

A Fast Computer Aided Design Method for Filters

HEAT TRANSFER IN ELECTRONIC PACKAGE SUBMITTED TO TIME-DEPENDANT CLIMATIC CONDITIONS

Assignment 5. Simulation for Logistics. Monti, N.E. Yunita, T.

2 Finite difference basics

Numerical Heat and Mass Transfer

THERMAL PERFORMANCE ENHANCEMENT OF SOLAR PARABOLIC TROUGH COLLECTOR USING SECONDARY REFLECTOR

Numerical Simulation and Experimental Study on a New Type of Variable-rate Fluidic Sprinkler

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity

Numerical Analysis of Heat Transfer and Pressure Drop in a Channel Equipped with Triangular Bodies in Side-By-Side Arrangement

Calculating the Quasi-static Pressures of Confined Explosions Considering Chemical Reactions under the Constant Entropy Assumption

PHYSICS - CLUTCH CH 28: INDUCTION AND INDUCTANCE.

The Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites

Universiti Malaysia Perlis, Pauh Putra Campus, Arau, Perlis, Malaysia. Jalan Permatang Pauh, Pauh, Pulau Pinang, Malaysia

Structure and Drive Paul A. Jensen Copyright July 20, 2003

A numerical study on the shell-side turbulent heat transfer enhancement of shell-and-tube heat exchanger with trefoilhole

Research & Reviews: Journal of Engineering and Technology

Lecture #06 Hotwire anemometry: Fundamentals and instrumentation

B and H sensors for 3-D magnetic property testing

Three-Phase Distillation in Packed Towers: Short-Cut Modelling and Parameter Tuning

Computer Simulations of Parallel-to-Series Conversion in Solid State Frame Transfer Image Sensors. J. Bisschop

Color Rendering Uncertainty

Adaptive sliding mode reliable excitation control design for power systems

Spatial Statistics and Analysis Methods (for GEOG 104 class).

Transcription:

5th Internatonal Conference on Cvl Engneerng and Transportaton (ICCET 015) Study on Effect of Dfferent Porosty on Thermally Drven Heat Transfer n a Centrfugal Force Feld 1 Xa Je1,a*, Chang Ha-png,b, and Tang Hou-xang3,c Department of Theory Tranng, Flght Drllmaster Tranng Base, Chna College of Energy and Power Engneerng, Nanng Unversty of Aeronautcs and Astronautcs, Chna 3 Department of Basc Course, Automoble Management Insttute, Chna a xx_chna@163.com, bchppe@nuaa.edu.cn, ctanghouxang@163.com Key words: porous medum; porosty; thermally drven Abstract: Based on the exploture of turbne blade super coolng technology, porous medum s nstalled n a new knd of coolng confguraton wth coolng tunnels. Experments and numercal smulatons are carred out to nvestgate the thermally drven heat transfer rules n the new knd of coolng confguraton wth dfferent porosty n a centrfugal force feld. The results of experments and numercal smulatons are consstent bascally. The results of study show that the heat transfer rules of the new coolng confguraton are dentcal under dfferent porosty. The thermally drven heat transfer of the coolng confguraton flled wth porous medum can be enhanced wth ncrease of the rotate speed, heat flux and coolng ar speed. At the same tme, the heat transfer effect can be weakened wth ncrease of the porosty under hgh porosty. Introducton That f ncreasng gas s nlet temperature effectvely and relably s the root of enhancng thermal effcency and thrust-weght rato accordng to gas turbne s thermal crculaton. Turbne blade s coolng desgn wll be the key to solve queston to reach future hgh-powered engne s technology ndex. Focusng on the stuaton, Professor Chang Ha-png and Academcan Guo Zeng-yuan brought up a new coolng technque based on the theory of thermal drve [1]. Slght coolng channels are the core of the technque. Coolng medum s put nto the slght coolng channels. Porous medum s nstalled n channels to enhance slght channels structure ntensty. At the same tme, flowng of the flud n the porous medum exst dsperson effect. It can make molecular groups mxed radally, make radal temperature dstrbuton unform and even, and eventually enhance heat exchange evdently [-4]. Along wth enhancng heat exchange, resstance to flow also augments a lot. But the ncreasng resstance can be overcome wth enough buoyancy. Based on the fundamental study[5-7], n ths paper, experments and numercal smulatons were carred out to nvestgate the thermally drven heat transfer rules n the new knd of coolng confguraton wth dfferent porosty and compare ther heat transfer effect. Expermental Research Expermental Faclty Fg1 s a scheme of expermental apparatus. The whole expermental apparatus contans power system, heatng system, coolng system and measurng system, etc. Coolng ar enters hollow shaft, passes gas collecton case nto expermental model, and streams out of the turntable. There are two layers hole boards n the gas collecton case to make gas enter the expermental model equably to smulate the turbne blade s work state. Fg and Fg3 are photos and sketch map of the testpeces. In order to obtan the closed channels llustrated n Fg3, one body shell, fve nner hulls and two end plates were manufactured. In the closed channels, foamed metals made wth colomony are used n the porous medum area(porosty s 98% and 93%), 015. The authors - Publshed by Atlants Press 683

whch s llustrated concretely n Fg. The external dmenson of the testpeces s 50mm*36mm*38mm. Except that the wdth of the porous medum channel s 3mm, the wdth of all the others s 1mm. Contract grant sponsors: Natonal Natural Scence Foundaton of Chna(507607) and the Foundatonal Research of the Natonal Defence (J1400D001) In the experments, the heatng theca posted on the body shell heats the cycle channels, and adustng electrc current can create dfferent heat. At the same tme, coolng ar goes through the nner hulls to form coolng channels, and takes away the heat dscharged by cycle flud. Adustng the ar flux can create dfferent coolng ar speed. The operatonal range of the experments ncludes: heat flux (5390-48500w/m ), coolng ar speed (3.1-96.m/s) and rotate speed (300-1380r/mn). In order to research the thermo-drvng heat transfer effectveness of flud n the closed channels nstalled wth porous medum, the author defnes a parameter (KH) of capacty of thermo-drvng. It characterzes the capacty of the flud n the closedchannels, whch takes heat to the coolng end wth the crcular flow of thermo-drvng. The concrete defnton ust as follows: q KH T max T mn 1.belt pulley 3.turntable skn 5.fexble couplng Where, q s heat flux of the heatng surface, T max and T mn respectvely s the hghest and lowest temperature measured from the surface of the walls along wth the orentaton of the crcular flow of thermo-drvng. Expermental Results and Dscusson Fg4, Fg5 and Fg6 llustrate the dagrams that the Fg. Photos of expermental mode temperature on the heatng surface changes along wth heat flux, rotatng speed and ntake velocty of the coolng ar, here the porosty s 93% n porous area. As llustrated n Fg4, wth the same heat flux densty, the temperature of the heatng end along wth the postve drecton of the Y axs consstently falls, because after absorbng heat, the temperature of flud ncreases and the densty decreases, and the flud affected by centrfugal force moves along wth the negatve drecton of the Y axs. That demonstrates the flow of thermo-drvng takes place n the enclosure after beng flled wth porous medum (If there s no thermo-drvng but only Fg.3 Schematc of expermental model thermo-conducton, the temperature on the heatng surface wll be unformty, and along wth the Y axs wll have no change. The temperature of the heatng end along wth the postve drecton of the Y axs consstently ncreases under rotaton. That s contrary to the stuaton under no rotaton.). Furthermore, the temperature at all the measurng ponts ncrease n varyng degrees along wth ncrease of heat flux. The hgher the heat flux, the greater the temperature dfference between the hghest and lowest measurng ponts. That demonstrates heat flux ncreases and thermo-drvng capacty mproves and then thermo-drvng heat transfer effectveness contnuously rses. Fg5 and Fg6 llustrate the temperature at all the measurng ponts decrease and the temperature dfference between the hghest and lowest measurng ponts decreases. That demonstrates 684 ar nlet ar outlet ar outlet 10 7.turntable 9.expermental model Fg.1 Schematc of expermental apparatus

rotatng speed (centrfugal acceleraton) or ntake velocty of the coolng ar ncreases, thermo-drvng heat transfer effectveness also rses. 360 350 n=900rpm 5390W/m 1100W/m 1600W/m 33700W/m 48500W/m 330 38 36 q=1600w/m 300rpm 540rpm 780rpm 100rpm 1380rpm 340 34 Temperature(k) 330 30 Temperature(K) 3 30 318 316 314 300 8 10 1 14 16 18 0 4 6 8 Y Dstance(mm) 31 8 10 1 14 16 18 0 4 6 8 Y Dstance(mm) Fg.4 Curves of temperature varyng wth heat fluxes Fg.5 Curves of temperature varyng wth rotatng speeds Temperature(K) 36 q=1600w/m 34 3 30 318 316 314 31 308 306 n=660rpm 8 10 1 14 16 18 0 4 6 8 Y Dstance(mm) 3.1m/s 37.4m/s 48.1m/s 64.1m/s 96.m/s Fg.6 Curves of temperature varyng wth coolng ar speeds 300 00 100 000 1900 1800 1700 1600 q=1600w/m ε =0.93% ε =0.98% 400 300 00 100 000 1900 1800 1700 1600 1500 1400 n=900rpm ε =0.93% ε =0.98% 1300 0 10000 0000 30000 40000 50000 350 q=1600w/m 300 50 00 150 100 050 000 1950 1900 1850 1800 1750 n=660rpm q(w/m ) Fg.7 Curves of KH varyng wth heat fluxes under two knds of porosty ε =93% ε =98% 1500 00 400 600 800 1000 100 1400 n(rpm) 30 40 50 60 70 80 90 100 v(m/s) Fg.8 Curves of KH varyng wth rotatng speeds Fg.9 Curves of KH varyng wth coolng ar speeds under two knds of porosty under two knds of porosty The paper researches the thermo-drvng heat transfer n the enclosure wth porous medum whose porosty s 93%, meanwhle, t makes a comparson between the results and those porosty s 98%. The 685

research demonstrates the enclosure has the same temperature dstrbuton rule under two knds of porosty. So, the temperature dstrbuton chart when the porosty s 98% s not gven n ths paper. Fg7, Fg8 and Fg9 respectvely demonstrate the dagrams of the comparson of the parameters of thermo-drvng heat transfer capacty between two knds of porosty when heat flux, rotatng speed and coolng ar speed change. The fgures llustrate that the parameter of thermo-drvng heat transfer capacty ncreases and thermo-drvng capacty mproves along wth ncrease of heat flux, rotatng speed and coolng ar speed under two knds of porosty. Although the heat transfer rule s the same, the enclosure wth 93% porosty has much better thermo-drvng heat transfer effectveness than the other. The reason s as follows. In porous medum area, sold medum s heat conducton coeffcent s hgher than the flud n ths examnaton. So, sold medum s heat conducton and convecton heat transfer between sold medum and flud have very mportant acton n porous medum area. The contact area between sold medum and flud ncreases along wth decrease of porosty. Consequently, t strengthens convecton heat transfer n porous medum area and mproves the enclosure s thermo-drvng heat transfer capacty and effectveness. Numercal Research Mathematcal Model In order to fnd out the flud s flow and heat transfer character, flow and temperature feld n the enclosure wth porous medum and the flud s thermo-drvng heat transfer capacty under dfferent porosty, the paper carres out numercal smulaton. The paper sets up mathematcal model accordng to expermental model s sze n numercal smulaton. Basc equatons are as follows. no porous medum u 0 x u u x P x x g u E P T k u x x x wth porous medum u 0 x 1 u u u x E P x P x x x k eff T x 1 ek e uk R x eff 1 u g ek e uk R x u eff s gyral angle velocty; s porosty; k eff s effectve heat conducton coeffcent(k eff =k f +(1-)k s, k f s flud heat conducton coeffcent, k s s sold heat conducton coeffcent); E s flud nner enerry; s the nfltraton rate of porous medum. The boundary condtons accordng to expermental parameters are as follows: Heatng surface s supposed changeless heat flux; coolng surface carres through heat transfer wth lqud n the enclosure and coolng ar, so t s supposed couplng surface; the nlet and outlet surface of the coolng ar channel are supposed velocty nlet and pressure outlet; all the other surfaces are supposed as perfect adabatc walls. In the numercal smulaton we adopt Gambt to buld grd and Fluent s Coupled Implct Steady Solver. Because there s shearng stran n flud feld under rotaton, we adopt RNG k- model and Two-Layer Model near wall. 686

Numercal Results and Analyss Fg10 s temperature dstrbutng of the expermental model. The temperature on the heatng surface along wth the postve drecton of the Y axs falls, whch s consstent wth expermental result. In order to verfy numercal result, the paper compares numercal heatng surface temperature wth expermental result. Fg11 s the comparson of numercal wth expermental result under the same condtons(q=5390w/m, n=900rpm, ). Fg11 llustrates that numercal and expermental result s change trand s the same. Numercal result s fve degrees hgher than expermental result, because the absolute adabatc condton s not realzed n expermentaton. 315 computaton experment Temperature(k) 305 300 95 0.000 0.005 0.010 0.015 0.00 0.05 0.030 Y Dstance(mm) Fg.10 Temperature dstrbutng of the expermental model Temperature(K) 317 316 315 314 313 31 311 309 308 307 306 305 304 303 0.000 0.005 0.010 0.015 0.00 0.05 0.030 Y Dstance(mm) ε =0.96 ε =0.93 ε =0.9 ε =0.86 ε =0.83 ε =0.8 1600 1500 1400 1300 100 1100 1000 Fg.11The comparson of numercal wth expermental result 900 0.78 0.80 0.8 0.84 0.86 0.88 0.90 0.9 0.94 0.96 0.98 ε Fg.1 Curves of temperature varyng wth porosty Fg.13 Curves of KH varyng wth porosty The research [9] ndcates that porous medum can strengthen heat transfer, but the flow drag ncreases more when the porosty s less than 0.8. So the desgn of the new coolng heat transfer confguraton takes bg porosty porous medum nto account. It can strengthen heat transfer, but the flow drag does not ncrease much. Therefore, the paper adopts bg porosty porous medum n numercal smulaton( 80%). Fg1 and Fg13 llustrate the dagrams that the temperature on the heatng surface and the parameters of thermo-drvng heat transfer capacty change along wth porosty under the condtons(q=5390w/m, n=900rpm, ). Fg1 llustrates heatng surface temperature and the temperature dfference between the hghest and lowest measurng ponts decreases along wth the decrease of porosty. That demonstrates porosty decreases, thermo-drvng heat transfer capacty rses. Fg13 also proves that. Fg1 also llustrates the declne extent of heatng surface temperature and the temperature dfference between the hghest and lowest measurng ponts decreases along wth the decrease of porosty. Fg13 llustrates the ncrease extent of the 687

parameter of thermo-drvng heat transfer capacty decreases along wth the decrease of porosty. These analyses demonstrate thermo-drvng heat transfer capacty ncreases but ts ncrease extent decreases along wth the decrease of porosty. The reason s as follows: when the porosty decreases, the addton of contact area of sold medum and flud strengthens convecton heat transfer, but the ncrease of flow drag weakens convecton heat transfer. It leads to the decrease of ncrease extent of thermo-drvng heat transfer capacty fnally. Concluson (1) The flow of thermo-drvng takes place n the new coolng confguraton after beng flled wth porous medum. () The new coolng confguraton has the same heat transfer rule under dfferent porosty: when rotatng speed, heat flux and ntake velocty of the coolng ar ncreases, ts thermo-drvng heat transfer capacty also rses. (3) The new coolng confguraton s thermo-drvng heat transfer capacty ncreases but ts ncrease extent decreases along wth the decrease of porosty. References [1] Chang H, Guo Z. A new coolng technque used n turbne blade. Inventon patent No. ZL01108, 001-04. [] M.Varahasamy,R.M.Fand.Heat Transfer by Forced Convecton n Ppes Packed wth Meda Whose Matrces Are Composed of Spheres.Int. [J].Heat Mass Transfer,1996, 39(18): 3931~3847. [3] Hsu.C.T.,Cheng P.. Closure Scheme of the Macroscopc Energy Equaton for Convecton Heat Transfer n Porous Meda.Int. [J].Commu.Heat Transfer, 1988, 15:689~703 [4] Morrson H L et al.[j]. Appl. Phys., 1949,0:107~109 [5] Mao Junku, Chang Hapng, Fang L. ExpermentalInvestgaton of Flow and Heat Transfer About thethermal Drvng of R1 n Rotatng Small-Sze CrcularEnclosures[J].Journal of Aerospace Power,003,18(4): 498~501. [6] Mao Junku, Chang Hapng, Fang L. Investgaton of Thermally Drven n Centrfugal Feld wth Lqud Crystal [J]. Journal of Engneerng Thermophyscs,004,5(): 9~31. [7]Yang Mn, Chang Hapng, Mao Junku. ExpermentalInvestgaton of Thermally Drvng Heat TransferCharacterstcs of HO and R1 n a Centrfugal Force Feld [J].Journal of Engneerng Thermophyscs,004,5(): 30~304. [8] Fluent Users Gude,1998 [9] Renken K.J,Poulkakous D. Experment and Analyss of Forced Convecton Heat Transport n a Packed Bed of Sphere, Int [J].Heat Mass Transfer 1998,31: 1399~1408 688