Unformatted text preview: Lecture 13 - Heat Transfer Applied Computational Fluid Dynamics
Instructor: André Bakker http://www.bakker.org © André Bakker (2002-2006) © Fluent Inc. (2002) 1 Introduction
• Typical design problems involve the determination of:
– Overall heat transfer coefficient, e.g. for a car radiator. – Highest (or lowest) temperature in a system, e.g. in a gas turbine, chemical reaction vessels, food ovens. – Temperature distribution (related to thermal stress), e.g. in the walls of a spacecraft. – Temperature response in time dependent heating/cooling problems, e.g. engine cooling, or how fast does a car heat up in the sun and how is it affected by the shape of the windshield? 2 Modes of heat transfer
• Conduction: diffusion of heat due to temperature gradients. A measure of the amount of conduction for a given gradient is the heat conductivity. • Convection: when heat is carried away by moving fluid. The flow can either be caused by external influences, forced convection; or by buoyancy forces, natural convection. Convective heat transfer is tightly coupled to the fluid flow solution. • Radiation: transfer of energy by electromagnetic waves between surfaces with different temperatures, separated by a medium that is at least partially transparent to the (infrared) radiation. Radiation is especially important at high temperatures, e.g. during combustion processes, but can also have a measurable effect at room temperatures. 3 Overview dimensionless numbers
• Nusselt number: Nu = hL / k f . Ratio between total heat transfer in a convection dominated system and the estimated conductive heat transfer. • Grashof number: Gr = L3 g∆ρ /ν 2 ρ w . Ratio between buoyancy forces and viscous forces. • Prandtl number: Pr = µ c p / k . Ratio between momentum diffusivity and thermal diffusivity. Typical values are Pr = 0.01 for liquid metals; Pr = 0.7 for most gases; Pr = 6 for water at room temperature. • Rayleigh number: Ra = Gr Pr = L3 ρ 2 gβ c p ∆T / µ k = L3 ρgβ∆T / µ α The Rayleigh number governs natural convection phenomena. • Reynolds number: Re = ρUL / µ . Ratio between inertial and viscous forces.
4 Enthalpy equation
• In CFD it is common to solve the enthalpy equation, subject to a wide range of thermal boundary conditions.
– Energy sources due to chemical reaction are included for reacting flows. – Energy sources due to species diffusion are included for multiple species flows. – The energy source due to viscous heating describes thermal energy created by viscous shear in the flow.This is important when the shear stress in the fluid is large (e.g. lubrication) and/or in high-velocity, compressible flows. Often, however, it is negligible. – In solid regions, a simple conduction equation is usually solved, although convective terms can also be included for moving solids. 5 Conjugate heat transfer
• “Conjugate heat transfer” refers to the ability to compute conduction of heat through solids, coupled with convective heat transfer in a fluid. • Coupled boundary conditions are available for wall zones that separate two cell zones. • Either the solid zone or the fluid zone, or both, may contain heat sources. • The example here shows the temperature profile for coolant flowing over fuel rods that generate heat. Grid Velocity vectors Temperature contours
Example: Cooling flow over fuel rods 6 Periodic heat transfer
• Also known as streamwise-periodic or fully-developed flow. • Used when flow and heat transfer patterns are repeated, e.g:
– Compact heat exchangers. – Flow across tube banks. • Geometry and boundary conditions repeat in streamwise direction. inflow outflow Outflow at one periodic boundary is inflow at the other
7 Heat conduction - Fourier’s law
• The heat flux is proportional to the temperature gradient: Q = q = − k ∇T A
where k(x,y,z,T) is the thermal conductivity. • In most practical situations conduction, convection, and radiation appear in combination. Also for convection, the heat transfer coefficient is important, because a flow can only carry heat away from a wall when that wall is conducting. temperature profile
dT dx hot wall cold wall x 8 Generalized heat diffusion equation
• If we perform a heat balance on a small volume of material…. heat conduction in T & q heat conduction out heat generation • … we get: rate of change of temperature ∂T & ρc = k∇2T + q ∂t heat cond. heat in/out generation α= k = thermal diffusivity ρc 9 Conduction example
• Compute the heat transfer through the wall of a home: Tout = 20° F Tout = 68° F Although slight, you can see the “thermal bridging” effect through the studs 2x6 stud k=0.15 W/m2-K sheetrock k=0.4 W/m2-K shingles k=0.15 W/m2-K fiberglass sheathing insulation 2-K k=0.15 W/m k=0.004 W/m2-K
10 Convection heat transfer
• Convection is movement of heat with a fluid. • E.g., when cold air sweeps past a warm body, it draws away warm air near the body and replaces it with cold air. flow over a heated block 11 Forced convection example
• Developing flow in a pipe (constant wall temperature). Tw
T Tw T∞ Tw T∞ Tw Tw bulk fluid temperature heat flux from wall x
12 Thermal boundary layer
• Just as there is a viscous boundary layer in the velocity distribution, there is also a thermal boundary layer. • Thermal boundary layer thickness is different from the thickness of the (momentum) viscous sublayer, and fluid dependent. The thickness of the thermal sublayer for a high-Prandtl-number fluid (e.g. water) is much less than the momentum sublayer thickness. For fluids of low Prandtl numbers (e.g. liquid metal), it is much larger than the momentum sublayer thickness. thermal boundary velocity boundary layer edge layer edge T∞ y T∞ , U
T ( y) δt δ δt ∝ L Tw Tw 13 Natural convection
• Natural convection (from a heated vertical plate). • As the fluid is warmed by the plate, its density decreases and a buoyant force arises which induces flow in the vertical direction. The force is proportional to ( ρ − ρ ∞ ) g . • The dimensionless group that governs natural convection is the Rayleigh number: T Tw
ρ = ρ (T ) u gravity Ra = Gr.Pr = gβ∆TL3 T∞ , ρ ∞ αυ
1 4 • Typically: Nu ∝ Ra x <x<1 3
14 Natural convection around a person
• Light weight warm air tends to move upward when surrounded by cooler air. • Thus, warm-blooded animals are surrounded by thermal plumes of rising warm air. • This plume is made visible by means of a Schlieren optical system that is based on the fact that the refraction of light through a gas is dependent on the density of the gas. • Although the velocity of the rising air is relatively small, the Reynolds number for this flow is on the order of 3000.
15 Natural convection - Boussinesq model
• Makes simplifying assumption that density is uniform.
– Except for the body force term in the momentum equation, which is replaced by: ( ρ − ρ 0 ) g = − ρ 0 β ( T − T0 ) g – Valid when density variations are small (i.e. small variations in T). • Provides faster convergence for many natural-convection flows than by using fluid density as function of temperature because the constant density assumptions reduces non-linearity. • Natural convection problems inside closed domains:
– For steady-state solver, Boussinesq model must be used. Constant density ρo allows mass in volume to be defined. – For unsteady solver, Boussinesq model or ideal gas law can be used. Initial conditions define mass in volume. 16 Newton’s law of cooling
• Newton described the cooling of objects with an arbitrary shape in a pragmatic way. He postulated that the heat transfer Q is proportional to the surface area A of the object and a temperature difference ∆T. • The proportionality constant is the heat transfer coefficient h(W/m2-K). This empirical constant lumps together all the information about the heat transfer process that we don’t know or don’t understand. T∞
Tbody q Q = q A = h A (Tbody − T∞ ) = h A ∆T
average heat transfer coefficient (W/m2-K)
17 Heat transfer coefficient
• h is not a constant, but h = h(∆T). • Three types of convection. • Natural convection. Fluid moves due to buoyancy. Typical values of h: Thot Tcold 4 - 4,000 W/m2-K h ∝ ∆T x 1 4 <x< 1 3 • Forced convection: flow is induced by external means. Tcold Thot 80 - 75,000 h = const
• Boiling convection: body is hot enough to boil liquid. Tcold Thot 300 - 900,000 h ∝ ∆T 2
18 Radiation heat transfer
• Thermal radiation is emission of energy as electromagnetic waves. • Intensity depends on body temperature and surface characteristics. • Important mode of heat transfer at high temperatures, e.g. combustion. • Can also be important in natural convection problems. • Radiation properties can be strong functions of chemical composition, especially CO2, H2O. • Radiation heat exchange is difficult solve (except for simple configurations). We must rely on computational methods. 19 Surface characteristics
(incident energy flux) ℜq (reflected) α q (absorbed) translucent slab τ q (transmitted) 1=α + ℜ +τ
20 transmittance Black body radiation
• A “black body”:
– Is a model of a perfect radiator. – Absorbs all energy that reaches it; reflects nothing. – Therefore α = 1, ℜ = τ = 0. • The energy emitted by a black body is the theoretical maximum: q =σ T4
• This is Stefan-Boltzmann law; σ is the Stefan-Boltzmann constant (5.6697E-8 W/m2K4). • The wavelength at which the maximum amount of radiation occurs is given by Wien’s law: λmax T = 2.898E − 3 [mK] • Typical wavelengths are λmax = 10 µm (far infrared) at room temperature and λmax = 0.5 µm (green) at 6000K.
21 Real bodies
• Real bodies will emit less radiation than a black body: q = εσ T • Here ε is the emissivity, which is a number between 0 and 1. Such a body would be called “gray” because the emissivity is the average over the spectrum. • Example: radiation from a small body to its surroundings.
– Both the body and its surroundings emit thermal radiation. – The net heat transfer will be from the hotter to the colder.
4 4 • The net heat transfer is then: Qnet = ε Aσ (Tw − T∞ ) • For small ∆T the term (Tw4-T∞4) can be 3 approximated as 4T (Tw − T∞ ) and Qnet = A hr ∆T with hr as an effective radiation heat transfer coefficient. 4 Qnet T∞ q∞ qw Tw A
22 The radiative heat transfer equation (RTE) 23 Radiation
• Radiation intensity transport equations (RTE) are solved.
– Local absorption by fluid and at boundaries links energy equation with RTE. • Radiation intensity is directionally and spatially dependent.
– Intensity along any direction can be reduced by:
• Local absorption. • Out-scattering (scattering away from the direction). – Intensity along any direction can be augmented by:
• Local emission. • In-scattering (scattering into the direction). • Four common radiation models are:
– – – – Discrete Ordinates Model (DOM). Discrete Transfer Radiation Model (DTRM). P-1 Radiation Model. Rosseland Model. 24 Discrete ordinates model
• The radiative transfer equation is solved for a discrete number of finite solid angles: ∂I s i σT 4 σ s + (a + σ s )I ( r , s ) = an 2 + ∂ xi π 4π
absorption 4π ∫ I (r , s ' )Φ ( s ⋅ s ' )dΩ '
0 • Advantages:
– – – – – emission scattering Solution method similar to that for the other conservation equations. Conservative method leads to heat balance for coarse discretization. Accuracy can be increased by using a finer discretization. Accounts for scattering, semi-transparent media, specular surfaces. Banded-gray option for wavelength-dependent transmission. • Limitations:
– Solving a problem with a large number of ordinates is CPU-intensive.
25 Discrete transfer radiation model (DTRM)
• Main assumption: radiation leaving surface element in a specific range of solid angles can be approximated by a single ray. • Uses ray-tracing technique to integrate radiant intensity along each ray: dI σT 4 = −α I + α ds π • Advantages:
– Relatively simple model. – Can increase accuracy by increasing number of rays. – Applies to wide range of optical thicknesses. • Limitations:
– Assumes all surfaces are diffuse (isotropic reflection). – Effect of scattering not included. – Solving a problem with a large number of rays is CPU-intensive.
26 P-1 model
• Main assumption: radiation intensity can be decomposed into series of spherical harmonics. Only first term in this (rapidly converging) series used in P-1 model. • Advantages:
– Radiative transfer equation easy to solve with little CPU demand. – Works reasonably well for combustion applications where optical thickness is large. – Easily applied to complicated geometries with curvilinear coordinates. – Effects of particles, droplets, and soot can be included. • Limitations:
– Assumes all surfaces are diffuse. – May result in loss of accuracy, depending on complexity of geometry, if optical thickness is small. – Tends to overpredict radiative fluxes from localized heat sources or sinks. • A further simplified (and faster but less accurate) version of this model is the Rosseland model.
27 Choosing a radiation model
• For certain problems, one radiation model may be more appropriate in general.
– Computational effort: P-1 gives reasonable accuracy with less effort. – Accuracy: DTRM and DOM more accurate. – Optical thickness: DTRM/DOM for optically thin media (optical thickness << 1); P-1 sufficiently accurate and faster for optically thick media. – Scattering: P-1 and DOM account for scattering. – Particulate effects: P-1 and DOM account for radiation exchange between gas and particulates. – Localized heat sources: DTRM/DOM with sufficiently large number of rays/ordinates is more appropriate. • DOM is the only model that works across the full range of conditions.
28 Wall heat flux calculation
heat flux : q = h f (Tw − T f ) + qrad ∆y for turbulent flows : h f follows from correlations describing the thermal boundary layer profile external wall radiation : qrad = ε extσ (T∞4 − Tw4 ) Tw = wall temperatur e T f = fluid cell temperatur e T∞ = user specified temperatur e for laminar flows : h f = kf ε ext = emissivity of external wall surface σ = Stefan − Boltzmann constant
h f = fluid side local heat transfer coefficient ∆y = normal distance fluid cell center to wall
29 Heat transfer optimization
• We have the following relations for heat transfer: – Conduction: Q = A ∆T k f / d – Convection: Q = A h ∆T – Radiation: Q = A hr ∆T • As a result, when equipment designers want to improve heat transfer rates, they focus on:
– Increasing the area A, e.g. by using profiled pipes and ribbed surfaces. – Increasing ∆T (which is not always controllable). – For conduction, increasing kf /d. – Increase h by not relying on natural convection, but introducing forced convection. – Increase hr, by using “black” surfaces.
30 Fluid properties
• Fluid properties such as heat capacity, conductivity, and viscosity can be defined as:
– – – – – Constant. Temperature-dependent. Composition-dependent. Computed by kinetic theory. Computed by user-defined functions. • Density can be computed by ideal gas law. • Alternately, density can be treated as:
– – – – Constant (with optional Boussinesq modeling). Temperature-dependent. Composition-dependent. User defined functions.
31 Phase change
heat capacity latent heat • Systems in which phase change occurs (e.g. melting, solidification, and sometimes evaporation) can be modeled as a single-phase flow with modified physical properties. • In that case, the medium gets the properties of one phase state below a certain critical temperature, and the properties of the other phase state above a second critical temperature. • Linear transitions for µ and ρ. • A “spike” in cp is added, the area of which corresponds to the latent heat. density conductivity viscosity Temperature • A second spike is added to the heat conductivity curve, to keep the ratio between heat capacity and thermal conductivity constant.
32 Thermal boundary conditions
• At flow inlets and exits.
– At flow inlets, must supply fluid temperature. – At flow exits, fluid temperature extrapolated from upstream value. – At pressure outlets, where flow reversal may occur, “backflow” temperature is required. • Thermal conditions for fluids and solids.
– Can specify energy source. • Thermal boundary conditions at walls.
– – – – – Specified heat flux. Specified temperature. Convective heat transfer. External radiation. Combined external radiation and external convective heat transfer.
33 Notes on convergence
• Heat transfer calculations often converge slowly. It is recommended to use underrelaxation factors of 0.9 or larger for enthalpy. If lower underrelaxation factors are used, obtaining a good solution may take prohibitively long. • If underrelaxation factors of 0.2 or lower have to be used to prevent divergence, it usually means that the model is ill-posed. • Deep convergence is usually required with scaled residuals having to be of the order 1E-6 or smaller. 34 Example: heat exchanger efficiency
• Problem: improve the efficiency of a tube-cooled reactor. • Non-standard design, i.e. traditional correlation based methods not applicable. • Solution: more uniform flow distribution through the shell that will result in a higher overall heat transfer coefficient and improved efficiency. 35 Heat exchanger - original design
Baffle “A” • Original design:
Injectors – – – – Bundle of tubes as shown. Repeated geometry. 3 different baffles, A, C, and D. Reactant injectors between baffles “A” and “D”. Fl ow on cti re Di Baffle “C” Baffle “D” 36 Heat exchanger - modeling approach
• 3-dimensional, steady, turbulent, incompressible, isothermal. • Bundle of tubes modeled as a non-isotropic porous medium. Two symmetry planes significantly reduce domain size. • Hybrid, unstructured mesh of 330,000 cells. • Zero thickness walls for baffles. • Leakage between baffles and shell wall (0.15” gap) modeled using thin prism cells. • Uniform inflow applied over a half-cylindrical surface upstream of the first baffle.
37 Heat exchanger - flow pattern
D A C A D A C A D Bundle of tubes Flow Direction Injection points (Shaded area) Compartments with low flow Compartment mainly served by nearest upstream injector Non-uniform flow distribution means low efficiency 38 Heat exchanger - modifications
• Design modifications:
Baffle “A” Relocated Injectors – Shorter baffles “ C’ ”. – Relocated and rotated injectors. Modified (Shorter) Baffle “ C’ ”
n tio ec r Di Fl ow Baffle “D” 39 Heat exchanger - improved flow pattern
• • • • Flow distribution after modifications. No recirculating fluid between baffles C’ and A. Almost uniform flow distribution. Problem has unique flow arrangement that does not allow traditional methods to be of any help. • A simplified CFD model leads to significantly improved performance of the heat exchanger/reactor.
D A C’ A D A C’ A D 40 Conclusion
• Heat transfer is the study of thermal energy (heat) flows: conduction, convection, and radiation. • The fluid flow and heat transfer problems can be tightly coupled through the convection term in the energy equation and when physical properties are temperature dependent. • Chemical reactions, such as combustion, can lead to source terms to be included in the enthalpy equation. • While analytical solutions exist for some simple problems, we must rely on computational methods to solve most industrially relevant applications. 41 ...
View Full Document
- Fall '10
- Heat Transfer