A Comprehensive Review of Nanofluid Heat Transfer in Porous Media

Abstract

In the present paper, recent advances in the application of nanofluids in heat transfer in porous materials are reviewed. Efforts have been made to take a positive step in this field by scrutinizing the top papers published between 2018 and 2020. For that purpose, the various analytical methods used to describe the flow and heat transfer in different types of porous media are first thoroughly reviewed. In addition, the various models used to model nanofluids are described in detail. After reviewing these analysis methods, papers concerned with the natural convection heat transfer of nanofluids in porous media are evaluated first, followed by papers on the subject of forced convection heat transfer. Finally, we discuss articles related to mixed convection. Statistical results from the reviewed research regarding the representation of various parameters, such as the nanofluid type and the flow domain geometry, are analyzed, and directions for future research are finally suggested. The results reveal some precious facts. For instance, a change in the height of the solid and porous medium results in a change in the flow regime within the chamber; as a dimensionless permeability, the effect of Darcy’s number on heat transfer is direct; and the effect of the porosity coefficient has a direct relationship with heat transfer: when the porosity coefficient is increased or decreased, the heat transfer will also increase or decrease. Additionally, a comprehensive review of nanofluid heat transfer in porous media and the relevant statical analysis are presented for the first time. The results show that Al2O3 nanoparticles in a base fluid of water with a proportion of 33.9% have the highest representation in the papers. Regarding the geometries studied, a square geometry accounted for 54% of the studies.

1. Introduction

A porous medium or a porous material is a solid material that contains pores. Depending on the ability of the porous medium to allow fluids to pass through it under the influence of external forces, it is classified as a permeable or a nonpermeable porous medium. Porous media may be non-dispersed or post-dispersed, homogeneous or heterogeneous, and multi-structure or the result of a combination of different structures [1]. Porosity, the main characteristic property of a porous medium, is a measure of the empty spaces in a material and is a fraction of the volume of the empty spaces over the total volume between 0 and 1. Most natural porous materials have a porosity of 0.6 (excluding hair). However, manufactured materials, such as metal foams, can have a porosity of up to 0.99. Table 1 shows the porosity value for different materials [2].

Investigating heat transfer characteristics in porous media is of great interest in various industries and engineering domains such as heat exchangers, heat storage, geothermal systems, and drying techniques. Including a porous material in a mechanical system can be considered a passive method for heat transfer enhancement. Indeed, the presence of a porous material alters the flow patterns and improves the overall thermal conductivity of the system [3,4,5,6,7].

A further improvement in heat transfer can be achieved by adding conductive nanoparticles to the base fluid. The working fluid is then considered a nanofluid. The added nanoparticles can be metallic and made of metals such as aluminum and copper, among others, or metal oxides, or non-metallic materials such as carbon. Hybrid nanofluids refer to the fluids in which two or more types of nanoparticles are dispersed. The main objective of dispersing nanoparticles in the fluid is to increase its thermal conductivity [8]. Nonetheless, adding nanoparticles to a fluid presents a drawback that should be avoided. Increasing the volume fraction of the nanoparticles above a certain level increases the fluid viscosity and may, as a result, hinder heat transfer.

Due to the double importance of nanofluids and porous materials in heat transfer enhancement, significant research has been performed to study the thermal characteristics of nanofluid flow in porous media. The present paper summarizes the recent papers dealing with this topic. Considering the large number of articles published since the last review article (2017 to the present) and their analytical complexities, the authors divide the articles into two parts: 2018–2020 and 2020 until the present. These parts are presented in two articles. The paper is organized as follows: in Section 2, the various models of fluid flow and heat transfer in porous media are recalled and the classifications of the nanofluids are presented. In Section 3, the works dealing with heat transfer in porous media are presented and are classified into free convection and forced convection. Finally, Section 4 presents a general conclusion.

2. Methods and Materials

2.1. Methods of Analysis of Porous Material

In this section, the various analytical methods used to describe porous media are presented based on the scale length of the medium.

2.1.1. Microporous Medium

Microporous media are at the size scale of the nanometer. Examples of such media include activated carbon, silica gels, carbon molecular sieves, and some crystalline structures such as zeolites.

2.1.2. Mesoporous Media

The pore size of mesoporous media is between 2 and 50 nm. This size applies to non-organic jellies such as alumina, silica powders, porous glass, and columnar or non-columnar bricks [2]. Since this review considers the heat transfer of nanofluids in porous media and nanofluids have particles with a diameter of more than 10 nm, the two methods described above cannot be used for this purpose. Therefore, we describe the macroscopic method in the following sections.

2.1.3. Macroporous Media

Macroporous media have a size scale of greater than 50 nanometers. Macroporous media are widely found in nature, such in soil, broken rocks, sandstones, wood materials, and various types of food. Foods are generally macroporous but have different porosity size properties in some cases. This feature is indicated by the distribution of pore size, which indicates two or more sizes. Manufactured materials such as thermal insulation materials, silicate, ceramics, cement, synthetic resins, and many other artificial materials are also macroporous materials. They have many applications, such as in fluid filters, electronic components, complex fiber structures, bioceramics, fluid chromatography, and biotechnology [2].

2.1.4. Macroscopic Governing Equations

The macroscopic equation is written as follows [2]:

(1)$\widehat{\nabla }.\mathrm{v}=0$

Considering the Boussinesq assumption, the macroscopic momentum is as follows:

(2)${\rho }_f\left[\frac{\partial v}{\partial t}+\widehat{\nabla }.\left(\frac{vv}{\varnothing }\right)\right]=-{\widehat{\nabla }}_p+{\mu }_f{\widehat{\nabla }}^{2}v+B-{\rho }_f\varnothing B\left({\widehat{T}}_f-{\widehat{T}}_{\infty }\right)g$

where

(3)$B=-\frac{1}{v}{\int }_{A_{fs}}{p}_{f}dS+\frac{{\mu }_f}{v}{\int }_{A_{fs}}\left(\nabla v_f\right).dS$

and

(4)$p=\varnothing {\widehat{p}}_f$

In these formulas $p,\varnothing ,\mathsf{\rho },\mathrm{and}\mathsf{\mu }$ are the pressure, porosity, density, and viscosity, respectively. The subscripts s and f are related to the solid and fluid.

The macroscopic energy equations for the fluid and solid are as follows:

(5)$\varnothing {\left(\rho C_p\right)}_f\left[\frac{\partial {\widehat{T}}_f}{\partial t}+\widehat{\nabla }.\left({\widehat{v}}_f{\widehat{T}}_f\right)\right]=\varnothing \widehat{\nabla }.\left[\left(k_f+k^{\prime }\right)\widehat{\nabla }{T}_f\right]+q_{sf}$

(6)$\left(1-\varnothing \right){\left(\rho C_p\right)}_s\frac{\partial {\widehat{T}}_s}{\partial t}=\left(1-\varnothing \right)\left[\widehat{\nabla }.\left(k_s\widehat{\nabla }{\widehat{T}}_s\right)\right]-q_{sf}$

In these formulas, $q_{sf}$ is the heat transfer between the solid matrix and fluid flow, and k represents the conductivity.

2.2. Nanofluid

Generally, a nanofluid is obtained by dispersing a certain volume fraction of nanoparticles in a base fluid. Nanoparticles are produced at different sizes from 10 nm to 100 nm, depending on their application. Since most of the fluids used for heat transfer have a low conductivity, heat transfer can be significantly improved with a uniform distribution of nanoparticles within the base fluid. Therefore, the typical size of nanofluid particles is greater than the typical pore size of microporous media and of the same order of magnitude as the typical pore size of mesoporous media. Therefore, nanofluids can only be used in macroporous media. Figure 1 shows a picture of the Titania nanoparticles [4].

Some of the exceptional properties of nanoparticles include their non-linear relationship between the conductivity and the concentration of solids, the low momentum of the particles, a higher mobility than microparticles, the strong dependence of conductivity on temperature, a strong increase in heat flux in the boiling region, and acceptable viscosity. These properties, provided to the fluids, are considered some of the most suitable and strongest choices.

• Thermal conductivity: the coefficient of the thermal conductivity of nanofluids depends on parameters such as the composition of the chemical percentage of nanoparticles, the volume percentage of nanoparticles, the surface-active substances, and the temperature. The coefficient of thermal conductivity is also influenced by mechanisms such as the Brownian motion of nanoparticles in the fluid, which increases mixing in the fluid, facilitates heat transfer and increases the coefficient of thermal conductivity;

• Size reduction: the small size of the nanoparticles reduces their motion and contact with the solid wall, reduces momentum, and ultimately reduces the possibility of erosion of parts such as heat exchangers, pipelines, and pumps;

• Stability: nanoparticles are less likely to be precipitated due to their low weight and small size, which prevents the problem of nanoparticle suspension caused by sedimentation. Presently, the instability of nanofluids hinders the application of nanofluids. The stability of the nanofluid means that the nanoparticles do not accumulate and precipitate at a significant rate and, as a result, the concentration of the floating nanoparticles is constant. Stokes’ law can be used to calculate the settling velocity of spherical particles in a quiescent fluid. This equation is obtained from the balance of gravity, buoyancy, and drag forces that act on particles. The stability of the nanofluid is a necessary condition for optimizing the properties of the nanofluid. Three general methods to increase nanofluid stability are:

  • Adding a surfactant;

  • Controlling the pH of the nanofluid;

  • Ultrasonic vibration.

2.2.1. Nanofluid Evaluation Methods

Researchers use different perspectives to analyze nanoscale behavior. Therefore, the modeling and formulation of nanoscale behavior also differs based on these perspectives. From one point of view, a certain volume fraction of nanoparticles is combined with the base fluid, but the nanofluid concentration does not change with respect to the initial concentration. The concentration of the nanofluid remains constant in different areas, but the newly formed fluid (nanofluid) has improved thermophysical properties compared to the base fluid. On the other hand, the movement of nanoparticles relative to the base fluid is not considered. Hence, this method or perspective is called a homogeneous or single-phase method. In the homogeneous method, three equations of mass conservation, momentum conservation, and energy conservation form the main structure of the model, and the thermophysical properties of the nano-fluid are substituted in conservation equations. The second method of analyzing nanoscale behavior is somewhat different from the homogeneous method. In this method, the movement of the nanoparticles changes the volume fraction of the nanoparticles (relative to the initial concentration) in different areas of the flow domain. Therefore, this method is called a non-homogeneous or pseudo-two-phase method. Several factors, such as gravitational force, Brownian forces, and thermophoresis forces, affect the heat transfer in this model. However, the Brownian and thermophoresis forces are dominant forces, according to ref. [8]. Because nanofluids consist of solid particles and a base fluid, one of the modeling approaches is to solve the conservation equations for the base fluid and the nanoparticles separately and is called the two-phase model.

2.2.2. Hybrid Nanofluid

Another type of nanofluid in which two or more nanoparticles exist within the base fluid is called a hybrid nanofluid [6]. When one nanoparticle type is dispersed in a combination of two or more fluids it is also called a hybrid nanofluid. The heat conductivity is improved significantly when using this type of nanofluid. For example, although the addition of aluminum oxide nanoparticles to water increases the conductivity, the addition of copper nanoparticles to the aluminum–water nanofluid can further increase the conductivity relative to the base fluid. Figure 2 shows the aluminum and copper hybrid nanofluid presented in ref. [9].

2.3. Fluid Flow Models

2.3.1. Darcy’s Equation

Based on the experimental data, Darcy’s equation provides a linear proportion between the fluid’s volumetric average velocity $\left(v\right)$ and its pressure difference $(∆p$) in addition to a porous media.

(7)$v=k\frac{∆p}{∆x}$

In this formula, $\frac{∆p}{∆x}$ is the pressure gradient. Hydraulic conductivity is an expression of how easily a fluid can flow throughout the hollows of a medium. Darcy’s equation is valid for creeping (viscous fluids having slow motion), isothermal incompressible flows [2].

2.3.2. Hazen–Darcy Equation

Darcy’s experiments were performed with a single-type fluid at a constant temperature; therefore, Darcy’s equation did not include the fluid’s viscosity,$\mu$. By providing a specific permeability relation, $K=k\mu$, Darcy’s equation is redefined in a viscosity-dependent form as follows:

(8)$v=\left(\frac{K}{\mu }\right)\frac{∆p}{∆x}$

The specific permeability $K$ is considered a hydraulic parameter and is independent of the fluid’s characteristics. The above-mentioned formula is the Hazen–Darcy equation, which is known as Darcy’s law [2].

2.3.3. Hazen–Dupuit–Darcy Equation

The quadratic Hazen–Dupuit–Darcy equation is derived from the analysis of steady-state, open-channel flows. The equation is defined based on an equilibrium state between the gravitational force and the shearing resistance as follows:

(9)$0=\frac{\partial p}{\partial x}+\frac{\mu }{K}v-C\rho v^2$

The second and third terms on the right-hand side express the viscous drag and geometrical drag, respectively. It is worthwhile to mention that the equation is only valid for steady-state and one-dimensional flows [2].

2.3.4. Brinkman–Hazen–Dupuit–Darcy Equation

Brinkman realized that when the permeability is low, the shearing tension of the fluid’s viscosity can be negligible in comparison to the viscous drag. Therefore, by adding the shearing tension term (a Laplacian one), the Brinkman–Hazen–Dupuit–Darcy equation is derived.

(10)$0=-\nabla p+\mu {\nabla }^{2}v-\frac{\mu }{K}\varphi v+\frac{C_F}{K^{1/2}}\rho {\varphi }^2\left|v\right.\left|v\right.$

where $\left|v\right|$ is the magnitude of the velocity. The fluid is capable of transforming the viscous shearing tension independent of the viscous drag. Both the K and C_F parameters have a wide variety in porous media in the engineering paradigms; for example, as packed beds, windowing environments, metal foams, and aerodynamic gels [2].

2.3.5. Brinkman–Forchheimer Equation

The modified equation is as follows.

(11)$\mathsf{\rho }\left[\frac{1}{\mathsf{\phi }}\frac{\partial \mathit{v}}{\partial t}+\frac{1}{{\phi }^2}\left(\mathit{v}.\nabla v\right)\right]=-\nabla p+{\mu }_e{\nabla }^2\mathit{v}-\frac{\mu }{K}\mathit{v}-\frac{c_F\rho }{K^{\raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{$2$}\right.}}\mathit{v}v$

The equation is applied for an incompressible fluid. The inertia term on the left-hand side is calculated by conventional averaging. The first viscosity term in the Brinkman equation and the last term in Forchheimer’s equation, the $C_F,$ is the Forchheimercoefficient. For more information about the mentioned formula, please see ref. [2].

2.4. Heat Transfer Models

Generally, there are two heat transfer models based on the volumetric averaging method: the single-equation model, which is based on the LTE (local thermal equilibrium) assumption, and the two-equation model, which is based on the LNTE (the local, non-thermal equilibrium) assumption. The two models will be discussed separately in terms of their assumptions and their application limitations [2].

2.4.1. LTE

By distinguishing the gradient operator in the microscopic and macroscopic coordinates, a simple form of the volumetric average is provided. By averaging the microscopic equations on a representative elementary volume (REV), the macroscopic equations for the mandatory relocation of an incompressible flow in a variant porous medium are derived.

(12)${\left(\rho C_p\right)}_m\frac{\partial \widehat{T}}{\partial t}+\overline{\nabla }.\left[v\widehat{T}\right]={\alpha }_f\left(\frac{k_m}{k_f}+\frac{k^{\prime }}{k_f}\right){\widehat{\mathit{\text{∇}}}}^2\widehat{T}$

According to the LTE assumptions, the single-equation model regarding the energy equation is highly capable of saving calculation time. The LTE assumption is that the temperature difference between the solid object and the fluid is negligible or is very small in comparison to the whole system’s temperature difference. This hypothesis can describe the inaccuracies in the energy transfer model.

2.4.2. LNTE

If the heat transfer between the solid object and the liquid is permissible, it leads to:

(13)$\left(1-\phi \right){\left(\rho c\right)}_s\frac{\partial T_s}{\partial t}=\left(1-\phi \right)\nabla .\left(k_s\nabla T_s\right)+\left(1-\phi \right)q_s^{‴}+h\left(T_f-T_s\right)$

(14)$\mathsf{\phi }{\left(\rho c_P\right)}_f\frac{\partial T_f}{\partial t}+\left(\rho c_P\right)v.\nabla T_f=\phi \nabla \left(k_f\nabla T_f\right)+\phi q_f^{‴}+h\left(T_s-T_f\right)$

In these formulas, $q^{‴}$ represents the volumetric heat generation, and h represents the convection heat transfer coefficient between the solid matrix and fluid flow.

2.4.3. Buongiorno’s Heterogeneous Model

Heat transfer among nanoparticles due to the effects of two phenomena, the thermophoresis sliding velocity and Brownian motion, results in a kind of heterogeneity in the nanofluid which is known for transfer or as Buongiorno’s model. In Buongiorno’s heterogeneous model, the thermophoresis and Brownian forces are the dominant forces exerted on the nanoparticles.

The thermophoresis force acts against the temperature gradient and tends to carry the nanoparticles from warm regions to cold ones. In contrast, the Brownian motion tends to transfer the nanoparticles from high-concentration regions to low-concentration regions. Due to the movement of nanoparticles in the base fluid, two important effects appear. The first is that the nanofluid in the low-concentration regions is lightweight and tends to move upward, while the nanofluid in the high-concentration regions is heavy and tends to move downward. The second is that the movements of the nanoparticles lead to energy transfer due to mass transfer [10].

The continuity equation, or the concentration of the particles using Buongiorno’s model in dimensional form, is as follows:

(15)$\frac{1}{\epsilon }V.\nabla \phi =\nabla .[D_B\nabla \phi +\frac{D_T}{T}\nabla \phi ]$

$D_B$ and $D_T$ are the Brownian motion and thermophoresis coefficient, respectively:

(16)$D_B=\frac{k_{B}T}{3\pi \mu d_{np}}$

(17)$D_T=(\frac{0.26k}{2k+k_{np}})(\frac{\mu }{\rho })\phi$

where ${\mathrm{k}}_{\mathrm{B}}{\mathrm{and}\mathrm{d}}_{\mathrm{np}}$ are the Boltzmann constant and the particle diameter, respectively.

For LTE, the energy equation in porous media is:

(18)$\rho c_{p}V.\nabla T=\nabla .k_m\nabla T+\epsilon {\rho }_{np}{c}_{p,np}[D_B\nabla \phi .\nabla T+D_T\frac{\nabla T.\nabla T}{T}]$

where $k_m$ is the effective thermal conductivity of the porous medium.

3. Results and Discussion

3.1. Free Convection Heat Transfer

When there is a temperature difference between a hot source and its surrounding fluid, the fluid’s density varies in terms of the temperature variation and the buoyancy forces come into play. Therefore, the fluid begins to flow on the solid’s surface in accordance with the temperature difference. In this case, heat transfer that occurs without any external stimulus and is completely autonomous is called free or natural convectional heat transfer.

Considering the importance of heat transfer, all its aspects and solutions should be taken into account to maximize the heat transfer speed. Natural heat transfer happens through the random motions of the molecules and the bulk motion of the fluid. Therefore, utilizing compounds with high heat conduction ratios in comparison to the fluids is useful. In this regard, using tiny particles with high heat conduction is considered. Studies showed that the smaller the size and dimension of the particles, the greater the effect on heat transfer rate [9,10,11,12,13].

3.2. Integrated Free Convectional Heat Transfer

Integrated free convectional heat transfer makes sense when two phenomena are combined: first, heat transfer due to the fluid’s displacement, and conduction heat transfer due to the contact between the fluid and the solid object. The type of medium that can host the integrated free convectional heat transfer phenomenon is a porous medium containing a fluid and nanoparticles. The integrated free convectional heat transfer phenomenon has several industrial applications, such as for particle storage, filtering, gas drying, underground pollution, maintaining cooling radioactive waste containers, soil cleaning by steam injection, heat insulation for buildings, solar collector technologies, electronic cooling, and many more [14].

Yekani-Motlagh et al. [15] investigated the free convection of a two-phase nanofluid in an inclined, porous, semi-annulus enclosure. They used Fe₃O₄-water as their magnetic nanofluid and used Darcy and Boungiorno’s models. They concluded that the Nusselt number increases when the nanoparticle volume fraction is increased.

The characteristics of the research papers related to free convection heat transfer and their results are summarized in Table 2.

3.3. Forced Convection Heat Transfer

Forced convection occurs when the fluid motion is generated by an external source and not only by the density difference inside the fluid.

Sheikholeslami et al. [43] performed analyses on the effects of Reynolds number, the volumetric quotient of the water–calcium oxide nanofluid, and the Hartmann and Darcy numbers on the forced convectional heat transfer in a container of a hot liquid. They realized that the heat profile decreases with the increase of the voltage, but heat conduction is improved by an increase in the Darcy and Reynolds numbers. In addition, the studies indicate that there is an inverse relationship between the temperature gradient and the Hartmann number.

Furthermore, Sheikholeslami et al. [44] evaluated the forced convectional heat transfer of water–aluminum oxide in the presence of a magnetic field and realized that the nanofluid velocity profile is in direct relation to the Reynolds number of the volumetric quotient of aluminum oxide but in an inverse relation to the Hartmann number. In addition, by the increase in the Lorentz force, the convectional heat transfer decreases. It was also revealed that the temperature gradient on the moving surface increases as the hot surface velocity and the volumetric aluminum oxide quotient increase. They reported that increasing the Reynold number resulted in the Nusselt number increasing. In addition, increasing the Hartmann number augmented heat transfer.

Ferdows and Alzahrani [45] numerically investigated the possibility of similar solutions as well as dual-branch solutions to evaluate the performance of nanoparticles associated with water: the base water. A steady-state condition was considered through a moving, flat, porous plate in the presence of magnetic fields. They concluded that the largest velocity profile and the smallest temperature distribution refer to the Cu-water nanofluid.

Table 3 summarizes several studies dealing with the forced convection of nanofluids.

3.4. Mixed Convection

Khademi et al. studied the mixed convection of nanofluids on a sloped, flat surface in a porous medium with the presence of a magnetic field. The boundary layer equations were solved numerically using the DQM method. This research identified a decrease in the Nusselt number [53].

Bondarenko et al. evaluated the mixed convection heat transfer of a nanofluid in a square container with insulated lateral walls and cold top and bottom walls. The boundary layer equations were solved using the FDM method [54]. Various research papers that considered the mixed convection of nanofluids are presented in Table 4.

3.5. Overall Review of Papers

Considering that the most important purpose of adding nanoparticles to a base fluid is the possibility of increasing the heat transfer, in the following section and in Table 5, the selected articles presented in Table 2, Table 3 and Table 4 in the previous section were examined and studied based on the increase or decrease of the Nusselt number, and a report of the increase or decrease in the Nusselt number is presented as a percentage. In addition, due to the importance of the size of nanoparticles, their sizes are presented in a separate column.

In most articles about the natural and mixed convection of a nanofluid in porous media, the Nusselt number is increased by adding the volume fraction of nanoparticles; however, some articles also reported a decrease. When nanoparticles are added to the base fluid, the conductivity definitely increases. Viscosity also increases; therefore, in low-velocity flows, such as free displacement heat transfer, the decrease in speed due to the increase in viscosity is quite evident and effective and it can dominate the increase in conductivity, causing a decrease in displacement heat transfer. The same applies to a low-speed combined heat transfer. In forced heat transfer, the increase in nanoparticles increases the heat transfer. Additionally, in some cases in which an increase in the nanoparticle diameter is observed, the natural convection heat transfer rate is increased by more than 5%. By increasing the nanoparticle diameter in the articles, the forced convection heat transfer rate was decreased.

3.6. Statistical Results

In this section, the statistical distribution of different parameters in the published papers on nanofluids in porous materials is presented. By reviewing the published research, it was revealed that Al₂O₃ nanoparticles in a base fluid of water with a proportion of 33.9% have the highest representation in the papers. Copper nanomaterials and Fe₃O₄-water, with a representation of 32.14% and 12.5%, respectively, occupy the second and third ranks. The share of each nanofluid is depicted in Table 6.

In general, it can be said that the most popular nanoparticle is the alumina nanoparticle due to its very good dispersion in the base fluid. Of course, the role of the stability of the nanofluid in the porous media is more important. Metal oxide nanoparticles have shown a dual behavior in free convection heat transfer which has been reported to increase in some cases and decrease in some cases. However, in forced convection heat transfer, the increase of any type of nanoparticle increases the heat transfer.

Regarding the studied geometries, a square geometry accounted for 54% of the studies. Hole and circular geometries were represented with a share of 17% and 13%, respectively, and occupy the next ranks. The details are depicted in Table 7.

Figure 3 illustrates the proportion of the models that were used for evaluation. Based on the figure, Darcy–Brinkman–Forchheimer has the most share at 16%, followed by the Darcy–Brinkman and Darcy models, with shares of 6% and 2%, respectively.

Figure 4 illustrates the share of the models used in the presence and absence of a magnetic field.

Figure 5 depicts the contribution of scientific publishers Science Direct, Springer, Tandf Online, John Wiley, and ASME in published papers related to this specific topic. According to this figure, the highest share of the studies was indexed in Science Direct at 56%.

In the statistical analysis section, the items that are important for heat transfer, according to the authors, were examined:

• The type of nanoparticle;

• The geometry under consideration.

Additionally, considering that the major contribution of the articles was related to numerical modeling, the model selected for investigating the flow in the porous medium was analyzed. It may be said that the last two forms are not directly related to the subject; however, as we know, external forces, including the magnetic force, have an effect on heat transfer, and approximately 30% of the subject articles dealt with it.

4. Conclusions

In this paper, works published between 2018 and 2020 that concern nanofluid heat transfer within porous media are reviewed. As the latest review in this field was published at 2017, it can be found at ref. [69]. Additionally, a good book was published in this field [70]. However, in this article, a statistical analysis of the published articles was performed. This has not been done before.

The results showed that heat transfer is improved by adding nanoparticles to the base fluid. On the contrary, by increasing the volume concentration of nanoparticles, the viscosity effect overcomes the thermal conductivity and hinders the heat transfer. The heat transfer is also reduced with an increase in the nanoscale conductivity relative to the porous matrix, and heat transfer is enhanced with the rise of the Darcy number and the porosity coefficient. Some important results emerged from this study. They are presented as follows:

• As a dimensionless permeability, hte effect of the Darcy number on heat transfer is direct. Therefore, as the Darcy number increases, the heat transfer increases. As the Darcy number decreases, the heat transfer decreases;

• The effect of the porosity coefficient has a direct relationship with the heat transfer, and the heat transfer will increase or decrease by increasing or decreasing the porosity coefficient;

• A change in the width of the solid and porous medium results in a change in the flow regime inside the chamber and has an increasing or decreasing effect on heat transfer;

• A change in the height of the solid and porous medium results in a change in the flow regime within the chamber and has an increasing or decreasing effect on heat transfer;

• Increasing the Ra number increases the heat transfer.

A detailed analysis of the share of each nanoparticle type and geometrical configuration revealed a disparity among the published papers. Future research should be directed towards the nanofluids that were less frequently investigated and should address new types of nanoparticles and hybrid nanofluids. Similarly, novel geometrical configurations should be subjected to more examination, such as circular and triangular enclosures. Finally, the economic aspect has not been widely considered for the various engineering applications. Studies addressing such aspects and evaluating the energy and long-term costs and gains of implementing nanofluids instead of conventional fluids in applications involving porous media should be emphasized in the future.

Finally, the main shortcomings inherent to the use of nanofluids in porous media are the increasing pressure drop, especially in forced convection, and the stability of the nanofluids (reduction due to the impact of nanoparticles on the solid matrix). It is proposed that research is performed on these problems so that they may be solved.

Author Contributions

Conceptualization, H.A.N., B.A., T.A., A.M.R., A.J.C.; Supervision, T.A. and H.A.N.; Investigation, H.A.N., B.A., T.A., A.M.R., A.J.C.; Methodology, H.A.N., B.A., T.A.; Writing—original draft, H.A.N., B.A., T.A.; Writing—review & editing, H.A.N., T.A., A.M.R., A.J.C.; Software, B.A., T.A.; Formal Analysis, H.A.N., B.A., T.A., A.M.R., A.J.C.; Project Adminstration H.A.N. All authors have read and agreed to the published version of the manuscript.

Data Availability Statement

Data are available upon request.

Acknowledgments

This study was supported by the Prince Sattam bin Abdulaziz University through project number (2022/RV/14).

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Boltzmann constant	$k_B$
Bejan Number	Be
Brownian motion	DB
Conduction factor for base fluid	$k_{bf}$
Conduction factor for nanoparticles	$k_p$
Darcy velocity	v
Diameter of the particle	$d_p$
Effective viscosity	${\mu }_e$
Fluid’s viscosity	$\mu$
Forchheimer’s dimensionless factor	CF
Hydraulic conduction	k
Magnetic field	M
Mass quotient of each phase	$C_k$
Nanofluid volume fraction	𝜑
Nanoparticle	np
Porosity	φ
Pressure	P
Pressure in saturated state	Ps
Specific permeability	K
Thermophoresis Coefficient	DT
Volumetric average velocity	v

Abbreviations

AR	Aspect Ratio
CFD	Computational Fluid Dynamics
CNT	Carbon nanotube
Da	Darcy
DQM	Differential Quadrature Method
FDM	Finite Difference Method
FEM	Finite Element Method
FVM	Finite Volume Method
Ha	Hartman Number
ISPH	Incompressible Smoothed Particle Hydrodynamics
KKL	Koo–Kleinstreuer–Li correlations
Le	Lewis Number
LNTE	Local Non-Thermal Equilibrium
LTE	Local Thermal Equilibrium
MFD	magnetic field dependent
MHD	Magnetohydrodynamics
NEPCM	Nano-Encapsulated Phase Change Materials
Nu	Nusselt Number
Pr	Prandtle Number
Ra	Raighly Number
Re	Reynolds Number
Ri	Richardson Number
SIMPLE	Semi-Implicit Method for Pressure Linked Equations
V. F.	Volume Fraction

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Figures and Tables

Figure 1. Titania nanoparticles (30 nm).

Figure 2. Aluminum and copper hybrid nanofluid.

Figure 3. The contribution of each of the models used to model the system in the published papers.

Figure 4. The presence or absence of magnetic fields to model the system in published papers.

Figure 5. The contribution of each publisher of the published papers related to nanofluids in a porous material.

Table 1

The porosity of various materials [2].

Material	Porosity $\varnothing$
Agar–agar	0.57–0.66
Black slate powder	0.12–0.34
Brick	0.45
Catalyst (Fischer–Tropsch, granules only)
Cigarette	0.17–0.49
Cigarette filters	0.02–0.12
Coal Concrete (ordinary mixes)	~0.10
Concrete (bituminous)
Copper powder (hot-compacted)	0.09–0.34
Corkboard
Fiberglass	0.88–0.93
Granular crushed rock	0.45
Hair (on mammals)	0.95–0.99
Hair felt
Leather	0.56–0.59
Hair felt limestone (dolomite)	0.04–0.10
Leather	0.37–0.50
Sand sandstone (“oil sand”)	0.08–0.38
Silica grains	0.65
Silica powder	0.37–0.49
Soil	0.43–0.54
Spherical packings (shaken well)	0.36–0.43
Wire crimps	0.68–0.76

Table 2

Research paper characteristics related to the free convection heat transfer.

Ref	Geometry Description	Nanofluid	Methodology	Results	Decision Variables
[15]	Inclined, porous, semi-annulus enclosure	Magnetic Fe3O4 -water	Free convection,Buongiorno and Darcy models,FVM,SIMPLE	- Adding nanoparticle volume fraction → Nu increases- Increase in porosity number → Nu increases	10 ≤ Ra ≤ 1000Porosity number = 0.4, 0.70 ≤ $\phi$ ≤ 0.040 ≤ inclination angle of cavity ≤ 90
[16]	Square enclosure and convection around a circular cylinder, different geometries of cylinders	Ag-water	Free convection,Darcy–Brinkman model	- Porous layer thickness increases (20% to 80%) → free convection performance decreases (up to 50%)	103 < Ra < 10610−5 < Da < 10−10% < thickness ofporous layer < 100%1 < thermal conductivity ratio 0 < $\phi$ < 0.1
[17]	Square enclosure	MWCNT–Fe3O4/water	Free convectiveMHD,MRT, Lattice–Boltzmann	- Increase in Ra → increase in heat transfer rate- Increase in Ha → decrease in Ra- Increase in Nu (+4.9%)	10−2 < Da < 10−1;103 < Ra < 105;0.4 < porosity < 0.9;0 < $\phi$ < 0.003;0 < Ha < 50;
[18]	Inclined square enclosure and exothermic chemical reaction administered by Arrhenius kinetics	Tilted nanofluid	Free convectiveBuongiorno nanofluid model,FEM	- Re increases → Nu decreases	Dissemination of streamlines;isotherms; iso-concentrations;and average Nusselt number
[19]	Square cavity and linearly heated left wall with composite nanofluid–porous layers	Cu-water	Free convection,Galerkin finite element method,Darcy–Brinkmann model	- Increase in Ra → intense streamlines	$\phi$ = 0.1;10−7 ≤ Da ≤ 1;103 ≤ Ra ≤ 107
[20]	Inverse T-shaped cavity	MWCNT–Fe3O4/water	Free convectionMHD,extended Darcy–Brinkman–Forchheimer model	- Lower inclination angle → higher Nu- Lower values of ratio of dimensionless convection coefficient and the magnetic field viscosity parameter → significant heat transfer enhancement	0 ≤ magnetic field viscosity parameter ≤ 1;0.7 ≤ porosity ratio ≤ 1.4;0 ≤ magnetic field inclination angle ≤ $\pi ;$0 ≤ ratio of dimensionless convection coefficient ≤ 10;Ha = 20;Ra = 105
[21]	Square cavity and two semicircular heat sources in the wall	MWCNT–Fe3O4/water	Free convection, FEM	- Ra = 1 × 104 → Nu increases with magnetic number	100 < Magnetic number < 5000;0.2 < Strength ratio of magnetic sources < 5;0 < Ha < 50;0.1 < porosity coefficient < 9
[22]	[Image omitted. Please see PDF.]	Nano-Encapsulated Phase ChangeMaterials (NEPCM)	Free convection.local thermal non-equilibrium(LTNE)	- Increase in thermal conductivity of porous medium → and increase in heat transfer	0 ≤ $\phi$ ≤ 0.05
[23]	Transient natural convection and a square cavity, considering nanoparticle sedimentation	Al2O3/water	Free convection	- Nu decreased- Reduction in convection heat transfer	104 < Ra < 107;10−5 < Da < 10−2
[24]	Square cavity	Ag–MgO/water	Free convection,LTNE,Darcy model,Galerkin FEM	- Increase in Ra→ increase in the vortex’s strength- Increase in heat transfer (5.85 times)	10 ≤ Ra ≤ 1000;0.1 ≤ ε≤0.9;0 ≤ φ ≤ 0.02;1 ≤ H ≤ 1000
[25]	Inclined enclosure with wavy walls and partially layered porous medium	Cu-Al2O3 water	Free convection,Galerkin FEM,Darcy–Brinkman model	- Increase in heat transfer	0 < inclination angle < 90;104 ≤ Ra ≤ 107;10−2 ≤ Da ≤ 10−5;0.2 ≤ porous layer width ≤ 0.8;1 ≤ number of undulations ≤ 4;0 ≤ $\phi$ ≤ 0.2
[26]	Eccentricity heat source and porous annulus	Cu-water	Free convection	- Increase in heat transfer	0 ≤ ϕ ≤ 0.04;103 ≤ Ra ≤ 106;10−4 ≤ Da ≤ 10−1;
[27]	Transient natural convection and non-Darcy porous cavity with an inner solid body	Al2O3 -Water	Free convection,Buongiorno model,Brinkman–Forchheimer extended Darcy formulation.FDM	- Higher Da→ uniform nanoparticle distribution- Increasing porosity → uniform nanoparticle distribution- Maximum Nu enhancement is approximately 30%	The porosity of the porous medium;Darcy number;The nanoparticles’ average volume fraction
[28]	Inner corrugated cylinders inside wavy enclosure and porous–nanofluid layers	Ag nanofluid	Free convection	- Increase in Ra and Da → increase in fluid flow strength and shear layer thickness- Increase in porous layer thickness→ decrease in heat transfer	106 ≥ Ra ≥ 103;0.1 ≥ Da ≥ 0.00001;0.2 ≥ vertical location (H) ≥ −0.2;6 ≥ number of sinusoidal inners;cylinders (N) ≥ 3
[29]	Inverse T-shaped cavity and trapezoidal heat source in the wall	Fe3O4-water	Free convection,magnetic fielddependent (MFD),FEM	- Local and average Nu increased	Darcy,Hartmann, andRayleigh numbers; inclination angle;cavity aspect ratio
[30]	Spherical electronic device	Cu-water	Free convection,SIMPLE algorithm	- Heat transfer increases- Average Nu increases	6.5 × 106 < Ra < 1.32 × 109;0 < $\phi$ < 10%;0 < thermal conductivity of the porous material’s matrix < 40
[31]	Tilted hemispherical enclosure	Water-ZnO	Free convectionexperiment	Increase in heat transfer	0 < inclination angle < 90;0 < $\phi$ < 8.22%
[32]	Wavy-walled porous cavity and inner solid cylinder	Al2O3/water	Free convection,FEM,Forchheimer–Brinkman extended Darcy model,Boussinesq approximation	- Higher values of Da → heat transfer enhancement	0 ≤ $\phi$ ≤ 0.04;10−6 < Da < 10−2;0.2 ≤ ε ≤ 0.8
[33]	Partitioned porous cavity for application in solar power plants	MWCNT–Fe3O4/water	Free convection,CFD method,volume averaging the microscopic equations	- Increase in Da, Ra → Nuave increases	103 < Ra < 106;0.5< porosity coefficient ratio < 1.8;0 < $\phi$ < 0.003;0.1 < Ri < 20;0.01 < Da < 100;Thermal conductivity ratio = 0.2, 0.4, 1, 5
[34]	Square cavity and inner sinusoidal vertical interface	Ag/water	Free convection,Galerkin FEM	- Increase in Da, Pr → Nuave increases	0.6 < power law index < 1.4;10−5 < Da < 10−1;0 < $\phi$ < 0.2;1 < undulation number (N) < 4;0.015 < Pr < 13.4;Ra = 105
[35]	Hot rectangular cylinder and cold circular cylinder	copper–water	Free convection,Brinkman-extendedDarcy model,Brinkman correlation	- Heat transfer enhanced	Rayleigh number;Hartmann number;Darcy number;magnetic field inclinationangle;nanoparticles volume fraction;nanoparticles shape factor;nanoparticles material;nanofluid thermal conductivity;dynamic viscosity models; nanofluid electrical conductivity correlation on streamlines;isotherms;local and average Nusselt numbers
[36]	Partially heated enclosure	Al2O3/water	Free convection,FEM,Brinkman equation	- Heat transfer rate augmented- Ra, Da increases → average velocity	103 < Ra < 106;0 < $\phi$ < 5%;0 < Ha < 100;0.001 < Da < 1
[37]	I -shaped cavity	Cu–water	Free convection,MHD,FDM	- Ha increases → Nu decreases- Ra increases → Nu increases- Maximum Nu occurs at B = 0.2- Minimum Nu occurs at B = 0.8	Ha;nanofluidvolume fraction;heat source size;location and angle of magnetic field on heat transfer;entropy generation;thermal performance
[38]	Porous enclosure	Cu, Al2O3 and TiO2/water	Free convection,MHD	- Increase in magnetic field intensity→ heat transfer deterioration- Enlarging nanoparticles, denser nanoparticles→ heat transfer deterioration	0 ≤ Ha ≤ 50;Nanoparticle volume fraction;Nanoparticle diameter
[39]	Inclined cavity	Al2O3-water	Free convectionEntropy generation	- Increase in chamber angle → increase in heat transfer- Adding nanoparticle volume fraction → increase in heat transfer	Rayleigh numberHartmann number;magnetic field angle changes;chamber angle changes;entropy parameter;radiationparameter;volume percent of nanoparticles
[40]	Cubical electronic component and hemispherical cavity	Water-ZnO	Free convection,control volume method	- Inclination increases → Nuava decreases- Nanofluid concentration increases → heat transfer increases	0 < volume fraction < 10%;Nuave
[41]	Inverted T-shape	MWCNT–Fe3O4/water	Free convection, thermal transmission	- Ha increases → Nuave decreases	Heat transfer performance;flow structures
[42]	Inverse T-shaped cavity and trapezoidal heat source in wall with wavy Wall	Magnetic Al2O3/water	Free convection,FEM,Koo–Kleinstreuer–Li (KKL) correlations	- Increase in Ra, decrease in Ha → increase in flow intensity-	Heatgeneration parameter;the shape factor of nanoparticles;Hartmann number;nanoparticle concentration;displacement of the trapezoidal heater wall; Rayleigh number; the amplitude of wavy wall

Table 3

Research paper characteristics related to the forced convection heat transfer.

Ref	Geometry Description	Nanofluid	Methodology	Results	Decision Variables
[45]	Moving surface	Water- based nanoparticles: copper (Cu), alumina (Al2O3), and titania (TiO2)	Forced convection, MHD	- Smallest temperature distribution: Cu- Largest velocity profile: Cu	Skin friction coefficient; Local Nu number
[46]	Channel, staggered, and in-line arrangements of square pillars	Al2O3-water	Forced convection,First and Second laws of thermodynamics,FVM	- Nu increases and decreases for the Re and nanofluidvolume fraction.- The Al2O3 nanoparticlesparticipation in the base fluid decreases the entropy generation.- The entropy generation and the Be decrease and increase with the nanofluid particle volume fraction.	Porosity: 0.84, 0.75, 0.91;Re = 10, 200, 300;Nanofluid V.F. = 4%
[47]	Multi-layered, U-shaped vented cavity and wall corrugation effects	CNT-water	Forced convection,FEM	- Heat transfer enhancement	100 < Re < 1000;0 < Ha < 50;10−4 < Da < 5 × 10−2
[48]	Lid-driven cavity and hot sphere obstacle	Al2O3-water	Forced convection,Lattice–Boltzmann method	- Rate of heat transfer enhances with the rise of permeability of porous media and velocityof lid wall. This is due to an enhanced temperature gradient with the increase ofDa and Re.	0.001 < Da < 100;0 < Ha < 40;30 < Re < 180
[49]	U-bend pipe	Al2O3-CuO-water	Forced convection,FEM,Darcy–Brinkman–Forchheimer equation	- Decrease in Da → increase in Nuave → Increase in pressure drop	10−4 < Da < 10−1
[50]	Cylinder	Al2O3–CuO–water	Forced convection,MHD,FVM	- Decrease in Da → increase in Nu → increase in Ha- Increase in Da and Ha → decrease in pressure drop- Adding metal nanoparticles → increase in Da and Ha	Nu number;0.0001 < Da < 0.1;0 < Ha < 40;Magnetic field orientation
[51]	Annulus with porous ribs	Al2O3-water	Forced convective	- Increase in porous ribs → increase in pressure drop	Nu
[52]	Horizontal plate	Water-based Cu/Alumina/Titania	Forced convection	- Heat transfer rate is higher for Cu than others- Increase in nanoparticles V.F. → increase in heat transfer rate	Nanoparticle volume fraction;porosity of porous media;thermal conductivity of porous media;effect of nanoparticle type on its heat transfer

Table 4

Research paper characteristics related to the mixed convection heat transfer.

Ref	Geometry Description	Nanofluid	Methodology	Results	Decision Variables
[53]	Inclined flat plate	Water-Cu	Mixed convection,MHD,DQM	Nu reduced	Ra = 105;Ha = 25
[54]	Lid-driven enclosure and two adherent porous blocks	Alumina/water	Mixed convection	-Ri < 1 → heat transfer enhancement- Ri ≥ 1 → reduction in heat transfer	0.01 ≤ Ri ≤ 10;0 ≤ $\phi$ ≤ 0.04
[55]	Rotating circular cylinder and trapezoidal enclosure	Cu-water	Mixed convection,MHD	- Decrease in stream function values → vertical magnetic field- Increase in Ha → increase in Nuave- Increase in Ha, thermal conductivity rate, cylinder radius, Da → increase in Nuave- Decrease in Ri → increase in Nuave	0 < Ha < 100;1 < Thermal conductivity ratio < 10;−5 < angular rotational velocity < 5;0.01 < Ri < 100;0 < Inclination angle < 90; 0.2 < Cylinder radius < 0.4;10−5 < Da < 10−1;0 < nanofluid concentration < 0.1
[56]	Square cavity with inlet and outlet ports	Water-based nanofluid	Mixed convection,Brownian diffusion,thermophoresis,FDM	- Increase in Re → cooling improvement- Ra = 10 → Nu = 1.071- Ra = 100 → Nu = 3.104- Ra = 1000 → Nu = 13.839- Ra = 10000 → Nu = 49.253$-\varnothing = 0.01\to Nu=31.6043$- $\varnothing = 0.02\to Nu=31.2538$- $\varnothing = 0.03\to Nu=30.829$	104 < Ra < 106;Pr = 6.82;10−5 < Da < 10−6;50 < Re < 300;ε = 0.5;Le = 1000
[57]	H-shaped cavity with cooler and heater cylinders	Cu-water	Mixed convection,Boussinesq approximation	- Increase in AR → decrease in heat transfer rateincrease in Da, decrease in Ri → increase in heat transfer rate	10−4 ≤ Da ≤ 10−2;1 ≤ Ri ≤ 100;1.4 ≤ AR ≤ 1.6;
[58]	Trapezoidal chamber	Cu-Al2O3/water	Mixed convection,FDM	- Increase in Re → increase in energy transport and convective circulation- Increase in Da → heat transfer enhancement	Reynolds number;Darcy number;nanoparticle volume fraction
[59]	Inclined cavity	Cu-water	Mixed convection,Darcy–Brinkman–Forchheimer model,SIMPLE algorithm	- Heat transfer rate increases with increasing Da.
[60]	Lid-driven square cavity	Al2O3/water	Mixed convection	- Decrease in Ri → increase in momentum- Ri = 100 → decrease in Darcy effects- Changing nanoparticles volume fraction and Da → significant changes in streamlined pattern- Higher Ri → more buoyancy effects- Increase in Da and Ri → less fluid resistance and more momentum penetration- Increase in Da → decrease in temperature, more uniformity in heat transfer	Ri = 0.01, 10 and 100;10−4 ≤ Da ≤ 10−2;0 ≤ $\phi$ ≤ 0.04
[61]	Stretching surface	---	Mixed convection, MHD	- For m < 1→ increase in velocity results in an increase in thermophoresis- For m > 1→ increase in velocity results in a decrease in thermophoresis.	Effects of buoyancy parameter;magnetic parameter;Brownian motion;thermophoresisparameter, etc., on velocity, temperature, and nanoparticle volume fraction
[62]	Square cavity and two rotating cylinders	Al2O3/water	Mixed convection	- Heat transfer enhancement (+ 20.4%)
[63]	Triangular shape, partitioned, lid-driven square cavity involving a porous compound	Ag–MgO/water	Mixed convectionMHD	- Nu enhancement (14.7%)	< Ri < 100;0 < Ha < 60;10−4 < Da < 5 × 10−2;0 < $\phi$ < 0.01
[64]	Vertical surface	Cu-water	Mixed convection,Laplacetransform techniqueCrank Nicolson method	- Increase in magnetic field strength → and decrease in fluid velocity- Porosity increases → fluid velocity decreases	Magneticparameter;porosity parameter;thermal and solute Grashof number;nanoparticle volumefraction parameter;time;Schmidt number;chemical reaction parameter;Prandtl number
[65]	Inclined cavity and porous layer	Cu-water	Mixed convection,incompressible smoothed particle hydrodynamics (ISPH)	- Ri increases → Nuave decreases-$\phi$ increases → overall heat transfer increases	0.001 < Ri < 100;10−5 < Da < 10−2;0 < $\phi$ < 0.05
[66]	[Image omitted. Please see PDF.]	CuO–Water	Mixed convection,entropy generation,Buongiorno’s two-phase model	- Increase in volume concentration → increase in Nuave- Maximum enhancement in cooling performance was 17.75%	- Volume concentration; development of a new predictive correlation
[67]	Gamma-shaped cavity	CuO–Water	Mixed convection,Entropy generation,FVM	- Increase in the Nusselt number with the volumefraction is more pronounced for the smallest heat source, a heat source placed at the lowest height from the bottomside, the lowest volumetric heat generation, the lowest imposed magnetic field, the lowest Darcy number, and for aporous media with the lowest solid to fluid thermal conductivity ratio. Increasing the nanoparticle volume fraction has a higher impact on the production of entropy than the enhancement in the heat transfer rate.	- Hartmann number; nanoparticle volume fraction; the length and location of a heat source
[68]	Rotating triangle chamber	Graphene Oxide generalized hybrid	Mixed convectionutilizing bvp4c solver	The velocity upsurges due to the dimensionless radius of the slender body parameter in case of the assisting flow.	0.025 ≤ $\phi$ ≤ 0.035

Table 5

Nu changes with nanofluid concentration and size.

Ref	Nu	Volume Fraction and Size
[15]	Nu increases 0.35%Heat transfer increases 0.48%	$0\le \phi$ ≤ 0.040.384 ≤ dp ≤ 2.5
[17]	Nu increases 4.9%Heat transfer increases 6.72%	$0<\phi$ < 0.0031.46 ≤ dp ≤ 2.8
[20]	Nu increases 3.6%Heat transfer increases 5.94%	$0<\phi$ < 0. 30.65 ≤ dp ≤ 1.41
[21]	Nu increases 6.01%Heat transfer increases 9.24%	$\phi$ < 0. 34.5 ≤ dp ≤ 7.6
[23]	Nu decreased 6.64%Heat transfer decreased 8.21%	$\phi$ < 0.050.51 ≤ dp ≤ 2.3
[26]	Nu increases 44.44%	0 ≤ ϕ ≤ 0.040.5 ≤ dp ≤ 20
[29]	Nu decreased 0.69%Heat transfer decreased 0.8%	0 ≤ ϕ ≤ 0.15 ≤ dp ≤ 15
[32]	Nu increases 5.69%Heat transfer increases 9.28%	$\phi$ ≤ 0.040.385 ≤ dp ≤ 33
[35]	Nu increases 24.98%	$\phi$ < 0.0038 ≤ dp ≤ 29
[50]	Nu increases 0.5%Heat transfer increases 2.3%	$\phi$ < 0.0050.52 ≤ dp ≤ 7.5
[53]	Nu decreased 30%	$\phi$ < 0.152 ≤ dp ≤ 24
[56]	Nu increases 66.6%	0 ≤ $\phi$ ≤ 0.044 ≤ dp ≤ 14
[57]	Nu decreased 12.28%	0 ≤ $\phi$ ≤ 0.010.89 ≤ dp ≤ 1.3
[60]	Nu decreased 26.08%	$\phi$ ≤ 0.040.384 ≤ dp ≤ 47
[62]	Nu increases 5.97%Heat transfer increases 20.4%	$\phi$ ≤ 23 ≤ dp ≤ 25
[63]	Nu increases 14.7%	$\phi$ ≤ 0.017 ≤ dp ≤ 44
[65]	Nu decreased 68.75%	$\phi$ ≤ 0.051.94 ≤ dp ≤ 6.29

Table 6

The share of each nanoparticle in published studies.

Nanoparticle	Share (%)
Al2O3-Water	33.9
Cu-Water	32.14
Fe3O4-Water	12.5
Ag-Water	5.35
Ag-MgO	3.75
Other nanoparticles	12.5

Table 7

The share of each of the studied geometries in published papers.

Geometry	Share (%)
Square	54
Hole	18
Circular	12
Other shapes	16

References

1. Lapwood, E.R. Convection of a fluid in a porous medium. Mathematical Proceedings of the Cambridge Philosophical Society; Cambridge University Press: Cambridge, UK, 1948; Volume 44, pp. 508-521. [DOI: https://dx.doi.org/10.1017/S030500410002452X]

2. Bejan, A.; Dincer, I.; Lorente, S.; Miguel, A.F.; Reis, A.H. Porous and Complex Flow Structures in Modern Technologies; Springer: New York, NY, USA, 2004.

3. Whitaker, S. Flow in porous media I: A theoretical derivation of Darcy’s law. Transp. Porous Media; 1986; 1, pp. 3-25. [DOI: https://dx.doi.org/10.1007/BF01036523]

4. Murshed, S.M.S.; Leong, K.C.; Yang, C. Thermophysical and electrokinetic properties of nanofluids—A critical review. Appl. Therm. Eng.; 2008; 28, pp. 2109-2125. [DOI: https://dx.doi.org/10.1016/j.applthermaleng.2008.01.005]

5. Meng, X.; Yang, D. Critical Review of Stabilized Nanoparticle Transport in Porous Media. J. Energy Resour. Technol.; 2019; 141, content [DOI: https://dx.doi.org/10.1115/1.4041929]

6. Boccardo, G.; Tosco, T.; Fujisaki, A.; Messina, F.; Raoof, A.; Aguilera, D.R.; Crevacore, E.; Marchisio, D.L.; Sethi, R. A review of transport of nanoparticles in porous media: From pore-to macroscale using computational methods. Nanomater. Detect. Remov. Wastewater Pollut.; 2020; pp. 351–381.-381. [DOI: https://dx.doi.org/10.1016/B978-0-12-818489-9.00013-X]

7. Ling, X.; Yan, Z.; Liu, Y.; Lu, G. Transport of nanoparticles in porous media and its effects on the co-existing pollutants. Environ. Pollut.; 2021; 283, 117098. [DOI: https://dx.doi.org/10.1016/j.envpol.2021.117098]

8. Hyun, J.M.; Choi, B.S. Transient natural convection in a parallelogram-shaped enclosure. Int. J. Heat Fluid Flow; 1990; 11, pp. 129-134. [DOI: https://dx.doi.org/10.1016/0142-727X(90)90006-W]

9. Takabi, B.; Shokouhmand, H. Effects of Al₂O₃-Cu/water hybrid nanofluid on heat transfer and flow characteristics in turbulent regime. Int. J. Mod. Phys. C; 2015; 26, 1550047. [DOI: https://dx.doi.org/10.1142/S0129183115500473]

10. Nield, D.A.; Bejan, A. Convection in Porous Media; Springer: New York, NY, USA, 2013.

11. Buongiorno, J. Convective transport in nanofluids. J. Heat Transfer.; 2006; 128, pp. 240-250. [DOI: https://dx.doi.org/10.1115/1.2150834]

12. Jou, R.Y.; Tzeng, S.C. Numerical research of natural convective heat transfer enhancement filled with nanofluids in rectangular enclosures. Int. Commun. Heat Mass Transf.; 2006; 33, pp. 727-736. [DOI: https://dx.doi.org/10.1016/j.icheatmasstransfer.2006.02.016]

13. Agarwal, D.K.; Vaidyanathan, A.; Sunil Kumar, S. Synthesis and characterization of kerosene-alumina nanofluids. Appl. Therm. Eng.; 2013; 60, pp. 275-284. [DOI: https://dx.doi.org/10.1016/j.applthermaleng.2013.06.049]

14. Motlagh, S.Y.; Golab, E.; Sadr, A.N. Two-phase modeling of the free convection of nanofluid inside the inclined porous semi-annulus enclosure. Int. J. Mech. Sci.; 2019; 164, 105183. [DOI: https://dx.doi.org/10.1016/j.ijmecsci.2019.105183]

15. Hussain, S.H.; Rahomey, M.S. Comparison of Natural Convection Around a Circular Cylinder with Different Geometries of Cylinders Inside a Square Enclosure Filled with Ag Nanofluid Superposed Porous Nanofluid Layers. J. Heat Transfer.; 2019; 141, 022501. [DOI: https://dx.doi.org/10.1115/1.4039642]

16. Sajjadi, H.; Amiri Delouei, A.; Izadi, M.; Mohebbi, R. Investigation of MHD natural convection in a porous media by double MRT lattice Boltzmann method utilizing MWCNT–Fe₃O₄/water hybrid nanofluid. Int. J. Heat Mass Transf.; 2019; 132, pp. 1087-1104. [DOI: https://dx.doi.org/10.1016/j.ijheatmasstransfer.2018.12.060]

17. Rahman, M.M.; Pop, I.; Saghir, M.Z. Steady free convection flow within a titled nanofluid saturated porous cavity in the presence of a sloping magnetic field energized by an exothermic chemical reaction administered by Arrhenius kinetics. Int. J. Heat Mass Transf.; 2019; 129, pp. 198-211. [DOI: https://dx.doi.org/10.1016/j.ijheatmasstransfer.2018.09.105]

18. Al-Srayyih, B.M.; Gao, S.; Hussain, S.H. Effects of linearly heated left wall on natural convection within a superposed cavity filled with composite nanofluid-porous layers. Adv. Powder Technol.; 2019; 30, pp. 55-72. [DOI: https://dx.doi.org/10.1016/j.apt.2018.10.007]

19. Mehryan, A.M.; Sheremet, M.A.; Soltani, M.; Izadi, M. Natural convection of magnetic hybrid nanofluid inside a double-porous medium using two-equation energy model. J. Mol. Liq.; 2019; 277, pp. 959-970. [DOI: https://dx.doi.org/10.1016/j.molliq.2018.12.147]

20. Izadi, M.; Mohebbi, R.; Delouei, A.A.; Sajjadi, H. Natural convection of a magnetizable hybrid nanofluid inside a porous enclosure subjected to two variable magnetic fields. Int. J. Mech. Sci.; 2019; 151, pp. 154-169. [DOI: https://dx.doi.org/10.1016/j.ijmecsci.2018.11.019]

21. Ghalambaz, M.; Hashem Zadeh, S.M.; Mehryan, S.A.M.; Haghparast, A.; Zargartalebi, H. Free convection of a suspension containing nano-encapsulated phase change material in a porous cavity; local thermal non-equilibrium model. Heliyon; 2020; 6, e03823. [DOI: https://dx.doi.org/10.1016/j.heliyon.2020.e03823]

22. Baghsaz, S.; Rezanejad, S.; Moghimi, M. Numerical investigation of transient natural convection and entropy generation analysis in a porous cavity filled with nanofluid considering nanoparticles sedimentation. J. Mol. Liq.; 2019; 279, pp. 327-341. [DOI: https://dx.doi.org/10.1016/j.molliq.2019.01.117]

23. Mehryan, S.A.M.; Ghalambaz, M.; Chamkha, A.J.; Izadi, M. Numerical study on natural convection of Ag–MgO hybrid/water nanofluid inside a porous enclosure: A local thermal non-equilibrium model. Powder Technol.; 2020; 367, pp. 443-455. [DOI: https://dx.doi.org/10.1016/j.powtec.2020.04.005]

24. Kadhim, H.T.; Jabbar, F.A.; Rona, A. Cu-Al₂O₃ hybrid nanofluid natural convection in an inclined enclosure with wavy walls partially layered by porous medium. Int. J. Mech. Sci.; 2020; 186, 105889. [DOI: https://dx.doi.org/10.1016/j.ijmecsci.2020.105889]

25. Gholamalipour, P.; Siavashi, M.; Doranehgard, M.H. Eccentricity effects of heat source inside a porous annulus on the natural convection heat transfer and entropy generation of Cu-water nanofluid. Int. Commun. Heat Mass Transf.; 2019; 109, 104367. [DOI: https://dx.doi.org/10.1016/j.icheatmasstransfer.2019.104367]

26. Alsabery, A.I.; Ismael, M.A.; Chamkha, A.J.; Hashim, I. Effect of nonhomogeneous nanofluid model on transient natural convection in a non-Darcy porous cavity containing an inner solid body. Int. Commun. Heat Mass Transf.; 2020; 110, 104442. [DOI: https://dx.doi.org/10.1016/j.icheatmasstransfer.2019.104442]

27. Abdulkadhim, A.; Hamzah, H.K.; Ali, F.H.; Abed, A.M.; Abed, I.M. Natural convection among inner corrugated cylinders inside a wavy enclosure filled with nanofluid superposed in porous–nanofluid layers. Int. Commun. Heat Mass Transf.; 2019; 109, 104350. [DOI: https://dx.doi.org/10.1016/j.icheatmasstransfer.2019.104350]

28. Molana, M.; Dogonchi, A.S.; Armaghani, T.; Chamkha, A.J.; Ganji, D.D.; Tlili, I. Investigation of Hydrothermal Behavior of Fe₃O₄-H₂O Nanofluid Natural Convection in a Novel Shape of Porous Cavity Subjected to Magnetic Field Dependent (MFD) Viscosity. J. Energy Storage; 2020; 30, 101395. [DOI: https://dx.doi.org/10.1016/j.est.2020.101395]

29. Baïri, A. Using nanofluid saturated porous media to enhance free convective heat transfer around a spherical electronic device. Chin. J. Phys.; 2020; 70, pp. 106-116. [DOI: https://dx.doi.org/10.1016/j.cjph.2020.03.023]

30. Baïri, A. Experimental study on enhancement of free convective heat transfer in a tilted hemispherical enclosure through Water-ZnO nanofluid saturated porous materials. Appl. Therm. Eng.; 2019; 148, pp. 992-998. [DOI: https://dx.doi.org/10.1016/j.applthermaleng.2018.11.115]

31. Alsabery, A.I.; Mohebbi, R.; Chamkha, A.J.; Hashim, I. Effect of the local thermal non-equilibrium model on natural convection in a nanofluid-filled wavy-walled porous cavity containing an inner solid cylinder. Chem. Eng. Sci.; 2019; 201, pp. 247-263. [DOI: https://dx.doi.org/10.1016/j.ces.2019.03.006]

32. Mohebbi, R.; Mehryan, S.A.M.; Izadi, M.; Mahian, O. Natural convection of hybrid nanofluids inside a partitioned porous cavity for application in solar power plants. J. Therm. Anal. Calorim.; 2019; 137, pp. 1719-1733. [DOI: https://dx.doi.org/10.1007/s10973-019-08019-9]

33. Al-Amir, Q.R.; Ahmed, S.Y.; Hamzah, H.K.; Ali, F.H. Effects of Prandtl Number on Natural Convection in a Cavity Filled with Silver/Water Nanofluid-Saturated Porous Medium and Non-Newtonian Fluid Layers Separated by Sinusoidal Vertical Interface. Arab. J. Sci. Eng.; 2019; 44, pp. 10339-10354. [DOI: https://dx.doi.org/10.1007/s13369-019-04115-y]

34. Dogonchi, S.; Sheremet, M.A.; Ganji, D.D.; Pop, I. Free convection of copper–water nanofluid in a porous gap between the hot rectangular cylinder and cold circular cylinder under the effect ofthe inclined magnetic field. J. Therm. Anal. Calorim.; 2019; 135, pp. 1171-1184. [DOI: https://dx.doi.org/10.1007/s10973-018-7396-3]

35. Akhter, R.; Ali, M.M.; Alim, M.A. Hydromagnetic natural convection heat transfer in a partially heated enclosure filled with porous medium saturated by a nanofluid. Int. J. Appl. Comput. Math.; 2019; 5, pp. 1-27. [DOI: https://dx.doi.org/10.1007/s40819-019-0638-7]

36. Armaghani, T.; Chamkha, A.; Rashad, A.M.; Mansour, M.A. Inclined magneto: Convection, internal heat, and entropy generation of nanofluid in an I-shaped cavity saturated with porous media. J. Therm. Anal. Calorim.; 2020; 142, pp. 2273-2285. [DOI: https://dx.doi.org/10.1007/s10973-020-09449-6]

37. Benos, L.T.; Polychronopoulos, N.D.; Mahabaleshwar, U.S.; Lorenzini, G.; Sarris, I.E. Thermal and flow investigation of MHD natural convection in a nanofluid-saturated porous enclosure: An asymptotic analysis. J. Therm. Anal. Calorim.; 2019; 143, pp. 751-765. [DOI: https://dx.doi.org/10.1007/s10973-019-09165-w]

38. Hajatzadeh Pordanjani, A.; Aghakhani, S.; Karimipour, A.; Afrand, M.; Goodarzi, M. Investigation of free convection heat transfer and entropy generation of nanofluid flow inside a cavity affected by the magnetic field and thermal radiation. J. Therm. Anal. Calorim.; 2019; 137, pp. 997-1019. [DOI: https://dx.doi.org/10.1007/s10973-018-7982-4]

39. Baïri, A. New correlations for free convection with water-ZnO nanofluid saturated porous medium around a cubical electronic component in the hemispherical cavity. Heat Transf. Eng.; 2020; 41, pp. 1275-1287. [DOI: https://dx.doi.org/10.1080/01457632.2019.1637142]

40. Izadi, M.; Oztop, H.F.; Sheremet, M.A.; Mehryan, S.A.M.; Abu-Hamdeh, N. Coupled FHD–MHD free convection of a hybrid nanoliquid in an inversed T-shaped enclosure occupied by partitioned porous media. Numer. Heat Transf. Part A Appl.; 2019; 76, pp. 479-498. [DOI: https://dx.doi.org/10.1080/10407782.2019.1637626]

41. Sadeghi, M.S.; Tayebi, T.; Dogonchi, A.S.; Armaghani, T.; Talebizadehsardari, P. Analysis of hydrothermal characteristics of magnetic Al₂O₃-H₂O nanofluid within a novel wavy enclosure during natural convection process considering internal heat generation. Math. Methods Appl. Sci.; 2020; mma.6520. [DOI: https://dx.doi.org/10.1002/mma.6520]

42. Sheikholeslami, M.; Hayat, T.; Muhammad, T.; Alsaedi, A. MHD forced convection flow of a nanofluid in a porous cavity with the hot elliptic obstaclethrough Lattice Boltzmann method. Int. J. Mech. Sci.; 2018; 135, pp. 532-540. [DOI: https://dx.doi.org/10.1016/j.ijmecsci.2017.12.005]

43. Sheikholeslami, M.; Hayat, T.; Alsaedi, A. Numerical simulation of nanofluid forced convection heat transfer improvement in the existence of magnetic field using lattice Boltzmann method. Int. J. Heat Mass Transf.; 2017; 108, pp. 1870-1883. [DOI: https://dx.doi.org/10.1016/j.ijheatmasstransfer.2017.01.044]

44. Ferdows, M.; Alzahrani, F. Dual solutions of nanofluid forced convective flow with heat transfer and porous media past a moving surface. Phys. A Stat. Mech. Its Appl.; 2020; 551, 124075. [DOI: https://dx.doi.org/10.1016/j.physa.2019.124075]

45. Torabi, M.; Torabi, M.; Ghiaasiaan, S.M.; Peterson, G.P. The effect of Al₂O₃-water nanofluid on the heat transfer and entropy generation of laminar forced convection through isotropic porous media. Int. J. Heat Mass Transf.; 2017; 111, pp. 804-816. [DOI: https://dx.doi.org/10.1016/j.ijheatmasstransfer.2017.04.041]

46. Selimefendigil, F.; Öztop, H.F. Magnetohydrodynamics forced convection of nanofluid in multi-layered U-shaped vented cavity with a porous region considering wall corrugation effects. Int. Commun. Heat Mass Transf.; 2020; 113, 104551. [DOI: https://dx.doi.org/10.1016/j.icheatmasstransfer.2020.104551]

47. Sheikholeslami, M. Influence of magnetic field on Al₂O₃-H₂O nanofluid forced convection heat transfer in a porous lid-driven cavity with hot sphere obstacle through LBM. J. Mol. Liq.; 2018; 263, pp. 472-488. [DOI: https://dx.doi.org/10.1016/j.molliq.2018.04.111]

48. Moghadasi, H.; Aminian, E.; Saffari, H.; Mahjoorghani, M.; Emamifar, A. Numerical analysis on laminar forced convection improvement of hybrid nanofluid within a U-bend pipe in porous media. Int. J. Mech. Sci.; 2020; 179, 105659. [DOI: https://dx.doi.org/10.1016/j.ijmecsci.2020.105659]

49. Aminian, E.; Moghadasi, H.; Saffari, H. Magnetic field effects on forced convection flow of a hybrid nanofluid in a cylinder filled with porous media: A numerical study. J. Therm. Anal. Calorim.; 2020; 141, pp. 2019-2031. [DOI: https://dx.doi.org/10.1007/s10973-020-09257-y]

50. Siavashi, M.; Talesh Bahrami, H.R.; Aminian, E.; Saffari, H. Numerical analysis on forced convection enhancement in an annulus using porous ribs and nanoparticle addition to base fluid. J. Cent. South Univ.; 2019; 26, pp. 1089-1098. [DOI: https://dx.doi.org/10.1007/s11771-019-4073-z]

51. Tlili, I.; Rabeti, M.; Safdari Shadloo, M.; Abdelmalek, Z. Forced convection heat transfer of nanofluids from a horizontal plate with convective boundary condition and a line heat source embedded in porous media. J. Therm. Anal. Calorim.; 2020; 141, pp. 2081-2094. [DOI: https://dx.doi.org/10.1007/s10973-020-09473-6]

52. Khademi, R.; Razminia, A.; Shiryaev, V.I. Conjugate-mixed convection of nanofluid flow over an inclined flat plate in porous media. Appl. Math. Comput.; 2020; 366, 124761. [DOI: https://dx.doi.org/10.1016/j.amc.2019.124761]

53. Bondarenko, D.S.; Sheremet, M.A.; Oztop, H.F.; Abu-Hamdeh, N. Mixed convection heat transfer of a nanofluid in a lid-driven enclosure with two adherent porous blocks. J. Therm. Anal. Calorim.; 2019; 135, pp. 1095-1105. [DOI: https://dx.doi.org/10.1007/s10973-018-7455-9]

54. Ali, F.H.; Hamzah, H.K.; Hussein, A.K.; Jabbar, M.Y.; Talebizadehsardari, P. MHD mixed convection due to a rotating circular cylinder in a trapezoidal enclosure filled with a nanofluid saturated with a porous media. Int. J. Mech. Sci.; 2020; 181, 105688. [DOI: https://dx.doi.org/10.1016/j.ijmecsci.2020.105688]

55. Sheremet, M.A.; Roşca, N.C.; Roşca, A.V.; Pop, I. Mixed convection heat transfer in a square porous cavity filled with a nanofluid with suction/injection effect. Comput. Math. Appl.; 2018; 76, pp. 2665-2677. [DOI: https://dx.doi.org/10.1016/j.camwa.2018.08.069]

56. Li, Z.; Barnoon, P.; Toghraie, D.; Balali Dehkordi, R.; Afrand, M. Mixed convection of non-Newtonian nanofluid in an H-shaped cavity with cooler and heater cylinders filled by a porous material: Two-phase approach. Adv. Powder Technol.; 2019; 30, pp. 2666-2685. [DOI: https://dx.doi.org/10.1016/j.apt.2019.08.014]

57. Cimpean, D.S.; Sheremet, M.A.; Pop, I. Mixed convection of hybrid nanofluid in a porous trapezoidal chamber. Int. Commun. Heat Mass Transf.; 2020; 116, 104627. [DOI: https://dx.doi.org/10.1016/j.icheatmasstransfer.2020.104627]

58. Rajarathinam, M.; Nithyadevi, N.; Chamkha, A.J. Heat transfer enhancement of mixed convection in an inclined porous cavity using Cu-water nanofluid. Adv. Powder Technol.; 2018; 29, pp. 590-605. [DOI: https://dx.doi.org/10.1016/j.apt.2017.11.032]

59. Nazari, S.; Ellahi, R.; Sarafraz, M.M.; Safaei, M.R.; Asgari, A.; Akbari, O.A. Numerical study on mixed convection of a non-Newtonian nanofluid with porous media in a two lid-driven square cavity. J. Therm. Anal. Calorim.; 2020; 140, pp. 1121-1145. [DOI: https://dx.doi.org/10.1007/s10973-019-08841-1]

60. Ghosh, S.; Mukhopadhyay, S. MHD mixed convection flow of a nanofluid past a stretching surface of variable thickness and vanishing nanoparticle flux. Pramana J. Phys.; 2020; 94, pp. 1-12. [DOI: https://dx.doi.org/10.1007/s12043-020-1924-y]

61. Tahmasbi, M.; Siavashi, M.; Abbasi, H.R.; Akhlaghi, M. Mixed convection enhancement by using optimized porous media and nanofluid in a cavity with two rotating cylinders. J. Therm. Anal. Calorim.; 2020; 141, pp. 1829-1846. [DOI: https://dx.doi.org/10.1007/s10973-020-09604-z]

62. Selimefendigil, F.; Chamkha, A.J. MHD mixed convection of Ag–MgO/water nanofluid in a triangular shape partitioned lid-driven square cavity involving a porous compound. J. Therm. Anal. Calorim.; 2020; 143, pp. 1467-1484. [DOI: https://dx.doi.org/10.1007/s10973-020-09472-7]

63. Madhura, K.R.; Iyengar, S.S. Impact of Heat and Mass Transfer on Mixed Convective Flow of Nanofluid Through Porous Medium. Int. J. Appl. Comput. Math.; 2017; 3, pp. 1361-1384. [DOI: https://dx.doi.org/10.1007/s40819-017-0424-3]

64. Aly, M.; Raizah, Z.A.S. Mixed Convection in an Inclined Nanofluid Filled-Cavity Saturated with a Partially Layered Porous Medium. J. Therm. Sci. Eng. Appl.; 2019; 11, 041002. [DOI: https://dx.doi.org/10.1115/1.4042352]

65. Asadi, A.; Molana, M.; Ghasemiasl, R.; Armaghani, T.; Pop, M.I.; Saffari Pour, M. A New Thermal Conductivity Model and Two-Phase Mixed Convection of CuO–Water Nanofluids in a Novel I-Shaped Porous Cavity Heated by Oriented Triangular Hot Block. Nanomaterials; 2020; 10, 2219. [DOI: https://dx.doi.org/10.3390/nano10112219] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33171855]

66. Chamkha, J.; Mansour, M.A.; Rashad, A.M.; Kargarsharifabad, H.; Armaghani, T. Magnetohydrodynamic Mixed Convection and Entropy Analysis of Nanofluid in Gamma-Shaped Porous Cavity. J. Thermophys. Heat Transf.; 2020; 34, pp. 836-847. [DOI: https://dx.doi.org/10.2514/1.T5983]

67. Khan, U.; Zaib, A.; Sheikholeslami, M.; Wakif, A.; Baleanu, D. Mixed convective radiative flow through a slender revolution bodies containing molybdenum-disulfide graphene oxide along with generalized hybrid nanoparticles in porous media. Crystals; 2020; 10, 771. [DOI: https://dx.doi.org/10.3390/cryst10090771]

68. Choi, S.U.S. Enhancing Thermal Conductivity of Fluids with Nanoparticles; No. ANL/MSD/CP-84938; CONF-951135-29 Argonne National Lab.(ANL): Argonne, IL, USA, 1995.

69. Kasaeian, A.; Daneshazarian, R.; Mahian, O.; Kolsi, L.; Chamkha, A.J.; Wongwises, S.; Pop, I. Nanofluid flow and heat transfer in porous media: A review of the latest developments. Int. J. Heat Mass Transf.; 2017; 107, pp. 778-791.

70. Kandelousi, M.; Ameen, S.; Akhtar, M.S.; Shin, H.S. Nanofluid Flow in Porous Media; IntechOpen: London, UK, 2020.