ISPH modeling of natural convection heat transfer with an analytical kernel renormalization factor

The objective of this study is to extend the attention of the incompressible smoothed particle hydrodynamics method (ISPH) in the heat transfer field. The ISPH method for the natural convection heat transfer under the Boussinesq approximation in various environments: pure-fluid, nanofluid, and non-Darcy porous medium is introduced. We adopted the improved analytical method for calculating the kernel renormalization factor and its gradient based on a quintic kernel function for the wall boundary treatment in the ISPH method. The proposed method requires no dummy particle layer to meet the impermeability condition and makes the heat flux over the wall boundary easy to implement. We performed four different numerical simulations of natural convection in cavities with increasing complexity in modeling and implementation: the natural convection in a square cavity with constant differentially heated wall temperature, natural convection with the heat flux from the bottom wall for a wide range of Rayleigh numbers, natural convection in a non-Darcy porous cavity fully filled with nanofluid in different flow regimes, and natural convection in a partially layered porous cavity. The results showed excellent agreement with results from literatures and the in-house P1–P1 finite element method code.     

Natural convection in an inclined rectangular porous cavity with uniform heat flux from the side

The problem of natural convection in an inclined rectangular porous layer enclosure is studied numerically. The enclosure is heated from one side and cooled from the other by a constant heat flux while the two other walls are insulated. The effect of aspect ratio, inclination angle and Rayleigh number on heat transfer is studied. It is found that the enclosure orientation has a considerable effect on the heat transfer. The negative orientation sharply inhibits the convection and consequently the heat transfer and a positive orientation maximizes the energy transfer. The maximum temperature within the porous medium can be considerably higher than that induced by pure conduction when the cavity is negatively oriented. The peak of the average Nusselt number depends on the Rayleigh number and the aspect ratio. The heat transfer between the two thermally active boundaries is sensitive to the effect of aspect ratio. For an enclosure at high or low aspect ratio, the convection is considerably decreased and the heat transfer depends mainly on conduction.     

Variable Permeability Effects on Natural Convection in a Vertical Porous Layer with Uniform Heat Flux from the Side

Transport in Porous Media - Natural convection heat transfer in a rectangular cavity filled with a saturated porous medium with variable permeability is investigated analytically and numerically....     

Stability analysis of soret-driven double-diffusive convection of Maxwell fluid in a porous medium

Stability analysis of double-diffusive convection for viscoelastic fluid with Soret effect in a porous medium is investigated using a modified-Maxwell-Darcy model. We use the linear stability analysis to investigate how the Soret parameter and the relaxation time of viscoelastic fluid effect the onset of convection and the selection of an unstable wavenumber. It is found that the Soret effect is to destabilize the system for oscillatory convection. The relaxation time also enhances the instability of the system. The effects of Soret coefficient and relaxation time on the heat transfer rate in a porous medium are studied using the nonlinear stability analysis, the variation of the Nusselt number with respect to the Rayleigh number is derived for stationary and oscillatory convection modes. Some previous results can be reduced as the special cases of the present paper.     

Vibrational convection in a saturated porous medium in the absence of thermal equilibrium between the solid and liquid phases

Thermal vibrational convection in a saturated porous medium is theoretically studied on the basis of a thermal nonequilibrium model, in accordance with which the temperatures of the porous medium and the saturating liquid can be different. In the high-frequency vibration approximation the averaged equations of convection are derived. The dependence of the vibration force direction on the interphase heat transfer coefficient and the vibration frequency is established. Vibrational convection in a cylindrical layer is studied. It is shown that, depending on the interphase heat transfer coefficient, the flows of two types differing in the liquid circulation direction can exist.     

正弦波温度边界下多孔介质方腔内非热平衡对流传热数值模拟

采用局部非热平衡模型,在方腔左侧壁面温度正弦波变化、右侧壁面温度均一的边界条件下,通过SIM-PLER算法数值研究了固体骨架发热多孔介质方腔内的稳态非达西自然对流,主要探讨了不同正弦波波动参数N及方腔的高宽比M/L对方腔内自然对流与传热的影响规律。计算结果表明:正弦波温度边界使得方腔内的流场出现了复杂的变化,流体及固体区域左侧壁面附近出现了周期性的正负变化的温度场分布,左侧壁面局部Nusselt数出现了周期性的震荡现象;存在一个最佳温度波动参数N=1,此时多孔介质方腔内的整体散热量达到最大值;增加方腔高宽比会显著地削弱方腔内的自然对流传热过程,小高宽比也会在一定的程度上削弱多孔介质方腔内的对流传热。     

Non-Darcy double-diffusive natural convection in axisymmetric fluid saturated porous cavities

Double-diffusive natural convection in a fluid saturated porous medium has been investigated using the finite element method. A generalised porous medium model is used to study both Darcy and non-Darcy flow regimes in an axisymmetric cavity. Results indicate that the Darcy number should be a separate parameter to understand flow characteristics in non-Darcy regime. The influence of porosity on heat and mass transfer is significant and the transport rates may differ by 25% or more, at higher Darcy and Rayleigh numbers. When compared with the Darcy and other specialised models of Brinkman and Forchheimer, the present generalised model predicts the least heat and mass transfer rates. It is also observed that an increase in radius ratio leads to higher Nusselt and Sherwood numbers along the inner wall.     

Natural Convection in a Cavity Filled with Porous Layers on the Top and Bottom Walls

Natural convection in a partially filled porous square cavity is numerically investigated using SIMPLEC method. The Brinkman-Forchheimer extended model was used to govern the flow in the porous medium region. At the porous-fluid interface, the flow boundary condition imposed is a shear stress jump, which includes both the viscous and inertial effects, together with a continuity of normal stress. The thermal boundary condition is continuity of temperature and heat flux. The results are presented with flow configurations and isotherms, local and average Nusselt number along the cold wall for different Darcy numbers from 10⁻¹ to 10⁻⁶, porosity values from 0.2 to 0.8, Rayleigh numbers from 10³ to 10⁷, and the ratio of porous layer thickness to cavity height from 0 to 0.50. The flow pattern inside the cavity is affected with these parameters and hence the local and global heat transfer. A modified Darcy–Rayleigh number is proposed for the heat convection intensity in porous/fluid filled domains. When its value is less than unit, global heat transfer keeps unchanged. The interfacial stress jump coefficients β ₁ and β ₂ were varied from  −1 to +1, and their effects on the local and average Nusselt numbers, velocity and temperature profiles in the mid-width of the cavity are investigated.     

基于蒙特卡罗随机有限元方法的随机多孔介质内流体自然对流不确定性研究

为分析孔隙率不确定性对多孔介质方腔内自然对流换热的影响,发展了一种基于KL（Karhunen-Loeve展开）-蒙特卡罗随机有限元算法的随机多孔介质内自然对流不确定性分析数理模型及有限元数值模拟程序框架。通过K-L展开及基于拉丁抽样法生成多孔介质孔隙率随机实现,并耦合多孔介质自然对流有限元程序,进行随机多孔介质内自然对流传热数值模拟,得出了多孔介质内流场与温度场平均值与标准偏差,并分析了孔隙率不确定性条件下Da数对Nu数的影响。结果表明,孔隙率不确定性对多孔介质方腔内自然对流有重要影响。随机多孔介质内流场及温度场与确定性条件下的流场及温度场存在一定偏差,Nu数标准偏差随着Da的增大先增大后减小。     

Effect of Local Thermal Non-equilibrium on the Onset of Convection in a Porous Medium Layer Saturated by a Nanofluid

The onset of convection in a horizontal layer of a porous medium saturated by a nanofluid is analytically studied. The model used for the nanofluid incorporates the effects of Brownian motion and thermophoresis. For the porous medium, the Darcy model is employed. The effect of local thermal non-equilibrium among the particle, fluid, and solid-matrix phases is investigated using a three-temperature model. The analysis reveals that in some circumstances the effect of LTNE can be significant, but for a typical dilute nanofluid (with large Lewis number and with small particle-to-fluid heat capacity ratio) the effect is small.     

Mixed Convection Boundary Layer Flow Past a Horizontal Circular Cylinder Embedded in a Bidisperse Porous Medium

Adopting a two-temperature and two-velocity model, appropriate to a bidisperse porous medium (BDPM) proposed by Nield and Kuznetsov (2008), the classical steady, mixed convection boundary layer flow about a horizontal, isothermal circular cylinder embedded in a porous medium has been theoretically studied in this article. It is shown that the boundary layer analysis leads to expressions for the flow and heat transfer characteristics in terms of an inter-phase momentum parameter, a thermal diffusivity ratio, a thermal conductivity ratio, a permeability ratio, a modified thermal capacity ratio, and a buoyancy or mixed convection parameter. The transformed partial differential equations governing the flow and heat transfer in the f-phase (the macro-pores) and the p-phase (the remainder of the structure) are solved numerically using a very efficient implicit finite-difference technique known as Keller-box method. A good agreement is observed between the present results and those known from the open literature in the special case of a traditional Darcy formulation (monodisperse system).     

Finite element analysis of convective heat transfer in porous media

The finite element method is used to analyse convective heat transfer in a porous medium. Convection past a vertical surface embedded in the medium and convection in a confined porous medium enclosure are analysed using the above method. The results are compared with those available in the literature and the agreement is found to be good. The method is applicable for two-dimensional analysis in a porous body of any arbitrary shape. The restriction of the boundary layer assumption is relaxed.     

The inertial effect on the natural convection flow within a fluid-saturated porous medium

The general momentum equation for fluid flow within a porous medium is supposedly valid for any fluid-porous medium configuration. One of the main concerns of using the general equations refers to the inclusion of both inertia terms, namely, the convective inertia term and the Forchheimer term. In this study, we go beyond the important discussion about the correctness of including both terms in the general momentum equations by focusing upon the effect of the convective inertia term on the heat transfer results. The fluid-porous medium system considered here is a cavity bounded by solid surfaces with vertical walls maintained at constant but different temperatures. The natural convection problem is solved numerically, and the results are compared with a general theory developed by using the method of scale analysis. It is demonstrated that the convective inertia term effect upon the heat transfer results is minor for 0.01 ≤ Pr ≤ 1, 10 ≤ Ra_D ≤ 10⁴, 10⁻⁸ ≤ Da ≤ 10⁻², and porosities 0.4 and 0.8. It is also shown that, contrary to the general belief, the convective inertial effect upon the heat transfer within the cavity is minimized when the Prandtl number is reduced.     

Non-darcy natural convection on a vertical cylinder in a saturated porous medium

The steady free convection boundary layer flow of non-Darcy fluid along an isothermal vertical cylinder embedded in a saturated porous medium using the Ergun model has been studied. The partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme developed by Keller. It is found that the heat transfer is strongly affected by the modified Grashof number which characterizes the non-Darcy fluid, and the curvature parameter. Also the heat transfer is found to be more than that of the flat plate.     

Effects of strong compressibility in natural convective flows through porous media with a near-critical fluid

The dynamics and heat transfer in a porous medium occupied by a liquid with parameters in the neighborhood of the critical point of “liquid-gas” transition are simulated numerically within the framework of the equations of dynamics of a porous medium with a compressible liquid phase and the Van-der-Waals equation of state. Adiabatic heating of the liquid phase in a porous layer initiated by a jump in temperature on one of the boundaries is investigated under microgravity conditions. Thermo-gravitational convection in the unsteady and steady-state regimes is simulated in rectangular domains and the effect of adiabatic heating on convection is studied. Calibration relations between the Rayleigh-Darcy and Prandtl numbers in the basic system of equations and their real analogs are obtained. A comparison is made with convection in a porous media occupied by a perfect gas.     

多孔介质中对流的周期性解与混沌

研究多孔介质内部有热源的对流传热．用高阶差分研究在不同的渗流瑞利（Ｒａｙｌｅｉｇｈ）数Ｒａ下对流随时间进展的演化情况（为比较起见也适当考察倾斜角的影响）．Ｒａ计算到大约１６０００．结果表明：在Ｒａ较小时对流是稳定的，Ｒａ增大到４６００出现了非稳定的但为周期性的解．随着Ｒａ进一步增大，出现一些混沌窗口．对于有侧斜角的情形，还出现阵发性     

Nonsimilar Boundary Layer Analysis of Free Convection Heat Transfer Over a Vertical Cylinder in Bidisperse Porous Media

This work presents a boundary layer analysis for the free convection heat transfer from a vertical cylinder in bidisperse porous media with constant wall temperature. A boundary layer analysis and the two-velocity two-temperature formulation are used to derive the nonsimilar governing equations. The transformed governing equations are solved by the cubic spline collocation method to yield computationally efficient numerical solutions. The effects of inter-phase heat transfer parameter, modified thermal conductivity ratio, and permeability ratio on the heat transfer and flow characteristics are studied. Results show that an increase in the modified thermal conductivity ratio and the permeability ratio can effectively enhance the free convection heat transfer of the vertical cylinder in a bidisperse porous medium. Moreover, the thermal nonequilibrium effects are strong for low values of the inter-phase heat transfer parameter.     

A two-phase,two-component model for natural convection in a porous medium

A numerical study of natural convection in a two-phase, two-component flow in a porous medium heated from below is presented. Interphase mass and energy transfer, latent heat and bouyancy effects are major physical features. This study extends earlier studies of natural convection based on single-phase, saturated porous medium models. The appearance of two-phase heat pipe zones in the flow has a marked effect on the fluid and heat flows as well as on the performance of the numerical methods. The numerical techniques for handling phase change, Jacobian construction and time step selection are discussed.     

Simultaneous heat and mass transfer under unsteady mixed convection along a vertical slender cylinder embedded in a porous medium

Unsteady laminar mixed convection flow (combined free and forced convection flow) along a vertical slender cylinder embedded in a porous medium under the combined buoyancy effect of thermal and species diffusion has been studied. The effect of the permeability of the medium as well as the magnetic field has been included in the analysis. The partial differential equations with three independent variables governing the flow have been solved numerically using a implicit finite difference scheme in combination with the quasilinearization technique. Computations have been carried out for accelerating, decelerating and oscillatory free stream velocity distributions. The effects of the permeability of the medium, buoyancy forces, transverse curvature and magnetic field on skin friction, heat transfer and mass transfer have been studied. It is found that the effect of free stream velocity distribution is more pronounced on the skin friction than on the heat and mass transfer. The permeability and magnetic parameters increase the skin friction, but reduce the heat and mass transfer. The skin friction, heat transfer and mass transfer are enhanced due to the buoyancy forces and curvature parameter. The heat transfer is strongly dependent on the viscous dissipation parameter and the Prandtl number, and the mass transfer on the Schmidt number.     

Triple-Diffusive Mixed Convection in a Porous Open Cavity

The triple-diffusive mixed convection heat and mass transfer of a mixture is analyzed in an enclosure filled with a Darcy porous medium. The mass transfer buoyancy effects due to concentration gradients of the dispersed components (pollutant components) are taken into account using the Boussinesq approximation model. The governing equations are transformed into a non-dimensional form, and six groups of non-dimensional parameters, including Darcy–Rayleigh number, Peclet number, two Lewis numbers for pollutant components 1 and 2 and two buoyancy ratio parameters for pollutant components 1 and 2, are introduced. The governing equations are numerically solved for various combinations of non-dimensional parameters using the finite element method. The effect of each group of non-dimensional parameters on the pollutant distribution and the heat transfer in the cavity is discussed. The results indicate that the presence of one pollutant component can significantly affect the pollutant distribution of the other component. When the Lewis number of a pollutant component is small, the increase in the bouncy ratio parameter of the proposed component always increases the Nusselt and Sherwood numbers in the cavity.