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.     

MHD natural convection flow with radiative heat transfer past an impulsively moving plate with ramped wall temperature

Influence of radiation on unsteady hydromagnetic natural convection transient flow near an impulsively moving vertical flat plate with ramped wall temperature in a porous medium is studied. Exact solution of momentum and energy equations, under Boussinesq approximation, is obtained in closed form by Laplace transform technique. The variations in fluid velocity and temperature are shown graphically whereas numerical values of skin friction and Nusselt number are presented in tabular form.     

Analytical and numerical investigation of double diffusion in thermally anisotropy multilayer porous medium

Double-diffusive natural convection within a multilayer anisotropic porous medium is studied numerically and analytically. The domain composed of two horizontal porous layers is subjected to a uniform horizontal heat flux and a vertical mass flux, where only the lower one is thermally anisotropic. Darcy model with classical Boussinesq approximation is used in formulating the mathematical model. The effect of thermal anisotropy and the relative width of the two layers on the flow and transfers is illustrated with characterising the transitions from the diffusive to the convective solution. Results were well compared with respect to a developed analytical approach, based on a parallel flow approximation for thermally anisotropic multilayer media.     

Vibration effect on convection onset in a system consisting of a horizontal pure liquid layer and a layer of liquid-saturated porous medium

The onset of convection in a system of two horizontal layers (a pure liquid and a porous medium saturated with the same liquid) heated from below under the action of vertical vibration is investigated. For describing the free thermal convection, in the liquid layer the Boussinesq approximation and in the porous layer the Darcy-Boussinesq approximation are used. In the limiting case of a thin liquid layer, effective boundary conditions on the upper boundary of the porous layer with account for convection in the liquid layer are obtained and it is shown that vibration has a stabilizing effect, whereas the presence of a liquid layer leads to destabilization. For an arbitrary liquid to porous layer thickness ratio the onset of convection is investigated numerically. In the case of a thin liquid layer there are two (short-and long-wave) unstable modes. In the case of thick layers the neutral curves are unimodal. Vibration has a stabilizing effect on perturbations with any wave number but affects short-wave perturbations much more strongly than long-wave ones.     

The Effect of Thermal Expansion on Porous Media Convection Part 1: Thermal Expansion Solution

The effect of thermal expansion on porous media convection is investigated by isolating first the solution of thermal expansion in the absence of convection which allows to evaluate the leading order effects that need to be included in the convection problem that is solved later. A relaxation of the Boussinesq approximation is applied and the relevant time scales for the formulated problem are identified from the equations as well as from the derived analytical solutions. Particular attention is paid to the problem of waves propagation in porous media and a significant conceptual difference between the isothermal compression problem in flows in porous media and its non-isothermal counterpart is established. The contrast between these two distinct problems, in terms of the different time scales involved, is evident from the results. While the thermal expansion is identified as a transient phenomenon, its impact on the post-transient solutions is found to be sensitive to the symmetry of the particular temperature initial conditions that are applied.     

Küppers-Lortz instability in rotating Rayleigh-Bénard convection in a porous medium

The Darcy-Lapwood-Brinkman model with the Boussinesq approximation is used to study Küppers-Lortz (KL) instability in the nonlinear regime of rotating Rayleigh-Bénard convection in a sparsely packed porous medium near the onset of stationary convection. The threshold Taylor numbers and critical angles for the onset of KL instability are obtained for different values of Λ, M for finite Prandtl numbers (1.5≤Pr≤100). Heat transfer is studied from Nusselt number at the onset of stationary convection.     

Investigation of convection in a horizontal cylindrical layer of a saturated porous medium

The structures of the convective motions and the nature of the heat transfer in a horizontal cylindrical layer are studied numerically for the Forchheimer model of a porous medium in the Boussinesq approximation. New asymmetric solutions of the equations of convection flow through a porous medium are found. Their development, domains of existence, and stability are investigated. One consists of a multivortex structure with asymmetric vortices in the near-polar region. Another asymmetric solution is realized at large Grashof numbers in the form of a convective plume deflected from the vertical. The threshold Grashof number of formation of the asymmetric motions depends on the Prandtl number and the cylindrical layer thickness.     

Influence of anisotropy of the permeability on the convection and heat transfer in a porous annular layer

Investigations into the convective transport of heat in porous materials are of interest for many applications in connection with the problem of increasing the efficiency of thermal insulation. In [1–5], convection in Isotropic porous media was considered. However, in many cases porous materials have an essential anisotropy of their permeability. Convective heat transfer has been inadequately studied for this case. In [6], the linearized equations were used to study the convection between infinite horizontal planes with a filling of an anisotropic material; the value of the critical Rayleigh number was found, and this agreed satisfactorily with experimental data. In the present paper, we investigate numerically convection between two infinite coaxial cylinders with an anisotropic porous filling, using the equations of convection in the Darcy—Boussinesq approximation [1–3]. The permeability tensor in the annular region is constructed from its principal values, which can be found experimentally. A method of calculation is developed and a parametric study made of the structure of the flow and of the local and averaged characteristics of the heat transfer, which are of interest for the design of thermal insulation.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 59–64, January–February, 1980.     

Influence of Property Variation on Natural Convection in a Gas Saturated Spherical Porous Annulus

The present study reports results of numerical investigation of natural convection in a spherical annulus filled with gas saturated porous medium. Boussinesq approximation has been relaxed considering variation of density, viscosity, and thermal conductivity of the gas with temperature. The momentum and energy equation along with boundary condition and property variation equations are solved by employing successive accelerated replacement scheme. There exists a reference parameter (q = 0.4) at which if thermo physical properties are evaluated, the average Nusselt number values for variable properties and constant properties are nearly same. Also, the effect of variable properties is negligible for temperature difference ratio, \(\theta ^{*}<0.1\) .     

Entropy Generation in Double-Diffusive Convection in a Square Porous Cavity using Darcy–Brinkman Formulation

The article reports a numerical study of entropy generation in double-diffusive convection through a square porous cavity saturated with a binary perfect gas mixture and submitted to horizontal thermal and concentration gradients. The analysis is performed using Darcy–Brinkman formulation with the Boussinesq approximation. The set of coupled equations of mass, momentum, energy and species conservation are solved using the control volume finite-element method. Effects of the Darcy number, the porosity and the thermal porous Rayleigh number on entropy generation are studied. It was found that entropy generation considerably depends on the Darcy number. Porosity induces the increase of entropy generation, especially at higher values of thermal porous Rayleigh number.     

Natural Convection Induced by a Centrifugal Force Field in a Horizontal Annular Porous Layer Saturated with a Binary Fluid

The Darcy Model with the Boussinesq approximation is used to study natural convection in a horizontal annular porous layer filled with a binary fluid, under the influence of a centrifugal force field. Neumann boundary conditions for temperature and concentration are applied on the inner and outer boundary of the enclosure. The governing parameters for the problem are the Rayleigh number, Ra, the Lewis number, Le, the buoyancy ratio, j{\varphi } , the radius ratio of the cavity, R, the normalized porosity, e{\varepsilon } , and parameter a defining double-diffusive convection (a = 0) or Soret induced convection (a = 1). For convection in a thin annular layer (R → 1), analytical solutions for the stream function, temperature and concentration fields are obtained using a concentric flow approximation and an integral form of the energy equation. The critical Rayleigh number for the onset of supercritical convection is predicted explicitly by the present model. Also, results are obtained from the analytical model for finite amplitude convection for which the flow and heat and mass transfer are presented in terms of the governing parameters of the problem. Numerical solutions of the full governing equations are obtained for a wide range of the governing parameters. A good agreement is observed between the analytical model and the numerical simulations.     

Effects of Viscous Dissipation on Unsteady Free Convection in a Fluid Past a Vertical Plate Immersed in a Porous Medium

The effects of viscous dissipation on unsteady free convection from an isothermal vertical flat plate in a fluid saturated porous medium are examined numerically. The Darcy–Brinkman–Forchheimer model is employed to describe the flow field. A new model of viscous dissipation is used for the Darcy–Brinkman–Forchheimer model of porous media. The simultaneous development of the momentum and thermal boundary layers are obtained by using a finite difference method. Boundary layer and Boussinesq approximation have been incorporated. Numerical calculations are carried out for various parameters entering into the problem. Velocity and temperature profiles as well as local friction factor and local Nusselt number are shown graphically. It is found that as time approaches infinity, the values of friction factor and heat transfer coefficient approach steady state.     

Anisotropy effect on the convection of a heat-conducting fluid in a porous medium and cosymmetry of the Darcy problem

The onset of convection in a porous rectangle is analyzed with account for the anisotropy of the thermal parameters and the permeability. For the Darcy–Boussinesq equations the conditions under which the problem pertains to the class of cosymmetric systems are established and explicit formulas for the critical Rayleigh numbers corresponding to the loss of stability of the mechanical equilibrium are derived. The critical numbers and the branching stationary convection regimes are calculated using a finite difference method conserving the problem cosymmetry.     

Numerical investigation of natural convective motions in a rotating cube

On the basis of a numerical integration of the Navier-Stokes equations in the Boussinesq approximation, the natural thermal convection in a rotating cube heated from below is investigated for three values of the Rayleigh number. The effect of the rotation velocity on the establishment of various forms of convection flows is studied. The existence of several types of convection flows is noted. Each of the motions is realized on a specific range of the rotation velocity of the cube. The domains of existence of the different types are determined. The results obtained are compared with experimental data.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 4, pp. 53–60, July–August, 1995.     

A coupled lattice BGK model for the Boussinesq equations

In this paper, a thermal lattice BGK model is developed for the Boussinesq incompressible fluids. The basic idea is to solve the velocity field and the temperature field using two independent lattice BGK equations, respectively, and then combine them into one coupled model for the whole system. The porous plate problem and the two‐dimensional natural convection flow in a square cavity with Pr=0.71 and various of Rayleigh numbers are simulated using the model. The numerical results are found to be in good agreement with the analytical solutions or those of previous studies. Copyright © 2002 John Wiley & Sons, Ltd.     

Numerical simulation of three-dimensional natural convection inside a heat generating anisotropic porous medium

Natural convection in anisotropic heat generating porous medium enclosed inside a rectangular cavity has been studied. A 3D finite volume based code is developed using the Darcy approximation and validated using experimental results of natural convection around an enclosed rod bundle. Subsequently, detailed simulation is carried out for a cavity, filled with orthotropic porous medium. The effects of heat generation, geometry and anisotropy are studied. Anisotropy is found to be of significant importance for both maximum value and distribution of temperature.     

Accuracy of the Boussinesq approximation for an incompressible fluid

The corrections of first order to the eigenvalues and critical Rayleigh numbers obtained in the Boussinesq approximation are determined for convection in a fluid with zero compressibility. The ratio of the equilibrium difference of the densities to a mean density of the fluid is taken as the small parameter. The corrections are found by the methods of perturbation theory for self-adjoint operators. It is shown that in the class of problems with symmetry with respect to a horizontal plane the first-order corrections vanish. The restrictions on the system needed if the Boussinesq approximation is to be meaningful in the problem of the occurrence of convective instability are established.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No, 2, pp. 19–26, March–April, 1981.     

An operator splitting numerical scheme for thermal/isothermal incompressible viscous flows

A numerical scheme for time‐dependent incompressible viscous fluid flow, thermally coupled under the Boussinesq approximation is presented. The scheme combines an operator splitting in the time discretization and linear finite elements in the space discretization, and is an extension of one previously applied for isothermal incompressible viscous flow governed by the Navier–Stokes equations. To show the efficiency of the scheme, numerical results are presented for mixed convection, and natural convection at high Rayleigh numbers. Restricting the scheme to the isothermal case, some numerical results at high Reynolds numbers are included, i.e. the scheme is tested for a small viscosity and a large force term, which are not trivial tasks to deal with. Copyright © 1999 John Wiley & Sons, Ltd.