Computer Modeling in Engineering & Sciences

DOI: 10.32604/cmes.2022.017391

ARTICLE

Bioconvection Cross Diffusion Effects on MHD Flow of Nanofluids over Three Different Geometries with Melting

Tahir Kamran, Muhammad Imran*, Muhammad N. Naeem and Mohsan Raza

Department of Mathematics, Government College University Faisalabad, Faisalabad, 38000, Pakistan
*Corresponding Author: Muhammad Imran. Email: drmimranranchaudhry@gmail.com
Received: 07 May 2021; Accepted: 09 June 2021

1  Introduction

Magnetohydrodynamics (MHD) has a vital role in several flow phenomena, including industry. The MHD is applied in various fields, including multidisciplinary technological areas like biochemical manufacturing heat transfer, automotive sector, ceramic technology, aerodynamic performance, metallurgical technology, mental operating techniques, and fluid dynamics. Magnetohydrodynamic initiation changes the flow field’s preferred direction by fluctuating boundary layer configuration flow. Electromagnetic fields play an important role in convection mechanisms like metallic casting, material processing fields, and nuclear reactor control plants. Zheng et al. [1] investigated the heat transformation in the two-dimensional flow of magnetohydrodynamic (MHD) by considering temperature and velocity slip effects through a shrinking porous surface. It was stated that the thermal boundary layer flow declined to owe to an increase in the shrinking parameters. The heat transmission and magnetohydrodynamic flow through the porous plate and shrinking surface were studied by Chauhan et al. [2]. The central concept was that the cooling rate was controlled by the electromagnetic number and suction parameter over channels bounded. Attia [3] examined the different physiological properties of the unsteady MHD Couette flowing heat exchange in a dusty fluid flow with variable viscosity. Kumar et al. [4] explored mass and heat transport across a vertical plate under the cross-diffusion impact in radiative mixed convective cases of magnetohydrodynamic fluid flow. The fluid movement increased with an increase in both the buoyancy force variables and the mixed convection parameter. Ibrahim et al. [5,6,7] observed heat exchange and MHD boundary layer nanofluid flow via vertical surfaces and a non–isothermal expanding surface using the similarity transformation technique. Ibrahim et al. [8,9] described nanoparticles’ MHD boundary layer flow through a stretched surface, influenced by double stratification and diffusion. Muhammad et al. [10] discussed the Eyring- Powell nanofluid with thermal radiation effects. Waqas et al. [11] checked the melting effects in a modified second-grade nanofluid. Danish et al. [12] explore the EMHD effects in Williamson nanofluid with activation energy and nonlinear thermal radiation. Waqas et al. [13] analyzed the flow and heat transfer improvement in radiative hybrid nanofluid flow across the disk. Mabood et al. [14] analyzed radiative nanofluid across the thin needle.

Several investigators studied heat transfer due to its various applications in applied technologies, manufacturing, and refrigeration. Thermally conductive fluid has subsequently become a common research topic for scientists. Nanotechnology is becoming a subject of significant interest in manufacturing and engineering for investigators as its requirement grows in industrial sectors. Non-Newtonian fluids are now a big subject for scientists and physiologists [15–17]. Moreover, the world is experiencing an energy problem because of global heating and pollution. As a result, researchers are working to develop new heat transfer techniques to make improvements that can be justified. Nanotechnology is one of the researchers' most recent techniques to help alleviate energy shortages. Choi et al. [18] introduced the nanofluids, and they have since achieved boundless attention. Nanofluids are the most using coolants in the industry. The thermal properties of a combination of nanometer-sized particles immersed in a fluid identified as nanofluids or base fluids improved due to the involvement of these nanomaterials. Buongiorno [19] explored how to augment the convection rate of heat transformation in two significantly slip situations: thermophoresis diffusion and Brownian motion behavior.

Bioconvection occurs because the microorganisms are denser than water; therefore, a motile gyrotactic microorganism moves upwardly in water. The production of the base fluid increased as the dynamically self-propelled motile gyrotactic microorganisms moved in a particular direction. Bioconvection is used in various areas, including fertilizer, biofuel, and manufacturing systems. Nanofluid bioconvection is considered for density stratified as well as patron formation due to nanomaterials apprentice, buoyancy forces, including microorganisms. Kuznetsov [20] reported that bioconvection in nanofluids is expected to occur if the concentration of nanoparticles is small so that nanoparticles do not cause any significant increase in the viscosity of the base fluid. Ali et al. [21] studied the bio-convectional of micropolar nanofluids in water-based concentrations containing microorganisms and nanoparticles. Nayak et al. [22] scrutinized the effects of velocity, temperature, solutal and motile gyrotactic microorganism slip lying on the biologically neutral Casson nano liquid flow across an exponential expanding electro- magnetic surface with chemical effects. Khan et al. [23] analyzed the bioconvection impact in the rheological nano liquid stress coupled with Arrhenius activation energy and Wu’s slip with magnetic properties. The bioconvection motion of nanofluids caused by gyrotactic motile microorganisms was mentioned by Kotha et al. [24]. The motile gyrotactic microorganisms and thermal radiation for the chemical bonded bio-convectional flow of tangent hyperbolic nanofluids were reported by Al-Khaled et al. [25]. For dispersing gyrotactic motile microorganisms and the density of gyrotactic microorganisms, including velocity profiles, are declining, the Schmidt number is performed by Makinde et al. [26]. The bioconvection Rayleigh number and the local thermal Grashof number are diminishing functions for motile gyrotactic microorganism diffusion is evaluated by Makinde et al. [27]. Saleem et al. [28] examined the density computed for motile gyrotactic microorganisms shrinks above the Peclet number under the bio convection effect. Waqas et al. [29] examined bioconvection effects on modified second-grade nanofluid microorganisms. Chu et al. [30] discussed the modified diffusion in tangent hyperbolic nanofluid under the bioconvection phenomenon. Many researchers [31–38] explore the bioconvection phenomenon in nanofluid flow.

Furthermore, affiliation with a related domain that is closely suitable in real life has a significant impact. Because of their significance in multiple fields of technology agriculture, including manufacturing applications, the flows generated through a wedge, cone, and vertical plate have gotten a lot of attention. Vajravelu et al. [39] developed a mathematical mechanism for transmitting heat reacting to flow via a cone and wedge. Chamkha [40] explained that the heat source parameter improved the local Nusselt number by considering the flow, including porosity. Many investigators in various fields, specifically fluid dynamics, used wedge as the geometry in their experiments because they implemented the above–mentioned work. Khan et al. [41] observed the flow of a nano liquid throughout a rotating wedge under the chemical process, including thermal radiation. Interestingly, fluid dynamics scientists are currently considering the action of heat transformation over wedge flows.

Whereas mass and heat transport occurs continuously in a moveable fluid. Compared to the Dufour effect, the Soret effect is negligible in fluids. Numerous researchers have reported cross-diffusion influences in boundary layer flow in recent decades. Abd El-Aziz [42] discussed the consequences of Dufour, radiation, Soret parameters, and Magnetohydrodynamic over a porous sheet that is often stretched. He discovered that the Dufour and Soret parameters have an important impact on temperature, momentum, and assemble distributions. Eckert et al. [43] and Mortimer et al. [44] have addressed a variety of flow events, including those in geophysics, where thermal diffusion and diffusion thermo consequences are highly suspected. Thermal–diffusion and Diffusion–thermo are imperative when density differences occur in the fluid flow, according to Devi et al. [45]. Thermal–diffusion and Diffusion–thermo influences have been most significant for interfering molecular weight gases in binary processes frequently confronted in chemical technology mechanisms.

MATLAB is an essential software tool for numerically solving and comparing the outcomes of various models, including physics, technology, and a variety of others. The variable flowing and thermal characteristics when a decelerating rotating disc is modified scrutinized by Rafiq et al. [46]. The consequence of thermal conductivity and the variable viscosity overheat transmission rates and drag was discovered in a two-way analysis for a numerical solution using Matlab-bvp4c and R-K shooting methods.

The study of bioconvection in nanofluid is beneficial in different fields like biofuels cells, biotechnology, bioengineering, oil refineries, biofuels, and agriculture engineering. To the author’s best knowledge, no studies have yet reported analyzing bioconvections for the heat and mass transfer of the MHD flow with cross-diffusion impacts by considering three different geometries by considering the thermal and solutal Robin’s conditions. So, we filled this gap with the assistance of the above-cited studies. The prime motive of this study is to investigate the consequences of the occurrence of motile microorganisms on nanofluid flow with thermal radiation, activation energy, and melting phenomena past the geometries cone, wedge, and the plate. Appropriate similarity transformations are used to convert the PDEs into ODEs, and these equations are solved by numerical technique through shooting methods using bvp4c solver. Outcomes are exposed graphically and conferred numerical values for various physical parameters.

2  Mathematical Development of the Problem

Assume an incompressible, convective, time independent, laminar MHD flow over geometries of three different shapes like cone, wedge and plate. The melting mechanism is measured. In the plane xy, we considered surface of the body is measured by x and normally outward to the surface measured by y. Fig. 1 shows that strength B0 of transverse electro–magnetic field is applying perpendicularly along the flow on surface of the geometry. The half angle of cone represented by α and full angle of wedge represented by Ω, the cone radius shown by r. Let ξ be the stream velocity function. The temperature of fluid, concentration and bioconvection near the plate is presented as Tw=T0+axs,Cw=C0+bxq and Nw=N0+cxn, here (a, b and c) are real constants, (s,q and n) are the energy, concentration, and bioconvection parameters for wall (the energy, concentration, and bioconvection over wall is constant when (s=q=n=0), (T0,C0 and N0) are the ambient temperature, concentration and bioconvection. Mixed convection heat flux, cross diffusion impacts are included in the given problem.

For the above considered problem the basic coupled nonlinear dimensional equations are given below [47]:

∂∂x(rm∂ξ∂y)−∂∂y(rm∂ξ∂x)=0,(1)

∂ξ∂y∂2ξ∂x∂y−∂ξ∂x∂2ξ∂y2=ν∂3ξ∂y3+(βt(T−T∞)+βc(C−C∞)+βn(N−N∞))gcos⁡α−σβ02ρ∂ξ∂y,(2)

∂ξ∂y∂T∂x−∂ξ∂x∂T∂y=kρcp∂2T∂y2+1(ρcp)f∂∂y(K(T)∂T∂y)+τ[DB∂C∂y∂T∂y+DTT∞(∂T∂y)2],(3)

∂ξ ∂y ∂C ∂x − ∂ξ ∂x ∂C ∂y = 1 (ρ c p ) f ∂ ∂y ( D(C) ∂C ∂y )+ D B ∂ 2 C ∂ y 2 + D T T ∞ ( ∂ 2 T ∂ y 2 ) − K 1 r 2 (C− C ∞ ) ( T T ∞ ) 2 exp( − E a κT ), (4)

∂ξ∂y∂N∂x+∂ξ∂x∂N∂y+[∂∂y(N∂C∂y)]bWc(Cw−C∞)=Dm∂∂y(∂N∂y).(5)

Figure 1: Schematic flow with coordinate system

The applied boundaries conditions are [48]

u(x,0)= u w =νx/ l 2 ,v(x,0)=0, k ∂T ∂y = h f ( T w −T), D B ∂C ∂y = h g ( C w −C), D m ∂N ∂y = h n ( N w −N), u→0,T→ T ∞ ,C→ C ∞ ,N→ N ∞ asy→∞. (6a)

−k ∂T ∂y | y=0 =ρ[ λ 0 +( T m − T 0 ) c s ] (6b)

1.    If m=0 and α≠0 then flow problem for wedge.

2.    If m=1 and α≠0 then flow problem for cone.

3.    If m=0 and α=0 then flow problem for plate.

Introduced the suitable similarities variables

ζ=yl−1,u=νxl−2∂f∂ζ,v=−ν(m+1)l−1f(ζ),θ(ζ)(Tw−T∞)+T∞=T,ϕ(ζ)(Cw−C∞)+C∞=C,χ(ζ)(Nw−N∞)+N∞=N.(7)

Transformation equations are as follows:

d3fdζ3+(m+1)fd2fdζ2−(dfdζ)2−Mdfdζ+λ(θ+Grϕ+Ncχ)cos⁡α=0,(8)

(1+ε1θ)d2θdζ2+ε1(dθdζ)2+Pr(Nbdθdζdϕdζ+Nt(dθdζ)2)+Pr((m+1)fdθdζ−sdfdζθ)=0,(9)

(1+ε2ϕ)d2ϕdζ2+ε2(dϕdζ)2+LePr((m+1)fdϕdζ−sdfdζϕ)+NtNbd2θdζ2−LePrσ∗∗(1+δθ)nϕexp⁡(−E1+δθ)=0,(10)

d2χdζ2+Lb((m+1)fdχdζ−sdfdζχ)−Pe[d2ϕdζ2(χ+δ1)+dχdζdϕdζ]=0,(11)

where s representing the temperature parameter.

The dimension-less conditions are as follows:

Meθ′(0)+Prf(0)=0,dfdζ=1,dθdζ=−Ω1(1−θ(0)),dϕdζ=−Ω2(1−ϕ(0))=0,dχdζ=−Ω3(1−χ(0))atζ=0(12)

dfdζ=0,θ=0,ϕ=0,χ=0asζ→0.(13)

In above equations M is the magnetic field, l is the characteristic length, λ is mixed convection parameter, Gr designed for Grashof number for concentration, Nc describe the Grashof number for microorganisms, Pr explain the Prandtl number, Nb demonstrate the Brownian motion, Nt illustrate the thermophoresis diffusion parameter, Lb confirm the bioconvection Lewis number, Pe represent the Peclet number, δ is the temperature difference parameter, microorganisms difference parameter denoted by δ1 and thermal Biot number characterized by Ω1, solutal Biot number indicated by Ω2 and the microorganisms Biot number symbolized by Ω3, ε1 is the temperature dependent thermal conductivity, ε2 denotes the Concentration diffusivity, σ∗∗ denotes the chemical reaction parameter and melting parameter can be defined as Me.

M=σ0β02l2ρν,λ=l2gβt(Tw−T∞)νuw,Gr=l2gβc(Cw−C∞)νuw,Nc=l2gBn(Nw−N∞)νuw,Pr=μcpk,Nb=τDB(Cw−C∞)ν,Nt=τDT(Tw−T∞)νT∞,Lb=νDm,Pe=bWcν,δ1=N∞Nw−N∞,Ω1=hfkl,Ω2=hgDBl,Ω3=hnDMl,Me=cp(T∞−Tm)λ0+cs(Tm−T0).(14)

The friction factor Cf, local Nusselt number Nu, local Sherwood number Sh and Sn defined as

Cf=Cf∗μuw=−d2fdζ2|ζ=0,Nu=hlk(Tw−T∞)=−dθdζ|ζ=o,Sh=hmlk(Cw−C∞)=−dϕdζ|ζ=0,Sn=hnlk(Nw−N∞)=−dχdζ|ζ=0.(15)

3  Numerical Scheme

In this portion the nonlinear dimensionless ordinary differential equations with appropriate boundary conditions are solved numerically through the solver bvp4c package in computational MATLAB. The higher order ODEs are converted into linear one by introducing the following variables such as

f=p,dfdζ=p1,d2fdζ2=p2,d3fdζ3=p′2,θ=p3,dθdζ=p4,d2θdζ2=p′4,ϕ=p5,dϕdζ=p6,d2ϕdζ2=p′6,χ=p7,dχdζ=p8,d2χdζ2=p′8.(16)

p2′=−(m+1)pp2+(p1)2+Mp1−λ(p3+Grp5+Ncp7)cos⁡γ,(17)

p4′=−ε1(p4)2−Pr(Nbp4p6+Nt(p4)2)−Pr((m+1)pp4−sp1p3)(1+ε1p3),(18)

p6′=−ε2(p6)2−LePr((m+1)pp6−sp1p5)−NtNbp′4+LePrσ∗∗(1+δp3)np5exp⁡(−E1+δp3)(1+ε2p5),(19)

p8′=−Lb((m+1)pp8−sp1p7)+Pe[p′6(p7+δ1)+p8p6],(20)

With

Mep4(0)+Prp(0)=0,p1=1,p4=−Ω1(1−p3(0)),p6=−Ω2(1−p5(0))=0,p8=−Ω3(1−p7(0))atζ=0(21)

p1=0,p3=0,p5=0,p7=0asζ→0.(22)

4  Description of the Results

In this section all the results of graphs are embellished and created by using the computational software MATLAB–16a (bvp4c) solver. For the understanding of figures, we must have to know about the color of curves demonstrated. The black curve demonstrates for wedge, red for cone and blue for sheet. The iterations of effective goods parameters increased or decreased with their individual behaviors. In this segment the graphic outcome interoperate for the different values of different parameters, like Pr=2.0, Pe=0.1, E=0.1, ε2=0.1, Lb=1.2, Le=1.2, Nb=0.1, Nt=0.1, Ω1=0.1, Ω2=0.2, Ω3=0.1, M=0.5, Nr=0.1, Nc=0.1, Me=0.5.

The influence of buoyancy ratio parameter Nr and velocity distribution f′ is illustrated in Fig. 2. It is noticed from this figure that for the highest values of Nr, f′ decreases (Fig. 3). Illustrates the fluctuation of the bioconvection Rayleigh number Nc with respect to the microorganism field with in the velocity profile f′. Such that the value of bioconvection Rayleigh number Nc increase due to the increased value of microorganism graph. The relationship between the buoyancy and viscosity is known as Rayleigh number within the fluids. The significant magnetic field parameter M against velocity distribution f′ has been examined in Fig. 4. The velocity of the fluid decreased with the increased value of M. Once more this behavior is defined as the Hartmann number is interlinked with the occurrence of Lorentz force. Figs. 5–8 are examines the variation of Me on motile density profile χ, concentration profile φ, temperature profile θ and velocity profile f′ respectively for three different geometries (wedge, cone and sheet). The fluctuation in Me reduces the motile density profile χ, concentration profile φ, temperature profile θ and velocity profile f′.

Figure 2: Pictogram of Nr over f′

Figure 3: Pictogram of Nc over f′

Figure 4: Pictogram of M over f′

Figure 5: Pictogram of Me over f′

Figure 6: Pictogram of Me over θ

Figure 7: Pictogram of Me over ϕ

Figure 8: Pictogram of Me over χ

Fig. 9 illustrates the impact of the Prandtl number over temperature distribution respectively for different value of m for cone, wedge and sheet. Fig. 10 demonstrates the impact of thermal Biot number via temperature distribution for different geometries for several values of m and α. For cone m=1.0 and for other two geometries its zero. Fig. 11 depicts that concentration is reduced via larger Prandtl number Pr. In Fig. 12 the activation energy E parameter is communicated for the nanoparticles volume distribution. In this graph when the value of E is increasing then the rate of reaction is decline due to higher temperature and greatest activation energy which develop the nanoparticles volume fraction. The concentration profile for ε2 is deeply observed in the Fig. 13. Concentration dependent thermal diffusivity enhances the concentration field of nanoparticles. Fig. 14 analyzed the effect of Lewis number Le over concentration ϕ. The decreasing in ϕ is noticed for demonstrative value of Le for three different geometries. The Le is defined to be the ratio of mass and thermal diffusivity. The ϕ is decreased for the higher value of Le in these geometries. Fig. 15 represented the variation of φ for the Brownian motion parameter Nb. The relationship between Nb and nanoparticles concentration is inverse. φ. From Fig. 16 the temperature distribution is upswing owing to the growing value of thermophoresis parameter Nt. For Nt the temperature profile improved, so the migration of particles from region to region exists as hot to cold region. Multi industrial processes use this concept now a days. The profile of concentration is increased for given value of Nt. The components of representational non–dimensional nanoparticles concentration ϕ is purposed through allotting unique values to Nt for three different geometries.

Figure 9: Pictogram of Pr over θ

Figure 10: Pictogram of Ω1 over θ

Figure 11: Pictogram of Pr over φ

Figure 12: Pictogram of E over φ

Figure 13: Pictogram of ε2 over φ

Figure 14: Pictogram of Le over φ

Figure 15: Pictogram of Nb over φ

Figs. 17 and 18 explain the variation of solutal Biot number Ω2, and microorganism Biot numbers Ω3 respectively for concentration and motile density profiles. The concentration and motile density profile depends upon the Biot number for all three conditions with respect to all three different geometries (wedge, cone and sheet). Here concentration and microorganism's field boosted for larger Biot numbers. Fig. 19 illustrates the impact of Peclet number on χ for three different geometries like wedge, cone and sheet. Largest values of Peclet number diminish the microorganism distribution. Fig. 20 illustrates that the bioconvection Lewis number Lb is a decreasing function for motile density profile χ of three geometries (wedge, cone and sheet). The χ diminished and the value of Lb are higher for the cause of increase in the diffusivity of the motile microorganism.

Figure 16: Pictogram of Nt over φ

Figure 17: Pictogram of Ω2 over φ

Figure 18: Pictogram of Ω3 over χ

Figure 19: Pictogram of Pe over χ

Figure 20: Pictogram of Lb over χ

5  Final Remarks

Heat and mass transformation of the nanofluid bioconvection flow under the melting mechanisms is scrutinized in the existence of cross-diffusion impacts. Mathematical modeling is fixed to scrutinize the flowing performance on vertically three different geometries (cone, wedge, and plate). Some significant remarks are synopsized below:

•   Fluid exposure, the fluid velocity is highest for the flow over the plate as compared with the other two geometries (cone and wedge).

•   The transformation rate of mass and heat is much better for the flow via cone than flow over the wedge.

•   It is investigated that both the Brownian motion diffusion parameter and thermophoretic diffusion parameter are directly proportional to the temperature of the fluid. At the same time, concentration is directly proportional to the thermophoretic parameter but inversely proportional to the Brownian diffusion parameter.

•   The temperature and concentration profile performance increased for the highest value in the thermophoretic parameter.

•   The microorganisms’ density is inversely proportional to the Peclet number and Lewis number.

•   This study analyzed the flow characteristics of non-Newtonian fluids extended for the present mathematical model.

•   Reducing the amount of melting enhanced the heat transfer rate and velocity of the fluid.

Funding Statement: This work is financially supported by the Government College University, Faisalabad and Higher Education Commission, Pakistan.

Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.

This work is licensed under a Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.