《高等传热学》课程教学资源（文献资料）Numerical simulation of thermomagnetic convection of air in a porous square enclusure under a magnetic quadrupole field using LTNE models

of Heat and(015)1 Contents lists available at ScenceDirect International Journal of Heat and Mass Transfer ELSEVIER journal homepage:www.elsevier.com/locate/ijhmt Numerical simulation of thermomagnetic convection of air in a porous square enclosure under a magnetic quadrupole field using LTNE models Changwei Jiang,Er Shi,Zhangmao Hu,Xianfeng Zhu,Nan Xie ARTICLE INFO ABSTRACT convection of air in a tw 23m March 2015 in this paper. ddle pla of decrease at first and then i 2015 Elsevier Ltd.All rights reserved. 1.Introduction etic field the heattransfer and uid fow has There are electro ic packaging purific al gro porousmedium using the homotopy analysis method. sfe dary con ditions have an important influence on the ngth and field have a strons porous square cavities.Chankim et al.analyzed theoretically porous medium4.Sathiyam 5analWzedtheconveCtw gnetic field.Magnetohydr in E-ma 7oo2aS2aa2&

Numerical simulation of thermomagnetic convection of air in a porous square enclosure under a magnetic quadrupole field using LTNE models Changwei Jiang ⇑ , Er Shi, Zhangmao Hu, Xianfeng Zhu, Nan Xie School of Energy and Power Engineering, Changsha University of Science and Technology, Changsha 410114, China Key Laboratory of Efficient and Clean Energy Utilization, College of Hunan Province, Changsha 410114, China article info Article history: Received 21 August 2014 Received in revised form 3 March 2015 Accepted 23 July 2015 Available online 6 August 2015 Keywords: Thermomagnetic convection Numerical simulation Porous media Magnetic quadrupole field Magnetic force abstract In this paper, thermomagnetic convection of air in a two-dimensional porous square enclosure under a magnetic quadrupole field has been numerically investigated. The scalar magnetic potential method is used to calculate the magnetic field. A generalized model, which includes a Brinkman term, a Forcheimmer term and a nonlinear convective term, is used to solve the momentum equations and the energy for fluid and solid are solved with the local thermal non-equilibrium (LTNE) models. The results are presented in the form of streamlines and isotherms and local and average Nusselt numbers. Numerical results are obtained for a range of the magnetic force parameter from 0 to 100, the Darcy num￾ber from 105 to 101 and dimensionless solid-to-fluid heat transfer coefficient from 1 to 1000. The results show that the magnetic force number, Darcy number, Rayleigh number and dimensionless solid-to-fluid heat transfer coefficient have significant effect on the flow field and heat transfer in a por￾ous square enclosure. The flow characteristics presents two cellular structures with horizontal symmetry about the middle plane of the enclosure and the Nusselt numbers are increased as the magnetic force number increases under the non-gravitational condition. The average Nusselt number respects the trend of decrease at first and then increases when the magnetic force number increases under gravitational condition. The non-equilibrium effect on fluid phase temperature and solid phase temperature gradually reduces with the increase of value of H. 2015 Elsevier Ltd. All rights reserved. 1. Introduction Natural convection heat transfer in porous enclosure is widely used in many industrial applications such as cooling of electronic devices, solar collectors, heat exchangers and so on. There are many open literature related to natural convection in porous enclosures [1–3]. Ellahi et al. [4,5] have analyzed the influence of variable viscosity and viscous dissipation on the non-Newtonian flow in porous medium using the homotopy analysis method. Numerical investigation of natural convection within porous square enclosures for various thermal boundary conditions has been done by Ramakrishna et al. [6]. It is found that thermal boundary conditions have an important influence on the flow and heat transfer characteristics during natural convection within porous square cavities. Chankim et al. [7] analyzed theoretically the onset of convection motion in an initially quiescent, horizontal isotropic porous layer. Magnetic field effect on the heat transfer and fluid flow has received much attention in recent years due to its importance in electronic packaging, purification of molten metals, crystal growth in liquids and many others [8–11]. Saleh et al. [12] analyzed the effect of a magnetic field on steady convection in a trapezoidal enclosure filled with a fluid-saturated porous medium by the finite difference method. Grosan et al. [13] examined the effects of a magnetic field and internal heat generation on natural convection heat transfer in an inclined square enclosure filled with a fluid-saturated porous medium. It was shown that both the strength and inclination angle of the magnetic field have a strong influence on convection modes. Nield studied MHD convection in porous medium [14]. Sathiyamoorthy [15] analyzed the convective heat transfer in a square cavity filled with porous medium under a magnetic field. Magnetohydrodynamic natural convection in a rectangular cavity under a uniform magnetic field at different angles with respect to horizontal plane has been investigated by Yu et al. [16]. They concluded that the heat transfer is not only determined by the strength of the magnetic field, but also http://dx.doi.org/10.1016/j.ijheatmasstransfer.2015.07.103 0017-9310/ 2015 Elsevier Ltd. All rights reserved. ⇑ Corresponding author at: School of Energy and Power Engineering, Changsha University of Science and Technology, Changsha 410114, China. Tel./fax: +86 731 85258409. E-mail address: cw_jiang@163.com (C. Jiang). International Journal of Heat and Mass Transfer 91 (2015) 98–109 Contents lists available at ScienceDirect International Journal of Heat and Mass Transfer journal homepage: www.elsevier.com/locate/ijhmt

99 Nomenclature density.bo=Br(T) n coordinates 4BB nanent magnets(T) dimensionless fluid temperature.= pefficient (Wm-2K-) solid temperature. netic nsit wm aticviscos average Nusselt number mass magnenss m e(K) reference value velocity vector solid Wakayama and cow orkers [26-28]studi d the ransfe uation for nagnetic conve ing a m milar to ic d app cubi 13-34 ed Finite Eleme t Method (C M). The results ind ases ile an opposite trend d that the magn e can be us un erical availabe from Neody mium-Iro out hat t f air or they found tha natural magne and depends on the relati rate of pa dients theoretically and experim ntally and found that th on of par s.How wa onal fo ces on the natura cated in n in n. e stability of p uid layer by the action im u diamagnetic field.Gray et al.125]obtained the nu erica number and magnet the esult ifcant effect on the flow field and heat transfer in a paramagnet

influenced by the inclination angle. Especially, when the aspect ratio is less or more than 1, it is found that the inclination angle plays a great role on flow and heat transfer. The natural convection heat transfer of Cu–water nanofluid in a cold outer circular enclosure containing a hot inner sinusoidal cir￾cular cylinder in the presence of horizontal magnetic field is inves￾tigated numerically by Sheikholeslami et al. [17] using the Control Volume based Finite Element Method (CVFEM). The results indi￾cate that in the absence of magnetic field, enhancement ratio decreases as Rayleigh number increases while an opposite trend is observed in the presence of magnetic field. In recent years, with the development of a superconducting magnet providing strong magnetic induction of 10 T or more, the natural convection of paramagnetic fluid like oxygen gas and air under magnetic field has been an interesting research topic inves￾tigated by many researchers [18,19]. The effect of the magnetic buoyancy force on the convection of paramagnetic fluids was first reported by Braithwaite et al. [20], they found that the effect of magnetic buoyancy force on convection depends on the relative orientation of the magnetic force and the temperature gradient. Carruthers and Wolfe [21] studied the thermal convection of oxy￾gen gas in a rectangular container with thermal and magnetic field gradients theoretically and experimentally, and found that the magnetic buoyancy force cancel out the influence of gravity buoy￾ancy force when rectangle enclosure heated from one vertical wall and cooled from opposing wall was located in horizontal magnetic field with vertical magnetic field gradient. Huang et al. [22–24] studied the stability of paramagnetic fluid layer by the action of a non-uniform magnetic field and the thermomagnetic convection of the diamagnetic field. Gray et al. [25] obtained the numerical similarity solutions for the two-dimensional plums driven by the interaction of a line heat source and a non-uniform magnetic field. Wakayama and coworkers [26–28] studied the magnetic control of thermal convection. Tagawa’s group [29–31] derived a model equation for magnetic convection using a method similar to the Boussinesq approximation and studied natural convection of para￾magnetic, diamagnetic and electrically conducting fluids in a cubic enclosure with thermal and magnetic field gradients at different thermal boundary. Ozoe and co-workers [18,32–34] studied natu￾ral convection of paramagnetic and diamagnetic fluids in a cylinder under gradient magnetic field at different thermal boundary and found that the magnetic force can be used to control heat transfer rate of paramagnetic and diamagnetic fluids. Yang et al. [35,36] numerically and experimentally investigated magnetothermal convection of air or oxygen gas in a enclosure by using the gradient magnetic field available from Neodymium–Iron–Boron permanent magnet systems, and pointed out that the enhancement or sup￾pression of magnetothermal convection of air or oxygen gas can be achieved by gradient magnetic field. Tomasz and co-workers [37–40] studied natural convection of paramagnetic fluids in a cubic enclosure under magnetic field by an electric coil and ana￾lyzed that the effect of inclined angle of electric coil, location of electric coil, Ra number on heat transfer rate of paramagnetic fluid. Above studies are concerned with the effect of magnetic force on natural convection of paramagnetic fluids. However, only a few studies are paid on the combined effects of both magnetic and gravitational forces on the natural convection of paramagnetic fluids in porous medium. Natural convection in an enclosure filled with a paramagnetic or diamagnetic fluid-saturated porous med￾ium under strong magnetic field was numerically investigated by Wang et al. [41–43]. Considering the effect of Darcy number, Rayleigh number and magnetic force number, the results of numerical investigation showed that the magnetic force had a sig￾nificant effect on the flow field and heat transfer in a paramagnetic Nomenclature b magnetic flux density (T) b0 reference magnetic flux density, b0 ¼ Br (T) B dimensionless magnetic flux Br magnetic flux density of permanent magnets (T) C C ¼ 1 þ 1 T0b cp fluid specific heat at constant pressure (J kg1 K1 ) Da Darcy number, j L2 fm magnetic force g gravitational acceleration (m s2 ) h solid-to-fluid heat transfer coefficient (W m2 K1 ) H dimensionless solid-to-fluid heat transfer coefficient, hL2 ekf H magnetic field intensity kf fluid thermal conductivity (W m1 K1 ) ks solid thermal conductivity (W m1 K1 ) L length of the enclosure, (m) Num average Nusselt number p pressure, Pa P dimensionless pressure p0 pressure at reference temperature, Pa p0 pressure difference due to the perturbed state, Pa Pr Prandtl number, Pr ¼ mf af Ra Rayleigh number, Ra ¼ gbðThTc ÞL3 af mf T0 T0 ¼ ThþTc 2 (K) Tc cold wall temperature (K) Tf fluid temperature (K) Th hot wall temperature (K) Ts solid temperature (K) u,v velocity components (m s1 ) U,V dimensionless velocity components U velocity vector x,y Cartesian coordinates X,Y dimensionless Cartesian coordinates Greek symbols af fluid thermal diffusivity (m s1 ) b thermal expansion coefficient (K1 ) c dimensionless magnetic strength parameter, c ¼ v0b2 0 lmgL hf dimensionless fluid temperature, hf ¼ Tf T0 ThTc hs dimensionless solid temperature, hs ¼ TsT0 ThTc e porosity l0 magnetic permeability of vacuum (H m1 ) lm magnetic permeability (H m1 ) lf fluid kinematic viscosity (kg m1 s 1 ) mf fluid dynamic viscosity (m2 s1 ) qf fluid density (kg m3 ) qs solid density (kg m3 ) v mass magnetic susceptibility (m3 kg1 ) v0 reference mass magnetic susceptibility (m3 kg1 ) vm volume magnetic susceptibility j Permeability (m2 ) K dimensionless thermal conductivity, K ¼ ekf ð1eÞks um scalar magnetic potential Subscripts 0 reference value c cold f fluid h hot s solid C. Jiang et al. / International Journal of Heat and Mass Transfer 91 (2015) 98–109 99

100 amr91(2015)9g-109 :Pyis the fluid density.kgm is the mass earfuture Thus the sty of on nat- ndudes the rce research and Even tho can be pres engineering app ough fluids in a rous cubic ncoimthanceantcalcolharet (2) n which.U is the velocity ye The on natural com en erical ing At the refere convectio ord rynumbedmeolid-t-heat 0=-+2+p (3) 2.Physical model where is the pressure at reference temperature.Pa:is the ility at reference temperature.mkg subtracting ge oesuareencos le the ot tw +22+-咖8 (4 the present study.the size of the rnethand 3.Mathematical formulation 4-4。+(）西-10+ (5) n饮pmeg 4=m+().-+ x-号 (7) x-9。=(器x-p)-+ N =w+)-+ 张(+)-+-8 1.Physical model and

or diamagnetic fluid-saturated porous medium. The application of strong magnetic field for porous medium may be found in the field of medical treatment, metallurgy, materials processing, combus￾tion. There may be plenty of applications in engineering field in the near future. Thus, the study of effect of magnetic force on nat￾ural convection in porous media is very important for both scien￾tific research and engineering application. Even though earlier studies on the natural convection of paramagnetic or diamagnetic fluids in a porous cubic enclosure with an electrical coil have been carried out by wang et al., a detailed investigation of natural con￾vection of air in a porous square enclosure under a magnetic quad￾rupole field using LTNE models has yet to appear in the literature. The present paper represents the results of a numerical investi￾gation on natural convection of air in a porous square enclosure under a magnetic quadrupole field using local thermal no-equilibrium (LTNE) models. The numerical investigation is car￾ried out for different governing parameters such as the magnetic force number, Darcy number, and dimensionless solid-to-fluid heat transfer coefficient etc. 2. Physical model The schematic of the system under consideration is shown in Fig. 1. The system consists of a porous square enclosure which is kept in a horizontal position and four permanent magnets which generate a magnetic field. The porous square enclosure filled with air is heated isothermally from left-hand side vertical wall and cooled isothermally from opposing wall while the other two walls are thermally insulated. The gravitational force acts in the Y direction. In the present study, the size of the enclosure L, the size of the permanent magnet L1 and the distance of permanent mag￾nets L2 are 0.024 m, 0.02 m and 0.03 m respectively. 3. Mathematical formulation 3.1. Governing equations The assumptions in the model are as follows: the fluid is consid￾ered to be steady, incompressible, and Newtonian fluid. Both the viscous heat dissipation and magnetic dissipation are assumed to be negligible. According to Braithwaite et al. [20], the magnetic force can be given as follows: fm ¼ vm 2lm rb2 ¼ qf v 2lm rb2 ð1Þ where, fm is the magnetic force; vm is the volumetric magnetic sus￾ceptibility; lm is the magnetic permeability, H m1 ; b is the mag￾netic flux density, T; qf is the fluid density, kg m3 ; v is the mass magnetic susceptibility, m3 kg1 . The Navier–Stokes equation which includes the magnetic force can be presented as: qf e2 U rU ¼ rp lf j U þ lf e r2 U qf e3=2 1:75 ffiffiffiffiffiffiffiffiffi 150 p jUjU ffiffiffi j p þ qf v 2lm rb2 þ qf g ð2Þ In which, U is the velocity vector; p is the pressure, Pa; lf is the fluid kinematic viscosity, kg m1 s1 ; g is the gravitational acceler￾ation, m s2 ; j is the permeability, m2 . At the reference state of the isothermal state, there will be no convection. Therefore, Eq. (2) becomes as follows. 0 ¼ rp0 þ qf0v0 2lm rb2 þ qf 0g ð3Þ where p0 is the pressure at reference temperature, Pa; qf0 is the fluid density at reference temperature, kg m3 ; v0 is the mass mag￾netic susceptibility at reference temperature, m3 kg1 , subtracting (3) from (2) gives: qf e2 U rU ¼ rp0 lf j U þ lf e r2 U qf e3=2 1:75 ffiffiffiffiffiffiffiffiffi 150 p jUjU ffiffiffi j p þ ðqf v qf 0v0Þ 2lm rb2 þ ðqf qf0Þg ð4Þ where: p ¼ p0 þ p0 , p0 is the pressure difference due to the perturbed state, Pa. Because qf and v are function of temperature, according to Taylor expansion method, qfv and qf can be respectively indicated as: qf v ¼ ðqf vÞ 0 þ @ðqf vÞ @Tf  0 ðTf T0Þþ ð5Þ qf ¼ qf 0 þ @qf @Tf  0 ðTf Tf 0Þþ ð6Þ For the paramagnetic fluid air, the mass magnetic susceptibility is inverse proportion to absolute temperature, according to the Curie’s law: v ¼ m Tf ð7Þ where m is the constant value; Tf is the fluid phase temperature, K; T0 = (Th + Tc)/2, K; Subscripts 0, h, c represent reference value, hot and cold respectively. So Eq. (5) can be written as: qf v ðqfvÞ 0 ¼ @qf @Tf v qf v Tf  0 ðTf T0Þþ ¼ qf bv qf v Tf  0 ðTf T0Þþ ¼ qf0v0b 1 þ 1 T0b ðTf T0Þþ ð8Þ The higher order small amount in Eq. (8) is omitted and gener￾ated into Eq. (4), Eq. (4) becomes as follows. 1 e2 U rU ¼ rp0 qf0 lf qf 0j U þ lf qf 0e r2 U 1 e3=2 1:75 ffiffiffiffiffiffiffiffiffi 150 p jUjU ffiffiffi j p v0b 2lm 1 þ 1 T0b ðTf T0Þrb2 þ bðTf T0Þg ð9Þ where b is thermal expansion coefficient, K1 . Fig. 1. Physical model and coordinate system. Accordingly, the governing equations can be written as: 100 C. Jiang et al. / International Journal of Heat and Mass Transfer 91 (2015) 98–109

C.Jiang et al.Int nd Mass Tranfer91(2015)98-109 verage Nusselt numbers with [45]at various values of Br. 25 Continuity equation: 贺+-0 Momentum equation △T 袋》贵货-源品 Present results Relative error +兰快+ (+制 0 02 Yang et al. Present work with 35](left:temperature field and right:velocity field)

Continuity equation: @u @x þ @v @y ¼ 0 ð10Þ Momentum equation: qf e2 u @u @x þ v @u @y  ¼ @p @x lf j u qf 1:75 ffiffiffiffiffiffiffiffiffi 150 p ðu2 þ v2Þ 1=2 ffiffiffi j p u e3=2 þ lf e @2 u @x2 þ @2 u @y2 ! 1 þ 1 T0b  v0bðTf T0Þ 2lm @ðb2 Þ @x ð11Þ Table 1 Comparison of the average Nusselt number Num for different grid resolution at e = 0.5, Pr = 0.71, H = 10, K = 10, cRa = 1  106 , Da = 1  103 under non-gravita￾tional condition. Grid dimension Num 30  30 2.0635 40  40 2.0519 50  50 2.0478 60  60 2.0460 70  70 2.0454 Table 2 Comparison of present results with [35]. DT Num Yang et al. [35] Present results Relative error/% 1 1.003 1.003 0 10 1.214 1.244 2.47 50 2.120 2.166 2.17 Fig. 2. The comparison of present results with [35] (left: temperature field and right: velocity field). Table 3 Comparison of average Nusselt numbers with [45] at various values of Br. Br Num Song et al. [45] Present results 0.0 T 4.520 4.525 0.5 T 4.506 4.524 1.5 T 4.352 4.360 2.5 T 3.898 3.923 C. Jiang et al. / International Journal of Heat and Mass Transfer 91 (2015) 98–109 101

(b) Fig.3.(a)Distribution of gradient of squ ion (VB"L (b)vec (-CRaPrB22)under the non-gravitational condito Ra=1x10 Ra=5×1 yRa=1×10 7Ra=1x10 rnthe(e)m iddle)and solid phase temperature (right)for Da-10-3,H-10 and 2-05 under the non (装+岁需柴-照”品 Fluid phase energy equation VK 偿 w要+-（德-功间 Solid phase energy equation: -1+) 0=1-9k(++- (14) +PygB(T;-To) (12)

qf e2 u @v @x þ v @v @y  ¼ @p @y lf j v qf 1:75 ffiffiffiffiffiffiffiffiffi 150 p ðu2 þ v2Þ 1=2 ffiffiffi j p v e3=2 þ lf e @2 v @x2 þ @2 v @y2 ! 1 þ 1 T0b  v0bðTf T0Þ 2lm @ðb2 Þ @y þ qf gbðTf T0Þ ð12Þ Fluid phase energy equation: ðqcpÞf u @Tf @x þ v @Tf @y  ¼ ekf @2 Tf @x2 þ @2 Tf @y2 ! þ hðTs TfÞ ð13Þ Solid phase energy equation: 0 ¼ ð1 eÞks @2 Ts @x2 þ @2 Ts @y2 ! þ hðTf TsÞ ð14Þ Fig. 3. (a) Distribution of gradient of square magnetic induction (rB2 ), (b) vectors of the magnetizing force (CcRaPrhfrB2 /2) under the non-gravitational condition. Fig. 4. Effect of cRa number on the streamlines (left), fluid phase temperature (middle) and solid phase temperature (right) for Da = 103 , H = 10 and e = 0.5 under the non￾gravitational convection. 102 C. Jiang et al. / International Journal of Heat and Mass Transfer 91 (2015) 98–109