《高等传热学》课程教学资源（文献资料）Numerical simulation of thermomagnetic convection of air in a porous cubic enclosure with an electric coil inclined in general orithentations using a LTNE model.pdf

Journal of Porous Media, 18 (11): 1059–1074 (2015) NUMERICAL SIMULATION OF THERMOMAGNETIC CONVECTION OF AIR IN A POROUS CUBIC ENCLOSURE WITH AN ELECTRIC COIL INCLINED IN GENERAL ORIENTATIONS USING A LTNE MODEL Changwei Jiang,∗ Er Shi, Xianfeng Zhu, & Nan Xie School of Energy and Power Engineering, Changsha University of Science and Technology, Changsha 410114, China; and Key Laboratory of Efficient and Clean Energy Utilization, College of Hunan Province, Changsha 410114, China ∗Address all correspondence to Changwei Jiang E-mail: cw jiang@163.com Original Manuscript Submitted: 11/9/2014; Final Draft Received: 2/25/2015 The natural convection heat transfer of paramagnetic fluid in porous media can be controlled by the gradient magnetic field. Thermomagnetic convection of air in a porous cubic enclosure with an inclined coil around the Y and Z axes is numerically investigated in the presence or absence of a gravity field. The Biot–Savart law is used to calculate the magnetic field. The Darcy model is used to solve the momentum, and the energies for fluid and solid are solved with the local thermal non-equilibrium (LTNE) models. The flow and temperature fields for the thermomagnetic convection of air are presented and the average Nusselt number on the hot wall is calculated and compared. The results show that both the magnetic force and the coil inclination have a significant effect on the flow field and heat transfer in a porous cubic enclosure, and the thermomagnetic convection heat transfer of air can be enhanced or controlled by applying a gradient magnetic field. KEY WORDS: thermomagnetic convection, numerical simulation, magnetic force, inclined electric coil, porous media 1. INTRODUCTION Natural convection in porous media is a subject of prime importance due to its use in recent engineering and scientific applications, such as in heat exchangers, solar receivers, geothermal exploitation, cooling of electronic devices, etc. Various researchers have carried out experimental and numerical investigations to study fluid flow and heat transfer in porous enclosures (Khalili and Shivakumara, 2003; Vafai, 2005; Pourshaghaghy et al., 2007; Bhardwaj and Dalal, 2013; Ellahi et al., 2013). Varol et al. (2006) analyzed the steady-state free-convection heat transfer in a porous medium–filled right-angled triangular enclosure. They found that heat transfer increases with the decreasing aspect ratio. Sankar et al. (2011) investigated natural convection flows in a vertical annulus filled with a fluid-saturated porous medium when the inner wall was subject to discrete heating. Numerical investigation of natural convection within porous square enclosures for various thermal boundary conditions was done by Ramakrishna et al. (2013), in which it was found that the thermal boundary conditions have an important influence on the flow and heat transfer characteristics during natural convection within porous square cavities. Kim and Choi (2013) theoretically analyzed the onset of convection motion in an initially quiescent, horizontal isotropic porous layer. Alloui and Vasseur (2013) investigated natural convection in a shallow porous cavity filled with a binary fluid. Natural convection under the effect of a magnetic field is of great importance in many industrial applications, 1091–028X/15/$35.00 ⃝c 2015 by Begell House, Inc. 1059

1060 Jiang et al. NOMENCLATURE B dimensionless magnetic flux b magnetic flux density (T) b0 reference magnetic flux density, (= µmi/L) (T) C = 1+ (1/T0β) cp fluid specific heat at constant pressure (J kg−1 K −1 ) Da Darcy number (= κ/L2 ) g gravitational acceleration (m s−2 ) H dimensionless solid-to-fluid heat transfer coefficient (= hL2/εkf) h solid-to-fluid heat transfer coefficient (W m−2 K −1 ) i electric current in a coil (A) kf fluid thermal conductivity (W m−1 K −1 ) ks solid thermal conductivity (W m−1 K −1 ) L length of a cubic enclosure (m) Num average Nusselt number P dimensionless pressure p pressure (Pa) Pr Prandtl number (= νf/αf) R dimensionless position vector r position vector (m) r radius of the coil (m) Ra Rayleigh number [= gβ (Th − Tc)L 3/αfνf] S dimensionless tangential element of a coil s tangential element of a coil (m) Tc cold wall temperature (K) Tf fluid temperature (K) Th hot wall temperature (K) Ts solid temperature (K) T0 reference temperature [= (Th + Tc)/2] (K) U,V ,W dimensionless velocity components u, v, w velocity components (ms−1 ) X,Y ,Z dimensionless Cartesian coordinates x, y, z Cartesian coordinates xEuler rotation angle around the X axis (◦ ) yEuler rotation angle around the Y axis (◦ ) zEuler rotation angle around the Z axis (◦ ) Greek Symbols αf fluid thermal diffusivity (m2 s −1 ) β thermal expansion coefficient (K−1 ) γ dimensionless magnetic strength parameter (= χ0b 2 0/µmgL) ε porosity θf dimensionless fluid temperature [= (Tf – T0)/(Th – Tc)] θs dimensionless solid temperature [= (Ts – T0)/(Th – Tc)] κ permeability (m2 ) Λ dimensionless thermal conductivity [= εkf /(1 – ε)ks] µf fluid kinematic viscosity (kg m−1 s −1 ) µm magnetic permeability (H m−1 ) µ0 magnetic permeability of vacuum (H m−1 ) νf fluid dynamic viscosity (m2 s −1 ) ρf fluid density (kg m−3 ) ρs solid density (kg m−3 ) χ mass magnetic susceptibility (m3 kg−1 ) χm volume magnetic susceptibility χm0 reference volume magnetic susceptibility χ0 reference mass magnetic susceptibility (m3 kg−1 ) Subscripts c cold f fluid h hot m magnetic s solid 0 reference value such as in solar energy technology, crystal growth, and liquid metal cooling blankets used in fusion reactors. During the past three decades, numerous theoretical and experimental studies on thermomagnetic convection in enclosures have been investigated (Bian et al., 1996; Sarris et al., 2005; Pirmohammadi et al., 2009; Grosan et al., 2009; Sathiyamoorthy and Chamkha, 2010). The natural convection heat transfer of Cu–water nanofluid in a cold outer circular enclosure containing a hot inner sinusoidal circular cylinder in the presence of a horizontal magnetic field was numerically investigated by Sheikholeslami et al. (2012) using the control volume-based finite-element method. The results indicated that in the absence of a magnetic field, the heat transfer enhancement ratio deJournal of Porous Media

Numerical Simulation of Thermomagnetic Convection of Air 1061 alues and found that the magnetic force an be used to on natural convection of non-Newtonian power-law fluids perimentally investigated magnetothermal convection of L2014 ed the fo me dition using the lattice Boltzmann method.They found or oxygen gas can be achieved by the gradient magnetic field.Bednarz et al. (2004,2005. 009)stu ed the ncrease wit r all phas In recent years,with the development of a supercon- alyzed the effect of the inclined angle of the electric co tic induc he location of the electric coil,and the ka number nd h been an interesting research topic investigated by many a paramagnetic field placed in a micro-gravity condition researchers(Qiet al,2001:Tagawa et al,2002).The ef and under a uniform vertical gradient magnetic field in an of the cold side buoyancy force on convection depends on the relative ori- of paramagnetic and diamagnetic fluids.However,only 50 the th ontainer with thermal and field Wang et al.(2007)and Zeng et al.(2007,2009)nu gradients and found that the magnetic buoyancy force onvection in an enclosure out mence of the buoyanc orc gnetic edraG cal wall and cooled from an ing wall ，was loc (Da Ra and mag with a vertical magnetic force numbers,the results of their numerical investig d gradient. Huang et a magpeticforc has a signincant e tion of the diamagnetic field.Gray et al.(2001)obtained plication of a strong magnetic field to a porous me numerical similarity solutions for two-dimensional plums may be found in the fields of medical treatment. met en by the i .200 and allurgy,ma ssing,and com and t and Wakavama (2002)studied the of engineering o erations.Thus.the study of the efTect 0)0) of magnetic force on natural convection in porous me. u et al.(2003)derive dia is impo ant n and enginee the Bor magn ngap pres ral convection of paramagnetic.diamagnetic.and electri ooled from the sidewalls with an electric coil inclined cally conducting fluids n a cubic nclosure with therma ound the Y and Z axes is numeric ally investigated.The and magnene ld gra now and heat magn (2005)studied the natural com ction of nation angle of an electric coil around the y axis paramagnetic and diamagnetic fluids in a cylinder under and Z axis ()Ra number.Da number.and magnetic a gradient magnetic field at different thermal boundary force parameter y. Volume 18,Number 11,2015

Numerical Simulation of Thermomagnetic Convection of Air 1061 creases as the Rayleigh (Ra) number increases, while an opposite trend is observed in the presence of a magnetic field. Kefayati (2014) studied the effect of a magnetic field on natural convection of non-Newtonian power-law fluids in a cavity with a linearly heated wall using the finitedifference lattice Boltzmann method (FDLBM). Mejri et al. (2014) analyzed the effects of a magnetic field on nanofluid flow in a cavity with a sinusoidal boundary condition using the lattice Boltzmann method. They found that the heat transfer and fluid flow decline with the Hartmann number and increase with Ra number for all phase deviations. In recent years, with the development of a superconducting magnet providing strong magnetic induction of 10 T or more, the natural convection of paramagnetic fluids such as oxygen gas or air under a magnetic field has been an interesting research topic investigated by many researchers (Qi et al., 2001; Tagawa et al., 2002). The effect of the magnetic buoyancy force on the convection of paramagnetic fluids was first reported by Braithwaite et al. (1991). They found that the effect of the magnetic buoyancy force on convection depends on the relative orientation of the magnetic force and the temperature gradient. Carruthers and Wolfe (1968) theoretically and experimentally studied the thermal convection of oxygen gas in a rectangular container with thermal and magnetic field gradients and found that the magnetic buoyancy force cancels out the influence of the gravity buoyancy force when the rectangular enclosure, heated from one vertical wall and cooled from an opposing wall, was located in a horizontal magnetic field with a vertical magnetic field gradient. Huang et al. (1998a,b,c); studied the stability of a paramagnetic fluid layer by the action of a nonuniform magnetic field and the thermomagnetic convection of the diamagnetic field. Gray et al. (2001) obtained numerical similarity solutions for two-dimensional plums driven by the interaction of a line heat source and a nonuniform magnetic field. Qi et al. (1999, 2001) and Wang and Wakayama (2002) studied the magnetic control of thermal convection. Kaneda et al. (2002), Shigemitsu et al. (2003), and Akamatsu et al. (2003) derived a model equation for magnetic convection using a method similar to the Boussinesq approximation and studied the natural convection of paramagnetic, diamagnetic, and electrically conducting fluids in a cubic enclosure with thermal and magnetic field gradients at different thermal boundary values. Akamatsu et al. (2004), Filar et al. (2005), and Fornalik et al. (2005) studied the natural convection of paramagnetic and diamagnetic fluids in a cylinder under a gradient magnetic field at different thermal boundary values and found that the magnetic force can be used to control the heat transfer rate of paramagnetic and diamagnetic fluids. Yang et al. (2003, 2004) numerically and experimentally investigated magnetothermal convection of air or oxygen gas in an enclosure by using the gradient magnetic field available from neodymium–iron–boron permanent magnet systems, and pointed out that enhancement or suppression of magnetothermal convection of air or oxygen gas can be achieved by the gradient magnetic field. Bednarz et al. (2004, 2005, 2008, 2009) studied the natural convection of paramagnetic fluids in a cubic enclosure under a magnetic field by an electric coil and analyzed the effect of the inclined angle of the electric coil, the location of the electric coil, and the Ra number on the heat transfer rate of the paramagnetic fluid. Saha (2013) studied thermomagnetic convection and heat transfer of a paramagnetic field placed in a micro-gravity condition and under a uniform vertical gradient magnetic field in an open square cavity with three cold sidewalls. The aforementioned studies were concerned with the effect of the magnetic force on the natural convection of paramagnetic and diamagnetic fluids. However, only a few studies have paid attention to the combined effects of both magnetic and gravitational forces on the natural convection of a paramagnetic fluid in a porous medium. Wang et al. (2007) and Zeng et al. (2007, 2009) numerically investigated natural convection in an enclosure filled with a paramagnetic or diamagnetic fluid–saturated porous medium under a strong magnetic field. Considering the effect of the Darcy (Da), Ra, and magnetic force numbers, the results of their numerical investigations showed that the magnetic force has a significant effect on the flow field and heat transfer in a paramagnetic or diamagnetic fluid–saturated porous medium. The application of a strong magnetic field to a porous medium may be found in the fields of medical treatment, metallurgy, materials processing, and combustion, and there may be plenty of applications in the near future in the field of engineering operations. Thus, the study of the effect of magnetic force on natural convection in porous media is important for both scientific research and engineering applications. In the present study, the thermomagnetic convection of air in a porous cubic enclosure heated and cooled from the sidewalls with an electric coil inclined around the Y and Z axes is numerically investigated. The flow and heat transfer behaviors of thermomagnetic convection are studied in relation to the effect of the inclination angle of an electric coil around the Y axis (yEuler) and Z axis (zEuler), Ra number, Da number, and magnetic force parameter γ. Volume 18, Number 11, 2015

1062 Jiang et al. 2. PHYSICAL MODEL A schematic of the system under consideration is shown in Fig. 1. The system consists of a porous cubic enclosure, which is kept in a horizontal position, and a one-turn electric circle coil, which generates a magnetic field. The porous cubic enclosure filled with air is heated isothermally from the left-hand-side vertical wall and cooled isothermally from the opposing wall, while the other four walls are thermally insulated. The electric coil is set centrally and coaxially with the enclosure. Inclination of the coil in general directions is given by the combination of two angles: yEuler and zEuler. The coil is inclined first around the Y axis (Coil 1) and then around the Z axis (Coil 2). The gravitational force acts in the minus Z direction. In the present study, the size of the cubic enclosure (L) is 0.03 m and the radius (r) of the coil is 0.05 m. 3. MATHEMATICAL FORMULATION 3.1 Governing Equations The hypotheses in this model are as follows: the fluid is considered to be steady, incompressible, and Newtonian, and both the viscous heat dissipation and magnetic dissipation are assumed to be negligible. According to Braithwaite et al. (1991), the magnetizing force can be given as follows: fm = χm 2µm ∇b 2 = ρfχ 2µm ∇b 2 (1) FIG. 1: Schematic diagram of the physical model and coordinate system where fm is the magnetic force; χm is the volumetric magnetic susceptibility; µm is the magnetic permeability (Hm−1 ); b is the magnetic flux density (T); ρf is the fluid density (kgm−3 ); and χ is the mass magnetic susceptibility (m3 kg−1 ). The Navier–Stokes equation which includes the magnetic force can be presented as ρfU·∇U=−∇p+µf∇2U− µf κ U+ ρfχ 2µm ∇b 2+ρfg (2) where U is the velocity vector (u, v, w); p is the pressure (Pa); µf is the fluid kinematic viscosity (kg m−1 s −1 ); g is the gravitational acceleration (ms−2 ); and κ is the permeability (m2 ). At the reference state of the isothermal state, there will be no convection. Therefore, Eq. (2) becomes the following equation: 0 = −∇p0 + ρf0χ0 2µm ∇b 2 + ρf0g (3) where p0 is the pressure at the reference temperature (Pa); ρ0 is the fluid density at the reference temperature (kgm−3 ); and χ0 is the mass magnetic susceptibility at the reference temperature (m3 kg−1 ). Subtracting Eq. (3) from Eq. (2) gives ρfU · ∇U = −∇p ′+µf∇2U − µf κ U+ (ρfχ−ρf0χ0) 2µm × ∇b 2 + (ρf − ρf0) g (4) where p = p0+p ′ , and p ′ is the pressure difference due to the perturbed state (Pa). Because ρf and χ are a function of temperature, ρfχ and ρf can be indicated according to the Taylor expansion method, respectively, as follows: ρfχ = (ρfχ) 0 + ( ∂ (ρfχ) ∂Tf ) 0 (Tf − T0) + · · · (5) ρf = ρf0 + ( ∂ρf ∂T ) 0 (Tf − T0) + · · · (6) For the paramagnetic fluid, the mass magnetic susceptibility is inversely proportional to absolute temperature, according to Curie’s law: χ = m Tf (7) where m is a constant value; Tf is the fluid temperature (K); T0 = (Th + Tc)/2 (K); and subscripts 0, h, and c represent the reference value, hot, and cold, respectively. Thus, Eq. (5) can be written as follows: Journal of Porous Media

Numerical Simulation of Thermomagnetic Convection of Air 1063 ρfχ − (ρfχ) 0 = ( ∂ρf ∂Tf χ − ρf χ Tf ) 0 (Tf − T0) + · · · = ( −ρfβχ−ρf χ Tf ) 0 (Tf −T0)+ · · · (8) = −ρf0χ0β ( 1+ 1 T0β ) (Tf −T0)+ · · · The small higher-order amount in Eq. (8) is omitted and generated into Eq. (4), and Eq. (4) then becomes U · ∇U = − ∇p ′ ρf0 + µf ρf0 ∇2U − ν κ U − χ0β 2µm ( 1+ 1 T0β ) × (Tf − T0) ∇b 2 − β (Tf − T0) g (9) where β is the thermal expansion coefficient (K−1 ). For the sake of convenience, the subscripts of the parameters and the superscripts of pressure are omitted, and the governing equations can be written as follows: • Continuity equation ∂u ∂x + ∂v ∂y + ∂w ∂z = 0 (10) • Momentum equations u ∂u ∂x +v ∂u ∂y +w ∂u ∂z =− 1 ρf ∂p ∂x +νf ( ∂ 2u ∂x 2 + ∂ 2u ∂y 2 + ∂ 2u ∂z 2 ) − νf κ u − ( 1 + 1 T0β ) χ0β (Tf − T0) 2µm ∂ ( b 2 ) ∂x (11) u ∂v ∂x + v ∂v ∂y + w ∂v ∂z = − 1 ρf ∂p ∂y + νf ( ∂ 2v ∂x 2 + ∂ 2v ∂y 2 + ∂ 2v ∂z 2 ) − νf κ v− ( 1+ 1 T0β ) χ0β (Tf −T0) 2µm ∂ ( b 2 ) ∂y (12) u ∂w ∂x + v ∂w ∂y + w ∂w ∂z = − 1 ρf ∂p ∂z + νf ( ∂ 2w ∂x 2 + ∂ 2w ∂y 2 + ∂ 2w ∂z 2 ) − νf κ w − ( 1 + 1 T0β ) χ0β (Tf − T0) 2µm ∂ ( b 2 ) ∂z + gβ (Tf − T0) (13) • Fluid phase energy equation (ρcp) f ( u ∂Tf ∂x + v ∂Tf ∂y + w ∂Tf ∂z ) = εkf ( ∂ 2Tf ∂x 2 + ∂ 2Tf ∂y 2 + ∂ 2Tf ∂z 2 ) + h (Ts − Tf) (14) • Solid phase energy equation 0 = (1−ε) ks ( ∂ 2Ts ∂x 2 + ∂ 2Ts ∂y 2 + ∂ 2Ts ∂z 2 ) +h (Tf −Ts) (15) where x, y, and z are the Cartesian coordinates; u, v, and w are the velocity components (ms−1 ); αf is the fluid thermal diffusivity (ms−1 ); νf is the fluid dynamic viscosity (m2 s −1 ); (cp) f is the fluid specific heat at constant pressure (J kg−1 K −1 ); ε is the porosity; kf is the fluid thermal conductivity (W m−1 K −1 ); ks is the solid thermal conductivity (W m−1 K −1 ); h is the solid-to-fluid heat transfer coefficient (W m−2 K −1 ); Tf is the fluid phase temperature (K); and Ts is the solid phase temperature (K). Equations (10)–(15) can be non-dimensionalized as follows: • Continuity equation ∂U ∂X + ∂V ∂Y + ∂W ∂Z = 0 (16) • Momentum equations U ∂U ∂X + V ∂U ∂Y + W ∂U ∂Z = − ∂P ∂X + Pr( ∂ 2U ∂X2 + ∂ 2U ∂Y 2 + ∂ 2U ∂Z2 ) − Pr Da U − γRaPrθf C 2 ∂ ( B2 ) ∂X (17) U ∂V ∂X + V ∂V ∂Y + W ∂W ∂Z = − ∂P ∂Y + Pr( ∂ 2V ∂X2 + ∂ 2V ∂Y 2 + ∂ 2V ∂Z2 ) − Pr Da V − γRaPrθf C 2 ∂ ( B2 ) ∂Y (18) U ∂W ∂X +V ∂W ∂Y +W ∂W ∂Z =− ∂P ∂Z +Pr( ∂ 2W ∂X2 + ∂ 2W ∂Y 2 + ∂ 2W ∂Z2 ) − Pr Da W −γRaPrθf C 2 ∂ ( B2 ) ∂Z +RaPrθf (19) • Fluid phase energy equation U ∂θf ∂X + V ∂θf ∂Y + W ∂θf ∂Z = ∂ 2θf ∂X2 + ∂ 2θf ∂Y 2 + ∂ 2θf ∂Z2 + H (θs − θf) (20) • Solid phase energy equation Volume 18, Number 11, 2015