Lecture 13 - Heat Transfer Applied Computational Fluid Dynamics Instructor: Andr Bakker Andr Bakker (2002-2004) 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 gT / 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 highvelocity, 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 Heat conduction - Fouriers 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 hot wall dT dx cold wall x

8 Generalized heat diffusion equation If we perform a heat balance on a small volume of material. heat conduction in T heat conduction out q heat generation we get: T c k2T q t rate of change of temperature 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 2x6 stud k=0.15 W/m2-K Although slight, you can see the thermal bridging effect through the studs sheetrock k=0.4 W/m2-K shingles k=0.15 W/m2-K fiberglass sheathing

insulation 2 k=0.15 W/m -K 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 T Tw T Tw T

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) Tw t t L 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: gTL3 Ra Gr.Pr

Typically: Nu Ra x 1 4 T (T ) Tw u gravity T , x 13 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 Newtons 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 dont know or dont understand. T q Tbody Q q A h A (Tbody T ) h A T h 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. 1 1 x h T 4 x3

Forced convection: flow is induced by external means. Typical values of h: Thot Tcold Tcold Thot 4 - 4,000 W/m2-K 80 - 75,000 h const Boiling convection: body is hot enough to boil liquid. 2 Tcold Thot 300 - 900,000 h T 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 q (W/m2) q (reflected) (incident energy flux) q (absorbed) translucent slab q (transmitted) 1 absorptance

transmittance reflectance 20 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 T 4 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 Wiens 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 4 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. The net heat transfer is then: Qnet A (Tw4 T4 ) For small T the term (Tw4-T4) can be 3 approximated as 4T (Tw T ) and Qnet A hr T with hr as an effective radiation heat transfer coefficient. Qnet T q qw Tw A 22 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 Wall heat flux calculation heat flux : q h f (Tw T f ) qrad for laminar flows : h f kf y for turbulent flows : h f follows from correlations describing the thermal boundary layer profile external wall radiation : qrad ext (T4 Tw4 ) Tw wall temperatur e T f fluid cell temperatur e T user specified temperatur e ext emissivity of external wall surface Stefan Boltzmann constant h f fluid side local heat transfer coefficien t

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 Q A hr T Radiation: 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 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. latent heat heat capacity 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

ow Fl D cti ire

Injectors Original design: Bundle of tubes as shown. Repeated geometry. 3 different baffles, A, C, and D. Reactant injectors between baffles A and D.

on

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 halfcylindrical surface upstream of the first baffle. 37 Heat exchanger - flow pattern Recirculation loops D

A C A D A C A D Bundle of tubes Flow Direction (Shaded area) Injection points Compartments with low flow Compartment mainly served by nearest upstream injector Non-uniform flow distribution means low efficiency

38 Heat exchanger - modifications

Baffle A Relocated Injectors Design modifications: Shorter baffles C . Relocated and rotated injectors. Modified (Shorter) Baffle C ow Fl

D on cti e ir

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