Next Article in Journal
PLGA Nanoparticles Decorated with Anti-HER2 Affibody for Targeted Delivery and Photoinduced Cell Death
Next Article in Special Issue
Lattice Model of Multilayer Adsorption of Particles with Orientation Dependent Interactions at Solid Surfaces
Previous Article in Journal
SARs for the Antiparasitic Plant Metabolite Pulchrol. 3. Combinations of New Substituents in A/B-Rings and A/C-Rings
Previous Article in Special Issue
Fractional Dynamics of HIV with Source Term for the Supply of New CD4+ T-Cells Depending on the Viral Load via Caputo–Fabrizio Derivative
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bioconvection Due to Gyrotactic Microorganisms in Couple Stress Hybrid Nanofluid Laminar Mixed Convection Incompressible Flow with Magnetic Nanoparticles and Chemical Reaction as Carrier for Targeted Drug Delivery through Porous Stretching Sheet

1
Deanship of Combined First Year, Umm Al-Qura University Makkah, Mecca P.O. Box 715, Saudi Arabia
2
Department of Mathematics, Abdul Wali Khan University, Mardan 23200, Pakistan
3
Department of Mathematics, University of Swabi, Swabi 94640, Pakistan
*
Author to whom correspondence should be addressed.
Molecules 2021, 26(13), 3954; https://doi.org/10.3390/molecules26133954
Submission received: 18 April 2021 / Revised: 11 June 2021 / Accepted: 23 June 2021 / Published: 28 June 2021

Abstract

:
In this paper, the steady electrically conducting hybrid nanofluid (CuO-Cu/blood) laminar-mixed convection incompressible flow at the stagnation-point with viscous and gyrotactic microorganisms is considered. Additionally, hybrid nanofluid flow over a horizontal porous stretching sheet along with an induced magnetic field and external magnetic field effects that can be used in biomedical fields, such as in drug delivery and the flow dynamics of the microcirculatory system. This investigation can also deliver a perfect view about the mass and heat transfer behavior of blood flow in a circulatory system and various hyperthermia treatments such as the treatment of cancer. The simple partial differential equations (PDEs) are converted into a series of dimensional ordinary differential equations (ODEs), which are determined using appropriate similarities variables (HAM). The influence of the suction or injection parameter, mixed convection, Prandtl number, buoyancy ratio parameter, permeability parameter, magnetic parameter, reciprocal magnetic prandtl number, bioconvection Rayleigh number, coupled stress parameter, thermophoretic parameter, Schmidt number, inertial parameter, heat source parameter, and Brownian motion parameter on the concentration, motile microorganisms, velocity, and temperature is outlined, and we study the physical importance of the present problem graphically.

1. Introduction

Bloodstream investigation in a human circulatory framework has developed amazing revenue in biotechnology and the world of medicine since most human diseases were caused by unsatisfactory supplies of blood to the lungs, veins, corridors, tissues, and systole stages. Numerous circulatory framework problems, including atherosclerosis embolism, aspiratory embolism and a blockage of blood supply in veins, nerves, and corridors, instigate a heart attack, stroke and ischemic thoracic inconvenience. By increased the blood temperature of 39 °C through 42 °C (hyperthermia) from average body temperature, the heart output doubles and the circulation increases. Nanomaterials have recently received a lot of attention in the field of biomedicine because of their useful applications, such as anticancer drug delivery, biosensing, antibacteria, and cell imaging, etc. Magnetic nanoparticles are extremely useful in magnetic drug targeting and magnetic resonance imaging agents, among other applications. Numerous studies have been conducted to determine the significance of nanoparticles in biological sciences. Chauhan and Tiwari [1] worked on non-Newtonian Herschel–Bulkley liquid and analyzed the heat transfer on blood movement in previous veins. They saw higher accuracy paper estimates generally decrease the blood velocity in veins. It is used for many positive therapies, such as cancers, cardiac drugs, and malignant development (Deussen [2] and Deniz [3]). With the aid of peristaltic and nanofluid, Bég and Tripathi [4] proposed the reconstruction of Mathematica’s bioengineering concept. Kothandapani and Prakash [5] found a heating source on an asymmetrical pointing channel on a non-Newtonian excessive digression nanofluid model. Akbar [6] examined delayed blood propagation of metal-based nanomaterial through the shaped stenotic route and explained nanomedicine applications. In Bhatti’s study [7], the properties and implementations of the vector viscosity blood clot model were investigated. The two-step model of peristalsis was considered in Dinarvand [8]. The constant laminar-blended mixture viscous and hybrid (CUO-Cu/blood) hybrid fluid flows near the plane stagnations on a level, permeable, linearly stretched board with an adjustable magnetic flux through a new nanoparticle and based liquid measurement [9]. Majee et al. conducted a methodical report on shaky blood progression with magnetic nanoparticles and is supposed to carry out the design of streams and nanoparticles in an infected blood vessel section that has atherosclerosis. Varshney [10] mathematically researched the pulsatile movement of blood going through a tightening vein, while the speed increase in the body is rambling. Jinga’s [11] research was conducted by combining a hybrid discrete component and a unit monitoring approach to calculate and stress the transmission behaviors of fractured crystalline rocks. The inspiration is the significance of insightful pressure impacts on the behavior, which are critical for the estimation of the environmental protection of many rock building projects, of the impurity transportation of broken crystalline rocks. Noorishad [12] presents a new means for the accurate testing of liquid stream lead in cracked permeable media, which is presented here. To do this, mechanical and liquid stream limits of both permeable and breakage media are used as part of an increase in Blot’s three-dimensional union hypothesis. Ellahi [13] introduced the peristaltic fluid stream between two coaxial cylinders of different forms and designs. The nanofluid consists of gold particles, while the pair of pressure fluids are filled as solvents. Choi and Eastman [14] created the term nanofluid, and this fluid is generated through a dilute inspection of solid particulate matter of 1–100 nm in constant fluids (oil, water, etc.). By the inclusion of ZnO, Cu, SiO2, TiO2, and Al2O3 nanopowders, the efficiency of the heat transfer of regular fluids has been greatly increased. In recent years, several investigators have discovered hypothetically and experimentally the characteristics of heat transfer from different nanoparticles in many manufacturing processes, development processes, and the application of renewable energy [15,16,17,18,19,20,21]. Researchers have developed various models for studying the Tiwari and Das model of nanofluids. Late in life, various researchers hypothetically and provisionally identified warmth motion attributes for several mechanical cycles, manufacturing, and environmentally friendly energy applications of different nanoparticles [15,16,17,18,19,20,21]. Specialists in nanofluids in which the Tiwari and Das model was presented with different models.
The macroscopic movement of fluid induces additional flexibility in swimming microorganisms, known as bioconvection, as a consequence of the 3-D variant in density over one region. The self-driving mobile microorganisms aim to boost the base fluid, creating a bioconvective stream in a specific direction. The travelling microorganisms are classified into different categories of chemical or oxytactical, gyrotactic, and negative gravitational characteristics. Nanoparticles are not self-regulated in comparison to mobile microorganisms, and the influence of the Brownian motion and the effect of thermophoresis is responsible for their motion. Nanofluid bioconvection is supposed to be feasible if the convergence of nanoparticles is low, and then, the choice for improved fluid thickness in the base is not sufficient. Basha [22] provided a mathematical response to the blood nanofluid fluid quality of a vehicle streaming across the plate, wedge, and stagnating stage. The effects of non-linear radiation, sticky dispersion, convinced magnetic field, and material reaction and the pertinence of the microorganism’s properties are evaluated. Bhatti [23] also studied the behavior of a changeable magnetic field and blood clot model using Jeffrey fluid nanoparticles and medication models. Ahmed [24] has regarded the magnetized, non-Darcy, permeable, laminar circulation of nanofluid and gyrotactic microorganisms. Kuznetsov [25] submitted nanoparticular suspension using the principle of Buongiorno for gyrotactic microorganisms. Raju [26] issued a mathematical study to determine Casson nanofluid stream liquid vehicles using gyrotactic microorganisms, which showed that adding the wedge point limits reduced the microorganism’s thickness. Raju and Rashad [27] addressed gyrotactic microorganisms’ impact on a chemical nanofluid flow over a vertical cylinder and noted that the amount of Rayleigh bioconvection increases the density of microorganisms. Liao (1992) [28] has observed that this technique is fast convergent to the approximate solution, and it is the best fit for the solution of non-linear problems. We considered an electrically guided coupled pressure crossover (CuO-Cu/blood) nanofluid stream comprised of gyrotactic microorganisms pushed near the plane stagnation-point over a level, permeable, extending layer alongside an outer attractive field and prompted attractive field impacts in the current study. This investigation will also provide a good image of the temperature and mass exchange activity of blood in a circulatory system, as well as various hyperthermia treatments, such as cancer care. The Buongiorno model was used to demonstrate thermophoresis and Brownian dispersion. The sheet is permeable, and the surface of the stretching sheet has an injection effect. Our model is mathematically formulated by deriving the governing equations and applying sufficient similarity transformations. Using the Homotopy Analysis process, we can solve our modeled problem.

1.1. Problem Mathematical Modeling

We are assuming the continuously laminar-varied convection compact viscous and the electrically guided coupling pressure Darcy–Forchheimer CuO-Cu/Blood hybrid nanofluid fluids and heat close to the stagnation point on a smooth, directly extending the plate below an exterior magnetic flux, as defined in Figure 1. The thermophoresis and Brownian diffusion effects were analyzed using the Buongiorno model. The sheet is permeable, and the stretching sheet’s surface has an injection effect. Moreover, the suction/injection velocity of the sheet is V 0 , while the linear stretching velocity is u w = c x . Furthermore, the sheet’s temperature depends on the stagnation point, i.e., T w = T + T 0 x l   . The control of the boundary layer and external magnetic flux towards x- direction can be represented by u e = a x and H e ( x ) = H 0 x l . We stress that this evaluation is a characteristic length l of the uniform magnetic field H 0 in the upstream infinity. Table 1 display the thermophysical properties of cp, ρ, k and β. The following principles can be expressed in the following terms: simple non-linear PDEs.
u x       +       v y = 0 ,
H 1 x       +       H 2 y = 0 ,
u u x       +       v u y μ e 4 π ρ h n f ( H 2 H 1 x       +     H 2   H 1 y ) = μ e d μ e d x   μ e 4 π ρ h n f H e d H e d x + μ h n f ρ h n f 2 u y 2 μ h n f ρ h n f u k * η 0 ρ h n f 4 u y 4 + ( ρ m ρ f ) ( n n ) ρ h n f g ( ρ p ρ f ) ( C C ) ρ h n f g + ( 1 C ) β ρ f ( T T ) g ρ h n f ,
u H 1 x       +     v H 1 y H 1 u x H 2 u y = η 0 2 H 1 y 2 ,
u T x       +       v T y = k h n f ( ρ c p ) h n f 2 T y 2 + τ [ D B ( C y T y ) + D T Τ ( T y ) 2 ] + Q ( ρ c p ) h n f ( T T ) ,
u C x       +       v C y = D B 2 C y 2 + ( D T Τ ) 2 T y 2 K r ( C C ) ,
u N x + v N y + b * W c ( C w C ) ( N C y ) y = D N N y y ,
Following to the boundary conditions:
u = u w = c x , v = V 0 , H 1 y = H 2 = 0 , T = T w = T + T 0 x l   ,   D B C y + D T Τ T y , N = N w     a t   y = 0 , u = μ e = a x , H 1 = H e = H 0 x l , T = T , C = C ,   N = N                                                                                           a t     y = .
where the position temperature T 0 , the atmospheric temperature Τ , magnetic permeability μ e , u and v , alongside x and y axes H 1 and H 2 are the velocity units with persuaded components of the magnetic field on each. The dimension k * is the magnetic absorption potential of the permeable intermediate. Diffusivity g is gravity acceleration, temperature is T , β is the thermal growth volumetric constant, and Q is the heat generation volumetric rate/absorption. ρ h n f , μ h n f , ( ρ c p ) h n f and k h n f are the density, viscosity, heat power volumetric, and hybrid nanofluid thermal conductivity calculated, respectively. These can be seen from Table 2.
As we note, the experimental type factor of the nanoparticles seen in Figure 2 is the classic Hamilton–Crosser approximate for real thermal conductivity.
It is worthwhile to mention here that, we suggest ϕ , ρ s and ( c p ) s as the equivalent nanoparticle volume fraction, the equivalent density of nanoparticles and the equivalents pecific heat at constant pressure of nanoparticles, respectively. Moreover, respectively, ϕ 1 and ϕ 2 are the 1st and 2nd nanoparticles’ volume fraction of these compatible formulae:
ρ s = ( ρ 1 × w 1 ) + ( ρ 2 × w 2 ) w 1 + w 2
( ρ c p ) s = ( ( c p ) 1 × w 1 ) + ( ( c p ) 2 × w 2 ) w 1 + w 2
Table 2. Adapted frameworks and thermophysical properties for the hybrid nanofluid.
Table 2. Adapted frameworks and thermophysical properties for the hybrid nanofluid.
Property Hybrid Nanofluid
Viscosity (μ) μ f ( 1 ϕ ) 2.5
Volumetric heat capacity (ρ) ( 1 ϕ ) ρ f + ϕ ρ s
Volumetric heat capacity (ρcp) ( 1 ϕ ) ( ρ c p ) f + ϕ ( ρ c p ) s
Thermal conductivity (k) k 2 + ( n 2 1 ) k f ( n 2 1 ) ϕ 1 ( k f k 2 ) k 2 + ( n 2 1 ) k n f + ϕ 2 ( k n f k 2 )
× k 1 + ( n 1 1 ) k f ( n 1 1 ) ϕ 1 ( k f k 1 ) k 1 + ( n 1 1 ) k n f + ϕ 2 ( k n f k 1 ) × k f ;
k n f = k 1 + ( n 1 1 ) k f ( n 1 1 ) ϕ 1 ( k f k 1 ) k 1 + ( n 1 1 ) k n f + ϕ 2 ( k n f k 1 ) × k f
ϕ 1 = w 1 ρ 1 w 1 ρ 1 + w 2 ρ 2 + w f ρ f ,      (11)
ϕ 2 = w 2 ρ 2 w 1 ρ 1 + w 2 ρ 2 + w f ρ f ,      (12)
ϕ = ϕ 1 + ϕ 2 ,        (13)
More, in Equations (12)–(15), a group of similar variables are introduced, and w 1 , w 2 and w f are the 1st and 2nd nanoparticles and the base fluid masses, respectively.
η = ( c ν f ) 1 / 2 y , ψ = ( c ν f ) 1 / 2 x f ( η ) , H 1 = H 0 x l g ( η ) , H 2 = ( ν f c l 2 ) 1 / 2 H 0 g ( η ) , θ ( η ) = Τ Τ Τ w Τ , ϕ ( η ) = C C C w C .
Putting Equation (17) into dimensional administering, Equations (1)–(4), (6), and (11) tips to accomplish these dimensionless non-linear overseeing ODEs.
K A 1 f v + f + 1 k 1 f λ ρ f ρ h n f [ θ N r ϕ R b ξ ] + A 1 ( ( a c ) 2 + f f f 2 ) + β A 2 ( g 2 g g 1 ) = 0 ,
Λ g + f g f g = 0 ,
θ + k f k h n f A 4 Pr θ + k f k h n f α θ + ( k c p ) f k h n f A 4 Pr [ N b θ ϕ + N t θ 2 ] = 0 ,
Φ + S c f Φ + N t N b θ A Φ = 0 ,
ξ P e ( ( ξ + σ 1 ) ϕ + ξ ϕ ) Pr L b f ξ = 0 .
A 1 = A 2 × ( 1 w 1 ρ 1 + w 2 ρ 2 w 1 ρ 1 + w 2 ρ 2 + w f ρ f + w 1 ρ 1 + w 2 ρ 2 w 1 ρ 1 + w 2 ρ 2 + w f ρ f ρ s ρ f ) , A 2 = ( 1 w 1 ρ 1 + w 2 ρ 2 w 1 ρ 1 + w 2 ρ 2 + w f ρ f ) A 3 = A 2 [ w 2 ρ 2 w 1 ρ 1 + w 2 ρ 2 + w f ρ f { ( 1 w 1 ρ 1 w 1 ρ 1 + w 2 ρ 2 + w f ρ f ) + w 1 ρ 1 w 1 ρ 1 + w 2 ρ 2 + w f ρ f ( ρ β ) 1 ( ρ β ) f + w 2 ρ 2 w 1 ρ 1 + w 2 ρ 2 + w f ρ f ( ρ β ) 2 ( ρ β ) f } ] A 4 = ( 1 w 1 ρ 1 + w 2 ρ 2 w 1 ρ 1 + w 2 ρ 2 + w f ρ f + w 1 ρ 1 + w 2 ρ 2 w 1 ρ 1 + w 2 ρ 2 + w f ρ f ( ρ c p ) s ( ρ c p ) f )
Comparably, replacing Equation (17) with Equation (7) provides one with the following dimensionless boundary conditions:
f ( 0 ) = s , f ( 0 ) = 1 , g ( 0 ) = g ( 0 ) = 0 , θ ( 0 ) = 1 , ϕ ( 0 ) = 1 , ξ ( 0 ) = 1 , f ( ) a c , g ( ) 1 , θ ( ) 0 , ϕ ( ) 0 , ξ ( ) 0 .
Obviously, primes signify the separation regarding η . In the current issue, administering parameters, for example, Prandtl number ( Pr ) , suction or injection parameter ( s ) , penetrability boundary ( k 1 ) , magnetic boundary ( β ) , blended convection or lightness boundary ( λ ) , corresponding magnetic Prandtl number ( Λ ) , inertial boundary ( F 1 ) , thermophoretic boundary ( N t ) , Brownian movement boundary ( N b ) , coupled pressure boundary ( K ) , Schmidt number ( S c ) , and heat source boundary ( α ) , are characterized as:
Pr = ν f α f , s = V 0 ( c ν f ) 1 / 2 , k 1 = c k * ν f , β = μ e 4 π ρ f ( H 0 l c ) 2 , α = Q ν f c k f , Λ = η 0 ν f , F 1 = υ f c b k f λ = G r x Re x 2 , G r x = g β f ( T T ) x 3 ( ν f ) 2 , Re x = u w x ν f = x 2 c ν f , R = 4 σ 1 T 3 k * k f , N b = τ D B ( C w C ) ν f , N t = τ D T ( Τ s Τ 0 ) ν f Τ 0 , S c = ν f D B N r = ( ρ p ρ f ) Δ C ρ f β ( 1 C ) Δ T , P e = b * W c D m , L b = α D N , σ 1 = N N w N , R b = ( ρ m ρ f ) g γ Δ n ρ f β ( 1 C ) Δ T .
where the local Grashof number, bioconvection Rayleigh number, bioconvection Peclet number, bioconvection Lewis number, buoyancy ratio parameter, concentration difference parameter, and the local Reynolds number, respectively, are denoted by the following symbols G r x   , R b , P e , L b , N r , σ 1 , and   Re x . The suction and injection should be noted and correlate to the suction.

1.2. Physical Quantities of Interest

For the above model, the local Nusselt number ( N u x ) , local Sherwood number ( S h x ) , and skin friction coefficient ( c f x ) are clear as follows:
c f x = τ w ρ f u w 2 , N u x = x q w k f ( T T ) , S h x = x q m D B ( C C )
where
τ w = μ h n f ( u y ) y = 0 , q w = k h n f ( ( T y ) q r ) | y = 0 , q m = D B ( C y ) | y = 0  
From above, we have the dimensional form as
[ R e x ] 1 2 C f = ( 1 w 1 ρ 1 + w 2 ρ 2 w 1 ρ 1 + w 2 ρ 2 + w f ρ f ) f ( 0 ) , [ R e x ] 1 2 N u x = k h n f k f θ ( 0 ) , [ R e x ] 1 2 S h = ϕ ( 0 ) , N n x Re x 1 / 2 = ξ ( 0 ) .

1.3. Solution by HAM

Boundary conditions (Equation (20)) with Equations (15)–(19) have been resolved via HAM. Mathematica programming is utilized for this objective.
L f ( f ) = f v , L g ( g ) = g , L θ ( θ ) = θ , L Φ ( Φ ) = Φ   , L ξ ( ξ ) = ξ ,
The linear operators are presented as:
L f ( e 1 + e 2 η + e 3 η 2 + e 4 η 3 + e 5 η 4 ) = 0 , L g ( e 6 + e 7 η + e 8 η 2 ) = 0 , L θ ( e 9 + e 10 η ) = 0 , L Φ ( e 11 + e 12 η ) = 0   , L ξ ( e 13 + e 14 η ) = 0
The non-linear operatives are chosen as N f , N g , N θ   a n d     N ϕ   and are identified in the following systems:
  N f   [ f ( η ;     ζ ) , g ( η ;     ζ ) , θ ( η ;     ζ ) , Φ ( η ;     ζ ) , ξ ( η ;     ζ ) ] =     K A 1 f η η η η η + f η η η + 1 k 1 f η λ ρ f ρ h n f [ θ N r Φ R b ξ ]   + A 1 ( ( a c ) 2 + f f η η f η 2 ) + β A 2 ( g η 2 g g η η 1 ) ,                      
N g   [ f ( η ;     ζ ) , g ( η ;     ζ ) ] =   Λ g η η η + f g η η f η η g ,
N θ [ θ ( η ; ζ ) , Φ ( η ; ζ ) ] = ( 1 + 4 3 R ) θ η η + k f k h n f A 4 Pr θ η + k f k h n f α θ + ( k c p ) f k h n f A 4 Pr ( N b θ η ϕ η + N t θ η 2 ) ,
N Φ [ Φ ( η ; ζ ) , f ( η ; ζ ) , θ ( η ; ζ ) ] = Φ η η S c f Φ η + N t N b θ η η A Φ ,
N ξ [ ξ ( η ; ζ ) , f ( η ; ζ ) , ξ ( η ; ζ ) ] = ξ η η P e ( ( ξ + σ 1 ) Φ η η + ξ η Φ η ) Pr L b f ξ η .
Moreover, BCs are:
f ( η ; ζ ) η | η = 0 = 1 ,   f ( η ; ζ ) | η = 0 = s , 2 g ( η ; ζ ) 2 η | η = 0 = g ( η ; ζ ) | η = 0 = 0 , θ ( η ; ζ ) | η = 0 = 1 , Φ ( η ; ζ ) | η = 0 = 1 , ξ ( η ; ζ ) | η = 0 f ( η ; ζ ) η | η = = a c , g ( η ; ζ ) η | η = = 1 , θ ( η ; ζ ) | η = = 0 , Φ ( η ; ζ ) | η = = 0 , ξ ( η ; ζ ) | η = = 0 .
Here, ζ is the embedding parameter. ζ [ 0 , 1 ] is used to standardize the convergence of the solution of f , g , θ and ϕ By choosing ζ = 0   and   ζ = 1 [27], we have:
f ( η ; 1 ) = f ( η ) , g ( η ; 1 ) = g ( η )         θ ( η ; 1 ) = θ ( η )   ,           Φ ( η ; 1 ) = Φ ( η ) , ξ ( η ; 1 ) = ξ ( η ) ,
Develop Taylor’s series for f ( η ; ζ ) , g ( η ; ζ ) , θ ( η ; ζ ) , Φ ( η ; ζ ) , and ξ ( η ; ζ ) about the point ζ = 0 :
f ( η ; ζ ) = f 0 ( η ) + n = 1 f n ( η ) ζ n g ( η ; ζ ) = g 0 ( η ) + n = 1 g n ( η ) ζ n θ ( η ; ζ ) = θ 0 ( η ) + n = 1 θ n ( η ) ζ n Φ ( η ; ζ ) = Φ 0 ( η ) + n = 1 Φ n ( η ) ζ n ξ ( η ; ζ ) = ξ 0 ( η ) + n = 1 ξ n ( η ) ζ n
f n ( η ) = 1 n ! f ( η ; ζ ) η | p = 0 , g n ( η ) = 1 n ! g ( η ; ζ ) η | p = 0 ,                   θ n ( η ) = 1 n ! θ ( η ; ζ ) η | p = 0 ,             Φ n ( η ) = 1 n ! Φ ( η ; ζ ) η | p = 0 , ξ n ( η ) = 1 n ! ξ ( η ; ζ ) η | p = 0 .
Moreover, BCs are:
f ( 0 ) = s , f ( 0 ) = 1 , g ( 0 ) = g ( 0 ) = 0 , θ ( 0 ) = 1 , Φ ( 0 ) = 1 , f ( ) a c , g ( ) 1 , θ ( ) 0 , Φ ( ) 0 .
Now:
n f ( η ) = μ h n f / μ f ρ h n f / ρ f f n 1 v + f n 1 + 1 k 1 f n 1 λ ρ f ρ h n f [ θ n 1 N r Φ n 1 R b ξ n 1 ] + A 1 ( ( a c ) 2 + j = 0 w 1 f w 1 j f j f n 1 2 ) + β A 2 ( g n 1 2 j = 0 w 1 g w 1 j g j 1 ) ,
n g ( η ) = Λ g n 1 + j = 0 w 1 f w 1 j g j 2 j = 0 w 1 f w 1 j g j ,
n θ ( η ) = ( 1 + 4 3 R ) ( θ n 1 ) + k f k h n f A 4 Pr θ n 1 + k f k h n f α θ n 1 + ( k c p ) f k h n f A 4 Pr ( N b j = 0 w 1 θ w 1 j ϕ j + N t θ n 1 2 ) ,
n Φ ( η ) = Φ n 1 S c j = 0 w 1 f w 1 j Φ j + N t N b θ n 1 A Φ n 1 ,
n ξ ( η ) = ξ n 1 P e ( j = 0 w 1 ( ξ w 1 j + σ 1 ) ϕ j + j = 0 w 1 ξ w 1 j ϕ j ) Pr L b j = 0 w 1 f w 1 j ξ j .
Additionally,
χ n = { 0 ,   if   n 1 1 ,   if   n > 1 .

2. Results and Discussion

We now deliberate the consequences of the current exploration from the relevant sketched graphical features on velocity, temperature, and concentration profiles.

2.1. Velocity Profile

The impacts of β , k 1 , λ , K on velocity fields f ( η ) are demonstrated in Figure 2, Figure 3, Figure 4 and Figure 5. The impacts of the magnetic parameter ( β ) on the dimensionless velocity field f ( η ) is seen in Figure 2. It is clear that as ( β ) increases, f ( η ) decreases. The velocity profile in the domain decreases as the magnetic parameter ( β ) is increased, as seen in the sketch. Physically, this happens when the Lorentz force increases as the magnetic parameter increases, resulting in a decrease in liquid velocity, as seen in Figure 2. Figure 3 illustrates how increasing the porosity factor increases the system’s resistance. Physically, this causes a decrease in fluid flow due to increased frictional force. The impacts of the couple-stress parameter K on the velocity field is also seen in Figure 4. The velocity profile increases as the value of the couple-stress parameter K is raised, as can be seen in this graph. For broad values, however, the rise in velocity would be negligible. That is, large values would result in a pure viscous fluid. The effect of λ on f ( η ) is depicted in Figure 5. It is reasonable to assume that f ( η ) has higher values of λ based on Figure 5. Enlarging λ induces an enrichment of pliable force, which causes the boundary layer to extend, as seen in Figure 5.

2.2. Dimensionless Induced Magnetic Field along x   Direction g ( η )

The impacts of β , s , a c on the dimensionless-made magnetic field along the x   direction g ( η ) have been shown in Figure 6 and Figure 7. The impacts of the magnetic parameter ( β ) and the ( a / c ) suction or injection parameter on g ( η ) are shown in Figure 6. It is very clear that by growing the values of ( β ) , g ( η ) reduces. We notice that an increment in ( β ) summons an attractive field upgrade (Lorentz power increase) and furthermore dimensionless g ( η ) . Figure 6 also shows the effect of ( a / c ) on g ( η ) . As a result, a decrease in ( a / c ) leads to reducing g ( η ) . Figure 7 portrays the impact of the suction/injection parameter over the velocity profile g ( η ) . Increasing the suction or injection parameter initially reduces the g ( η ) profile, but at that point, as we grow in the direction of the center of the hydrodynamic edge stream, its pattern changes to the reverse direction.

2.3. Temperature Profile

Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12 show the effects of physical factors Pr , N b , N t , α on the temperature distribution θ ( η ) . The results of the Prandtl number on temperature and concentration profiles are seen in Figure 8. The temperature profile indicates a decrease with rising Pr values, as can be seen in Figure 8. Physically, raising the Prandtl number creates a decrease in thermal diffusivity, which is caused by a decrease in the temperature profile. The investigation of the heat source parameter α on θ ( η ) is depicted in Figure 9. When more heat is applied to the layer, the temperature of fluid particles in the whole domain increases, while in the case of a heat sink, the reverse result occurs, as seen in Figure 10. In the sink situation, the layer loses a lot of heat, dropping the temperature of the fluids. The effect of the thermophoresis boundary ( N t ) on the heat profile is shown in Figure 11. The temperature dispersion and the warm limit layer shows expanded conduct for the expansion of the thermophoresis boundary. Truly, the reformist idea of the thermophoretic boundary brings about an expansion of the thermophoretic power inside the liquid system, following improvement in the temperature profile and related limit layer. The effect of the Brownian dispersion boundary ( N b ) on the heat profile is shown in Figure 12. For the Brownian movement boundary, truly, the reformist nature warms the actual arrangement. This warming winds up moving nanoparticles from the colder extending sheet district to the quiet liquid locale.

2.4. Concentration Profile

The outcome of physical factors S c , N t , N b has been examined for concentration distribution Φ ( η ) in Figure 13, Figure 14 and Figure 15. The graphical results for S c are displayed in Figure 13. This is very clear that S c increases as the concentration profile Φ ( η ) declines. Therefore, the concentration decreases, S c and kinematic viscosity physically increase due to the reduction of molecular diffusion. Figure 14 shows that the growth in the parameter of the Brownian motion N b produces the reduction in the concentration profile Φ ( η ) of the fluid because nanoparticles move from the high concentration region to the region with less concentration. It is obvious that the increase in motion of the nanoparticles produces a high Brownian motion, and this irregular Brownian motion from high concentration areas to lower concentration regions reduces the momentum of the fluid. Figure 15 demonstrates that the escalation in the thermophoresis boundary N t produces the high concentration profile Φ ( η ) of the stream. These meet with zero at the boundary stream layer. The small variation in the thermophoresis boundary prompts fast movement in the liquid particles, making an abundance of heat energy and prompting a gigantic expansion in the focus dispersion. With an improvement in the calculation of N t , Figure 15 indicates a huge expansion in the focus dissemination.

2.5. Microorganism Distribution

The portrayed behavior of several estimations of the bioconvection L b and Peclet number in Figure 16 and Figure 17 indicates that improving the bioconvection and Peclet number induces a fast reduction in the thickness for motile microorganisms. That is, the thickness of motile microorganisms was decreased, and to be sure, by reinforcing the bioconvection Lewis number and Peclet number, the decrease in microorganisms’ dispersion was deciphered. This produces the thickness and limit layer thickness slumped for motile microorganisms by raising the worth in L b and P e . Figure 18 depicts the effect of the σ 1 on the rising parameters, which slows down the density of motile microorganisms. Figure 19, Figure 20, Figure 21 and Figure 22 show h-curves graphs of f ( 0 ) , θ ( 0 ) , Φ ( 0 ) , and ξ ( 0 ) , respectively.

2.6. Tables Discussion

Table 3 shows that c f x increased when the values of F 1 , k 1 , β increased. The c f x is decreased when the value of K increased. Table 4 shows that N u x increased when the values of R , α increased. The N u x decreased when the values of Pr , N b , N t increased. Table 5 shows that S h increased when the value of S c increased. The S h decreased when the values of N b , N t increased. Table 6 shows that N n x increased when the values of L b , σ 1 , P e increased. Table 7 is the comparison between the ND solves method and HAM method.

3. Conclusions

In this exploration, we analyzed the persistent laminar-blended convection of thick viscous and electrically leading sets of nanofluid cross breed Darcy–Forchheimer CuO-Cu/Blood stress stream close to the stagnation-point in the plane past a level permeable extending load up. This is utilized in biomedical fields, for example, the miniature roundabout framework’s stream elements and particularly in the inventory of medications. The fundamental partial differential equations (PDEs) are modified to a bunch of dimensionless ordinary differential equations (ODEs) with the assistance of reasonable comparability variables. These coupled ODEs are then solved by utilizing the Homotopy Analysis Method (HAM).
After detail study of the work, that the following conclusions were observed:
  • When increasing the value of the magnetic field parameter, the porosity factor velocity profiles decrease.
  • The velocity profile rises with a rise in the value of couple-stress parameter K .
  • The velocity profile displays a rising feature for greater values of λ .
  • With the enhancement of the strength of the thermophoresis parameter and the Brownian diffusion parameter, the temperature profile increases.
  • With the enhancement of the strength of the heat source ( s > 0 ) , the fluid temperature increases; on the other hand, an increase in the heat sink strength ( s < 0 ) decreases the temperature.
  • By increasing the value of the Prandtl factor, the fluid temperature decreases.
  • With an increase in the strength of the thermophoresis parameter and Brownian diffusion parameter, both have reverse impact on the concentration profile.
  • With an increase in the value of S c , the concentration profile decreases.
  • The Φ ( η ) portrayed a decreasing tendency with the rising number of L e .
  • The density of the moving microorganisms inside the fluid reduces for large values of P e .
  • The density number reduces to increase L b .

Author Contributions

Formal analysis, M.J., M.N., W.U.J., and R.J.; methodology, F.M.A. and M.J.; software, M.Z.; writing—original draft, M.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

There is no data availability regarding this research.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Not applicable.

Abbreviations

SymbolDiscription
uwstretching velocity (m/s)
u and vvelocity components along x- and y-axes (m/s)
H1 and H2magnetic field components along x- and y-axes
η0Magnetic diffusivity
βMagnetic field parameter
βvolumetric thermal expansion coefficient
μhnfviscosity of hybrid nanofluid
khnfthermal conductivity of hybrid nanofluid
ShSherwood number
ϕnanoparticle volume fraction
k*mean absorption coefficient
ssuction or injection parameter
w1, w2the first nanoparticle, the second nanoparticle
λmixed convection or buoyancy parameter
μeMagnetic permeability
αheat source parameter
ScSchmidt number
CfSkin friction coefficient
ρSdensity of nanoparticles
CwSurface concentration
TTemperature of the fluid (K)
CConcentration of the fluid
NnConcentration of the fluid
V0suction/injection velocity
F1Inertia coefficient
θDimentionless temperature (-)
ρfDensity (kg × m−3)
f’Dimentionless velocity (-)
RbBịoconvectịon Rayleigh number
μfDynamic viscosity (kg × m−1 × s−1)
PeBioconvection Peclet number
T∞Ambient temperature (K)
T0Reference temperature (K)
kDimensional permeability
gacceleration due to gravity (m/s2)
λ1diffusive constant parameter
Qvolumetric rate of heat generation/absorption
ρhnfDensity of hybrid nanofluid (kg × m−3)
(ρcp)hnfvolumetric heat capacity of hybrid nanofluid (m2 × s−2 × K−1)
RexReynolds number
ϕ1, ϕ2volume fraction of first and second nanoparticles
(cp)sspecific heat at constant pressure of nanoparticles (m2 × s−2 × K−1)
wfbase fluid masses
k1Non-dimensional permeability parameter
Ntthermophoretic parameter
ΛReciprocal magnetic prandtl number
Kcoupled stress parameter
NuxNusselt number
PrPrandtl number
TwSurface temperature (K)
NbBrownian motion parameter
lcharacteristic length (m)
σconcentration difference parameter
GrGrashof number
kThermal conductivity (W × m−1 × K−1)
NrBuoyancy ratio parameter
ΦDimentionless concentration (-)
ηSimilarity variable
νfKinematic viscosity (m2 × s−1)
(ρc)pEffective heat capacity of nanoparticles (m2 × s−2 × K−1)
LbBioconvection Lewis number

References

  1. Tiwari, A.; Chauhan, S.S. Effect of varying viscosity on two-fluid model of pulsatile blood flow through porous blood vessels: A comparative study. Microvasc. Res. 2019, 123, 99–110. [Google Scholar] [CrossRef]
  2. Deussen, A. Hyperthermie and hypothermie. Anaesthesist 2007, 56, 907–911. [Google Scholar] [CrossRef]
  3. Deniz, T.; Agalar, C.; Agalar, F.; Comu, F.M.; Caglayan, O.; Alpay, Y.; Saygun, O. The effect of hypothermia on splanchnic flows and lung in a two-hit hemorrhagic shock model. J. Surg. Res. 2010, 158, 121–126. [Google Scholar] [CrossRef]
  4. Bég, O.A.; Tripathi, D. Mathematica simulation of peristaltic pumping with double-diffusive convection in nanofluids: A bio-nano-engineering model. Proc. Inst. Mech. Eng. N 2011, 225, 99–114. [Google Scholar] [CrossRef]
  5. Kothandapani, M.; Prakash, J. Influence of heat source, thermal radiation, and inclined magnetic field on peristaltic flow of a hyperbolic tangent nanofluid in a tapered asymmetric channel. IEEE Trans. Nanobiosci. 2014, 14, 385–392. [Google Scholar] [CrossRef]
  6. Akbar, N.S. Metallic nanoparticles analysis for the blood flow in tapered stenosed arteries: Application in nanomedicines. Int. J. Biomath. 2016, 9, 1650002. [Google Scholar] [CrossRef]
  7. Bhatti, M.M.; Zeeshan, A.; Ellahi, R. Heat transfer analysis on peristaltically induced motion of particle-fluid suspension with variable viscosity: Clot blood model. Comput. Methods Prog. Biomed. 2016, 137, 115–124. [Google Scholar] [CrossRef] [PubMed]
  8. Dinarvand, S.; Rostami, M.N.; Dinarvand, R.; Pop, I. Improvement of drug delivery micro-circulatory system with a novel pattern of CuO-Cu/blood hybrid nanofluid flow towards a porous stretching sheet. Int. J. Numer. Methods Heat Fluid Flow 2019, 29, 4408–4429. [Google Scholar] [CrossRef]
  9. Majee, S.; Shit, G.C. Modeling and simulation of blood flow with magnetic nanoparticles as carrier for targeted drug delivery in the stenosed artery. Eur. J. Mech. B/Fluids 2020, 83, 42–57. [Google Scholar] [CrossRef]
  10. Varshney, G.; Katiyar, V.K.; Kumar, S. Numerical modeling of pulsatile flow of blood through a stenosed tapered artery under periodic body acceleration. J. Mech. Med. Biol. 2010, 10, 251–272. [Google Scholar] [CrossRef]
  11. Jinga, L.; Minb, K.B.; Baghbananc, A.; Zhaod, Z. Understanding coupled stress, flow and Jinga transport processes in fractured rocks. Geosyst. Eng. 2013, 16, 2–25. [Google Scholar] [CrossRef]
  12. Noorishad, J.; Ayatollahi, M.S.; Witherspoon, P.A. A Finite-Element Method for Coupled Stress and Fluid Flow Analysis in Fractured Rock Masses. Int. J. Rock Mech. Min. Sci. Geomech. 1982, 19, 185–193. [Google Scholar] [CrossRef]
  13. Ellahi, R.; Zeeshan, A.; Hussain, F.; Asadollahi, A. Peristaltic Blood Flow of Couple Stress Fluid Suspended with Nanoparticles under the Influence Chemical Reaction and Activation Energy. Symmetry 2019, 11, 276. [Google Scholar] [CrossRef] [Green Version]
  14. Choi, S.U.S.; Eastman, J.A. 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. [Google Scholar]
  15. Ghalambaz, M.; Behseresht, A.; Behseresht, J.; Chamkha, A.J. Effects of nanoparticles diameter and concentration on natural convection of the Al2O3-water nanofluids considering variable thermal conductivity around a vertical cone in porous media. Adv. Powder Technol. 2015, 26, 224–235. [Google Scholar] [CrossRef]
  16. Zaraki, A.; Ghalambaz, M.; Chamkha, A.J.; Ghalambaz, M.; De Rossi, D. Theoretical analysis of natural convection boundary layer heat and mass transfer of nanofluids: Effects of size, shape and type of nanoparticles, type of base fluid and working temperature. Adv. Powder Technol. 2015, 26, 935–946. [Google Scholar] [CrossRef]
  17. Reddy, P.S.; Chamkha, A.J. Dufour effects on MHD convective flow of Al2O3-water and TiO2-water nanofluids past a stretching sheet in porous media with heat generation/absorption. Adv. Powder Technol. 2016, 27, 1207–1218. [Google Scholar] [CrossRef]
  18. Reddy, P.S.; Sreedevi, P.; Chamkha, A.J. MHD boundary layer flow, heat and mass transfer analysis over a rotating disk through porous medium saturated by Cuwater and Ag-water nanofluid with chemical reaction. Adv. Powder Technol. 2017, 307, 46–55. [Google Scholar] [CrossRef]
  19. Mohebbi, R.; Rashidi, M.M. Numerical simulation of natural convection heat transfer of a nanofluid in an Lshaped enclosure with a heating obstacle. J. Taiwan Inst. Chem. Eng. 2017, 72, 70–84. [Google Scholar] [CrossRef]
  20. Ganesh, N.V.; Kameswaran, P.K.; Al-Mdallal, Q.M. Non-Linear thermal radiative Marangoni boundary layer flow of cAl2O3 nanofluids past a stretching sheet. J. Nanofluids 2018, 7, 1–7. [Google Scholar] [CrossRef]
  21. Tiwari, R.J.; Das, M.K. Heat transfer augmentation in a two-sided lid-driven differentially heated square cavity using nanofluids. Int. J. Heat Mass Transf. 2007, 50, 2002–2018. [Google Scholar] [CrossRef]
  22. Thameem Basha, H.; Sivaraj, R. Numerical simulation of blood nanofluid flow over three different geometries by means of gyrotactic microorganisms: Applications to the flow in a circulatory system. J. Mech. Eng. Sci. 2020, 1–20. [Google Scholar] [CrossRef]
  23. Bhatti, M.M.; Zeeshan, A.; Ellahi, R. Simultaneous effects of coagulation and variable magnetic field on peristaltically induced motion of Jeffrey nanofluid containing gyrotactic microorganism. Microvasc. Res. 2017, 110, 32–42. [Google Scholar] [CrossRef]
  24. Ahmed, S.E.; Mahdy, A. Laminar MHD natural convection of nanofluid containing gyrotactic microorganisms over vertical wavy surface saturated non-Darcian porous media. Appl. Math. Mech. 2016, 37, 471–484. [Google Scholar] [CrossRef]
  25. Kuznetsov, A.V. Non-oscillatory and oscillatory nanofluid bio-thermal convection in a horizontal layer of finite depth. Eur. J. Mech. B Fluids 2011, 30, 56–165. [Google Scholar] [CrossRef]
  26. Raju, C.S.K.; Mainul, M.; Sivasankar, T. The radiative flux of Casson fluid through a travelling coil packed with gyrotatic microorganisms. Adv. Powder Technol. 2017, 28, 575–583. [Google Scholar] [CrossRef]
  27. Rashad, A.M.; Chamkha, A.J.; Mallikarjuna, B.; Abdou, M.M.M. A nano-fluid-containing mixed bioconvection flow past the slender vertical cylinder. Front. Heat Mass Transf. 2018, 10, 21. [Google Scholar]
  28. Liao, S.J. An Explicit, Totally Analytic Approximate Solution for Blasius Viscous Flow Problems. Int. J. Non Linear Mech. 1999, 34, 759–778. [Google Scholar] [CrossRef]
Figure 1. Two-dimensional problem and coordinate system of the Schematic diagram.
Figure 1. Two-dimensional problem and coordinate system of the Schematic diagram.
Molecules 26 03954 g001
Figure 2. The effect of β on f ( η ) when k 1 = 0.3 , K = 0.7 , R b = 0.4 , N r = 2.0 , λ = 1.1 .
Figure 2. The effect of β on f ( η ) when k 1 = 0.3 , K = 0.7 , R b = 0.4 , N r = 2.0 , λ = 1.1 .
Molecules 26 03954 g002
Figure 3. The effect of k 1 on f ( η ) when β = 1.0 , K = 1.2 , R b = 0.4 , N r = 2.0 , λ = 1.1 .
Figure 3. The effect of k 1 on f ( η ) when β = 1.0 , K = 1.2 , R b = 0.4 , N r = 2.0 , λ = 1.1 .
Molecules 26 03954 g003
Figure 4. The effect of K on f ( η ) when β = 1.0 , k 1 = 1.2 , R b = 0.4 , N r = 2.0 , λ = 1.1 .
Figure 4. The effect of K on f ( η ) when β = 1.0 , k 1 = 1.2 , R b = 0.4 , N r = 2.0 , λ = 1.1 .
Molecules 26 03954 g004
Figure 5. The effect of λ on f ( η ) when β = 1.0 , k 1 = 1.2 , R b = 0.4 , N r = 2.0 .
Figure 5. The effect of λ on f ( η ) when β = 1.0 , k 1 = 1.2 , R b = 0.4 , N r = 2.0 .
Molecules 26 03954 g005
Figure 6. The effect of β , a c on g ( η ) when S = 1.0 , Λ = 1.0 .
Figure 6. The effect of β , a c on g ( η ) when S = 1.0 , Λ = 1.0 .
Molecules 26 03954 g006
Figure 7. The effect of s > 0 , s < 0 on g ( η )   when β = 1.5 , a c = 4.0 , Λ = 1.0 .
Figure 7. The effect of s > 0 , s < 0 on g ( η )   when β = 1.5 , a c = 4.0 , Λ = 1.0 .
Molecules 26 03954 g007
Figure 8. The effect of Pr on θ ( η ) when N b = 0.3 , α = 1.0 , N t = 0.2 .
Figure 8. The effect of Pr on θ ( η ) when N b = 0.3 , α = 1.0 , N t = 0.2 .
Molecules 26 03954 g008
Figure 9. The effect of α on θ ( η ) when N t = 0.2 , Pr = 12 , N b = 0.5 .
Figure 9. The effect of α on θ ( η ) when N t = 0.2 , Pr = 12 , N b = 0.5 .
Molecules 26 03954 g009
Figure 10. The effect of α on θ ( η ) when N t = 0.2 , Pr = 12 , N b = 0.5 .
Figure 10. The effect of α on θ ( η ) when N t = 0.2 , Pr = 12 , N b = 0.5 .
Molecules 26 03954 g010
Figure 11. The effect of N t on θ ( η ) when N b = 0.3 , α = 1.0 , Pr = 12 .
Figure 11. The effect of N t on θ ( η ) when N b = 0.3 , α = 1.0 , Pr = 12 .
Molecules 26 03954 g011
Figure 12. The effect of N b on θ ( η ) when N t = 0.2 , α = 1.0 , Pr = 7.5 , R = 0.8 .
Figure 12. The effect of N b on θ ( η ) when N t = 0.2 , α = 1.0 , Pr = 7.5 , R = 0.8 .
Molecules 26 03954 g012
Figure 13. The effect of S c on Φ ( η ) when N t = 0.3 , N b = 0.2 .
Figure 13. The effect of S c on Φ ( η ) when N t = 0.3 , N b = 0.2 .
Molecules 26 03954 g013
Figure 14. The effect of N b on Φ ( η ) when N t = 0.3 , S c = 1.0 .
Figure 14. The effect of N b on Φ ( η ) when N t = 0.3 , S c = 1.0 .
Molecules 26 03954 g014
Figure 15. The impact of N t on Φ ( η ) when N b = 0.3 , S c = 1.0 .
Figure 15. The impact of N t on Φ ( η ) when N b = 0.3 , S c = 1.0 .
Molecules 26 03954 g015
Figure 16. The effect of P e on ξ ( η ) when L b = 1.5 , σ 1 = 0.5 .
Figure 16. The effect of P e on ξ ( η ) when L b = 1.5 , σ 1 = 0.5 .
Molecules 26 03954 g016
Figure 17. The effect of L b on ξ ( η ) when P e = 0.4 , σ 1 = 0.5 .
Figure 17. The effect of L b on ξ ( η ) when P e = 0.4 , σ 1 = 0.5 .
Molecules 26 03954 g017
Figure 18. The effect of σ 1 on ξ ( η ) when P e = 0.4 , L b = 1.5 .
Figure 18. The effect of σ 1 on ξ ( η ) when P e = 0.4 , L b = 1.5 .
Molecules 26 03954 g018
Figure 19. The h-curve graph for f ( 0 ) .
Figure 19. The h-curve graph for f ( 0 ) .
Molecules 26 03954 g019
Figure 20. The h-curve graph for θ ( 0 ) .
Figure 20. The h-curve graph for θ ( 0 ) .
Molecules 26 03954 g020
Figure 21. The h-curve graph for Φ ( 0 ) .
Figure 21. The h-curve graph for Φ ( 0 ) .
Molecules 26 03954 g021
Figure 22. The h-curve graph for ξ ( 0 ) .
Figure 22. The h-curve graph for ξ ( 0 ) .
Molecules 26 03954 g022
Table 1. Nanoparticles and primary base fluid have the following thermophysical properties [8].
Table 1. Nanoparticles and primary base fluid have the following thermophysical properties [8].
Thermophysical PropertiesFluid Phase (Blood)Copper Oxide (CuO)Copper (Cu)
cp (J/kgK)3594533385
ρ (kg/m3)106365008933
k (w/mK)0.49217.65400
β × 105 (L/K)0.181.81.67
Particles size (nm)-295–25
Table 3. The effect on the Skin friction Re s 1 2 C f = ( 1 w 1 ρ 1 + w 2 ρ 2 w 1 ρ 1 + w 2 ρ 2 + w f ρ f ) f ( 0 ) , with the microorganism of different physical parameters. The impact of various physical parameters over.
Table 3. The effect on the Skin friction Re s 1 2 C f = ( 1 w 1 ρ 1 + w 2 ρ 2 w 1 ρ 1 + w 2 ρ 2 + w f ρ f ) f ( 0 ) , with the microorganism of different physical parameters. The impact of various physical parameters over.
F1k1βKNrλ ( 1 w 1 ρ 1 + w 2 ρ 2 w 1 ρ 1 + w 2 ρ 2 + w f ρ f ) f ( 0 ) ,
0.3 0.30.50.41.00.6 2.5313646
0.5 2.6768415
0.7 2.8147213
0.3 1.2310745
0.4 1.3324613
0.6 1.4345015
0.5 0.6156063
1.0 0.8367109
1.5 1.1348550
0.4 3.5157305
0.8 3.1516109
1.0 3.0357126
1.0 2.9461262
2.0 3.5240782
3.0 3.7180281
0.6 2.0146849
1.0 1.6130745
1.5 1.2154047
Table 4. The effect on the Nusselt number k h n f k f ( 1 + 4 3 R ) θ ( 0 ) with the microorganism of different physical parameters.
Table 4. The effect on the Nusselt number k h n f k f ( 1 + 4 3 R ) θ ( 0 ) with the microorganism of different physical parameters.
αPrNtNb k h n f k f ( 1 + 4 3 R ) θ ( 0 )
−0.58.50.70.5 3.6374583
−0.5 3.4584563
−1.0 3.1524204
8.5 2.3482139
9.5 2.6434058
10.5 2.8464381
0.2 2.7125861
0.4 2.5439308
0.6 2.1643615
0.5 0.9614792
0.7 1.4739615
0.9 1.7905246
Table 5. The effect on the Sherwood number Re x 1 2 S h x = ϕ ( 0 ) . with the microorganism of different physical parameters.
Table 5. The effect on the Sherwood number Re x 1 2 S h x = ϕ ( 0 ) . with the microorganism of different physical parameters.
NtNbSc ϕ ( 0 )
0.20.31.0 1.9352426
0.4 2.2568761
0.6 2.5917425
0.3 1.2518701
0.7 1.4507817
0.9 1.7591572
1.0 1.4523497
1.3 0.5248979
1.5 0.1935468
Table 6. The effect on the motile N n x Re x 1 / 2 = ξ ( 0 ) with the microorganism of different physical parameters.
Table 6. The effect on the motile N n x Re x 1 / 2 = ξ ( 0 ) with the microorganism of different physical parameters.
Peσ1Lb ξ ( 0 )
0.20.30.5 3.3467521
0.3 3.2674302
0.4 3.1459609
0.3 4.2670293
0.4 4.1523549
0.5 4.0795164
0.5 2.3895186
1.5 2.1384597
2.0 2.0634175
Table 7. A comparison between ND solves and HAM methods.
Table 7. A comparison between ND solves and HAM methods.
ηf′(η) θ(η)ϕ(η)ξ(η)
NDSolveHAMNDSolveHAMNDSolveHAMNDSolveHAM
0 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000
0.5 1.267861 1.267872 1.034338 1.034340 0.583358 0.583360 0.8434886 0.8434888
1.0 0.838512 0.838518 1.021331 1.021333 0.307213 0.307215 0.7330826 0.7330829
1.5 0.03184 0.03188 0.965956 0.965958 0.143647 0.143650 0.6444098 0.6444100
2.0 0.82835 0.82837 0.874652 0.874654 0.057064 0.057066 0.5616697 0.5616700
2.5 1.08151 1.08153 0.754801 0.754804 0.016656 0.016658 0.4770605 0.4770600
3.0 0.64058 0.64060 0.614284 0.614286 0.001988 0.001990 0.3896790 0.3896793
3.5 0.2365287 0.2365289 0.461086 0.461088 0.000735 0.000736 0.3012565 0.3012567
4.0 1.0290494 1.0290496 0.302960 0.302963 0.004461 0.004460 0.2111275 0.2111277
4.5 1.2528515 1.2528517 0.147122 0.147125 0.006107 0.006110 0.1139330 0.1139333
5.0 0.7499995 0.7499997 0 0 0 0 0 0
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Alharbi, F.M.; Naeem, M.; Zubair, M.; Jawad, M.; Jan, W.U.; Jan, R. Bioconvection Due to Gyrotactic Microorganisms in Couple Stress Hybrid Nanofluid Laminar Mixed Convection Incompressible Flow with Magnetic Nanoparticles and Chemical Reaction as Carrier for Targeted Drug Delivery through Porous Stretching Sheet. Molecules 2021, 26, 3954. https://doi.org/10.3390/molecules26133954

AMA Style

Alharbi FM, Naeem M, Zubair M, Jawad M, Jan WU, Jan R. Bioconvection Due to Gyrotactic Microorganisms in Couple Stress Hybrid Nanofluid Laminar Mixed Convection Incompressible Flow with Magnetic Nanoparticles and Chemical Reaction as Carrier for Targeted Drug Delivery through Porous Stretching Sheet. Molecules. 2021; 26(13):3954. https://doi.org/10.3390/molecules26133954

Chicago/Turabian Style

Alharbi, F. M., Muhammad Naeem, Muhammad Zubair, Muhammad Jawad, Wajid Ullah Jan, and Rashid Jan. 2021. "Bioconvection Due to Gyrotactic Microorganisms in Couple Stress Hybrid Nanofluid Laminar Mixed Convection Incompressible Flow with Magnetic Nanoparticles and Chemical Reaction as Carrier for Targeted Drug Delivery through Porous Stretching Sheet" Molecules 26, no. 13: 3954. https://doi.org/10.3390/molecules26133954

APA Style

Alharbi, F. M., Naeem, M., Zubair, M., Jawad, M., Jan, W. U., & Jan, R. (2021). Bioconvection Due to Gyrotactic Microorganisms in Couple Stress Hybrid Nanofluid Laminar Mixed Convection Incompressible Flow with Magnetic Nanoparticles and Chemical Reaction as Carrier for Targeted Drug Delivery through Porous Stretching Sheet. Molecules, 26(13), 3954. https://doi.org/10.3390/molecules26133954

Article Metrics

Back to TopTop