Next Article in Journal
Regularity, Asymptotic Solutions and Travelling Waves Analysis in a Porous Medium System to Model the Interaction between Invasive and Invaded Species
Next Article in Special Issue
Simulation of Dissipative Hybrid Nanofluid (PEG-Water + ZrO2 + MgO) Flow by a Curved Shrinking Sheet with Thermal Radiation and Higher Order Chemical Reaction
Previous Article in Journal
Optimization of the Cognitive Processes in a Virtual Classroom: A Multi-objective Integer Linear Programming Approach
Previous Article in Special Issue
Peristaltic Transport of Carreau Nanofluid in Presence of Triple Diffusion in an Asymmetric Channel by Multi-Step Differential Transformation Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

3D Flow of Hybrid Nanomaterial through a Circular Cylinder: Saddle and Nodal Point Aspects

by
Javali K. Madhukesh
1,†,
Gosikere K. Ramesh
2,
Govinakovi S. Roopa
3,
Ballajja C. Prasannakumara
1,
Nehad Ali Shah
4,† and
Se-Jin Yook
5,*
1
Department of Mathematics, Davangere University, Davangere 577007, India
2
Department of Mathematics, K.L.E. Society’s J.T. College, Gadag 582101, India
3
Department of Mathematics, Malnad College of Engineering, Hassan 573202, India
4
Department of Mechanical Engineering, Sejong University, Seoul 05006, Korea
5
School of Mechanical Engineering, Hanyang University, 222 Wangsimni-ro, Seongdong-gu, Seoul 04763, Korea
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work and are co-first authors.
Mathematics 2022, 10(7), 1185; https://doi.org/10.3390/math10071185
Submission received: 9 March 2022 / Revised: 1 April 2022 / Accepted: 2 April 2022 / Published: 5 April 2022
(This article belongs to the Special Issue Mathematical Problems in Mechanical Engineering)

Abstract

:
This mathematical model explains the behavior of sinusoidal radius activity in stagnation point three-dimensional flow of hybrid nanoparticles through a circular cylinder. The energy equation of heat source/sink effect and the mass equation of Arrhenius energy of activation and chemical reaction effects are incorporated. Self-relation transformations are adopted to reduce the PDEs to ODEs, then the RKF-45 method is solved with shooting proficiency. The nodal and saddle point action is studied in pertinent parameters for thermal, mass, and velocity curves. Further statistical values of skin friction, Nusselt number, and Sherwood number of both nodal and saddle points are portrayed in tables format. It is ascertained that higher values of activation energy and reaction rate enhance the concentration curve. In addition, the nodal point curves are always less than saddle point curves.

1. Introduction

Nanofluids are very fine engineering colloidal particles created by combining nanoparticles and a base liquid. The diameter of nanoparticles typically ranges from 1–100 nm, and the base liquid is typically water or ethyl glycerol. Choi [1] in 1995 was the first to suggest the effects of thermal conductivity and heat transmission across nanofluids. In addition, based on his assessment, he discovered that nanofluids are more efficient than base fluids. Because the nanofluid is one of the developing liquids, it has resulted in numerous scientific inventions in the nano industries, according to Choi’s survey [2]. Many of academics have examined a range of analyses of fluid models across various platforms (see refs. [3,4,5,6,7,8,9]). Hybrid nanofluids have recently taken precedence over nanofluids due to increased thermal conductivity and a more efficient heat transference rate. Hybrid nanofluids are created by combining two or more nanoparticles with a base liquid. The various applications of hybrid nanofluids were studied by Sarkar et al. [10]. Stability and various characteristics of CuO–TiO2/water-based hybrid nanofluid were examined by Asadi et al. [11]. Lund et al. [12] studied the dual solution of the impact of thermal radiation effect over magnetohydrodynamic flow of a Cu–Al2O3/H2O-based hybrid nanofluid. Ramesh et al. [13] investigated the Darcy–Forchheimer model of hybrid nanofluid in a stretchable convergent/divergent channel. They found variations of velocity and solid volume fraction in the channel. Ramesh and Madhukesh [14] investigated the effect of hybrid carbon nanotubes in the presence of activation energy and heat source/sink. They found that hybrid CNTs nanomaterials have a greater rate of heating/cooling than single CNTs nanomaterials.
Many aspects, such as activation energy, heat source/sink, and slip parameters, all play a significant part in thermal and concentration analyses of fluid flow. These variables enable us to examine the behavior of fluids over a range of parameters. Alghamdi et al. [15] studied the role of dyadic chemical response for the MHD flow of viscous fluid with suspended nanoparticles through the gyrating disk. Rekha et al. [16] considered AA7072-AA7075/water-based hybrid nanofluid flowing through a cone, wedge, and plate in the presence of activation energy. They discussed the impact of activation energy and mass transfer phenomena in the presence of hybrid nanofluid. Ramesh [17] used Buongiorno’s model to analyze the significance of chemical reaction and energy of activation over a viscoelastic nanomaterial through a surface with stretching. Further, this model was taken and work was extended by Alsadi et al. [18], who, over stretching sheet, studied the entropy generation impact, and Asma et al. [19], who examined the influence of Arrhenius energy of activation on Darcy–Forchheimer nano liquid in three dimension.
A three-dimensional Williamson nanofluid with Darcy–Forchheimer permeable media was emphatically deliberated upon by Ramzan et al. [20]. They found the behavior of heat flux under the Fourier and Fick’s laws and examined the Williamson and Hartman number over a velocity profile. Jagan et al. [21] deliberated the 3D unsteady MHD convective stream of nano liquid with double stratification using the homotopy analysis method. Nayak et al. [22] explored the three-dimensional free convective moment of a radioactive nanofluid under a variable magnetic field. They observed the influence of the Hartmann number on the temperature and thickness of the related boundary layer. Irfan et al. [23] scrutinized the nonlinear mixed convection on radioactive flow and activation energy on the 3D flow of a Carreau nanofluid. Alwawi et al. [24] investigated the ethyl-glycerol-based nanofluid flow around a circular cylinder by utilizing the Tiwari and Das [25] nanofluid model. In this study, they considered the magnetohydrodynamic effect and analyzed the heat transfer effects.
Inspired by all the above insights, we developed a 3D flow of hybrid nanoparticles ( Fe 3 O 4 + Go ) through a circular cylinder in the stagnation point region. From the literature point of view, the problem has not yet been studied by considering activation energy. Normal partial coupled equations are regenerated ordinary coupled equations. Furthermore, equations are solved numerically with the prominent RKF-45 method. A good discussion is made with the help of graphs.

2. Problem Formulation

Consider a time-independent, three-dimensional flow of a hybrid nanofluid ( Fe 3 O 4 + Go ) with water as a base fluid in the direction of the circular cylinder (see Figure 1). The radius of the cylinder variegates sinusoidally. The positions B, A, and C on the cylinder represent the minimum and maximum stagnation points. A connecting line connects A to B and B to C, splitting the flow that passes on either side of the cylinder. Additionally, the flow is taken with u e * & v e * free stream speed along the x & y direction, with c = b * a * representing the special streamline. The flow lines of the nodal stagnation point are ranging from 0 < c < 1 , The flow lines of the saddle stagnation point are ranging from 1 < c < 0 . The flow is said to plane when c = 0 . The Arrhenius energy of activation, chemical reaction, and uniform heat source/sink effects are examined. The equations which represent the above-assumed flow are given as follows (see [26,27]).
Continuity equation
u x + v y + w z = 0
Momentum equation
u u x + v u y + w u z = a * 2 x + ν h n f 2 u z 2
u v x + v v y + w v z = b * 2 y + ν h n f 2 v z 2
Temperature equation
u T x + v T y + w T z = α h n f 2 T z 2 + Q 0 ( ρ c p ) h n f ( T T )
Concentration equation
u C x + v C y + w C z = D h n f 2 C z 2 κ r 2 ( T T ) n e E a K T ( C C )
The interrelated boundary constraints are
u = 0 , v = 0 , w = 0 , γ k h n f T z = T T w , C = C w a t z = 0 ,
u u e * , v v e * , T T , C C   a s z
where ( u ,   v ,   w   and   T ) are the components of velocity and temperature along x, y, and z directions. a * 2 and b * 2 represent the free stream-dependent constants. The term ν ( = μ ρ ) denotes kinematic viscosity, and μ and ρ represent dynamic viscosity and density of the hybrid nanofluid, respectively. α ( = k ( ρ C p ) ) denotes hybrid nanofluid’s thermal diffusivity, and the term k represents thermal conductivity. T w & C w represent the temperature and concentration at the wall, and T & C represent ambient temperature and concentration. Q 0 is the uniform heat source/sink coefficient, κ r 2 is the reaction rate, n denotes the fitted rate constant, K is the Boltzmann constant, E a is activation energy, and ( T T ) n e E a K T is modified Arrhenius function. γ is thermal slip, and the subscript h n f denotes hybrid nanofluid. Thermophysical features of water and nanoparticles are listed in Table 1 (see [28]).
The efficient characteristics of hybrid nanofluid are stated below.
D h n f = ( 1 ϕ 1 ) 2.5 ( 1 ϕ 2 ) 2.5 D f , ( ρ C p ) h n f = ( 1 ϕ 1 ) 2.5 ( 1 ϕ 2 ) 2.5 ( ρ C p ) f + ϕ 1 ( ρ C p ) S 1 + ϕ 2 ( ρ C p ) S 2 , μ h n f = μ f ( 1 ϕ 1 ) 2.5 ( 1 ϕ 2 ) 2.5 , ρ h n f = ( 1 ϕ 2 ) [ ( 1 ϕ 1 ) ρ f + ϕ 1 ρ s 1 ] + ϕ 2 ρ s 2 , ( ρ c p ) h n f = ( 1 ϕ 2 ) [ ( 1 ϕ 1 ) ( ρ c p ) f + ϕ 1 ( ρ c p ) s 1 ] + ϕ 2 ( ρ c p ) s 2 , k h n f = k s 2 + ( m 1 1 ) k b f ( m 1 1 ) ϕ 2 ( k b f k s 2 ) k s 2 + ( m 1 1 ) k b f + ϕ 2 ( k b f k s 2 ) k b f , where   k b f = k s 1 + ( m 1 1 ) k f ( m 1 1 ) ϕ 1 ( k f k s 1 ) k s 1 + ( m 1 1 ) k f + ϕ 1 ( k f k s 1 ) k f }
Here, ϕ 1 , ϕ 2 indicates the volume fraction of individual nanoparticles. Here, suffixes h n f denote hybrid nanofluid, f denote base fluid, and s denote solid particle, and m 1 is shape factor parameter. In the current study, we take a sphere shape ( m 1 = 3 ) [29].
The variables are used to obtain the similarity equations:
u = a * x f ( η ) , v = b * y g ( η ) , w = a * ν f ( f ( η ) + c g ( η ) ) , T = T + ( T w T ) θ ( η ) , C = C + ( C w C ) χ ( η ) , η = z a * ν f . }
The transformed ordinary differential equations along with the conditions are
1 ( 1 ϕ 1 ) 2.5 ( 1 ϕ 2 ) 2.5 [ { ( 1 ϕ 2 ) ( 1 ϕ 1 ) + ϕ 1 ρ s 1 ρ f } + ϕ 2 ρ s 2 ρ f ] f + f f + c f g f 2 + 1 = 0
1 ( 1 ϕ 1 ) 2.5 ( 1 ϕ 2 ) 2.5 [ { ( 1 ϕ 2 ) ( 1 ϕ 1 ) + ϕ 1 ρ s 1 ρ f } + ϕ 2 ρ s 2 ρ f ] g + f g + c g g c g 2 + c = 0
k h n f k f Pr [ { ( 1 ϕ 2 ) ( 1 ϕ 1 ) + ϕ 1 ( ρ c p ) s 1 ( ρ c p ) f } + ϕ 2 ( ρ c p ) s 2 ( ρ c p ) f ] θ + f θ + c θ g + S θ = 0
( 1 ϕ 1 ) 2.5 ( 1 ϕ 2 ) 2.5 χ + S c f χ + c   S c χ g β   S c ( 1 + δ θ ) n E x p ( E 1 + δ θ ) χ = 0
The reduced boundary constraints are
f ( 0 ) = 0 , g ( 0 ) = 0 , f ( 0 ) = 0 , g ( 0 ) = 0 , θ ( 0 ) = 1 + λ k h n f k f θ ( 0 ) ,   χ ( 0 ) = 1 }   at   η 0
f ( ) 1 , g ( ) 1 , θ ( ) 0 ,   χ ( ) 0   as   η
The variables which are present in the Equations (9)–(13) are as follows:
c = b * a * is a constant that represents a streamline, λ = γ k f a * ν f is the thermal slip parameter, Pr = μ f ( C p ) f k f is a Prandtl number, S c = ν f D f   is a Schmidt number, E = E a K T is the activation energy, δ = T w T T is temperature difference parameter, β = κ r 2 a * is reaction rate, and S = Q 0 a * ( ρ C p ) h n f is the heat source/sink parameter.
Skin friction, Nusselt, and Sherwood numbers are given by
C f x = τ w x ρ f u w 2 , C f y = τ w y ρ f u w 2 , N u x = x q w k f ( T w T ) , S h x = D f x j w ( C w C )
where τ w x and τ w y are the shear stresses surface in the x and y directions.
The surface heat flux is denoted as q w and mass transfer j w is defined as
τ w x = [ μ h n f u z ] z = 0 ,   τ w y = [ μ h n f v z ] z = 0 ,   q w = k h n f ( T z ) z = 0 ,   j w = 1 D h n f ( C z ) z = 0
From the above equations, the skin friction, Sherwood, and Nusselt numbers are
Re x C f x = f ( 0 ) ( 1 ϕ 1 ) 2.5 ( 1 ϕ 2 ) 2.5 , ( x / y ) Re x C f y = c g ( 0 ) ( 1 ϕ 1 ) 2.5 ( 1 ϕ 2 ) 2.5 ,   N u x Re x = k h n f k f θ ( 0 ) ,   S h x Re x = χ ( 0 ) ( 1 ϕ 1 ) 2.5 ( 1 ϕ 2 ) 2.5   .
where Re = x 2 a * ν f is Reynolds number.

3. Numerical Procedure and Validation of Code

The reduced equations as stated in (10)–(13) and boundary conditions (14) and (15) are numerically tackled with the help of a well-known numerical approach i.e., Runge–Kutta–Fehlberg-45 order. To obtain the numerical solution, the following steps are followed.
  • Convert the BVP into IVP of the first order.
  • Apply the shooting procedure to guess the missing boundary values.
  • Apply the RKF-45 method to obtain the solution to IVP.
  • Find the residuals for all the boundary conditions.
  • If the residual error is greater than the error tolerance, adjust the initial guesses.
  • If the residual error is less than error tolerance, numerical results are obtained.
Let us consider
Δ 1 = ( 1 ϕ 1 ) 2.5 ( 1 ϕ 2 ) 2.5 , Δ 2 = { ( 1 ϕ 2 ) ( 1 ϕ 1 ) + ϕ 1 ρ s 1 ρ f } + ϕ 2 ρ s 2 ρ f , Δ 3 = { ( 1 ϕ 2 ) ( 1 ϕ 1 ) + ϕ 1 ( ρ c p ) s 1 ( ρ c p ) f } + ϕ 2 ( ρ c p ) s 2 ( ρ c p ) f , f = ϑ 1 , f = ϑ 2 , f = ϑ 3 g = ϑ 4 , g = ϑ 5 , g = ϑ 6 θ = ϑ 7 , θ = ϑ 8 χ = ϑ 9 , χ = ϑ 10
ϑ 3 = Δ 1 Δ 2 ( ϑ 1 ϑ 3 + c ϑ 3 ϑ 4 ϑ 2 2 + 1 )
ϑ 6 = Δ 1 Δ 2 ( ϑ 1 ϑ 6 + c ϑ 6 ϑ 4 c ϑ 5 2 + c )
ϑ 8 = Pr Δ 3 k f ( ϑ 1 ϑ 8 + c ϑ 8 ϑ 4 + S ϑ 7 ) / k h n f
ϑ 10 = ( S c ϑ 1 ϑ 10 + c   S c ϑ 10 ϑ 4 β   S c ( 1 + δ ϑ 7 ) n E x p ( E 1 + δ ϑ 7 ) ϑ 9 ) / Δ 1
and
ϑ 2 ( 0 ) = 0 , ϑ 5 ( 0 ) = 0 , ϑ 1 ( 0 ) = 0 , ϑ 4 ( 0 ) = 0 , ϑ 7 ( 0 ) = 1 + λ k h n f k f ϑ 8 ( 0 ) ,   ϑ 9 ( 0 ) = 1 &   ϑ 3 ( 0 ) = ξ 1 , ϑ 6 ( 0 ) = ξ 2 , ϑ 10 ( 0 ) = ξ 3 }
The algorithm for the numerical procedure (Figure 2) and RKF-45 scheme ([30,31,32]) is given as follows.
In the RKF-45 method, each step uses the following values:
s 1 = f ( x i , y i ) h 1 s 2 = f ( x i + 1 4 h , y i + 1 4 s 1 ) h 1 s 3 = f ( x i + 3 8 h , y i + 3 32 s 1 + 9 32 s 2 ) h 1 s 4 = f ( x i + 12 13 h , y i + 1932 2147 s 1 7200 2147 s 2 + 7296 2147 s 3 ) h 1 s 5 = f ( x i + h , y i + 439 216 s 1 845 4104 s 4 8 s 2 + 3680 513 s 3 ) h 1 s 6 = f ( x i + 1 2 h , y i 8 27 s 1 + 2 s 2 11 40 s 5 3544 2565 s 3 1859 4104 s 4 ) h 1
Then, using a Runge–Kutta technique of order 4, an estimate of the I.V.P. result is generated.
y i + 1 = y i + 25 216 s 1 1 5 s 5 + 2197 4104 s 4 + 1408 2565 s 3
The Runge–Kutta technique of order 5 is given by
z i + 1 = y i + 16 135 s 1 + 2 55 s 6 9 50 s 5 + 28561 56430 s 4 + 6656 12825 s 3
The numerical results are validated with available previously published works (see Gangadar et al. [33], Bhattacharyya and Gupta [34], Dinarvand [35], Bachok et al. [36]) and those are found that best match each other in Table 2.

4. Results and Discussion

The main aspect of the current section is to frame up the research, the influence of solid volume fraction, heat source/sink parameter, activation energy, temperature difference parameter, reaction rate, and Schmidt number. To understand this, the ODEs (9)–(12) with the boundary constraints (13) are solved by applying the Runge–Kutta–Felberg-45 order technique along with the shooting method. The thermophysical characteristics of water and nanoparticles are represented in Table 1. Numerical validation of the present study with existing works is provided in the Table 2. Computational values of f ( 0 ) , g ( 0 ) , θ ( 0 ) and χ ( 0 ) are represented in Table 3 for diverse values of c . The computational values of f ( 0 ) , g ( 0 ) , θ ( 0 ) and χ ( 0 ) with Pr = 6.2 , n = 0.2 , α = 0.5 and ϕ 1 = 0.1 for saddle point are presented in Table 4 and for nodal point are presented in Table 5.
From Table 4, it is detected that rise in value of ϕ 2 upsurges the profiles f ( 0 ) , g ( 0 ) , θ ( 0 ) and χ ( 0 ) . Nodal point region is an inclined function of ϕ 2 . It is noted that skin friction, Nusselt, and Sherwood numbers are also increased by increase in the solid volume fraction but the in saddle point region, skin friction for f ( 0 ) profile increases. Whereas Nusselt and Sherwood numbers increase with a decrease in g ( 0 ) profile, for all increasing values, the heat source/sink parameter S , Nusselt, and Sherwood numbers decrease and there is no variation observed in f ( 0 ) and g ( 0 ) profiles. The same observations are seen in the saddle point region also. For the temperature difference parameter δ , the Nusselt number increases but the Sherwood number diminishes gradually. For saddle point, Nusselt number decreases but variations are observed in Sherwood number. For Schmidt number S c , Sherwood number increases in both nodal and saddle point regions. For reaction rate parameter β , Sherwood number increases but decreases in activation energy E parameter in both nodal and saddle point, flow can be seen in Table 5.
Figure 3, Figure 4, Figure 5 and Figure 6 illustrate the impact of special streamline parameter c over f ( η ) ,   g ( η ) ,   θ ( η )   and   χ ( η ) profiles. From Figure 3 and Figure 4 it is detected that at saddle point, a rise in the value of c slightly increases the hybrid nanofluid velocity in both the profiles f ( η )   &   g ( η )   , and the same behavior is observed for nodal point also, but temperature and concentration profiles show the opposite behavior, which is depicted in Figure 5 and Figure 6.
Figure 7 displays the impact of solid volume fraction ϕ 2 over f ( η ) profile. It is evident from the figure that a rise in the value of the solid fraction ϕ 2 increases the velocity f ( η ) in both the nodal and saddle points. From the physical point of interpretation, the addition of Go nanoparticles increases the volume fraction, and velocity increases. The same behavior is seen in the g ( η ) profile (see Figure 8). Moreover, it is noted that the velocity of the nodal point is comparatively more than the saddle point. Figure 9 displays the impact of solid volume fraction ϕ 2 over temperature profile θ ( η ) ; as the volume of the solid fraction increases, the thermal distribution continues increasing from the nodal stagnation point to the saddle stagnation point. Figure 10 draws the impact of solid volume fraction ϕ 2 over concentration profile χ ( η ) . As the volume of the solid fraction increases, the thickness of boundary layer also upsurges, leading to a rise in the concentration of the hybrid nanofluid from the nodal stagnation point to the saddle stagnation point.
Figure 11 validates the impact of heat source/sink parameter ( S ) over the thermal gradient θ ( η ) . S = (<0,0,>0) denotes heat sink, absence of heat source/sink, and heat source, respectively. The existence of a heat sink acts as an exchanger, transmitting heat from the surface into the hybrid nanofluid. Temperature is produced by the system in the presence of a heat source. One can see from Figure 11 that as the heat source/sink parameter values continue increasing, the thermal distribution also increases, because of an upsurge in the thickness of the boundary layer. In the case of the heat sink, the thermal dispersion is low, and the temperature is generated by the interface in the case of the heat source. In this case, a heat source outperforms a heat sink in terms of thermal performance. In addition, it is evident from the diagram that thermal distribution is more in the saddle point than the nodal point. Figure 12 displays the effect of temperature difference parameter ( δ ) on thermal profile θ ( η ) . Thermal distribution continues diminishing as the temperature difference parameter values rise. The thermal difference in the saddle point is better than a nodal point of flow.
Figure 13 provides the influence of S c over χ ( η ) . The concentration curve continues decreasing with the slight upsurge in the values of Sc. A similar behavior is observed in both the saddle and nodal points of the flow. The same observation is seen for the reaction rate parameter β , which is plotted in Figure 14. Figure 15 gives the impact of the activation energy parameter E over the concentration profile χ ( η ) . One can observe from the figure that the concentration curve continues increasing with the rise in value of E . The saddle point flow exhibits more energy distribution than the nodal point flow.

5. Final Remarks

The key intention of this article is to draw out the behavior of sinusoidal radius activity in stagnation point three-dimensional flow through a circular cylinder with the suspension of hybrid nanoparticles. The best results are obtained by applying the RKF-45 method and shooting method. Considerable outcomes of this study are given below.
  • A rise in the value of the streamline parameter upsurges the flow velocity and decreases the thermal distribution and concentration.
  • Better thermal gradient and concentration are seen in the enhancement of volume fraction.
  • Boundary layer thickness and concentration are decreased with increase in the Schmidt number.
  • Thermal distribution and concentration are more in the saddle point than in the nodal point of a hybrid nanofluid.

Author Contributions

Conceptualization, J.K.M. and N.A.S.; methodology, G.K.R.; software, G.S.R. and N.A.S.; validation, B.C.P. and S.-J.Y.; resources, N.A.S. and J.K.M.; data curation, G.S.R.; writing—original draft preparation, G.K.R. and B.C.P.; writing—review and editing, N.A.S. and S.-J.Y.; visualization, G.S.R.; project administration, G.K.R.; funding acquisition, S.-J.Y., J.K.M. and N.A.S. contributed equally to this work and are co-first authors. All authors have read and agreed to the published version of the manuscript.

Funding

This work was conducted under the Technology Innovation Program (or Industrial Strategic Technology Development Program—material part package type) as “Development of fire suppression-type high safety module and demonstration of safety for future eco-friendly medium and large secondary battery (No. 20015986)”, funded by the Ministry of Trade, Industry & Energy (MOTIE), Republic of Korea.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

a * & b * Free stream dependent constants
c Gradient of streamline
C Concentration
C w Wall concentration
C Ambient concentration
C p Specific heat
C f x &   C f y    Skin friction along x and y direction
D Diffusivity
EActivation energy parameter
E a Activation energy
f Fluid
f ( η ) Dimensionless velocity
g ( η ) Dimensionless velocity
h n f Hybrid nanofluid
j w Mass transfer
k Thermal conductivity
κ r 2 Reaction rate
K Boltzmann constant
n Fitted rate constant
N u x Nusselt number
PrPrandtl number
q w Surface heat flux
Q 0 Uniform heat source/sink coefficient
Re Reynolds number
s 1 Solid particle of Fe 3 O 4
s 2 Solid particle of Go
S Heat source/sink parameter
S c Schmidt number
S h x Sherwood number
T Temperature
T w Wall temperature
T Ambient temperature
u e *   & v e * Free stream velocity
u , v & w Velocity components
x , y & z Coordinate axis
Greek symbols
α Thermal diffusivity
μ Dynamic viscosity
ρ Density
ν Kinematic viscosity
ρ C p Heat capacitance
  γ Thermal slip
λ Thermal slip parameter
δ Temperature difference parameter
β Reaction rate
τ w x Shear stresses surface in the x-direction
τ w y Shear stresses surface in the y-direction
ϕ 1 The solid volume fraction of Fe 3 O 4
ϕ 2 The solid volume fraction of Go
θ ( η ) Dimensionless temperature
χ ( η ) Dimensionless concentration

References

  1. Choi, S.U.S.; Eastman, J. Enhancing thermal conductivity of fluids with nanoparticles. In Proceedings of the ASME International Mechanical Engineering Congress and Exposition, San Francisco, CA, USA, 12–17 November 1995; Volume 66. [Google Scholar]
  2. Choi, S.U.S. Nanofluids: From Vision to Reality Through Research. J. Heat Transf. 2009, 131, 033106. [Google Scholar] [CrossRef]
  3. 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] [PubMed]
  4. Rehman, F.U.; Nadeem, S.; Rehman, H.U.; Haq, R.U. Thermophysical analysis for three-dimensional MHD stagnation-point flow of nano-material influenced by an exponential stretching surface. Results Phys. 2018, 8, 316–323. [Google Scholar] [CrossRef]
  5. Fetecau, C.; Shah, N.A.; Vieru, D. General Solutions for Hydromagnetic Free Convection Flow over an Infinite Plate with Newtonian Heating, Mass Diffusion and Chemical Reaction. Commun. Theor. Phys. 2017, 68, 768. [Google Scholar] [CrossRef]
  6. Sheikholeslami, M. Magnetic field influence on nanofluid thermal radiation in a cavity with tilted elliptic inner cylinder. J. Mol. Liq. 2017, 229, 137–147. [Google Scholar] [CrossRef]
  7. Chamkha, A.J.; Abbasbandy, S.; Rashad, A.M.; Vajravelu, K. Radiation Effects on Mixed Convection over a Wedge Embedded in a Porous Medium Filled with a Nanofluid. Transp. Porous Media 2011, 91, 261–279. [Google Scholar] [CrossRef]
  8. Animasaun, I.L.; Shah, N.A.; Wakif, A.; Mahanthesh, B.; Sivaraj, R.; Koriko, O.K. Ratio of Momentum Diffusivity to Thermal Diffusivity: Introduction, Meta-Analysis, and Scrutinization, 1st ed.; Chapman and Hall/CRC: London, UK, 2022. [Google Scholar] [CrossRef]
  9. Saeed, M.; Kim, M.-H. Heat transfer enhancement using nanofluids (Al2O3-H2O) in mini-channel heatsinks. Int. J. Heat Mass Transf. 2018, 120, 671–682. [Google Scholar] [CrossRef]
  10. Sarkar, J.; Ghosh, P.; Adil, A. A review on hybrid nanofluids: Recent research, development and applications. Renew. Sustain. Energy Rev. 2015, 43, 164–177. [Google Scholar] [CrossRef]
  11. Asadi, A.; Alarifi, I.M.; Foong, L.K. An experimental study on characterization, stability and dynamic viscosity of CuO-TiO2/water hybrid nanofluid. J. Mol. Liq. 2020, 307, 112987. [Google Scholar] [CrossRef]
  12. Lund, L.A.; Omar, Z.; Khan, I.; Sherif, E.-S.M. Dual Solutions and Stability Analysis of a Hybrid Nanofluid over a Stretching/Shrinking Sheet Executing MHD Flow. Symmetry 2020, 12, 276. [Google Scholar] [CrossRef] [Green Version]
  13. Ramesh, G.K.; Shehzad, S.A.; Tlili, I. Hybrid nanomaterial flow and heat transport in a stretchable convergent/divergent channel: A Darcy-Forchheimer model. Appl. Math. Mech. 2020, 41, 699–710. [Google Scholar] [CrossRef]
  14. Ramesh, G.; Madhukesh, J. Activation energy process in hybrid CNTs and induced magnetic slip flow with heat source/sink. Chin. J. Phys. 2021, 73, 375–390. [Google Scholar] [CrossRef]
  15. Alghamdi, M. Significance of Arrhenius Activation Energy and Binary Chemical Reaction in Mixed Convection Flow of Nanofluid Due to a Rotating Disk. Coatings 2020, 10, 86. [Google Scholar] [CrossRef] [Green Version]
  16. Rekha, M.B.; Sarris, I.E.; Madhukesh, J.K.; Raghunatha, K.R.; Prasannakumara, B.C. Activation Energy Impact on Flow of AA7072-AA7075/Water-Based Hybrid Nanofluid through a Cone, Wedge and Plate. Micromachines 2022, 13, 302. [Google Scholar] [CrossRef] [PubMed]
  17. Ramesh, G. Analysis of active and passive control of nanoparticles in viscoelastic nanomaterial inspired by activation energy and chemical reaction. Phys. A Stat. Mech. Its Appl. 2020, 550, 123964. [Google Scholar] [CrossRef]
  18. Alsaadi, F.E.; Ullah, I.; Hayat, T.; Alsaadi, F.E. Entropy generation in nonlinear mixed convective flow of nanofluid in porous space influenced by Arrhenius activation energy and thermal radiation. J. Therm. Anal. 2019, 140, 799–809. [Google Scholar] [CrossRef]
  19. Asma, M.; Othman, W.; Muhammad, T. Numerical study for Darcy–Forchheimer flow of nanofluid due to a rotating disk with binary chemical reaction and Arrhenius activation energy. Mathematics 2019, 7, 921. [Google Scholar] [CrossRef] [Green Version]
  20. Ramzan, M.; Gul, H.; Zahri, M. Darcy-Forchheimer 3D Williamson nanofluid flow with generalized Fourier and Fick’s laws in a stratified medium. Bull. Pol. Acad. Sci. Tech. Sci. 2020, 68, 327–335. [Google Scholar]
  21. Jagan, K.; Sivasankaran, S.; Bhuvaneswari, M.; Rajan, S. Effect of Non-linear Radiation on 3D Unsteady MHD Nanoliquid Flow over a Stretching Surface with Double Stratification. Trends Math. 2019, 109–116. [Google Scholar] [CrossRef]
  22. Nayak, M.K.; Shaw, S.; Chamkha, A.J. 3D MHD Free Convective Stretched Flow of a Radiative Nanofluid Inspired by Variable Magnetic Field. Arab. J. Sci. Eng. 2019, 44, 1269–1282. [Google Scholar] [CrossRef]
  23. Irfan, M.; Khan, W.A.; Khan, M.; Gulzar, M.M. Influence of Arrhenius activation energy in chemically reactive radiative flow of 3D Carreau nanofluid with nonlinear mixed convection. J. Phys. Chem. Solids 2019, 125, 141–152. [Google Scholar] [CrossRef]
  24. A Alwawi, F.; Alkasasbeh, H.T.; Rashad, A.; Idris, R. Heat transfer analysis of ethylene glycol-based Casson nanofluid around a horizontal circular cylinder with MHD effect. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2020, 234, 2569–2580. [Google Scholar] [CrossRef]
  25. Tiwari, R.K.; Das, M.K. Heat transfer augmentation in a two-sided lid-driven differentially heated square cavity utilizing nanofluids. Int. J. Heat Mass Transf. 2007, 50, 2002–2018. [Google Scholar] [CrossRef]
  26. Nadeem, S.; Abbas, N. On both MHD and slip effect in micropolar hybrid nanofluid past a circular cylinder under stagnation point region. Can. J. Phys. 2019, 97, 392–399. [Google Scholar] [CrossRef]
  27. Tlili, I.; Waqas, H.; Almaneea, A.; Khan, S.U.; Imran, M. Activation Energy and Second Order Slip in Bioconvection of Oldroyd-B Nanofluid over a Stretching Cylinder: A Proposed Mathematical Model. Processes 2019, 7, 914. [Google Scholar] [CrossRef] [Green Version]
  28. Shah, N.A.; Wakif, A.; El-Zahar, E.R.; Thumma, T.; Yook, S.-J. Heat transfers thermodynamic activity of a second-grade ternary nanofluid flow over a vertical plate with Atangana-Baleanu time-fractional integral. Alex. Eng. J. 2022, 61, 10045–10053. [Google Scholar] [CrossRef]
  29. Rashid, U.; Baleanu, D.; Iqbal, A.; Abbas, M. Shape Effect of Nanosize Particles on Magnetohydrodynamic Nanofluid Flow and Heat Transfer over a Stretching Sheet with Entropy Generation. Entropy 2020, 22, 1171. [Google Scholar] [CrossRef]
  30. Manohar, G.R.; Venkatesh, P.; Gireesha, B.J.; Madhukesh, J.K.; Ramesh, G.K. Performance of water, ethylene glycol, engine oil conveying SWCNT-MWCNT nanoparticles over a cylindrical fin subject to magnetic field and heat generation. Int. J. Model. Simul. 2021, 1–10. [Google Scholar] [CrossRef]
  31. Mathews, J.H.; Fink, K.D. Numerical Methods Using MATLAB, 3rd ed.; Prentice Hall: Hoboken, NJ, USA, 1998. [Google Scholar]
  32. Ramesh, G.; Madhukesh, J.; Das, R.; Shah, N.A.; Yook, S.-J. Thermodynamic activity of a ternary nanofluid flow passing through a permeable slipped surface with heat source and sink. Waves Random Complex Media 2022, 1–21. [Google Scholar] [CrossRef]
  33. Gangadhar, K.; Nayak, R.E.; Rao, M.V.S.; Kannan, T. Nodal/Saddle Stagnation Point Slip Flow of an Aqueous Convectional Magnesium Oxide–Gold Hybrid Nanofluid with Viscous Dissipation. Arab. J. Sci. Eng. 2021, 46, 2701–2710. [Google Scholar] [CrossRef]
  34. Bhattacharyya, S.; Gupta, A. MHD flow and heat transfer at a general three-dimensional stagnation point. Int. J. Non-Linear Mech. 1998, 33, 125–134. [Google Scholar] [CrossRef]
  35. Dinarvand, S. Nodal/saddle stagnation-point boundary layer flow of CuO–Ag/water hybrid nanofluid: A novel hybridity model. Microsyst. Technol. 2019, 25, 2609–2623. [Google Scholar] [CrossRef]
  36. Bachok, N.; Ishak, A.; Nazar, R.; Pop, I. Flow and heat transfer at a general three-dimensional stagnation point in a nanofluid. Phys. B Condens. Matter 2010, 405, 4914–4918. [Google Scholar] [CrossRef]
Figure 1. Geometry of the problem.
Figure 1. Geometry of the problem.
Mathematics 10 01185 g001
Figure 2. Numerical flow chart.
Figure 2. Numerical flow chart.
Mathematics 10 01185 g002
Figure 3. The curve f ( η ) on c .
Figure 3. The curve f ( η ) on c .
Mathematics 10 01185 g003
Figure 4. The curve g ( η ) on c .
Figure 4. The curve g ( η ) on c .
Mathematics 10 01185 g004
Figure 5. The curve θ ( η ) on c .
Figure 5. The curve θ ( η ) on c .
Mathematics 10 01185 g005
Figure 6. The curve χ ( η ) on c .
Figure 6. The curve χ ( η ) on c .
Mathematics 10 01185 g006
Figure 7. The curve f ( η ) on ϕ 2 .
Figure 7. The curve f ( η ) on ϕ 2 .
Mathematics 10 01185 g007
Figure 8. The curve g ( η ) on ϕ 2 .
Figure 8. The curve g ( η ) on ϕ 2 .
Mathematics 10 01185 g008
Figure 9. The curve θ ( η ) on ϕ 2 .
Figure 9. The curve θ ( η ) on ϕ 2 .
Mathematics 10 01185 g009
Figure 10. The curve χ ( η ) on ϕ 2 .
Figure 10. The curve χ ( η ) on ϕ 2 .
Mathematics 10 01185 g010
Figure 11. The curve θ ( η ) on S .
Figure 11. The curve θ ( η ) on S .
Mathematics 10 01185 g011
Figure 12. The curve θ ( η ) on δ .
Figure 12. The curve θ ( η ) on δ .
Mathematics 10 01185 g012
Figure 13. The curve χ ( η ) on S c .
Figure 13. The curve χ ( η ) on S c .
Mathematics 10 01185 g013
Figure 14. The curve χ ( η ) on β .
Figure 14. The curve χ ( η ) on β .
Mathematics 10 01185 g014
Figure 15. The curve χ ( η ) on E .
Figure 15. The curve χ ( η ) on E .
Mathematics 10 01185 g015
Table 1. Thermophysical properties of nanoparticles and water.
Table 1. Thermophysical properties of nanoparticles and water.
Particles ρ ( k g / m 3 ) c p ( J / k g   K ) k ( W / m K )
Fe 3 O 4 51806709.7
Go 18007175000
Water997.141790.613
Table 2. Numerical validation for various values of c in the absence of S , λ , ϕ 1 , and ϕ 2 .
Table 2. Numerical validation for various values of c in the absence of S , λ , ϕ 1 , and ϕ 2 .
Gangadar et al. [33]Bhattacharyya and
Gupta [34]
Dinarvand [35]Bachok et al. [36]Present Study
Parameter c c c c c
−0.50.5−0.50.5−0.50.5−0.50.5−0.50.5
Re x C f x 1.23021.26691.23121.26791.23251.2681-1.26811.23081.2675
Re x C f y 0.05580.49910.05570.49930.05570.4993-0.49940.05550.4990
N u x Re x 1.12271.29381.12351.33021.12371.3301-1.33021.12311.2979
Table 3. Computational values of f ( 0 ) , g ( 0 ) , θ ( 0 ) and χ ( 0 ) for different values of c .
Table 3. Computational values of f ( 0 ) , g ( 0 ) , θ ( 0 ) and χ ( 0 ) for different values of c .
c f ( 0 ) g ( 0 ) θ ( 0 ) χ ( 0 )
−0.51.127855−0.1022240.5236660.551725
−0.21.1238000.3074320.5213010.550452
0.01.1300490.5230080.5246050.565968
0.21.1406880.6987890.5301970.587699
0.51.1614750.9150780.5407860.624976
Table 4. For nodal point, computational values of f ( 0 ) , g ( 0 ) , θ ( 0 ) and χ ( 0 ) with Pr = 6.2 ,   n = 0.2 ,   α = 0.5 and ϕ 1 = 0.1 .
Table 4. For nodal point, computational values of f ( 0 ) , g ( 0 ) , θ ( 0 ) and χ ( 0 ) with Pr = 6.2 ,   n = 0.2 ,   α = 0.5 and ϕ 1 = 0.1 .
ϕ 2 S δ S c β E f ( 0 ) g ( 0 ) θ ( 0 ) χ ( 0 )
0.010.50.50.50.50.51.1614750.9150780.5407860.624976
0.02 1.1730170.9241710.5275880.631374
0.03 1.1812750.9306770.5143900.637650
0.01−0.50.50.50.50.51.1614750.9150780.6345730.627710
0.0 1.1614750.9150780.5918760.626489
0.5 1.1614750.9150780.5407860.624976
0.010.50.00.50.50.51.1614750.9150780.8570720.613593
0.5 1.1614750.9150780.5407860.624976
1.0 1.1614750.9150780.3950140.629594
0.010.50.50.10.50.51.1614750.9150780.5407860.315530
0.3 1.1614750.9150780.5407860.505093
0.5 1.1614750.9150780.5407860.624976
0.010.50.50.50.00.51.1614750.9150780.5407860.521808
0.5 1.1614750.9150780.5407860.624976
1.0 1.1614750.9150780.5407860.716755
0.010.50.50.50.50.01.1614750.9150780.5407860.703487
0.51.1614750.9150780.5407860.624976
1.01.1614750.9150780.5407860.579943
Table 5. For saddle point, computational values of f ( 0 ) ,   g ( 0 ) ,   θ ( 0 ) and χ ( 0 ) with Pr = 6.2 ,   n = 0.2 ,   α = 0.5 and ϕ 1 = 0.1 .
Table 5. For saddle point, computational values of f ( 0 ) ,   g ( 0 ) ,   θ ( 0 ) and χ ( 0 ) with Pr = 6.2 ,   n = 0.2 ,   α = 0.5 and ϕ 1 = 0.1 .
ϕ 2 S δ S c β E f ( 0 ) g ( 0 ) θ ( 0 ) χ ( 0 )
0.010.50.50.50.50.51.127855−0.1022240.5236660.551725
0.02 1.139063−0.1032400.5110810.557363
0.03 1.147081−0.1039670.4984760.563019
0.01−0.50.50.50.50.51.127855−0.1022240.6243210.554988
0.0 1.127855−0.102224Not Converging
0.5 1.127855−0.1022240.5236660.551725
0.010.50.00.50.50.51.127855−0.1022240.8148500.539718
0.5 1.127855−0.1022240.5236660.551725
1.0 1.127855−0.1022240.3804860.511333
0.010.50.50.10.50.51.127855−0.1022240.5236660.261420
0.3 1.127855−0.1022240.5236660.437491
0.5 1.127855−0.1022240.5236660.551725
0.010.50.50.50.00.51.127855−0.1022240.5236660.422603
0.5 1.127855−0.1022240.5236660.551725
1.0 1.127855−0.1022240.5236660.659084
0.010.50.50.50.50.01.127855−0.1022240.5236660.642921
0.51.127855−0.1022240.5236660.551725
1.01.127855−0.1022240.5236660.497428
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Madhukesh, J.K.; Ramesh, G.K.; Roopa, G.S.; Prasannakumara, B.C.; Shah, N.A.; Yook, S.-J. 3D Flow of Hybrid Nanomaterial through a Circular Cylinder: Saddle and Nodal Point Aspects. Mathematics 2022, 10, 1185. https://doi.org/10.3390/math10071185

AMA Style

Madhukesh JK, Ramesh GK, Roopa GS, Prasannakumara BC, Shah NA, Yook S-J. 3D Flow of Hybrid Nanomaterial through a Circular Cylinder: Saddle and Nodal Point Aspects. Mathematics. 2022; 10(7):1185. https://doi.org/10.3390/math10071185

Chicago/Turabian Style

Madhukesh, Javali K., Gosikere K. Ramesh, Govinakovi S. Roopa, Ballajja C. Prasannakumara, Nehad Ali Shah, and Se-Jin Yook. 2022. "3D Flow of Hybrid Nanomaterial through a Circular Cylinder: Saddle and Nodal Point Aspects" Mathematics 10, no. 7: 1185. https://doi.org/10.3390/math10071185

APA Style

Madhukesh, J. K., Ramesh, G. K., Roopa, G. S., Prasannakumara, B. C., Shah, N. A., & Yook, S. -J. (2022). 3D Flow of Hybrid Nanomaterial through a Circular Cylinder: Saddle and Nodal Point Aspects. Mathematics, 10(7), 1185. https://doi.org/10.3390/math10071185

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