Next Article in Journal
A New Incommensurate Fractional-Order COVID 19: Modelling and Dynamical Analysis
Next Article in Special Issue
Taylor DQN: An Optimization Method for Aircraft Engine Cleaning Schedule
Previous Article in Journal
Change Point Analysis for Kumaraswamy Distribution
Previous Article in Special Issue
Study of Steam-Induced Convection in a Rotating Vertical Flow Channel
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Simultaneous Features of CC Heat Flux on Dusty Ternary Nanofluid (Graphene + Tungsten Oxide + Zirconium Oxide) through a Magnetic Field with Slippery Condition

1
Department of Physics, College of Science, King Faisal University, P.O. Box 400, Al-Ahsa 31982, Saudi Arabia
2
Laboratory of Fluid Mechanics, Physics Department, Faculty of Sciences of Tunis, University of Tunis El Manar, Tunis 2092, Tunisia
Mathematics 2023, 11(3), 554; https://doi.org/10.3390/math11030554
Submission received: 20 December 2022 / Revised: 15 January 2023 / Accepted: 18 January 2023 / Published: 20 January 2023
(This article belongs to the Special Issue Mathematical Methods on Mechanical Engineering)

Abstract

:
The purpose of this work is to offer a unique theoretical ternary nanofluid (graphene/tungsten oxide/zirconium oxide) framework for better heat transfer. This model describes how to create better heat conduction than a hybrid nanofluid. Three different nanostructures with different chemical and physical bonds are suspended in water to create the ternary nanofluid (graphene/tungsten oxide/zirconium oxide). Toxic substances are broken down, the air is purified, and other devices are cooled thanks to the synergy of these nanoparticles. The properties of ternary nanofluids are discussed in this article, including their thermal conductivity, specific heat capacitance, viscosity, and density. In addition, heat transport phenomena are explained by the Cattaneo–Christov (CC) heat flow theory. In the modeling of the physical phenomena under investigation, the impacts of thermal nonlinear radiation and velocity slip are considered. By using the right transformations, flow-generating PDEs are converted into nonlinear ordinary differential equations. The parameters’ impacts on the velocity and temperature fields are analyzed in detail. The modeled problem is graphically handled in MATLAB using a numerical technique (BVP4c). Graphical representations of the important factors affecting temperature and velocity fields are illustrated through graphs. The findings disclose that the performance of ternary nanofluid phase heat transfer is improved compared to dusty phase performance. Furthermore, the magnetic parameter and the velocity slip parameter both experience a slowing-down effect of their respective velocities.

1. Introduction

Nanofluids are liquids that include ultrafine particles contained in a dilute solution. These particles are suspended in the liquid. These fluids have thermophysical qualities that are far better than those of pure fluids. It has been established that increasing the amount of carbon, copper, or other nanoparticles with high thermal conductivity that is added to water, ethylene glycol, or oil may enhance the materials’ ability to transport heat. The thermal conductivity of the nanofluids is significantly raised, and the particle size and volume fraction of suspended particles have an impact on this phenomenon. Choi [1] invented the term “nanofluids”, which defines a liquid suspension of ultrasmall particles (less than 50 nm in diameter). With the fast development of nanofabrication, various affordable combinations of fluid/particles are instantly accessible. A complete assessment by Buongiorno [2] investigated convective transport in nanofluids, who argued that an appropriate explanation for the phenomenon exists, although unexpected growth in heat conductivity and viscosity is yet to be found. Recently, Reddy et al. [3] investigated hybrid dusty fluid flow through a Cattaneo–Christov heat flux model. Waqas et al. [4] discussed the heat transfer analysis of hybrid nanofluid flow with thermal radiation through a stretching sheet. Gurdal et al. [5] studied the turbulent flow and heat transfer characteristics of Ferro-nanofluid flowing in a dimpled tube under a magnetic field effect. Many researchers [6,7,8,9] have addressed the convective heat transfer during the flow of nanoliquid.
The passage of heat energy through the particles of a liquid may be described as a sort of heat transfer known as thermal radiation. It is a condition that must be met in a wide variety of applications, such as the planning of propulsion devices for space spacecraft, nuclear power plants, and gas turbines. In space applications, where devices are required to operate at high temperatures to obtain the required thermal efficiency, the impact of thermal radiation on controlling the heat transfer in certain applications and on estimating the thermal effects for processes with high temperatures is unavoidable. This is because thermal radiation has an impact on both controlling the heat transfer in certain applications and on estimating the thermal effects for processes with high temperatures. Duwairi and Duwairi [10] studied gray viscous fluid flow with magnetic and radiation parameters influences. Cortell [11] analyzed radiation and viscous dissipation impacts for the flow of the thermal boundary layer along a nonlinear stretching surface. Recently, Bakar and Soid [12] studied the MHD stagnation-point flow and heat transfer over an exponentially stretching/shrinking vertical sheet in a micropolar fluid with a buoyancy effect. Ali et al. [13] analyzed the flow and heat transfer over stretching/shrinking and porous surfaces. Azam et al. [14] discussed the transient bioconvection and activation energy impacts on Casson nanofluid with gyrotactic microorganisms and nonlinear radiation. Reddy et al. [15] studied the unsteady absorption flow and dissipation heat transfer over a non-Newtonian fluid. Gnaneswara et al. [16] investigated the effect of thermal conductivity on Blasius–Rayleigh–Stokes flow and heat transfer over a moving plate by considering the magnetic dipole moment. Some recent studies related to the consideration of heat transfer theory can be found in [17,18,19,20,21,22].
Fourier’s law of heat conduction describes heat transfer in a given medium due to temperature differences [23]. It is considered the basis for the investigation of thermal phenomena. Using this law, the temperature field equation is developed into a parabolic-type equation stating that heat transfer with initial disruption having infinite speed propagates throughout the medium under certain conditions. This heat conduction paradox requires modification of Fourier’s law. Cattaneo [24] presented a modified form of Fourier’s law known as the Maxwell–Cattaneo (MC) model by introducing the thermal relaxation time parameter. This modification converted the temperature field equation from a parabolic- to a hyperbolic-type equation, describing that heat transfer has finite speed in the entire medium. Christov [25] investigated the frame invariance of the MC model and replaced the time derivative with Oldroyd’s upper-convected time derivative. Recently, Reddy et al. [26] studied zero-mass flux and Cattaneo–Christov heat flux through a Prandtl non-Newtonian nanofluid in Darcy–Forchheimer porous space. Machireddy et al. [27] investigated the impact of Cattaneo–Christov heat flux on the hydromagnetic flow of non-Newtonian fluids filled with a Darcy–Forchheimer porous medium. Khan et al. [28] discussed Cattaneo–Christov (CC) heat and mass fluxes in the stagnation-point flow of Jeffrey nanoliquids by a stretching surface. Tausif et al. [29] investigated the modified homogeneous and heterogeneous chemical reaction and flow performance of Maxwell nanofluid with the Cattaneo–Christov heat flux law. Some recent studies related to the consideration of CC theory can be found in [30,31,32].
In this study, we investigated dusty fluid flow subjected to a magnetic field over a stretching sheet. Impacts of thermal radiation, velocity slip, and conditions are considered in the modeling of the physical phenomena. The use of CC heat flux theory is carried out in the energy equation for analyzing the heat transfer phenomena. Considering the appropriate variables, the governing PDEs are converted to nonlinear ODEs. MATLAB software is used for investigating the numerical solution. Impacts of fluid velocity and temperature are analyzed through different pertinent parameters in the considered problem.
Specifically, this work will address the research questions mentioned below:
  • What will be the influence of the magnetic parameter on fluid heat transfer rate, temperature, and velocity?
  • What is the significance of the CC theory on the dusty ternary fluid model?
  • What will be the impact of thermal radiation on heat transfer rate and temperature?
  • What will be the influence of the velocity slip parameter on fluid velocity, temperature, and heat transfer rate?
We organized this article as follows. In Section 2, we formulate the physical phenomena of the problem and transformations of PDEs. Result analyses are provided in Section 3. Concluding remarks of the present article are given in Section 4.

2. Mathematical Formulation

Consider a continuous 2D boundary layer flow and the heat transfer of a dusty ternary nanofluid across a stretching sheet moving at a constant speed U w = b x . With the slot serving as the point of origin, the x -axis is drawn along the surface that is being stretched in the direction that the motion is going, and the y -axis is drawn perpendicular to the sheet in the direction that is facing away from the fluid. It is presumable that the flow will not extend beyond the area where y is greater than 0. Convective heat transfer is responsible for keeping the temperature of the sheet’s surface at a constant value known as T f .
Some of the assumptions taken into consideration are listed below:
  • The dust particles are taken to be small enough and of sufficient number to be treated as a continuum and allow concepts such as density and velocity to have physical meaning.
  • The dust particles are assumed to be spherical in shape, having the same radius and mass, and are undeformable.
  • The Cartesian coordinate system is located in such a way that the x -axis and y -axis are taken along (and normal to) the surface, respectively, while the origin of the system is located at the leading edge.
  • The dust particles are assumed to be uniform in size, and the density number of the dust particle is taken as constant throughout the flow.
The equations that regulate the flow of dusty ternary nanofluid can be written as follows, according to the standard boundary layer approximations (Reddy et al. [3]):
u x + v y = 0
u u x + v u y = μ t h n f ρ t h n f 2 u y 2 σ B 0 2 ρ t h n f u
( ρ c p ) t h n f ( u T x + v T y ) + λ 1 ( u 2 2 T x 2 + v 2 2 T y 2 + 2 u v 2 T x y + ( u u x + v u y ) T x + ( u v x + v v y ) T y ) = k t h n f 2 T y 2 + N c p f τ t ( T p T ) + N τ v ( u p u ) 2 + 1 ρ c p q r y + μ t h n f ( u y ) 2
Dust Phase
u p x + v p y = 0
u p u p x + v p u p y = K m ( u u p )
( u p T p x + v p T p y ) = c t h p f c t h m f τ T ( T T P )
By using the Rosseland approximation to the radiative heat flow expression found in Equation (3), one obtains the following:
q r = 4 σ * 3 k * T 4 y = 16 σ * 3 k * T 3 T y
where x and y , respectively, represent coordinate axes along the continuous surface in the direction of motion and perpendicular to it. ( u ,   v ) and ( u p ,   v p ) denote the velocity components of the nanofluid and dust phases along the x - and y -directions, respectively. N is the number density of dust particles. λ 1 is the thermal relaxation time. m is the mass concentration of dust particles. K = 6 π μ n f   r is the Stokes drag constant. r is the radius of dust particles. A is the slip constant. σ * is the Stefan–Boltzmann constant. k * is the mean absorption coefficient.
Boundary Conditions:
u = U w + A ν f u y ,       v = 0 ,     k f T y = h f ( T f T )         at     y = 0 , u 0 ,             u p 0 ,         v p v ,     T T ,     T p T         a s               y ,
The expression for viscosity ( μ t h n f ), density ( ρ t h n f ), specific heat ( ( ρ C p ) t h n f ), and thermal conductivity ( k t h n f ) of the ternary nanofluid is as follows (Prakash et al. [9]):
μ t h n f = μ f ( 1 ϕ 1   ) 2.5 ( 1 ϕ 2   ) 2.5 ( 1 ϕ 3   ) 2.5 ,   ( ρ C p ) t h n f = [ ( 1 ϕ 2 ) [ ( 1 ϕ 1 ) ( ρ C p ) f + ϕ 1 ( ρ C p ) s 1 ] + ϕ 2 ( ρ C p ) s 2 ] ,   ρ t h n f = ( 1 ϕ 1 ) [ ( 1 ϕ 2 ) { ( 1 ϕ 3 ) ρ f + ϕ 3 ρ 3 } + ϕ 2 ρ 2 ] + ϕ 1 ρ 1 ,   k t h n f k h n f = k 1 + 2 k h n f 2 ϕ 1 ( k h n f k 1 ) k 1 + 2 k h n f + ϕ 1 ( k h n f k 1 ) , k h n f k n f = k 2 + 2 k n f 2 ϕ 2 ( k n f k 2 ) k 2 + 2 k n f + ϕ 2 ( k h n f k 2 ) , k n f k n f = k 3 + 2 k f 2 ϕ 3 ( k f k 3 ) k 3 + 2 k n f + ϕ 3 ( k f k 3 )
where ϕ 1 ,     ϕ 2 , and ϕ 3 are the nanoparticle volume fraction of copper/alumina/zirconium oxide, respectively. C p is the specific heat. k f denotes the thermal conductivity of the regular fluid.
The nonlinear PDEs in the model problem are converted to nonlinear ODEs by way of a similarity transformation, as shown below.
u = b x f ( η ) , v = b v f f ( η ) , u p = b x F ( η ) , v p = b v f F ( η ) η = b v f y , θ ( η ) = T T T f T ,       θ p ( η ) = T p T T f T
with T = T f ( 1 + ( θ w 1 ) θ )   and   θ w = T f T , θ w > 1 the temperature ratio parameter.
Equations (1) and (4) are completely satisfied by these transformations, and the other equations in the model, [2,3,5,6], together with the boundary conditions, are translated as follows:
Fluid Phase:
μ t h n f ρ t h n f f + ( f f f 2 ) + μ t h n f [ l β v ( F f ) Q F ] = 0 .
k t h n f k f ( θ + R [ ( 1 + ( θ w 1 ) θ ) 3 θ + 3 ( θ w 1 ) θ 2 ( 1 + ( θ w 1 ) θ ) 2 ] ) + Pr ( ρ c p ) f ( ρ c p ) t h n f [ Γ ( f f θ f 2 θ ) + f θ + l β T [ θ p θ ] + l E c β v [ F f ] 2 ] + P r E c f 2 = 0 .
Dust Phase
F F F 2 + β υ ( f F ) = 0
θ p F + α β T ( θ p θ ) = 0
f ( 0 ) = 0 ,     f ( 0 ) = 1 + δ ( 1 + ϕ 1 ) 2.5 ( 1 + ϕ 2 ) 2.5 ( 1 + ϕ 2 ) 2.5 f ( 0 ) ,     θ ( 0 ) = B i ( 1 θ ( 0 ) ) ,   f ( ) 0 ,     F ( ) 0 ,     F ( ) f ( ) ,     θ ( ) 0 ,     θ p ( ) 0 .
where Q = σ B 0 2 ρ f , Pr = ( μ c p ) f k f , E c = U w 2 ( T f T ) c p f , R = 4 σ * T 3 3 k t h n f k * ,   β T = 1 τ T c , β v = 1 τ v c , τ v = m K , l = m N ρ f , γ = c p f c m f   Γ = b λ 1 , δ = A U w ν f and B i = ν f   c h f k f .
The physical quantities of C f and N u are;
C f = τ w ρ f U w 2 ,   N u = x q w k f ( T f T )   at   y = 0 ,
where τ w and q w are given by
τ w = μ t h n f ( u y ) q w = k t h n f ( T y )   at   y = 0 .
In terms of the nondimensional variables, we express them as
C f R e ( 1 2 ) = ( μ t h n f μ f ) f ( 0 )   and     N u ( R e ) ( 1 2 ) = k t h n f k f [ 1 + R θ w 3 ] θ ( 0 ) .
where Re = U w 2 b v f is the local Reynolds number.

3. Numerical Method

Solving the nonlinear system of ODEs generated by a numerical technique (BVP4c) in MATLAB (see Equations (10)–(13)) may be achieved numerically with the aid of boundary constraints using the approach described in Equation (14). Table 1 shows the physical properties of the ternary nanofluid. Equations of a higher level of order differential type are produced. A variety of physical influences on the velocity and temperature profiles for both dusty and ternary nanofluid cases are analyzed graphically. Calculated values for C f ( R e ) 1 2 , and N u ( R e ) 1 2 are tabulated below (Tables 3 and 4). At the outset, we reduce these ODEs to the first-order form by simplifying the higher-order terms.
Let f = y 1 ,   f = y 2 ,   f = y 3 ,   f = y 4 ,   θ = y 5 ,   θ = y 6 ,   θ = y 7   F = y 8 ,   F = y 9 ,   F = y 10 ,   θ p = y 11 ,   θ p = y 12 .
Fluid Phase:
y 4 = ρ t h n f μ t h n f [ ( y 1 y 3 y 2 2 ) + μ t h n f [ l β v ( y 9 y 2 ) Q y 2 ] ]   ,
y 7 = k f k t h n f [ ( 1 + R [ ( 1 + ( θ w 1 ) y 5 ) 3 + 3 ( θ w 1 ) y 6 2 ( 1 + ( θ w 1 ) y 5 ) 2 ] ) + Pr ( ρ c p ) f ( ρ c p ) t h n f [ Γ ( y 1 y 2 y 6 y 1 2 y 6 ) + y 1 y 6 + l β T [ y 11 y 5 ] + l E c β [ y 9 y 2 ] 2 ] + P r E c y 3 2 . ] ,
Dust Phase:
y 10 = 1 y 8 [ y 9 2 β v ( y 2 y 9 ) ]
y 12 = 1 y 8 [ α β T ( y 11 y 5 ) ]
Boundary Conditions:
( 0 ) = 0 ,     y 2 ( 0 ) = 1 + δ ( 1 + ) 2.5 y 3 ( 0 ) ,     y 6 ( 0 ) = B i ( 1 y 5 ( 0 ) ) ,     y 2 ( ) 0 ,     y 9 ( ) 0 ,     y 8 ( ) y 1 ( ) ,     y   5 ( ) 0 ,     y 11 ( ) 0

4. Results and Discussion

In this paper, our study is focused on the dusty ternary nanofluid flow with the Cattaneo–Christov model over a stretching sheet in the presence of a magnetic field, nonlinear radiation, and convective heat transfer. We analyzed the effects of the physical parameters on the velocity and thermal fields through figures and tables by using the fixed parametric values Pr = 6.9 ,   E c = 0.5 ,   R = 0.6 ,     θ w = 1.2 ,   δ = 0.7 ,   Q = 0.5 ,   B i = 0.8 ,   ϕ 1 = ϕ 2 = ϕ 3 = 0.005   , and Γ = 0.4 . Additionally, the results of the comparison of the Prandtl number to previously published studies are shown in Table 2. We can conclude that the current results are in excellent agreement with previous findings in the literature. Furthermore, C f ( R e ) 1 2 and N u ( R e ) 1 2 values for the present analysis are computed and given in Table 3 and Table 4.
Figure 1 shows how Q affects the velocity field ( f ( η ) and F ( η ) ) for both ternary nanofluid and dusty fluid cases. A decreasing trend can be observed in f ( η ) and F ( η ) scenarios for boosting values of Q . Physically, the magnetic field depends on Lorentz force, which is stronger for a larger magnetic field. Therefore, the velocity profile declines with more magnetic fields. In addition, the influence of the ternary nanofluid is examined more thoroughly than in the case of dusty fluid. In addition, the appropriate layer thickness is scaled back in both scenarios when the Q value is increased.
Figure 2, Figure 3 and Figure 4 depict the instances of dusty fluid and ternary nanofluid and the influence of ϕ 1 , ϕ 2 , and ϕ 3 on the F ( η ) and f ( η ) profiles, respectively. A decrease in the velocity curve that occurs as a consequence of the incorporation of solid nanoparticles ( ϕ 1 , ϕ 2 , and ϕ 3 ) into the nanomaterial was seen in both dusty and ternary nanofluid cases. This outcome is due to collisions between nanoparticles that have been extensively dispersed across the environment.
In order to explore the impact of the δ parameter on the velocity distributions of the F ( η ) and f ( η ) profiles, Figure 5 is displayed. It is clear from these results that changing the value of the parameter results has a significant reduction of the velocity distribution for both the F ( η ) and f ( η ) situations. This is due to the fact that the fluid’s near-surface velocity does not sustain correctly for the extended surface velocity when taking slip effect into consideration. Therefore, if the value of the slip velocity parameter is increased, the slide speed will be increased accordingly. Following this, there is a decrease in the liquid velocity as a result of the stretching surface deformation, which can only be communicated to the fluid under slip conditions.
In Figure 6, we can see how R affects the temperature field for both ternary nanofluid phase and dusty phase ( θ ( η ) and θ p ( η ) ) profiles. It seems to be enhanced when R improves. In accordance with the physical principle, as the radiation parameter increases, heat is radiated at the fluid at a much higher rate than when the mean absorption coefficient decreases. Physically, the effect of thermal radiation is to boost the transfer, as by increasing thermal radiation, the thermal boundary layer increases. Therefore, it has been reported that the process of thermal radiation reduction should proceed at a faster rate. The use of radiation can regulate the distribution of temperature and flow, and such applications can be used in pseudoscience to monitor blood pressure via the magnetotherapy process.
Figure 7 depicts the analysis of the characteristics of the function θ ( η ) and θ p ( η ) in relation to rising values of θ w for a ternary nanofluid and dusty fluid condition, respectively. It is witnessed that the thermal field and thermal layer thickness of the ternary nanofluid and dusty fluid condition both increase at an exponential rate when the θ w parameter is increased. In addition, for high values of θ w , a rise in fluid temperature enriches the whole fluid case. This is the situation when the temperature of the fluid is increased. In addition, the enhancement of fluid temperature is optimal in the case of the ternary nanofluid, followed by the case of dusty fluid.
Figure 8 displays the changes in the thermal field for the ternary nanofluid and dusty fluid instances as a function of the E c parameter. To improve the E c parameter, it is evident that the thermal field rises in both ternary nanofluid and dusty fluid phases. Physically, the Eckert number is the ratio of kinetic energy to the specific enthalpy difference between wall and fluid. Therefore, an increase in the Eckert number causes the transformation of kinetic energy into internal energy by work that is done against the viscous fluid stresses. Due to this, increasing E c enhances dramatically the temperature of the fluid.
In Figure 9, we see how varying the Γ parameter affects the thermal distribution for both θ ( η ) and θ p ( η ) cases. It is inferred that when the Γ parameter increases, the fluid temperature decreases. An increase in the Γ parameter decreases the temperature field as it decelerates the relaxation of heat flow.
The thermal field for the θ ( η ) and θ p ( η ) instances is shown to be affected by the ϕ 1 , ϕ 2 , and ϕ 3 parameters in Figure 10, Figure 11 and Figure 12. The fact that the nanoparticle is so tiny may be deduced from the amount of heat that it generates when it releases the accumulated energy. It is possible that in both cases ( θ ( η ) and θ p ( η ) ), extra energy will be needed for the mixing of additional nanoparticles. In turn, this will cause an increment in the temperature and the thickness of the associated layer.
Figure 13 shows the effect of the B i parameter on thermal field for both θ ( η ) and θ p ( η ) instances. It is noted that the thermal field for both θ ( η ) and θ p ( η ) instances increases with an increase in the B i parameter. Furthermore, the interrelated thickness of the thermal layer also enhances by rising the values of the B i parameter.
In Table 3 and Table 4, the impact of various substantial factors over skin friction and the Nusselt number for the fluid flow system has been portrayed in numerical manners. The fluid motion decreases with expansion in magnetic effects and slip factor. The friction factor decreases with augmenting values of the nanoparticle volume friction parameter ( ϕ 1 ,   ϕ 2 , and ϕ 3 ), as shown in Table 3. The Nusselt number is maximized for higher values of the nanoparticle volume friction parameter ( ϕ 1 ,     ϕ 2 , and ϕ 3 ), as shown in Table 4. An increment of the values of the Eckert number, Prandtl number, and Biot number enhances the Nusselt number. Additionally, the Nusselt number rises with the rising values of thermal relaxation. One can conclude that the Nusselt number is enhanced for higher values of radiation and temperature parameters.

5. Conclusions

By using a CC heat flux model, we investigated the movement of a dusty ternary nanofluid over a stretching sheet. The shown partial differential equations are converted to ordinary differential equations and numerically solved via the BVP4c technique. Graphs depict the concentration, the velocity, and the temperature. Overall, we are able to draw the following conclusions:
  • The magnetic field depends essentially on Lorentz force, which is predominant for a larger magnetic field. Due to this, the velocity profiles decrease by further increasing the values of the magnetic parameter.
  • The velocity profiles are decreasing by increasing the values of the slip parameter. However, once the slip velocity parameter is increased, the slide speed will also be decreased accordingly.
  • The Eckert number is defined as the ratio between kinetic energy and the specific enthalpy difference between the wall and the fluid. Thus, increasing the Eckert number causes an increment in the temperature profile for both dusty and ternary phases.
  • The thermal radiation role is to boost heat transfer by enlarging the thermal radiation and thermal boundary layers for both dusty and ternary phases.
  • The tiny nanoparticle effect may be deduced from the amount of heat generated by the accumulated energy. Due to this, temperature profiles are enhanced by increasing values of the nanoparticle volume friction parameter.
  • An increment in the value of thermal relaxation generates a decrement in the temperature distribution.
  • Temperature profiles increase by increasing the values of the temperature ratio parameter.
  • The performance of the dusty phase heat transfer has significantly enhanced compared to the performance of the ternary nanofluid phase.

Funding

This work was supported by the Deanship of Scientific Research, Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia (Grant No. 2256).

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

This work was supported by the Deanship of Scientific Research, Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia (Grant No. 2256).

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

B 0 2 Magnetic field
c p f Specific heat of fluid particle
c m f Specific heat of dust particle
C f Friction factor
E c Eckert number
K Stokes drag coefficient
k t h n f Thermal conductivity of the nano fluid
k f Thermal conductivity of the fluid
k * Mean absorption coefficient
l Mass concentration of particles
m Mass concentration of dust particles
Q Magnetic parameter
N Number density of dust particle
N u Nusselt number
Pr Prandtl number
q r Radiative heat flux
R Radiation parameter
u , u p Velocity   along   x   axis   ( m . s 1 )
v , v p Velocity   along   y   axis   ( m . s 1 )
B i Biot number
δ Slip parameter
T Ambient fluid temperature
T w Wall temperature
T , T p Temperature of fluid and dust particle
Greek letters
α Ratio of specific heat
β v Fluid particle interaction parameter for velocity
β t Fluid particle interaction parameter for temperature
γ Thermal relaxation parameter
τ t Relaxation time of dust particle for temperature
τ v Relaxation time of dust particle for velocity
λ 1 Relaxation time for heat flux
σ Electrical conductivity
ϕ 1 ,   ϕ 2 , ϕ 3 Volume fraction of nano particle
σ * Stefan Boltzmann constant
( ρ c p ) Heat capacity of the fluid
( ρ C p ) t h n f Heat capacity of the nanofluid
k t h n f Thermal conductivity of ternary hybrid fluid
k f Thermal conductivity of fluid

References

  1. Choi, U.; Eastman, J. 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. Buongiorno, J. Convective transport in nano fluids. J. Heat Transf. 2006, 128, 240–250. [Google Scholar] [CrossRef]
  3. Reddy, M.; Rani, M.; Kumar, K.; Prasannakumar, B.; Lokesh, H. Hybrid dusty fluid flow through a Cattaneo–Christov heat flux model. Phys. A Stat. Mech. Its Appl. 2020, 551, 123975. [Google Scholar] [CrossRef]
  4. Waqas, H.; Farooq, U.; Liu, D.; Abid, M.; Imran, M.; Muhammad, T. Heat transfer analysis of hybrid nanofluid flow with thermal radiation through a stretching sheet: A comparative study. Int. Commun. Heat Mass Transf. 2022, 138, 106303. [Google Scholar] [CrossRef]
  5. Gürdal, M.; Pazarlıoğlu, H.K.; Tekir, M.; Arslan, K.; Gedik, E. Numerical investigation on turbulent flow and heat transfer characteristics of ferro-nanofluid flowing in dimpled tube under magnetic field effect. Appl. Therm. Eng. 2021, 200, 117655. [Google Scholar] [CrossRef]
  6. Ashrafi, H.; Pourmahmoud, N.; Mirzaee, I.; Ahmadi, N. Performance improvement of proton-exchange membrane fuel cells through different gas injection channel geometries. Int. J. Energy Res. 2022, 46, 8781–8792. [Google Scholar]
  7. Hussain, I.; Prakash, D.; Abdalla, B.; Muthtamilselvan, M. Analysis of Arrhenius activation energy and chemical reaction in nanofluid flow and heat transfer over a thin moving needle. Curr. Nanosci. 2023, 19, 39–48. [Google Scholar]
  8. Sadiq, M.N.; Sajid, M.; Abbas, T.; Mahmood, K.; Bashir, S.; Mahmood, W. Numerical Simulation for flow and heat transfer of a nanofluid over lubricated stretchable surface. J. Nanofluids 2023, 12, 506–513. [Google Scholar] [CrossRef]
  9. Prakasha, D.; Sudharani, M.; Kumar, K.G.; Chamkha, A.J. Comparative study of hybrid (graphene/magnesium oxide) and ternary hybrid (graphene/zirconium oxide/magnesium oxide) nanomaterials over a moving plate. Int. Commun. Heat Mass Transf. 2023, 140, 106557. [Google Scholar] [CrossRef]
  10. Tiwari, R.K.; Das, M.K. Heat transfer augmentation in a two-sided liddriven differentially heated square cavity utilizing Nano fluids. Int. J. Heat Mass Transf. 2007, 52, 2002–2018. [Google Scholar]
  11. Khilili, S.; Dinarvand, S.; Hosseini, R.; Saber, M.; Pop, I. Magnetohydrodynamic stagnation-point flow towards stretching/shrinking permeable plate in porous medium filled with a nano fluid. J. Process Mech. Eng. 2014, 228, 309–319. [Google Scholar]
  12. Bakar, F.N.A.; Soid, S.K. MHD Stagnation-Point Flow and Heat Transfer Over an Exponentially Stretching/Shrinking Vertical Sheet in a Micropolar Fluid with a Buoyancy Effect. J. Adv. Res. Numer. Heat Transf. 2022, 8, 50–55. [Google Scholar]
  13. Ali, A.; Marwat, D.N.K.; Ali, A. Analysis of flow and heat transfer over stretching/shrinking and porous surfaces. J. Plast. Film. Sheeting 2022, 38, 21–45. [Google Scholar] [CrossRef]
  14. Azam, M.; Abbas, N.; Kumar, K.G.; Wali, S. Transient bioconvection and activation energy impacts on Casson nanofluid with gyrotactic microorganisms and nonlinear radiation. Waves Random Complex Media 2022, 1–20. [Google Scholar] [CrossRef]
  15. Reddy, M.G.; Krishnamurthy, M.R.; Praveena, M.M.; Naik, L.S.; Prakasha, D.G.; Kumar, K.G. Unsteady absorption flow and dissipation heat transfer over a non-Newtonian fluid. Waves Random Complex Media 2022, 1–15. [Google Scholar] [CrossRef]
  16. Reddy, M.G.; Sudharani, M.V.V.N.L.; Praveena, M.M.; Kumar, K.G. Effect of thermal conductivity on Blasius–Rayleigh–Stokes flow and heat transfer over a moving plate by considering magnetic dipole moment. Eur. Phys. J. Plus 2022, 137, 1–13. [Google Scholar]
  17. Gil, P. Flow and heat transfer characteristics of single and multiple synthetic jets impingement cooling. Int. J. Heat Mass Transf. 2023, 201, 123590. [Google Scholar] [CrossRef]
  18. Şahin, Y.S.; Toprak, B.I.; Solmaz, I.; Bayer, O. Investigation of flow and heat transfer behavior of integrated pin fin-aluminum foam heat sink. Appl. Therm. Eng. 2023, 219, 119504. [Google Scholar] [CrossRef]
  19. Bejawada, S.G.; Nandeppanavar, M.M. Effect of thermal radiation on magnetohydrodynamics heat transfer micropolar fluid flow over a vertical moving porous plate. Exp. Comput. Multiph. Flow 2023, 5, 149–158. [Google Scholar] [CrossRef]
  20. Gaur, P.K.; Jha, A.K. Heat Transfer for MHD Flow in an Inclined Channel with Heat Generation/Absorption. In Advances in Mathematical Modelling, Applied Analysis and Computation; Springer: Singapore, 2023; pp. 273–280. [Google Scholar]
  21. Aljaloud, A.S.M. Hybrid nanofluid mixed convection in a cavity under the impact of the magnetic field by lattice Boltzmann method: Effects of barrier temperature on heat transfer and entropy. Eng. Anal. Bound. Elem. 2023, 147, 276–291. [Google Scholar] [CrossRef]
  22. Yang, D.; Chen, L.; Feng, Y.; Chen, H. Comparisons of Supercritical Loop Flow and Heat Transfer Behavior Under Uniform and Nonuniform High-Flux Heat Inputs. Nucl. Sci. Eng. 2023, 197, 74–91. [Google Scholar] [CrossRef]
  23. Fourier, J. Théorie Analytique de Chaleur. chez Firmin Didot. Pere Et Fils; Editions Jacques Gabay: Paris, France, 1822. [Google Scholar]
  24. Cattaneo, C. Sulla conduzione del calore. Atti Sem. Mat. Fis. Univ. Modena 1948, 3, 83–101. [Google Scholar]
  25. Christov, C.I. On frame indifferent formulation of the Maxwell-Cattaneo model of finite-speed heat conduction. Mech. Res. Commun. 2009, 36, 481–486. [Google Scholar] [CrossRef]
  26. Reddy, M.; Vijayakumari, P.; Kumar, K.; Shehzad, S. Zero-mass flux and Cattaneo–Christov heat flux through a Prandtl non-Newtonian nanofluid in Darcy–Forchheimer porous space. Heat Transf. 2021, 50, 220–233. [Google Scholar] [CrossRef]
  27. Machireddy, G.R.; Praveena, M.; Rudraswamy, N.G.; Kumar, G.K. Impact of Cattaneo–Christov heat flux on hydromagnetic flow of non-Newtonian fluids filled with Darcy–Forchheimer porous medium. Waves Random Complex Media 2021, 1–18. [Google Scholar] [CrossRef]
  28. Khan, S.A.; Hayat, T.; Alsaedi, A. Cattaneo Christov (CC) heat and mass fluxes in Stagnation point flow of Jeffrey nanoliquids by a stretched surface. Chin. J. Phys. 2022, 76, 205–216. [Google Scholar] [CrossRef]
  29. Tausif, S.M.; Das, K.; Kundu, P.K. Modified Homogeneous and Heterogeneous Chemical Reaction and Flow Performance of Maxwell Nanofluid with Cattaneo–Christov Heat Flux Law. J. Eng. Thermophys. 2022, 31, 64–77. [Google Scholar] [CrossRef]
  30. Hayat, T.; Fatima, A.; Muhammad, K.; Alsaedi, A. Heat transfer and entropy analysis in squeezing flow of hybrid nanofluid (Au-CuO/NaAlg) with DF (Darcy-Forchheimer) and CC (Cattaneo-Christov) heat flux. Mater. Sci. Eng. B 2023, 288, 116150. [Google Scholar] [CrossRef]
  31. Ullah, H.; Ullah, K.; Raja, M.A.Z.; Shoaib, M.; Nisar, K.S.; Islam, S.; Weera, W.; Al-Harbi, N. Numerical treatment of squeezed MHD Jeffrey fluid flow with Cattaneo Chrisstov heat flux in a rotating frame using Levnberg-Marquard method. Alex. Eng. J. 2023, 66, 1031–1050. [Google Scholar] [CrossRef]
  32. Shahzad, A.; Imran, M.; Tahir, M.; Khan, S.; Akgül; Abdullaev, S.; Yahia, I.S. Brownian motion and thermophoretic diffusion impact on Darcy-Forchheimer flow of bioconvective micropolar nanofluid between double disks with Cattaneo-Christov heat flux. Alex. Eng. J. 2023, 62, 1–15. [Google Scholar] [CrossRef]
  33. Muhammad, K.; Hayat, T.; Alsaedi, A.; Momani, S. Mixed convective slip flow of hybrid nanofluid (MWCNTs + Cu + Water), nanofluid (MWCNTs + Water) and base fluid (Water): A comparative investigation. J. Therm. Anal. Calorim. 2021, 143, 1523–1536. [Google Scholar] [CrossRef]
  34. Ghadikolaei, S.S.; Hosseinzadeh, K.; Hatami, M.; Ganji, D.D. MHD boundary layer analysis for micropolar dusty fluid containing Hybrid nanoparticles (Cu’Al2O3) over a porous medium. J. Mol. Liq. 2018, 268, 813–823. [Google Scholar] [CrossRef]
  35. Hosseinzadeh, K.; Alizadeh, M.; Ganji, D.D. Hydrothermal analysis on MHD squeezing nanofluid flow in parallel plates by analytical method. Int. J. Mech. Mater. Eng. 2018, 13, 4. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The f ( η ) and F ( η ) illustration against Q .
Figure 1. The f ( η ) and F ( η ) illustration against Q .
Mathematics 11 00554 g001
Figure 2. The f ( η ) and F ( η ) illustration against ϕ 1 .
Figure 2. The f ( η ) and F ( η ) illustration against ϕ 1 .
Mathematics 11 00554 g002
Figure 3. The f ( η ) and F ( η ) illustration against ϕ 2 .
Figure 3. The f ( η ) and F ( η ) illustration against ϕ 2 .
Mathematics 11 00554 g003
Figure 4. The f ( η ) and F ( η ) illustration against ϕ 3 .
Figure 4. The f ( η ) and F ( η ) illustration against ϕ 3 .
Mathematics 11 00554 g004
Figure 5. The f ( η ) and F ( η ) illustration against δ .
Figure 5. The f ( η ) and F ( η ) illustration against δ .
Mathematics 11 00554 g005
Figure 6. The θ ( η ) and θ p ( η ) illustration against R .
Figure 6. The θ ( η ) and θ p ( η ) illustration against R .
Mathematics 11 00554 g006
Figure 7. The θ ( η ) and θ p ( η ) illustration against θ w .
Figure 7. The θ ( η ) and θ p ( η ) illustration against θ w .
Mathematics 11 00554 g007
Figure 8. The θ ( η ) and θ p ( η ) illustration against E c .
Figure 8. The θ ( η ) and θ p ( η ) illustration against E c .
Mathematics 11 00554 g008
Figure 9. The θ ( η ) and θ p ( η ) illustration against Γ .
Figure 9. The θ ( η ) and θ p ( η ) illustration against Γ .
Mathematics 11 00554 g009
Figure 10. The θ ( η ) and θ p ( η ) illustration against ϕ 1 .
Figure 10. The θ ( η ) and θ p ( η ) illustration against ϕ 1 .
Mathematics 11 00554 g010
Figure 11. The θ ( η ) and θ p ( η ) illustration against ϕ 2 .
Figure 11. The θ ( η ) and θ p ( η ) illustration against ϕ 2 .
Mathematics 11 00554 g011
Figure 12. The θ ( η ) and θ p ( η ) illustration against ϕ 3 .
Figure 12. The θ ( η ) and θ p ( η ) illustration against ϕ 3 .
Mathematics 11 00554 g012
Figure 13. The θ ( η ) and θ p ( η ) illustration against B i .
Figure 13. The θ ( η ) and θ p ( η ) illustration against B i .
Mathematics 11 00554 g013
Table 1. Thermophysical properties of a ternary nanofluid (graphene + sirconium oxide + tungsten oxide) (Prakasha et al. [9]) and Muhammad et al. [33]).
Table 1. Thermophysical properties of a ternary nanofluid (graphene + sirconium oxide + tungsten oxide) (Prakasha et al. [9]) and Muhammad et al. [33]).
Thermophysical PropertiesBase FluidTernary Nanofluid
H 2 O GrapheneZirconium OxideTungsten OXIDE
ρ ( kg / m 2 )997.122005680 7160
C p   ( j / kgK ) 41795000502 96.15
k ( w / mk ) 0.6137901.7 1.63
Table 2. Validation of the results in comparison with those found in the literature.
Table 2. Validation of the results in comparison with those found in the literature.
Pr Ghadikolaei et al. [34]Hosseinzadeh et al. [35]Reddy et al. [3]Present Results
0.70.45380.45410.45390.45415
2.00.91130.91140.91130.91133
7.01.89541.89541.89541.89545
Table 3. The friction factor values for various pertinent parameters.
Table 3. The friction factor values for various pertinent parameters.
Q δ ϕ 1 ϕ 2 ϕ 3 Friction Factor
Dusty PhaseTernary Phase
0.5 0.7852391.511657
1.0 0.6130421.447579
1.5 0.5080941.406297
0.5 0.6202071.470950
1.0 0.5411871.464362
1.5 0.4460561.453310
0.005 0.7405011.477960
0.001 0.6882521.475486
0.015 0.6408701.472762
0.005 0.9674951.480191
0.001 0.8283021.475561
0.015 0.5494981.460959
0.0050.7405011.477961
0.0010.6882521.475486
0.0150.5494981.460959
Table 4. The Nusselt number values for various pertinent parameters.
Table 4. The Nusselt number values for various pertinent parameters.
Pr E c B i Γ R θ w ϕ 1 ϕ 2 ϕ 3 Nusselt Number
Dust PhaseTernary Phase
0.5 1.446971.56233
1.0 1.461932.69770
1.5 1.482743.51973
0.1 1.502431.94796
0.3 1.437311.93508
0.5 1.394141.92532
0.1 1.511662.23259
0.2 1.447581.74300
0.3 1.406341.44461
0.1 1.571421.97687
0.2 1.849962.01510
0.3 2.083912.01747
0.5 1.458661.91198
1.0 1.536131.88219
1.5 1.579341.86686
0.8 1.286821.48919
1.2 1.328621.64629
1.4 1.368651.79380
0.005 1.062391.88550
0.001 1.164551.89615
0.015 1.278721.91104
0.005 1.394141.92532
0.001 1.437311.93508
0.015 1.511662.23259
0.0051.174551.87615
0.0011.228721.89104
0.0151.304141.92532
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Souayeh, B. Simultaneous Features of CC Heat Flux on Dusty Ternary Nanofluid (Graphene + Tungsten Oxide + Zirconium Oxide) through a Magnetic Field with Slippery Condition. Mathematics 2023, 11, 554. https://doi.org/10.3390/math11030554

AMA Style

Souayeh B. Simultaneous Features of CC Heat Flux on Dusty Ternary Nanofluid (Graphene + Tungsten Oxide + Zirconium Oxide) through a Magnetic Field with Slippery Condition. Mathematics. 2023; 11(3):554. https://doi.org/10.3390/math11030554

Chicago/Turabian Style

Souayeh, Basma. 2023. "Simultaneous Features of CC Heat Flux on Dusty Ternary Nanofluid (Graphene + Tungsten Oxide + Zirconium Oxide) through a Magnetic Field with Slippery Condition" Mathematics 11, no. 3: 554. https://doi.org/10.3390/math11030554

APA Style

Souayeh, B. (2023). Simultaneous Features of CC Heat Flux on Dusty Ternary Nanofluid (Graphene + Tungsten Oxide + Zirconium Oxide) through a Magnetic Field with Slippery Condition. Mathematics, 11(3), 554. https://doi.org/10.3390/math11030554

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