Boger nanofluid: significance of Coriolis and Lorentz forces on dynamics of rotating fluid subject to suction/injection via finite element simulation.

This study briefings the roles of Coriolis, and Lorentz forces on the dynamics of rotating nanofluids flow toward a continuously stretching sheet. The nanoparticles are incorporated because of their unusual qualities like upgrade the thermal transportation, which are very important in heat exchangers, modern nanotechnology, electronics, and material sciences. The primary goal of this study is to improve heat transportation. Appropriate similarity transformations are applied for the principal PDEs to transform into nonlinear dimensionless PDEs. A widely recognized Numerical scheme known as the Finite Element Method is employed to solve the resultant convective boundary layer balances. Higher input in the solvent fraction parameter has a rising effect on the primary velocity and secondary velocity magnitude, and decreasing impact on the distributions of temperature. It is seen that growing contributions of the Coriolis, and Lorentz forces cause to moderate the primary and secondary velocities, but the temperature and concentration functions show opposite trend. The concentration, temperature, and velocities distributions for suction case is prominently than that of injection case, but inverse trend is observed for local Nusselt and Sherwood numbers. These examinations are relevant to the field of plastic films, crystal growing, paper production, heat exchanger, and bio-medicine.

Introduction

Background study

The fluids which are engineered by the homogeneous dispersion of metal and metallic oxide particles at the nanoscale are known as nanofluids[1]. The theoretical and experimental investigations confirm from the literature that the existence of nanoparticles in host fluid greatly affects the thermophysical properties of base fluids which exhibits poor conductivity properties[2,3]. The convective nanofluid heat transfer flow attracts numerous researchers due to their fascinating applications in every field of Science and Engineering. To name a few, the diamond and ceramic nanoparticles are used to upsurge the dielectric properties of the mineral oil, the fluid with nanoparticles can be used for direct absorption of the sunlight in solar collectors, the Zink and Titanium Oxide particles have antibacterial behavior and therefore which can be used for biomedical applications such as drug delivery and cancer therapy, and many more[4-6]. Ali et al.[7] discussed convective unsteady hybrid nanofluid flow within an upright channel by using Laplace Transform (LT) technique, in which authors developed fractional Maxwell fluid model is developed by using Caputo fractional differential operator. Awais et al.[8] considered nanofluid assisting and opposing flow for exploring the Lie group analysis numerically with Adams-Bash forth technique (ABT). Abdelmalek et al.[9] considered Cauchy’s stress tensor for third-grade fluid to investigate the dissipative nanofluid flow past an elongating sheet by adopting built-in Matlab routine bvp4c. Farooq et al.[10] adopted a modified Cattaneo–Christov mass and heat flu model for the numerical exploration of the Carreau nanofluid flow by horizontal stretching cylinder subject to the convective conditions. Subba Rao et al.[11] studied the nanofluid (Buongiorno model) flow above the isothermal upright cone. There are some articles related to boundary layer flow along a nanoparticles[12-14]. The non-Newtonian fluids are technological importance. In nature, Boger fluids are also non-Newtonian fluids[15,16].

There have been rising fascination in the investigation of the impact of Lorentz force produced due to magnetic field because of its modern applications in fluid engineering, assembling of plastic materials, hot rolling, production of glass fiber, and medical treatment[17,18]. Magnetic doubly stratified mixed convective flow with heat generation/absorption effects reported by Abbasi et al.[19] and Alhussain et al.[20] securitized numerically on Cattaneo-Christov flux model for the magneto nanofluid flow past a spinning cone embedded in an anisotropic permeable medium in the presence of cross-diffusion (Soret and Dufour), Navier slip, and Stefan blowing effects. The insight into significance of particles aggregation and carbon nanotube subject to rotating environment studied by Acharya et al.[21,22]. Jawad et al.[23] obtained approximate analytical solutions using HAM for the energy irreversibility analysis on hydromagnetic nanofluid flow past a nonlinear elongating permeable sheet. Jamshed[24] employed Keller Box Method (KBM) to study the MHD nanofluid flow past a nonlinearly stretching sheet with entropy generation and viscous dissipation effects.

In the recent years, the growing interest in the examination of Coriolis force effects on the dynamics of fluid flow over a stretching sheet because of several real applications in astrophysical and geophysical problems, hot rolling, fiber production, and centrifugal bio-reactor[25]. Chu et al.[26] explored numerically on nanofluid flow past a bidirectional periodically moving surface in the presence of nonlinear radiation and heat source/sink effects by employing (HAM) Homotopy Analysis Method. Ali et al.[18] employed GFEM (Galerkin Finite Element Method) to explore the buoyant driven transient bio convective Maxwell nanofluid rotating 3D flow above the Riga plate with binary chemical reaction and activation energy. Ahmed et al.[27] numerically investigated heat source/sink effects on stagnation point flow over stretching/shrinking rotating disk, Bilal et al.[28] studied chemical reaction effects on magnetized nanofluid past a rotating cone with HAM, and khan et al.[29] reported on the impact of heat source/sink effects subject to revised mass flu conditions for flow past rotating elongating cylinder. Ali et al.[30] employed FEM to explore rotating Casson Carreau nanofluid flow over a sheet with Darcy–Forchheimer medium.

The scrupulous literature assessment cited above confirms that a little attention is made towards the Lorentz force embedded in Bogar nanofluid rotating flow past stretching surface under the influence of the Coriolis force and suction/injection, and present analysis is unique because of solvent fraction, Lorentz and Coriolis impact on dynamic of Bogar nanofluid. To the best of the author’s knowledge, all the cited reports are not deliberated on the elaborated problem. Therefore, the prime aim of the current investigation is to investigate the heat and mass transfer effects on transient hydromagnetic rotating nanofluid 3D flow with suction/injection effects. Very recently, the hydromagnetic boundary layer flow for both Newtonian and non-Newtonian nanofluid flow fields was investigated by numerous authors[31-33] by employing GFEM with a weighted residual approach. The flow governing coupled highly nonlinear PDEs are solved by adopting a control volume approach by employing a Galerkin Finite Element Method[34,35] with a weighted residual approach. The various flow field characteristics for numerous significant parameter changes are discussed and depicted pictorially. The computational results obtained through Matlab code blocks are corroborated with the existing literature and found to be a tolerable correlation. The friction factor, Nusselt number, and Sherwood number values are computed and documented in tabular form. This numerical study is relevant to the polymers, heat exchangers, bio-medicine, and astrophysical and geophysical problems.

Research questions

This study provide answers to below related scientific research questions:

To observe the influence of Coriolis and Lorentz forces on the Boger fluid dynamics of temperature, fluid velocity, and concentration distributions with suction/injection?

What is the effect of combine affect of Brownian motion, thermophoresis, Coriolis and Lorentz forces on the mass transfer rate, skin friction tensors, and heat transfer rate with suction/injection?

Assess the Brownian motion, and thermophoresis parameters effect the temperature and tiny particles concentration profile?

Geometry and mathematical formulation

The transient magnetohydrodynamic 3D rotating flow of nanofluid past a bidirectional stretching sheet is considered. The physical flow configuration and coordinate system of elaborated problem is shown in Fig. 1, and flow is confined to . The nanofluid flow rotating about the z-direction with angular constant velocity . At , the surface is elongated in the direction of x along with velocity . In axial direction (z-direction), a static and uniform field of magnetic is applied. Less magnetic filed is induced due to low magnetic Reynolds number, so that Hall current and Ohmic dissipation are insignificant[30]. The ambient concentration and temperature is and , while surface concentration and temperature is represented by , and , respectively. The field of velocity for the current elaborate problem is considered by . In view of above assumptions and boundary layer approximation, the equations of conservation of mass, linear-momentums, temperature, and concentration are expressed as[36,37]:

Figure 1

Physical configuration and coordinate system.

In the above equations, C and T) are the tiny particles concentration, and temperature of fluid, () are the thermophoresis, and Brownian diffusion coefficients, are respectively the thermal diffusivity, density, dynamic viscosity of fluid, solvent fraction parameter, and ratio of relaxation time, respectively. The boundary conditions are[38,39]:

The complexity of elaborated problem is reduced by introducing the below similarity transformations (see Refs.[36,38]):In view of Eq. (10), the Eq. (1) is valid and the Eqs. (2–9) becomes below mention non-linear PDEs into transformed coordinate systems ().where primes () denote the d.r.t . The rotating , magnetic (M), Prandtl number (Pr), Lewis number (Le), Brownian motion , thermophoresis , and suction/injection () parameters are defined as:

Expressions of local Nusselt number, Sherwood number, and coefficients of skin friction are defined as:

Here, tensor of skin friction at wall are (x-direction) and (y-direction), the mass, and heat flux at the sheet are , and . Using Eq. (10):

Numerical solution

Finite element method

This numerical technique is very powerful tool reported by the experts to solve the modern engineering and industrial problems due to its accuracy and computability[18,40]. In this technique, below steps are involved:

Finite-element discretisation

First, the whole computational domain is distributed into sub-interval of a finite number and these sub-intervals are called elements, and set of these elements represent the finite element mesh.

Generation of the element equations

On the base of finite element mesh, a variational formulation of the problem is performed.

In above (i), the element equations are generated with the aid of approximation solution.

Construct a stiffness matrix by using interpolation functions.

Assembly of element equations

After generation of the elements equation, we assembled all the generated elements equations through inter element continuity condition, and obtained a model which called global finite element model over a entire computational domain.

Imposition of boundary conditions

The elaborated boundary conditions of the current problem are utilized on the obtained assembled equations.

Solution of assembled equations

Finally, we solved the assembled equations iteratively. To justify the validity of finite element computations, a study for grid independence is conducted. The whole domain is divided into different mesh density of grids size, and observe no more variation after , so we fixed the all computations on grid size (see Table 1). The comparison with earlier studies are documented in Tables 2 and 3 in particular cases are carried out to get the conformism of the precision of the solution methodology. It is witnessed that the current numerical computations are in close agreement with the existing literature in some limiting cases. The Finite Element computations for friction factor along with axial and transverse directions respectively tabulated in Table 2 for different values of rotation parameter at . It is observed from Table 2 that the numerical results so obtained are in good correlation with the results reported by Ali et al.[41], and Wang[42]. Further, in Table 3, the Nusselt number values are corroborated between Ali et al.[43] and Shafique et al.[44], and present FEM results for diverse values , and found that they are in good agreement. Hence, the confidence in the numerical computations and confirm that the Finite Element Computations obtained through Matlab code are attained the good convergence rate.

Table 1

Study of Grid independence for different grid sizes at .

Grid size	-F^′′(ζ,0)	-G^′(ζ,0)	-θ^′(ζ,0)	-Φ^′(ζ,0)
10 × 10	2.7911	1.9890	0.7145	3.2679
40 × 40	2.7611	1.0827	0.6711	3.1553
75 × 75	2.7553	1.0747	0.6865	3.0912
100 × 100	2.7550	1.0743	0.6861	3.0908
120 × 120	2.7549	1.0740	0.6860	3.0905

Table 2

Comparison of skin friction and for various values of at when ignore other involved parameters.

λ	Ali et al.[41]		Wang[42]		Present
	-F^′′(0)	-G^′(0)	-F^′′(0)	-G^′(0)	-F^′′(0)	-G^′(0)
0.0	1.00000	0.00000	1.0000	0.0000	1.00000	0.00000
1.0	1.32501	0.83715	1.3250	0.8371	1.32501	0.83715
2.0	1.65232	1.28732	1.6523	1.2873	1.65232	1.28732
5.0	2.39026	2.15024	–	–	2.39026	2.15024

Table 3

Comparison of Nusselt number for and different values of rotating parameters when other parameters are zeros.

λ	β	Pr	Ali et al.[43]	Shafique et al.[44]	FEM (our results)
0.2	0.2	1.0	0.546683	0.54670	0.5466828
–	0.4	–	0.528090	0.52809	0.5280903
–	0.6	–	–	0.51009	0.5100870
–	0.8	–	0.492547	0.49255	0.4925468

To solve the Eqs. (11) to (14) along boundary conditions (15), first consider:The set of Eqs. (11)–(15) reduced to lesser order:

For numerical computation, the plate length fixed at and thickness at . The variational form connected along with Eqs. (18)–(22) can be written as:

Here , , , for trial functions. Divide the domain () into 4-nodded elements (see Fig. 2). The related approximations of finite element are:

Figure 2

Rectangular grid.

Here, , , , and are the linear interpolation functions for are given by:

The develop finite-elements model of the equations is below:where and (m,n=1,2,3,4) are defined as:andwhere, , , , , and supposed to be the known values. Evaluate 6 functions at each node, acquired 61206 equations after assembly are nonlinear, so linearize through an iterative scheme with required precision.

Results and discussion

In this section, we show the significance of suction/injection on dynamics of Maxwell rotating fluid subject to Coriolis force, magnetohydrodynamic, and gyrotactic microorganism via finite element simulation. In each of the figures for these quantities, three sets of curves are drawn for varying values of suction/injection parameter: (i) (suction), (ii) (static), and (iii) (for injection and Table 1). The default values of the parameters in the study were: , , , , , .

The primary and secondary velocity distributions for different values of solvent fraction parameter, magnetic field parameter, rotation parameter, and unsteady parameter are portrayed in Figs. 3, 4, 5 and 6, respectively. The velocity distributions for various parameters is portrayed in Fig. 3a,b and figure exhibits that the velocity distributions is upgraded by solvent fraction parameter . Actually, at low Reynolds number, the flow rate of Bogar nanofluid increases so at constant viscosity, Bogar fluid is highly elastic. Figure 4a,b displays that and for different values of magnetic field parameter. The existence of the resistive force in the form of Lorentz force is due to the inclusion of enhancing external magnetic field and leads to deceleration of the axial momentum in Fig. 4a while an opposite behavior is perceived for transverse momentum in Fig. 4b. The axial momentum and Transverse momentum for several values of rotation parameter portrayed in Fig. 5a,b respectively. It is witnessed from Fig. 5a that the declining of the axial momentum for growing values of rotation parameter due to Coriolis force whereas an opposed action is alleged for transverse momentum in Fig. 5b. For growing values of time-dependent parameter, the increase in the thickness of momentum boundary layer in the transverse direction and reduction of the thickness of the momentum boundary layer in the axial direction is depicted in Fig. 6a,b respectively. Thus, it confirms that the unsteadiness parameter is vital in controlling the momentum in the axial direction. Moreover from these graphs, it is perceived that the reduces against the strength of (injection), but it is slightly boosted up directly with (suction) parameter.

Figure 3

Impact of on in x-direction, and in y-direction at .

Figure 4

Impact of M on in x-direction, and in y-direction at .

Figure 5

Impact of on in x-direction, and in y-direction at .

Figure 6

Impact of on in x-direction, and in y-direction at .

The sketches of (friction factor) along axial and (transverse directions) near the surface for the range of is represented in Fig. 7a,b for . Figure 7a shows that for increasing , the distribution is increased gradually up to a constant rate and then after no significant change is observed. Whereas for enhancing , a significant reduction in near the sheet surface is observed. For increasing , the distribution of is reduced gradually up to a constant rate and then after no significant change is observed as shown Fig. 7b, while enhancing . The significant difference in the values near the sheet surface is observed. Moreover from these graphs, it is perceived that the values of and for the case of (injection) are lesser than the case of (suction).

Figure 7

Impact of on in x-direction, and in y-direction.

The distribution for various parameters is portrayed in Figs. 8, 9 and 10. The influence of effect of solvent fraction parameter and time-dependent on temperature is shown in Fig. 8a,b, respectively. The profiles is decreased along with , and increased against growing strength of unsteady parameter. Figure 9 exhibits that the temperature distributions is upgraded by magnetic field parameter. The net resulting net force usually known as resistive Lorentz force between the external magnetic field and internal electric field controls the flow momentum, which is depicted in Fig. 9a, while the thickness of the thermal boundary layer is increased for enhancing as shown in Fig. 9b. The impact of thermophoresis and Brownian motion parameters on temperature is shown in Fig. 10a,b respectively. The distribution of temperature profiles is increased along with and . Further, from these graphs, it is perceived that the reduces against the strength of (injection), but it is boosted up directly with (suction) parameter. The sketches of local Nusselt number at is represented in Fig. 11a,b for respectively. For increasing M and , the distribution of is decreased gradually. A significant reduction in near the sheet surface is observed for growing values of . It is also perceived from the figure that the for the case of (injection) has a larger amount in .

Figure 8

Impact of and on at .

Figure 9

Impact of M and on at .

Figure 10

Impact of and on at .

Figure 11

Fluctuation of against , , M, and .

Figure 12a,b displays that for different values of magnetic field M, and rotation parameters. The concentration profiles are enhanced for increasing values of the magnetic field and rotating parameters as shown in Fig. 12a,b. From Fig. 13a,b revealed that the distribution profile is decrease again growing input of , and it notably rise against the . It is reported notably fall against the Lewis number Le, and time dependent parameter parameters (see Fig. 14a,b). The behavior of local Sherwood number is portrayed in Fig. 15a,b for respectively at . For increasing M and , the distribution of is decreased. Whereas for increasing , an opposite trend is observed and it is also witnessed that the (injection) case has larger than that of (suction) case.

Figure 12

Impact of M and on at .

Figure 13

Impact of and on at .

Figure 14

Impact of Le and on at .

Figure 15

Fluctuation of against , , M, and .

Conclusions

In this report, the finite element analysis on the transient magnetohydrodynamic three-dimensional rotating flow of nanofluid flow past a bidirectional stretching sheet with suction/injection, Coriolis, and Lorentz forces has been explored numerically. Based on the observed analysis of outcomes, it is worth concluding that: After successful computational struggle of parametric influence on the fluid dynamics, this study can be extend in future for Maxwell-nanofluid, Oldroyd-B nanofluid, comparative study between tangent-hyperbolic and Maxwell nanofluid, and viscoelastic Jeffrey’s nanofluids.

Higher input in the Coriolis and Lorentz has a decreasing effect on the primary velocity and secondary velocity magnitude, and

an increasing impact on the distributions of temperature and concentration.

increase the magnitude of (skin friction factor).

the velocity, temperature, and concentrations components reduce against the strength of injection.

Higher input in the solvent fraction parameter has a rising effect on the primary velocity and secondary velocity magnitude, and

decreasing impact on the distributions of temperature.

decrease the magnitude of (skin friction factor).

Simultaneous increase in Brownian motion and thermophoresis parameters have an increasing impact on the distribution of temperature, and

a negative effects on and positive effects on .

injection has a larger amount in .

injection case has larger and than that of suction case.

Tiny particles concentration diminishes against incremented Lewis number, unsteady, and Brownian motion parameters and exhibits rise against thermophoresis.