Abstract

In this study, an attempt is made to explore the two-phase Casson nanofluid passing through a stretching sheet along a permeable surface with the effects of chemical reactions and gyrotactic microorganisms. By utilizing the strength of similarity transforms the governing PDEs are transformed into set of ODEs. The resulting equations are handled by using a proficient numerical scheme known as the shooting technique. Authenticity of numerical outcomes is established by comparing the achieved results with the MATLAB built-in solver bvp4c. The numerical outcomes for the reduced Nusselt number and Sherwood number are exhibited in the tabular form, while the variations of some crucial physical parameters on the velocity, temperature, and concentration profiles are demonstrated graphically. It is observed that Local Nusselt number rises with the enhancement in the magnetic field parameter, the porous media parameter, and the chemical reactions, while magnetic field parameter along with porous media parameter retards the velocity profile.

1. Introduction

The study of two-phase nanofluids has taken considerable hype due to its vast applications in the domain of science and technology, industrial studies, engineering, and biomedical engineering [1]. Due to utmost importance of the Casson nanofluid problems, several methodologies have been employed to enhance the thermal appearances of the regular fluids. Das [2] used the convective slip, for analyzing heat properties of nanofluids fluid through a porous stretched surface. Khan et al. [3] discovered the dynamics of a nanofluid through boundary-layer region over a stretching surface. Gorla and Chamkha [4] solved the nanofluid flow problem over a horizontal plate using a permeable surface along natural convective boundary-layer conditions through a horizontal plate. Noghrehabadi et al. [5] explored the physics of heat transfer by employing slip in the boundary layer region for the flow of a nanofluid together with dynamic nanoparticles. In another study Noghrehabadi et al. [6] observed decline in the Nusselt number when slip factor is enhanced. Ibrahim and Shanker [7] investigated the flow of nanofluid in boundary layer region along with heat transfer effect past a porous stretching surface by incorporating magnetic field properties, slip at boundary, and thermal radiation. Goyal and Bhargava [8] reported the heat transfer flow of viscoelastic nanofluid with velocity slip at boundary past a stretched surface. Recently, nanofluid flow along stretching sheet with Lorentz’s impact has taken huge inspiration for the researchers’ community because of its significant involvement in MHD generators, nuclear reactor studies, plasma physics, wire drawing, manufacturing of glass fiber and hot rolling, etc.

Recently the flow involving heat and mass conversion investigations on fluids along with chemical effects past through stretching surface has proven significant importance in metallurgy and industrial studies of chemical engineering like polymer production and food processing. Some possible applications are the study of drying, temperature and moisture distribution, decline of temperature in a cooling tower, and desert cooler flow [9]. The traditional fluids, namely, water, oil, and ethylene glycol compound, are considered very weak fluids for the transfer of heat. The heat involvement through the surface and the medium is strongly affected by the thermal conductivity. Sheikholeslami et al. [10] used the radiation effects for free convective heat transfer and they recovered that the intensity of heat transfer reduces due to the Lorentz forces. Garoosi et al. [11] reported the combined convection and naturally convective heat transfer between heated cylinders of differential length filled with nanofluid. Turkyilmazoglu [12] investigated mathematically MHD fluid flow through porous surface with heat involvement persuaded by stretching or shrinking sheet for the two and three dimensional deformed objects. Boungiorno [13] worked on convective dynamics of nanofluids by taking Brownian diffusion and established the unnatural convective heat transfer viewed in nanofluids. He also explained the importance of two nanoparticles, Brownian diffusion and thermophoresis. Congedo et al. [14] with the aid of finite volume method calculated the numerical results for naturally convective heat transfer nanofluids flow. Rashidi et al. [15] made investigations to explore the flow of nanofluid within permeable rotating disk in the presence of magnetic field by employing the law of thermodynamics. They also explained the importance of the magnetic disk, in the application of heat expansion development for long lasting energy systems.

The problems of porous medium have attained much consideration of the researchers due to its vast utilities in heat exchange, water-flow in reservoirs, rock mechanics, acoustics study, soil mechanics, production of oil, and catalytic reactors. The porous media play a crucial role in increasing the overall efficiency of the solar systems, improving the heating process of green houses and in storing the solar heat in packed beds. Vafai [16] investigated the forced convection flow in a permeable channel and achieved the expressions for velocity and temperature profiles with the help of exact solution. They used the porous medium properties and considered the fluid to be homogeneous and isotropic. Nakayama et al. [17] used the non-Darcy forced convective heat transfer in a medium restricted by channel boundaries. Nield et al. [18] presented an analytical technique to formulate the forced convection in a plane working on a saturated bidispersed permeable channel. They observed that with the finite thermal resistance effects reduction happens in the intensity of unstable local thermal and transfer of heat to the porous medium. Bakier [19] used the Runge-Kutta method while studying the impacts of thermal radiation on the mixed convection in a permeable geometry.

The bioconvection flow is the large scale convection of fluid flow produced by the density gradient formed by swimming of dynamic microorganisms. Oberbeck-Boussinesq calculation with aid of the similarity transformation for MHD boundary layer flow with transfer of heat and mass in a water-oriented nanofluid having gyrotactic microorganisms passing through a transverse plate was reported by Khan et al. [20]. Xu et al. [21] used the mixed convection flow passing on a stretching sheet with both gyrotactic microorganisms and nanoparticles. They found that the thermophoresis parameter rises as the Nusselt number declines. Kuznetsov [22] presented the commencement of bioconvection having a nanofluid and gyrotactic microorganisms. To find an analytical solution, Galerkin and a linear instability method were used for the nonoscillatory condition. Kuznetsov and Nield [23] proposed the model of mixed bioconvection and later on this model was studied by Xu and Pop [24]. Mahdy [25] discussed a transverse cone in a permeable media, when bioconvection in a peristaltic flow takes place in an asymmetric channel that is comprised of nanofluid enclosing the gyrotactic microorganisms. Akbar [26] reported the bio-nano-engineering model by using bioconvection peristaltic flow in a channel of asymmetric geometry using nanofluid comprising gyrotactic microorganisms.

The electrically conducting two-dimensional steady incompressible viscous nanofluid flow passing a stretching sheet has diverse applications in the industry. Some of the key applications of viscous flow through a stretching sheet [2731] can be narrated as aerodynamic extrusion of plastic sheets, condensation mechanism of metallic plate through a cooling bath, and manufacturing of a polymer sheet through a dye, processing of food stuffs, glass fiber, metallic packaging, rope rolling, and aluminum jug industrial advancements, etc.

The prime features of current study are described as follows:(i)A novel mathematical model is developed for two-phase Casson nanofluid through a stretching sheet on a permeable surface with the variations of chemical reactions and gyrotactic microorganisms.(ii)By utilizing the strength of similarity transforms the resulting PDEs are changed into set of ODEs.(iii)The resulting ODEs are tackled by using a well established proficient numerical scheme known as the shooting technique.(iv)Numerical results are validated by comparing the achieved results with the MATLAB built-in solver bvp4c.(v)Numerical outcomes for the reduced Nusselt number and Sherwood number are exhibited in the tabular form while the impacts of crucial parameters are investigated graphically.

2. Mathematical Interpretation of the Problem

A novel electrically conducting two-dimensional steady incompressible viscous nanofluid flow is investigated through a stretching surface having plate temperature and species concentration , respectively. The stretching sheet is taken along -axis in the direction of the motion and y-axis is taken in the transverse direction of the plate. The elementary equations for momentum, energy, concentration, and gyrotactic microorganism [32, 33] are described as in [32] while the coordinate setting and physical configuration are exhibited in Figure 1.The suggested boundary conditions are written asThe nonuniform permeability of the source is taken as whereas the dimensional heat and chemical reaction parameters were taken as and , respectively. The similarity transformations are given asUsing the above transformations, the ordinary differential equations (ODEs) are given below:The boundary conditions may be written in the dimensionless form asIn the above equations, different parameters have the following formulations:

3. Solution Methodology

A brief description on the numerical approach used for the solution of the above nonlinear, nondimensional, and coupled ODEs (8)-(11) along with the boundary conditions (12) has been presented in this section. The nonlinear systems of equations are tackled with aid of well known numerical technique, namely, the shooting method. For comparison purpose, the same system is solved by using the MATLAB built-in solver bvp4c. The domain has been taken as . Due to the observed ignorable variations in the solution for , has been chosen as 10 throughout the article. To convert the system of higher order ODEs given in Eqs. (8)-(11) to set of first order ODEs, is denoted by , by , by , and by . The resulting equations are written asWith the new notations, the set of initial conditions from (12) gets the following form: By choosing the missing initial conditions and , the resulting initial value problem is obtained by the shooting method. Newton’s method has been opted to update the missing initial conditions. The numerical solution is gathered by defining following accuracy criteria:where is very small and positive number. The numerical outcomes of this work are obtained with .

A critical analysis of present investigation with already obtained result available in the literature for against the Prandtl number and Eckert number with those of Ferdows et al. [34] and Bidin and Nazar [35] is shown in Table 1. The computed result exhibited in Table 1 shows impressive match with the already available results.

4. Results and Discussions

It is quite clear from Table 2 that the amount of heat transfer reduces with the rise in magnetic parameter , nonuniform permeability , dimensionless inertia factor , mass convective parameter , Eckert number , thermophoresis factor , and Brownian motion factor ; on other hand quite different behavior is noticed for the thermal convective parameter and Prandtl parameter . However Sherwood number is found to have a decreasing impact due to the magnetic parameter , nonuniform permeability , dimensionless inertia factor , and Prandtl parameter , while an opposite behavior is seen for the thermal convective parameter , which is the mass convective parameter, Eckert number , thermophoresis factor , and Brownian motion factor .

Table 3 demonstrates the behavior of local Nusselt number values and Sherwood number values for varying various flow indicators with , .1, and =1. The rate of heat transfer increases for the heat generation/absorption and viscosity ratio factor while an opposite behavior is noticed for Lewis number , chemical reaction factor , and Casson fluid factor . Furthermore, Sherwood number increases for the Lewis number , chemical reaction factor , and viscosity ratio factor and decreases for the heat generation/absorption parameter and Casson fluid factor .

In order to have a deep insight of the solution, graphical illustrations of different profiles are presented. Figure 2 reflects the behavior of velocity profile for various values of the magnetic field parameter . The magnetic field along with the electric current generates a force which is reputed as Lorentz force creates obstacle in the fluid motion, which results in decrements in the velocity profile. Since the magnetic force is applied against the direction of the flow which raises stationary fluid absorption on the stretching sheet also retards the flow momentum. As a result, the flow becomes slow. Figure 3 speculates the conduct of the velocity profile against the porous media parameter . The porosity parameter escalates when the fluid becomes more viscous. Due to an enhancement in the viscosity, the velocity of the fluid is declined. The graphical impact of the dimensionless inertia parameter is demonstrated in Figure 4. It is quite evident that, by enhancing the inertia parameter, the thermal boundary layer thickness reduces, which creates a decrement in the velocity profile. Analysis of increasing Prandtl number value on velocity profile is exhibited in Figure 5. The Prandtl number exhibits the ratio of kinematic viscosity to thermal diffusivity [3639], and the value of Prandtl number is considered about 0.7 for gasses while for water it is about 7. As temperature is increased, the value of Prandtl number declines; therefore we have considered the values of Prandtl number between 0.1 and 1.5. In heat transfer mechanism, the Prandtl number maintains the momentum and thermal boundary layers thickness. If we take the value of Pr small, it reflects that the heat diffusivity can be compared to the momentum of the flow in a short span of time. This indicates that, for the case of liquid metals, the thermal boundary layer thickness is much higher than to the boundary layer velocity. So, by enhancing Prandtl number Pr temperature profile declines; on the other hand thermal boundary layer rises. By escalating the thermal diffusivity values, the boundary layer thickness declines because reduction is seen in the Prandtl number. Figure 6 indicates the variations of the thermal convective parameter on the velocity. When increase in the thermal convective parameter values is made, a progressive increment in the flow velocity is seen. The impact of mass convective parameter on the velocity profile is demonstrated in Figure 7. The principle factor is the Buoyancy forces that are responsible to reduce the velocity profile. By enhancing the value of , the velocity increases within the boundary layer which brings a decrement in the velocity profile. Figure 8 reflects the impact of Brownian motion parameter on the temperature profile. To enhance the values, the intermolecular collision between the atoms is increased, which escalates the temperature of the fluid at the end. The variations of the thermophoresis parameter on the temperature profile plots are highlighted in Figure 9. By escalating the thermophoresis parameter, the molecules move from hot surface to the cold one, which results in an enhancement in the thermal layer boundary, which decreases the temperature profile. The graphical effects of Lewis number on the concentration profile are exhibited in Figure 10. Chen et al. [40] experimentally used the values of Lewis number between 0 and 2. Various researchers [27, 4042] have studied stretching sheet phenomena theoretically and experimentally and used similar ranges in their studies. Decrease in the velocity profile is recorded close to the boundary wall and an enhancement is seen by moving away from the boundary. The evidence behind this process is the occurrence of the boundary layer thickness that increases the nanoparticle volume fraction and its boundary layer. Figure 11 highlights the impacts of the chemical reaction parameter on the nanoparticle concentration profile. It depicts that the concentration profile decreases near the boundary wall and the concentration of the species escalates by moving away from the boundary wall. The effects of the thermophoresis parameter on the nanoparticles concentration profile are depicted in Figure 12. It is concluded that huge gap in temperature drives the fluid more quickly from hot surface to the cold. And finally, the depth of the concentration boundary layer increases which ultimately augments the concentration profile. Figure 13 represents the behavior of the mass friction field versus the Brownian motion parameter . By enhancing the Brownian motion parameter, the thickness of thermal boundary layer increases due to this concentration profile. Figure 14 displays the behavior of the bioconvection Lewis number parameter against the density of motile microorganisms profile. An increment in the value of leads to a decrement of the microorganisms diffusion. Ultimately, thickness of the boundary layer along with density of the motile microorganisms declines.

5. Conclusion

The current investigation is designed to solve the two-phase Casson nanofluid passing over a stretching sheet using the permeable surface along with the effects of the chemical reactions and gyrotactic microorganisms. The systems of nonlinear PDEs are transformed to ODEs by employing the similarity transformations. The well known shooting numerical technique along with the fourth order Runge-Kutta numerical approach is deployed to solve these ODEs. The software used for solving this system of nonlinear equations is MATLAB R2016(a) package. Some novel findings of this investigation are presented as follows:(i)It is evident from this study that Local Nusselt number rises with the enhancement in magnetic field parameter, the porous media parameter, and the chemical reactions.(ii)The magnetic field parameter along with porous media parameter shows a decreasing effect on the velocity profile .(iii)The inertia parameter and the Prandtl number demonstrate a decreasing effect on the velocity profile .

Nomenclature

:Uniform magnetic field
:Velocity components
:Density of the fluid
:Kinematic viscosity
:Reference kinematic viscosity
:Uneven permeability
:Heat at fixed pressure
:Temperature of the fluid
:Concentration of the fluid
:Inertia value
:Brownian coefficients
:Thermophoresis diffusion coefficient
:Electrical conductivity
:Reference velocity
:Reference temperature
:Concentration
:Length
:Gyrotactic microorganisms
:Heat generation coefficients
:Chemical reaction coefficient
:Casson fluid factor
:Viscosity ratio factor
:Inertia factor
:Thermal convective factor
:Mass convective factor
:Bioconvection Peclet number
:Magnetic field factor
:Prandtl number
:Brownian motion factor
:Thermophoresis factor
:Lewis number
:Heat generation/absorption factor
:Chemical reaction factor
:Eckert number
:Bioconvection Lewis number
:Microorganism concentration difference parameter.

Data Availability

All the data used in this research work is available within the manuscript.

Conflicts of Interest

The authors declare that they have no conflicts of interest.