Next Article in Journal
Cohesive Zone Model Modification Techniques According to the Mesh Size in Finite Element Models of Stiffened Panels with Debonding
Next Article in Special Issue
New Design and Optimization of a Jet Pump to Boost Heavy Oil Production
Previous Article in Journal / Special Issue
On Numerical Modeling of Thermal Performance Enhancementof a Heat Thermal Energy Storage System Using a Phase Change Material and a Porous Foam
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Thermal Performances in a Ventilated Room Using LBM-MRT: Effect of a Porous Separation

1
Laboratory of Multiphase Transport and Porous Media (LTPMP), Faculty of Mechanical and Proceeding Engineering, University of Sciences and Technology Houari Boumediene, Algiers 16111, Algeria
2
Laboratory of Energy and Mechanical Engineering (LEMI), Faculty of Technology, University of M’Hamed Bougara Boumerdes, Boumerdes 35000, Algeria
*
Author to whom correspondence should be addressed.
Computation 2022, 10(1), 4; https://doi.org/10.3390/computation10010004
Submission received: 30 November 2021 / Revised: 29 December 2021 / Accepted: 4 January 2022 / Published: 10 January 2022
(This article belongs to the Special Issue Computational Heat, Mass, and Momentum Transfer—III)

Abstract

:
This article demonstrates the feasibility of porous separation on the performance of displacement ventilation in a rectangular enclosure. A jet of fresh air enters the cavity through an opening at the bottom of the left wall and exits through an opening at the top of the right wall. The porous separation is placed in the center of the cavity and its height varies between 0.2 and 0.8 with three values of thickness, 0.1, 0.2, and 0.3. The heat transfer rate was calculated for different intervals of Darcy (10−6 ≤ Da ≤ 10), Rayleigh (10 ≤ Ra ≤ 106), and Reynolds (50 ≤ Re ≤ 500) numbers. The momentum and the energy equations were solved by the lattice Boltzmann method with multiple relaxation times (LB-MRT). Schemes D2Q9 and D2Q5 were chosen for the velocity and temperature fields, respectively. For porous separation, the generalized Darcy–Brinkman–Forchheimer model was adopted. It is represented by a term added in the standard LB equations. For the dynamic domain, numerical simulations revealed complex flow structures depending on all control parameters. The results showed that the thermal field, mainly in the second compartment, is very dependent on the size and permeability of the porous separation. However, they have no influence on the transfer rate.

1. Introduction

Heat transfer in enclosures has been of great interest to researchers due to the many applications arising from such geometries. Among the numerous works and papers published in the field of thermal flow in closed cavities, we cite the work of de Vahl Davies [1]. The latter considered as the reference work concerning heat transfer in a closed enclosure. The author considered a driving thermal gradient on the vertical walls, while the horizontal walls were considered adiabatic. The results have been validated several times and are considered as a benchmark for the case of closed cavities. Thereafter, Oztop [2] numerically analyzed the convective heat transfer in a porous rectangular enclosure inclined by a variable angle α. One of the side walls is kept at a constant temperature and another one is partially cooled, while the other walls are adiabatic. After solving by the finite volume method, he found that the tilt angle and the aspect ratio dominate the heat transfer evolution. In the same year, Sathiyamoorthy et al. [3] analyzed the natural convection in closed cavity using the finite volume method. The bottom wall is uniformly heated, but the vertical walls are linearly heated while the top wall is adiabatic. The results were presented for a Rayleigh number (103 < Ra < 105) and a Prandtl number (0.7 ≤ Pr ≤ 10). Varol et al. [4] studied the steady convection in a rectangular enclosure filled with a porous medium. The temperature profile in the base wall varies according to sinusoidal function, while the other walls are isolated. The problem was solved using a finite volume method, and the controls parameters are: Rayleigh number (10 ≤ Ra ≤ 1000), aspect ratio (0.25 ≤ AR ≤ 1.0), and amplitude of the temperature (0.25 ≤ λ ≤ 1.0). As the main result, they found that the increase in amplitude λ causes the increase in convective heat transfer, whereas decreases with the diminution of the aspect ratio. Ameziani et al. [5] added the effect of a moving wall on the thermal exchange by considering the cooperating and opposing forces of natural and forced convection. The authors found that a heat transfer minimum designated by the Nusselt number is observed for a specific pair of (Ra-Re). Huelsz and Rechtman [6] treated the heat transfer by the natural convection of air in a tilted, differentially heated square cavity. The simulation was carried out by the lattice Boltzmann method for a laminar regime with variation in the inclination angle and the Rayleigh number. The results show that for a given value of the Rayleigh number, a hysteresis is observed, whatever the value of the tilt angle.
In the case of ventilated cavities, the number of studies is rather abundant; in the last few years, we can cite the work of Omri and Nasrallah [7] on the transient mixed convection in the laminar regime. The cavity is cooled by a fresh fluid at a temperature lower than the initial temperature of the cavity. The governing equation was solved with a finite volume method. The dynamic and thermal fields are numerically resolute for a wide range of Re and Ri, and for different fluid inlet–outlet positions. In configuration (A), the results show that a single vortex turning in a clockwise direction is generated when Re increases with Ri < 1, or with an increase in Ri at moderate values of Re. However, the flow becomes multi-cellular when inertial and buoyancy effects are significant. In configuration B, increasing Ri reverses the fluid clockwise, but increasing Re reverses the fluid counterclockwise, and the dynamic field has several structures, depending on the values of the control parameters. Najam et al. [8] numerically studied the mixed convection in a T-shaped cavity heated with a constant heat flow and cooled with an air jet from the bottom. The heater blocks are identical, and the system is symmetrical about a vertical axis through the centers of the openings. A finite difference method was used to solve the governing equations. The results obtained for H/L = 1 indicate the existence of multiple solutions, and the heat transfer depends substantially on the control parameter values. Later, Gan [9] used CFD for the simulation of buoyancy-driven natural ventilation in a vertical enclosure for different heat flow intensities and wall thermal distributions using two computational domains. The simulation results show that the ventilation rate and the heat transfer coefficient depend on the inlet position as well as on the cavity size. For the case of ventilated cavities with porous media, studies are rather rare; we quote the work of Tong and Subramaniant [10], who studied natural considered convection in two rectangular cavities filled with porous media, with the goal to establish the heat transfer characteristics. The flow in the porous region was modelled by a modified revised Darcy’s law with the no-slip condition. Through the finite element solution, the authors concluded that there are regions wherein the heat transfer can be reduced by partially filling the cavity with the porous medium. A few years later, Moraga et al. [11] used the finite volume method (MVF) to study mixed convection in a ventilated rectangular cavity, vertically divided into two different porous media; they studied the effect of several control parameters such as the Darcy number, Reynolds number, Richardson number, aspect ratio, and flow direction. They found that the friction coefficient of the walls is larger when Re = 500 and Ri = 10 in both flow directions. Additionally, when the Darcy number increases, the velocity gradients increase near the walls, resulting in an increase in the friction coefficient (0 < Cf Re < 47). Thereafter, Mehrizi et al. [12] carried out a numerical analysis of heat transfer by forced convection in a ventilated cavity with an obstacle heated at a constant temperature. The walls of this cavity and the inlets/outlets are adiabatic. They used LBM for incompressible flow in a porous medium and the Darcy–Brinkman–Forchheimer model to consider the effect of the Reynolds number and Prandtl on the heat flow and temperature distribution. They found that increasing the Reynolds and Prandtl numbers enhanced the heat transfer characterized by the average Nusselt number. Oztop et al. [13] considered the natural convection heat transfer in a partially open cavity filled with porous media. They used the Darcy–Brinkman–Forcheimer model in the case where the left vertical wall had a constant temperature. They examined the different control parameters (porosity, Darcy number, Grashof number, and length of the heated wall). The authors found that the Nusselt number increases when the Grashof number increases, due to strengthening buoyancy-driven flows. In partially porous cavities also, Liu and He [14] are interested LBM-MRT modelling of incompressible flow. Porous medium is modeled by the Brinkman and Forchheimer models as extensions of the Darcy model. Additionally, through the Chapman–Enskog analysis, they were able to establish the generalized Navier–Stokes equations and found that their numerical results agreed with those findings with the analytical method. Shuja et al. [15] numerically studied heat transfer in a ventilated cavity with two porous blocks. The authors considered a moderate Reynolds number (Re = 100) over a large porosity range for two aspect ratios. The obtained results showed that the aspect ratio has little influence on the transfers which are directly dependent on the porosity. In second numerical research with the same configuration, the same authors [16] studied the flow field and heat transfer (Nusselt number). They showed that, for the first block (inlet fresh air side), the increase in Gr enhances the heat transfer and favors the natural convection phenomenon. Hireche et al. [17] performed a numerical analysis to determine the thermal convection characteristics in a ventilated rectangular cavity with an aspect ratio of (A = 2) and a porous partition in the middle. They concluded that heat transfer increases with an increase in Re and Ra values. Additionally, they concluded that the effect of the Darcy number is very significant for small permeability, while the effect of the porous layer height is hardly significant with a maximum deviation of less than 7%. Recently, several papers have been published on ventilated cavities, with different geometries and opening positions, even changing the physical properties of fluid [18,19,20,21,22,23,24,25,26].
From this synthesis, it can be seen that much of the numerical work has been carried out on ventilated cavities. For the case of displacement ventilation, few works have treated heat transfer in a rectangular room with a porous separation. However, most of the previous studies do not take into account the type of separation (solid or porous) and its dimensions (thickness and height), parameters on which this work is based—in particular, the separation thickness. In this paper, we explore the possibility of modeling the heat removal process by displacement ventilation in professional premises, large rooms, open plan offices, or general hospital divisions, providing air exchange and being among the most efficient systems. It presents an interesting solution for enhancing air quality, energy saving, and thermal comfort in buildings.

2. Mathematical Formulation

The present study focuses on convection in a rectangular open cavity of height H and length L = 2H, as shown in Figure 1, partitioned by a vertical porous wall of thickness “e” varying between H/10 and 3H/10 as well as a height that varies between H/5 and 4H/5. The porous wall is saturated by a single fluid. The air in the cavity is renewed through two diametrically opposed openings located at the lower and upper corners of the left and right walls, respectively. The walls of the cavity are impermeable, so the left wall is heated at a constant temperature and the other walls are adiabatic.
Some assumptions have been considered to simplify the mathematical formulation of our model. The fluid is Newtonian and incompressible. The fluid flow into the cavity is laminar and two-dimensional, while the Boussinesq approximation has been applied to model the density variations. The porous bloc is assumed to be homogeneous, isotropic, and in thermodynamic equilibrium with the fluid. The thermo-physical properties of fluid are considered constant over the temperature scale. It is also assumed that radiative transfer is neglected.

2.1. Governing Equations

2.1.1. Fluid Zone

The equation of continuity expresses the theory of mass conservation of a moving fluid. In an elementary volume, the mass of the fluid removed or injected is equal to the sum of the variation fluid density in an elementary time interval, and the mass flows through the surface of this volume. The transport equations express the theory of momentum and energy conservation. They are defined as:
  • Continuity equation:
u = 0
  • Equation of momentum conservation:
    ρ f [ u t + u ( u ) ] 1 = P 2 + μ 2 u 3 + ρ f G 4 ,
    where:
    (1)
    Macroscopic inertial forces;
    (2)
    Pressure gradient;
    (3)
    Diffusion term;
    (4)
    Volume force, where the density is expressed according to the Boussinesq approximation.
  • Equation of energy conservation:
ρ f c p [ T t + u ( T ) ] = ( k f T )
where:
cp: Specific heat.
kf: Fluid conductivity.

2.1.2. Porous Zone

Based on the generalized Darcy–Brinkman–Forchheimer model, the dimensional governing equations for convection heat transfer in a homogeneous, isotropic, and fluid-saturated porous medium can be written as [27,28,29]:
  • Continuity equation:
ε ρ f t + ( ρ f u ) = 0
  • Equation of momentum conservation:
    ρ f [ u t + u ( u ) ] 1 = P 2 + μ e f f 2 u 3 υ K u 4 1.75 ε 150 ε K | u | u 5 + ρ f G 6
    (1)
    Macroscopic inertial forces;
    (2)
    Pressure gradient;
    (3)
    Brinkman’s viscous term;
    (4)
    Darcy term;
    (5)
    Forchheimer term (microscopic inertial forces [30]);
    (6)
    Volume forces.
  • Equation of energy conservation:
( ρ 0 c p ) m T t + ( ρ 0 c p ) f ( T u ) = ( k e f f T )
where ( ρ 0 c p ) m and ( ρ 0 c p ) f are, respectively, the thermal heat capacity of the porous medium and the fluid, respectively, and keff is the effective thermal conductivity.

2.2. Dimensionless Equations

Governing equations are written in dimensionless form using the reduced variables:
X = x H , Y = y H , U = u U 0 , V = v U 0 , P = p ρ U 0 2 , θ = T T c T , T = T h T c , t = t H / U 0
The governing equations system is given in dimensionless form in the next section.

2.2.1. Fluid Zone

  • Continuity equation:
U X + V Y = 0
  • Equation of momentum conservation following X
U t + [ U U X + V U Y ] = P X + 1 Re ( 2 U X 2 + 2 U Y 2 )
  • Equation of momentum conservation following Y :
V t + [ U V X + V V Y ] = P Y + 1 Re ( 2 V X 2 + 2 V Y 2 ) + R a Pr Re θ
  • Equation of energy conservation:
θ t + ( U θ ) X + ( V θ ) Y = 1 Re Pr ( 2 θ X 2 + 2 θ Y 2 )

2.2.2. Porous Zone

  • Continuity equation:
U X + V Y = 0
  • Equation of momentum conservation following X [30]:
U t + [ U ε U X + V ε U Y ] = ε P X 1 D a Re U 1.75 150 ( U 2 + V 2 ) 1 / 2 D a V + 1 Re ( 2 U X 2 + 2 U Y 2 )
  • Equation of momentum conservation following Y [30]:
V t + [ U ε V X + V ε V Y ] = ε P Y 1 D a Re V 1.75 150 ( U 2 + V 2 ) 1 / 2 D a V + 1 Re ( 2 V X 2 + 2 V Y 2 ) + R a Pr Re θ  
  • Equation of energy conservation:
γ θ t + ( U θ ) X + ( V θ ) Y = R k Re Pr ( 2 θ X 2 + 2 θ Y 2 )
where γ and Rk are, respectively, the heat capacity ratio and the thermal conductivity ratio. In our work, they have been assumed unit values.
The non-dimensional form of the latter equations gives the appearance of the following non-dimensional parameters mentioned in Table 1.
The hydrodynamic boundary conditions, characterized by the non-slip condition of the fluid particles on the solid walls and their impermeability, are given:
Left   vertical   wall :   U ( 0 , Y ) = V ( 0 , Y ) = 0 ; 0.1 Y 1
Right   vertical   wall :   U ( 2 , Y ) = V ( 2 , Y ) = 0 ; 0 Y 0.9
Horizontal   surfaces :   U ( X , 0 ) = V ( X , 0 ) = U ( X , 1 ) = V ( X , 1 ) = 0 ; 0 X 2
The   inlet :   { X = 0 0 Y 0.1 U = 1 , V = 0 .
The   outlet :   { X = 2 0.9 Y 1 U X = 0 , V = 0 .
Thermal boundary conditions:
Left   vertical   wall :   X = 0 ; 0.1 Y 1 , θ = 1
The   inlet :   θ = 0 , 0 Y 0.1
Right   vertical   wall :   X = 2 , 0 Y 0.9 , θ X | X = 2 = 0
The   outlet :   θ X = 0 , 0.9 Y 1
Horizontal   surfaces :   θ Y | Y = 0 = 0 , θ Y | Y = 1 = 0 ; 0 X 2
Finally, we express the heat transfer rate on the active surface by the averaged and local Nusselt numbers, defined, respectively, as:
N u = a H 1 N u X Y  
with:
N u X = θ X | X = 0

3. Lattice Boltzmann Method

The Boltzmann method on lattice, commonly called the LBM (lattice Boltzmann method), is another approach for numerical simulation of physical phenomena. It is a relatively new method that differs from traditional methods. The construction of the Boltzmann method on a lattice is represented by two independent steps: the development of statistical physics and the appearance of cellular automata [31,32]. In 1872, Boltzmann proposed his famous equation describing the spatio-temporal evolution of a function f representing the distribution of particles with a given velocity at a given place and time. This function is often named the distribution function and thus depends on time, space, and velocity:
f t + c i f x i + F i m f c i = ( f t ) c o l l
The right side of the Boltzmann equation, called the collision operator, represents the influence of collisions between particles. If this term is equal to zero, the particles are simply affected and subjected to the action of the force F present on the left-hand side. If the particles collide, the evolution of the system depends on the form of the collision operator.
The work of Chapman [33] and Enskog [34], validated experimentally by the chemist FW Dootson, led to direct relations between the Boltzmann equation and those of Navier–Stokes (1823), thus extending the work of Boltzmann. However, the collision operator chosen by Chapman and Enskog remained very complex and typic. It was then necessary to wait many decades until the development of a simple collision model; it was the mathematicians Bhatnagar, Gross, and Krook who gave birth to it in 1954 [35]. This model, established from complex theoretical developments, is based on the idea that a particle collision can be described as the relaxation in a given time of the particles, towards a given equilibrium state. This fundamental idea makes it possible to write the collision operator in a very simple form that is named the BGK operator, from the name of its creators. Statistical physics showed, as early as 1954, that the Boltzmann equation provided with the BGK operator could be used to describe fluid mechanics flows governed by the Navier–Stokes equations. The Boltzmann method on lattices could have been born at that time, but the notion of discretization and networking, specific to computer use, only emerged much later with the intensive development of numerical simulation.
In what follows, we will not consider any external force (Fi), and the effects of gravity will be neglected. We notice that the density, momentum, and kinetic energy of the fluid are found by considering all possible speeds:
ρ = ƒ d c
ρ u = c ƒ d c
ρ e = 1 2 | c u | 2 m ƒ d c
Considering only binary collisions transforms the velocities [c1; c2] (pre-collision) into [c1′; c2′] (post-collision). In the absence of a collision, the particles are in thermodynamic equilibrium. Such a state is described by an equilibrium distribution function feq, a solution to the equation:
( f t ) c o l l = Ω
Such a function therefore verifies:
ƒ i e q = w i ρ [ 1 + 3 ( c i u ) + 9 2 ( c i u ) 2 3 2 u 2 ]  
where:
wi: network constants;
C: speed of sound.

3.1. Discretization of the Lattice Boltzmann Equation

The lattice Boltzmann method can be applied for the simulation of fluid flow, heat, and mass transfer. The fluid is considered as a set of particles, moving at a discrete velocity in specified directions, according to the lattice structure, and colliding (interact) with each other at the nodes of the lattice. The distribution function is the probability of the presence of a particle at a certain node with a certain velocity. The collision and propagation of particles is governed by the Boltzmann equation of the lattice, which reflects the evolution of these functions. The external force implemented by a source term:
f i ( x + e i δ t , t + δ t ) f i ( x , t ) = Ω ( f i ) + δ t F i
where f i is the distribution function with the velocity e i at the lattice node x at time t, δ t is the discrete time step, Ω ( f i ) is the discrete collision operator, and F i is the external forces term.

Collision Operator

There are several approaches to calculating the collision operator, such as:
(a)
The BGK approximation:
Bhatnagar, Gross, and Krook (BGK) developed a simplified model of the collision operator in 1954 [35]. The collision operator can be written as follows:
Ω = ω ( ƒ e q ƒ ) = 1 τ ( ƒ e q ƒ )
ω = 1 τ
τ: relaxation factor;
ω: collision frequency.
The local distribution function is denoted by feq, which is the Maxwell–Boltzmann distribution function.
The Boltzmann equation (without external forces) can be written after the application of the BGK approximation, as follows:
f t + c f = 1 τ ( ƒ e q ƒ )
The discrete velocity of Boltzmann equation can be written as follows:
f i t + c i f i = 1 τ ( ƒ i e q ƒ i )
The discretization of the previous equation gives us:
f i ( x + c i Δ t , t + Δ t ) = f i ( x , t ) + Δ t τ [ f i e q ( x , t ) f i ( x , t ) ]
(b)
Multiple Relaxation time MRT:
In order to overcome these deficiencies of the LBGK model (numerical instability, applicable only for Pr = 1), d’Humières [36] proposes the approach of moment, known as the multiple relaxation time model. To simulate the evolution of macroscopic quantities, different relaxation times are used. By running the collision in an orthogonal space–time, an efficient and flexible scheme is obtained. This model has optimal stability compared to the BGK model. In the MRT scheme, the collision term is written as [37]:
Ω i = M 1 S i [ m i ( x , t ) m i e q ( x , t ) ]
Therefore, the Boltzmann equation is written as:
f i ( x + c i Δ t , t + Δ t ) f i ( x , t ) = M 1 S i [ m i ( x , t ) m i e q ( x , t ) ]
The velocity is equal to the ratio between displacement dx and the unit time dt. m ( x , t ) = e t   m i e q ( x , t ) are moment vectors; m = ( m 0 , m 1 , m 2 ,   , m n ) T . The relaxation matrix is a diagonal matrix, and the linkage between velocities and moment spaces is given from linear transformation as follows:
m = M f
f = M 1 m

3.2. Dynamic Lattice Boltzmann Model

The terminology used in the LBM method to describe the dimension and number of model directions is as follows: DnQm, where: n and m designate the problem dimension (1D, 2D, and 3D) and the number of the function directions (discrete velocities).

3.2.1. The Two-Dimensional Nine-Speed Model (D2Q9)

This is the most widely used model, especially for the study of flow problems. The D2Q9 model is based on a square lattice of step δ x = δ y = 1 , and each lattice pattern is characterized by nine discrete velocities ci(i = 0, …, 8) [38].
For the hydrodynamic field, the two-dimensional nine-velocity D2Q9 model is used in this work as well as the Boltzmann equation of the multi-relaxation time lattice (MRT-LBM) [39,40]:
f i ( x + e i δ t , t + δ t ) f i ( x , t ) = M 1 S i [ m i ( x , t ) m i e q ( x , t ) ]  
The nine discrete velocities are given by:
e i = { ( 0 , 0 ) i = 0 [ cos ϕ , sin ϕ ] c ϕ = ( i 1 ) π 2 i = 1 , 2 , 3 , 4 [ cos ϕ , sin ϕ ] 2 c ϕ = ( 2 i 9 ) π 4     i = 5 , 6 , 7 , 8
where c = δ x / δ t = δ y / δ t is the lattice velocity, in which δ x / δ t and δ y / δ t are the lattice cell widths; these quantities are chosen equal to the unit δ x = δ y = δ t = 1 .
The distribution function of local equilibrium f i e q for the hydrodynamic field in porous media is calculated by equation:
f i e q = w i ρ ( 1 + 3 e i v + 9 ( e i v ) 2 2 ε + 3 v 2 2 ε )  
where the weighting factors wi are given in the literature as follows: w 0 = 4 9 ; w 1 4 = 1 9 ; w 5 8 = 1 36 .
M as the transformation matrix, which projects fi and fieq into space–time with m = MF and meq = MFeq, where m and meq are vectors of macroscopic and macroscopic equilibrium variables, respectively. The distribution functions in space–time are given by:
m = ( ρ e φ j x δ t 2 ρ F x q x j y δ t 2 ρ F y q y p x x p x y ) = [ 1 1 1 1 1 1 1 1 1 4 1 1 1 1 2 2 2 2 4 2 2 2 2 1 1 1 1 0 1 0 1 0 1 1 1 1 0 2 0 2 0 1 1 1 1 0 0 1 0 1 1 1 1 1 0 0 2 0 2 1 1 1 1 0 1 1 1 1 0 0 0 0 0 0 0 0 0 1 1 1 1 ] ( f 0 f 1 f 2 f 3 f 4 f 5 f 6 f 7 f 8 ) = M f
where ρ is the fluid density, e is related to the energy, φ is related to the square of the energy, jx,y are the kinetic momentum components j = (jx,jy) = (ρvx,ρvy), qx,y correspond to the energy flow in two directions, and pxx,xy correspond to the diagonal and off-diagonal components of the stress rate tensor.
The density ρ and the momentum jx,y are conserved quantities, while the other six velocity moments are not conserved. The equilibrium moments for the non-conserved moments meq can be defined as [41]:
e e q = 2 ρ + 3 ρ 0 | v | 2 ε ; ϕ e q = ρ 3 ρ 0 | v | 2 ε ; q x e q = ρ 0 v x ; q y e q = ρ 0 v y ; p x x e q = ρ 0 ( v x 2 v y 2 ) ε ; p x y e q = ρ 0 v x v y ε
In the above equilibrium moments, the incompressibility approximation has been used, the fluid density is written as ρ = ρ0 + δρδ0 (δρ is the density fluctuation), and (jx,jyρ0v). The average density of the fluid ρ0 is taken to be 1 to simplify the model. I is the identity matrix and S is a diagonal relaxation matrix given by:
S = d i a g ( s 0 , s 1 , s 2 , s 3 , s 4 , s 5 , s 6 , s 7 , s 8 )
Si are the relaxation rates. These values must be between 0 and 2 to keep the model stable. In our simulations, they are considered as follows:
S = d i a g ( 1 , 1.2 , 1.4 , 1 , 1.2 , 1 , 1.2 , 1 τ , 1 τ )  
In the LBM-MRT method, the fluid kinematic viscosity is related to the relaxation parameter of the flow field by the following relation:
υ = 1 3 ( τ 1 2 )
The macroscopic fluid variables, density and velocity, are calculated as follows:
Mass conservation can be described as the summation of the distribution functions in each the lattice node, which represent the macroscopic density of the fluid [42]:
ρ = i = 0 8 f i
The conservation of momentum can be described as the summation of the product of the microscopic lattice velocities by the distribution function:
ρ u = i = 0 8 f i e i
Thus:
u = 1 ρ i = 0 8 f i e i

3.2.2. D2Q5-MRT-LB Equation for the Thermal Field

For the temperature field, the two-dimensional model D2Q5 with five discrete velocities is used in our work.
The LBM-MRT equation can be written as follows:
g i ( x + e i δ t , t + δ t ) g i ( x , t ) = N 1 E i [ n i ( x , t ) n i e q ( x , t ) ]  
where gi(x,t) is the temperature distribution function. In the D2Q5 model, the five discrete speeds are given by the following relation:
e i = { ( 0 , 0 ) i = 0 [ cos ϕ , sin ϕ ] c ϕ = ( i 1 ) π 2 i = 1 , 2 , 3 , 4
The equilibrium distribution function is given by:
g i e q = w i T [ 1 + 5 e i v ]  
Here, the weighting factors wi are given as follows; w 0 = 3 / 5 ;   w 1 4 = 1 / 10 .
N is an orthogonal transformation matrix 5 × 5, which projects gi and gieq into space–time with n = Ng and neq = Ngeq, where n and neq are, respectively, the variables of macroscopic and macroscopic equilibrium temperature.
For temperature, the distribution functions in space–time are given by [43]:
n = ( n 0 n 1 n 2 n 3 n 4 ) = | 1 1 1 1 1 0 1 0 1 0 0 0 1 0 1 4 1 1 1 1 0 1 1 1 1 | ( g 0 g 1 g 2 g 3 g 4 ) = N g  
E is the diagonal relaxation matrix:
E = d i a g ( r 0 , r 1 , r 2 , r 3 , r 4 )
To ensure good stability of computation, the relaxation rates are chosen in our work as follows:
E = d i a g ( 1 , 1 / τ T , 1 / τ T , 1.25 , 1.25 )
where τ T M is the relaxation parameter of temperature field.
The collision process D2Q5 LBM-MRT proceeds as follows:
n i * ( x , t ) = n i ( x , t ) E [ n i ( x , t ) n i e q ( x , t ) ]
Between the five moments, only the temperature is kept:
n 0 = T = i = 0 8 g i  
In equilibrium, the non-conserved moments n e q = N g e q are written as follows [44]:
n 1 e q = v x T ; n 2 e q = v y T ; n 3 e q = ϖ T ; n 4 e q = 0

3.3. Boundary Condition

The hydrodynamic boundary conditions that govern our problem are mainly the adherence of the fluid to all walls as well as a horizontal imposed inlet velocity and a velocity establishment on the outlet opening. For the thermal boundary conditions, we have an active left wall (hot temperature) with an insulation of all other walls. We also consider a fresh air source (cold temperature) at the inlet with a thermal establishment at the outlet. One of the most important parts of a Boltzmann lattice simulation is the establishment of boundary conditions. The establishment of boundary conditions for classical computational methods used in hydrodynamics (CFD) is carried out by obtaining the flow variables. This is not the case for the LBM method, where the distributions functions entered in the domain must be determined. Therefore, the appropriate equations for calculating of distribution functions at the domain boundary for each boundary condition.
In our work, we used the rebound boundary condition [45,46,47,48] to model the adhesion to all walls. We also used the boundary conditions for the imposed velocity type (Zu and He [42]) to model the imposed velocities and temperatures.

4. Validation of the Calculation Code

To validate our computational code, we compared our results with previous work (those of Rahman et al. [49]). The latter have studied mixed convection in a ventilated square cavity with a horizontal heat square block; the Richardson number for this study varies from 0.0 to 5.0, the Reynolds number has been fixed at Re = 100, and the fluid in the cavity is the air with Pr = 0.71. The comparison of results (Nusselt number at the hot wall) obtained from our computational code and those of the author considered was performed for the zero-dimensional limiting case of the obstacle, for different Richardson numbers (Figure 2). The confrontation shows a good agreement between our numerical simulations and those of the previous author.

5. Results and Interpretations

The required number of lattices is obtained through a grid sensitivity analysis to ensure accurate results. Due to the multiple parameters controlling the flow pattern, temperature field, and heat transfer, all calculations were performed by setting a constant Prandtl number (Pr = 0.71; the air) and a finite aspect ratio A = 2. The Reynolds number varies up to 500 and the Rayleigh number up to 106. The ratio of conductivities and the ratio of viscosities between the fluid and porous media are equal to unity.
Figure 3 shows the stream functions (tangent to velocity vector) for different Darcy numbers, Rayleigh numbers, and for three Reynolds values. Note that the negative values of the stream functions are shown as dotted lines and that the porous medium thickness is taken to be constant at Ep = 0.2. In order to cover all the ranges considered, three values of Darcy are taken (10−2; fluid regime 10−4 low-permeability medium and Da = 10−6 for the Darcian regime). The Rayleigh number is assumed in the range of 102–106 in order to sweep from the diffusive to the convective regime. When Da = 10−2 (fluid regime), the porous medium does not offer resistance to the fluid. As a result, the stream lines are continuous, and no distortion is observed. In this situation, the fluid enters through the inlet (bottom left) as a jet which blossoms and takes the form of a fan; the fluid is then sucked through the outlet. When the Rayleigh number increases, we notice the formation of a natural convection cell then a second cell (for Ra = 106) in the vicinity of the parietal heating. These natural convection cells disappear for large Reynolds values. It should be noted that for low Reynolds values, the increase in thermal draft leads to a recirculation cell (dotted line because of its negative values) in the vicinity of the bottom right corner. The increase in Reynolds number drives the fluid (like a jet) towards the right wall and it is aspirated to the outlet in the vicinity of the second compartment of the cavity (Re = 250). For a higher Reynolds number (Re = 500), the jet projected on the left wall comes back to the right; it is then aspirated by the ascending flow all along the two upper walls (hot wall and the one on the top), creating a vacuum at the bottom of the cavity, quickly filled by a recirculation cell. When the Darcy number decreases (i.e., Da = 10−4 and Da = 10−6), the porous wall behaves like a solid and the moving fluid bypasses the porous medium following the path of least resistance. It should also be mentioned that when the Darcy number is very low, the fluid is directly projected towards the outlet, creating a stagnant zone downstream of the porous wall, filled by a negative recirculation. In this situation, the thermal comfort is affected in this second compartment totally or partially filled by this recirculation.
For high Reynolds values (Re = 500), the main observation is that the fluid follows a back-and-forth path along the cavity and is then aspirated towards the hot wall; it is then deviated along the two upper walls to the exit orifice. By following the path of the top wall, the flow becomes more complex with a recirculation on the whole upper part (between the fluid jet and the porous medium). Note that for lower permeabilities (Da = 10−6), the effect of the inlet jet is concentrated on the first compartment with the appearance of the recirculation cell between the top wall and the porous medium, which leads to another cell of positive values in the second compartment. For reasons of illustration space, Re = 50 has been omitted, considered as predictable in Figure 3. The fluid regime, Da = 10−2, has been omitted since it considers the fluid medium extensively studied in the literature.
In Figure 4, the stream functions are represented for two thicknesses of porous medium (Ep = 0.1 and 0.3) with a scan of two Reynolds numbers (Re = 250 and Re = 500). In this figure also, we find the illustrations of the figures related to two Darcy numbers (Da = 10−4 and 10−6) as well as three Rayleigh numbers, which give us the sweep of the range Ra [102–106]. We observe in this figure that the thickness of the porous wall preserves the same principle of general flow movement, mainly with three combined movements. With a jet of fresh fluid coming from the inlet port, the last bypasses the porous medium, sometimes to the outlet and other times to the active wall for large Reynolds values. Note that the porous wall thickness has very little influence on the flow, and its effect is only visible for intermediate permeabilities (Da = 10−4) wherein the dimension of the recirculation zone of the second compartment is affected.
The influence of the porous medium height is illustrated in Figure 5 for the two heights considered (Hp = 0.2 and Hp = 0.8). This figure confirms the relative independence of the flow structures with the porous wall thickness in the considered thickness range Ep ∈ [0.1–0.3]. Note that when the height of the porous wall is important, the second compartment is characterized by at least one recirculation zone for low Rayleigh values. When the thermal draft increases (Ra = 106), a recirculation zone is also obtained in the first compartment of the cavity. When the porous wall is one fifth of the chamber height (Hp = 0.2), the stream functions follow the same path as the fluid regime (Figure 3 for Da = 10−2), the only difference being the size of the swirls.
Figure 6, Figure 7 and Figure 8 illustrate isotherms as a function of the control parameters for the same conditions chosen in Figure 3, Figure 4 and Figure 5, respectively. It is noted that the interval between two successive isotherms is constant and is equal to 0.05. It is also mentioned that the dimensionless temperature varies from the minimum value at the inlet (θ = 0) to the maximum value near the active wall (θ = 1). It can be seen from these figures that the porous medium does not offer any resistance to the isotherms, and no distortion is observed. This is probably due to the fact that a unit conductivity ratio has been considered. The most pertinent remark is that the influence of the porous wall thickness is not very visible on the isotherms (confirmed in Figure 8), which may indicate that the Nusselt numbers (on the heat transfer figures) will be independent of the porous wall thickness.
In Figure 6, we notice that for large values of Darcy number (i.e., Da = 10−2), the increase in Rayleigh number changes isotherms from a distributed form (fan-shaped) to vertical lines stuck at the hot wall. When Re increases, we show a flattening of the thermal boundary layer, probably due to the large flow rate of the parietal fluid, which carries more energy. In the left-hand compartment of the cavity, a vertical stratification zone is formed which clearly illustrates the thermal boundary layer. Additionally, in the right compartment, horizontal thermal stratification is observed as the Reynolds number increases. So, isotherms pass from an almost constant temperature between 0.15 and 0.25 to a vertical temperature stratification along a range depending on the Rayleigh number.
For a low-permeability porous partition, Da = 10−6, it can be seen that the isotherms are not affected except for large Reynolds values where the flow structures are affected by the presence of the porous obstacle.
The effect of the thickness and height of the porous wall is illustrated in Figure 8, wherein it can be seen that the partition influences the isotherms in the cavity. This influence is more visible in the second compartment when the thermal draft is the maximum for the maximum height considered (Hp = 0.8). When the height is minimal, no influence is observed.
For average values of the Darcy number (Da = 10−2), the increase in Ra gives isotherms passing from a distributed form (aerated) to vertical lines grouped near the hot wall in the form of a thermal boundary layer. With increasing Re, this thermal boundary layer flattens out due to the high fluid flow. In the rest of the cavity (right section), a stratification zone is formed. Subsequently, when the Reynolds number increases, a horizontal thermal stratification is observed, which will give rise to a natural convection motion represented in the graphs of the flow functions. For the porous partition with low permeability (Da = 10−6), its presence is remarkable because at its level the isotherms lose their horizontal shape going towards Re > 250 and Ra > 104 and the stratification is visible on both sides of the partition with a slight difference in the range of variation that represents larger values in the left space of the cavity.
Figure 9 shows the transfer rate, an expression of the Nusselt number, as a function of the control parameters. On the abscissa axis, we find the Darcy number, which gives us the relative permeability of the porous wall. The different lines on the graph represent the variation in the driving parameters (i.e., Reynolds and Rayleigh numbers) of the problem for two heights of porous media considered (Ep = 0.3, Hp = 0.2 and Hp = 0.8). We can see very well on this figure that when the height of the porous wall is small (Hp = 0.2), its effect is negligible on the heat transfer and the curves have almost horizontal linear lines depending only on the intensities of the two natural and mechanical motors (i.e., Rayleigh and Reynolds numbers). This is certainly due to the fact that when the height of the porous wall is equal to 20% of the height of the chamber, its effect remains local and influences the transfers very little. It is shown in this figure that the heat transfer is directly dependent on the Rayleigh and Reynolds numbers.
When the height is the maximum (Hp = 0.8, on Figure 8), the heat transfer increases with a reduction in the permeability, but its influence remains weak (less than 10% increase). This can be explained by the fact that for low Darcy numbers, the flow is disturbed by a solid obstacle generating vortices that swirl the inner fluid, thus favoring heat transfer. Note that for the case of a porous partition height of 80%, we are in the presence of an optimization situation wherein an intermediate Reynolds of Re = 250 gives us the maximum heat transfer.
Figure 10 shows the results of numerical simulations of the heat transfer rate (Nusselt number) as a function of the geometric dimensions of the porous wall (Hp and Ep, respectively, height and thickness) and for the driving parameters of the ventilation system (i.e., Reynolds and Rayleigh numbers). Note that, in this figure, we have considered the case of a porous medium with low permeability (Da = 10−6). The shape of the curves is slightly increasing linearly, for the Rayleigh of large values, depending on the thickness of the porous wall for the height Hp = 0.8. This is due to the fact that, in this situation, the fluid is constricted in the first compartment and then directed towards the active wall, thus increasing the heat transfer. For the other Rayleigh cases, the increase in transfer is not visible and the plots show an almost horizontal linear trend. Note that optimization cases are found wherein a larger heat transfer is obtained for intermediate Reynolds numbers; (Re,Ra) = (250,10+6) gives a larger heat transfer than (500,10+6).
For the natural convection case (Ra = 106), the evolution of the Nusselt number is increasing linearly with the Darcy number on a log–log scale. Note that the effect of the porous medium (permeability, height, and thickness) is not remarkable, with the same result for the behavior of the curves, which tend to join.
The analysis of the heat transfer rate (Nusselt number) in the case of forced convection shows that it increases as the Reynolds number increases. This is obvious since the fluid is attracted to the heated wall (see graphs of the stream function), thus increasing the flow through the hot wall and subsequently causing an intensive extraction of hot fluid.

6. Conclusions

The flow structures and heat transfer have been numerically analyzed in a displacement-ventilated chamber. The conservation equations have been written by adopting some simplifying assumptions and then solved by the lattice Boltzmann method. The numerical results of the developed computational code have been compared with those of the literature. The findings of our research are quite convincing; thus, the following conclusions can be drawn:
  • The flow structures are complex and dependent on the control parameters where we find the spatial competition between a fluid jet of forced convection and the acceleration by natural convection at the hot wall;
  • When the permeability of the porous wall is high (Da = 10−2), the porous wall offers no resistance, which can be defined as a fluid regime;
  • When the permeability of the porous wall is very small (Da = 10−6), the porous wall behaves like a solid in this limit of the Darcien regime;
  • The porous wall has very little influence on the heat transfer, but its dimensions and permeability influence the temperature field, mainly in the second compartment;
  • The elevation of the separation thickness does not significantly change the curves evolution of the Nusselt number according to Ra and Re. The differences between these curves for different porous separation thickness are reduced, and the maximum values obtained for Re = 500 and Ra = 106.

Author Contributions

Conceptualization, Z.H. and D.E.A.; methodology, Z.H., N.H., L.N. and D.E.A.; software, L.N.; validation, Z.H. and N.H.; formal analysis, L.N., N.H. and D.E.A.; investigation, Z.H.; data curation, L.N.; writing—original draft preparation, Z.H. and D.E.A.; writing—review and editing, Y.H., Z.H. and D.E.A.; visualization, L.N.; supervision, D.E.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data sharing is not applicable to this article.

Acknowledgments

The authors acknowledge the national program of research PRFU.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

CLattice speed
csSpeed of sound (cs = 1/ 3 )
ciDiscrete velocities
DaDarcy number Da = K/H2
EpWidth of porous wall
FiExternal forces in Boltzmann equation
ƒ First population: distribution function of momentum equation
feqEquilibrium distribution function
GGravity (m/s2)
g Second population: distribution function of temperature
hCoefficient of convective heat exchange
HCavity height (m)
HpPorous separation height
JMomentum components
KPermeability of the porous medium
Fluid inlet width (m)
LCavity length (m)
m Moments of the mathematical model
meqEquilibrium moments of the mathematical model
MTransformation matrix of the mathematical model
N Thermal transformation matrix of the mathematical model
NuThe averaged Nusselt number
NuxThe local Nusselt number
PrThe Prandtl number Pr = υ/α
rjPosition node
RaThe Rayleigh number Ra = (gβ∆TrefH3)/υα
ReThe Reynolds number Re = U0H/υα
RiThe Richardson number Ri = Ra/(Pr·Re2)
SMatrix of the relaxation rates elements
tDimensionless time
TDimensional temperature
uDimensional velocity
UDimensionless horizontal velocity
VDimensionless vertical velocity
xDimensional longitudinal coordinate
XDimensionless longitudinal coordinate
yDimensional vertical coordinate
YDimensionless vertical coordinate
Subscript
CCold
EEffective
HHot
Greek symbols
αThermal diffusivity (m2/s)
β Thermal expansion coefficient (K−1)
ε Porosity of the porous medium (ε = 0,98)
ε ¯ The second-order energy
μ Dynamic viscosity (kg/ms−1)
υ Kinematic viscosity (m2/s)
δ t Temporal step
τ Relaxation time
ρDensity (Kg/m3)
ΘDiagonal relaxation matrix of σi
ϕ Energy flux
θDimensionless temperature

References

  1. De Vahl Davis, G. Natural convection of air in a square cavity: A benchmark numerical solution. Int. J. Numer. Methods Fluids 1983, 3, 249–264. [Google Scholar] [CrossRef]
  2. Oztop, H.F. Natural convection in partially cooled and inclined porous rectangular enclosures. Int. J. Therm. Sci. 2007, 46, 149–156. [Google Scholar] [CrossRef]
  3. Sathiyamoorthy, M.; Basak, T.; Roy, S.; Pop, I. Steady natural convection flows in a square cavity with linearly heated side wall. Int. J. Heat Mass Transf. 2007, 50, 766–775. [Google Scholar] [CrossRef]
  4. Varol, Y.; Oztop, H.F.; Pop, I. Numerical analysis of natural convection for a porous rectangular enclosure with sinusoidally varying temperature profile on the bottom wall. Int. Commun. Heat Mass Transf. 2008, 35, 56–64. [Google Scholar] [CrossRef]
  5. Ameziani, D.E.; Guo, Y.; Bennacer, R.; Ganaoui, M.; Bouzidi, M. Competition Between Lid-Driven And Natural Convection In Square Cavities Investigated With A Lattice Boltzmann Method. Comput. Therm. Sci. 2010, 2, 269–282. [Google Scholar] [CrossRef]
  6. Huelsz, G.; Rechtman, R. Heat transfer due to natural convection in an inclined square cavity using the lattice Boltzmann equation method. Int. J. Therm. Sci. 2013, 65, 111–119. [Google Scholar] [CrossRef]
  7. Omri, A.; Nasrallah, S.B. Control volume finite element numerical simulation of mixed convection in an air-cooled cavity. Numer. Heat Transf. Part A Appl. 1999, 36, 615–637. [Google Scholar]
  8. Najam, M.; El Alami, M.; Hasnaoui, M.; Amahmid, A. Étude numérique de la convection mixte dans une cavité en forme de «T» soumise à un flux de chaleur constant et ventilée par le bas à l’aide d’un jet d’air vertical. Comptes Rendus Mécanique 2002, 330, 461–467. [Google Scholar] [CrossRef]
  9. Gan, G. Simulation of buoyancy-driven natural ventilation of buildings—Impact of computational domain. Energy Build. 2010, 42, 1290–1300. [Google Scholar] [CrossRef]
  10. Tong, T.W.; Subramaniant, E. Natural convection in rectangular enclosures partially filled with a porous medium. Int. J. Heat Fluid Flow 1986, 7, 3–10. [Google Scholar] [CrossRef]
  11. Moraga, N.O.; Sánchez, G.C.; Riquelme, J.A. Unsteady mixed convection in a vented enclosure partially filled with two non-Darcian porous layers. Numer. Heat Transf. Part A Appl. 2010, 57, 473–495. [Google Scholar] [CrossRef]
  12. Mehrizi, A.A.; Sedighi, K.; Afrouzi, H.H.; Aghili, A.L. Lattice Boltzmann Simulation of Forced Convection in Vented Cavity Filled by Porous Medium with Obstruction. World Appl. Sci. J. 2012, 16, 31–36. [Google Scholar]
  13. Oztop, H.F.; Al-Salem, K.; Varol, Y.; Pop, I. Natural convection heat transfer in a partially opened cavity filled with porous media. Int. J. Heat Mass Transf. 2011, 54, 2253–2261. [Google Scholar] [CrossRef]
  14. Liu, Q.; He, Y.L. Multiple-relaxation-time lattice Boltzmann modeling of incompressible flows in porous media. Physica A 2015, 429, 215–230. [Google Scholar] [CrossRef] [Green Version]
  15. Shuja, S.Z.; Yilbas, B.S.; Khan, S.M.A. Flow subjected to porous blocks in the cavity: Consideration of block aspect ratio and porosity. Chem. Eng. J. 2008, 139, 84–92. [Google Scholar] [CrossRef]
  16. Shuja, S.Z.; Yilbas, B.S.; Kassas, M. Flow over porous blocks in a square cavity: Influence of heat flux and porosity on heat transfer rates. Int. J. Therm. Sci. 2009, 48, 1564–1573. [Google Scholar] [CrossRef]
  17. Hireche, Z.; Nasseri, L.; Ameziani, D.E. Heat transfer analysis of a ventilated room with a porous partition: LB-MRT simulations. Eur. Phys. J. Appl. Phys. 2020, 91, 20904. [Google Scholar] [CrossRef]
  18. Lage, J.L.; Bejan, A.; Anderson, R. Efficiency of transient contaminant removal from a slot ventilated enclosure. Int. J. Heat Mass Transf. 1991, 34, 2603–2615. [Google Scholar] [CrossRef]
  19. Lage, J.L.; Bejan, A.; Anderson, R. Removal of contaminant generated by a discrete source in a slot ventilated enclosure. Int. J. Heat Mass Transf. 1992, 35, 1169–1180. [Google Scholar] [CrossRef]
  20. Deng, Q.-H.; Zhou, J.; Mei, C.; Shen, Y.-M. Fluid, heat and contaminant transport structures of laminar double-diffusive mixed convection in a two-dimensional ventilated enclosure. Int. J. Heat Mass Transf. 2004, 47, 5257–5269. [Google Scholar] [CrossRef]
  21. Girault, M.; Maillet, D.; Bonthoux, F.; Galland, B.; Martin, P.; Braconnier, R.; Fontaine, J.R. Estimation of time-varying pollutant emission rates in a ventilated enclosure: Inversion of a reduced model obtained by experimental application of the modal identification method. Inverse Probl. 2008, 24, 015021. [Google Scholar] [CrossRef]
  22. Saidi, M.H.; Sajadi, B.; Molaeimanesh, G.R. The effect of source motion on contaminant distribution in the cleanrooms. Energy Build. 2011, 43, 966–970. [Google Scholar] [CrossRef]
  23. Xamán, J.; Ortiz, A.; Álvarez, G.; Chávez, Y. Effect of a contaminant source (CO2) on the air quality in a ventilated room. Energy 2011, 36, 3302–3318. [Google Scholar] [CrossRef]
  24. Serrano-Arellano, J.; Xamán, J.; Álvarez, G. Optimum ventilation based on the ventilation effectiveness for temperature and CO2 distribution in ventilated cavities. Int. J. Heat Mass Transf. 2013, 62, 9–21. [Google Scholar] [CrossRef]
  25. Du, R.; Gokulavani, P.; Muthtamilselvan, M.; Al-Amri, F.; Abdalla, B. Influence of the Lorentz force on the ventilation cavity having a centrally placed heated baffle filled with the Cu—Al2O3—H2O hybrid nanofluid. Int. Commun. Heat Mass Transf. 2020, 116, 104676. [Google Scholar] [CrossRef]
  26. Palaniappan, G.; Murugan, M.; Mdalla, Q.; Doh, B. Numerical investigation of open cavities with parallel insulated baffles. Int. J. Heat Technol. 2020, 38, 611–621. [Google Scholar] [CrossRef]
  27. Hsu, C.-T.; Cheng, P. Thermal dispersion in a porous medium. Int. J. Heat Mass Transf. 1990, 33, 1587–1597. [Google Scholar] [CrossRef]
  28. Seta, T.; Takegoshi, E.; Okui, K. Lattice Boltzmann simulation of natural convection in porous media. Math. Comput. Simul. 2006, 72, 195–200. [Google Scholar] [CrossRef]
  29. Hassan, S.; Himika, T.A.; Molla, M.M.; Hasan, F. Lattice Boltzmann simulation of fluid flow and heat transfer through partially filled porous media. Comput. Eng. Phys. Modeling 2019, 2, 38–57. [Google Scholar]
  30. Ergun, S. Fluid flow through packed columns. Chem. Eng. Prog. 1952, 48, 89–94. [Google Scholar]
  31. Cercignani, C.; Berman, A.S. Theory and Application of the Boltzmann Equation. J. Appl. Mech. 1976, 43, 521. [Google Scholar] [CrossRef]
  32. Frisch, U.; Hasslacher, B.; Pomeau, Y. Lattice-Gas Automata for the Navier-Stokes Equation. Phys. Rev. Lett. 1986, 56, 1505–1508. [Google Scholar] [CrossRef] [Green Version]
  33. Chapman, S.; Cowling, T.G. The Mathematical Theory of Non-Uniform Gases. Math. Gaz. 1939, 23, 488. [Google Scholar] [CrossRef]
  34. Enskog, D. Kinetische Theorie der Vorgënge in Mässigverdunntën Gasen. Ph.D. Thesis, Upsala University, Upsala, Sweden, 1917. [Google Scholar]
  35. Bhatnagar, P.L.; Gross, E.P.; Krook, M. A Model for Collision Processes in Gases. I. Small Amplitude Processes in Charged and Neutral One-Component Systems. Phys. Rev. 1954, 94, 511–525. [Google Scholar] [CrossRef]
  36. d’Humières, D. Generalized lattice Boltzmann equations. Prog. Astraunotics Aeronotics 1992, 159, 450–458. [Google Scholar]
  37. Shah, N.; Dhar, P.; Chinige, S.K.; Geier, M.; Pattamatta, A. Cascaded collision lattice Boltzmann model (CLBM) for simulating fluid and heat transport in porous media. Numer. Heat Transf. Part B Fundam. 2017, 72, 211–232. [Google Scholar] [CrossRef]
  38. Luo, L.-S. Lattice-Gaz Automata and Lattice Boltzmann Equations for Two-Dimensional Hydrodynamics. Ph.D. Thesis, Georgia Tech, Atlanta, GA, USA, 1993. [Google Scholar]
  39. Higuera, F.J.; Succi, S.; Benzi, R. Lattice Gas Dynamics with Enhanced Collisions. Europhys. Lett. (EPL) 1989, 9, 345–349. [Google Scholar] [CrossRef]
  40. Li, Q.; He, Y.L.; Tang, G.H.; Tao, W.Q. Improved axisymmetric lattice Boltzmann scheme. Phys. Rev. E 2010, 81, 056707. [Google Scholar] [CrossRef] [Green Version]
  41. Liu, Q.; He, Y.L.; Li, Q.; Tao, W.Q. A multiple-relaxation-time lattice Boltzmann model for convection heat transfer in porous media. Int. J. Heat Mass Transf. 2014, 73, 761–775. [Google Scholar] [CrossRef] [Green Version]
  42. Zou, Q.; He, X. On pressure and velocity boundary conditions for the lattice Boltzmann BGK model. Phys. Fluids 1997, 9, 1591–1598. [Google Scholar] [CrossRef] [Green Version]
  43. Mezrhab, A.; Amine Moussaoui, M.; Jami, M.; Naji, H.; Bouzidi, M. Double MRT thermal lattice Boltzmann method for simulating convective flows. Phys. Lett. A 2010, 374, 3499–3507. [Google Scholar] [CrossRef]
  44. Lallemand, P.; Luo, L.S. Theory of the lattice Boltzmann method: Acoustic and thermal properties in two and three dimensions. Phys. Rev. E 2003, 68, 036706. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Himrane, N.; Ameziani, D.E.; Nasseri, L. Study of thermal comfort: Numerical simulation in a closed cavity using the lattice Boltzmann method. SN Appl. Sci. 2020, 2, 785. [Google Scholar] [CrossRef] [Green Version]
  46. Nasseri, L.; Himrane, N.; Ameziani, D.E.; Bourada, A.; Bennacer, R. Time-periodic cooling of Rayleigh–Bénard convection. Fluids 2021, 6, 87. [Google Scholar] [CrossRef]
  47. Hireche, Z.; Nasseri, L.; Ameziani, D.E. Study of Periodic Thermal Exchange in a Cavity Ventilated by Displacement. Arab. J. Sci. Eng. 2020, 45, 5751–5768. [Google Scholar] [CrossRef]
  48. Nasseri, L.; Rahli, O.; Ameziani, D.E.; Bennacer, R. Study of mixed convection in closed enclosure with a ceiling fan. Eur. Phys. J. Appl. Phys. 2019, 86, 20902. [Google Scholar] [CrossRef]
  49. Rahman, M.M.; Alim, M.A.; Saha, S.; Chowdhury, M.K. Effect of the presence of a heat conducting horizontal square block on mixed convection inside a vented square cavity. Nonlinear Anal. Model. Control. 2009, 14, 531–548. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Physical problem.
Figure 1. Physical problem.
Computation 10 00004 g001
Figure 2. Validation of the evolution of the average Nusselt number as a function of Ri.
Figure 2. Validation of the evolution of the average Nusselt number as a function of Ri.
Computation 10 00004 g002
Figure 3. Stream function as function of Ra and Re for Da = 10−2, 10−4 and 10−6 with Ep = 0.2 and Hp = 0.4.
Figure 3. Stream function as function of Ra and Re for Da = 10−2, 10−4 and 10−6 with Ep = 0.2 and Hp = 0.4.
Computation 10 00004 g003
Figure 4. Stream function as function of Ra and Re for Da = 10−4 and 10−6 with Ep = 0.1 and 0.3 and Hp = 0.4.
Figure 4. Stream function as function of Ra and Re for Da = 10−4 and 10−6 with Ep = 0.1 and 0.3 and Hp = 0.4.
Computation 10 00004 g004
Figure 5. Stream function as function of Ra and Re for Da = 10−6 with Ep = 0.1–0.15 and Hp = 0.2 and 0.8.
Figure 5. Stream function as function of Ra and Re for Da = 10−6 with Ep = 0.1–0.15 and Hp = 0.2 and 0.8.
Computation 10 00004 g005
Figure 6. Evolution of isotherms as function of Ra and Re for Da = 10−2, 10−4 and 10−6 with Ep = 0.2 and Hp = 0.4.
Figure 6. Evolution of isotherms as function of Ra and Re for Da = 10−2, 10−4 and 10−6 with Ep = 0.2 and Hp = 0.4.
Computation 10 00004 g006
Figure 7. Evolution of isotherms as function of Ra and Re for Da = 10−4 and 10−6 with Ep = 0.1 and 0.3 and Hp = 0.4.
Figure 7. Evolution of isotherms as function of Ra and Re for Da = 10−4 and 10−6 with Ep = 0.1 and 0.3 and Hp = 0.4.
Computation 10 00004 g007
Figure 8. Evolution of isotherms as function of Ra and Re for Da = 10−6 with Ep = 0.1–0.3 and Hp = 0.2 and 0.8.
Figure 8. Evolution of isotherms as function of Ra and Re for Da = 10−6 with Ep = 0.1–0.3 and Hp = 0.2 and 0.8.
Computation 10 00004 g008
Figure 9. Evolution de Nusselt number as function of Da for different values of Re, Ra and Hp (Ep = 0.3).
Figure 9. Evolution de Nusselt number as function of Da for different values of Re, Ra and Hp (Ep = 0.3).
Computation 10 00004 g009
Figure 10. Evolution of Nu as function of porous dimensions (Hp and Ep) for different values of Re and Ra (Da = 10−6).
Figure 10. Evolution of Nu as function of porous dimensions (Hp and Ep) for different values of Re and Ra (Da = 10−6).
Computation 10 00004 g010
Table 1. Dimensionless number expressions.
Table 1. Dimensionless number expressions.
TypesExpression
Prandtl number Pr = υ / α
Reynolds number Pr = U 0 H / υ
Rayleigh number Ra = ( g β T r e f H 3 ) / υ α
Darcy number D a = K / H 2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hireche, Z.; Himrane, N.; Nasseri, L.; Hamrioui, Y.; Ameziani, D.E. Analysis of Thermal Performances in a Ventilated Room Using LBM-MRT: Effect of a Porous Separation. Computation 2022, 10, 4. https://doi.org/10.3390/computation10010004

AMA Style

Hireche Z, Himrane N, Nasseri L, Hamrioui Y, Ameziani DE. Analysis of Thermal Performances in a Ventilated Room Using LBM-MRT: Effect of a Porous Separation. Computation. 2022; 10(1):4. https://doi.org/10.3390/computation10010004

Chicago/Turabian Style

Hireche, Zouhira, Nabil Himrane, Lyes Nasseri, Yasmine Hamrioui, and Djamel Eddine Ameziani. 2022. "Analysis of Thermal Performances in a Ventilated Room Using LBM-MRT: Effect of a Porous Separation" Computation 10, no. 1: 4. https://doi.org/10.3390/computation10010004

APA Style

Hireche, Z., Himrane, N., Nasseri, L., Hamrioui, Y., & Ameziani, D. E. (2022). Analysis of Thermal Performances in a Ventilated Room Using LBM-MRT: Effect of a Porous Separation. Computation, 10(1), 4. https://doi.org/10.3390/computation10010004

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop