# Unsteady Bioconvection Squeezing Flow in a Horizontal Channel with Chemical Reaction and Magnetic Field Effects

Authors: Qingkai Zhao; Hang Xu; Longbin Tao
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

The time-dependent mixed bioconvection flow of an electrically conducting fluid between two infinite parallel plates in the presence of a magnetic field and a first-order chemical reaction is investigated. The fully coupled nonlinear systems describing the total mass, momentum, thermal energy, mass diffusion, and microorganisms equations are reduced to a set of ordinary differential equations via a set of new similarity transformations. The detailed analysis illustrating the influences of various physical parameters such as the magnetic, squeezing, and chemical reaction parameters and the Schmidt and Prandtl numbers on the distributions of temperature and microorganisms as well as the skin friction and the Nusselt number is presented. The conclusion is drawn that the flow field, temperature, and chemical reaction profiles are significantly influenced by magnetic parameter, heat generation/absorption parameter, and chemical parameter. Some examples of potential applications of such bioconvection could be found in pharmaceutical industry, microfluidic devices, microbial enhanced oil recovery, modeling oil, and gas-bearing sedimentary basins.


---

## Body

## 1. Introduction

The unsteady squeezing channel flow, caused by the moving boundary under the influence of external normal stresses or vertical velocities, is often encountered in fabrication of hydrodynamical machines, accelerators, compression and injection moulding, lubrication equipment, and polymer processing and so on. Since Stefan [1] initiated the investigation of squeezing flow by using the lubrication approximation, many researchers successively considered squeezing flow problems for various geometrical configurations in many different ways. Moore [2] concluded that the Stefan equation [1] is inadequate to describe the molecular mechanism of viscosity of squeezing flows; instead, he indicated that the influences such as surface finish, viscoelastic liquids, elastomeric surfaces, and molecular effects also need to be considered partially or entirely based on the degree of complexity of problems. P. S. Gupta and A. S. Gupta [3] noticed that the unsteady squeezing channel flow problem could be simplified significantly via similarity variables in the case that the distance between the paralleled plates varies as the square root of a linear function of time. Duwairi et al. [4] considered the heat transfer effects on the unsteady squeezing channel flow, in which case they assumed that the paralleled walls are heated uniformly by a constant temperature. The unsteady squeezing flow problems were further investigated by different researchers with consideration of various physical problems and computational techniques, such as Verma [5], Mustafa et al. [6], and Hayat et al. [7].It is known that a large amount of microorganisms are suspended in appropriate aqueous environments, such as oceans and rivers, puddles, and droplets. Their behaviours govern the dynamics of many of these interactions with fluid mechanics. The hydrodynamic characteristics of single-celled microorganisms which are motile and self-propelled were firstly studied by Pedley et al. [8] and Pedley and Kessler [9]. Pedley’s theory [8] was then adopted by Kuznetsov [10] to consider the thermal bioconvection, in which a macroscopic convective motion was induced in a fluid layer by the combined effect of density stratification caused by the upswimming of oxytactic microorganisms and heating from below. Kuznetsov [11] then investigated the influences of nanoparticles on the behaviours of gyrotactic microorganisms and concluded that nanoparticles can be either stabilizing or destabilizing the density of microorganisms, depending on whether the basic nanoparticle distribution is bottom-heavy or top-heavy. His idea was confirmed by Xu and Pop [12] via an analysis of a mixed bioconvection nanofluid flow over a stretching flat sheet with uniform free stream with an improved nanofluid model proposed by Kuznetsov and Nield [13]. The novel behaviours of thermobioconvection in nanofluids were further investigated by Xu and Pop [14] and Raees et al. [15].Physically, magnetic fields can induce currents in a moving electrically conductive fluid, which in turn polarizes the fluid and reciprocally changes the magnetic field itself. Such fluid was named as magnetohydrodynamics (MHD) by Alfvén [16]. Since MHD flows have been found to be very useful in heat and mass transfer progresses in many industrial procedures such as polymer extrusion, drawing of copper wires, continuous stretching of plastic films and artificial fibers, hot rolling, wire drawing, glass-fiber, metal extrusion, and metal spinning, many studies [17–27] have been done towards understandings of their transport mechanisms and novel applications. It is worth mentioning that Khan and Makinde [28] conducted an investigation about the magnetic effect on bioconvection of a nanofluid along a vertical stretching sheet. They concluded that the ratio of convective to conductive heat transfer across (normal to) the boundary decreases with the magnetic parameter increasing, while the density number of the motile microorganisms decreases with magnetic parameter enlarging. Their work indicates that the MHD flow could be used to alter the distributions of microorganisms in potential applications such as removal of microorganisms attached in a water tank.In this paper, we shall consider an unsteady bioconvection squeezing flow in a horizontal channel in the presence of both chemical reaction and magnetic effect. In a flowing fluid that contains heat transfer, mass transfer, and chemical reactions, the flow velocity is always to affect partially the transport of reactants, products, and heat and hence influences the overall mass and heat transfer rates. An analysis for accurate estimates of the transport effects upon the performance of processes with combined chemical reaction and magnetic effect in a combined flow regime is thus very desirable. As far as we know, no such study has been available in the literature. This study attempts to meet such a need.


## 2. Mathematical Descriptions

Consider the unsteady MHD flow of an electrically conducting fluid containing motile microorganisms between two infinite parallel plates in the presence of a first-order chemical reaction. As shown in Figure1, the coordinate system is chosen such that the x-axis is along the lower plate and y-axis is perpendicular to the main flow direction, respectively. It is assumed that the two plates are placed at the distance y=ht=[ν1-at/b]1/2 while the upper plate is moving towards or away from the lower stationary plate with the velocity vt=dh/dt. Here, a and b are positive numbers, ν is the kinetic viscosity, and t denotes the time. Obviously, it is known that 1-at must be greater than zero, which indicates that a must be less than 1/t. Physically, when a=0, it means that the two plates are fixed. As 0<a<1/t, it shows that the upper plate is squeezed to the lower one. While as a<0, it indicates the upper plate moves away from the lower plate. It is also assumed that the lower and upper plates are maintained, respectively, at the constant temperatures T1 and T2, constant concentration of microorganisms N1 and N2, and constant chemical reaction concentrations C1 and C2. The time-varying magnetic field of strength Bt is imposed along the y-axis. Based on the above assumptions, the conservation equations for total mass, momentum, thermal energy, mass diffusion, and microorganisms equations are expressed by(1)∇·V=0,∂V∂t+V·∇V=-1ρ∇p+ν∇2V+1ρJ×B,∂T∂t+V·∇T=α∇2T+QstρcpT-T0,∂C∂t+V·∇C=D∇2C-KtC-C0,∂N∂t=-∇·j,where V is the velocity vector, J is the electric current density in the fluid, B is the magnetic induction, ρ is the fluid density, p is the pressure, ν is the kinematic viscosity of the fluid, T is the temperature, T0 is the reference temperature, α is the thermal diffusivity, Qst is the volumetric rate of heat generation/absorption, C is the concentration of the chemical reaction, C0 is the reference concentration, D is the mass diffusivity, Kt is the reaction rate, N is the concentration of motile microorganism, and j is the flux of microorganisms.Figure 1
Physical model and coordinate system.It is assumed that the magnetic Reynolds number is small (≪1) so that the induced magnetic field can be omitted. It is also assumed that the heating due to the plates and the chemical reaction is weak so that microorganisms can survive and their gyrotactic behaviour can keep unchanged. With those assumptions, we can expand (1) into the following forms:(2)∂ξ∂t+u∂ξ∂x+v∂ξ∂y=ν∂2ξ∂x2+∂2ξ∂y2+σB2tuρ,∂T∂t+u∂T∂x+v∂T∂y=α∂2T∂x2+∂2T∂y2+QstρcpT-T0,∂C∂t+u∂C∂x+v∂C∂y=D∂2C∂x2+∂2C∂y2-KtC-C0,∂N∂t+u∂N∂x+v∂N∂y+∂∂xNv~+∂∂yNv~=Dn∂2N∂x2+∂2N∂y2,subject to the boundary conditions(3)u=0,v=0,T=T1,C=C1,N=N1aty=0,u=0,v=dhdt,T=T2,C=C2,N=N2aty=ht,where u and v are the velocity components in the x- and y-directions defined by u=∂ψ/∂y and v=-∂ψ/∂x with ψ being the stream function, ξ=-∇2ψ is the vorticity function, N is the density of motile microorganisms, v~=[bcWc/(C1-C0)]∂C/∂y is the average swimming velocity vector of the oxytactic microorganisms with bc being the chemotaxis constant and Wc being the maximum cell swimming speed, and Dn is the diffusivity of microorganisms.We then introduce the following similarity transformations and nondimensional quantities:(4)ψx,y=bν1-at1/2xfη,u=bx1-atf′η,v=-bν1-at1/2fη,η=bν1-at1/2y,θη=T-T0T1-T0,ϕη=C-C0C1-C0,wη=NN1.In view of (4), (2) can be reduced to the following nondimensional form:(5)f′′′′η-Mf′η-3βf′′η-βηf′′′η-f′ηf′′η+fηf′′′η=0,(6)1Prθ′′η+Qθη-βηθ′η+fηθ′η=0,(7)1LePrϕ′′η-K0ϕη-βηϕ′η+fηϕ′η=0,(8)w′′η-Scβηw′η+Scfηw′η-Pewηϕ′′η-Peϕ′ηw′η=0.subject to the following boundary conditions:(9)f′0=0,f0=0,θ0=1,ϕ0=1,w0=1,f′1=0,f1=β,θ1=δθ,ϕ1=δϕ,w1=δw,where M is the magnetic field parameter, β is the squeezing parameter, Pr is the Prandtl number, Q is the heat generation/absorption parameter, Le is the Lewis number, K0 is the chemical reaction parameter, Sc is the Schmidt number, Pe is the bioconvection Péclet number, and δθ, δϕ, and δw are constants, which are defined by(10)M=σνρB02,β=a2b,Pr=να,Q=Q0ρcp,Le=αD,K0=b1-at-1Kt,Sc=νDn,Pe=bcWcDn,δθ=T2-T0T1-T0,δϕ=C2-C0C1-C0,δw=N2N1.In above equationB0=b/ν1-at-3/4Bt is a reference magnetic field parameter, Q0=b/1-at-1Qst is a reference heat generation/absorption parameter, and K0 is the chemical reaction parameter.Practical interests of this problem are the local skin frictionCfx, the local Nusselt number Nux, the local wall mass flux Shx, and the local wall motile microorganisms flux Qx, which are defined by(11)Cfx=τw1/2ρUx2,Nux=xqwkT1-T0,Shx=xjwDC1-C0,Qx=xhwDnN1-N0,where(12)τw=μ∂u∂yy=0,qw=-k∂T∂yy=0,jw=-D∂C∂yy=0,hw=-Dn∂N∂yy=0.Without loss of generality, we setN0=0. Substituting (12) into (11), we obtain(13)Rex1/2Cfx=2f′′0,Rex-1/2Nux=-θ′0,Rex-1/2Shx=-ϕ′0,Rex-1/2Qx=-w′0,where Rex=Uwx/ν is the local Reynolds number.


## 3. Results and Discussion

Equations (5)–(8) and their boundary conditions (9) are solved numerically for various values of the parameters using the shooting method along with fourth-order Runge-Kutta scheme. In the computation, the value of the step size Δη between 0.001 and 0.01 was used, in order that the numerical values obtained are mesh independent. In addition, we choose δθ=1/2, δϕ=0, and δw=1 for all considered cases, which are not against physical requirements. Further, to check the accuracy of the results, we compare them with those given by the BVPh2.0 programme (which is based on the homotopy analysis methods technique). Here, the initial approximations f0η, θ0η, ϕ0η, and w0η are chosen based on (9) as below:(14)f0η=3βη2-2βη3,θ0η=1+δθ-1η,ϕ0η=1+δϕ-1η,w0η=1+δw-1η.The auxiliary linear operator are chosen as:(15)Lfη=∂4fη∂η4,Lθη=∂2θη∂η2,Lϕη=∂2ϕη∂η2,Lw=∂2wη∂η2.Through calculation and comparison, excellent agreement is found, as shown in Figure 2.Figure 2
Comparison of the BVPh2.0 programme results and the numerical results forβ=1.5, Q=1, K0=1, and Pr=Le=Sc=Pe=1, where symbols are numerical results and lines are BVPh2.0 programme results.We first consider the influence of the magnetic parameterM on variations of various profiles. It is shown in Figure 3 that f′η develops rapidly with M evolving, the enlargement of M causes the enhancement of the flow velocity in the vicinity of the lower plate but the reduction of the flow velocity on the other side. This is due to the fact that the enhanced Lorentz force imparts additional momentum into the squeezing channel, which accelerates the flow motion in the lower part of the channel. Synchronously the flow near the upper plate is retarded owing to the nonslip boundary. It is found in Figure 4 that the temperature profiles decrease continuously as M increases. The reason is that the enhanced Lorentz force increases the thermal boundary-layer thickness in the channel. Particularly, we notice that the temperature profiles near the lower plate reduce a bit more quickly than that near the upper plate, which is due to the fact that the flow velocity in the vicinity of the lower plate is higher than that near the upper plate. As shown in Figure 5, the influence of the magnetic parameter M on the concentration of chemical reaction ϕη is similar to that on the temperature profiles θη; that is to say, the larger the value of M, the lower the profile ϕη. Physically, this variation is mainly caused by the mass transportation induced by the enhanced Lorentz force. It is illustrated in Figure 6 that M has obvious effect on variation of the motile density of the microorganisms wη. wη decreases as M enlarges, since the increase of M accelerates the fluid motion, which is helpful to alter the distributions of microorganisms.Figure 3
Variation off′η with η for various M in the case of β=1.5, Q=1, K0=1, and Pr=Le=Sc=Pe=1.Figure 4
Variation ofθη with η for various M in the case of β=1.5, Q=1, K0=1, and Pr=Le=Sc=Pe=1.Figure 5
Variation ofϕη with η for various M in the case of β=1.5, Q=1, K0=1, and Pr=Le=Sc=Pe=1.Figure 6
Variation ofwη with η for various values of parameter M in the case of β=1.5, Q=1, K0=1, and Pr=Le=Sc=Pe=1.We next consider the effects of the chemical reaction parameterK0 and the heat generation/absorption parameter Q on various profiles. In view of (5), (6), and (7), it is readily found that fη and θη are uncoupled with ϕη; that is to say, the variation of ϕη due to the evolvement of K0 does not affect the velocity and temperature distributions at all. Therefore, we need only to consider the influence of K0 on ϕη and wη. As shown in Figure 7, the increase of K0 causes the decrease of ϕη. Physically, K0 is related to the chemical reaction rate. The larger this value is, the faster the reaction rate is. As a result, the concentration of the chemical species is reduced as K0 evolves. It is illustrated in Figure 8 that the density of motile microorganisms wη reduces as K0 enlarges. This is due to the fact that the enhancement of K0 reduces the concentration of chemical species, which leads to the decrease of the density of motile microorganisms. Similarly, it is found from (5)–(8) that the heat generation/absorption parameter Q is only associated with θη via the reduced energy equation (8) and θη has no relation with any other functions. We therefore only need to consider its influence on the temperature distribution. As shown in Figures 9 and 10, for Q>0, corresponding to the heat source generation case, θη increase rapidly with Q increasing. On the other hand, for Q<0, corresponding to the heat source absorption case, θη reduces rapidly with Q decreasing. In other words, the heat generation/absorption parameter Q has a great influence on the temperature.Figure 7
Variation ofϕη with η for various values of K0 in the case of β=1.5, Q=1, M=1, and Pr=Le=Sc=Pe=1.Figure 8
Variation ofwη with η for various values of K0 in the case of β=1.5, Q=1, M=1, and Pr=Le=Sc=Pe=1.Figure 9
Variation ofθ(η) with η for various positive values of Q in the case of β=1.5, M=1, K0=1, and Pr=Le=Sc=Pe=1.Figure 10
Variation ofθη with η for various negative values of Q in the case of β=1.5, M=1, K0=1, and Pr=Le=Sc=Pe=1.We now discuss the effects of various physical parameters on physically important quantities. It is known thatM and δθ can affect the local Nusselt number Rex-1/2Nux, according to (6) and its boundary condition (9) for temperature distribution on the upper plate. As shown in Table 1, the local Nusselt number increases gradually as the magnetic parameter M enlarges. This is due to the fact that the increase of M means more additional momentum caused by the enhanced Lorentz force is imparted into the sneezing channel. As a result, the fluid flow is accelerated and the temperature difference between this plate and the reference value is increased. Similarly, we know that M and δϕ take effects on the local wall mass flux Rex-1/2Shx. It is seen in Table 2 that the local wall mass flux increases continuously as M enlarges. This trend indicates that the magnetic parameter plays a role on variation of the nanoparticles distribution. When M is fixed, it is found in this table that the local wall mass flux reduces gradually as δϕ evolves.Table 1
The local Nusselt numberRex-1/2Nux for various values of parameters δθ and M in the case of δϕ=0, δw=1, Q=K0=Pr=Le=Sc=Pe=1, and β=1.5.
δθ M = 0 M = 50 M = 100 M = 200 0 0.62504 0.68927 0.75182 0.85656 0.2 0.39363 0.44398 0.49285 0.57423 0.4 0.16221 0.19868 0.23388 0.29191 0.6 −0.06921 −0.04661 −0.02509 0.00959 0.8 −0.30062 −0.29191 −0.28406 −0.27273 1 −0.53204 −0.53721 −0.54303 −0.55505Table 2
For differentM, the values of the local wall flux Rex-1/2Shx for various values of parameter δϕ in the case of δθ=0.5, δw=1, Q=K0=Pr=Le=Sc=Pe=1, and β=1.5.
δϕ M=0 M=50 M=100 M=200 0 1.28437 1.33808 1.39126 1.48271 0.2 1.11792 1.16167 1.20508 1.28000 0.4 0.95147 0.98525 1.01890 1.07729 0.6 0.78501 0.80883 0.83271 0.87459 0.8 0.61856 0.63242 0.64653 0.67188 1 0.45210 0.45600 0.46035 0.46918In practical application, It is interesting to give some approximate formulae for calculating those physical quantities. For instance, by regression analysis, we obtain the estimated formulae for various physical quantities as functions ofM and/or K0 in the case of Q=Pr=Le=Sc=Pe=1 and β=1.5, (16)Rex1/2Cfx=9.035322+0.177944M,Rex-1/2Nux=0.04904+0.00524M,Rex-1/2Shx=1.1097864+0.0081446M+0.20924432K0,Rex-1/2Qx=0.1403712+0.00089064M+0.22188504K0.Similarly, we are able to give the estimated formulae by linear regression for other physical parameters. For instance, we are able to obtain the approximate expressions regarding δθ, δϕ and δw by defining(17)PhQ=B0+B1δθ+B2δϕ+B3δw,where PhQ can be replaced by various quantities including the local Nusselt number, the local wall mass flux, and the local wall motile microorganisms. When the value of δθ, δϕ, and δw are, respectively, taken as 0, 0.4, 0.8, 1.2, 1.6, an 2, the coefficients B1, B2, and B3 can be computed by linear regression, as shown in Tables 3–5. Note that here we set Q=Pr=Le=Sc=Pe=K0=1 and β=1.5.Table 3
Linear regression coefficients and the maximum relative errorε for the local Nusselt number Rex-1/2Nux.
M B 0 (y-intercept) B1 B2 B3 ε 0 0.6250 −1.1571 0 0 5.64E-07 50 0.6893 −1.2265 0 0 5E-07 100 0.7518 −1.2949 0 0 1.57E-16 200 0.8566 −1.4116 0 0 5E-07 500 0.9942 −1.5765 0 0 5E-07Table 4
Linear regression coefficients and the maximum relative errorε for the local wall flux Rex-1/2Shx.
M B 0 (y-intercept) B 1 B 2 B 3 ε 0 1.2844 0 −0.8323 0 5.39 E - 07 50 1.3381 0 −0.8821 0 4.59 E - 07 100 1.3913 0 −0.9309 0 5.82 E - 07 200 1.4827 0 −1.0135 0 5.39 E - 07 500 1.6144 0 −1.1258 0 5.12 E - 07Table 5
Linear regression coefficients and the maximum relative errorε for the local wall motile microorganisms flux Rex-1/2Qx.
M B 0 (y-intercept) B 1 B 2 B 3 ε 0 1.3019 0 0.1535 −1.0013 0.0401 50 1.4114 0 0.1082 −1.0534 0.0407 100 1.5201 0 0.0628 −1.1044 0.0411 200 1.7075 0 −0.0158 −1.1907 0.0419 500 1.9784 0 −0.1263 −1.3083 0.0431


## 4. Conclusions

The unsteady squeezing flow of an incompressible viscous fluid with suspension of small motile microorganisms between two infinite parallel plates in the presence of both chemical reaction and magnetic effect has been studied. The governing equations embody the total mass, momentum, thermal energy, and mass diffusion and microorganisms are reduced to a set of nonlinear ordinary differential equations via a set of similarity transformations. The effects of the various parameters on the distributions of velocity, temperature, chemical reaction concentration, and density of motile microorganisms, as well as the important physical quantities, have been examined in detail. The major contributions can be summarized as follows:(i)
This problem solved by two methods and the results obtained by numerical methods have good agreement with the analytical solution.(ii)
The magnetic parameterM takes actions on each profile since it imparts additional momentum into the squeezing channel. The increase of M leads to the enhancement of the velocity profile, but causes the decrease of the profiles of the temperature, the chemical reaction concentration, and the motile density of microorganisms. In short, the magnetic field can be used to control the flow characteristics.(iii)
The chemical parameterK0 plays a role in the chemical reaction concentration and motile density of microorganisms profiles. The increase of K0 causes the decrease of both profiles.(iv)
The heat generation/absorption parameterQ is critical for the profiles of the temperature. The temperature is found to increase as we go from strong absorption to strong generation of heat.(v)
The estimated expressions for physical quantities of practical interests are given by linear regression for a range of parameters.



---

*Source: 2541413-2017-01-24.xml*