Abstract

This study presents numerical simulations on double-diffusive flow of a nanofluid in two cavities connected with four vertical gates. Novel shape of an outer square shape mounted on a square cavity by four gates was used. Heterogeneous porous media and -water nanofluid are filled in an inner cavity. Outer rectangle shape is filled with a nanofluid only, and its left walls carry high temperature and high concentration . The right walls of a rectangle shape carry low temperature and low concentration and the other walls are adiabatic. An incompressible smoothed particle hydrodynamics (ISPH) method is applied for solving the governing equations of velocities, temperature, and concentration. Results are introduced for the effects of a buoyancy ratio , Darcy parameter , solid volume fraction , and porous levels. Main results are indicated in which the buoyancy ratio parameter adjusts the directions of double-diffusive convection flow in an outer shape and inner cavity. Adding more concentration of nanoparticles reduces the flow speed and maximum of the velocity field. Due to the presence of a porous medium layer in an inner cavity, the Darcy parameter has slight changes inside the rectangle shape.

1. Introduction

Double-diffusive convection is a form of convection flow resulting from the variations of the diffusion rates [1]. Double-diffusive has several applications such as fuel cells, diffusion of chemical pollutants in ocean, and nuclear waste storage. In addition, heat and mass transfer in porous media are surviving in numerous engineering applications including oceanography, geophysics, chemical engineering, and thermal engineering [2]. Double-diffusive convection in an annulus has been studied intensively due to its applications in dying and cleaning operations and energy in solar ponds [3].

The adjoint influences of thermal and solutal buoyancy forces generate complex flow in a porous annulus. Beji et al. [4] utilized the Darcy model to investigate double-diffusive convection in a porous annulus. Double-diffusive convection in a porous annulus subjected to mass and heat fluxes was introduced by Marcoux et al. [5]. Nithiarasu et al. [6] generalized a model for non-Darcy flow of double-diffusive convection in a porous annulus. Lee et al. [7] examined the multilayer flow in double-diffusive convection of salt water filled in a rotating annulus. Chen et al. [8] adopted lattice Boltzmann model to investigate the double-diffusive in a vertical annulus by considering opposing gradients of the temperature and concentration. Goyeau [9] studied the double-diffusive convection in a porous cavity below Darcy number and porosity impacts. Chamkha and Al-Mudhaf [10] used the finite difference method to check the effects of an inclination angle and buoyancy ratio on double-diffusive in a tilted porous cavity.

In the recent years, nanofluid has many applications in solar energy systems, nuclear reactors, heat exchangers, and so on [1119]. Bhatti et al. [20] discussed the impacts of chemical reaction and nonlinear thermal radiation on unsteady 3D boundary layer flow of a viscous nanofluid having gyrotactic microorganisms through a stretching porous cylinder. Mohebbi and Rashidi [21] used the Lattice Boltzmann Method (LBM) to study numerically the natural convection flow of nanofluid in shaped enclosure containing a heating obstacle. Chowdhury et al. [22] studied the double-diffusive convection in an enclosure filled by porous media and nanofluid with considering heat generation effects. The ISPH method was adopted by Aly and Raizah [23] to simulate double-diffusive convection in an enclosure filled with a nanofluid.

The topic of the fluid flow and heat transfer in heterogeneous porous media has many applications in remediation, enhanced oil recovery, and geological CO2 reservation. The porous media are heterogeneous porous media when the permeability varies from point to point in a medium [2426].

Since the last decade, mesh-free methods became a well-alternative tool for overcoming disadvantage of mesh methods in simulating free-surface flows and high deformation problems. Lucy [27] and Gingold and Monaghan [28] firstly introduced smoothed particle hydrodynamics to simulate astrophysical problems. The smoothed particle hydrodynamics (SPH) method is employed in several fields [2935]. An incompressible SPH (ISPH) method has the ability for solving several complex problems by showing a good efficiency [3642].

The SPH method is considered an excellent numerical method for the problems of a highly deforming surface/interface such as free-surface flow, dam break, wave dynamics, impact flows, and bubble dynamics. In addition, the SPH method deals with interfaces naturally without any special requirements. The incompressible version of the SPH method entitled the ISPH method and has mesh-free advantages. The ISPH method can handle the double-diffusive convection flow in a complex cavity by an easy way.

In the current research, the ISPH method is used for the simulation of thermo-solutal convection in an outer square shape connected with a center square cavity. The outer square shape and inner square cavity are filled with -water nanofluid, and the inner cavity is saturated by a heterogeneous porous medium. It is found that an augmentation in buoyancy ratio parameter increases the thermo-solutal convection in both of outer square shape and inner cavity. Moreover, the directions of thermo-solutal convection are strongly depending on the values of buoyancy ratio parameter. Presence of a porous medium in an inner square cavity only makes slight effects of Darcy parameter and porous medium levels in maximum velocity field. An incremental in solid volume fraction augments viscosity of water, and consequently, the flow speed decreases.

2. Mathematical Formulations

Figure 1 presents the initial physical model and initial particle distributions. In the current physical model, the outer square shape has a height and the width of a square shape equals 0.1. The inner center cavity is a square cavity with a height , and there are four normal gates between an outer square shape and an inner square cavity. Only the left walls of an outer square shape are carrying and and right walls have and , and the other walls are adiabatic. The nanofluid occupied an outer square shape and inner square cavity, and the nanofluid is modeled by a one-phase model. Heterogeneous porous media are saturated inner square cavity only and the porous media are in thermal equilibrium with nanofluid. Thermophysical properties of water (base fluid) and at are introduced in Table 1.

Time dependent of double-diffusive convection for nanofluid flow according to [43, 44] is

The porous matrix expressions arewhere

The expressions of a nanofluid according to [4547] arewhere is Brownian velocity, is a freezing point of a base fluid, and is Boltzmann’s coefficient.

2.1. Dimensionless Boundary Conditions

The dimensionless boundary conditions arewhere is a normal vector.

3. ISPH Method

This section summarizes the following steps of the ISPH method.

Predicted velocities:

Pressure Poisson equation (PPE):

Corrected velocities:

Thermal energy:

Concentration equation:

Positions are updated as

Shifting technique is

3.1. SPH Formulation

Any hydrodynamic function can be defined aswhere is a smoothing length and is a kernel function. Equation (13) is defined by SPH approximation as

Quintic spline kernel function iswhere . First derivative is approximated by SPH:

First derivatives of velocity and pressure are defined as

Here, first derivative of a velocity is corrected by a kernel gradient normalization as

SPH approximations of second derivatives for the velocity and pressure are

The calculations of FORTRAN-90 code of the ISPH method were performed by Shaheen II supercomputer owned by King Abdullah University of Science and Technology, Thuwal, Saudi Arabia (project number K1467) (https://www.hpc.kaust.edu.sa/content/shaheen-ii).

4. Validation Test

In this section, the validation test of the ISPH method for the problem of natural convection in a partitioned porous cavity has been conducted. Figure 2 shows the temperature profiles in a partitioned porous cavity among the results from the current results of the ISPH method and Beckermann et al. [48]. The temperature profiles were evaluated at three different locations in a cavity: , , and . It is clear that the results of the ISPH method have a well agreement with the experimental and numerical results of Beckermann et al. [48].

5. Results and Discussion

The numerical results are presented in terms of concentration, temperature, and velocity field characteristics in order to get a clear insight on the problem of double-diffusive of a nanofluid in an outer square shape mounted on the four vertical gates in a porous cavity. Figures 35 show the impacts of a buoyancy ratio parameter on the characteristic of the temperature, concentration, and velocity field. It is revealed that an increase in a buoyancy ratio parameter rises the temperature and concentration distributions inside an outer square shape and inner square cavity. Moreover, at opposing flow mode , the double-diffusive convection occurs in the bottom area of an outer square shape and in an inner cavity. While at aiding flow mode , the double-diffusive convection occurs in the top area of an outer square shape and in a cavity. The maximum of the velocity field increases by as the buoyancy ratio parameter increases from to . Moreover, the maximum of the velocity field appears in the left-top area of an outer square shape at opposing flow mode , and it appears in the bottom-left part of an outer square shape at aiding flow mode .

The effects of the Darcy parameter on the characteristic of the temperature, concentration, and velocity field have been shown in Figures 68. As a Darcy parameter represents a useful factor for explaining the fluid flow in the porous materials, a little decrease in the temperature and concentration distributions is occurring when the Darcy parameter decreases from to . Maximum of the velocity field is reducing as a Darcy parameter decreases from to . Moreover, a lower velocity field in the center of a square cavity appears at a lower Darcy parameter . From these figures, it is noted there are slight changes in the characteristics of the temperature, concentration, and velocity field below the effects of a Darcy parameter. The physical reason returns to the presence of a porous medium layer in the inner cavity only, while an outer square shape contains nanofluid only.

Figures 911 introduce the influences of the solid volume fraction on the characteristics of temperature, concentration, and velocity field. It is seen that the temperature and concentration distributions have slight decrease according to adding more concentration of nanoparticles until . The concentration of the nanoparticles is limited to to avoid solidification between the nanofluid and porous media. Adding more concentration of the nanoparticles increases the nanofluid viscosity. As a result, the maximum of the velocity field decreases by as the solid volume fraction increases from 0 to 0.05.

The variations of the porous levels according to different and on the temperature, concentration, and velocity field profiles have been shown in Figures 1214. In Figure 12, the highest distributions of the temperature are obtained at a heterogeneous porous medium and lowest temperature distributions are obtained at a homogeneous porous medium . In Figure 13, it is seen that the concentration distributions are affected strongly by changing the porous medium levels. A heterogeneous porous medium gives the highest concentration distributions in an inner cavity and lowest concentration distributions in an outer square shape. In contrast, a homogeneous porous medium gives the lowest concentration distributions in an inner cavity and the highest concentration distributions in an outer square shape. In Figure 12, there are slight changes in the velocity field below different levels of porous media. The physical reason returns to the presence of porous media in an inner cavity only.

6. Conclusion

The novelty of this work is to simulate double-diffusive convection of nanofluid inside a new novel shape of two square cavities connected by four vertical gates. Three different levels of porous media related to the values of and were conducted. An inner cavity is filled by a heterogeneous porous medium and -water nanofluid. The outer square shape with four gates is filled by a nanofluid only. An incompressible version of the SPH method was utilized to solve the nondimensional governing equations of the current problem. The main finding of the current results can be summarized as follows:(i)An augmentation in a buoyancy ratio parameter increases the temperature and concentration distributions in an outer square shape and in an inner cavity(ii)At opposing flow mode, , and the direction of buoyancy force from upwards to downwards and the maximum of velocity field occur at a top-left part of an outer square shape(iii)At aiding flow mode, , and the direction of buoyancy force from downwards into upwards and the maximum of velocity field occur at a bottom-left part of an outer square shape(iv)A decrease in the Darcy parameter reduces the flow speed in an inner porous cavity, and it has slight effects in an outer square shape(v)Adding more concentration of nanoparticles increases the viscosity of a host fluid, and consequently, the flow speed inside an outer square shape and in an inner cavity decreases(vi)Variation in porous medium levels changes the concentration distributions inside an outer square shape and in an inner cavity

Nomenclature

:Concentration
:Heat capacity
:Darcy parameter
:Effective diffusivity
:Gravitational acceleration,
:Porous height
:Smoothing length
:Permeability
:Boltamann’s coefficient,
:Thermal conductivity,
:Length of an inner square cavity
:Lewis number
:Buoyancy ratio parameter
:Dimensional pressure,
:Dimensionless pressure
:Prandtl number
:Water freezing point
:Dimensional temperature,
:Rayleigh number:
:Kernel gradient normalization
:Dimensional time,
:Mass,
:Distance between particles
:Dimensionless velocity components
:Dimensional velocity components,
:Brownian velocity
:Smoothing function
:Dimensional Cartesian coordinates,
:Dimensionless Cartesian coordinates
Greek symbols
:Effective thermal diffusivity,
:Coefficient of a kernel function
:Thermal expansion coefficient,
:Constant to avoid zero dominator
:Change rate of in direction
:Change rate of in direction
:Dimensionless concentration
:Solid volume fraction
:Porosity
:Relaxation coefficient
:Dynamic viscosity,
:Density,
:Numerical density
:Any hydrodynamic function
:Dimensionless time
:Dimensionless temperature.
Subscripts
:Cold
:Hot
:Fluid
:Nanofluid
:Porous medium.

Data Availability

Data used to support the findings of the study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest associated with this publication.

Acknowledgments

The first author extends his appreciation to the Deanship of Scientific Research at King Khalid University, Abha, Saudi Arabia, for funding this work through the Research Group Project under Grant no. RGP.2/17/42. This work was supported by the National Natural Science Foundation of China (no. 71601072), Key Scientific Research Project of Higher Education Institutions in Henan Province of China (no. 20B110006), and the Fundamental Research Funds for the Universities of Henan Province.