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 (10
3 < Ra < 10
5) 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.
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:
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:
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:
Such a function therefore verifies:
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:
where
is the distribution function with the velocity
at the lattice node
x at time
t,
is the discrete time step,
is the discrete collision operator, and
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:
τ: 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:
The discrete velocity of Boltzmann equation can be written as follows:
The discretization of the previous equation gives us:
- (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]:
Therefore, the Boltzmann equation is written as:
The velocity is equal to the ratio between displacement dx and the unit time dt.
are moment vectors;
. The relaxation matrix is a diagonal matrix, and the linkage between velocities and moment spaces is given from linear transformation as follows:
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
, 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]:
The nine discrete velocities are given by:
where
is the lattice velocity, in which
and
are the lattice cell widths; these quantities are chosen equal to the unit
.
The distribution function of local equilibrium
for the hydrodynamic field in porous media is calculated by equation:
where the weighting factors
wi are given in the literature as follows:
.
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:
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]:
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:
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:
In the LBM-MRT method, the fluid kinematic viscosity is related to the relaxation parameter of the flow field by the following relation:
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]:
The conservation of momentum can be described as the summation of the product of the microscopic lattice velocities by the distribution function:
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:
where
gi(
x,
t) is the temperature distribution function. In the D2Q5 model, the five discrete speeds are given by the following relation:
The equilibrium distribution function is given by:
Here, the weighting factors wi are given as follows; .
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]:
E is the diagonal relaxation matrix:
To ensure good stability of computation, the relaxation rates are chosen in our work as follows:
where
M is the relaxation parameter of temperature field.
The collision process D2Q5 LBM-MRT proceeds as follows:
Between the five moments, only the temperature is kept:
In equilibrium, the non-conserved moments
are written as follows [
44]:
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.
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 10
2–10
6 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 = 10
6) 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 [10
2–10
6]. 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 = 10
6), 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.