Next Article in Journal
A Collision Reduction Adaptive Data Rate Algorithm Based on the FSVM for a Low-Cost LoRa Gateway
Next Article in Special Issue
Investigation of Efficient Optimization Approach to the Modernization of Francis Turbine Draft Tube Geometry
Previous Article in Journal
Functional Ergodic Time Series Analysis Using Expectile Regression
Previous Article in Special Issue
Steady-State Solutions for Two Mixed Initial-Boundary Value Problems Which Describe Isothermal Motions of Burgers’ Fluids: Application
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Mixed Convection on Two-Phase Nanofluid Flow Past a Vertical Plate in Brinkman-Extended Darcy Porous Medium with Nield Conditions

1
Department of Civil Engineering, College of Engineering, University of Ha’il, Hail 81451, Saudi Arabia
2
Laboratoire d’Ingénierie Géotechnique, Ecole Nationale d’Ingénieurs de Tunis, Université de Tunis El Manar, Tunis 1007, Tunisia
3
Department of Mathematical Sciences, Faculty of Science and Technology, Universiti Kebangsaan Malaysia, UKM, Bangi 43600, Malaysia
4
Department of Mathematics and Social Sciences, Sukkur IBA University, Sukkur 65200, Pakistan
5
Department of Mathematical Sciences, Federal Urdu University of Arts, Science & Technology, Gulshan-e-Iqbal, Karachi 75300, Pakistan
6
Center of Research, Faculty of Engineering, Future University in Egypt, New Cairo 11835, Egypt
7
Department of Mathematics, College of Science, King Khalid University, Abha 62529, Saudi Arabia
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(20), 3918; https://doi.org/10.3390/math10203918
Submission received: 23 September 2022 / Revised: 19 October 2022 / Accepted: 20 October 2022 / Published: 21 October 2022
(This article belongs to the Special Issue Computational Fluid Dynamics II)

Abstract

:
The rapid advancement in technology in recent years has shown that nanofluids are very vital to further development in science and technology. Moreover, many industrial specifications cannot be met by allowing natural convection only, hence the need to incorporate forced convection and natural convection into a single flow regime. The research aims to quantify the mixed convective two-phase flow past a vertical permeable surface in a Brinkman-Extended Darcy porous medium (BEDPM) induced by nanofluid, with heat and mass transfer. In addition, the Nield condition is also incorporated. The model of the problem was initially constructed in the vital form of leading governing equations (LGEs). These LGEs are specifically called partial differential equations (PDEs) (because of two or more independent variables) which were later converted into a set of the single independent variable of ordinary differential equations (ODEs) by implementing the similarity transformations. The set of single independent ODEs was numerically solved via the boundary value problem of fourth-order (bvp4c) technique. The bvp4c is one of the most frequently recommended built-in MATLAB subroutines based on the three-stage Labatto formula. The impact of several physically embedded influential parameters on the fluid flow, along with mass and thermal properties of the nanofluid in a Brinkman-Extended Darcy porous medium for the cases of buoyancy assisting flow (BAF) and buoyancy opposing flow (BOF), were investigated and argued. The numerical outcomes clarify that the porosity parameter reduces the velocity, whereas the concentration and the temperature enhance in the case of the buoyancy assisting and buoyancy opposing flows. In addition, the wall drag force elevates for the larger value of the dimensionless permeability parameter K 1 and the buoyancy ratio parameter N , while it declines for the modified porosity parameter ε 1 .
MSC:
76D05; 76D10

1. Introduction

Heat transfer improvement in industrial and engineering applications has received a lot of interest from researchers in recent years. This is because the reliability of most equipment in several fields, such as heat exchangers and electronic devices, is highly dependent on the phenomenon of heat transport. Due to their very small thermal conductivity, base fluids such as ethylene glycol, water, and oil limit the rate of heat transfer. The aforementioned fluids’ deficiencies are overcome by incorporating a single kind of nano-sized particle into these fluids. Choi and Eastman [1] were the first to address this process and referred to the resulting mixture as a nanofluid. Researchers have considered a variety of nanoparticle combinations, including metal oxides (Al2O3, CuO), semiconductors (SiO2, TiO2), as well as metals (Cu, Fe, Al). Therefore, distinct works on nanoparticles can be found in the books of Das et al. [2], Minkowycz et al. [3], and Minea [4] as well as recent research papers by Zaib et al. [5], Usman et al. [6], Waini et al. [7], Kotresh et al. [8], Khan et al. [9], Abderrahmane et al. [10,11,12].
In contrast, employing nanoliquid saturated in a porous media is a well-organized method to develop the characteristics of convection heat transfer for industrial processes. As a result, researchers have paid close attention to the porous media procedure. This is because in several engineering areas, including groundwater, oil flow filtration, thermal insulation, and entire types of heat exchangers, a similar type of scheme was encountered. It can also be found in combustion schemes, cooling of electronic equipment, and turbine blades cooling, where heating or cooling is needed. The mixing of high and low-energies fluids has a great impact on the ability of these devices. Merkin [13] investigated the mixed convective flow through a vertical plate embedded in a porous media and discovered that the presence of double solutions occurs in the opposing flow. Later, Ahmad and Pop [14] stretched this problem to nanofluid flow. The mixed convective flow through a horizontal circular cylinder induced by nanofluid saturated in a porous media was inspected by Nazar et al. [15]. They claimed that as the solid volume fraction increases, the magnitude of the friction factor decreases. Chamkha et al. [16] investigated the radiation impact on the fluid flow of mixed convective packed with a nanofluid past a wedge embedded in a porous medium. Roşca et al. [17] developed a mixed convective flow past a vertical flat plate through a porous medium filled with a nanofluid. According to their numerical results, the nanoparticle volume fraction reduces the velocity for the buoyancy assisting flow at the plate. Sheremet and Pop [18] investigated mixed convective flow within a square lid-driven cavity induced by a nanofluid, whereas Bakar et al. [19] scrutinized the radiation effect on the mixed convective flow through a nanofluid past a vertical permeable cylinder saturated in a porous medium, concluding that the nanoparticle of Al2O3 travels quicker in separating the thicknesses of the boundary layer, followed by Cu and TiO2. Waini et al. [20] found double solutions of mixed convective flow with hybrid nanofluid from a vertical flat surface immersed in a porous media.
Several realistic applications rely on heat and mass transfer, including food processing, automotive radiators, refrigeration, air conditioning, petrochemical refining, and microwave ovens. Dash and Das [21] examined the Hall current impact on magneto flow via a porous accelerated plate, along with the transfer of heat and mass characteristics induced by heat generation. The time-dependent free convective flow induced by rotating elastic viscous fluid past a porous vertical plate in a porous media with heat and mass transfer (HT and MT) is premeditated by Panda et al. [22]. Das et al. [23] explored the time-dependent flow with MT via a porous accelerated plate and obtained the solution numerically as well as analytically. Sanjayanand and Khan [24] examined the combined effects of heat and mass transfer on the flow induced by viscoelastic fluid past a stretched exponential sheet with viscous dissipation. The impact of mass transfer on the time-dependent flow past a stretchable sheet subjected to a porous medium and chemically reactive species was examined by Elbashbeshy et al. [25]. The agglomeration effects of mass and heat transfer on thin fluid film induced by nanofluid through a time-dependent stretchable sheet with a viscous dissipation effect were considered by Qasim et al. [26]. They observed that the heat transfer declines due to viscous dissipation as well as for nanofluids parameters. The exploration of Brownian motion on magneto flow of a Casson nanofluid from a stretchable sheet with irregular radiation was scrutinized by Ramudu et al. [27]. Recently, Alam and Marwat [28] inspected MT and HT in a viscous flow past an erratic permeable plate with irregular thickness.
The use of combined forced and natural convection is trending in industries and manufacturing sectors for heat and mass transmission in a fluid. The applications range from the installation of grain storage and electronic circuits to the extraction of oil and heat exchangers. In many of these applications, buoyancy forces are generated due to the diffusion of masses of different chemical species and pressure gradients at different spatial positions. The effect of buoyancy on heat and mass transfer is very significant and has been studied by many authors. Roşca et al. [29] analyzed the effects of buoyancy in a flow under a combined presence of forced and natural convection in a porous medium. The model is developed with the notion that the free stream velocity u is not equal to zero. A numerical analysis was conducted on the governing equations, and the results indicate a dual solution exists when the flow is opposing. Mahanthesh et al. [30] studied the effects of heat and mass transfer on natural convection flow. With the aid of Laplace transform, a closed-form solution was obtained for the proposed model. They found out that the flow and heat transfer can be controlled by adjusting the nanoparticle volume fraction. Waini et al. [31] administered the flow over a traveling needle-like surface and also established that the solution is non-unique for a certain range of values of the flow parameters. Similar results were observed, in that the two solution branches respond oppositely to the flow parameters. Meanwhile, Bakar et al. [32] administered the flow over a shrinking surface in the presence of slip factor and thermal radiation. Dual solutions were also observed and the results remained similar for the two solutions. The stability analysis shows that the lower branch is unstable while the branch of the upper solution is stable.
Based on a comprehensive assessment of the previously stated review of the literature, it is worth noting that no comprehensive exploration of nanofluid caused by mixed convection instabilities developed through a vertical plate subjected to a porous medium with infinite extensions under the Nield condition has been conducted. As a result of the numerous applications of nanofluids in modern nanotechnology and the cooling of sophisticated electronic and industrial equipment, this research not only provides the interactions of nanofluids over a vertical plate immersed in a Brinkman-extended Darcy porous medium, but also investigates the impact of Nield conditions on the transport phenomenon. The boundary layer equation is scaled down into the non-dimensional form by using some similarity transformation and the resulting non-dimensional equation is solved.

2. Mathematical Scenario of the Problem

Consider a steady flow near an impermeable vertical flat plate subject to mixed convection embedded in a porous medium filled with nanofluid flow flowing with non-uniform free-stream velocity u ( x ) = A x and zero mass flux conditions. The analysis of HT and MT with the Brinkman-extended Darcy model was also incorporated in the given formulation, where we shall assume the velocity square term has been neglected. As portrayed in Figure 1, where the Cartesian coordinates x and y measured, respectively, along with the plate and normal to it. The plate is assumed to be kept at a constant temperature T w , while T signifies the temperature of the surrounding fluid, where T w > T corresponds to the heated plate (buoyancy assisting flow) while T w < T corresponds to the cooling plate (buoyancy opposing flow). Moreover, the concentration C at a specific constituent on the fluid side of the plate varies from C w to C in the solution through porous media with C w > C . The following fundamental governing equations, which we can write down, were developed using the aforementioned conditions with the help of boundary-layer and Boussinesq approximations.
u x + v y = 0 ,
u u x + v u y = u d u d x + ε 2 υ e f f 2 u y 2 ε 2 υ f K ( u u ) + ε 2 g [ ( β T ) f ( T T ) + ( β C ) f ( C C ) ] ,
u T x + v T y τ [ D B T y C y + D T T ( T y ) 2 ] = α f 2 T y 2 ,
u C x + v C y D T T ( 2 T y 2 ) = D B 2 C y 2 .
For the current physical model, the subjected boundary conditions (BCs) are:
{ u = 0 , v = 0 , T = T w , D B C y + D T T T y = 0 ,   at   y = 0 , u u ( x ) , T T , C C   as   y .
Here, ( u , v ) denote the components of velocity in ( x , y ) axes, T signifies the temperature, g signifies the acceleration because of gravity, μ f ,   ρ f signify the dynamic viscosity and the density, respectively, α f thermal diffusivity, ( ( β T ) f , ( β C ) f ) signify the corresponding coefficients of thermal and concentration volumetric expansions, ( D B , D T ) signify the coefficients of Brownian and the thermophoretic diffusions, τ signifies the ratio of the effective heat capacity of the solid nanoparticles to the heat capacity of the working base fluid’s and K signifies the permeability of the porous medium.
To further simplify the examination of the model, substantial similarity transformations are introduced [13]:
η = y u x α f , ψ = x u α f F ( η ) , G ( η ) = T T T w T , S ( η ) = C C C w C ,
where ψ indicates the stream function, characterized by v = ψ / x and u = ψ / y . Further, differentiating the stream function to get the following simplest form of the velocity component can take place as follows:
u = u F ( η )   and   v = u α f x ( F η F ) .
For making the governing similarity equations, we have to assume the values of ( β T ) f and ( β C ) f in the form ( β T ) f = x ( β T ) f * and ( β C ) f = x ( β C ) f * (see Mohamed et al. [33]), where ( β T ) f * and ( β C ) f * are arbitrary constants. Therefore, the aforementioned Equation (6) implies two buoyancy or mixed convection parameters λ λ T and N λ λ C , respectively, defined as:
λ T = ε 2 g ( β T ) f * ( T w T ) A 2 = R a x T P e x 2 Pr   and   λ C = ε 2 g ( β C ) f * ( C w C ) A 2 = R a x C P e x 2 Pr ,
where R a x T = ε 2 g ( β T ) f ( T w T ) x 3 υ e f f α f and R a x C = ε 2 g ( β C ) f ( C w C ) x 3 υ e f f α f are the respective Rayleigh numbers, P e x = u x α f is the Peclet number, Pr = υ e f f α f is the Prandtl number, and the buoyancy or mixed convection ratio factor, N , which clarifies the rivalry between the concentration and thermal mixed convective forces, is demarcated as;
N = ( β C ) f * ( C w C ) ( β T ) f * ( T w T ) .
Furthermore, the choice of negative or positive selection of ( β C ) f * (see Bejan and Khair, [34]) suggests that the buoyancy ratio parameter will be negative or positive (see Equation (9)). On the other hand, the extraordinary condition λ C = N λ excludes the character of the term mixed convective ratio constraint from the set of governing non-dimensional equations. Consequently, we examine the current model further by focusing solely on the rest of the comprised influential parameters. Hence, by substituting Equation (6) into the governing equations, we obtain:
ε 1 F + F F F 2 + 1 K 1 ( F 1 ) + ( G + N S ) λ T = 0 ,
G + F G + N b S G + N t G 2 = 0 ,
S + L e F S + N t N b G = 0 ,
along with boundary conditions (BCs) are:
{ F ( 0 ) = 0 , F ( 0 ) = 0 , G ( 0 ) = 1 , N b S ( 0 ) + N t G ( 0 ) = 0   at   η = 0 , F ( η ) 1 , G ( η ) 0 , S ( η ) 0   as   η .
In the above equations, the dimensionless control parameters are the Brownian diffusion parameter, the thermophoresis parameter, the modified porosity parameter, the dimensionless permeability parameter, and the Lewis number. These are defined in the mathematical notations as:
N b = τ ( C w C ) D B / α f , N t = τ ( T w T ) D T / α f T , ε 1 = ε 2 υ e f f / α f = ε 2 Pr , K 1 = ε 2 υ f / K A , L e = α f / D B .

2.1. Skin Friction Coefficient

The skin friction coefficients or friction factor is defined as follows:
( P e x 1 / 2 / Pr ) C f = μ f ( u y ) y = 0 ρ f u 2 = F ( 0 ) .

2.2. Heat Transfer

The heat transfer is demarcated as follows:
P e x 1 / 2 N u x = x ( k f T y ) y = 0 k f ( T w T ) = G ( 0 ) ,
where P e x = u x / α f corresponds to the Peclet number. It is appealing to observe that due to modified conditions such as N b S ( 0 ) + N t G ( 0 ) = 0 , the expression of mass transfer identically takes place to zero specified by Kuznetsov and Nield [35].

3. Numerical Methodology and Validation of the Code

In the present research work, the governing equations for a steady flow with heat and mass transfer near an impermeable vertical flat plate comprising buoyancy impact saturated in a porous media occupied with a two-phase nanofluid flow model along with subjected non-uniform free-stream velocity and zero mass flux conditions are mathematically modeled. These modeled sets of PDEs were changed into (ordinary) similarity equations by implementing the similarity variables. Mathematically, the composed form of the set of differential Equations (10)–(12) along with appropriate BCs (13) was coupled and highly nonlinear, therefore, the solution in the closed explicit form of these ODEs was complicated to find. For this reason, it is possible to use a well-known numerical process to which we may quickly revert to its effective strategies for achieving the emphasized goals or aims. In terms of computation, the bvp4c (boundary value problem of fourth-order) finite–difference approach is one of the most suggested built-in codes available in MATLAB sub-routines (Shampine et al. [36]), which may be used successfully to offer fourth-order correct results in this phenomenon. Therefore, the given method can be used to solve first-order initial value problems. To continue our scheme, here, we present the new-fangled notations which can be read as follows:
F = E 1 , F = E 2 , F = E 3 , G = E 4 , G = E 5 , S = E 6 , S = E 7 .
The nonlinear first-order set of ODEs can be documented due to the assessment of these new mathematical notations in the acquired similarity equations, which can take place as:
( E 1 E 2 E 3 E 4 E 5 E 6 E 7 ) = ( E 2 E 3 ( E 2 2 E 1 E 3 1 + K 1 ( E 2 1 ) λ T ( E 4 + N E 6 ) ) ε 1 E 5 E 1 E 5 N b E 5 E 7 N t E 5 2 E 7 L e E 1 E 7 N t N b E 5 ) ,
and to measure the aforementioned Equation (17) to be a single-point boundary value problem (IVP). Therefore, the subject ICs (known and unknown) are premeditated by:
{ E 1 ( 0 ) = 0 ,     E 2 ( 0 ) = 0 ,     E 3 ( 0 ) = A 1 , E 4 ( 0 ) = 1 ,     E 5 ( 0 ) = N b N t E 7 ( 0 ) ,   E 6 ( 0 ) = A 2 , E 7 ( 0 ) = N t N b E 5 ( 0 ) .
Additionally, the values of unknown arbitrary constants A 1 and A 2 are assumed to be simple for problems with a single solution. Moreover, the procedure of the scheme had been illustrated practically and graphically by Chu et al. [37] and Khan et al. [38], respectively, in a clear-cut way: that the practice of the code is constantly repeated until the far-field BCs are held at once:
{ E 2 ( η ) = 0 ,   E 4 ( η ) = 0 ,   E 6 ( η ) = 0 .
The mesh selection and error process are bolstered in this case by the continuous solution of residual. The error tolerance was fixed to 10−6. Moreover, the edge of the boundary layer thickness η = η max is taken to be a fixed finite value of seven, which can hold the border ailments (13) asymptotically for the variations of the distinct parameters to maintain the necessary accuracy, while the mesh-size selection value has been adjusted between 0.001 and 0.1. Based on these observations, several numerical values of the skin friction coefficient for varying values of   ε 1 and   λ T when   K 1 = 1.0 (without the influence of the other comprised controlling parameters) are quantitatively compared with the available former work of Rosali et al. [39] and Khan et al. [40] as elucidated completely in Table 1. This tabular demonstration shows that the current outcomes and the comparative ones are in superb settlement, indicating that the results supplied by the bvp4c numerical approach are unquestionably real and accurate.

4. Analysis of Results

This fragment of the research work is divulged particularly to demonstrate the detailed description of the results in the form of distinct tables as well as several graphs. The acquired tabular form and graphical behavior of the outcomes were constructed in the requisite form of the dimensionless velocity F ( η ) , the concentration distribution S ( η ) , the temperature distribution G ( η ) , the shear stress or friction factor coefficient C f   , and the heat transfer N u x , as clarified in Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13 as well as in Table 2 and Table 3. For the computational simulations, we have fixed the range, such as ( η < η = 7.0 ) , to develop the asymptotically converged behavior of the curves owing to growing values of the embedded control flow parameters such as ε 1 (modified porosity parameter), λ T (mixed convection parameter), K 1 (dimensionless permeability parameter), N t (thermophoresis parameter), N b (Brownian motion parameter), L e (Lewis number), and N (buoyancy ratio parameter), whose default or fixed values are ε 1 = 0.5 , K 1 = 1.2 , N t = 0.1 , N b = 0.1 , L e = 10 , and N = 0.5 . Meanwhile, the mixed convection parameter λ T value has been taken to be 3 and −3 for the corresponding buoyancy assisting and buoyancy opposing flows, respectively, where the Prandtl number is taken to be fixed, such as 6.2 for the considered working base fluid (water).
The influence of the modified porosity parameter  ε 1 on the dimensionless velocity, F ( η ) , the temperature distribution, G ( η ) , and the concentration distribution, S ( η ) , for the buoyancy assisting and opposing flows against the pseudo-similarity variable η is schematically revealed in Figure 2, Figure 3 and Figure 4, respectively. From Figure 2, it is seen that F ( η ) solution curves shrinkages for the phenomenon of BOF as well as BAF owing to the escalating values of ε 1 . Furthermore, the momentum boundary layer thickness for the case of opposing flow is lower as compared to the phenomenon of the buoyancy assisting flow owing to the larger value of the modified porosity parameter ε 1 . Graphically, the gap in the outcome curves for the buoyancy assisting flow (BAF) is more relative to the case of the buoyancy opposing flow (BOF) when we boosted up the significant impacts of the modified porosity parameter. Physically, it is observed that the effective kinematic viscosity of the nanofluidic medium upsurges directly with the larger value of the modified porosity parameter; as a result, the velocity of the nanofluid and the thickness of the momentum boundary layer decelerate as we move forward against the pseudo-similarity variable. On the other hand, the profiles of temperature and concentration distributions are continuously enriched for both cases of BAF and BOF, due to the augmentation of modified porosity parameter ε 1 , as clarified amazingly in the image form (see Figure 3 and Figure 4), respectively. Further, it is noted from the outcomes that the thickness of the concentration and thermal boundary layer (TBL) escalates with higher ε 1 . Meanwhile, the gap for the solutions of temperature in both cases of mixed convection parameter looks completely similar, whereas for the solutions of concentration the gap in curves of BAF is more relative to the BOF, since the resistance to the nanofluid flow elevates with superior consequences of the modified porosity parameter, which ultimately upgraded the solution domain possibility for concentration and temperature profiles.
Figure 2. F ( η ) against ε 1 .
Figure 2. F ( η ) against ε 1 .
Mathematics 10 03918 g002
Figure 3. G ( η ) against ε 1 .
Figure 3. G ( η ) against ε 1 .
Mathematics 10 03918 g003
Figure 4. S ( η ) against ε 1 .
Figure 4. S ( η ) against ε 1 .
Mathematics 10 03918 g004
Figure 5, Figure 6 and Figure 7 elucidate the influence of the dimensionless permeability parameter K 1 on profiles of F ( η ) , G ( η ) and S ( η ) of the nanofluidic medium flow for the cases of BAF and BOF against the pseudo-similarity variable η , respectively. From the outcome of these images, it is detected that the velocity curves decline for BAF and incline for BOF due to the augmentation in the value of the permeability parameter while the behavior of the temperature and concentration distribution profiles are completely changed, and hence, ultimately, it is upsurges and shrinkages for both cases of BAF and BOF, respectively. In addition, the momentum BLT diminished but the thermal and concentration boundary layer elevates for the BAF due to the growing value of K 1 while the entire profile tendency was opposite for the case of BOF owing to the augmentation in K 1 . In general, a big dimensionless permeability parameter K 1 is produced by a small permeability of a porous medium K . Because the fluid may travel more freely due to the high permeability of a porous medium, a decrease in K 1 , raises the velocity profile for the case of BAF. Since as velocity increases the temperature falls, and vice versa, as a result, as K 1 decreases the temperature decreases. In addition, regarding the improving values of K 1 , the gap between the solution curves of the concentration profiles for both cases of BAF and BOF are comparatively larger than the gap available in the solution curves of the temperature distribution and velocity profile.
Figure 5. F ( η ) against K 1 .
Figure 5. F ( η ) against K 1 .
Mathematics 10 03918 g005
Figure 6. G ( η ) against K 1 .
Figure 6. G ( η ) against K 1 .
Mathematics 10 03918 g006
Figure 7. S ( η ) against K 1 .
Figure 7. S ( η ) against K 1 .
Mathematics 10 03918 g007
The consequence of the thermophoresis diffusion parameter N t on the temperature distribution G ( η ) and concentration distribution S ( η ) for the cases of BAF and BOF are emphasized in Figure 8 and Figure 9, respectively. From Figure 8, it seems to be worth mentioning that for both cases of buoyancy or mixed convection parameter, the ascending motion of nanomaterials helps to simplify heat transport filled with a nanofluidic medium by increasing G ( η ) and thickening the related region of the boundary layer. Alternatively, the concentration profile upsurges for the case of BAF due to the larger N t , while for the case of BOF, the tendency of concentration solution curves initially declines and then is ultimately boosted (see Figure 9). In addition, the gap in the solution curves of the temperature profile is slightly lesser than the gap available in the concentration profile curves because of the enhancing value of the thermophoresis parameter N t .
Figure 8. G ( η ) against N t .
Figure 8. G ( η ) against N t .
Mathematics 10 03918 g008
Figure 9. S ( η ) against N t .
Figure 9. S ( η ) against N t .
Mathematics 10 03918 g009
Figure 10 and Figure 11 clarify the behavior of the temperature distribution G ( η ) and concentration distribution S ( η ) for the BAF and BOF, owing to the larger impacts of the Brownian motion parameter N b , respectively. Outcomes reveal that the occurrence of Brownian mass diffusion term encourages the thermo-migration of fluid particles confines of a nanofluidic medium. Owing to this well-known statement, the temperature distribution G ( η ) elevates for the case of BAF and decelerates for BOF due to the escalating values of N b , as emphasized in Figure 10. On the other hand, the concentration profile declines for the BAF and BOF if we increase the influences of N b , as shown graphically in Figure 11. Further, this behavior of the outcomes was generally acceptable due to the superior values of N b , which abruptly lessen the contribution of N t in the concentration equation because of the diffusive term   ( N t / N b ) ( d 2 S / d η 2 ) . Therefore, N t and N b divulge the opposite tendency towards ϕ ( ξ ) .
Figure 10. G ( η ) against N b .
Figure 10. G ( η ) against N b .
Mathematics 10 03918 g010
Figure 11. S ( η ) against N b .
Figure 11. S ( η ) against N b .
Mathematics 10 03918 g011
From an engineering perspective, the behavior of the outcomes of the shear stress C f   and heat transfer N u x for sundry values of the influential control parameters due to the cases of BAF and BOF are illustrated graphically and quantitatively in Figure 12 and Figure 13 and Table 2 and Table 3, respectively. Graphically, Figure 12 and Figure 13 exemplify the influence of the modified porosity factor   ε 1 on the shear stress F ( 0 ) and heat transfer G ( 0 ) profile for the BAF and BOF against the dimensionless permeability parameter K 1 , respectively. Results demonstrate the shear stress F ( 0 ) and heat transfer G ( 0 ) shrinkages for the case of BAF in the entire domain of K 1 , escalating the modified porosity parameter   ε 1 , while for the case of BOF the F ( 0 ) and G ( 0 ) initially grows and then declines for the rest of the domain against K 1 . In the quantitative approach, Table 2 marks a significant improvement in shear stress C f   due to the boosted value of N (buoyancy ratio parameter), and K 1 (dimensionless permeability parameter) for the cases of the buoyancy or mixed convection factor. At the same time, a diminishing behavior of the magnitude of friction factor is depicted for the case of BAF as well as BOF at the vertical plate surface of the porous medium owing to larger impacts of   ε 1 (modified porosity parameter). In the thermal properties, Table 3 reveals an improvement in the local Nusselt number N u x with the escalating value of the factors   ε 1 (modified porosity factor), N t (thermophoresis diffusion parameter), L e (Lewis number), K 1 (dimensionless permeability parameter), and N (buoyancy ratio parameter) for the case of BOF, while it is declined due to N b (Brownian motion parameter). However, a weakening pattern was noticed in the heat transfer N u x for the larger support in the value of the constraints   ε 1 (modified porosity parameter), N t (thermophoresis parameter), N b (Brownian motion parameter), N b (Lewis number), and K 1 (dimensionless permeability parameter) for the case of BOF, while it upsurges due to N (buoyancy ratio parameter).
Figure 12. Variations of the shear stress for several values of ε 1 .
Figure 12. Variations of the shear stress for several values of ε 1 .
Mathematics 10 03918 g012
Figure 13. Variations of the heat transfer for several values of ε 1 .
Figure 13. Variations of the heat transfer for several values of ε 1 .
Mathematics 10 03918 g013

5. Conclusions

The present research work was carried out to investigate the steady flow with heat and mass transfer near an impermeable vertical plate comprising the impact of mixed convection embedded in a porous media filled with a two-phase nanofluid flow model. The non-uniform free-stream velocity and Nield conditions are also incorporated in the current mathematical modeling of the problem. Therefore, the notable findings derived from the considered research work are as follows:
  • Outcomes of the given model are discussed throughout the research work for the cases of buoyancy assisting flow (BAF) and buoyancy opposing flow (BOF) with several varying embedded control parameters.
  • The profile of velocity declines for the case of BAF as well as BOF due to the larger impacts of ε 1 , while the trend of velocity is the same for the case of BAF and the opposite for the case of BOF with the larger value of the dimensionless permeability parameter K 1 .
  • The dimensionless permeability parameter K 1 and the Brownian motion parameter N b raise the temperature for the case of BAF and decline for the case of BOF. On the other hand, the temperature abruptly improves for both cases of mixed convection parameter due to N t and ε 1 .
  • The concentration profile decelerates for the case of BAF as well as BOF owing to higher impacts of parameter N b , whilst the modified porosity parameter ε 1 provides a substantial increment for both cases of buoyancy parameter.
  • The local skin friction coefficient elevates for the larger value of K 1 and N while it is declined for ε 1 .
  • The higher impacts of the buoyancy ratio parameter N enrich the wall mass transfer for the cases of BAF as well as for BOF.
It is worth mentioning that the work addressed above can be further extended by involving the time-dependent flow or involving some different aspects such asnon-linear radiation, the application of solar energy, viscous dissipation, or chemical reaction. Moreover, the geometry has been changed from the vertical plate to a static wedge or cylinder.

Author Contributions

Conceptualization, U.K. and A.Z.; methodology, U.K. and A.Z.; software, U.K. and A.Z.; validation, U.K., H.G. and A.Z.; formal analysis, A.I., U.K., A.Z. and H.G.; investigation, A.I., Z.R. and S.M.E.; resources, Z.R.; data curation, A.I., H.G. and S.M.E.; writing—original draft preparation, A.I., H.G., S.M.E. and Z.R.; writing—review and editing, U.K., H.G., S.M.E. and Z.R.; visualization, Z.R.; supervision, A.I.; project administration, Z.R. and S.M.E.; funding acquisition, Z.R. and S.M.E. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the research center of the Future University in Egypt, 2022.

Data Availability Statement

Not applicable.

Acknowledgments

The authors are thankful for the partially support of the research center of the Future University in Egypt, 2022.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Choi, S.U.S.; Eastman, J.A. Enhancing Thermal Conductivity of Fluids with Nanoparticles. In Proceedings of the 1995 International Mechanical Engineering Congress and Exhibition, San Francisco, CA, USA, 12–17 November 1995. [Google Scholar]
  2. Das, S.K.; Choi, S.U.S.; Yu, W.; Pradeep, T. Nanofluids: Science and Technology; Wiley-Interscience: Hoboken, NJ, USA, 2007. [Google Scholar]
  3. Minkowycz, W.J.; Sparrow, E.M.; Abraham, J.P. Nanoparticle Heat Transfer and Fluid Flow; CRC Press: New York, NY, USA; Taylor & Fracis Group: Abingdon, UK, 2013. [Google Scholar]
  4. Minea, A.A. Advances in New Heat Transfer Fluids: From Numerical to Experimental Techniques; CRC Press: New York, NY, USA; Taylor & Fracis Group: Abingdon, UK, 2017. [Google Scholar]
  5. Zaib, A.; Haq, R.U.; Sheikholeslami, M.; Khan, U. Numerical analysis of effective Prandtl model on mixed convection flow of γAl2O3–H2O nanoliquids with micropolar liquid driven through wedge. Phys. Scr. 2020, 95, 035005. [Google Scholar] [CrossRef]
  6. Usman, M.; Zubair, T.; Hamid, M.; Haq, R.U. Novel modification in wavelets method to analyze unsteady flow of nanofluid between two infinitely parallel plates. Chin. J. Phys. 2020, 66, 222–236. [Google Scholar] [CrossRef]
  7. Waini, I.; Ishak, A.; Pop, I. Nanofluid flow on a shrinking cylinder with Al2O3 nanoparticles. Mathematics 2021, 9, 1612. [Google Scholar] [CrossRef]
  8. Kotresh, M.J.; Ramesh, G.K.; Shashikala, V.K.R.; Prasannakumara, B.C. Assessment of Arrhenius activation energy in stretched flow of nanofluid over a rotating disc. Heat Transf. 2021, 50, 2807–2828. [Google Scholar] [CrossRef]
  9. Khan, U.; Zaib, A.; Bakar, S.A.; Ishak, A. Stagnation-point flow of a hybrid nanoliquid over a non-isothermal stretching/shrinking sheet with characteristics of inertial and microstructure. Case Stud. Therm. Eng. 2021, 26, 101150. [Google Scholar] [CrossRef]
  10. Abderrahmane, A.; Al-Khaleel, M.; Mourad, A.; Laidoudi, H.; Driss, Z.; Younis, O.; Guedri, K.; Marzouki, R. Natural convection within inversed T-shaped enclosure filled by nano-enhanced phase change material: Numerical investigation. Nanomaterials 2022, 12, 2917. [Google Scholar] [CrossRef]
  11. Abderrahmane, A.; Qasem, N.A.A.; Mourad, A.; Al-Khaleel, M.; Said, Z.; Guedri, K.; Younis, O.; Marzouki, R. Enhancing the melting process of shell-and-tube PCM thermal energy storage unit using modified tube design. Nanomaterials 2022, 12, 3078. [Google Scholar] [CrossRef]
  12. Abderrahmane, A.; Younis, O.; Al-Khaleel, M.; Laidoudi, H.; Akkurt, N.; Guedri, K.; Marzouki, R. 2D MHD mixed convection in a zigzag trapezoidal thermal energy storage system using NEPCM. Nanomaterials 2022, 12, 3270. [Google Scholar] [CrossRef]
  13. Merkin, J.H. Mixed convection boundary layer flow on a vertical surface in a saturated porous medium. J. Eng. Math. 1980, 14, 301–313. [Google Scholar] [CrossRef]
  14. Ahmad, S.; Pop, I. Mixed convection boundary layer flow from a vertical flat plate embedded in a porous medium filled with nanofluids. Int. Commun. Heat Mass Transf. 2010, 37, 987–991. [Google Scholar] [CrossRef]
  15. Nazar, R.; Tham, L.; Pop, I.; Ingham, D.B. Mixed convection boundary layer flow from a horizontal circular cylinder embedded in a porous medium filled with a nanofluid. Transp. Porous Media 2011, 86, 517–536. [Google Scholar] [CrossRef]
  16. Chamkha, A.J.; Abbasbandy, S.; Rashad, A.M.; Vajravelu, K. Radiation effects on mixed convection over a wedge embedded in a porous medium filled with a nanofluid. Transp. Porous Media 2012, 91, 261–279. [Google Scholar] [CrossRef]
  17. Roşca, N.C.; Roşca, A.V.; Groşan, T.; Pop, I. Mixed convection boundary layer flow past a vertical flat plate embedded in a non-Darcy porous medium saturated by a nanofluid. Int. J. Numer. Methods Heat Fluid Flow 2014, 24, 970–987. [Google Scholar] [CrossRef]
  18. Sheremet, M.A.; Pop, I. Mixed convection in a lid-driven square cavity filled by a nanofluid: Buongiorno’s mathematical model. Appl. Math. Comput. 2015, 266, 792–808. [Google Scholar] [CrossRef]
  19. Bakar, S.A.; Arifin, N.M.; Md Ali, F.; Bachok, N.; Nazar, R.; Pop, I. A stability analysis on mixed convection boundary layer flow along a permeable vertical cylinder in a porous medium filled with a nanofluid and thermal radiation. Appl. Sci. 2018, 8, 483. [Google Scholar] [CrossRef] [Green Version]
  20. Waini, I.; Ishak, A.; Groşan, T.; Pop, I. Mixed convection of a hybrid nanofluid flow along a vertical surface embedded in a porous medium. Int. Comm. Heat Mass Transf. 2020, 114, 104565. [Google Scholar] [CrossRef]
  21. Dash, G.C.; Das, S.S. Hall effects on MHD flow along an accelerated porous flat plate with mass transfer and internal heat generation. Math. Eng. Indust. 1999, 7, 389–404. [Google Scholar]
  22. Panda, J.P.; Dash, G.C.; Das, S.S. Unsteady free convective flow and mass transfer of a rotating elastico-viscous liquid through porous media past a vertical porous plate. AMSE J. Mod. Meas. Cont. B 2003, 72, 47–59. [Google Scholar]
  23. Das, S.S.; Sahoo, S.K.; Dash, G.C. Numerical solution of mass transfer effects on unsteady flow past an accelerated vertical porous plate with suction. Bull. Malays. Math. Sci. Soc. 2006, 29, 33–42. [Google Scholar]
  24. Sanjayanand, E.; Khan, S.K. On heat and mass transfer in a viscoelastic boundary layer flow over an exponentially stretching sheet. Int. J. Therm. Sci. 2006, 45, 819–828. [Google Scholar] [CrossRef]
  25. Elbashbeshy, E.M.A.; Emam, T.G.; Abdel-Wahed, M.S. Mass transfer over unsteady stretching surface embedded in porous medium in the presence of variable chemical reaction and suction/injection. Appl. Math. Sci. 2011, 5, 557–571. [Google Scholar]
  26. Qasim, M.; Khan, Z.H.; Lopez, R.J.; Khan, W.A. Heat and mass transfer in nanofluid thin film over an unsteady stretching sheet using Buongiorno’s model. Eur. Phys. J. Plus 2016, 131, 16. [Google Scholar] [CrossRef]
  27. Ramudu, A.C.V.; Kumar, K.A.; Sugunamma, V.; Sandeep, N. Heat and mass transfer in MHD Casson nanofluid flow past a stretching sheet with thermophoresis and Brownian motion. Heat Transf. 2020, 49, 5020–5037. [Google Scholar] [CrossRef]
  28. Alam, A.; Marwat, D.N.K. Heat and mass transfer on a stretching/ shrinking and porous sheet of variable thickness with suction and injection. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2021, 235, 5297–5308. [Google Scholar] [CrossRef]
  29. Roşca, A.V.; Roşca, N.C.; Pop, I. Mixed convection heat and mass from a vertical surface embedded in a porous medium. Transp. Porous Media 2015, 109, 279–295. [Google Scholar] [CrossRef]
  30. Mahanthesh, B.; Gireesha, B.J.; Gorla, R.S.R. Heat and mass transfer effects on the mixed convective flow of chemically reacting nanofluid past a moving/stationary vertical plate. Alex. Eng. J. 2016, 55, 569–581. [Google Scholar] [CrossRef] [Green Version]
  31. Waini, I.; Ishak, A.; Pop, I. Dufour and Soret effects on Al2O3-water nanofluid flow over a moving thin needle: Tiwari and Das model. Int. J. Num. Meth. Heat Fluid Flow 2020, 31, 766–782. [Google Scholar] [CrossRef]
  32. Bakar, S.A.; Md Arifin, N.; Khashi’ie, N.S.; Bachok, N. Hybrid nanofluid flow over a permeable shrinking sheet embedded in a porous medium with radiation and slip impacts. Mathematics 2021, 9, 878. [Google Scholar] [CrossRef]
  33. Mohamed, M.K.A.; Salleh, M.Z.; Ishak, A. Effects of viscous dissipation on mixed convection boundary layer flow past a vertical moving plate in a nanofluid. J. Adv. Res. Fluid Mech. Therm. Sci. 2020, 69, 1–18. [Google Scholar] [CrossRef]
  34. Bejan, A.; Khair, K.R. Heat and mass transfer by natural convection in a porous medium. Int. J. Heat Mass Transf. 1985, 28, 909–918. [Google Scholar] [CrossRef]
  35. Kuznetsov, A.V.; Nield, D.A. Natural convective boundary-layer flow of a nanofluid past a vertical plate. Int. J. Therm. Sci. 2010, 49, 243–247. [Google Scholar] [CrossRef]
  36. Shampine, L.F.; Gladwell, I.; Thompson, S. Solving ODEs with Matlab; Cambridge University Press: Cambridge, UK, 2003. [Google Scholar]
  37. Chu, Y.M.; Khan, U.; Shafiq, A.; Zaib, A. Numerical simulations of time-dependent micro-rotation blood flow induced by a curved moving surface through conduction of gold particles with non-uniform heat sink/source. Arab. J. Sci. Eng. 2021, 46, 2413–2427. [Google Scholar] [CrossRef]
  38. Khan, U.; Zaib, A.; Ishak, A.; Abu Bakar, S.; Animasaun, I.L.; Yook, S.-J. Insights into the dynamics of blood conveying gold nanoparticles on a curved surface when suction, thermal radiation, and Lorentz force are significant: The case of Non-Newtonian Williamson fluid. Math. Comput. Simul. 2022, 193, 250–268. [Google Scholar] [CrossRef]
  39. Rosali, H.; Ishak, A.; Pop, I. Mixed convection stagnation-point flow over a vertical plate with prescribed heat flux embedded in a porous medium: Brinkman-extended Darcy formulation. Trans. Porous Media 2011, 90, 709–719. [Google Scholar] [CrossRef]
  40. Khan, U.; Zaib, A.; Khan, I.; Nisar, K.S.; Baleanu, D. Insights into the stability of mixed convective Darcy–Forchheimer flows of cross liquids from a vertical plate with consideration of the significant impact of velocity and thermal slip conditions. Mathematics 2020, 8, 31. [Google Scholar] [CrossRef]
Figure 1. Physical model of the problem.
Figure 1. Physical model of the problem.
Mathematics 10 03918 g001
Table 1. Validations of   C f for varying values of   ε 1 and   λ T when   K 1 = 1.0 and the rest of the controlling parameters are absent.
Table 1. Validations of   C f for varying values of   ε 1 and   λ T when   K 1 = 1.0 and the rest of the controlling parameters are absent.
  ε 1   λ T   C f
Rosali et al. [39]Khan et al. [40]Present Results
0.1−0.54.15084.13894.13883456
1.06.48746.48646.48737453
2.07.76117.76147.76132675
0.5−0.51.88211.88381.88374657
1.02.85972.84532.84523589
2.03.39443.39443.39431567
0.7−0.51.59671.60081.60073245
1.02.40742.41242.41232897
2.02.85142.84992.84984563
1.0−0.51.34181.34881.34371245
1.02.00502.00502.00501345
2.02.36902.36202.36201965
Table 2. Numerical values of shear stress C f for several values of influential parameters.
Table 2. Numerical values of shear stress C f for several values of influential parameters.
ε 1 K 1 N Assisting Flow Opposing Flow
0.250.750.506.1487376−0.98035728
0.50 4.1754402−0.75291722
0.75 3.3192327−0.59739661
1.00 2.8170001−0.47427638
0.500.750.504.1754402−0.75291722
0.80 4.1854688−0.74086792
0.85 4.1955118−0.72308883
0.90 4.2055694−0.69391372
0.500.750.404.1623374−0.76864007
0.434.1662691−0.76368910
0.464.1702001−0.75894785
0.504.1754402−0.75291722
Table 3. Numerical values of heat transfer N u x for several values of influential parameters.
Table 3. Numerical values of heat transfer N u x for several values of influential parameters.
ε 1 N t N b L e K 1 N Assisting FlowOpposing Flow
0.250.100.10100.750.500.754349840.21760305
0.500.707140960.22880305
0.750.676747550.23700267
1.000.654252120.25756980
0.500.100.10100.750.500.707140960.22880305
0.300.650153450.23965033
0.500.600021870.23109820
0.700.555590800.21769050
0.500.100.10100.750.500.707140960.22880305
0.110.706819430.22631427
0.120.706551170.22514745
0.130.706323940.22402853
0.500.100.102.00.750.500.717763620.15705795
4.00.711194580.18919402
10.00.707140960.22880305
15.00.706252850.24466209
0.50.10.1100.750.50.717763620.22880305
0.800.706977650.23609696
0.850.706818680.24595134
0.900.706664440.26005883
0.50.10.1100.750.400.706433040.22295484
0.430.706645630.22485683
0.460.706858040.22662686
0.500.707140960.22880305
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gasmi, H.; Khan, U.; Zaib, A.; Ishak, A.; Eldin, S.M.; Raizah, Z. Analysis of Mixed Convection on Two-Phase Nanofluid Flow Past a Vertical Plate in Brinkman-Extended Darcy Porous Medium with Nield Conditions. Mathematics 2022, 10, 3918. https://doi.org/10.3390/math10203918

AMA Style

Gasmi H, Khan U, Zaib A, Ishak A, Eldin SM, Raizah Z. Analysis of Mixed Convection on Two-Phase Nanofluid Flow Past a Vertical Plate in Brinkman-Extended Darcy Porous Medium with Nield Conditions. Mathematics. 2022; 10(20):3918. https://doi.org/10.3390/math10203918

Chicago/Turabian Style

Gasmi, Hatem, Umair Khan, Aurang Zaib, Anuar Ishak, Sayed M. Eldin, and Zehba Raizah. 2022. "Analysis of Mixed Convection on Two-Phase Nanofluid Flow Past a Vertical Plate in Brinkman-Extended Darcy Porous Medium with Nield Conditions" Mathematics 10, no. 20: 3918. https://doi.org/10.3390/math10203918

APA Style

Gasmi, H., Khan, U., Zaib, A., Ishak, A., Eldin, S. M., & Raizah, Z. (2022). Analysis of Mixed Convection on Two-Phase Nanofluid Flow Past a Vertical Plate in Brinkman-Extended Darcy Porous Medium with Nield Conditions. Mathematics, 10(20), 3918. https://doi.org/10.3390/math10203918

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

Article Metrics

Back to TopTop