Next Article in Journal
Identification of Linear Time-Invariant Systems: A Least Squares of Orthogonal Distances Approach
Next Article in Special Issue
Advanced Computational Framework to Analyze the Stability of Non-Newtonian Fluid Flow through a Wedge with Non-Linear Thermal Radiation and Chemical Reactions
Previous Article in Journal
An Augmented Social Network Search Algorithm for Optimal Reactive Power Dispatch Problem
Previous Article in Special Issue
Significance of Darcy–Forchheimer Law, Activation Energy, and Brownian Motion of Tiny Particles on the Dynamics of Rotating MHD Micropolar Nanofluid
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ternary Hybrid Nanofluid Flow Containing Gyrotactic Microorganisms over Three Different Geometries with Cattaneo–Christov Model

1
Department of Applied Science, Meerut Institute of Engineering and Technology, Meerut 250 005, Uttar Pradesh, India
2
Department of Mathematics, Graphic Era Deemed to be University, Dehradun 248 002, Uttarakhand, India
3
Department of Mechanical Engineering, Sejong University, Seoul 05006, Republic of Korea
4
Department of Mathematics, Statistics and Computer Science, G. B. Pant University of Agriculture and Technology, Pantnagar 263 145, Uttarakhand, India
5
Center of Research, Faculty of Engineering, Future University in Egypt, New Cairo 11835, Egypt
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work and are co-first authors.
Mathematics 2023, 11(5), 1237; https://doi.org/10.3390/math11051237
Submission received: 1 February 2023 / Revised: 27 February 2023 / Accepted: 1 March 2023 / Published: 3 March 2023
(This article belongs to the Special Issue Computational Fluid Dynamics II)

Abstract

:
The movement of microorganism cells in fluid influences various biotic processes, including septicity and marine life ecology. Many organic and medicinal applications need to look into the insight of mechanism in nanofluids containing a microbial suspension. The current paper concerns the bioconvection of a ternary hybrid nanofluid (Al2O3-Cu-CNT/water) flow containing motile gyrotactic microorganisms toward three different geometries (a flat plate, a wedge, and a cone) in the occurrence of natural convection, radiation, and heat source/sink. The Cattaneo–Christov theory is employed to develop the model. The equations are solved by using the “bvp4c function in MATLAB”. The influence of the crucial significant factors on the motile microorganisms’ density, velocity, temperature, nanoparticles’ concentration, microbe density gradient, and transmission rates of heat and mass is discussed. The results depict that the heat transmission rate is highest for the flow toward the cone, whereas the mass transmission rate and microbe density gradient are highest for the flow toward the wedge. In addition, the higher estimates of the thermal relaxation parameter corresponding to the Cattaneo–Christov theory act to enhance the rate of heat transmission. The results of the current study will be useful to many microbial-enhanced oil recovery systems, carriage processes, architectural design systems, medicinal fields that utilize nanofluids, and so on.

1. Introduction

The importance of NF and HNF flow over a wedge, a cone, and a flat plate is due to their vast applications. These applications are seen in, but are not restricted to, fiber technology, the design of vessels for nuclear waste retention, steam generators, spacecraft design, and solar power collectors. Initially, for the aforementioned applications, manufacturers utilized lubricants, air, and water to accomplish heat transmission (HT). However, such natural materials do not fulfill the standards at the industrial scale. Many working fluids have poor thermal conductivity (TC), which is considered a hurdle in HT, and it slows down the working of thermal flow models. In the present world, the production of several gadgets and apparatuses utilized in industrial and technical arenas has evolved significantly. For instance, in various devices, temperature rises over time, owing to electrical resistance. The heat-transport capacity of such devices is diminished due to electric resistance, resulting in a technical defect. To lessen the possibility of a technical malfunction, many devices and components must be heat-intolerant. Therefore, these applications require a coolant that has superior TC and heat transfer qualities. Previous research (Choi [1]) has demonstrated that dispersing nanoparticles (NPs) in a base fluid dramatically changes the characteristics and heat transfer qualities of a convectional base fluid (BF). Furthermore, if necessary, the characteristics of nanofluids (NFs) can be tailored to a specific purpose. The main limitation with an NF containing a single type of NPs is that the NPs have better thermal linkages or good rheological features, but not both. The properties of NFs may be enhanced by altering the nanoparticles volume fraction (NVF); however, this has a restriction because of the difficulties arising owing to viscosity escalation. To erase the limitation of NFs, researchers studied hybrid nanofluids (HNFs) that aim to elevate the properties of nanoparticles of different materials. Early HNFs studies employed NFs with two different particles; however, in the last three years, a distinct research trend has evolved that consists of mixing several NPs to generate an ideal nanofluid: three NPs in a working fluid. Such NFs are mentioned as “ternary nanofluids”, “ternary hybrid nanofluids” (THNF), or “trihybrid nanofluids.” The authors of the present study have studied the comparative behavior of bioconvective THNF flow over three different geometries (a flat plate, a wedge, and a cone). The main aim of the authors is to investigate the HT rate of THNF (Al2O3-Cu-CNT/water) over the surface of three different geometries (a flat plate, a wedge, and a cone).
In their review article, Sahoo [2] and Adun et al. [3] emphasized the shapes of NPs, the material of NPs, the size of NPs, and the proportion of NPs in optimizing the properties of THNF. Furthermore, selecting the appropriate NPs mixture, as well as optimal NPs mixing and stability, are critical variables that may influence THNF development, notably for HT applications (Xuan et al. [4]). The invention of THNFs and the experimental validation of correlations (Sahoo and Kumar [5] and Sahoo [6]) have paved the way for researchers working in numerical flow analysis to investigate the possibility of THNFs as a viable coolant in industrial HT applications. Animasaun et al. [7] expounded the significance of convective heating on the transport of THNF over a horizontal surface stretching with a linear rate. They investigated the flow with water as a BF and aluminum oxide, aluminum, and silver NPs, all with different shapes. They concluded that convective heating could be utilized for achieving a high HT rate. Alanazi et al. [8] recently studied the dynamics of ternary-hybrid nanofluid subject to magnetic flux density and heat source or sink on a convectively heated surface. They deduced that increasing thermo-migration of NPs causes HT to reduce. Raju et al. [9] explicated the nonlinear movements of axisymmetric ternary hybrid nanofluids in a thermally radiated expanding or contracting permeable Darcy Walls with different shapes and densities by using simple linear regression. Ramzan et al. [10] expounded on the influence of ion slip and Hall current on kerosene-based THNF flow towards a surface rotating about an axis with silver, copper, and graphene oxide NPs.
Improvement of thermal control is a concern of manufacturers of the electrical division, thermal systems, power sector, and medicinal bids. To study the HT rate in any medium, it is of utmost importance to see the regulations of HT in that particular medium. For a long period, HT was analyzed with the Fourier law [11], which says that HT occurs at an infinite speed. Later, Cattaneo [12] modified the law of Fourier by generalizing it and stating that HT occurs at a restricted speed in any medium. He incorporated a time relaxation factor to generalize the Fourier model, and it reduced the nature of the parabolic thermal equation to a hyperbolic equation. Christov [13] achieved the material invariant form of Cattaneo’s law through the Oldroyd derivative. After the successful formulation by Cattaneo and Christov, the model came to be known as the Cattaneo–Christov heat flux model (CCHM). Numerous researchers have calculated the HT rate of different flows with a time lag factor using CCHM. Venkateswarlu et al. [14] expounded on the significance of CCHM in nanofluid flow past a surface with a linear stretching rate. They discussed the nanofluid flow with BF as water, discussed the flow pattern, and presented a comparison of MoS2/water and MgO/water nanofluid. Rawat and Kumar [15] analyzed Cu-water NF flow over a stretching surface with CCHM and inferred that time lag presence reduces the HT rate. Ramzan et al. [16] expounded on the flow behavior of dusty micropolar NF flow over a sheet with CCHM. Lv et al. [17] studied the significance of CCHM in Reiner–Rivlin NF flow over a disk rotating about an axis with bioconvection phenomena. Abderrahmane et al. [18] elucidated the impact of MHD Hybrid Nanofluid Mixed Convection Heat Transfer and Entropy Generation in a 3-D Triangular Porous Cavity with Zigzag Wall and Rotating Cylinder.
Another intriguing study field is the bioconvection phenomenon, which covers a range of applications in the real world. The directed movement of diverse microorganisms functions as the base for many bioconvection models. Different species are distinguished by the direction in which these microorganisms swim. Gyrotactic microbes (GM) are the ones that travel upwards against gravity in stagnant water, resulting in the higher section of the fluid being denser than the lower. When the higher portion gets too dense, owing to microbe accumulation, it becomes unstable, and microbes fall, causing bioconvection. This bioconvection pattern is maintained by returning swimming microorganisms. Applications of this phenomenon include biological tissues, biofuels, biotechnology, microsystems, transportation processes, enzyme biosensors, microbial-enhanced oil recovery, and so on. Furthermore, bioconvective nanofluid has vast applications, including the development of nanoparticle processing, architectural design, nanoparticle stability, polymer coating, automobile coolants, and sterilizing applications in thermal sciences, among others. The novel idea of nanofluid with different nanoparticles are presented in [19,20,21]. Shi et al. [22] explicated the idea of a magnetic bioconvective stream of a Cross nanoliquid over a stretched sheet. They determined that the growing Lewis number causes the density of gyrotactic microorganisms (GM) to fall. Al Khaled and Khan et al. [23] expounded on the thermal properties of a Casson nanoliquid flow containing GM over a surface in motion. Bhatti et al. [24] explicated the applications of GM in blood-based non-Newtonian nanoliquid flow through a narrow artery. They utilized the Sutterby fluid model to study blood flow. Kairi et al. [25] expounded on the significance of the suspension of GM in a Casson nanoliquid flow past an inclined, stretched sheet with Marangoni convection. Ali et al. [26] explicated the magnetic Cross nanomaterial flow over an extending horizontal plate and horizontal cylinder and studied the significance of GM and CCHM in their flow model.
Of late, several researchers have published papers on the flow of NFs or HNFs over a cone, a wedge, and a flat plate (Mishra and Kumar [27] and Yaseen et al. [28], Mahanthesh and Mackolil [29], and Gumber et al. [30]). In this research communication, authors have modeled and investigated the bioconvective THNF flow (Al2O3-Cu-CNT/water flow) over three different geometries (a flat plate, a wedge, and a cone) with motile gyrotactic microorganisms (GM). Moreover, the significance of CCHM, heat generation/absorption, natural convection, and thermal radiation are also studied. Some past research (see, for instance, Rawat et al. [31], Sandeep and Reddy [32], and Reddy et al. [33]) has conducted such comparative studies. However, until now, no author(s) have attempted or published a study, which shows the comparative behavior of bioconvective THNF (Al2O3-Cu-CNT/water) flow over three different geometries (a flat plate, a wedge, and a cone). The results of the current study will be useful to many microbial-enhanced oil recovery systems, carriage processes, architectural design systems, medicinal fields that utilize nanofluids, and so on. The novelty and aim of this research communication and the aspects this paper attempts to explore are as mentioned in the following points:
  • The mathematical flow model of THNF (Al2O3-Cu-CNT/water) over three different geometries (a flat plate, a wedge, and a cone).
  • The application of CCHM in the THNF flow.
  • The significance of motile GM in THNF (Al2O3-Cu-CNT/water).
  • Investigating the HT rate of THNF flow over three different geometries (a flat plate, a wedge, and a cone) and finding the condition/geometries under present modeling for which THNF has the maximum HT rate.
  • Comparison of the HT rate of THNF flow with the HT rate of HNF and NF in the case of all three geometries.

2. Flow Model and Governing Equations

2.1. Flow Assumptions and Mathematical Model

This research communication explores the flow of THNF (Al2O3-Cu-CNT/water) toward three different geometries (a cone, a wedge, and a flat plate) containing motile GM. This section of the paper deals with the modeling of the flow of THNF in the boundary layer region (BLR) and the assumptions of the model. To study the THNF flow, the surface of a cone, wedge, and flat plate are assumed coincidental to the x-axis, and the y-axis is normal to it (refer to Figure 1). The flow is instigated by the manifestation of the convection effect. The half angle of the wedge and cone are denoted by α, and the radius of the cone is denoted by r. To overview the influence of the presence of time lag in the heat transfer (HT) and mass transfer (MT) rate, the authors modeled the flow with the Cattaneo–Christov theory. Furthermore, heat source/sink, thermal radiation, and injection/suction effects are also deliberated together to study the flow of THNF in the present problem. In the last term on the right hand side in the energy equation (Equation (3)), q r denotes the thermal radiative heat flux, and it is defined later, after governing equations.
Based on the assumptions and aforesaid discussion, the bioconvective flow of THNF (Al2O3-Cu-CNT/water) toward three different geometries (a flat plate, a wedge, and a cone) is represented by the following PDEs [34,35]:
Continuity equation:
( r n u ) x + ( r n v ) y = 0
Since the thermal BLR thickness is very minuscule, the r = x s i n ( α ) .
Momentum equations:
u u x + v u y = ν t h n f ( 2 u y 2 ) + ( ρ β ) t h n f ρ t h n f g ( T T ) cos α
Energy equation:
( ρ C p ) t h n f { u T x + v T y + τ t ( u 2 2 T x 2 + ( u u x + v u y ) T x + ( u v x + v v y ) T y + 2 u v 2 T x y + v 2 2 T y 2 Q o ( p C p ) t h n f ( u T x + v T y ) ) } = k t h n f 2 T y 2 + Q 0 ( T T ) q r y
Concentration equation:
u C x   +   v C y   + τ c [ u 2 2 C x 2 + ( u u x + v u y ) C x + ( u v x + v v y ) C y + 2 u v 2 C x y + v 2 2 C y 2 ] = D n 2 C y 2
Microorganisms’ equation:
u N x + v N y + ( b W c C w C ) ( y ( N C y ) ) = D m 2 N y 2
where “ ( u , v ) are the velocity components along ( x , y ) directions”, respectively. In addition, the variables used in Equations (1)–(6) have the following meaning: “T—temperature”, “C—concentration”, “N—motile microorganisms concentration”, “Q0—heat source/sink coefficient”, “b—chemotaxis constant”, “Wc—swimming cells speed”, “ τ t —thermal relaxation time”, “ τ c —solutal relaxation time”. Furthermore, “k—thermal conductivity, μ —dynamic viscosity, ν —kinematic viscosity, C p —heat capacity, ρ —density, σ —electrical conductivity, Dn—Brownian diffusion coefficient, and Dm—microorganisms diffusion”. Moreover, subscript w—wall condition, —free stream condition, thnf—ternary hybrid nanofluid, nf—nanofluid, and bf—fluid.
Furthermore, at the surface, suction/injection is applied with vw as velocity of mass transfer, and surface is moving with initial velocity, as u w = ν f x / l 2 . The temperature at the surface T w is taken as a variable T w = T + a x s , where T is the free stream temperature of THNF, and s is the parameter related to surface temperature. In addition, C w and N w are concentrations of NPs and GM at the surface, whereas C and N are concentrations of NPs and GM at the free stream.
The Boundary Conditions (BCs) are:
u = u w = ν f   x l 2 ,         v = v w ,         T = T w = T + a x s ,     C = C w ,       N = N w   at         y   =     0 u 0 ,                           T T ,   C C ,           N N   as         y }
The Rosseland theory is used to approximate the thermal radiative heat flux q r , and it is demarcated as follows [15]:
q r = 4 σ * 3 k * T 4 y
where “k* is the mean absorption coefficient”, and “σ* is the Stefan-Boltzmann constant”. In Equation (7), the value of q r is estimated by assuming T 4 4 T 3 T 3 T 4 with the help of expansion via Taylor’s series. This assumption is valid for flow that has minuscule temperature differences in the layers. To arrive at such an assumption, higher-order terms are neglected in Taylor’s series expansion about T . The energy equation (Equation (3)) is further solved with the help of Equation (7) and the aforementioned assumption of T 4 .
The flow of THNF (Al2O3-Cu-CNT/water) toward three different geometries (a cone, a wedge, and a flat plate) is characterized by the following conditions [32,35]:
(a)
THNF flow toward cone: n = 1 and α ≠ 0.
(b)
THNF flow toward wedge: n = 0 and α ≠ 0.
(c)
THNF flow toward flat plate: n = 0 and α ≠ 0.

2.2. Properties of Ternary Hybrid Nanofluid

The current subsection refers to the modeling and investigation of THNF. The characteristics of THNF are subjected to the base fluid (water), NPs, and the NVF of Al2O3, Cu, and CNT NPs, as presented in Table 1 [15,36]. Moreover, the authors investigated the effect of differently shaped nanoparticles (CNT—cylindrical, Cu—platelet, and Al2O3—spherical) on the thermal behavior of the THNF in this study.
Table 2 lists the models used to develop thermophysical correlations of the properties of THNF [36,37]. These correlations are valid and provide precise results when the NVF is taken between the range of 0.01 to 0.1 and temperature range is between 35 °C to 50 °C. In this paper, NVF of Al2O3, Cu, and CNT NPs is denoted by φ 1 , φ 2 , and φ 3 . Furthermore, the following subscripts are used for: sp1—Al2O3, sp2—Cu, and sp3—CNT NPs.

2.3. Conversion of the Model Equations with Similarity Transformation

The flow governing PDEs (Equations (1)–(6)) are first modified using the similarity transformation [32,38]:
η = y l ,       u = ν f x l 2 f ( η ) ,           v =   ( n + 1 ) ν f f ( η ) l ,     θ ( η ) = T T T w T ,   ϕ ( η ) = C C C w C ,     ζ = N N N w N
The continuity equation (Equation (1)) is identically satisfied, and the same can be verified by utilizing the similarity transformations from Equation (7) in the continuity equation (Equation (1)). The reduced form of Equations (2)–(5) and BCs (Equation (6)) after using similarity transformations are as follows:
Reduced Momentum equation:
f 2   ( n + 1 ) f f = 1 2 f + 3 2 λ θ cos α
Reduced Energy equation:
s f θ   ( n + 1 ) f θ + γ t [ s ( s 1 ) f 2 θ + s ( θ f 2   ( n + 1 ) θ f f ) + ( n + 1 ) 2 θ f f 2 ( n + 1 ) s f f θ +   ( n + 1 ) 2 f 2 θ Q 4 ( s θ f   ( n + 1 ) f θ ) ] = 1 4 Pr ( 5 + 4 3 R d ) θ + Q 4 θ
Reduced Concentration equation:
ϕ + L e ( n + 1 ) f ϕ   L e ( n + 1 ) 2 γ c { f f ϕ + f 2 ϕ } = 0  
Reduced Microorganisms’ equation:
ζ P e ( ζ ϕ + ( ϖ + ζ ) ϕ ) + ( n + 1 ) L b f ζ = 0
Reduced BCs:
f ( 0 ) = 1 ,         f ( 0 ) = S u ( n + 1 ) ,         θ ( 0 ) = 1 ,     ϕ ( 0 ) = 1 ,       ζ ( 0 ) = 1     at               η   =   0 f ( 0 ) = 0 ,                           θ ( 0 ) = 0 ,           ϕ ( 0 ) = 0 ,                     ζ ( 0 ) = 0                   as         η }
In the aforementioned Equations (9)–(12) and BCs (Equation (13)), the nondimensional parameters are encountered, and they are named and defined as follows:
λ ( = G r x / R e x 2 ) —Natural convection parameter, where G r x ( = g β f ( T w T ) l 4 / x ν f 2 ) —Local Grashof number, R e x = ( u w x / v f ) —Local Reynolds number. Furthermore, P r ( = ( ρ c p ) f ν f / k f ) —Prandtl number, R d = ( 4 σ * T 3 / k f k * ) —Thermal radiation parameter, Q ( = Q 0 l 2 / ( ρ C p ) f ν f   ) —Heat source/sink parameter, and γ t = ( τ t ν f / l 2 ) —Thermal relaxation parameter. In addition, L e ( = v f / D n ) —Lewis number, L b ( = v f / D m ) —bioconvection Lewis number, P e ( = b W c / D m ) —bioconvection Peclet number, ϖ ( = N 2 / ( N 1 N 2 ) ) —bioconvection constant, γ c = ( τ c ν f / l 2 ) —Solutal relaxation parameter, and S u = ( v w l / v f ) —Suction/Injection parameter. Furthermore, 1 = μ t h n f / μ f ,   2 = ρ t h n f / ρ f ,   3 = ( ρ β ) t h n f / ( ρ β ) f ,   4 = ( ρ C p ) t h n f / ( ρ C p ) f ,   5 = k t h n f / k f are constants based on THNF correlations and properties (see Table 2)”.
To have similar solutions of the Equations (9)–(12) with the BCs (Equation (13)), the term of x must vanish from these equations. It is noticed that the natural convection parameter λ ( = G r x / R e x 2 ) appearing in Equation (9) when computed in the simplified form will be written as λ = g β f ( T w T ) l 4 / ν f 2 x . Further, the BCs (Equation (6)) imply that surface temperature is in a variable form, as T w = T + a x s . So, further simplification will lead to the value of the natural convection parameter as λ = g β f a x s 1 l 4 / ν f 2 . Therefore, a similar solution will exist only when s = 1 because the term of x will vanish. For a similar solution to exist and for the computations of the results of this study, the value of the parameter related to surface temperature is strictly restricted to s = 1 .

3. Engineering Parameters

The analysis of the concerned and related engineering parameters is necessary while studying a flow model because it helps in applying the bids of the flow model in real-life applications. The “Nusselt number Nu” (representing heat transfer rate), “Sherwood number Sh” (representing mass transfer rate), and motile microorganism density gradient (Nm), at the surface of the flat plate, wedge, and cone, respectively, are [24,39]:
N u = ( 5 + R d ) θ ( 0 )
S h = φ ( 0 )
N m = ζ ( 0 )

4. Methodology of Numerical Approach

The current section demonstrates the approach taken to solve the model depicted by similarity equations (Equations (9)–(12)) and BCs (13). The transfigured equations are nonlinear and coupled. Deducing an analytical solution to such a system is not an easy task. Therefore, the authors attempted to solve the system (Equations (9)–(12)) and BCs (13) numerically. The numerical solution of the Equations is derived with the “bvp4c function (a built-in package in MATLAB)”. The “bvp4c function” uses a finite difference scheme together with a precision of fourth-order with the help of the “3-stage Lobatto IIIA formula”. The step-by-step working of bvp4c function is provided by Shampine et al. [40]. To derive the solution of the present problem, the Equations (9)–(12) are transformed into ODEs of first-order with the help of following transformation:
( 1 ) = f ,   ( 2 ) = f ,   ( 3 ) = f ,   ( 4 ) = θ ,   ( 5 ) = θ ,   ( 6 ) = φ ,   ( 7 ) = φ , ( 8 ) = ζ , ( 9 ) = ζ
Utilizing the new variables, the following MATLAB syntax is used:
1 = 2 1 ( ( 2 ) 2   ( n + 1 ) ( 1 ) ( 3 ) 3 2 λ ( 4 ) cos α ) ; 2 = Q 4 ( 4 ) ( 2 ) ( 4 ) +   ( n + 1 ) ( 1 ) ( 5 ) γ [ ( 4 ) ( 2 ) 2   ( n + 1 ) ( 4 ) ( 1 ) ( 3 ) + ( n + 1 ) 2 ( 5 ) ( 1 ) ( 2 ) 2 ( n + 1 ) ( 1 ) ( 2 ) ( 5 ) Q 4 ( ( 4 ) ( 2 )   ( n + 1 ) ( 1 ) ( 5 ) ) ] 1 4 Pr ( 5 + 4 3 R d γ   ( n + 1 ) 2 ( 1 ) 2 ) ; 3 = L e ( n + 1 ) ( 1 ) ( 7 )   L e ( n + 1 ) 2 γ c ( 1 ) ( 2 ) ( 7 ) ( 1 + L e ( n + 1 ) 2 γ c ( 1 ) 2 ) ; 4 = P e ( ( 9 ) ( 7 ) + ( ϖ + ( 8 ) ) ( 3 ) ) ( n + 1 ) L b ( 1 ) ( 9 ) ;
The following MATLAB syntax is used for BCs at the surface and far-field:
( 2 ) ( 0 ) = 1 ,         ( 1 ) ( 0 ) = S u ( n + 1 ) ,         ( 4 ) ( 0 ) = 1 ,   ( 6 ) ( 0 ) = 1 ,       ( 8 ) ( 0 ) = 1   at         η   =   0 ( 2 ) ( η ) = 0 ,                         ( 4 ) ( η ) = 0 , ( 6 ) ( η ) = 0 ,               ( 8 ) ( η ) = 0   as         η }
In addition, the values of slopes not known are guessed to begin the process of solving the equations. The iteration process is repeated, and the solution is acknowledged only when BCs are met (see Figure 2). In this model, ηmax was taken as 6, as the conditions at this point were satisfied asymptotically, and step size was taken as 0.01. To authenticate the results of the present model, a comparison is outlined with the published results of Vajravelu and Nayfeh [38] (see Table 3). The results are reasonably consistent, confirming that the present conclusions are effective.

5. Results and Discussion

This communication explores the flow of THNF (Al2O3-Cu-CNT/water) toward three different geometries (a flat plate, a wedge, and a cone) containing motile GM. This section of the paper deals with the analysis of results and their discussion of THNF flow toward the three different geometries (a flat plate, a wedge, and a cone). The discussion will be focused on the velocity f ( η ) , temperature θ ( η ) , NPs concentration φ ( η ) , and microorganism concentration profile ζ ( η ) inside the BLR. The “bvp4c function in MATLAB” was operated to deduce the numerical solution of this model, and various properties of the flow were analyzed with the solution deduced. During the derivation of the solution, the following general values of the parameters (present in Equations (18) and (19)) were used for computations: Pr = 6.2 (water), Rd = 10, Su = 0.1, Pe = 0.1, Ω = 0.2, γ t = 1.5 , γ c = 0.01 , Le = 1.5, Lb = 0.5, λ = 0.2, Q = −2, s = 1 , and α = π 4 (for cone and wedge).

5.1. Discussion of Velocity Profiles, Velocity Boundary Layer Patterns, and Streamlines

Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7 show the influence of the natural convection parameter λ , NVF ϕ i , and suction/injection parameter Su on the velocity f ( η ) , respectively. Figure 3 displays that the velocity f ( η ) of THNF rises with an increase in the natural convection parameter λ . The natural convection phenomenon is the mechanism of transportation of heat. This phenomenon causes THNF motion because of variation in fluid density, which occurs due to temperature gradients. The applications of natural convection can be seen in food and chemical industries and also in solar applications. In the present model, the increasing value of the parameter λ causes the buoyancy forces to increase and helps the velocity increase in all three cases (i.e., THNF flow toward a cone, a wedge, and a flat plate). Figure 4, Figure 5 and Figure 6 display the velocity f ( η ) of THNF with a rise in the NVF ϕ i of the NPs. The increasing NVF ϕ 1 of Al2O3 nanoparticles in the working fluid causes the velocity f ( η ) to decrease. In addition, the increasing NVF ϕ 2 and ϕ 3 of Cu and CNT nanoparticles results in an increase in velocity f ( η ) . The decrease in velocity due to an increase in NVF of Al2O3 nanoparticles is due to the spherical shape of these NPs, as the spherical shape has a larger surface area for interacting with other NPs, and due to this, the effective viscosity of THNF increases. Consequently, the velocity f ( η ) of THNF decreases with an increment in the NVF of Al2O3 nanoparticles. Figure 7 displays that the velocity f ( η ) of THNF drops with a surge in the numerical value of the suction/injection parameter Su. The application of suction/injection at the BLR has a central role in the boundary layer separation and is often used to delay the boundary layer separation.
The increasing suction at the surface causes the thickness of the momentum boundary layer to reduce; consequently, the layers of the flow are strongly held toward the surface destruct the momentum, and hence, velocity decreases.
In Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7, it is perceived that the velocity of THNF is highest for the flow toward the plate and is lowest for the flow toward the cone. The flat plate provides a plane and a smooth surface for the motion of the THNF, whereas the curved surface of the cone provides obstruction to the flow, and hence, for the cone, the velocity is lowest. The above results are also validated by the velocity boundary layer patterns and streamlines drawn in Figure 8 and Figure 9. Figure 8 and Figure 9 show the velocity boundary layer configuration and streamlines for the THNF flow toward the plate, wedge, and cone. The velocity boundary layer configuration in Figure 8 show that the flow toward the cone has lower velocity in comparison to the flow toward the plate and wedge. In addition, streamlines in Figure 9 also depict the same. It can be seen that for THNF flow toward the plate and wedge, the streamlines are narrow, whereas for the THNF flow toward the cone, the streamlines are wider. The wider streamlines depict that the velocity is low for the particular case.

5.2. Discussion of Temperature Profiles

Figure 10, Figure 11, Figure 12, Figure 13, Figure 14 and Figure 15 show the influence of heat source/sink parameter Q, NVF ϕ , thermal radiation parameter Rd, suction/injection parameter Su, thermal relaxation parameter γ t , and natural convection parameter λ on the temperature θ , respectively. Figure 10 displays that the temperature θ of the THNF is ominously affected by the heat source/sink parameter Q values. The temperature θ increases on increasing the value of parameter Q. The negative (Q < 0) and positive (Q > 0) values of parameter Q depict the mechanism of a heat sink and source, respectively.
The results depict that on amplifying the magnitude of the heat sink in the model, the temperature θ of THNF flow falls. The reason for the falling temperature is that the heat sink absorbs the heat produced in the BLR during the flow, which results in a falling temperature. Figure 11, Figure 12 and Figure 13 display the temperature θ of THNF with a rise in the NVF ϕ i of the NPs. The increasing NVF ϕ 1 of Al2O3 nanoparticles in the working fluid causes the temperature θ to increase. In addition, the increasing NVF ϕ 2 and ϕ 3 of Cu and CNT nanoparticles results in a decrease in temperature θ . The increase in temperature due to an increase in NVF of Al2O3 nanoparticles is due to the spherical shape of these NPs, as the spherical shape has a larger surface area for interacting with other NPs, and due to this, the effective viscosity of THNF increases and a greater amount of heat is released. Consequently, the temperature θ of THNF increases with an increment in the NVF of Al2O3 nanoparticles. Figure 14 displays that the temperature θ of THNF surges with intensification in the radiation parameter Rd. This means that the radiation provides energy to the fluid particles, which surges their movement, and hence, the temperature rise. Furthermore, it is also seen that the temperature θ of the THNF is considerably affected by the radiation parameter Rd, which confirms its crucial role in the heat transmission process.
Figure 15 displays that the temperature θ of THNF falls with a rise in the numerical value of the suction/injection parameter Su. The increasing suction at the surface causes the width of the momentum boundary layer to reduce; consequently, the layers of the flow with zero momentum are provided energy, and they gain momentum and start their motion again. This process destroys the flow and reduces the friction at the boundary layer due to layers with zero momentum. Consequently, the temperature falls. Figure 16 displays that the temperature θ of THNF falls for higher estimates of the thermal relaxation parameter γ t . The nonzero parameter γ t denotes lag in the heat conduction. The zero value of the parameter γ t resembles the traditional Fourier’s law for heat conduction. The higher values of the parameter γ t resemble the increased amount of time during the heat conduction. This revelation character of the Cattaneo–Christov model differentiates it from the “Fourier law of heat conduction”. The results indicate that increasing time lag during the heat transfer leads to lower temperature. Figure 17 displays that the temperature θ of THNF decreases with an increment in the natural convection parameter λ . The natural convection phenomenon is the mechanism of transportation of heat. This phenomenon causes THNF motion because of variation in fluid density, which occurs due to temperature gradients. The increasing temperature gradients cause the temperature to fall.
In Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15, Figure 16 and Figure 17, it is perceived that the temperature θ of THNF flow is highest for the flat plate and is lowest for the flow toward the cone. The flat surface of the plate allows the smooth motion of the flow and the increasing interaction of the NPs and working fluid particles, causing the system to interact properly, and consequently, the temperature rises. However, the curved surface of the cone obstructs the movement and the interaction of particles, and hence, the temperature is comparatively low for the flow toward the cone.

5.3. Discussion of Nanoparticles Concentration and Microorganisms Concentration

Figure 18, Figure 19 and Figure 20 show the influence of the Lewis number Le, suction/injection parameter Su, and solutal relaxation parameter γ c on the concentration ϕ of NPs, respectively. Figure 18 shows the fall in the concentration of NPs with a rise in the Lewis number Le. The Lewis number Le is negatively correlated with molecular diffusivity. The increment in the Lewis number Le is seen as the lower molecular diffusivity. The lower diffusivity of the NPs in the BLR causes a fall in the concentration of NPs. Figure 19 displays that the NPs concentration is negatively correlated to the suction/injection parameter. Stronger suction at the surface destroys the flow and momentum, reducing the diffusivity of the NPs. Hence, the concentration falls due to the stronger suction.
Figure 20 depicts the influence of the solutal relaxation parameter γ c on concentration of NPs, and γ c 0 signifies the manifestation of lag in mass transfer. The zero value of parameter γ c corresponds to the traditional Fick’s law for mass diffusion. The higher values of the parameter γ c resemble the increased amount of time during the mass diffusion. Results depict that the higher estimates of parameter γ c lead NPs’ concentration to rise.
Figure 21, Figure 22, Figure 23 and Figure 24 depict the fluctuation of the microorganism concentration profile ζ ( η ) with the bioconvection Lewis number Lb, bioconvection Peclet number Pe, bioconvection constant ϖ , and suction/injection parameter Su. Figure 21, Figure 22 and Figure 23 show that the microorganism concentration profile ζ ( η ) falls with higher estimates of the bioconvection Lewis number Lb, Peclet number Pe, and bioconvection constant. The bioconvection Lewis number Lb is in charge of the system’s bioconvection. It increases heat diffusivity at the surface, lowering the density of microbes. It also diminishes the width of the boundary layer, which is proportional to the density of the motile microbe. Figure 24 shows that a rise in the suction/injection parameter Su causes the microorganism concentration profile ζ ( η ) to decrease. A stronger effect of the suction/injection parameter reduces the width of momentum and the NPs’ concentration boundary layer; as a result, the movement of layers and particles reduces, and it further causes the obstruction in the movement of the microorganism. Consequently, the microorganism concentration profile ζ ( η ) falls with increasing application of suction at the surface. In Figure 18, Figure 19, Figure 20, Figure 21, Figure 22, Figure 23 and Figure 24, it is observed that the concentration of NPs and microorganism concentration is highest for the flow toward the flat plate and is lowest for the flow toward the cone. The flat surface of the plate allows the smooth motion of the flow, and more and more diffusion of the NPs occurs; similarly, it also provides the smooth functioning of the microorganisms’ movement. However, the curved surface of the cone obstructs the diffusion of NPs and microorganisms’ movement; hence, these two profiles are comparatively low for the flow toward the cone.

5.4. Discussion of Nusselt Number, Sherwood Number, and Motile Microorganisms’ Density Gradient

Table 4 and Table 5 display the influence of flow parameters on the Nusselt (Nu) number depicting HT rate, Sherwood (Sh) number depicting mass transmission rate (MT rate), and microorganisms’ density gradient (Nm) at the surface. Figure 25 is drawn to present the variation in Nusselt number (Nu). The higher estimates of the natural convection parameter λ cause the HT rate to rise, whereas the MT rate and microorganisms’ density gradient fall. The rising natural convection phenomenon causes the HT rate to rise due to the transportation of heat. This phenomenon causes THNF motion because of variation in fluid density, which occurs due to temperature gradients. The increasing magnitude of the suction parameter Su and heat sink parameter Q (Q < 0) acts to increase the HT rate. The magnification in the HT rate is due to a reduction in thermal BLR thickness. The reduced thickness allows more transmission of heat through the surface. On the contrary, the increasing suction parameter causes the MT rate and microorganisms’ density gradient (Nm) to decrease. Higher estimates of the radiation parameter Rd increase the HT rate. This means that the radiation provides energy to the particles, which surges their movement and, hence, the interaction with other particles, causing the HT rate to rise.
On increasing the NVF φ 2 and φ 3 of Cu and CNT nanoparticles, the HT rate rises, but on the contrary, on increasing the NVF φ 1 of Al2O3 nanoparticles, the HT rate decreases. This is because of the fact that when increasing the NVF of Cu and CNT nanoparticles, thermal boundary layer thickness decreases, whereas on increasing the NVF of Al2O3 nanoparticles, thermal boundary layer thickness increases. The reduction in thermal boundary layer thickness causes a greater amount of heat to transfer through the surface. Consequently, the HT rate increases. The amplification of the thermal relaxation parameter γ t causes the HT rate to increase at the surface. This result implies that the higher time lag present during the transmission causes an enhanced HT rate. Furthermore, higher solutal relaxation parameter γ c acts to increase the MT rate and microorganisms’ density gradient (Nm). This implies that a higher time lag during the mass transmission causes an enhanced MT rate, and it also helps the microorganisms’ density to increase. The rising Lewis number acts to decrease the MT rate and microorganisms’ density gradient (Nm). In addition, bioconvection Peclet number Pe, bioconvection constant ϖ , and bioconvection Lewis number Lb are negatively correlated with the microorganisms’ density gradient (Nm). The results depict that the heat transmission rate is highest for the flow toward the cone, whereas the mass transmission rate and microbe density gradient are highest for the flow toward the wedge.

6. Conclusions

This paper deals with the bioconvective 2D THNF flow (Al2O3-Cu-CNT/water flow) containing a microbial suspension toward three different geometries (a flat plate, a wedge, and a cone) with the suction/injection effect. The Cattaneo–Christov theory is used to analyze HT rate and MT rate. The THNF flow is modeled bearing in mind the effects of natural convection, radiation, and heat source/sink. The in-built function “bvp4c in MATLAB” is used to solve the equations, and the solutions are shown in graphs. The key findings of the study are:
A stronger application of suction causes the thickness of the momentum boundary layer to reduce.
The temperature of THNF increases with higher radiation parameter and heat source/sink parameter.
The increasing value of the thermal relaxation parameter corresponding to the Cattaneo–Christov theory acts to enhance the heat transmission rate.
The microorganism concentration profile decreases with higher bioconvection Lewis number.
The heat transmission rate is highest for the flow toward the cone.
Mass transmission rate and microbe density gradient are highest for the flow toward the wedge.

Future Scope of Research

The present study discusses the flow of ternary hybrid nanofluid without the nanoparticle aggregation effect. The study can be extended for nanoparticle aggregation effect with proper modeling of fluid flow by utilization of validated thermophysical correlations.

Author Contributions

Conceptualization: M.Y., S.K.R., N.A.S., M.K. and S.M.E.; Methodology: M.Y., S.K.R. and N.A.S.; Writing of first draft: S.K.R. and N.A.S.; Review and editing of draft: M.K. and S.M.E. All authors have read and agreed to the published version of the manuscript.

Funding

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

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

Thermal conductivityTC
Heat transmissionHT
NanoparticlesNPs
Base fluidBF
NanofluidsNFs
Nanoparticles volume fractionNVF
Hybrid nanofluidsHNFs
Ternary hybrid nanofluidTHNF
Cattaneo–Christov heat flux modelCCHM
Gyrotactic microorganismGM
Aluminum oxideAl2O3
CopperCu
Carbon nanotubeCNTs
Boundary layer regionBLR

References

  1. Choi, S.U.S. Enhancing Thermal Conductivity of Fluids with Nanoparticles. Am. Soc. Mech. Eng. Fluids Eng. Div. FED 1995, 231, 99–105. [Google Scholar]
  2. Sahoo, R.R. Thermo-Hydraulic Characteristics of Radiator with Various Shape Nanoparticle-Based Ternary Hybrid Nanofluid. Powder Technol. 2020, 370, 19–28. [Google Scholar] [CrossRef]
  3. Adun, H.; Kavaz, D.; Dagbasi, M. Review of Ternary Hybrid Nanofluid: Synthesis, Stability, Thermophysical Properties, Heat Transfer Applications, and Environmental Effects. J. Clean. Prod. 2021, 328, 129525. [Google Scholar] [CrossRef]
  4. Xuan, Z.; Zhai, Y.; Ma, M.; Li, Y.; Wang, H. Thermo-Economic Performance and Sensitivity Analysis of Ternary Hybrid Nanofluids. J. Mol. Liq. 2021, 323, 114889. [Google Scholar] [CrossRef]
  5. Sahoo, R.R.; Kumar, V. Development of a New Correlation to Determine the Viscosity of Ternary Hybrid Nanofluid. Int. Commun. Heat Mass Transf. 2020, 111, 104451. [Google Scholar] [CrossRef]
  6. Sahoo, R.R. Experimental Study on the Viscosity of Hybrid Nanofluid and Development of a New Correlation. Heat Mass Transf. 2020, 56, 3023–3033. [Google Scholar] [CrossRef]
  7. Animasaun, I.L.; Yook, S.J.; Muhammad, T.; Mathew, A. Dynamics of Ternary-Hybrid Nanofluid Subject to Magnetic Flux Density and Heat Source or Sink on a Convectively Heated Surface. Surf. Interfaces 2022, 28, 101654. [Google Scholar] [CrossRef]
  8. Alanazi, M.M.; Ahmed Hendi, A.; Ahammad, N.A.; Ali, B.; Majeed, S.; Shah, N.A. Significance of Ternary Hybrid Nanoparticles on the Dynamics of Nanofluids over a Stretched Surface Subject to Gravity Modulation. Mathematics 2023, 11, 809. [Google Scholar] [CrossRef]
  9. Raju, C.S.K.; Ameer, N.A.; Kiran, S.; Nehad, A.S.; Se-jin, Y.; Dinesh, M.K. Nonlinear linear movements of axisymmetric ternary hybrid nanofluids in a thermally radiated expanding or contracting permeable Darcy Walls with different shapes and densities: Simple linear regression. Int. Commun. Heat Mass Transfer 2022, 135, 106110. [Google Scholar] [CrossRef]
  10. Ramzan, M.; Dawar, A.; Saeed, A.; Kumam, P.; Sitthithakerngkiet, K.; Lone, S.A. Analysis of the Partially Ionized Kerosene Oil-Based Ternary Nanofluid Flow over a Convectively Heated Rotating Surface. Open Phys. 2022, 20, 507–525. [Google Scholar] [CrossRef]
  11. Fourier, J. Théorie Analytique de La Chaleur; Chez Firmin Didot: Paris, France, 1822. [Google Scholar]
  12. Cattaneo, C. Sulla Conduzione Del Calore. Atti Sem. Mat. Fis. Univ. Modena 1948, 3, 83–101. [Google Scholar]
  13. 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]
  14. Venkateswarlu, S.; Varma, S.V.K.; Durga Prasad, P. MHD Flow of MoS2 and MgO Water-Based Nanofluid through Porous Medium over a Stretching Surface with Cattaneo-Christov Heat Flux Model and Convective Boundary Condition. Int. J. Ambient Energy 2020, 43, 2940–2949. [Google Scholar] [CrossRef]
  15. Rawat, S.K.; Kumar, M. Cattaneo–Christov Heat Flux Model in Flow of Copper Water Nanofluid Through a Stretching/Shrinking Sheet on Stagnation Point in Presence of Heat Generation/Absorption and Activation Energy. Int. J. Appl. Comput. Math. 2020, 6, 112. [Google Scholar] [CrossRef]
  16. Ramzan, M.; Gul, H.; Baleanu, D.; Nisar, K.S.; Malik, M.Y. Role of Cattaneo–Christov Heat Flux in an MHD Micropolar Dusty Nanofluid Flow with Zero Mass Flux Condition. Sci. Rep. 2021, 11, 19528. [Google Scholar] [CrossRef]
  17. Lv, Y.P.; Gul, H.; Ramzan, M.; Chung, J.D.; Bilal, M. Bioconvective Reiner–Rivlin Nanofluid Flow over a Rotating Disk with Cattaneo–Christov Flow Heat Flux and Entropy Generation Analysis. Sci. Rep. 2021, 11, 15859. [Google Scholar] [CrossRef]
  18. Abderrahmane, A.; Qasem, N.A.A.; Younis, O.; Marzouki, R.; Mourad, A.; Shah, N.A.; Chung, J.D. MHD Hybrid Nanofluid Mixed Convection Heat Transfer and Entropy Generation in a 3-D Triangular Porous Cavity with Zigzag Wall and Rotating Cylinder. Mathematics 2022, 10, 769. [Google Scholar] [CrossRef]
  19. Kumar, M.D.; Raju, C.S.K.; Sajjan, K.; El-Zahar, E.R.; Shah, N.A. Linear and quadratic convection on 3D flow with transpiration and hybrid nanoparticles. Int. Commun. Heat Mass Transfer 2022, 134, 105995. [Google Scholar] [CrossRef]
  20. Sajjan, K.; Shah, N.A.; Ahammad, N.A.; Raju, C.S.K.; Kumar, M.D.; Weera, W. Nonlinear Boussinesq and Rosseland approximations on 3D flow in an interruption of Ternary nanoparticles with various shapes of densities and conductivity properties. AIMS Math. 2022, 7, 18416–18449. [Google Scholar] [CrossRef]
  21. Kuznetsov, A.V. Bio-Thermal Convection Induced by Two Different Species of Microorganisms. Int. Commun. Heat Mass Transf. 2011, 38, 548–553. [Google Scholar] [CrossRef]
  22. Shi, Q.-H.H.; Hamid, A.; Khan, M.I.; Kumar, R.N.; Gowda, R.J.P.P.; Prasannakumara, B.C.; Shah, N.A.; Khan, S.U.; Chung, J.D.; Naveen Kumar, R.; et al. Numerical Study of Bio-Convection Flow of Magneto-Cross Nanofluid Containing Gyrotactic Microorganisms with Activation Energy. Sci. Rep. 2021, 11, 16030. [Google Scholar] [CrossRef] [PubMed]
  23. Al-Khaled, K.; Ullah Khan, S. Thermal Aspects of Casson Nanoliquid with Gyrotactic Microorganisms, Temperature-Dependent Viscosity, and Variable Thermal Conductivity: Bio-Technology and Thermal Applications. Inventions 2020, 5, 39. [Google Scholar] [CrossRef]
  24. Bhatti, M.M.; Marin, M.; Zeeshan, A.; Ellahi, R.; Abdelsalam, S.I. Swimming of Motile Gyrotactic Microorganisms and Nanoparticles in Blood Flow Through Anisotropically Tapered Arteries. Front. Phys. 2020, 8, 95. [Google Scholar] [CrossRef] [Green Version]
  25. Kairi, R.R.; Shaw, S.; Roy, S.; Raut, S. Thermosolutal Marangoni Impact on Bioconvection in Suspension of Gyrotactic Microorganisms over an Inclined Stretching Sheet. J. Heat Transf. 2021, 143, 031201. [Google Scholar] [CrossRef]
  26. Ali, A.; Sarkar, S.; Das, S.; Jana, R.N. Investigation of Cattaneo–Christov Double Diffusions Theory in Bioconvective Slip Flow of Radiated Magneto-Cross-Nanomaterial Over Stretching Cylinder/Plate with Activation Energy. Int. J. Appl. Comput. Math. 2021, 7, 208. [Google Scholar] [CrossRef]
  27. Mishra, A.; Kumar, M. Numerical Analysis of MHD Nanofluid Flow over a Wedge, Including Effects of Viscous Dissipation and Heat Generation/Absorption, Using Buongiorno Model. Heat Transf. 2021, 50, 8453–8474. [Google Scholar] [CrossRef]
  28. Yaseen, M.; Rawat, S.K.; Kumar, M. Falkner–Skan Problem for a Stretching or Shrinking Wedge With Nanoparticle Aggregation. J. Heat Transf. 2022, 144, 102501. [Google Scholar] [CrossRef]
  29. Mahanthesh, B.; Mackolil, J. Flow of Nanoliquid Past a Vertical Plate with Novel Quadratic Thermal Radiation and Quadratic Boussinesq Approximation: Sensitivity Analysis. Int. Commun. Heat Mass Transf. 2021, 120, 105040. [Google Scholar] [CrossRef]
  30. Gumber, P.; Yaseen, M.; Kumar, S.; Kumar, M. Heat Transfer in Micropolar Hybrid Nanofluid Flow Past a Vertical Plate in the Presence of Thermal Radiation and Suction/Injection Effects. Partial Differ. Equ. Appl. Math. 2022, 5, 100240. [Google Scholar] [CrossRef]
  31. Rawat, S.K.; Upreti, H.; Kumar, M. Comparative Study of Mixed Convective MHD Cu-Water Nanofluid Flow over a Cone and Wedge Using Modified Buongiorno’s Model in Presence of Thermal Radiation and Chemical Reaction via Cattaneo-Christov Double Diffusion Model. J. Appl. Comput. Mech. 2021, 7, 1383–1402. [Google Scholar] [CrossRef]
  32. Sandeep, N.; Reddy, M.G. Heat Transfer of Nonlinear Radiative Magnetohydrodynamic Cu-Water Nanofluid Flow over Two Different Geometries. J. Mol. Liq. 2017, 225, 87–94. [Google Scholar] [CrossRef]
  33. Reddy, M.G.; Rani, M.V.V.N.L.S.; Kumar, K.G.; Prasannakumara, B.C. Cattaneo–Christov Heat Flux and Non-Uniform Heat-Source/Sink Impacts on Radiative Oldroyd-B Two-Phase Flow across a Cone/Wedge. J. Braz. Soc. Mech. Sci. Eng. 2018, 40, 95. [Google Scholar] [CrossRef]
  34. Jayachandra Babu, M.; Sandeep, N.; Saleem, S. Free Convective MHD Cattaneo-Christov Flow over Three Different Geometries with Thermophoresis and Brownian Motion. Alex. Eng. J. 2017, 56, 659–669. [Google Scholar] [CrossRef]
  35. Garia, R.; Rawat, S.K.; Kumar, M.; Yaseen, M. Hybrid Nanofluid Flow over Two Different Geometries with Cattaneo–Christov Heat Flux Model and Heat Generation: A Model with Correlation Coefficient and Probable Error. Chin. J. Phys. 2021, 74, 421–439. [Google Scholar] [CrossRef]
  36. Arif, M.; Kumam, P.; Kumam, W.; Mostafa, Z. Heat Transfer Analysis of Radiator Using Different Shaped Nanoparticles Water-Based Ternary Hybrid Nanofluid with Applications: A Fractional Model. Case Stud. Therm. Eng. 2022, 31, 101837. [Google Scholar] [CrossRef]
  37. Elnaqeeb, T.; Animasaun, I.L.; Shah, N.A. Ternary-Hybrid Nanofluids: Significance of Suction and Dual-Stretching on Three-Dimensional Flow of Water Conveying Nanoparticles with Various Shapes and Densities. Z. Naturforsch.—Sect. A J. Phys. Sci. 2021, 76, 231–243. [Google Scholar] [CrossRef]
  38. Vajravelu, K.; Nayfeh, J. Hydromagnetic Convection at a Cone and a Wedge. Int. Commun. Heat Mass Transf. 1992, 19, 701–710. [Google Scholar] [CrossRef]
  39. Yaseen, M.; Rawat, S.K.; Kumar, M. Cattaneo–Christov Heat Flux Model in Darcy–Forchheimer Radiative Flow of MoS2–SiO2/Kerosene Oil between Two Parallel Rotating Disks. J. Therm. Anal. Calorim. 2022, 147, 10865–10887. [Google Scholar] [CrossRef]
  40. Shampine, L.F.; Gladwell, I.; Thompson, S. Solving ODEs with MATLAB; Cambridge University Press: Cambridge, UK, 2003; ISBN 9780511615542. [Google Scholar]
Figure 1. Physical model.
Figure 1. Physical model.
Mathematics 11 01237 g001
Figure 2. Flow chart.
Figure 2. Flow chart.
Mathematics 11 01237 g002
Figure 3. Effect of λ on f ( η ) .
Figure 3. Effect of λ on f ( η ) .
Mathematics 11 01237 g003
Figure 4. Effect of φ 1 on f ( η ) .
Figure 4. Effect of φ 1 on f ( η ) .
Mathematics 11 01237 g004
Figure 5. Effect of φ 2 on f ( η ) .
Figure 5. Effect of φ 2 on f ( η ) .
Mathematics 11 01237 g005
Figure 6. Effect of φ 3 on f ( η ) .
Figure 6. Effect of φ 3 on f ( η ) .
Mathematics 11 01237 g006
Figure 7. Effect of Su on f ( η ) .
Figure 7. Effect of Su on f ( η ) .
Mathematics 11 01237 g007
Figure 8. Velocity boundary layer patterns for (a) flat plate, (b) wedge, and (c) cone.
Figure 8. Velocity boundary layer patterns for (a) flat plate, (b) wedge, and (c) cone.
Mathematics 11 01237 g008aMathematics 11 01237 g008b
Figure 9. Streamlines for (a) flat plate, (b) wedge, and (c) cone.
Figure 9. Streamlines for (a) flat plate, (b) wedge, and (c) cone.
Mathematics 11 01237 g009
Figure 10. Effect of Q on θ ( η ) .
Figure 10. Effect of Q on θ ( η ) .
Mathematics 11 01237 g010
Figure 11. Effect of φ 1 on θ ( η ) .
Figure 11. Effect of φ 1 on θ ( η ) .
Mathematics 11 01237 g011
Figure 12. Effect of φ 2 on θ ( η ) .
Figure 12. Effect of φ 2 on θ ( η ) .
Mathematics 11 01237 g012
Figure 13. Effect of φ 3 on θ ( η ) .
Figure 13. Effect of φ 3 on θ ( η ) .
Mathematics 11 01237 g013
Figure 14. Effect of Rd on θ ( η ) .
Figure 14. Effect of Rd on θ ( η ) .
Mathematics 11 01237 g014
Figure 15. Effect of Su on θ ( η ) .
Figure 15. Effect of Su on θ ( η ) .
Mathematics 11 01237 g015
Figure 16. Effect of γ t on θ ( η ) .
Figure 16. Effect of γ t on θ ( η ) .
Mathematics 11 01237 g016
Figure 17. Effect of λ on θ ( η ) .
Figure 17. Effect of λ on θ ( η ) .
Mathematics 11 01237 g017
Figure 18. Effect of Le on ϕ ( η ) .
Figure 18. Effect of Le on ϕ ( η ) .
Mathematics 11 01237 g018
Figure 19. Effect of Su on ϕ ( η ) .
Figure 19. Effect of Su on ϕ ( η ) .
Mathematics 11 01237 g019
Figure 20. Effect of γ c on ϕ ( η ) .
Figure 20. Effect of γ c on ϕ ( η ) .
Mathematics 11 01237 g020
Figure 21. Effect of Lb on ζ ( η ) .
Figure 21. Effect of Lb on ζ ( η ) .
Mathematics 11 01237 g021
Figure 22. Effect of Pe on ζ ( η ) .
Figure 22. Effect of Pe on ζ ( η ) .
Mathematics 11 01237 g022
Figure 23. Effect of Su on ζ ( η ) .
Figure 23. Effect of Su on ζ ( η ) .
Mathematics 11 01237 g023
Figure 24. Effect of ϖ on ζ ( η ) .
Figure 24. Effect of ϖ on ζ ( η ) .
Mathematics 11 01237 g024
Figure 25. Nusselt number (a) flat plate, (b) wedge, and (c) cone.
Figure 25. Nusselt number (a) flat plate, (b) wedge, and (c) cone.
Mathematics 11 01237 g025aMathematics 11 01237 g025b
Table 1. Thermo-physical properties of water, Al2O3, Cu, and CNT nanoparticles [15,36].
Table 1. Thermo-physical properties of water, Al2O3, Cu, and CNT nanoparticles [15,36].
ρ (kg/m3)Cp (J/kgK)k (W/mK)β (K−1) ShapeSphericity
Water997.141790.61321
Al2O33970765400.85Sphericalψ = 1
Cu89333854011.67Plateletψ = 0.612
CNT260042566001.6 × 10−6Cylindricalψ = 0.52
Table 2. Thermophysical properties of ternary hybrid nanofluid [36,37].
Table 2. Thermophysical properties of ternary hybrid nanofluid [36,37].
Density ρ t h n f = ( 1 ϕ 1 ϕ 2 ϕ 3 ) ρ b f + ϕ 1 ρ s p 1 + ϕ 2 ρ s p 2 + ϕ 3 ρ s p 3
Heat capacitance ( ρ C p ) t h n f = ( 1 ϕ 1 ϕ 2 ϕ 3 ) ( ρ C p ) b f + ϕ 1 ( ρ C p ) s p 1 + ϕ 2 ( ρ C p ) s p 2 + ϕ 3 ( ρ C p ) s p 3 ,
Modified Maxwell model: k n f i k b f = k i + ( n 1 ) k b f + ( n 1 ) ϕ i ( k i k b f ) k i + ( n 1 ) k b f ϕ i ( k i k b f )
where, n = ( 3 ψ ) is shape factor.
Nanoparticle—1(Spherical) { μ n f 1 μ b f = 1 + 2.5 ϕ 1 + 6.2 ϕ 1 2 k n f 1 k b f = k 1 + 2 k b f + 2 ϕ 1 ( k 1 k b f ) k 1 + 2 k b f ϕ 1 ( k 1 k b f )
Nanoparticle—2(Platelet) { μ n f 2 μ b f = 1 + 13.5 ϕ 2 + 904.4 ϕ 2 2 k n f 2 k b f = k 2 + 3.9 k b f + 3.9 ϕ 2 ( k 2 k b f ) k 2 + 3.9 k b f ϕ 2 ( k 2 k b f )
Nanoparticle—3(Cylindrical) { μ n f 3 μ b f = 1 + 37.1 ϕ 3 + 612.6 ϕ 3 2 k n f 3 k b f = k 3 + 4.7 k b f + 4.7 ϕ 3 ( k 3 k b f ) k 3 + 4.7 k b f ϕ 3 ( k 3 k b f )
Viscosity μ t h n f = μ n f 1 ϕ 1 + μ n f 2 ϕ 2 + μ n f 3 ϕ 3 ϕ t h n f
Thermal conductivity k t h n f = k n f 1 ϕ 1 + k n f 2 ϕ 2 + k n f 3 ϕ 3 ϕ t h n f Where ϕ t h n f = ϕ 1 + ϕ 2 + ϕ 3
Table 3. Comparison of values of f ( 0 )   and θ ( 0 ) with results of Vajravelu and Nayfeh [38] for the flow over cone when Rd = 0, γ t = γ c = 0 , Su = 0; Le = Lb = Pe = ϖ = 0.
Table 3. Comparison of values of f ( 0 )   and θ ( 0 ) with results of Vajravelu and Nayfeh [38] for the flow over cone when Rd = 0, γ t = γ c = 0 , Su = 0; Le = Lb = Pe = ϖ = 0.
QPrGrMs
f ( 0 )   ( See   [ 38 ] )
f ( 0 )   [ Present ]
θ ( 0 )   ( See   [ 38 ] )
θ ( 0 )   [ Present ]
−50.3−0.51−2.1−0.155592−0.15570252−2.237475−2.23538986
−50.3−0.512.1−0.156001−0.15588966−2.232780−2.23418090
−50.3−0.532.1−0.126400−0.12634921−2.233732−2.23472524
−51.00.532.10.1252600.12541914−2.245321−2.24047000
Table 4. Numerical values of the Nusselt number (Nu), Sherwood number (Sh), and motile microorganisms’ density gradient (Nm).
Table 4. Numerical values of the Nusselt number (Nu), Sherwood number (Sh), and motile microorganisms’ density gradient (Nm).
Cone ( n = 1 ,   α = π 4 ) Wedge ( n = 0 ,   α = π 4 ) Plate ( n = 0 ,   α = 0 )
Pe ϖ γ t γ C LeLb λ QRdSNuShNmNuShNmNuShNm
0.10.21.50.011.50.50.2−2100.124.92767−1.51382−0.8794723.00199−1.09132−0.633223.02106−1.09284−0.63419
0.2 --−1.03107--−0.74255--−0.74368
0.4 --−1.33868--−0.96453--−0.96592
0.10.7 --−0.93948--−0.67622--−0.67727
1.7 --−1.0595--−0.76227--−0.76342
0.20.9 20.81458−1.51501−0.8805619.29351−1.09256−0.6341419.31295−1.09457−0.63551
1.3 23.60737−1.51414−0.8797621.80354−1.09167−0.6334621.82288−1.09332−0.63456
1.50.06 -−1.38953−0.86634-−1.04645−0.62838-−1.04774−0.62935
0.16 -−1.16796−0.84311-−0.96177−0.6193-−0.96266−0.62024
0.012 -−1.81643−0.91215-−1.3191−0.65807-−1.32073−0.65908
3 -−2.3356−0.96938-−1.71464−0.70207-−1.71641−0.70309
1.50.6 --−0.97675--−0.69807--−0.69919
0.7 --−1.06895--−0.76131--−0.76253
0.5−1.8 24.55437−1.47761−0.8522722.49837−1.04949−0.6057522.28255−1.03021−0.59317
2.2 25.26119−1.5448−0.9025923.4335−1.12445−0.6550323.61331−1.13766−0.66375
0.2−1.4 21.88554−1.51457−0.8801320.27092−1.09219−0.6338520.294−1.09404−0.6351
−0.8 18.43186−1.51581−0.8813117.15315−1.09365−0.6350217.18372−1.09606−0.63672
0.4 7.816058−1.53274−0.902687.173743−1.11078−0.651757.408173−1.11788−0.65803
−212 26.66715−1.51439−0.8799824.58693−1.09197−0.6336824.61154−1.09374−0.63487
14 28.25102−1.51495−0.8804926.03504−1.0926−0.6341726.06564−1.09462−0.63555
16 29.70916−1.51549−0.8810127.37275−1.09321−0.6346527.40974−1.09547−0.63623
10−0.122.27847−1.287−0.7790519.99349−0.85639−0.5317620.0157−0.85824−0.53294
0.124.92767−1.51382−0.8794723.00199−1.09132−0.633223.02106−1.09284−0.63419
0.327.94773−1.7543−0.9859226.47821−1.3467−0.7431526.49433−1.3479−0.74397
Table 5. Numerical values of the Nusselt number (Nu), Sherwood number (Sh), and motile microorganisms’ density gradient (Nm) for volume fraction of nanoparticles.
Table 5. Numerical values of the Nusselt number (Nu), Sherwood number (Sh), and motile microorganisms’ density gradient (Nm) for volume fraction of nanoparticles.
. Cone ( n = 1 ,   α = π 4 ) Wedge ( n = 0 ,   α = π 4 ) Plate ( n = 0 ,   α = 0 )
φ 1 φ 2 φ 3 NuShNmNuShNmNuShNm
0.020.020.0224.92767−1.51382−0.8794723.00199−1.09132−0.633223.02106−1.09284−0.63419
0.03 24.84579−1.41239−0.817122.92604−1.01785−0.5891822.94572−1.01933−0.59014
0.04 24.78749−1.33885−0.7720322.87163−0.96459−0.5573722.89167−0.96603−0.5583
0.020.03 25.06144−1.67807−0.9809323.12018−1.21024−0.704723.13655−1.21166−0.70564
0.04 25.26947−1.98846−1.1740423.30667−1.43523−0.8407923.31994−1.43655−0.84169
0.020.0325.15922−1.86527−1.0972423.20391−1.346−0.7867223.21942−1.34746−0.78771
0.0425.4086−2.36089−1.4076223.42386−1.70561−1.0055123.43616−1.70703−1.00649
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

Yaseen, M.; Rawat, S.K.; Shah, N.A.; Kumar, M.; Eldin, S.M. Ternary Hybrid Nanofluid Flow Containing Gyrotactic Microorganisms over Three Different Geometries with Cattaneo–Christov Model. Mathematics 2023, 11, 1237. https://doi.org/10.3390/math11051237

AMA Style

Yaseen M, Rawat SK, Shah NA, Kumar M, Eldin SM. Ternary Hybrid Nanofluid Flow Containing Gyrotactic Microorganisms over Three Different Geometries with Cattaneo–Christov Model. Mathematics. 2023; 11(5):1237. https://doi.org/10.3390/math11051237

Chicago/Turabian Style

Yaseen, Moh, Sawan Kumar Rawat, Nehad Ali Shah, Manoj Kumar, and Sayed M. Eldin. 2023. "Ternary Hybrid Nanofluid Flow Containing Gyrotactic Microorganisms over Three Different Geometries with Cattaneo–Christov Model" Mathematics 11, no. 5: 1237. https://doi.org/10.3390/math11051237

APA Style

Yaseen, M., Rawat, S. K., Shah, N. A., Kumar, M., & Eldin, S. M. (2023). Ternary Hybrid Nanofluid Flow Containing Gyrotactic Microorganisms over Three Different Geometries with Cattaneo–Christov Model. Mathematics, 11(5), 1237. https://doi.org/10.3390/math11051237

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