Skip to content
BY 4.0 license Open Access Published by De Gruyter September 24, 2021

Gyrotactic microorganism and bio-convection during flow of Prandtl-Eyring nanomaterial

  • Tasawar Hayat , Inayat Ullah , Khursheed Muhammad EMAIL logo and Ahmed Alsaedi
From the journal Nonlinear Engineering

Abstract

Our main intension behind this work is to investigate Prandtl-Eyring nanomaterial in presence of gyrotactic microorganisms. Flow is generated via stretching sheet and is subject to melting heat effect. Radiation and dissipation are addressed. Entropy rate is also reported. Nanofluid effects are explored through Buongiorno model for nanofluid by considering Brownian motion and thermophoresis impacts. Problem related modelling is done by obtaining PDEs and these PDEs are then transmitted into ODEs by using appropriate similarity variables. Homotopic technique has been employed to obtain a convergent series solution of the considered problem. Graphical results have been presented to investigate the impact of different prominent variables over fluid velocity, temperature distribution, nanofluid concentration and on microorganism concentration. Entropy analysis has been discussed and the physical quantities such as surface drag force, Nusselt number, local Sherwood number and microorganism density number for the current problem is obtained. Velocity boost against higher melting and fluid parameters. Temperature of the fluid reduces with an increment in melting and radiation parameters while it intensifies through Prandtl and Eckert number, Brownian motion and thermophoresis parameters. Decay in concentration is noticed against higher values of melting and thermophoresis parameters while it increases for higher Schmidt number and Brownian motion parameter. Microorganism field boosts with higher values of Peclet number and microorganism concentration difference parameter. Moreover entropy generation rate intensifies against higher radiation parameter and Brickman number.

1 Introduction

Recently advancements in nanotechnology got much attention from scientists and researchers due to its enormous applications in industrial and technological zones. The introduction of nanofluid brings a lot of its applications to the real world. In the beginning Choi [1] observed that the heat transportation capacity (thermal conductivity) of normal fluids can be improved by adding small nanosize particles with the characteristics of high thermal conductivity. He presented the idea of replacing normal/traditional fluids with these more advance fluids with high thermal efficiency. These advance fluids are known as nanofluids. Nanofluid plays a vital role in technological devices, heat-related equipment, nuclear reactor, space technology, processing of generic drugs, microelectronics, radiators, cooling system, bio-sensors and many more. Buongiorno [2] explained the transportation of mass and heat transfer in nanofluids. He presented the seven slip mechanism and concluded that out of all these thermophoresis and Brownian diffusion is of particular interest. Ahmad et al. [3] investigated the nanomaterial flow of magnetohydrodynamic non-Newtonian fluid. Turkyilmazoglu [4] studied the mass and heat transport of unsteady natural convective flow of nanoliquid with radiation past a vertical surface. MHD forced convection in nanofluid is reported to Sheikholeslami and bhatti [5]. Hayat et al. [6] considered melting effect in CNTs flow with stagnation point by variable thickness surface. Some recent investigations on nanofluid can be see through Refs. [7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17].

Bio-convection appears due to the up swimming of microorganisms in a particular direction since the density of microorganisms are higher than the base fluid. Therefore, the upper surface of the fluid becomes denser due to the gathering of microorganism. Thus upper surface of fluid is unstable and eventually the microorganism falls down causing bioconvection. The up swimming of microorganisms back to the top carry on bioconvection process. The utilization of gyrotactic microorganism in different nanoliquid enhances mass transfer, microscale mixing and also enhances the stability of nanofluids. Khan et al. [18] studied bioconvection nanofluid flow considering gyrotactic microorganism. Mutuku and Makinde [19] explored MHD bioconvection flow of nanoliquid past a porous vertical surface with gyrotactic microorganism. Hayat et al. [20] examined bioconvection in flow of Walter-B nanomaterials. Kuznetsov [21] examined the bioconvection nanofluid flow with the suspension of both gyrotactic microorganism and nanoparticles.

There is no doubt that the phenomenon of melting heat transfer has been studied extensively due to its ever rising usages in industrial and technological processes. Melting heat transfer plays an essential role in coil exchanger, magma solidification, semiconductor material preparation, welding process, crystal growth, permafrost melting, thermocouple and heat engines. Krishnamurthy et al. [22] considered chemical reactions and radiation during nanofluid flow with melting heat. Prasannakumara et al. [23] investigated melting heat, thermal radiation, MHD and heat source effects doing flow in region of stagnation point. Both melting heat transfer phenomenon and thermal radiation in magnetohydrodynamic flow past a moving surface is addressed by Das [24]. Further explorations on melting heat transfer phenomenon can be seen through Refs. [25, 26, 27, 28, 29].

Very little attention is paid by the researchers towards flow of non-Newtonian fluids in presence of microorganisms. This consideration narrow down when melting heat is considered. Here we have considered melting effect in Prandtl-Eyring nanofluid flow containing gyrotactic microorganism past a stretching sheet. Viscous dissipation and thermal radiation effects are also addressed. Entropy generation is accounted. Appropriate variables are used to transform governing expressions (PDEs) into ODEs. HAM technique is employed for solutions development. Entropy, velocity, concentration, microorganism and temperature under involved physical variables is evaluated graphically.

2 Mathematical formulation

Two dimensional Prandtl-Eyring nanofluid with motile gyrotactic microorganisms is considered by a stretched sheet. Sheet is subject to melting phenomena. Thermal radiation and viscous dissipation elaborates heat transport phenomena. Surface is stretched with the velocity (u = ax), with constant (a). Flow geometry is given in Figure 1.

Figure 1 Flow diagram.
Figure 1

Flow diagram.

After executing aforementioned assumptions and boundary layer approximations we obtain [30, 31]:

(1) ux+vy=0,

(2) uux+vuy=AρfC*2uy2A2ρfC*3(uy)22uy2,

(3) uTx+vTy=kf(ρcp)f2Ty2+1(ρcp)f[AC*(uy)2A6C*3(uy)4]+16σ*T33k*(ρcp)f2Ty2+τ(DBCyTy+DTT(Ty)2),}

(4) uCx+vCy=DB(2Cy2)+DTT(2Ty2),

(5) uNx+vNy+bWc(CCm)[y(NCy)]=Dm2Ny2,

With boundary conditions

(6) u=ax,T=Tm,C=Cm,N=Nm,v=0,aty=0,u0,TT,CC,NNasy,k(Ty)y=0=ρnf[λ1+cs(TmT0)]v(x,0).}

Consider the transformations [6, 12]:

(7) u=axf(η),v=aνf(η),η=aνy,θ(η)=TTmTTm,ϕ(η)=CCmCCm,χ(η)=NNmNNm.}

Making use of these transformations in Eqs. (16) we get

(8) αf′″(η)αβf2(η)f′″(η)f2(η)+f(η)f(η)=0,

(9) θ(η)+Rdθ(η)+PrEc(αf2(η)αβf4(η))+Prf(η)θ(η)+PrNbθ(η)φ(η)+PrNtθ2(η)=0,

(10) ϕ(η)+Scf(η)ϕ(η)+NtNbθ(η)=0,

(11) χ(η)+Lbf(η)χ(η)Pe[ϕ(η)(χ(η)+Ω)+ϕ(η)χ(η)]=0,

with transformed boundary conditions

(12) f(0)=1,θ(0)=0,χ(0)=0,ϕ(0)=0,f()=0,ϕ()=1,θ()=1,χ()=1,Prf(0)+Meθ(0)=0.}

Related parameters are

(13) Rd=16σ*T33kk*,Nb=τDB(CCm)υ,α=AμC*,β=a3x26C*2ν,Nt=τDT(TTm)Tν,Pe=bWcDm,Pr=υα,Me=cp(TTm)λ1+cs(TmT0),Sc=υDB,Lb=υDm,Ec=(ax)2cp(TTm),Ω=NmNNm.}

3 Quantities of interest

Skin friction, Nusselt, Sherwood and Microorganism density numbers (Cfx, Nux, Shx, Nnx) are expressed as

(14) Cfx=τw12ρUw2,Nux=xqwk(TTm),Shx=xJwDB(CCm),Nnx=xJnDm(NNm).}

Here (τw,qw,Jw,Jn) are given by

(15) τw=(AC*(uy)A6C*3(uy)3)|y=0,qw=kf(1+16σ*T33kk*)(Ty)|y=0,Jw=DB(Cy)|y=0,Jn=Dm(Ny)|y=0.}

Making use in Eq. (13) we get

(16) 12Rex1/2Cfx=αf(0)13αβ[f(0)]3,Rex1/2Nux=(1+Rd)θ'(0),ShxRex1/2=ϕ(0),NnxRex1/2=χ(0).}

3.1 Entropy analysis

Entropy is given by [29]

(17) EG=kT2(1+16σ*3kk*T3)(Ty)2+RDC(Cy)2+RDN(Ny)2+RDT(TyCy)+1T(AC*(uy)2A6C*(uy)4)+RDT(NyTy).}

In non-dimensionalized form

(18) SG=α1(1+Rd)θ2(η)+L1α2α1ϕ2(η)+L2α3α1χ2(η)+L1θ(η)ϕ(η)+Br[αf2(η)αβf4(η)]+L2χ(η)θ(η).}

Bejan number is given as

(19) Bejannumber=EntropyduetomassandheattransferTotalentropyofthesystem

or

(20) Be=α1(1+Rd)θ2(η)+L1α2α1ϕ2(η)+L2α3α1χ2(η)+L1θ(η)ϕ(η)+L2χ(η)θ(η)α1(1+Rd)θ2(η)+L1α2α1ϕ2(η)+L2α3α1χ2(η)+L1θ(η)ϕ(η)+Br[αf2(η)αβf4(η)]+L2χ(η)θ(η)

Where

(21) L1=RD(CCm)k,Br=μ(ax)2k(TTm),L2=RD(NNm)k,SG=νEGkα1a,α1=TTmT,α2=CCmC,α3=NNmN.}

4 Solution procedure

Obtained non-linear ODEs with BCs are tackled analytically through homotopy analysis method (HAM). This technique is very helpful for solving highly non-linear ODEs. This method works on basis of defining appropriate initial guesses and linear operators. Thus initial guesses and linear operators are [32, 33, 34]:

(22) f0(η)=(1eη)MePr,θ0(η)=1eη,ϕ0(η)=1eη,χ0(η)=1eη,}

(23) Lf=3η3η,Lθ=2η21,Lφ=2η21,Lχ=2η21,}

which has the property

(24) Lf=[c0+c1eη+c2eη],Lθ=[c3eη+c4eη],Lφ=[c5eη+c6eη],Lχ=[c7eη+c8eη].}

Here ci (i = 1, 2, 3..., 8) are all positive constants. where (c1 = c3 = c5 = c7 = 0)

Zeroth order deformation equation is defined as

(25) Lf[F(η;p)f0(η)](1p)=pfNf[F(η;p)],

(26) Lθ[θ(η;p)θ0(η)](1p)=pθNθ[θ(η;p)],

(27) Lϕ[ϕ(η;p)ϕ0(η)](1p)=pϕNϕ[ϕ(η;p)],

(28) Lχ[χ(η;p)χ0(η)](1p)=pχNχ[χ(η;p)],

(29) F(0;p)η=1,θ(0;p)=0,χ(;p)=1,χ(0;p)=0,θ(;p)=1,ϕ(;p)=1,F(;p)η=0,ϕ(0;p)=0,PrF(0;p)+Meθ(0;p)η=0.}

Where p ∈ [0, 1] is embedding parameter while (ħf, ħθ, ħϕ, ħχ) are convergence control parameters.

The non-linear operators are written as

(30) Nf=α3F(η;p)η3(F(η;p)η)2+F(η;p)2F(η;p)η2αβ(2F(η;p)η2)23F(η;p)η3,

(31) Nθ=PrNb(θ(η;p)ηφ(η;p)η)+2θ(η;p)η2+PrNt(θ(η;p)η)2+Pr(F(η;p)θ(η;p)η)+Rd2θ(η;p)η2+αPrEc(2F(η;p)η2)2αβPrEc(2F(η;p)η2)4,}

(32) Nϕ=2ϕ(η;p)η2+NtNb2θ(η;p)η2+ScF(η;p)ϕ(η;p)η,

(33) Nχ=2χ(η;p)η2+LbF(η;p)χ(η;p)ηPe[2ϕ(η;p)η2(χ(η;p)+Ω)+ϕ(η;p)ηχ(η;p)η].

The solution in series form is written as

(34) f(η)=f0(η)+fm(η),θ(η)=θ0(η)+θm(η),ϕ(η)=ϕ0(η)+ϕm(η),χ(η)=χ0(η)+χm(η).}

5 Convergence analysis

Convergence of solution is highly dependent on auxiliary variables (ħf, ħθ, ħϕ, ħχ). These variables have great role is obtaining a convergent series solution. For this sketches of ħ− curves have been presented in Figure 2. Convergence regions are −1. 4 ≤ ħf ≤ −0.2, −1.5 ≤ ħθ ≤ −0.5, −1.5 ≤ ħϕ ≤ −0.1, and −1.5 ≤ ħχ ≤ −0.2.

Figure 2 Combine ħ− curves.
Figure 2

Combine ħ− curves.

6 Discussion

The investigations of different influential parameters over velocity field, temperature distribution, mass concentration, microorganism concentration and on entropy rate is the main aim behind this section

6.1 Velocity

Variations of melting and fluid parameters (Me, α) for velocity field (f′(η)) are depicted in Figures 34. Figure 3 displays outcomes of melting parameter for velocity. Velocity enhances against higher melting variable (Me). Physically with an increment in melting (Me) the heat transfer enhances due to high convection from hot fluid to cold melting surface and as a result velocity boosts. Figure 4 elucidate impact of material parameter (α) for velocity (f′(η)). Clearly an augmentation in velocity occurs against higher material parameter (α). It is quite clear that absolute viscosity and fluid parameter (α) are related inversely. Hence higher fluid parameter (α) leads to intensification in velocity of fluid (f′(η)).

Figure 3 f′(η) versus Me.
Figure 3

f′(η) versus Me.

Figure 4 f′(η) versus α.
Figure 4

f′(η) versus α.

6.2 Temperature field

Influences of different prominent variables on temperature are displayed in Figures 510). Figure 5 depicts melting parameter (Me) variations for temperature distribution. Higher melting parameter decays temperature distribution. Physically for higher melting (Me) convective flow enhances which cause more heat transfer from hot fluid towards melting surface and as a results the temperature decays. Figure 6 illustrates variation of Brownian motion parameter against temperature. Temperature rises for higher (Nb). Outcomes of radiation (Rd) for temperature is displayed in Figure 7. A reduction in temperature occurs against larger radiation parameter. Figure 8 elucidate temperature variations for thermophoresis parameter (Nt). An augmentation in thermophoresis variable results in temperature enhancement. Figure 9 presents Prandtl number (Pr) outcomes for temperature. Temperature boosts against higher Prandtl number. Decline in thermal diffusivity occurs for higher Prandtl number and temperature enhances against higher (Pr). Figure 10 depict Schmidt number (Sc) variations for temperature. Clearly temperature augmentation occurs for higher Schmidt number.

Figure 5 θ(η) versus Me.
Figure 5

θ(η) versus Me.

Figure 6 θ(η)versus Nb.
Figure 6

θ(η)versus Nb.

Figure 7 θ(η) versus Rd.
Figure 7

θ(η) versus Rd.

Figure 8 θ(η) versus Nt.
Figure 8

θ(η) versus Nt.

Figure 9 θ(η) versus Pr.
Figure 9

θ(η) versus Pr.

Figure 10 θ(η) versus Ec.
Figure 10

θ(η) versus Ec.

Table 1

Convergence of f″(0), θ′(0), ϕ′(0), χ′(0) when α = β = Pr = Pe = 0.5, Nt = Nb = Rd = 0.01, Me = Sc = 0.1, Lb = 0.2 at various approximations.

Order of approximation f″(0) θ′(0) ϕ′(0) χ′(0)

1 0.8125 0.6664 0.4833 0.667
5 0.8284 0.2883 0.2039 0.3045
10 0.8329 0.2063 0.1691 0.2161
15 0.8304 0.1839 0.1504 0.1909
20 0.8304 0.1839 0.1504 0.1803
22 0.8304 0.1839 0.1504 0.1803
Table 2

Nomenclature for the considered problem.

Expression Name Unit
Pe Peclet number Dimensionless
Sc Schmidt number Dimensionless
Nt Thermophoresis parameter Dimensionless
Wc Swimming cell speed m.s−1
T Temperature K
Br Brinkman number Dimensionless
Uw Stretching velocity m.s−1
AC* Fluid material coefficient kg.m−1.s−1
k Thermal conductivity s−3.kg.m.K−1
α1, α2 Temperature and concentration difference parameters Dimensionless
Me Melting parameter Dimensionless
Pr Prandtl number Dimensionless
Cfx Skin friction coefficient Dimensionless
cs Heat capacity m2.s−1.K−1
b Chemotaxis constant M
μf Dynamic viscosity kg.m−1.s−1
Rex Local Reynold number Dimensionless
Tm Reference temperature K
cp Specific heat m2.s−1.K−1
Shx Sherwood number Dimensionless
Nux Nusselt number Dimensionless
(ρcp)f Fluid heat capacity kg.m−1.s−2.K−1
ρf Density of fluid kg.m−3
Ec Eckert number Dimensionless
f Dimensionless velocity profile Dimensionless
DT Thermophoresis diffusion coefficient m2 s−1
C, N (Concentration, Microorganisms) fluid mol.m−3
ϕ Dimensionless concentration profile Dimensionless
ϕ′(0) (Concentration, Microorganisms) surface mol.m−3
a Dimensional constant (positive) s−1
θ Dimensionless temperature profile Dimensionless
(νf) Kinematic viscosity m2.s−1
u, v Components of velocity along x, y direction m.s−1
(ρcp)nf Nanofluid heat capacity kg.m−1.s−2.K−1
Dm Microorganism diffusion coefficient m2 s−1
Nb Brownian motion parameter Dimensionless
Nnx Microorganism density number Dimensionless
k* Mean absorption coefficient m−1
hf Coefficient of heat transfer W.m−2 k−1
SG Entropy rate Dimensionless
Ω Microorganism concentration difference parameter Dimensionless
Rd Radiation parameter Dimensionless
L1 Concentration diffusion parameter Dimensionless
L2 Microorganism diffusion parameter Dimensionless
DB Brownian motion coefficient m2 s−1
α, β Material parameters Dimensionless
C, N (Concentration, Microorganisms) ambient mol.m−3
τ Ratio of heat capacities Dimensionless
Tw, T Surface and ambient temperature K
χ Dimensionless microorganism profile Dimensionless
θ′(0) Stefan-Boltzmann constant kg.s−3.K−4
λ1 Latent heat kJ.kg−1
Lb Bioconvection Lewis number Dimensionless
α3 Microorganism difference parameter Dimensionless

6.3 Concentration

Variations of concentration for different variables are presented in Figures 1114. Figure 11 display outcome of melting parameter for concentration (ϕ(η)). An augmentation in melting parameter result in concentration decay. Figure 12 depict variations of Brownian motion parameter on concentration. Concentration boosts against higher (Nb). Physically molecular diffusivity decreases with augmentation in Brownian motion parameter which results in concentration enhancement. Figure 13 display thermophoresis outcome on concentration (ϕ(η)). A decline in concentration is observed for larger values of thermophoresis (Nt). Figure 14 elucidate outcome of Schmidt number for concentration (ϕ(η)). With an increment in Schmidt number concentration boosts. Basically mass transfer decays form hot fluid to cold surface which is responsible for concentration enhancement.

Figure 11 ϕ(η) versus Me.
Figure 11

ϕ(η) versus Me.

Figure 12 ϕ(η) versus Nb.
Figure 12

ϕ(η) versus Nb.

Figure 13 ϕ(η) versus Nt.
Figure 13

ϕ(η) versus Nt.

Figure 14 ϕ(η) versus Sc.
Figure 14

ϕ(η) versus Sc.

6.4 Microorganism

Influences of bio convection lewis number (Lb), microorganism concentration difference (Ω), Peclet number (Pe) and melting parameter (Me) on microorganism field are displayed in Figures 1518. Figure 15 displays impact of bio convection Lewis number on microorganism field. A decline in microorganism field occurs for higher values bio convection Lewis number. Influence of Peclet number on microorganism field is plotted in Figure 16. Higher Peclet number results in microorganism field enhancement. Figure 17 illustrates the impact of (Ω) on microorganism field. An augmentation in microorganism field occurs against higher (Ω). Figure 18 depict outcomes of melting parameter for microorganism field. Microorganism field reduces against higher melting variable (Me).

Figure 15 xi(η) versus Lb.
Figure 15

xi(η) versus Lb.

Figure 16 xi(η) versus Pe.
Figure 16

xi(η) versus Pe.

Figure 17 xi(η) versus Ω.
Figure 17

xi(η) versus Ω.

Figure 18 χ(η) versus Me.
Figure 18

χ(η) versus Me.

6.5 Entropy rate

Variations of radiation (Rd) and Brinkman number (Br) on entropy rate are presented in Figures 1920. Figure 19 depicts outcomes of Brinkman number for entropy rate (SG(η)). An augmentation in entropy rate occurs against higher Brinkman number. Since viscous forces increases for higher Brinkman number which causes disturbance in the system and consequently entropy rate boosts against Brinkman number. Figure 20 shows outcome of radiation for entropy rate (SG(η)). Clearly entropy rate boosts against higher radiation variable (Rd).

Figure 19 SG(η) versus Br.
Figure 19

SG(η) versus Br.

Figure 20 SG(η) versus Rd.
Figure 20

SG(η) versus Rd.

7 Conclusion

Presented work deals with flow of Prandtl fluid in the presence of microorganisms. Surface is subjected to melting heat. Boungiorno model for nanofluid describes nanofluid characteristics while heat transfer characteristics are explored via thermal radiation and viscous dissipation. It is concluded that velocity increases with an increase in fluid and melting parameters. Decline in temperature of the fluid is noticed for both radiation and melting parameters while it intensifies for larger Brownian motion and thermophoresis parameters. Higher melting parameter causes reduction in concentration while it increases with higher Brownian motion parameter and Schmidt number. Higher Peclet number leads to intensification in microorganism field while it decays against bio-concentration Lewis number. Entropy production rate is higher for Brinkman number and radiation parameter. In future this work can be extended for cavity problems and can be handled by finite element method.

  1. Funding information: The authors state no funding involved.

  2. Author contributions: All authors have accepted responsibility for the entire content of this manuscript and approved its submission.

  3. Conflict of interest: The authors state no conflict of interest.

References

[1] Choi SU. Enhancing thermal conductivity of fluids with nanoparticles. ASME-Publications-Fed. 1995;231:99–106.Search in Google Scholar

[2] Buongiorno J. Convective transport in nanofluids. J Heat Transfer. 2006;128(3):240–50.10.1115/1.2150834Search in Google Scholar

[3] Ahmad S, Khan MI, Hayat T, Alsaedi A. Inspection of Coriolis and Lorentz forces in nanomaterial flow of non-Newtonian fluid with activation energy. Phys A Statt Mech App. 2020;54015:123057. https://doi.org/10.1016/j.physa.2019.123057.10.1016/j.physa.2019.123057Search in Google Scholar

[4] Turkyilmazoglu M, Pop I. Heat and mass transfer of unsteady natural convection flow of some nanofluids past a vertical infinite flat plate with radiation effect. Int J Heat Mass Transf. 2013;59:167–71.10.1016/j.ijheatmasstransfer.2012.12.009Search in Google Scholar

[5] Sheikholeslami M, Bhatti MM. Forced convection of nanofluid in presence of constant magnetic field considering shape effects of nanoparticles. Int J Heat Mass Transf. 2017;111:1039–49.10.1016/j.ijheatmasstransfer.2017.04.070Search in Google Scholar

[6] Hayat T, Muhammad K, Farooq M, Alsaedi A. Melting heat transfer in stagnation point flow of carbon nanotubes towards variable thickness surface. AIP Adv. 2016;6(1):015214.10.1063/1.4940932Search in Google Scholar

[7] Abbas N, Nadeem S, Saleem S, Issakhov A. Transportation of modified nanofluid flow with time dependent viscosity over a Riga plate: exponentially stretching. Ain Shams Eng J. 2018;8:552–8.10.1016/j.asej.2021.01.034Search in Google Scholar

[8] Alharbi SO. Influence of wall slip and jump in wall temperature on transport of heat energy in hybrid nanofluid. J Therm Anal Calorim. 2021;144:847–54.10.1007/s10973-020-09428-xSearch in Google Scholar

[9] Rabbi KM, Sheikholeslami M, Karim A, Shafee A, Li Z, Tlili I. Prediction of MHD flow and entropy generation by artificial neural network in square cavity with heater-sink for nanomaterial. Physica A. 2020;541:123520.10.1016/j.physa.2019.123520Search in Google Scholar

[10] Khan M, Ahmad L, Yasir M, Ahmed J. Numerical analysis in thermally radiative stagnation point flow of Cross nanofluid due to shrinking surface: dual solutions. Appl Nanosci. 2021. https://doi.org/10.1007/s13204-021-01861-0.10.1007/s13204-021-01861-0Search in Google Scholar

[11] Ghahremanian A, Abbassi A, Mansoori Z, Toghraie D. Investigation the nanofluid flow through a nanochannel to study the effect of nanoparticles on the condensation phenomena. J Mol Liq. 2020;3111:113310.10.1016/j.molliq.2020.113310Search in Google Scholar

[12] Hayat T, Khan SA, Ijaz Khan M, Alsaedi A. Theoretical investigation of Ree-Eyring nanofluid flow with entropy optimization and Arrhenius activation energy between two rotating disks. Comput Methods Programs Biomed. 2019 Aug;177:57–68.10.1016/j.cmpb.2019.05.012Search in Google Scholar PubMed

[13] Hayat T, Shah F, Khan MI, Alsaedi A. Numerical simulation for aspects of homogeneous and heterogeneous reactions in forced convection flow of nanofluid. Results Phys. 2018;8:206–12.10.1016/j.rinp.2017.11.035Search in Google Scholar

[14] Rasool G, Shafiq A, Alqarni MS, Wakif A, Khan I, Bhutta MS. Numerical Scrutinization of Darcy-Forchheimer Relation in Convective Magnetohydrodynamic Nanofluid Flow Bounded by Nonlinear Stretching Surface in the Perspective of Heat and Mass Transfer. micromach. 2021;12: https://doi.org/10.3390/mi12040374.10.3390/mi12040374Search in Google Scholar PubMed PubMed Central

[15] Khan WA, Rashad AM, Abdou MM, Tlili I. Natural bioconvection flow of a nanofluid containing gyrotactic microorganisms about a truncated cone. Eur J Mech BFluids. 2019;75:133–42.10.1016/j.euromechflu.2019.01.002Search in Google Scholar

[16] Mutuku WN, Makinde OD. Hydromagnetic bioconvection of nanofluid over a permeable vertical plate due to gyrotactic microorganisms. Comput Fluids. 2014;95(22):88–97.10.1016/j.compfluid.2014.02.026Search in Google Scholar

[17] Khan MI, Waqas M, Hayat T, Khan MI. Alsaedi, Behavior of stratification phenomenon in flow of Maxwell nanomaterial with motile gyrotactic microorganisms in the presence of magnetic field. Int J Mech Sci. 2017;131–132:426–34.10.1016/j.ijmecsci.2017.07.009Search in Google Scholar

[18] Kuznetsov AV. The onset of nanofluid bioconvection in a suspension containing both nanoparticles and gyrotactic microorganisms. Int Commun Heat Mass Transf. 2010;37(10):1421–5.10.1016/j.icheatmasstransfer.2010.08.015Search in Google Scholar

[19] Farooq S, Hayat T, Alsaedi A, Ahmad B. Numerically framing the features of second order velocity slip in mixed convective flow of Sisko nanomaterial considering gyrotactic microorganisms. Int J Heat Mass Transf. 2017;112:521–32.10.1016/j.ijheatmasstransfer.2017.05.005Search in Google Scholar

[20] Krishnamurthy MR, Gireesha BJ, Prasannakumara BC, Gorla RS. Thermal radiation and chemical reaction effects on boundary layer slip flow and melting heat transfer of nanofluid induced by a nonlinear stretching sheet. Nonlinear Eng. 2016;5(3):147–59.10.1515/nleng-2016-0013Search in Google Scholar

[21] Prasannakumara BC, Gireesha BJ, Manjunatha PT. Melting phenomenon in MHD stagnation point flow of dusty fluid over a stretching sheet in the presence of thermal radiation and non-uniform heat source/sink. Int J Comput Methods Eng Sci Mech. 2015;16(5):265–74.10.1080/15502287.2015.1047056Search in Google Scholar

[22] Das K. Radiation and melting effects on MHD boundary layer flow over a moving surface. Ain Shams Eng J. 2014;5(4):1207–14.10.1016/j.asej.2014.04.008Search in Google Scholar

[23] Siddiqui BK, Batool S, Hassan QM, Malik MY. Irreversibility analysis in the boundary layer MHD two dimensional flow of Maxwell nanofluid over a melting surface. Ain Shams Eng J. 2021;3: https://doi.org/10.1016/j.asej.2021.01.017.10.1016/j.asej.2021.01.017Search in Google Scholar

[24] Qi R, Wang Z, Ren J, Wu Y. Numerical investigation on heat transfer characteristics during melting of lauric acid in a slender rectangular cavity with flow boundary condition. Int J Heat Mass Transf. 2020;157:119927.10.1016/j.ijheatmasstransfer.2020.119927Search in Google Scholar

[25] Hayat T, Shah F, Alsaedi A, Ahmad B. Entropy optimized dissipative flow of effective Prandtl number with melting heat transport and Joule heating. Int Commun Heat Mass Transf. 2020;111:104454.10.1016/j.icheatmasstransfer.2019.104454Search in Google Scholar

[26] Duan J, Li F. Transient heat transfer analysis of phase change material melting in metal foam by experimental study and artificial neural network. J Energy Storage. 2021;33:102160.10.1016/j.est.2020.102160Search in Google Scholar

[27] Hayat T, Inayatullah, Alsaedi A. Inayatullah, Alsaedi A. Development of bioconvection flow of nanomaterial with melting effects. Chaos Solitons Fractals. 2021;148:111015.10.1016/j.chaos.2021.111015Search in Google Scholar

[28] Rehman KU, Malik AA, Malik MY, Tahir M, Zehra I. On new scaling group of transformation for Prandtl-Eyring fluid model with both heat and mass transfer. Results Phys. 2018;8:552–8.10.1016/j.rinp.2017.12.071Search in Google Scholar

[29] Samal SK, Moharana MK. Second law analysis of recharging microchannel using entropy generation minimization method. Int J Mech Sci. 2021;193(1):106174.10.1016/j.ijmecsci.2020.106174Search in Google Scholar

[30] Bilal M, Sagheer M, Hussain S. Numerical study of magnetohydrodynamics and thermal radiation on Williamson nanofluid flow over a stretching cylinder with variable thermal conductivity. Alex Eng J. 2018;57(4):3281–9.10.1016/j.aej.2017.12.006Search in Google Scholar

[31] Farooq U, Zhao YL, Hayat T, Liao SJ. Application of the HAMbased mathematica package BVPh 2.0 on MHD Falkner-Skan flow of nanofluid. Comput Fluids. 2015;111:64–75.10.1016/j.compfluid.2015.01.005Search in Google Scholar

[32] Khan NA, Khan H. A Boundary layer flows of non-Newtonian Williamson fluid. Nonlinear Eng. 2014;3(2):107–15.10.1515/nleng-2014-0002Search in Google Scholar

[33] Ali N, Khan SU, Abbas Z. Unsteady flow of third grade fluid over an osscillatory stretching sheet with thermal radiation and heat source/sink. Nonlinear Eng. 2015;4(4):223–36.10.1515/nleng-2015-0019Search in Google Scholar

[34] Naveed M, Abbas Z, Sajid M. Hydromagnetic flow over an unsteady curved stretching surface. Eng Sci Technol Int J. 2019;19:841–45.10.1016/j.jestch.2015.11.009Search in Google Scholar

Received: 2021-03-21
Accepted: 2021-06-17
Published Online: 2021-09-24

© 2021 Tasawar Hayat et al., published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Downloaded on 29.5.2024 from https://www.degruyter.com/document/doi/10.1515/nleng-2021-0015/html
Scroll to top button