Electromagnetohydrodynamic bioconvective flow of binary fluid containing nanoparticles and gyrotactic microorganisms through a stratified stretching sheet.

Bioconvection has recently been the subject of dispute in a number of biotechnological fields that depend on fluids and their physical properties. When mixed nanofluids are subjected to heat and mass transmission, the process of bioconvection occurs. This attempt conveys the theoretical analysis of two-dimensional electrically conducting and magnetically susceptible binary fluid containing nanoparticles and gyrotactic microorganisms past a stratified stretching surface. Furthermore binary chemical reaction, thermal radiation, and activation energy are taken into assumptions. The analytical solution based on HAM has been performed. The convergence of HAM is presented with the help of figures. The present study is compared with previously published results and has established an excessive agreement which validate the present study. It is perceived that the presence and absence of an electric field influences the variations in fluid velocities due to presence of magnetic field. The micropolar constant heightens the velocity and microrotation of the fluid flow. The buoyancy parameter and bioconvection Rayleigh number diminish the velocity function while these parameters show dual impact on microrotation function. The skin friction and couple stress escalates with the increasing buoyancy ratio parameter and bioconvection Rayleigh number.

Introduction

The expression "nanofluid" refers to the suspended nanoparticles that maximize the combined heat and mass transfer phenomena within a typical fluid. Nanoparticles have fascinated the researchers’ interest in today’s modern era due to their significant importance in the fields of electronics, food science, biosensors, biomedicine, and mechanical engineering. Additionally, the movement of respective nanoparticles in designated structures is strongly dependent on the elementary concepts of cancer treatment, selective drug delivery, chemotherapy, fermentation science, and nano-medicine. It is a well-established fact that the fluids flowing through microchannels in cooling and heating systems are entirely dependent on the heat transfer particles produced by nanofluids. As a result, nanofluid dynamics is a critical term to grasp for all fields concerned with nonmaterial suspensions in some manner in order to achieve optimum productivity. Choi[1] pioneered the concept of nanoparticles with improved thermophysical properties, which was subsequently expanded by a number of scientists. Buongiorno[2] defined the seven slip mechanisms in nanoparticle movement, focusing on Brownian motion and thermophoresis effects. Hsiao[3] examined mixed convection and slip flow in flow of nanofluid configured by a stretched surface in the existence of both electrical and magnetic field aspects. Turkyilmazoglu[4] presented the heat transformer of several nanofluids containing Ag, Al2O3, CuO, Cu, and TiO2 nanoparticles through a plane wall jet. Hsiao[5] probed the hydromagnetic nanofluids flow heat and mass transfer through a stretching sheet with magnetic and viscous dissipation effects. In order to investigate the thermophoresis and Brownian movement of the nanoparticles, Ahmed et al.[6] used the Buongiorno's model for the stagnation point Maxwell nanofluid flow past a rotating disk. Sandeep and Animasaun[7] probed the enhanced thermal transmissions of electrically conducting water based nanofluids containing aluminum alloy nanoparticles with magnetic field impact. The effect of radiation thermal transmission on nanofluid between two pipes with horizontal magnetic field is investigated by Sheikholeslami et al.[8]. Shahzadi and Nadeem[9] presented an effective mathematical model for the blood based peristaltic flow containing two different types of nanoparticles past a porous material with velocity slip conditions and magnetic impact. Sandeep et al.[10] addressed the water based nanofluids containing magnetite nanoparticles. Sheikholeslami and Bhatti[11] investigated the heat transmission in a nanofluid over a porous semi-annulus with magnetic field. Raza et al.[12] deliberated the influence magnetic field on Casson fluid flow containing suspended nanoparticles past a nonlinear permeable surface with velocity slip condition. Siavashi et al.[13] offered the mixed convection flow of power-law nanofluid containing CuO nanoparticles in a porous enclosure. In addition, the related studies are mentioned in[14-19].

The activation energy is described as the least energy needed to initiate a chemical reaction. Spontaneous reactions are chemical reactions that require lower activation energy. Nuclear reactions accompanying both fusion and fission of nuclei are crucial, but they necessitate greater activation energy. The importance of activation energy has a significant impact on nanoparticle movement in simple carrier fluids. Chemical engineering, food manufacturing, and the mechanics of oil water emulsions all have a strong demand for activation energy[20-22]. Makinde et al.[23] analyzed the chemically reactive and thermally radiative unsteady fluid flow over a porous plate. Khan et al.[24] presented the features of activation energy and entropy generation in Carreau-Yasuda fluid past an extending sheet with magnetic field. The significance of activation energy on an electrically conducting magnetized third grade fluid containing gyrotactic microorganisms over a stretching sheet was probed by Chu et al.[25]. Moreover, the related studies can be found in[26-35].

Bioconvection has recently been the subject of dispute in a number of biotechnological fields that depend on fluids and their physical properties. When mixed nanofluids are subjected to heat and mass transmission, the process of bioconvection occurs. Kuznetsov[36] proposed the bioconvection in a nanofluid having gyrotactic microorganisms. Kuznetsov[37] developed the same idea by adding nanoparticles in order to the stability of nanofluid. Xun et al.[38] analyzed the bioconvective fluid flow containing gyrotactic microorganisms between two rotating plates. Further studies based on nanofluids containing gyrotactic microorganisms are mentioned in[39-45].

Based on the literature review, the theoretical analysis of two-dimensional electrically conducting and magnetically susceptible viscoelastic micropolar nanofluid containing nanoparticles and gyrotactic microorganisms through a stratified stretching sheet has not been performed yet. Thus, the authors have presented the viscoelastic micropolar nanofluid containing nanoparticles and gyrotactic microorganisms through a stratified stretching sheet. Furthermore bioconvection, binary chemical reaction, thermal radiation, and activation energy influences are taken into assumptions.

Problem formulation

Let us assume the incompressible and electrically conducting MHD two-dimensional bioconvective viscoelastic micropolar nanofluid flow containing gyrotactic microorganisms which propagate over a stratified stretching sheet. Microorganisms are brought to become the nanoparticles stable. The nanoparticles do not affect the microorganisms’ velocity and swimming direction. The stretching velocity of the sheet is assumed as along direction, whereas direction is normal to the nanofluid flow. Magnetic and electric fields are applied normal to the nanofluid flow. Furthermore, binary chemical reaction, activation energy, thermal radiation, mixed convection, and Joule heating influences are taken into consideration. Figure 1 indicates the physical representation of the flow problem. Following the above assumptions, the leading equations are formulated as[33,46-48]:

Figure 1

Physical illustration of the flow problem.

Continuity equation

Momentum equation

with boundary conditions:

The correspondence transformations are defined as[5,49-51]:

Using (4), (1) is identically contented, (2) and (3) are transformed as: where is the micropolar constant, is the viscoelastic parameter, is the magnetic parameter, is the electric parameter, is the mixed convection parameter, is the buoyancy ratio parameter, and is the bioconvection Rayleigh number.

Angular momentum equation

with boundary conditions:

Using (4), (7) and (8) are transformed as:

Temperature equation

with boundary conditions:

Here, is the radiative heat flux and is defined as:

Expending by mean of Taylor series, we have:

Using (4), (11) and (12) are transformed as:where is the thermal radiation parameter, is the Prandtl number, is the thermal stratification parameter, is the Eckert number, is the heat generation parameter.

Concentration equation

with boundary conditions:

Using (4), (17) and (18) are transformed as:where is the Schmidt number, is the reaction rate parameter, is the activation energy parameter, is the mass stratification parameter, and is the temperature difference parameter.

Motile density equation

with boundary conditions:

Using (4), (21) and (22) are transformed as:with transformed boundary conditions:where is the bioconvection Lewis number, is the bioconvection Peclet number, microorganisms’ concentration difference parameter, and is the motile density stratification parameter.

The expressions of skin friction, couple stress, Nusselt number, Sherwood number, and density number can be written as:and

Using (4), (25)-(29) are transformed as:where .

HAM solution

The linear operators and initial guesses are taken as:withwhere are called constants.

The zeroth order deformation can be written aswhere is the embedded parameter and , , , , and are auxiliary factors.

The th order deformation can be written as:

For and , we can write:

When varies from 0 to 1, the solutions varies from initial to final. Using Taylor’s series to expand the solutions: i.e.

The series (73–77) converges by choosing , i.e.

HAM convergence

HAM solution is operated to investigate the analytical solution of the present analysis. HAM is associated with the axillary parameters , , , , and . These parameters are responsible to utilize and control the convergence area of the series solutions. The convergence areas of the velocity, microrotation, temperature, concentration, and motile density functions are , , , and (see Fig. 2a–e).

Figure 2

(a) curve for . (b) curve for . (c) curve for . (d) curve for . (e) curve for .

Validation of the present analysis

An analytical scheme called HAM has been focused for the solution of the flow problem. The present analysis has been compared with previously published results by Wakif et al.[50], and Eldabe and Ouaf[52] and found a great agreement (see Tables 1 and 2).

Table 1

Numerical comparison of the present values with Wakif et al.[50] when .

α_1	M	-f^″0
		[50]	Present values
0.0	0.0	1.000000000	1.000000000
0.0	16.0	4.123105625	4.123105625
0.1	16.0	4.346134936	4.346134936
0.2	16.0	4.609772228	4.609772228
0.3	16.0	4.928053803	4.928053803
0.1	4.0	1.054092553	1.054092553
0.1	9.0	3.333333333	3.333333333
0.1	16.0	4.346134936	4.346134936

Table 2

Numerical comparison of the present values with Hsiao[5] when .

M	K	g^′0
		[52]	Present values
0.0	0.1	0.09500	0.095000
0.5	0.1	0.10509	0.105091
1.0	0.1	0.11212	0.112131
0.1	0.0	0.00000	0.000000
0.1	0.5	0.21116	0.211169
0.1	2.0	0.35855	0.358562

Results and discussion

In this section, variations in the flow profiles of binary fluid due to embedded parameters defined in Table 3, are displayed with the help of Figures, and discussed in detail. Figure 3a shows the flow profiles of the modeled problem. Here, we have verified the boundary conditions for , , , , and . Figure 3b–f show the variation in , , , , and due to . Here, it is found that the increasing micropolar constant increases the velocity and microrotation functions (see Fig. 3b,c). Physically the greater change in micropolar constant is associated with the fluid viscosity. The higher change in micropolar constants diminishes the binary fluid viscosity due to which the velocity function heightens. Consequently, an increasing impact in velocity function due to micropolar constants is depicted (see Fig. 3b). A similar impact of micropolar constant against microrotation function is observed as shown in Fig. 3c. In addition, the binary fluid has maximum viscosity when the micropolar constant converges to zero (i.e.). The decreasing impact in temperature, concentration, and motile density functions due to micropolar constant is shown in Figs. 3d–f. Actually, the higher change in micropolar constants reduces the viscosity of the nanoparticles of due to which the temperature and concentration of the nanoparticles reduces. Thus, decreasing impacts in temperature and concentration functions are depicted. A similar impact of micropolar constant against motile density function is depicted (see Fig. 3f). Figures 4a–e) show the variation in , , , , and due to . A decreasing impact in velocity profile is depicted here. This impact is due to the buoyancy forces which lead the velocity function to decrease. However, the microrotation function increases from while a decreasing impact is depicted as . Furthermore , , and functions increase with greater . Figure 5a–e show the variation in , , , , and due to . The decreasing impact of on and is depicted (see Fig. 5a, b). Physically, higher values of heighten the buoyancy forces which conclude the decreasing impact in . However, increases from while a decreasing impact is depicted as . Furthermore, , , and increase with greater . Figure 6a, b show the variation in via in the absence and presence of . reduces with greater in the absence of . The Lorenz Force which is greater with a larger magnetic field is obviously dependent on the magnetic field. The Lorenz force augments the frictional force in the absence of , acting as a retarding force that opposes the flow of binary fluid. In the presence of , the magnetic parameter has increasing impact on velocity profile (see Fig. 6b). Physically, increases the body force which accelerates the flow velocity. Thus, an increasing conduct is observed here. A similar impact is depicted in[53]. Figures 6c,d show the variation in via in the absence and presence of . In the absence of , a decreasing impact of on is depicted (see Fig. 6c). However in the presence of , magnetic field parameter has dual impact on (see Fig. 6d). Figure 7a, b display the variation in and via . The greater increases the velocity function whereas the microrotation function reduces with greater electric field parameter. With larger values of , the velocity function boosts. Physically, the electric field acts as a decreasing force and lowering the fluid's frictional strength and thereby increasing the fluid velocity. However, electric field has the opposite effect on microrotation function. Figure 8a–c denote the variation in , , and via thermal, mass, and motile density stratification parameters, respectively. Physically, the thermal stratification parameter increases the fluid density which results reduction in the temperature function. A decreasing impact of mass, and motile density stratification parameters are also depicted on concentration and motile density functions. Figure 9a–c show the variation in due to , , and respectively. The influence of on is displayed in Fig. 9a. The greater thermal radiation increases the temperature function. Physically, the greater thermal radiation means production of more heat to the fluid flow system. Thus, the greater thermal radiation parameter heightens the temperature function. Figure 9b displays the inspiration of on . The greater increases the temperature function. Physically, the increasing dissipation improves the thermal conductivity of the binary fluid which consequently enhances the temperature boundary layer thickness. Figure 9c shows the influence of on . The greater increases . Physically, increasing releases energy to the fluid flow which accordingly heightens the thermal function. Thus, increases with greater . Figure 10a–c represent the variation in due to , , and respectively. The impact of on is displayed in Fig. 10a. The greater reduces . Physically, the increasing chemical reaction results the thickening in a concentration boundary layer. With greater the term increases significantly. Thus, the greater chemical reaction parameter declines the concentration function. Figure 10b shows the impact of on . increases with greater . Actually, the greater activation energy heightens the thickness of mass transport boundary layer which consequently increases . Figure 10c specifies the deviation in due to . The greater reduces . The Schmidt number has opposite relation with mass diffusivity which consequently reduces . Figure 11a,b show the variation in due to , and respectively. The higher reduces . Physically, the diffusivity of microorganisms declines with greater which consequently reduces . Figure 11b shows the variation in due to . It is found that reduces with greater . Figures 12a–d specify the impacts of , , , and on skin friction. The increasing , and increases the skin friction whereas the increasing reduces the skin friction. Figures 13a,d specify the impacts of , , , and on couple stress. The increasing and increases the couple stress whereas the increasing and reduces the couple stress. Figures 14a–f specify the impacts of , , , , and on Nusselt number. The increasing , , and increases the Nusselt number whereas the increasing , and reduces the Nusselt number. Figure 15a–c specify the impacts of , , and on Sherwood number. The increasing and increases the Sherwood number whereas the increasing reduces the Sherwood number. Figure 16a–d specify the impacts of , , and on density number. The greater increases the density number, while , and has reverse impacts on density number. Tables 4, 5, 6, 7 and 8 show the numerical values of , , , , and via different embedded parameters. The outcomes are discussed in Figs. 12, 13, 14, 15, 16.

Table 3

Embedded parameters and their default values.

Parameter	Default value	Parameter	Default value	Parameter	Default value
K	1.0	Nr	0.8	Ec	0.4
α_1	0.5	Rb	0.7	δ	0.6
M	0.5	Rd	0.3	Sc	0.5
E¯	0.2	Pr	1.0	ϖ	0.3
λ	0.9	S	0.2	E	0.1
Q	0.2	Lb	0.7	Ω	0.2
ε	0.1	Pe	0.4	B	0.2

Figure 3

(a) Flow profiles for , , , , and , when . (b) Alteration in via . (c) Alteration in via . (d) Alteration in via . (e) Alteration in via . (f) Alteration in via .

Figure 4

(a) Alteration in via . (b) Alteration in via . (c) Alteration in via . (d) Alteration in via . (e) Alteration in via .

Figure 5

(a) Alteration in via . (b) Alteration in via . (c) Alteration in via . (d) Alteration in via . (e) Alteration in via .

Figure 6

(a) Alteration in via when . (b) Alteration in via when . (c) Alteration in via when . (d) Alteration in via when .

Figure 7

(a) Alteration in via . (b) Alteration in via .

Figure 8

(a) Alteration in via . (b) Alteration in via . (c) Alteration in via .

Figure 9

(a) Alteration in via . (b) Alteration in via . (c) Alteration in via .

Figure 10

(a) Alteration in via . (b) Alteration in via . (c) Alteration in via .

Figure 11

(a) Alteration in via . (b) Alteration in via .

Figure 12

(a) versus . (b) versus . (c) versus . (d) versus .

Figure 13

(a) versus . (b) versus . (c) versus . (d) versus .

Figure 14

(a) versus . (b) versus . (c) versus . (d) versus . (e) versus . (f) versus .

Figure 15

(a) versus . (b) versus . (b) versus .

Figure 16

(a) versus . (b) versus . (c) versus . (d) versus .

Table 4

Numerical values of via different embedded parameters.

λ	Nr	Rb	E¯	-f^″0
0.1	0.1	0.1	0.1	1.076865
0.2	0.1	0.1	0.1	1.038653
0.3	0.1	0.1	0.1	0.986367
0.1	0.2	0.1	0.1	1.084852
0.1	0.4	0.1	0.1	1.106473
0.1	0.6	0.1	0.1	1.129064
0.1	0.1	0.2	0.1	1.079995
0.1	0.1	0.4	0.1	1.087640
0.1	0.1	0.6	0.1	1.099853
0.1	0.1	0.1	0.2	1.346806
0.1	0.1	0.1	0.4	1.357537
0.1	0.1	0.1	0.6	1.368079

Table 5

Numerical values of via different embedded parameters.

λ	Nr	Rb	E¯	g^′0
0.1	0.1	0.1	0.1	0.106432
0.2	0.1	0.1	0.1	0.100175
0.3	0.1	0.1	0.1	0.094562
0.1	0.2	0.1	0.1	0.107527
0.1	0.4	0.1	0.1	0.109753
0.1	0.6	0.1	0.1	0.116474
0.1	0.1	0.2	0.1	0.106436
0.1	0.1	0.4	0.1	0.109292
0.1	0.1	0.6	0.1	0.116468
0.1	0.1	0.1	0.2	0.340543
0.1	0.1	0.1	0.4	0.335789
0.1	0.1	0.1	0.6	0.329064

Table 6

Numerical values of via different embedded parameters.

K	M	Rd	S	δ	Ec	θ^′0
0.1	0.1	0.1	0.1	0.1	0.1	1.225757
0.2	0.1	0.1	0.1	0.1	0.1	1.236580
0.3	0.1	0.1	0.1	0.1	0.1	1.246725
0.1	0.2	0.1	0.1	0.1	0.1	1.358975
0.1	0.4	0.1	0.1	0.1	0.1	1.347853
0.1	0.6	0.1	0.1	0.1	0.1	1.339841
0.1	0.1	0.2	0.1	0.1	0.1	1.546853
0.1	0.1	0.4	0.1	0.1	0.1	1.580637
0.1	0.1	0.6	0.1	0.1	0.1	1.625721
0.1	0.1	0.1	0.1	0.1	0.1	1.431547
0.1	0.1	0.1	0.2	0.1	0.1	1.414789
0.1	0.1	0.1	0.3	0.1	0.1	1.397825
0.1	0.1	0.1	0.1	0.1	0.1	1.653368
0.1	0.1	0.1	0.1	0.2	0.1	1.636795
0.1	0.1	0.1	0.1	0.3	0.1	1.613404
0.1	0.1	0.1	0.1	0.1	0.1	1.357437
0.1	0.1	0.1	0.1	0.1	0.2	1.347533
0.1	0.1	0.1	0.1	0.1	0.3	1.335367

Table 7

Numerical values of via different embedded parameters.

Sc	Q	Ec	ϕ^′0
0.1	0.1	0.1	0.357446
0.2	0.1	0.1	0.359881
0.3	0.1	0.1	0.360563
0.1	0.1	0.1	0.247463
0.1	0.2	0.1	0.233853
0.1	0.3	0.1	0.225164
0.1	0.1	0.1	0.115788
0.1	0.1	0.2	0.125367
0.1	0.1	0.3	0.136352

Table 8

Numerical values of via different embedded parameters.

Lb	Pe	B	Ω	χ^′0
0.1	0.1	0.1	0.1	0.368757
0.2	0.1	0.1	0.1	0.379642
0.3	0.1	0.1	0.1	0.386325
0.1	0.1	0.1	0.1	0.224678
0.1	0.2	0.1	0.1	0.214795
0.1	0.3	0.1	0.1	0.204526
0.1	0.1	0.1	0.1	0.904324
0.1	0.1	0.2	0.1	0.896537
0.1	0.1	0.3	0.1	0.885735
0.1	0.1	0.1	0.1	0.886474
0.1	0.1	0.1	0.2	0.868522
0.1	0.1	0.1	0.3	0.846723

Conclusion

This analysis has been performed for electrically conducting MHD binary fluid containing nanoparticles and gyrotactic microorganisms through a stratified stretching sheet. The present analysis has been performed in the presence of electric and magnetic fields. An analytical scheme called HAM has been dedicated for the solution of the flow problem. The present analysis is compared with previously published results and has found a great agreement. The final comments are listed as:

The presence and absence of an electric field has affected the binary fluid velocities due to the magnetic field.

The velocity and microrotation functions have escalated, while the temperature, concentration, and motile density functions have reduced via microrotation constant.

The buoyancy parameter and bioconvection Rayleigh number have reduced the velocity function, while a dual impression of buoyancy parameter and bioconvection Rayleigh number on microrotation function are depicted.

The thermal, mass, and density stratification parameters have reduced the temperature, concentration, and motile density functions.

The greater chemical reaction parameter and Schmidt number have reduced the concentration function while an opposite behavior is observed via activation energy parameter.

The higher bioconvection Peclet and Lewis number have reduced the motile density function.

The increasing buoyancy ratio parameter, bioconvection Rayleigh number, and electric parameter have increased the skin friction, whereas the increasing mixed convection parameter have reduced the skin friction.

The increasing buoyancy ratio parameter and bioconvection Rayleigh number have increased the coupe stress, whereas the increasing mixed convection parameter and electric parameter have reduced the couple stress.