Next Article in Journal
Interlayer Microstructure Analysis of the Transition Zone in the Silicon/Perovskite Tandem Solar Cell
Next Article in Special Issue
Using Non-Fourier’s Heat Flux and Non-Fick’s Mass Flux Theory in the Radiative and Chemically Reactive Flow of Powell–Eyring Fluid
Previous Article in Journal
Cost Assessment of a Tokamak Fusion Reactor with an Inventive Method for Optimum Build Determination
Previous Article in Special Issue
Details on the Hydrothermal Characteristics within a Solar-Channel Heat-Exchanger Provided with Staggered T-Shaped Baffles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mechanism of Solute and Thermal Characteristics in a Casson Hybrid Nanofluid Based with Ethylene Glycol Influenced by Soret and Dufour Effects

1
Faculty of Mechanical Engineering and Ship Technology, Gdańsk University of Technology, Narutowicza 11/12, 80-233 Gdańsk, Poland
2
Institute of Structural Analysis, Poznan University of Technology, Piotrowo 5 Street, 60-965 Poznan, Poland
3
Department of Applied Mathematics and Statistics, Institute of Space Technology, P.O. Box 2750, Islamabad 44000, Pakistan
4
Section of Mathematics, International Telematic University Uninettuno, Corso Vittorio Emanuele II, 39, 00186 Roma, Italy
5
Department of Statistics and Operations Research, College of Science, King Saud University, P.O. Box 2455, Riyadh 11451, Saudi Arabia
*
Authors to whom correspondence should be addressed.
Energies 2021, 14(20), 6818; https://doi.org/10.3390/en14206818
Submission received: 8 September 2021 / Revised: 14 October 2021 / Accepted: 15 October 2021 / Published: 18 October 2021
(This article belongs to the Special Issue Recent Advances in Solar Energy Collectors: Models and Applications)

Abstract

:
This article models a system of partial differential equations (PDEs) for the thermal and solute characteristics under gradients (concentration and temperature) in the magnetohydrodynamic flow of Casson liquid in a Darcy porous medium. The modelled problems are highly non-linear with convective boundary conditions. These problems are solved numerically with a finite element approach under a tolerance of 10 8 . A numerical algorithm (finite element approach) is provided and a numerical procedure is discussed. Convergence is also observed via 300 elements. Simulations are run to explore the dynamics of flow and the transport of heat and mass under parametric variation. To examine the impact of a temperature gradient on the transport of mass and the role of a concentration gradient on the transport of heat energy, simulations are recorded. Remarkable changes in temperature and concentration are noted when Dufour and Soret numbers are varied.

1. Introduction

Several techniques for the optimisation of thermal transfer have been used in practice. The recent development in nanotechnology makes it possible to synthesize solid structures at a nano-size. Furthermore, investigators have shown that the suspension of nano-sized solid structures is important in fluids, observing a rise in the thermal conductivity of liquids and noting that fluids with nanostructures behave as suitable coolants. These fluids with nano-sized particles are called nanofluids. The thermal performance of such coolants is optimized due to the dispersion of such nanostructures. Consequently, energy losses can be minimized. Based on this fact, researchers have studied heat transfer in nanofluids and observed a remarkable enhancement in that heat transfer. For example, Masuda et al. [1] presented in the early 1990s the concept of nanometer-sized particles. Choi [2] worked on nano-sized particles in the Argonne National Laboratory of USA and introduced the term nanofluid, a fluid with nano-sized solid structures. Phelan et al. [3] published on the applications related to nanoparticles in fluid during the heat transport phenomenon. Lee et al. [4] devised the strategy for measuring the impact of oxide nanoparticles on the thermal conductivity of fluids. Eastman et al. [5] also conducted a study on enhancing a fluid’s thermal conductivity due to the immersion of nanoparticles in it. Huaqing et al. [6] analyzed the effects of the suspension of alumina particles on the increase in thermal conductivity of a fluid. Yimin and Li [7] published significant outcomes related to the suspension of nanoparticles in fluid and their impact on the effectiveness of thermal conductivity. Phillbot et al. [8] discussed how the dispersion of nanoparticles causes thermal enhancement. Naseem et al. [9] discussed an enhancement in temperature performance, inserting impacts of nanoparticles past a stretching surface via a shooting approach. Nazir et al. [10] investigated the comparative performance of thermal aspects among hybrid nanoparticles and nanoparticles using the finite element method. Koriko et al. [11] performed a bio-convection flow in a thixotropic nanofluid under a magnetic field considering the stretching surface. Ali et al. [12] investigated the consequences of peristaltic flow in a Jeffrey liquid using hybrid nanoparticles along with the slip conditions under the action of an MHD flow. Tian et al. [13] studied 3D configurations of hybrid nanoparticles using the impacts of magnetic fields in fins. Mumraiz et al. [14] performed a study of entropy generation in the presence of hybrid nanoparticles along with a variable heat flux. Awais et al. [15] discussed the impacts of pressure drop in thermal energy including nanofluids. Nazir et al. [16] used the finite element method to consider the consequences of thermal energy in a Carreau—Yasuda liquid when inserting hybrid nanoparticles past a heated cone.
The transport of heat during industrial processing plays an essential role in a product being manufactured. For example, food processing, the welding of joints of mechanical tools, the fabrication of materials, casting processes, and coating processes, etc., are engineering applications where the process of heat transport is an essential factor. Theoretical and experimental investigations have been published regarding industrial applications in various emerging applied fields where heat transfer is encountered. It is for this reason, for instance, that Bibi et al. [17] developed the boundary value problems associated with thermal transport in a Williamson fluid exposed to a magnetic field and thermal radiation. They observed that thermal transport can be increased by dispersing nanomaterials in the fluid. Ijaz et al. [18] studied the entropy generation in a Sisko fluid when joule heating and activation energy are significant. Majeed et al. [19] discussed the transfer of heat in a fluid between two coaxially rotating disks. They used non-Fourier’s law to study the impact of thermal relaxation on the transport of heat energy. Ali et al. [20] analyzed the effect of thermal relaxation time on the transfer of heat energy in a Jeffery fluid. Tanveer et al. [21] studied the peristaltic phenomenon in fluid theoretically. They examined the role of the contraction and expansion of muscles on food transport and the transfer of heat. Tanveer et al. [22] examined the dynamics of the flow of blood in a micro-channel. Khan et al. [23] modelled the movement of fluid and heat transfer in MHD (magnetohydrodynamic) flow over the parabolic surface. Abbas et al. [24] discussed the extension of the Yamada—Ota and Xue models for a micropolar fluid subject to a stagnation point. Rahman et al. [25] considered implementing the finite element method (FEM) for simulations associated with heat transfer under hydrodynamic forces in a grooved channel having two partially heated circular cylinders. Zahri et al. [26] performed numerical simulations related to the thermally magnetized rectangular chamber optimization of a practically heated continuous stream. Hayat et al. [27] discussed the thermal radiation emitted by the MHD (magnetohydrodynamic) flow of a Maxwell fluid during thermal energy transport. Hayat et al. [28] examined the impact of the role of nanoparticles on the transfer of heat energy in a couple stress fluid squeezed by two parallel plates. Saif et al. [29] also studied the impact of nanomaterials in optimizing the transport of heat in a second-grade fluid. Hayat et al. [30] discussed the effect of the diffusion of wall temperature into a fluid with nanoparticles over a non-uniformly moving surface.
A porous medium makes a vital impact in adjusting the thickness of the MBL (momentum boundary layer). Therefore, several researchers have considered a porous medium in the modelling of thermal processes in flowing fluids. Some recent references are given in [31,32,33].
Solute and thermal characteristics have been shown to have significant effects in various studies in view of engineering and industrial processes. Several researchers have studied heat transfer phenomena caused by the concentration and temperature gradients. The mechanism of heat transfer that occurs due to the concentration gradient is called the Dufour effect, whereas the mechanism of heat transfer that happens due to the temperature gradient is called the Soret effect. Therefore, heat transfer phenomena make a vital impact on the Dufour and Soret theories [34,35,36,37,38]. Naseem et al. [39] studied the impacts of thermal energy, including the variable properties under a magnetic field, considering the Soret and Dufour impacts over a melting surface. Naseem et al. [40] discussed the characterizations of thermal energy in the presence of the Soret and Dufour effects.
The literature reveals that the Lorentz force, porous mediums, convective boundary conditions, and buoyant force in the presence of a dispersion of nano-metallic structures have not been studied yet. Therefore, a constant magnetic field and the Soret and Dufour effects are herein considered. This article consists of five sections. Comprehensive literature examples are given in the introduction section. The development system of non-linear partial differential equations in a porous medium exposed to a magnetic field is given in Section 2 where boundary conditions are studied as well. The problems are solved numerically using the shooting technique. This is an appropriate technique as concluded in Refs. [41,42,43,44]. This numerical procedure is explained in Section 3. The graphical consequences are addressed in Section 4. Section 5 visualizes the prime findings of the current model.

2. Physical Statement of the Problem

A phenomenon of heat and mass transfer in a Casson liquid is observed towards the vertical melting surface. A constant magnetic field is exerted along the y-axis of the surface. The fluid runs due to the stretching of the surface moving along with velocity ( u = a x ). The Nano and hybrid nanoparticles have an important role when inserted into the base fluid for simulating the enhancement of heat and solute in fluid particles. The external heat source/sink is also considered in the heat equation carrying Joule heating, viscous dissipation, and the Dufour and Soret effects. The geometry of fluid behavior is considered in Figure 1. The following assumptions are considered:
A constant magnetic field is considered,
Dufour and Soret effects are addressed,
A porous vertical surface is considered,
Casson fluid particles are inserted,
A phenomenon of Joule heating is noticed,
The finite element method is used,
Nano and hybrid nanoparticles are observed.
The following PDEs (partial differential equations) are used utilizing BLA (boundary layer thickness):
u x + v y = 0 ,
u u x + v u y = 1 ρ h n f ( μ h n f + p y 2 π c ) 2 u 2 y + ( β h n f ) T g ( T T ) + ( β h n f ) c g ( C C ) σ h n f B 0 2 u ρ h n f μ h n f u K 1 ,
T x + v T y = K h n f ( ρ c p ) h n f 2 u 2 y + 1 ( ρ c p ) h n f ( μ h n f + p y 2 π c ) ( u y ) 2 + σ h n f B 0 2 u 2 ( ρ c p ) h n f + Q 0 ( ρ c p ) h n f ( T T ) + D K T u C s C p 2 C 2 y ,
u C x + v C y = ( D B ) h n f 2 C 2 y + D T T 2 T 2 y ,
Equation (1) is called the 2D continuity equation for incompressible and steady flow. It is noticed that boundary layer approximations and assumptions are used in the law of conservation mass. Equation (2) is known as the momentum equation in x and y directions for steady and incompressible flow. Equation (2) is derived from momentum conservation law considering the magnetic field, Casson fluid and the Bouncy effect. Equations (3) and (4) are energy and concentration equations considering effects (heat generation, Joule heating, chemical reaction, Soret and Dufour) over a vertical surface.
Slip conditions are considered and derived using the slip theory concept. Slip theory is defined as the velocity of the fluid at the boundary that is not equal to the velocity of the boundary. The mathematical form of slip conditions near the surface of the boundary is defined as:
a x = u ( x , 0 ) , v ( x , 0 ) = 0 , γ k h n f T y ( x , 0 ) = h f ( T T ( x , 0 ) , T ( x , ) = T γ D h n f C y ( x , 0 ) = h c ( C C ( x , 0 ) ) , 0 = u ( x , 0 ) , C ( x , ) = C , } ,
Composite between hybrid nanoparticles and nanomaterials [45,46] are:
ρ h n f ρ f = ( 1 φ 2 ) { ( 1 φ 1 ) + φ 1 ρ s 1 ρ f } + φ 2 ρ s 2 ρ f ,
( ρ c p ) h n f ( ρ c p ) f = ( 1 φ 2 ) { ( 1 φ 1 ) + φ 1 ( ρ c p ) s 1 ( ρ c p ) f } + φ 2 ( ρ c p ) s 2 ( ρ c p ) f ,
μ h n f μ f = 1 ( 1 φ 1 ) 2.5 ( 1 φ 2 ) 2.5 ,
σ h n f σ b f = σ s 2 + ( n 1 ) σ b f ( n 1 ) φ 2 ( σ b f σ s 2 ) σ s 2 + 2 σ b f + φ 2 ( σ b f σ s 2 ) ,   σ b f σ f = σ s 1 + 2 σ f ( n 1 ) φ 1 ( σ f σ s 1 ) σ s 1 + 2 σ f + φ 1 ( σ f σ s 1 ) ,
k h n f k b f = k s 2 + 2 k b f ( n 1 ) φ 2 ( k b f k s 2 ) k s 2 + 2 k b f + φ 2 ( k b f k s 2 ) ,   k b f k = k s 1 + 2 k f ( n 1 ) φ 1 ( k f k s 1 ) k s 1 + 2 k f + φ 1 ( k f k s 1 ) ,
D h n f D f = 1 ( 1 φ 1 ) ( 1 φ 2 ) ,
Required transformations are:
u = a x f ( η ) ,   v = a υ f f ( η ) , η = a υ f y , θ ( η ) = T T T w T , ψ ( η ) = a υ f x f ( η ) ,   ϕ ( η ) = C C C w C
Hence, one gets:
υ h n f ν f ( 1 + 1 β ) f [ ( f ) 2 f f ] + ( G r ) T θ ( η ) + ( G r ) C ϕ ( η ) M ( σ h n f σ f ) ( ρ f ρ h n f ) f K f = 0 , ( K h n f K f ) ( ρ f ρ h n f ) ( ( C p ) f ( C p ) h n f ) θ + ( υ h n f ν f ) ( ( C p ) f ( C p ) h n f ) ( 1 + 1 β ) P r E c ( f ) 2
+ ( σ h n f σ f ) ( ( C p ) f ( C p ) h n f ) ( ρ f ρ h n f ) M E c (   f ) 2 + f P r θ + P r β * θ + P r D f ϕ = 0 ,
ϕ + S c S r θ + S c f ϕ =   0 ,
Associated BCs are formulated as:
f ( 0 ) = 0 ,       θ ( 0 ) = β i γ ( 1 θ ( 0 ) ) ,     ϕ ( 0 ) = E i α ( 1 ϕ ( 0 ) ) ,   f ( 0 ) = 1   f ( ) = θ ( ) = ϕ ( ) = 0     .   }
( G r ) t = g ( T w T ) ( β h n f ) T a 2 x , ( G r ) c = g ( C w C ) ( β h n f ) C a 2 x , M = σ f B 0 2 ρ f a , K = μ f a K 1 E c = a 2 x 2 c p ( T w T ) ,   P r = μ f ( c p ) f k f , β * = Q 0 a ( ρ c p ) f , D f = D k T ( C w C ) C s C p ( T w T ) , S c = D T ( T w T ) ( 1 φ 1 ) ( 1 φ 2 ) υ f T D f ( C w C ) , S c = υ f ( 1 φ 1 ) ( 1 φ 2 ) D f   , B i = h f k h n f υ f a . }
It is noticed that the physical magnitude reveals a dynamically mutual fluid-structure interaction at a rigid interface in which the strength of viscous frictional forces are evaluated as the trough skin friction coefficient called factor C f . Surface force is also called the skin friction coefficient. The skin friction coefficient, Nusselt number and concentration gradient are:
C f = τ x y | y = 0 ρ f U 0 2 = f ( 0 ) R e x 1 2 ( 1 φ 1 ) 2.5 ( 1 φ 2 ) 2.5 , N u = x k h n f T y | y = 0 k f ( T w T ) = R e x 1 2 k h n f k f θ ( 0 ) , S h = x D h n f C y | y = 0 D f ( C w C ) = R e x 1 2 ( 1 φ 1 ) ( 1 φ 2 ) ϕ ( 0 ) , }
where R e x = a x 2 ν f is the Reynolds number. Thermo-physical properties are shown in Table 1. Moreover, τ is shear stress (force per unit area) applied perpendicular to the surface of the sheet.

3. Numerical Procedure

The concept of the FEM is generated by structure mechanics. The vital role of the FEM is to handle complex geometries, unstructured grids and curved cells. An important phenomenon of the FEM is the discretization of the domain into elements. The FEM is a strong approach in terms of convergence, accuracy and stability rather than other numerical methods. Detailed discussion of the Galerkin finite element method can be found in [10,44,45].
The discretization of the problem domain into finite elements: Firstly, the domain of the current model is divided into the finite number of elements and a polynomial kind of solution is assumed over each element. Such approximation is then applied in the weak form (with the weighted residual scheme) which is derived from the strong form (Equations (13)–(16)). Therefore, the shape functions are multiplied with the residuals and integrated over the whole domain. The residuals of the current model are:
η e η e + 1 W 1 [ f S ] d η = 0 ,
η e η e + 1 W 2 [ υ h n f ν f ( 1 + 1 β ) S [ ( S ) 2 f S ] + ( G r ) T θ ( η ) + ( G r ) C ϕ ( η ) M ( σ h n f σ f ) ( ρ f ρ h n f ) S K S ] d η = 0 ,
η e η e + 1 W 3 [ ( K h n f K f ) ( ρ f ρ h n f ) ( ( C p ) f ( C p ) h n f ) θ ' ' + P r β * θ + ( σ h n f σ f ) ( ( C p ) f ( C p ) h n f ) ( ρ f ρ h n f ) M E c ( S ' ) 2 + f P r θ ' P r D f ϕ ' ' + ( υ h n f ν f ) ( ( C p ) f ( C p ) h n f ) ( 1 + 1 β ) P r E c ( S ' ) 2 ] d η = 0 ,
η e η e + 1 W 4 [ ϕ + S c S r θ + S c ( 1 φ 2 ) 2.5 ( 1 φ 1 ) 2.5 f ϕ ] d η = 0 ,
where W 1 ,   W 2 ,   W 3 and W 4 are called the shape functions, f is considered as S for the Galerkin approach. Shape functions and variables used in the problem are:
χ j = ( 1 ) l 1 1 η η l 1 1 η l η l 1 ,
where l is considered as 1 and 2, respectively:
f = l = 1 2 ( f l χ j ) ,   F = l = 1 2 ( F l χ j ) ,  
θ = l = 1 2 ( θ l χ j ) ,   ϕ = l = 1 2 ( ϕ l χ j ) .
Problem domain: The computational domain is made with the help of the problem domain while the problem domain is considered [ 0 ,   ) . It is observed that asymptotic BCs are satisfied by assuming   η = 6 . Therefore, in this current model, the computation domain is [ 0 ,   6 ] .
Assembly process: This procedure plays a vital role to obtain the stiffness matrix, force vector, and integral boundary vector over the whole domain of the problem. Finally, the global stiffness matrix is simulated with the help of the assembly process. The Galerkin approach is used to find the stiffness elements that are:
K i j 11 = η e η e + 1 χ i ( d χ j d η ) d η , K i j 12 = η e η e + 1 χ i ( χ j ) d η ,
K i j 13 = 0 ,   K i j 14 = 0 , K i j 21 = 0 , K i j 42 = 0 ,
K i j 23 = η e η e + 1 ( G r ) T χ i ( χ j ) d η , B i 1 = 0 , B i 4 = 0 ,
K i j 24 = η e η e + 1 ( G r ) C χ i ( χ j ) d η , B i 2 = 0 , B i 3 = 0 ,
K i j 33 = η e η e + 1 [ ( K h n f K f ) ( ρ f ρ h n f ) ( ( C p ) f ( C p ) h n f ) d χ j d η d χ i d η + P r β * χ i ( χ j ) + F ¯ P r χ i d χ j d η ] d η ,
K i j 22 = η e η e + 1 [ ν h n f ν f ( 1 + 1 β ) d χ j d η d χ i d η S ¯ χ i d χ j d η f ¯ χ i d χ j d η ] M ( σ h n f σ f ) ( ρ f ρ h n f ) χ i ( χ j ) K χ i ( χ j ) ] d η ,
K i j 31 = η e η e + 1 [ ( υ h n f ν f ) ( ( C p ) f ( C p ) h n f ) ( 1 + 1 β ) P r E c S ¯ χ i d χ j d η ( σ h n f σ f ) ( ( C p ) f ( C p ) h n f ) ( ρ f ρ h n f ) M E c S ¯ χ i d χ j d η ] d η ,
K i j 34 = η e η e + 1 [ P r D f d χ j d η d χ i d η ] d η , K i j 32 = 0 , K i j 34 = 0 ,
K i j 43 = η e η e + 1 ( d χ j d η d χ i d η ) d η , K i j 51 = 0 , K i j 41 = 0 ,
K i j 43 = η e η e + 1 [ d χ j d η d χ i d η + f ¯ χ i d χ j d η S c ( 1 φ 2 ) 2.5 ( 1 φ 1 ) 2.5 ] d η ,
Picard linearization approach: This approach is used to convert non-linear equations to linear equations. This approach is called the linearization approach. In the current model, the Picard linearization approach is utilized to obtain the linear system of equations:
S ¯ = j = 1 2 χ j S j ¯ ,   f ¯ = j = 1 2 χ j f j ¯ .
Here S j ¯ and f j ¯ are nodal (unknown) values. The flow chart of this Picard linearisation approach is captured in Figure 2.
The Picard approach is used to simulate the heat conduction problems related to the local fractional. This approach is also used to determine the non-differential solution of the heat equation in view of the fractal media and non-differential solution. One of the most vital advantages is that it converts highly non-linear equations into linear equations. The engineering and heat transfer problems are associated with PDEs (partial differential equations). The Picard approach is used to simulate such equations in terms of the engineering and industrial fields.
Error analysis: The error analysis is simulated by the following equation that is:
E r r o r = | τ j τ j 1 | ,
and convergence criteria is mentioned as:
M a x | τ j τ j 1 | < 10 5 .
Meshfree analysis: In this analysis, convergence analysis is tested in terms of a grid-independent study. This analysis is carried out in Table 2. The simulations are recorded at the mid of each element. The solution converges after the simulations of 300 elements. Hence, 300 elements are required for convergence of the problem. The mesh-free analysis is known as a grid-independent study. The outcomes are recorded of f ( η m a x 2 ) ,   θ ( η m a x 2 ) and ϕ ( η m a x 2 ) at the mid of each 300 elements. Hence, the solution of the problem is converged at 300 elements. Table 3 reveals the comparative results of the Nusselt number with published work [40].
Comparative study: It is observed that the results of the present problem are simulated with the help of the finite element method. Table 3 reveals a validation of results in view of the Nusselt number for the case of nanofluid. In Table 3, published results [47] are simulated by a shooting approach whereas these results are verified with the finite element method. Table 4. Shows the comparative study of FEM and shooting method.

4. Outcomes and Discussion

Flow analysis: The effects of pertinent physical parameters are displayed in Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15 and Figure 16. This parameter is defined as the ratio between the viscous force and buoyancy force.
Variation in the velocity of nanofluids under the influences of ( G r ) c is presented in Figure 3. Figure 3 reveals that the velocity of the nanofluids is an increasing function of ( G r ) c and the reason for this is that when ( G r ) c increases, it loosens up the intermolecular forces within the nanofluids and thereby enhances the velocity of the nanofluids. Figure 4 displays the behavior of ( G r ) t on the nanofluids’ temperature. It is observed in Figure 4 that the temperature of the nanofluids’ has a direct relation with ( G r ) t , which is due to an enhancement in the nanofluids’ volumetric expansion. The effect of K   on the velocity of nanofluids is displayed in Figure 5. Figure 5 reports that the velocity of the nanofluids is a decreasing function of K . Physically, K enhances the resistance to nanofluids causing a reduced velocity of the nanofluids. Figure 6 demonstrates the impact of β (this parameter appears due to using a Casson liquid in the flow analysis) on the velocity profile of the nanofluids. It is observed in Figure 6 that the velocity of the nanofluids is a decreasing function of β . Physically, when β is added to the nanofluids, the friction between particles slows. An inverse proportional relation is addressed among the Casson number and the motion of fluid particles. In Figure 1, Figure 2, Figure 3 and Figure 4, it is noted that the velocity profiles of the hybrid nanofluid ( M O S 2 S i O 2 ) is higher than the nanofluid ( M O S 2 ). Figure 5 and Figure 6 elucidate the impacts of the magnetic field (the magnetic field is called the vector field near the region of the electric field in forces related to where the magnetic field is observed) on the velocity and temperature of nanofluids. It is noted in Figure 5 and Figure 6 that velocity is enhanced with an increase in the magnetic field, while an opposite trend is found in the temperature profile. Physically, the Lorentz force produced due to the magnetic field is in the same direction of the nanofluids’ motion, so the flow boundary layer thickness increases and the thermal boundary layer decreases. Figure 5 also shows that the velocity of the nanofluid ( M O S 2 ) is higher when compared to the hybrid nanofluid ( M O S 2 S i O 2 ), while an opposite trend is noted in the temperature profile (Figure 6).
Temperature analysis: Figure 7 demonstrates the effect of E c (this parameter expresses ration among the enthalpy difference and the flow’s kinetic energy) on the temperature of nanofluids. It is observed from Figure 7, that E c has an inverse relation with the temperature profile. The reason for this is that the heat transfer rate at the surface increases, thus, the thermal boundary layer thickness decreases with increases in the numerical value of E c . Figure 8 depicts the impact of β on the temperature of nanofluids. The temperature of nanofluids decreases with the increase of β . The external heat source has an important impact on the temperature gradient of nanofluids, resulting in a decrement in both the thermal state of the nanofluid and temperature distribution. The variation effect of Bi on the nanofluids’ temperature is plotted in Figure 9. Figure 9 shows that the temperature of the nanofluids’ increases with an increase in the Bi. Physically, the temperature gradient near the surface increases when increasing the numerical values of Bi and hence the thickness of the thermal boundary layer increases. Figure 10 reveals the variation in temperature with the effect of P r (fractional between momentum to thermal boundary layers). It is noted from Figure 10 that the temperature of the nanofluids’ decreases as the numerical value of P r is increased. This occurred because the thermal conductivity of the nanofluids’ decreased with an increase in the P r . Figure 7, Figure 8, Figure 9 and Figure 10 report that the temperature profile is greater in the case of the hybrid nanofluid ( M O S 2 S i O 2 ). The impact of S r on the concentration profile is depicted in Figure 11. It is observed from Figure 11 that   S r has a direct relation with the concentration profile. The reason for this is that the S r impact refers to the mass flux from a lower to greater solute concentration produced by the temperature gradient. The variation in the concentration profile under the influences of S c is presented in Figure 12. It is observed in Figure 12 that the concentration profile decreased with an increased numerical value of S c (ratio between mass diffusivity and momentum diffusivity). Physically, the kinematic viscosity increases with an increase in the S c which turns the molecular diffusion, therefore, the velocity of the nanofluids’ decreases.
Concentration analysis: Figure 13 shows that the concentration profile increases versus the higher values of S r .  Figure 14 reveals a decreasing trend into solute particles when S c is increased. The role of E i versus concentration is observed in Figure 15 while the concentration increases. The reason for this is that as E i increases, mass diffusivity decreases and as a result of the concentration of nanofluids increases. Figure 16 reveals the variation of the concentration profile under the influences of α . It is noted from Figure 16 that the concentration profile has an inverse relation with α .
Analysis of the Sherwood, skin friction coefficient and Nusselt number: Table 5 presents the surface force, Nusselt number and rate of mass diffusion for ( M O S 2 S i O 2 ) and nanofluid ( M O S 2 ). The following remarks have been noted from Table 5. Under the effect of E i The local skin friction coefficient, Nusselt number and Sherwood number is higher in the case of the hybrid nanofluid ( M O S 2 S i O 2 ). Surface force is more significant in the case of the hybrid nanofluid ( M O S 2 S i O 2 ) with the effect of B i and ( G r ) t . Under the effects of B i and ( G r ) t , the rate of heat transfer is higher in the nanofluid ( M O S 2 ) as compared to the hybrid nanofluid ( M O S 2 S i O 2 ).

5. Key Points of the Problem

The presented model, related to thermo-physical correlations, is developed for selected characteristics of heat and mass diffusions under the influence of viscous and Joule heating, and temperature and concentration gradients. The theoretical concepts are implemented utilizing the finite element method, with a comprehensive parametric study also presented. The following outcomes are crucial:
  • Momentum boundary thickness is decreased against a variation in a magnetic field.
  • The adjusting of BLT is controlled by varying the magnitude of the magnetic number.
  • The fluid magnetic field interactions with hybrid nanofluid particles are more significant than the magnetic fluid interaction in the case of the nanofluid. Therefore, Joule heating in the hybrid nanofluid is more significant than that in the nanofluid.
  • The role of dissipation of the thermal energy and heat generation is helpful for an enhancement in thermal performance. Further, BLT is inclined versus the variation of dissipation of thermal energy and heat generation.
  • The rate of solute particles is inclined versus enhancement in the Soret number.
  • In the porous medium, drag force exists due to the flow end; hence the convective transfer of heat and mass is compromised.
  • Maximum production of heat energy is achieved for the case of hybrid nanoparticles rather than the production of heat energy for nanoparticles.
  • Maximum acceleration is produced in the motion of particles for hybrid nanoparticles rather than the case of nanoparticles.
  • Temperature and concentration gradients are significantly boosted for hybrid nanoparticles rather than nanoparticles.
  • Convergence of the problem is obtained for 300 elements.

Author Contributions

Conceptualization, methodology and software, M.B.H.; validation, and formal analysis, U.N.; data curation, writing—original draft preparation, W.S.; project administration, H.A. and funding acquisition, S.A. All authors have read and agreed to the published version of the manuscript.

Funding

Research Supporting Project number (RSP-2021/167), King Saud University, Riyadh, Saudi Arabia.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

SymbolsUsed for
v ,   u Velocities
μ Fluid viscosity
p y Yield stress
β Casson fluid number
T ,   T Temperature and ambient temperature
σ Electrical conductivity
C p Specific heat capacitance
K 1 Porosity number
D Mass diffusion
b f Base fluid
γ Thermal slip
φ 1 ,   φ 2 ,   φ Volume fractions
C Ambient concentration
( G r ) c Grashof number
Infinity
P r Prandtl number
M Magnetic number
D f Dufour number
S r Soret number
a Stretching rate in x-direction
N u Nusselt number
τ Shear stress
y ,   x Space coordinates
T w Wall temperature
PDEsPartial differential equations
g Gravitational force
θ ,   ϕ Dimensionless temperature and concentration
B 0 Magnetic field
k Thermal conductivity
ODEsOrdinary differential equations
C s Concentration susceptibility
D m
h n f ,   n f Hybrid nanofluid and nanofluid
s 1 ,   s 2 Solid particles
( G r ) t Grashof number
η Independent variable
E c Eckert number
β * Heat generation number
S h Sherwood number
C f Skin friction coefficient
R e Reynolds number
S c Schmidt number
E G Ethylene glycol
B L A Boundary layer approximation

References

  1. Masuda, H.; Ebata, A.; Teramae, K. Alteration of thermal conductivity and viscosity of liquid by dispersing ultra-fine particles. Dispersion of Al2O3, SiO2 and TiO2 ultra-fine particles. Netuse Bussei 1993, 7, 227–233. [Google Scholar] [CrossRef]
  2. Choi, S.U.; Eastman, J.A. Enhancing Thermal Conductivity of Fluids with Nanoparticles; Argonne National Lab: Lemont, IL, USA, 1995. [Google Scholar]
  3. Phelan, P.E.; Bhattacharya, P.; Prasher, R.S. Nanofluids for heat transfer applications. Annu. Rev. Heat Transf. 2005, 14, 255–275. [Google Scholar] [CrossRef]
  4. Lee, S.; Choi, S.U.-S.; Li, S.; Eastman, J. Measuring thermal conductivity of fluids containing oxide nanoparticles. J. Heat Transf. 1999, 121, 280–289. [Google Scholar] [CrossRef]
  5. Eastman, J.A.; Jeffery, A.; Choi, S.U.; Shaoping, L.; Thompson, L.J.; Lee, S. Enhanced thermal conductivity through the development of nanofluids. MRS Online Proc. Libr. Arch. 1996, 457, 220–236. [Google Scholar] [CrossRef] [Green Version]
  6. Huaqing, X.; Wang, J.; Tonggeng, X.; Liu, Y.; Wu, Q. Thermal conductivity enhancement of suspensions containing nano-sized alumina particles. J. Appl. Phys. 2002, 91, 4568–4572. [Google Scholar]
  7. Yimin, X.; Li, Q. Heat transfer enhancement of nanofluids. Int. J. Heat Fluid Flow 2000, 21, 58–64. [Google Scholar]
  8. Keblinski, P.; Phillpot, S.; Choi, S.; Eastman, J. Mechanisms of heat flow in suspensions of nano-sized particles (nanofluids). Int. J. Heat Mass Transfer. 2002, 45, 855–863. [Google Scholar] [CrossRef]
  9. Naseem, T.; Nazir, U.; Sohail, M.; Alrabaiah, H.; Sherif, E.-S.M.; Park, C. Numerical exploration of thermal transport in water-based nanoparticles: A computational strategy. Case Stud. Therm. Eng. 2021, 45, 101334. [Google Scholar] [CrossRef]
  10. Nazir, U.; Nawaz, M.; Alharbi, S.O. Thermal performance of magnetohydrodynamic complex fluid using nano and hybrid nanoparticles. Phys. A Stat. Mech. Its Appl. 2020, 553, 124345. [Google Scholar] [CrossRef]
  11. Koriko, O.K.; Shah, N.A.; Saleem, S.; Chung, J.D.; Omowaye, A.J.; Oreyeni, T. Exploration of bioconvection flow of MHD thixotropic nanofluid past a vertical surface coexisting with both nanoparticles and gyrotactic microorganisms. Sci. Rep. 2021, 11, 16627. [Google Scholar] [CrossRef]
  12. Ali, A.; Saleem, S.; Mumraiz, S.; Saleem, A.; Awais, M.; Marwat, D.N.K. Investigation on TiO2—Cu/H2O hybrid nanofluid with slip conditions in MHD peristaltic flow of Jeffrey material. J. Therm. Anal. Calorim. 2021, 143, 1985–1996. [Google Scholar] [CrossRef]
  13. Tian, M.-W.; Rostami, S.; Aghakhani, S.; Goldanlou, A.S.; Qi, C. A techno-economic investigation of 2D and 3D configurations of fins and their effects on heat sink efficiency of MHD hybrid nanofluid with slip and non-slip flow. Int. J. Mech. Sci. 2021, 189, 105975. [Google Scholar] [CrossRef]
  14. Mumraiz, S.; Ali, A.; Awais, M.; Shutaywi, M.; Shah, Z. Entropy generation in electrical magnetohydrodynamic flow of Al2O3—Cu/H2O hybrid nanofluid with non-uniform heat flux. J. Therm. Anal. Calorim. 2021, 143, 2135–2148. [Google Scholar] [CrossRef]
  15. Awais, M.; Ullah, N.; Ahmad, J.; Sikandar, F.; Ehsan, M.M.; Salehin, S.; Bhuiyan, A.A. Heat transfer and pressure drop performance of Nanofluid: A state-of-the-art review. Int. J. 2021, 9, 100065. [Google Scholar]
  16. Nazir, U.; Sohail, M.; Selim, M.M.; Alrabaiah, H.; Kumam, P. Finite element simulations of hybrid nano-Carreau Yasuda fluid with hall and ion slip forces over rotating heated porous cone. Sci. Rep. 2021, 11, 19604. [Google Scholar] [CrossRef]
  17. Manoj, C.; Kumar, S.; Bhandari, D.R.; Das, K.P.; Mann, I. Development and characterisation of Al2Cu and Al2Al nanoparticle dispersed water and ethylene glycol based nanofluid. Mat. Sci. Eng. 2007, 4, 141–148. [Google Scholar]
  18. Ijaz, M.; Ayub, M.; Malik, M.Y.; Khan, H.; A Alderremy, A.; Aly, S. Entropy analysis in nonlinearly convective flow of the Sisko model in the presence of Joule heating and activation energy: The Buongiorno model. Phys. Scr. 2020, 95, 025402. [Google Scholar] [CrossRef]
  19. Majeed, A.H.; Bilal, S.; Mahmood, R.; Malik, M.Y. Heat transfer analysis of viscous fluid flow between two coaxially rotated disks embedded in permeable media by capitalising non-Fourier heat flux model. Phys. A Stat. Mech. Its Appl. 2020, 540, 123182. [Google Scholar] [CrossRef]
  20. Ali, U.; Rehman, K.U.; Malik, M.Y. Thermal energy statistics for Jeffery fluid flow regime: A generalised Fourier’s law outcomes. Phys. A Stat. Mech. Its Appl. 2020, 542, 123428. [Google Scholar] [CrossRef]
  21. Tanveer, A.; Khan, M.; Salahuddin, T.; Malik, M.; Khan, F. Theoretical investigation of peristaltic activity in MHD based blood flow of non-Newtonian material. Comput. Methods Programs Biomed. 2020, 187, 105225. [Google Scholar] [CrossRef]
  22. Tanveer, A.; Salahuddin, T.; Khan, M.; Malik, M.; Alqarni, M. Theoretical analysis of non-Newtonian blood flow in a microchannel. Comput. Methods Programs Biomed. 2020, 191, 105280. [Google Scholar] [CrossRef]
  23. Khan, M.; Salahuddin, T.; Malik, M.; Alqarni, M.; Alqahtani, A. Numerical modeling and analysis of bioconvection on MHD flow due to an upper paraboloid surface of revolution. Phys. A Stat. Mech. Its Appl. 2020, 553, 124231. [Google Scholar] [CrossRef]
  24. Abbas, N.; Nadeem, S.; Malik, M. On extended version of Yamada–Ota and Xue models in micropolar fluid flow under the region of stagnation point. Phys. A Stat. Mech. Its Appl. 2020, 542, 123512. [Google Scholar] [CrossRef]
  25. Rehman, K.U.; Al-Mdallal, Q.M.; Qaiser, A.; Malik, M.; Ahmed, M. Finite element examination of hydrodynamic forces in grooved channel having two partially heated circular cylinders. Case Stud. Therm. Eng. 2020, 18, 100600. [Google Scholar] [CrossRef]
  26. Zahri, M.; Al-Kouz, W.; Rehman, K.U.; Malik, M.Y. Thermally Magnetised Rectangular Chamber Optimization (TMRCO) of Partially Heated Continuous Stream: Hybrid Meshed Case Study. Case Stud. Therm. Eng. 2020, 22, 100770. [Google Scholar] [CrossRef]
  27. Hayat, T.; Sajjad, R.; Abbas, Z.; Sajid, M.; Hendi, A.A. Radiation effects on MHD flow of Maxwell fluid in a channel with porous medium. Int. J. Heat Mass Transf. 2011, 54, 854–862. [Google Scholar] [CrossRef]
  28. Hayat, T.; Sajjad, R.; Alsaedi, A.; Muhammad, T.; Ellahi, R. On squeezed flow of couple stress nanofluid between two parallel plates. Results Phys. 2017, 7, 553–561. [Google Scholar] [CrossRef]
  29. Saif, R.S.; Hayat, T.; Ellahi, R.; Muhammad, T.; Alsaedi, A. Stagnation-point flow of second grade nanofluid towards a non-linear stretching surface with variable thickness. Results Phys. 2017, 7, 2821–2830. [Google Scholar] [CrossRef]
  30. Hayat, T.; Sajjad, R.; Muhammad, T.; Alsaedi, A.; Ellahi, R. On MHD non-linear stretching flow of Powell–Eyring nanomaterial. Results Phys. 2017, 7, 535–543. [Google Scholar] [CrossRef]
  31. Hayat, T.; Haider, F.; Muhammad, T.; Alsaedi, A. Darcy-Forchheimer flow due to a curved stretching surface with Cattaneo-Christov double diffusion: A numerical study. Results Phys. 2017, 7, 2663–2670. [Google Scholar] [CrossRef]
  32. Hayat, T.; Saif, R.S.; Ellahi, R.; Muhammad, T.; Ahmad, B. Numerical study for Darcy-Forchheimer flow due to a curved stretching surface with Cattaneo-Christov heat flux and homogeneous-heterogeneous reactions. Results Phys. 2017, 7, 2886–2892. [Google Scholar] [CrossRef]
  33. Saif, R.S.; Hayat, T.; Ellahi, R.; Muhammad, T.; Alsaedi, A. Darcy–Forchheimer flow of nanofluid due to a curved stretching surface. Int. J. Numer. Methods Heat Fluid Flow 2019, 29, 2–20. [Google Scholar] [CrossRef]
  34. Hayat, T.; Nawaz, M. Soret and Dufour effects on the mixed convection flow of a second-grade fluid subject to Hall and ion-slip currents. Int. J. Numer. Methods Fluids 2011, 67, 1073–1099. [Google Scholar] [CrossRef]
  35. Nawaz, M.; Alsaedi, A.; Hayat, T.; Alhothauli, M.S. Dufour and Soret effects in an axisymmetric stagnation point flow of second grade fluid with newtonian heating. J. Mech. 2013, 29, 27–34. [Google Scholar] [CrossRef]
  36. Nawaz, M.; Hayat, T.; Alsaedi, A. Dufour and Soret effects on MHD flow of viscous fluid between radially stretching sheets in porous medium. Appl. Math. Mech. 2012, 33, 1403–1418. [Google Scholar] [CrossRef]
  37. Hayat, T.; Nawaz, S.; Alsaedi, A.; Rafiq, M. Mixed convective peristaltic flow of water based nanofluids with Joule heating and convective boundary conditions. PLoS ONE 2016, 11, e0153537. [Google Scholar] [CrossRef] [Green Version]
  38. Naseem, T.; Nazir, U.; Sohail, M. Contribution of Dufour and Soret effects on hydromagnetized material comprising temperature-dependent thermal conductivity. Heat Transf. 2021, 50, 7157–7175. [Google Scholar] [CrossRef]
  39. Naseem, T.; Nazir, U.; El-Zahar, E.R.; Algelany, A.M.; Sohail, M. Numerical Computation of Dufour and Soret Effects on Radiated Material on a Porous Stretching Surface with Temperature-Dependent Thermal Conductivity. Fluids 2021, 6, 196. [Google Scholar] [CrossRef]
  40. Al-Mdallal QMSyam, M.I.; Anwar, M.N. A collocation-shooting method for solving fractional boundary value problems. Commun. Nonlinear Sci. Numer. Simul. 2010, 15, 3814–3822. [Google Scholar] [CrossRef]
  41. Chang, S.H. Numerical solution of Troesch’s problem by simple shooting method. Appl. Math. Comput. 2010, 216, 3303–3306. [Google Scholar] [CrossRef]
  42. Attili, B.S.; Syam, M.I. Efficient shooting method for solving two-point boundary value problems. Chaos Solitons Fractals 2008, 35, 895–903. [Google Scholar] [CrossRef]
  43. Lee, J.; Kim, D.H. An improved shooting method for computation of effectiveness factors in porous catalysts. Chem. Eng. Sci. 2005, 60, 5569–5573. [Google Scholar] [CrossRef]
  44. Nazir, U.; Sohail, M.; Alrabaiah, H.; Selim, M.M.; Thounthong, P.; Park, C. Inclusion of hybrid nanoparticles in hyperbolic tangent material to explore thermal transportation via finite element approach engaging Cattaneo-Christov heat flux. PLoS ONE 2021, 16, e0256302. [Google Scholar] [CrossRef] [PubMed]
  45. Qureshi, I.H.; Nawaz, M.; Abdel-Sattar, M.A.; Aly, S.; Awais, M. Numerical study of heat and mass transfer in MHD flow ofnanofluid in a porous medium with Soret and Dufour effects. Heat Transf. 2021, 50, 4501–4515. [Google Scholar] [CrossRef]
  46. Rana, S.; Nawaz, M.; Alharbi, S.O.; Malik, M.Y. Thermal enhancement in coolant using novel hybrid nanoparticles with mass transport. Case Stud. Therm. Eng. 2021, 28, 101467. [Google Scholar] [CrossRef]
  47. Hafeez, M.B.; Amin, R.; Nisar, K.S.; Jamshed, W.; Abdel-Aty, A.H.; Khashan, M.M. Heat transfer enhancement through nanofluids with applications in automobile radiator. Case Stud. Therm. Eng. 2021, 27, 01192. [Google Scholar] [CrossRef]
Figure 1. The physical arrangement and coordinates system.
Figure 1. The physical arrangement and coordinates system.
Energies 14 06818 g001
Figure 2. Working principle related to finite element method.
Figure 2. Working principle related to finite element method.
Energies 14 06818 g002
Figure 3. Velocity profile for ( G r ) c   when H a 2 = 0.2 ,   P r = 2.36 ,   E c = 0.01 ,   S c = 0.8 ,   S r = 0.1 ,   K = 0.1 ,   B i = 0.05 ,   and   β * = 0.1 .
Figure 3. Velocity profile for ( G r ) c   when H a 2 = 0.2 ,   P r = 2.36 ,   E c = 0.01 ,   S c = 0.8 ,   S r = 0.1 ,   K = 0.1 ,   B i = 0.05 ,   and   β * = 0.1 .
Energies 14 06818 g003
Figure 4. Variation in velocity for ( G r ) t   when M = 0.2 ,   P r = 2.36 ,   E c = 0.01 ,   S c = 0.8 ,   S r = 0.1 ,   K = 0.1 ,   B i = 0.05 , and   β * = 0.1 .
Figure 4. Variation in velocity for ( G r ) t   when M = 0.2 ,   P r = 2.36 ,   E c = 0.01 ,   S c = 0.8 ,   S r = 0.1 ,   K = 0.1 ,   B i = 0.05 , and   β * = 0.1 .
Energies 14 06818 g004
Figure 5. Variation in velocity for K   when H a 2 = 0.2 ,   P r = 2.36 ,   E c = 0.01 ,   S c = 0.8 ,   S r = 0.1 ,   B i = 0.05 ,   and   β * = 0.1 .
Figure 5. Variation in velocity for K   when H a 2 = 0.2 ,   P r = 2.36 ,   E c = 0.01 ,   S c = 0.8 ,   S r = 0.1 ,   B i = 0.05 ,   and   β * = 0.1 .
Energies 14 06818 g005
Figure 6. Variation in velocity for β   when P r = 2.36 ,   E c = 0.01 ,   S c = 0.8 ,   S r = 0.1 ,   B i = 0.05 ,   and   β * = 0.1 .
Figure 6. Variation in velocity for β   when P r = 2.36 ,   E c = 0.01 ,   S c = 0.8 ,   S r = 0.1 ,   B i = 0.05 ,   and   β * = 0.1 .
Energies 14 06818 g006
Figure 7. Variation in velocity for M when P r = 2.36 ,   E c = 0.01 ,   S c = 0.8 ,   S r = 0.1 ,   B i = 0.05 ,   and   β * = 0.1 .
Figure 7. Variation in velocity for M when P r = 2.36 ,   E c = 0.01 ,   S c = 0.8 ,   S r = 0.1 ,   B i = 0.05 ,   and   β * = 0.1 .
Energies 14 06818 g007
Figure 8. Variation in temperature for M when P r = 2.36 ,   E c = 0.01 ,   S c = 0.8 ,   ( G r ) c = 0.1 ,   K = 0.1 ,   B i = 0.05 ,   β * = 0.1 ,   and   S r = 0.1 .
Figure 8. Variation in temperature for M when P r = 2.36 ,   E c = 0.01 ,   S c = 0.8 ,   ( G r ) c = 0.1 ,   K = 0.1 ,   B i = 0.05 ,   β * = 0.1 ,   and   S r = 0.1 .
Energies 14 06818 g008
Figure 9. Variation in temperature for E c when M = 0.2 ,   P r = 2.36 ,   S r = 0.1 ,   S c = 0.8 ,   ( G r ) c = 0.1 ,   K = 0.1 ,   B i = 0.05 ,   and   β * = 0.1 .
Figure 9. Variation in temperature for E c when M = 0.2 ,   P r = 2.36 ,   S r = 0.1 ,   S c = 0.8 ,   ( G r ) c = 0.1 ,   K = 0.1 ,   B i = 0.05 ,   and   β * = 0.1 .
Energies 14 06818 g009
Figure 10. Variation in temperature for β when M = 0.2 ,   P r = 2.36 ,   S r = 0.1 ,   S c = 0.8 ,   ( G r ) c = 0.1 ,   K = 0.1 ,   B i = 0.05 ,   and   β * = 0.1 .
Figure 10. Variation in temperature for β when M = 0.2 ,   P r = 2.36 ,   S r = 0.1 ,   S c = 0.8 ,   ( G r ) c = 0.1 ,   K = 0.1 ,   B i = 0.05 ,   and   β * = 0.1 .
Energies 14 06818 g010
Figure 11. Variation in temperature for B i when M = 0.2 ,   P r = 2.36 ,   S r = 0.1 ,   S c = 0.8 ,   ( G r ) c = 0.1 ,   K = 0.1 ,   and   β * = 0.1 .
Figure 11. Variation in temperature for B i when M = 0.2 ,   P r = 2.36 ,   S r = 0.1 ,   S c = 0.8 ,   ( G r ) c = 0.1 ,   K = 0.1 ,   and   β * = 0.1 .
Energies 14 06818 g011
Figure 12. Variation in temperature for P r when M = 0.2 ,   S r = 0.1 ,   S c = 0.8 ,   ( G r ) c = 0.1 ,   K = 0.1 ,   and   β * = 0.1 .
Figure 12. Variation in temperature for P r when M = 0.2 ,   S r = 0.1 ,   S c = 0.8 ,   ( G r ) c = 0.1 ,   K = 0.1 ,   and   β * = 0.1 .
Energies 14 06818 g012
Figure 13. Variation in concentration for S r when M = 0.3 ,   P r = 2.5 ,   E c = 0.3 ,   S c = 0.3 ,   ( G r ) c = 0.1 ,   K = 0.1 ,   B i = 0.3 ,   and   β * = 0.5 .
Figure 13. Variation in concentration for S r when M = 0.3 ,   P r = 2.5 ,   E c = 0.3 ,   S c = 0.3 ,   ( G r ) c = 0.1 ,   K = 0.1 ,   B i = 0.3 ,   and   β * = 0.5 .
Energies 14 06818 g013
Figure 14. Variation in concentration for S c when M = 0.2 ,   P r = 2.36 ,   S r = 0.1 ,   ( G r ) c = 0.1 ,   K = 0.1 ,   B i = 0.05 ,   and   β * = 0.1 .
Figure 14. Variation in concentration for S c when M = 0.2 ,   P r = 2.36 ,   S r = 0.1 ,   ( G r ) c = 0.1 ,   K = 0.1 ,   B i = 0.05 ,   and   β * = 0.1 .
Energies 14 06818 g014
Figure 15. Variation in concentration for E i when M = 0.2 ,   P r = 2.36 ,   S r = 0.1 ,   S c = 0.8 ,   ( G r ) c = 0.1 ,   K = 0.1 ,   and   β * = 0.1 .
Figure 15. Variation in concentration for E i when M = 0.2 ,   P r = 2.36 ,   S r = 0.1 ,   S c = 0.8 ,   ( G r ) c = 0.1 ,   K = 0.1 ,   and   β * = 0.1 .
Energies 14 06818 g015
Figure 16. Variation in concentration for α when M = 0.2 ,   P r = 2.36 ,   S r = 0.1 ,   S c = 0.8 ,   ( G r ) c = 0.1 ,   K = 0.1 ,   and   β * = 0.1 .
Figure 16. Variation in concentration for α when M = 0.2 ,   P r = 2.36 ,   S r = 0.1 ,   S c = 0.8 ,   ( G r ) c = 0.1 ,   K = 0.1 ,   and   β * = 0.1 .
Energies 14 06818 g016
Table 1. Values of properties of nano-structures and hybrid nano-structures.
Table 1. Values of properties of nano-structures and hybrid nano-structures.
Physical PropertyEthylene Glycol (EG)MOS2MOS2/SiO2
ρ ( K g / m 3 ) 1113.526505060
k ( W / ( m · K ) ) 2430730397.746
c p ( J / ( k g · K ) ) 0.2531.534.5
σ / ( Ω / m ) 4.3 × 10 5 0.0005 1 × 10 18
Table 2. Mesh-free analysis of f ( η m a x 2 ) ,   θ ( η m a x 2 ) and ϕ ( η m a x 2 ) at the mid of the computational domain [ 0 ,   6 ] .
Table 2. Mesh-free analysis of f ( η m a x 2 ) ,   θ ( η m a x 2 ) and ϕ ( η m a x 2 ) at the mid of the computational domain [ 0 ,   6 ] .
Number of Elements f ( η m a x 2 ) θ ( η m a x 2 ) ϕ ( η m a x 2 )
300.0080898667630.015037047580.02481424434
600.0086925213440.0087312239850.01864316747
900.0088678524790.0063777687770.01586843723
1200.0089521091390.0051110484920.01419676870
1500.0090019590020.0043076053870.01304491414
1800.0090349383700.0037473986440.01218677926
2100.0090584985100.0033318981050.01151399434
2400.0090764118660.0030099495250.01096716142
2700.0090343201360.0030099495250.01096716142
3000.0090943050310.0030103482130.01092535030
Table 3. Comparative study of a temperature gradient with published results of Qureshi et al. [45] in the case of nanofluid.
Table 3. Comparative study of a temperature gradient with published results of Qureshi et al. [45] in the case of nanofluid.
[45]present results
Nusselt numberNusselt number
0.680.681052103137
k n f 0.721410.723331807103
0.824580.824720819103
Table 4. Comparative numerical values of Nusselt number with shooting method for nanofluid.
Table 4. Comparative numerical values of Nusselt number with shooting method for nanofluid.
Shooting Method [47]Finite Element Method
N u N u
0.93910.93890213031
k n f 1.066961.06691013039
Table 5. Local skin friction coefficient, Nusselt number and Sherwood number for hybrid nanofluid ( M O S 2 S i O 2 ) and nanofluid ( M O S 2 ) when M = 0.2 ,     P r = 2.36 ,     S r = 0.1 ,     S c = 0.8 ,     ( G r ) c = 0.1 ,     K = 0.1 ,   and   β * = 0.1 .
Table 5. Local skin friction coefficient, Nusselt number and Sherwood number for hybrid nanofluid ( M O S 2 S i O 2 ) and nanofluid ( M O S 2 ) when M = 0.2 ,     P r = 2.36 ,     S r = 0.1 ,     S c = 0.8 ,     ( G r ) c = 0.1 ,     K = 0.1 ,   and   β * = 0.1 .
Hybrid   Nanofluid   ( M O S 2 S i O 2 ) Nanofluid   ( M O S 2 )
R e 1 2 C f R e 1 2 N u R e 1 2 S h R e 1 2 C f R e 1 2 N u R e 1 2 S h
0.20.7972750.5848090.4716310.3662150.5155420.392980
E i 0.30.7814520.5829140.5539710.3614980.5136250.437901
0.40.7750780.5821430.5873030.3588240.5125340.463449
0.60.8130440.6960150.7028770.3594940.6121300.591962
0.40.8313210.5656370.3600100.7476940.6233370.346218
B i 0.50.8293790.7824240.3828620.7370700.6458480.344631
0.60.8106320.6028740.3574820.7295930.6618060.343501
0.80.7993290.6234480.3560740.7464630.8322620.372349
0.01.0976110.5643460.3500030.9950860.6189570.336153
( G r ) t 0.10.9543710.5767800.3546130.8621830.6335710.340746
0.20.8191970.5873800.3585360.7370700.6458480.344631
0.30.6903220.5966710.3619610.6179280.6564700.348014
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hafeez, M.B.; Sumelka, W.; Nazir, U.; Ahmad, H.; Askar, S. Mechanism of Solute and Thermal Characteristics in a Casson Hybrid Nanofluid Based with Ethylene Glycol Influenced by Soret and Dufour Effects. Energies 2021, 14, 6818. https://doi.org/10.3390/en14206818

AMA Style

Hafeez MB, Sumelka W, Nazir U, Ahmad H, Askar S. Mechanism of Solute and Thermal Characteristics in a Casson Hybrid Nanofluid Based with Ethylene Glycol Influenced by Soret and Dufour Effects. Energies. 2021; 14(20):6818. https://doi.org/10.3390/en14206818

Chicago/Turabian Style

Hafeez, Muhammad Bilal, Wojciech Sumelka, Umar Nazir, Hijaz Ahmad, and Sameh Askar. 2021. "Mechanism of Solute and Thermal Characteristics in a Casson Hybrid Nanofluid Based with Ethylene Glycol Influenced by Soret and Dufour Effects" Energies 14, no. 20: 6818. https://doi.org/10.3390/en14206818

APA Style

Hafeez, M. B., Sumelka, W., Nazir, U., Ahmad, H., & Askar, S. (2021). Mechanism of Solute and Thermal Characteristics in a Casson Hybrid Nanofluid Based with Ethylene Glycol Influenced by Soret and Dufour Effects. Energies, 14(20), 6818. https://doi.org/10.3390/en14206818

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