# Bioconvection Flow of MHD Viscous Nanofluid in the Presence of Chemical Reaction and Activation Energy

Authors: Muhammad Imran Asjad; Muhammad Zahid; Fahd Jarad; Abdullah M. Alsharif
Journal: Mathematical Problems in Engineering
 (2022)

---

## Abstract

Enhancement of heat transfer due to stretching sheets can be appropriately controlled by the movement of the nanofluids. The concentration and settling of the nanoparticles in the viscous MHD fluid and bioconvection are addressed. In this scenario, the fluid flow occurring in the presence of a normal and uniform magnetic field, thermal radiation, and chemical reaction is taken into account. For the two-dimensional flow with heat and mass transfer, five dependent variables and three independent variables constitute the system of partial differential equations. For this purpose, similarity functions are involved to convert these equations to corresponding ODEs. Then, the Runge–Kutta method with shooting technique is used to evaluate the required findings with the utilization of MATLAB script. The fluid velocity becomes slow against the strength of the magnetic parameter. The temperature rises with the parameter of Brownian motion and thermophoresis. The bioconvection Lewis number diminishes the velocity field. Compared with the existing literature, the results show satisfactory congruences.


---

## Body

## 1. Introduction

The convoluted and quick process in massive machinery and little gadgets has produced a significant problem of thermal imbalance. Varied extraneous techniques like fins and fans are used; however, their utility is restricted because of their giant size. In 1995, the scientist Choi and Eastman [1] introduced that the nanosized particles mixed in the fluid called nanofluid have more capacity of heat transfer as compared with fluid without nanosized particles. Das et al. [2] explained the recent and future applications of fluid involving nanosized particles. Khan et al. [3] using the shooting method analyzed flow features of Williamson nanofluid influenced by variable viscosity depending on temperature and Lorentz force past an inclined nonlinear extending surface. Koo and Kleinstreuer [4] described influences of convection, conduction, viscous dissipation, and thermal transportation on nanofluid flow in a microchannel. Sui et al. [5] introduced the CattaneoeChristov model with double diffusion to analyze the influence of slip velocity, Brownian motion, and variable viscosity on the transportation of an upper convected Maxwell nanofluid through stretching sheet. Imran et al. [6] determined an unsteady stream of Maxwell fluid through an accelerated exponentially vertical surface with influences of radiation, Newtonian heating, MHD, and slip condition taken into account. Khan et al. [7] investigated the flow of micropolar base nanofluid through stretching sheet with thermal radiation and magnetic dipole. Sheikholeslami and Rokni [8] scrutinized magnetic field impacts on the thermal transport rate in a nanofluid. Seyyedi et al. [9] analyzed the entropy generation for Cu-water nanofluid having a semi-annulus porous wavy cavity in the presence of a magnetic field. Molana et al. [10] discussed the characteristics of heat transfer and natural convection of nanofluid past a porous cavity with a constant inclined magnetic field. Dogonchi et al. [11] explained the characteristics of natural convection and magnetic field on nanofluid flow through porous medium with effects of Hartmann number, Rayleigh number, and Darcy number taken into account. Shaw et al. [12] scrutinized the impact of nonlinear thermal and entropy generation on Casson nanofluid flow with rotating disk and also described the brain function. Chamkha et al. [13] explained MHD nanofluid flow through cavity using the control-volume-based finite element method with effects of natural convection, thermal radiation, and shape factor of nanoparticles taken into account. Dogonchi et al. [14] numerically introduced the importance of the Cattaneo–Christov theory of heat conduction through triangular semicircular heater with viscosity dependent on the magnetic field. Seyyedi et al. [15] described the entropy generation and natural convection heat transfer of Cu-water nanoliquid through the hexagonal cavity. Sadeghi et al. [16] analyzed the thermal behavior of magnetic buoyancy-driven flow in ferrofluid-filled wavy enclosure furnished with two circular cylinders.Stretching sheets because of their wide applications like plastic film, metal drawing and spinning, glass fiber, paper processing, and heat moving have become an important topic in the past decades. Recently, some researchers investigated the magnetohydrodynamic flow with the different effects such as viscous dissipation and chemical reaction using stretching sheet (see Ismail et al. [17], Rajput et al. [18], Swain et al. [19], Reddy et al. [20], Jat and Chand [21], Sajid and Hayat [22], Ishak [23], Abd El-Aziz [24], Makinde [25], and Goud et al. [26]).Bioconvection described the phenomena in which living microorganisms denser than water swim upward in suspensions. These microorganisms pile up in the layer of the upper surface and lower surface becomes less dense. The microorganisims fall down due to unstability of density distribution. Bioconvection has applications in biological systems and biotechnology such as purifying cultures, enzyme biosensors, and separating dead and living cells [27]. Raees et al. [28] examined the unsteady stream of bioconvection mixed nanofluid having gyrostatic motile microorganisms through a horizontal channel. Siddiqa et al. [29] numerically studied the bioconvection flow of nanofluid having mass and thermal transportation along with gyrotactic microorganisms through a curved vertical cone. Abbasi et al. [30] introduced the bioconvection stream of viscoelastic nanofluid because of gyrotactic microorganisms past a rotating extending disk having zero mass flux and convective boundary condition and also described the relatable parameters influences on velocity, temperature, local density, Sherwood number, and Nusselt number in detail. Chu et al. [31] analyzed the stream of bioconvection MHD fluid through extending sheet with the significance of motile microorganisms, activation energy, thermophoresis diffusion, Brownian motion, and chemical reaction taken into account. Henda et al. [32] examined the magnetized bioconvection flow of fluid past an extending cylinder with thermal radiation, activation energy, and heat source. Khan et al. [33] scrutinized bioconvection stream of viscous nanofluid through (cone, wedge, and plate) multiple geometries with effects of heat flux, cross-diffusion, and Cattaneo–Christov.Inspired by the above literature survey, our interest pertains to extending the results of Goud et al. [26] to investigate a more general problem, including bioconvection of nanofluid transportation with the effects of chemical reaction and radiation to avoid probable settling of nanoentities. The connotation of such meaningful attributes can be a useful extension, and the results can be utilized for desired effective thermal transportation in the heat exchanger of various technological processes.


## 2. Problem Formulation

Here, we considered steady incompressible magnetohydrodynamic nanofluid flow through exponentially stretching sheet along thex-axis and y-axis taken to be normal with velocity U˜w=a0ex/l as shown in Figure 1. A magnetic field is applied to the flow region and acts in the y-direction. A mild diffusion of microorganisms and nanoparticles is set in the fluid. Thermal radiation is considered, and bioconvection takes place because of microorganisms’ movement. The fluid velocity for two-dimensional flow is u˜,v˜.Figure 1
Geometry of the problem.Under the above conditions, the governing equations are as follows [20, 26]. Continuity equation is as follows:(1)u˜x+v˜y=0,momentum equation is as follows:(2)u˜u˜x+v˜u˜y=νu˜yy−σρB02u˜+1ρgβρ1−C˜∞T˜−T˜∞−gρp−ρfC˜−C˜∞−gγρm−ρfN˜−N˜∞,energy equation is as follows:(3)u˜T˜x+v˜T˜y=αT˜yy−1ρCp∂qr∂y+τDBT˜yC˜y+DTT˜∞T˜y2,concentration equation is as follows:(4)u˜C˜x+v˜C˜y=DC˜yy−KrC˜−C˜∞−Kr2C˜−C˜∞T˜T˜∞nexp−EakT˜+DTT˜∞T˜yy,bioconvection equation is as follows:(5)u˜N˜x+v˜N˜y+dWc∂∂yN˜ΔC˜C˜y=N˜yyDn,with constraints(6)u˜=U˜w,v˜=0,T˜=T˜w,C˜=C˜w,N˜=N˜waty⟶0,u˜⟶0,T˜⟶T˜∞,C˜⟶C˜∞,N˜⟶N˜∞asy⟶∞.Now, introducing(7)U˜w=a0ex/l,T˜w=T˜∞+T˜0ex/2l,C˜w=C˜∞+C˜0ex/2l,N˜w=N˜∞+N˜0ex/2l,under the Rosseland approximation qr [26], equation (3) can be written as(8)u˜T˜x+v˜T˜y=T˜yyα+16σ∗T∞3ρCp3k1+τDBT˜yC˜y+DTT˜∞T˜y2.Introducing similarity transformation,(9)η=ya02νℓex/2l,u˜=a0ex/lf′η,v˜=−a0ν2ℓex/2lfη+ηf′η,ψ=2νLa0fηex/2l,T˜=T˜∞+T˜0ex/2l,C˜=C˜∞+C˜0ex/2l,N˜=N˜∞+N˜0ex/2l.In view of the above appropriate relations, equation (1) is satisfied identically and equations (2)–(5), respectively, become(10)f″′−Mf′−2f′2+ff″+λθ−Nrϕ−Rbχ=0,(11)1+43Kθ″+Prfθ′−Prθf′+θ′Nbϕ′+Ntθ′=0,(12)ϕ″+fϕ′−Crϕ−ϕf′−σmϕ1+δθnexp−E1+δθSc+NtNbθ″=0,(13)χ″ξ+LbPrfξχ′ξ−LbPrf′ξχξ−Peσ1ϕ″ξ+χξϕ″ξ+χ′ξϕ′ξ=0,and the constraints reduce to(14)f′0=1,f0=0,ϕ0=1,θ0=1,χ0=1,atη=0,(15)f′∞⟶0,ϕ∞⟶0,θ∞⟶0,χ∞⟶0.asη⟶∞.The associated parameters are(16)M=2σB02ℓρU˜w,Pr=να,λ=1−C∞βgTw−T∞2ℓU˜w2,Nr=ρp−ρfCw−C∞β1−C∞ρTw−T∞,Rb=ρm−ρfγNw−N∞ρ1−C∞βTw−T∞,Nt=τDTTw−T∞νT∞,Nb=τDBCw−C∞ν,σm=2Kr2ℓU˜w,δ=Tw−T∞T∞,K=4σ∗T∞3k∗K,σ1=N∞Nw−N∞,E=EakT∞,Sc=νD,Lb=αDn,Pe=dWcDn,Cr=2ℓKrU˜w,where M is the magnetic field parameter, Pr is the Prandtl number, λ is the mixed convection parameter, Nr is the buoyancy ratio number, Rb is the bioconvection Rayleigh number, Nt is the thermopherosis diffusion factor, Nb is the Brownian factor, σm denotes the dimensionless reaction rate, δ is used as the temperature distinction parameter, K is the radiation parameter, σ1 is the bioconvective difference parameter, E means the nondimensional energy activation, Sc is the Schmidt number, Lb is the bioconvection Lewis number, Pe is the peclet number, and Cr is the chemical reaction parameter.The wall shear stress, thermal flux, and mass flux, respectively, are given as(17)τw−μ∂u˜∂yy=0=0,qw+k∂T˜∂yy=0=0,jw+D∂c˜∂yy=0=0.Cf (skin friction), Nux (Nusselt number), and Shx (sherwood number) in dimensionless form are(18)Cf=f′′02Rex,Nux=−Rexθ′0,ShxRe−1/2=−ϕ′0.


## 3. Results and Discussion

Physical meanings of the final nondimensional formulation of time-independent MHD flow of nanofluid due to stretch of an exponential sheet in the presence of chemical reaction along the boundary constraints are solved numerically. Table1 contains results for −θ′0 (Nusselt number). Comparison of the results indicates acceptable agreement to validate this numeric procedure. In Figure 2, the velocity of the flow seems to be reduced significantly when magnetic parameter M0.0≤M≤2.5 is increased because high values of magnetic field parameter improve the contradictory force known as Lorentz force. The improvement of mixed convection parameter λ causes to boost the flow velocity f′η as shown in Figure 3. From Figures 4 and 5, significant rising behavior of θη is noticed with an enhanced value of Brownian motion parameter Nb and thermophoresis parameter Nt. The fast random motion of nanoparticles characterized by larger Nb is responsible for enhanced heat transfer to raise θη. Similarly, the higher Nt means a greater thermophoretic effect which moves the nanoparticles hotter regime to the colder one and increases the thermal distribution. The similarly larger value of E provides strength to ϕη as depicted in Figure 6. Figure 7 displays the decrement in ϕη due to the larger value of chemical reaction, and the chemical reaction becomes faster to recede nanoparticles concentration. The bioconvection Rayleigh Rb and parameter are responsible for given direct increment to χη as demonstrated in Figure 8.Table 1
Comparison ofθ′0 with changed values of K, M, and Pr.
KMPrIshak [23]Goud et al. [26]Bidin and Nazar [34]Our results0.001.00.954 80.954 7840.954 70.954 810 62.01.471 51.471 4621.471 41.471 454 03.01.869 11.869 0731.869 11.869 068 85.02.500 12.500 1112.500 128 010.03.660 43.660 3463.660 369 31.01.00.861 10.861 0970.861 508 61.00.00.531 20.531 1700.531 311 21.00.450 50.450 68700.450 695 5Figure 2
Influences ofM on f′.Figure 3
Influences ofλ on f′.Figure 4
Influences of Nb onθ.Figure 5
Influence of Nt onϕ.Figure 6
Influences ofE on ϕ.Figure 7
Influences of Cr onϕ.Figure 8
Influences of Rb onχ.


## 4. Conclusions

Theoretical and numeric analysis for magnetohydrodynamic of nanofluid owing to sudden stretched in an exponential sheet has been made in this communication. Effects of the emerging parameters are enumerated on the physical field, namely, velocity, temperature, and microorganisms distribution. Significant outcomes are summarized as follows:(i)
The velocity reduces withM and boosts with λ(ii)
The conclusion of nanoparticles characterized by parameters Nb and Nt shows an increment in the temperature profile(iii)
Concentration recurses with Cr and is enhanced withE(iv)
Bioconvection parameter is increased with Rb



---

*Source: 1707894-2022-02-07.xml*