# Implementing the Galerkin Method Associated with the Shifted Vieta-Lucas Polynomials for Studying Numerically the Bionanofluid Flow Which Is Saturated by Gyrotactic Microorganisms over a Slippery Stretching Sheet

Authors: M. M. Khader; M. M. Babatin; Ahmed M. Megahed; A. Eid
Journal: Journal of Mathematics
 (2022)

---

## Abstract

Heat transfer is a critical function in many technical, industrial, home, and commercial structures. As a result, the purpose of this study is to investigate the effects of slip velocity and variable fluid characteristics on Casson bionanofluid flow across a stretching sheet that has been saturated by gyrotactic microorganisms. The suggested system will be converted to a computationally tractable form using the Galerkin method. The shifted Vieta-Lucas polynomials are then used as basis functions on the provided domain to solve the nonlinear system of ordinary differential equations that has been constructed (ODEs). The results are presented in the form of graphs and tables to assess the impact of the problem’s governing parameters. The estimated solutions produced by using the proposed techniques were physically acceptable and accurate. The current outcomes are confirmed by comparing them to the available literature. It appears that the temperature distribution is enhanced whereas the velocity distribution declines, caused by rising values of the magnetic parameter, slip parameter, and Casson parameter. Also, the local Nusselt number escalates with the strength of the viscosity parameter while the friction drag decays with the same parameter. In addition, the effectiveness and accuracy of the proposed method are satisfied by computing and the residual error function.


---

## Body

## 1. Introduction

Many technological and industrial operations, such as the manufacture of glass fibers, polymer refining, and the extrusion process of aerodynamics, rely heavily on nanofluid flow through stretched sheets [1, 2]. Nanofluids are also very useful in biomedical research. Cancer therapies, nanocryosurgery, and magnetic resonance imaging are just a few examples [3]. The suspended stability of nanoparticles is commonly found to be greatly increased when gyrotactic bacteria are present. As a result of their numerous practical applications in assorted disciplines of biotechnology and applicable science, the innovative research of nanofluid flow encompassing gyrotactic microorganisms as well as heat (mass) transmission across a stretching surface has gotten a lot of attention from the research community and scientists. Various forms of microorganisms exist, including gravitaxis, oxytocic, chemotaxis, and gyrotactic. Microorganisms have been employed to make a variety of industrial and commercial products, including fertilizers, medicine delivery, and waste-derived biofuel [4]. In recent years, many researchers have written about these remarkable phenomena [5–8]. On a stretched surface, recently, Hosseinzadeh et al. [9] investigated the MHD nanofluid flow in conjunction with thermal radiation, microorganisms, and nanoparticles. Sankad et al. [10] very recently examined a semi-intuitive and numerical solution for the system governing the fluid flow problem containing gyrotactic microorganisms across a linearly expanding sheet.In fact, most nonlinear ordinary differential equations (ODEs) have no exact solutions, so numerical and approximate techniques are the semiunique way to solve these types of ODEs. One of the most useful tools to simulate the ODEs is the spectral method [11]. The most famous advantage of these methods is their capability to generate accurate outcomes with a very small degree of error of freedom. The orthogonality property of some important polynomials as the shifted Vieta-Lucas polynomials is used to approximate functions on their domain. These polynomials have a main and important role in these methods for ODEs [12]. Many researchers used and implemented these polynomials to solve numerically many problems such as in [13] to solve the 1D and 2D nonlinear generalized Benjamin-Bona-Mahony-Burgers equation, in [14] to solve the 1D and 2D sinh-Gordon equation, and in [15] to get the approximate solution for certain types of fourth-order BVPs. In the same field of study, many kinds of research are also introduced [16].The major objective of this study is to use the Galerkin method to numerically investigate the results acquired while investigating a non-Newtonian Casson dissipative nanofluid flow containing gyrotactic microorganisms and influenced by a magnetic field and slip velocity phenomenon. The novelty of the proposed method is that it gets a more accurate numerical solution than methods that appear in the literature. For the first time, the Galerkin method and the shifted Vieta-Lucas polynomials are combined to take advantage of the properties of each of them in approximation on one hand and in increasing the applicability and accuracy of the proposed technique on the other hand, which in turn gives good and closest solutions for the real solutions for the problem under study. As a result of applying the proposed technique to the Casson bionanofluid flow caused by a stretched sheet saturated by gyrotactic microorganisms in the presence of slip velocity, the provided study is novel according to the literature review.


## 2. Description of the Problem

A two-dimensional non-Newtonian Casson bionanofluid saturated by gyrotactic microorganisms is addressed with slip velocity, magnetic field, and viscous dissipation. The fluid flow is created by the stretching of an elastic impermeable sheet at a velocity ofuw=ax. During this examination, the x−axis is focused on the direction of flow, but the y−axis is perpendicular to it. Furthermore, the microbe concentration, fluid concentration, and temperature at the sheet surface are represented by Nw, Cw, and Tw (Figure 1).Figure 1
Flow diagram of the problem.In the case of varying properties, the equations that regulate the physical model can be written as [10](1)∂u∂x+∂v∂y=0,(2)u∂u∂x+v∂u∂y=1ρ∞∂∂yμT1+1β∂u∂y−σ0B02ρ∞u,(3)u∂T∂x+v∂T∂y=1ρ∞cp∂∂yκT∂T∂y+τDBC∂C∂y∂T∂y+DTT∞∂T∂y2+1ρ∞cpμT1+1β∂u∂y2,(4)u∂C∂x+v∂C∂y=∂∂yDBC∂C∂y+DTT∞∂2T∂y2,(5)u∂N∂x+v∂N∂y+b˜WcCw−C∞∂∂yN∂C∂y=∂∂yDmC∂N∂y,where the velocity components in the x and y directions are u and v, respectively. Likewise, the magnetic field’s strength is denoted by B0, T is the temperature of the nanofluid, DT is the thermophoretic diffusion coefficient, β is the Casson parameter, ρ∞ is the density of the nanofluid at the ambient, and Dm is the diffusivity of microorganisms. Also, DB is the Brownian diffusion coefficient, b˜ is the chemotaxis constant, τ is the ratio of the nanoparticle’s effective heat capacitance to that of the base fluid, and Wc represents the maximum speed at which a cell can swim. Further, the microbe concentration, fluid concentration, and the temperature all at the ambient of the surface are represented by N∞,C∞, and T∞. On the other hand, the following are the pertinent physical boundary conditions:(6)u=ax+λ1μTμ∞1+1β∂u∂y,T=Tw,v=0,C=Cw,N=Nwaty=0,(7)u⟶0,T⟶T∞,C⟶C∞,N⟶N∞,asy⟶∞,wherea>0 is the rate of stretching for the sheet, λ1 is the velocity slip coefficient, and μ∞ is the Casson nanofluid viscosity at the ambient. Furthermore, the fluctuation in viscosity, thermal conductivity, and nanoparticle and microorganism diffusivities are represented with the corresponding formulae functions, based on works by Noghrehabadi and Behseresht [17] and Megahed [18].(8)μT=μ∞Exp−αT−T∞Tw−T∞,κT=κ∞1+ε1T−T∞Tw−T∞,(9)DBC=d11+ε2C−C∞Cw−C∞,DmC=d21+ε3C−C∞Cw−C∞,(10)η=yaρ∞μ∞,u=axf′η,v=−aρ∞μ∞fη,where α is the viscosity parameter, ε1 is the thermal conductivity parameter, ε2 is the mass diffusivity parameter, and ε3 is the microorganism diffusivity parameter. Now we will use the dimensionless variables listed below [10]:(11)θη=T−T∞Tw−T∞,ϕη=C−C∞Cw−C∞,χη=N−N∞Nw−N∞.To get the nonlinear ordinary differential equations given below, we plug (8)–(11) into (2)–(7).(12)1+1βf‴−αθ′f″e−αθ+ff″−f′2−Mf′=0,(13)1Pr1+ε1θθ″+ε1θ′2+fθ′+Ntθ′2+Nb1+ε2ϕθ′ϕ′+Ec1+1βf″2e−αθ=0,(14)1+ε2ϕϕ″+Scfϕ′+ε2ϕ′2+NtNbθ″=0,(15)1+ε3ϕχ″+ε3ϕ′χ′+Sbfχ′−Peϕ′χ′+σ+χϕ″=0.The converted boundary conditions are as follows:(16)f′0=1+λ1+1βf″0e−αθ0,f0=0,θ0=1,ϕ0=1,χ0=1,(17)f′∞⟶0,θ∞⟶0,ϕ∞⟶0,χ∞⟶0,where the parameters present in equations (12)–(17) can be defined as follows: M=σ0B02/aρ∞ is the magnetic parameter, λ=λ1aρ∞/μ∞ is the slip velocity parameter, Pr=μ∞cp/κ∞ is the Prandtl number, Nt=τDTTw−T∞/ν∞T∞ is the thermophoresis parameter, Nb=τd1Cw−C∞/ν∞ is the Brownian motion parameter, Ec=uw2/cpTw−T∞ is the Eckert number, Sc=ν∞/d1 is the Schmidt number, Sb=ν∞/d2 is the bioconvection Schmidt number, Pe=b˜Wc/d2 is the bioconvection Peclet number, and σ=N∞/Nw−N∞ is the bioconvection constant. Further, the local skin friction factor Cfx, local Nusselt number Nux, local Sherwood number Shx, and local density number of motile microorganism Nnx are all important engineering factors. They are formulated as (18)Rex1/2Cfx=−1+1βf″0e−α,Rex−1/2Nux=−1+ε1θ′0,Rex−1/2Shx=−ϕ′0,Rex−1/2Nnx=−χ′0,where Rex=uwx/ν∞ is the local Reynolds number.


## 3. Basic Concepts of the Proposed Methods

Based on Sturm-Liouville’s theory, some special polynomials are created through well-known linear algebra methods. These polynomials are often found to use the power series method to solve the ODEs. It is usually demonstrated that they may be constructed using a generating function, the Rodrigues formula, for each specific polynomial and, ultimately, a contour integral. These polynomials are of great importance in the theory of approximation, mathematical physics, engineering, the theory of mechanical quadrature, and so forth. This section shows a brief description of the most important applications for one special polynomial, namely, shifted Vieta-Lucas polynomials. In addition, we will introduce some concepts about the Galerkin method.


### 3.1. The Galerkin Method

This subsection introduces the Galerkin method which is based on the approximation concept of the exact solution utilizing a set of basis polynomials. The mathematical theory behind the Galerkin technique proves to be immensely useful in laying the groundwork for both qualitative and quantitative analysis of the ODEs, which have been studied extensively since the mid-nineteenth century. Many researchers used and implemented these polynomials to solve numerically many problems such as in [19] to get the approximate solution for linear second-order ordinary differential equations with mixed boundary conditions, and in [20] they were used to solve the wave propagation problems. To clarify this method, we consider the following operator equation [20]:(19)ℒu=g,inΩ.ℒ is usually a differential operator acting on u. The approximate solution for (19) is as follows:(20)ux≈uNx=∑i=0Nγiψix,where the coordinate function ψix is the approximating function and N is a natural number. Note that the weights γi should satisfy the criteria that the difference (the residual of the approximation which is denoted by R) is minimal (i.e., nonzero):(21)R=ℒuNx−g=ℒ∑i=0Nγiψix−g≠0.The parametersγi’s in the weighted residual technique are derived by requiring the residual R to disappear in the weighted-integral sense; that is:(22)∫ΩψixRx,γidxdy=0,i=0,1,…,N.For more details about the Galerkin method, the reader is urged to examine [19–21].


### 3.2. Vieta-Lucas and the Shifted Vieta-Lucas Polynomials

In this part of the paper, we are presenting the basic definitions of the Vieta-Lucas and the shifted Vieta-Lucas polynomials and their notations and properties that we will use in our study and they are necessary to reach our goal [22].We are researching for a class of orthogonal polynomials which lies at the heart of our research. Using the recurrence relations and analytical formula of these polynomials, they can be generated to construct a new family of orthogonal polynomials that will be known as Vieta-Lucas polynomials.The Vieta-Lucas polynomialsVLmz of degree m∈ℕ0 are defined by the following relation [22]:(23)VLmz=2cosmξ,ξ=arccosz2,ξ∈0,π,z≤2.It is easy like other famous functions; we can prove that these Vieta-Lucas polynomials satisfy the following recurrence formula:(24)VLmz=zVLm−1z−VLm−2z,m=2,3,…,VL0z=2,VL1z=z.Using the transformationz=4/ℓx−2, we can generate from the family of Vieta-Lucas polynomials a new class of orthogonal polynomials on the interval 0,ℓ, which are the orthogonal family of the shifted Vieta-Lucas polynomials, and it will be denoted by VLmsx and can be obtained as follows:(25)VLmsx=VLm4ℓx−2.The shifted Vieta-Lucas polynomialsVLmsx satisfy the following recurrence relation:(26)VLm+1sx=4ℓx−2VLmsx−VLm−1sx,m=1,2,…,where VL0sx=2,VL1sx=4/ℓx−2. Also, we find VLms0=2−1m and VLmsℓ=2,m=0,1,2,….For more details about the shifted Vieta-Lucas polynomials, see [22, 23].


## 3.1. The Galerkin Method

This subsection introduces the Galerkin method which is based on the approximation concept of the exact solution utilizing a set of basis polynomials. The mathematical theory behind the Galerkin technique proves to be immensely useful in laying the groundwork for both qualitative and quantitative analysis of the ODEs, which have been studied extensively since the mid-nineteenth century. Many researchers used and implemented these polynomials to solve numerically many problems such as in [19] to get the approximate solution for linear second-order ordinary differential equations with mixed boundary conditions, and in [20] they were used to solve the wave propagation problems. To clarify this method, we consider the following operator equation [20]:(19)ℒu=g,inΩ.ℒ is usually a differential operator acting on u. The approximate solution for (19) is as follows:(20)ux≈uNx=∑i=0Nγiψix,where the coordinate function ψix is the approximating function and N is a natural number. Note that the weights γi should satisfy the criteria that the difference (the residual of the approximation which is denoted by R) is minimal (i.e., nonzero):(21)R=ℒuNx−g=ℒ∑i=0Nγiψix−g≠0.The parametersγi’s in the weighted residual technique are derived by requiring the residual R to disappear in the weighted-integral sense; that is:(22)∫ΩψixRx,γidxdy=0,i=0,1,…,N.For more details about the Galerkin method, the reader is urged to examine [19–21].


## 3.2. Vieta-Lucas and the Shifted Vieta-Lucas Polynomials

In this part of the paper, we are presenting the basic definitions of the Vieta-Lucas and the shifted Vieta-Lucas polynomials and their notations and properties that we will use in our study and they are necessary to reach our goal [22].We are researching for a class of orthogonal polynomials which lies at the heart of our research. Using the recurrence relations and analytical formula of these polynomials, they can be generated to construct a new family of orthogonal polynomials that will be known as Vieta-Lucas polynomials.The Vieta-Lucas polynomialsVLmz of degree m∈ℕ0 are defined by the following relation [22]:(23)VLmz=2cosmξ,ξ=arccosz2,ξ∈0,π,z≤2.It is easy like other famous functions; we can prove that these Vieta-Lucas polynomials satisfy the following recurrence formula:(24)VLmz=zVLm−1z−VLm−2z,m=2,3,…,VL0z=2,VL1z=z.Using the transformationz=4/ℓx−2, we can generate from the family of Vieta-Lucas polynomials a new class of orthogonal polynomials on the interval 0,ℓ, which are the orthogonal family of the shifted Vieta-Lucas polynomials, and it will be denoted by VLmsx and can be obtained as follows:(25)VLmsx=VLm4ℓx−2.The shifted Vieta-Lucas polynomialsVLmsx satisfy the following recurrence relation:(26)VLm+1sx=4ℓx−2VLmsx−VLm−1sx,m=1,2,…,where VL0sx=2,VL1sx=4/ℓx−2. Also, we find VLms0=2−1m and VLmsℓ=2,m=0,1,2,….For more details about the shifted Vieta-Lucas polynomials, see [22, 23].


## 4. The Proposed Numerical Scheme

In this section, we are going to implement the Galerkin method associated with the shifted Vieta-Lucas polynomials to solve numerically the proposed model (12)–(15) with its boundary conditions (16) and (17).In the Galerkin method, we take a finite space solution’s domain0,η∞ [24].We will apply the proposed numerical method with the following steps:Step 1.
We use the auxiliary parameterf′=Υ to reduce (12) of the second order. So we can rewrite system (12)–(15) as follows:(27)f′=Υ,(28)Υ″−αθ′Υ′+β1+βeαθfΥ′−Υ2−MΥ=0,(29)1+ε1θθ″+ε1θ′2+Prfθ′+Ntθ′2+Nb1+ε2ϕθ′ϕ′+Ec1+1βΥ′2e−αθ=0,(30)1+ε2ϕϕ″+Scfϕ′+ε2ϕ′2+NtNbθ″=0,(31)1+ε3ϕχ″+ε3ϕ′χ′+Sbfχ′−Peϕ′χ′+σ+χϕ″=0.
The associated transformed boundary conditions become(32)Υ0=1+λ1+1βΥ′0e−αθ0,f0=0,θ0=ϕ0=χ0=1,(33)Υη∞=θη∞=ϕη∞=χη∞=0.Step 2.
The idea behind the proposed method assumes that the approximate solution for (21)–(33) can be provided as a linear combination of shifted the Vieta-Lucas polynomials as follows:(34)f¯η=∑i=0mci1VLisη,Υ¯η=∑i=0mci2VLisη,θ¯η=∑i=0mci3VLisη,ϕ¯η=∑i=0mci4VLisη,χ¯η=∑i=0mci5VLisη.Step 3.
The values of the constants (weights)ci1,ci2,ci3,ci4,ci5,i=0,1,…,m can be given by applying the boundary conditions (32) and (33) and the Galerkin approach to solve system (27)–(31) as follows:(35)∫0η∞f′−ΥVLjsηdη=0,∫0η∞Υ″−αθ′Υ′+β1+βeαθfΥ′−Υ2−MΥVLjsηdη=0,∫0η∞1+ε1θθ″+ε1θ′2+Prfθ′+Ntθ′2+Nb1+ε2ϕθ′ϕ′+Ec1+1βΥ′2e−αθVLjsηdη=0,∫0η∞1+ε2ϕϕ″+Scfϕ′+ε2ϕ′2+NtNbθ″VLjsηdη=0,∫0η∞1+ε3ϕχ″+ε3ϕ′χ′+Sbfχ′−Peϕ′χ′+σ+χϕ″VLjsηdη=0.
This is equivalent to(36)∫0η∞f¯′−Υ¯VLjsηdη=0,(37)∫0η∞−Υ¯′VLj′sη+−αθ¯′Υ¯′+β1+βeαθ¯f¯Υ¯′−Υ¯2−MΥ¯VLjsηdη=−Υ¯′η∞VLjsη∞+Υ¯′0VLjs0,(38)∫0η∞−1+ε1θ¯θ¯′VLj′sη+Prfθ¯′+Nb1+ε2ϕθ¯′ϕ′+Ec1+1βΥ¯′2e−αθ¯VLjsηdη=−1+ε1θ¯η∞θ¯′η∞VLjsη∞+1+ε1θ¯0θ¯′0VLjs0,(39)∫0η∞−1+ε2ϕ¯ϕ¯′+NtNbθ¯′VLj′sη+Scf¯ϕ¯′VLjsηdη=−1+ε2ϕ¯η∞ϕ¯′η∞VLjsη∞+1+ε2ϕ¯0ϕ¯′0VLjs0−NtNbθ¯′η∞VLjsη∞+NtNbθ¯′0VLjs0,(40)∫0η∞−1+ε3ϕ¯χ¯′+σ+χ¯ϕ¯′VLj′sη+Sbf¯χ¯′−P+1ϕ¯′χ¯′VLjsηdη−1+ε3ϕ¯η∞χ¯′η∞VLjsη∞+1+ε3ϕ¯0χ¯′0VLjs0−σ+χ¯η∞ϕ¯′η∞VLjsη∞+σ+χ¯0ϕ¯′0VLjs0.Step 4.
After considerable rearrangement, with the result of substituting (36) with m=5 (as a special case to show the method) into the system of equations (36)–(38) and in the associated boundary conditions (32) and (33), we can write the nonlinear equations (36)–(40) in a matrix notation by using the Newton iteration method (NIM) as follows:(41)Ψp+1=Ψp−J−1ΨpFΨp,where Ψp=c01p,c11p,…,c51p;c02p,c12p,…,c52p;c03p,c13p,…,c53p;c04p,c14p,…,c54p;c05p,c15p,…,c55pT, FΨs is the vector which represents the nonlinear equations, and J−1Ψs is the inverse of the Jacobian matrix [25].
Note that, here in constructing the nonlinear system of equations, we take into account the factsVLms0=2−1m and VLmsη∞=2,m=0,1,2,…, in the right hand side of the system of equations (36)–(40).Step 5.
The values of all coefficientscjk,k=1,2,3,4,5;j=0,1,2,…,m are acquired using the direct approach of solving the obtained nonlinear system of equations by the NIM and by substituting these values into formulae (34), the approximate solution for the proposed model (12)–(17) in the specified range 0,8 is generated.


## 5. Results and Discussion

Table1 compares the skin friction values to those previously published. It should be emphasized that the obtained findings are in line with the previously reported outcomes. This lends credibility to the current mathematical formulation’s authenticity and precision, as well as the numerical solutions.Table 1
Comparison ofCfxRe1/2 with the results of Hayat et al. [26] when α=λ=0 and β⟶∞.
MHayat et al. [26]Present work0.01.000001.0000000001.01.414211.4142079815.02.449482.44947879410.03.316623.31661692750.07.141427.141419850Numerical values are sketched in Figures2–11 to illustrate the findings, and this part includes a full explanation of the impacts of the pertinent parameters on the velocity, temperature, concentration, and motile bacteria density profiles. The skin friction coefficient, local Nusselt number, local Sherwood number, and local density number of motile microorganisms are also calculated in this part with the same output responses. The implications of the dimensionless velocity f′η, dimensionless temperature θη, dimensionless concentration ϕη, and dimensionless microorganism χη due to the magnetic number M are depicted in Figure 2. Clearly, when the magnetic number is enhanced, the velocity declines rapidly, whereas for dimensionless temperature, dimensionless concentration, and dimensionless microorganisms, the opposite pattern is observed. Physically, when the magnetic number is present, a reduction in the boundary layer thickness and a diminished behavior for the velocity of the nanofluid along the sheet occur. Physically, because the Lorentz force is proportional to the magnetic field, the larger the Lorentz force, the lower the fluid velocity and the higher the temperature.Figure 2
(a) Impact ofM on f′η. (b) Impact of M on θη. (c) Impact of M on ϕη. (d) Impact of M on χη.
(a)(b)(c)(d)Figure 3
(a) Impact ofλ on f′η. (b) Impact of λ on θη. (c) Impact of λ on ϕη. (d) Impact of λ on χη.
(a)(b)(c)(d)Figure 4
(a) Impact ofα on f′η. (b) Impact of α on θη. (c) Impact of α on ϕη. (d) Impact of α on χη.
(a)(b)(c)(d)Figure 5
(a) Impact ofβ on f′η. (b) Impact of β on θη. (c) Impact of β on ϕη. (d) Impact of β on χη.
(a)(b)(c)(d)Figure 6
(a) Impact ofε1 on θη. (b) Impact of Ec on θη.
(a)(b)Figure 7
(a) Impact ofNt on ϕη. (b) Impact of Nt on χη.
(a)(b)Figure 8
(a) Impact ofNt on θη. (b) Impact of Nb on θη.
(a)(b)Figure 9
(a) Impact ofNb on ϕη. (b) Impact of Nb on χη.
(a)(b)Figure 10
(a) Impact ofε2 on ϕη. (b) Impact of ε2 on χη.
(a)(b)Figure 11
(a) Impact ofε2 on θη. (b) Impact of ε3 on χη.
(a)(b)Figure3 conveys the impacts of slip velocity parameter λ on velocity f′η, temperature θη, concentration ϕη, and dimensionless microorganism χη distributions, respectively. The dimensionless microorganism, dimensionless temperature, and dimensionless concentration of the Casson nanofluid rise as the slip velocity parameter is increased, whereas the contrary trend occurs for the dimensionless velocity with the same parameter. The resistance force that was created to act towards raising the fluid temperature distribution, but acts against the velocity field, demonstrates the physical importance of slip velocity parameter on the fluid velocity characteristics.Figure4 depicts the behavior of the viscosity parameter α for velocities, temperature, concentration, and microbe field in Casson nanofluid motion. The viscosity parameter has a strong connection with the shear stress, and as a result, with greater values of the viscosity parameter, velocity increases besides the sheet, while the contrary trend happens away from the sheet.Furthermore, when the same parameter is improved, the temperature, concentration, and microbe field all increase because a boost in viscosity creates higher resistance in the nanofluid layers which raises the temperature field, as well as concentration and microbe field. To assess the impact of the Casson parameterβ on the flow and heat mass transfer as well as the microbe field, Figure 5 is depicted. The velocity field drops as the Casson parameter is increased, whereas the temperature, concentration, and microbe fields show the opposite pattern. Physically, when the fluid reduces to the Newtonian type β⟶∞, it is possible to accomplish a robust warming process throughout the thermal boundary layer.Figure6 shows the temperature distribution as a function of the Eckert number Ec and the thermal conductivity parameter ε1. It is obvious from this diagram that increasing the value of both parameters enhances the thickness of the thermal boundary layer as well as the temperature distribution. Physically, when the dissipation impact is addressed, the distribution of the temperature is at a higher level due to internal heating related to the dissipation phenomenon. Physically, increasing the Eckert number reduces the temperature difference between the surface Tw and the ambient temperature T∞, improving temperature dispersion. Figure 7 depicts the distributions of dimensionless concentration ϕη and dimensionless microorganisms χη for a horizontal sheet with various thermophoresis parameter values Nt.As can be observed in this graph, as the thermophoresis parameter is increased, both the dimensionless concentration and the dimensionless microorganisms rise. Figure8 also depicts the temperature profiles for the thermophoresis Nt and Brownian motion Nb parameters, respectively. When both parameters are increased, the fluid temperature θη in the boundary layer medium is observed to rise. Physically, as the fluid dynamic viscosity coefficient ν∞ grows, the magnitude of both the Brownian motion and thermophoresis parameters rises, and therefore the fluid temperature rises.When the Brownian motion parameterNb is adjusted for Casson nanofluid flow, Figure 9 shows the variation in concentration function ϕη and dimensionless microorganisms function χη plotted against independent variable η. Increasing Nb suppresses both the concentration distribution and the spread of dimensionless microbes.Figure10 depicts the behavior of the concentration function ϕη and the dimensionless microorganisms function χη against η for the mass diffusivity parameter ε2, respectively. It may be concluded that the mass diffusivity parameter has a greater impact on the concentration function than on the function of dimensionless microorganisms. As can be seen, raising the mass diffusivity parameter ε2 resulted in more bacterial spreading and a more even distribution of concentrations.Figure11(a) depicts the changing of the dimensionless temperature θη for various values of the mass diffusivity parameter ε2. For a constant amount of η, the dimensionless temperature θη increases slightly with the mass diffusivity parameter. This occurs since boosting the mass diffusivity parameter causes inductive heat to be generated within the boundary layer, resulting in a higher temperature profile. Figure 11(b), on the other hand, shows how the microbe diffusivity parameter ε3 affects the dimensionless microorganisms profile χη. It is worth noting that the dimensionless microorganisms profile in the Casson nanofluid is heavily influenced by the microorganism diffusivity parameter ε3. We discovered that raising ε3 raises the profile of nanofluid bacteria.To validate the accuracy of the approximate method, we evaluated through Table2 the residual error function (REF) [27] of the current method with the values of parameters M=1.0,λ=α=ε1=ε2=ε3=σ=Ec=Nt=0.2,Nb=0.5, β=Sb=0.7, with Pr=Sc=2.0,Pe=0.6. These values show the thoroughness of the proposed method in this article and confirm that the current method gives better accuracy. In addition, we studied the effect of the order of the approximation by the Galerkin method associated with the shifted Vieta-Lucas polynomials on the skin friction and the Nusselt number through Table 3. From this table, we can get more accurate values of these parameters by increasing the order of the approximation.Table 2
Values for the REF in the current method.
ηRef offηRef ofθηRef ofϕηRef ofχη0.04.753951E − 076.951753E − 075.156874E − 071.963741E − 071.04.456123E − 064.789123E − 061.654123E − 062.654123E − 052.05.852741E − 061.002140E − 052.741258E − 064.021345E − 063.01.951321E − 073.852123E − 063.963742E − 076.852789E − 064.06.853214E − 062.456951E − 067.021541E − 079.964578E − 065.07.021340E − 060.953214E − 079.115500E − 065.123045E − 056.00.654789E − 076.012548E − 062.551423E − 052.015975E − 077.08.963741E − 057.654789E − 050.852123E − 063.975320E − 078.06.954512E − 062.654123E − 067.789357E − 078.012548E − 06Table 3
The order of approximationN vs. the skin friction −f″0 and the Nusselt number −θ′0 when β=0.7,M=1.0, Nb=0.5,Pe=0.6,Pr=Sc=2.0, and Nt=Ec=α=λ=σ=ε1=ε2=ε3=0.2.
The order of approximation−f″0−θ′050.661485200.25098510100.662299810.25100987150.662365290.25115874200.662371980.25117981250.662372110.25118215300.662372110.25118215The drag force, heat transfer coefficient, mass transfer coefficient, and microorganism transfer coefficient are now represented in nondimensional form asCfxRex1/2, NuxRex−1/2, ShxRex−1/2, and NnxRex−1/2, respectively, versus various values of the controlling parameters in Table 4. The heat transfer coefficient improves as the magnetic number, slip velocity parameter, and viscosity parameter increase; however the mass transfer coefficient and microorganism transfer coefficient exhibit the reverse trend for the same parameters. Both the thermal conductivity and mass diffusivity parameters appear to have a little effect on the local skin friction coefficient, as well as a gradual decrease in the local Sherwood number, local Nusselt number, and microorganism transfer rate at the wall. Furthermore, it is undeniable that raising the Brownian motion parameter of Eckert number lowers the local Nusselt number, whereas the local skin friction coefficient, the local Sherwood number, and the local microbe transfer rate at the wall show the opposite trend. What is more, in the same table, it is indicated that the local skin friction coefficient, the local Sherwood number, the local Nusselt, and the microorganism transfer coefficient all are rapidly dropping for the increment of the Casson parameter.Table 4
Values forCfxRex1/2, NuxRex−1/2, ShxRex−1/2, and NnxRex−1/2 for various values of M,λ,α, β,ε1,Ec,Nt,Nb,ε2, and ε3 with Pr=2.0,Sc=2.0,Pe=0.6,σ=0.2, and Sb=0.7.
Mλαβε1EcNtNbε2ε3CfxRex1/2NuxRex−1/2ShxRex−1/2NnxRex−1/20.00.20.20.70.20.20.20.50.20.20.9647560.4063880.6979430.6574971.00.20.20.70.20.20.20.50.20.21.3170400.3013940.6105250.5619282.00.20.20.70.20.20.20.50.20.21.5520210.2275610.5468690.4964671.00.00.20.70.20.20.20.50.20.21.9019810.2674030.7510620.6810021.00.20.20.70.20.20.20.50.20.21.3170400.3013940.6105250.5619281.00.30.20.70.20.20.20.50.20.21.1481060.3128810.5647970.5224321.00.20.00.70.20.20.20.50.20.21.4479310.2937120.6067530.5582941.00.21.00.70.20.20.20.50.20.20.8928760.3212130.6025380.5549281.00.21.50.70.20.20.20.50.20.20.7016440.3226710.5760410.5304531.00.20.20.50.20.20.20.50.20.21.4137510.3054970.6202850.5726981.00.20.21.00.20.20.20.50.20.21.2317900.2966910.5990640.5499371.00.20.22.00.20.20.20.50.20.21.1120510.2879030.5776870.5286251.00.20.20.70.00.20.20.50.20.21.3150710.2589270.6118150.5630611.00.20.20.70.20.20.20.50.20.21.3170400.3013940.6105250.5619281.00.20.20.70.50.20.20.50.20.21.3195850.3592160.6092840.5606461.00.20.20.70.20.00.20.50.20.21.3122110.4331940.5817580.5508271.00.20.20.70.20.20.20.50.20.21.3170400.3013940.6105250.5619281.00.20.20.70.20.50.20.50.20.21.3243240.1031640.6537990.5782841.00.20.20.70.20.20.10.50.20.21.3154120.3257050.6154210.5629651.00.20.20.70.20.20.20.50.20.21.3170400.3013940.6105250.5619281.00.20.20.70.20.20.50.50.20.21.3215110.2384650.6030820.5653591.00.20.20.70.20.20.20.20.20.21.3093820.4563790.5087590.5272581.00.20.20.70.20.20.20.50.20.21.3170400.3013940.6105250.5619281.00.20.20.70.20.20.21.00.20.21.3276240.1254730.6307760.5668841.00.20.20.70.20.20.20.50.00.21.3159920.3149110.7494590.6170731.00.20.20.70.20.20.20.50.20.21.3170400.3013940.6105250.5619281.00.20.20.70.20.20.20.50.50.21.3184420.2840450.4807110.511491.00.20.20.70.20.20.20.50.20.01.3170400.3013940.6105250.6391921.00.20.20.70.20.20.20.50.20.21.3170400.3013940.6105250.5619281.00.20.20.70.20.20.20.50.20.51.3170400.3013940.6105250.480001


## 6. Conclusions

The nondimensionality governing equations of MHD non-Newtonian Casson dissipative bionanofluid flow saturated by gyrotactic microorganisms with variable transport characteristics and slip velocity phenomenon are solved numerically using the Galerkin technique and shifted Vieta-Lucas polynomials. Under the controller parameters, the temperature function, velocity function, concentration function, and microbe concentration function are evaluated using graphical representations. The efficiency and accuracy of the results were ascertained by comparing the approximate solutions with the solution obtained from the special case of the problem or by evaluating the REF. The numerical simulation work is performed with the aid of the Mathematica software package.The following are the most crucial decisive points based on the numerical results:(1)
As the magnetic parameter, slip velocity parameter, viscosity parameter, and Casson parameter are increased, the dimensionless microorganisms increase; however, as the Brownian motion parameter is increased, the opposite trend is observed.(2)
When the magnetic parameter, slip velocity parameter, viscosity parameter, and Eckert number are increased, the temperature function displays an increasing trend.(3)
For increasing levels of the magnetic number, slip velocity parameter, Casson parameter, and viscosity parameter, the nanofluid’s velocity collapses.(4)
The skin friction factor shrinks as the viscosity, slip velocity, and Casson parameters swell, whereas the magnetic parameter swells in the opposite direction.(5)
Nusselt number values rise with cumulative increases in thermal conductivity quantity, slip, and viscosity parameters but fall with increases in magnetic factor and Casson parameter.



---

*Source: 5236196-2022-08-27.xml*