Next Article in Journal
Entropy Generation in a Mass-Spring-Damper System Using a Conformable Model
Next Article in Special Issue
Application of Fractional Derivative Without Singular and Local Kernel to Enhanced Heat Transfer in CNTs Nanofluid Over an Inclined Plate
Previous Article in Journal
Influence of Ramped Wall Temperature and Ramped Wall Velocity on Unsteady Magnetohydrodynamic Convective Maxwell Fluid Flow
Previous Article in Special Issue
Utilization of Second Order Slip, Activation Energy and Viscous Dissipation Consequences in Thermally Developed Flow of Third Grade Nanofluid with Gyrotactic Microorganisms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Numerical Exploration of Modified Second-Grade Nanofluid with Motile Microorganisms, Thermal Radiation, and Wu’s Slip

1
School of Mathematics, Xijing University, Xi’an 710123, China
2
Department of Mathematics, Government College University Faisalabad, Punjab 38000, Pakistan
3
Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia
4
Department for Management of Science and Technology Development, Ton Duc Thang University, Ho Chi Minh City 758307, Vietnam
5
Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City 758307, Vietnam
*
Author to whom correspondence should be addressed.
Submission received: 6 January 2020 / Revised: 19 February 2020 / Accepted: 21 February 2020 / Published: 3 March 2020
(This article belongs to the Special Issue Future and Prospects in Non-Newtonian and Nanofluids)

Abstract

:
This study is carried out to scrutinize the gyrotactic bioconvection effects on modified second-grade nanofluid with motile microorganisms and Wu’s slip (second-order slip) features. The activation energy and thermal radiation are also incorporated. The suspended nanoparticles in a host fluid are practically utilized in numerous technological and industrial products such as metallic strips, energy enhancement, production processes, automobile engines, laptops, and accessories. Nanoparticles with high thermal characteristics and low volume fraction may improve the thermal performance of the base fluid. By employing the appropriate self-similar transformations, the governing set of partial differential equations (PDEs) are reduced into the ordinary differential equations (ODEs). A zero mass flux boundary condition is proposed for nanoparticle diffusion. Then, the transmuted set of ODEs is solved numerically with the help of the well-known shooting technique. The numerical and graphical illustrations are developed by using a collocation finite difference scheme and three-stage Lobatto III as the built-in function of the bvp4c solver via MATLAB. Behaviors of the different proficient physical parameters on the velocity field, temperature distribution, volumetric nanoparticles concentration profile, and the density of motile microorganism field are deliberated numerically as well as graphically.

1. Introduction

In recent years, the tremendous progress in research and wide-ranging solicitations of functional nanoparticles have been noticed. Nanoparticles have numerous potential applications in the many fields such as the biomedical field, microbiology, supramolecular and colloid chemistry, material sciences, petroleum sciences, social sciences, etc. Nanoparticles have further physical features that must be determined for an undivided depiction, for example surface characteristics, sizes, and shapes. Nanoparticles belong to a broad interdisciplinary field in pharmaceutical medicines, thermal systems, electronics, nuclear reactors, the chemical industry, etc. Nanoparticles are also expressed in ink, polymer, and natural fluids such as blood. In general, nanoparticles-based technologies are focused on opportunities for improving the efficiency, sustainability, and speed of already existing processes. A nanofluid is an advanced type of fluid that contains nanometer-sized particles. As a result of their sub-nanoscopic size, they have unique material features and manufactured nanomaterials that may find practical applications in a mixture of areas, including catalysis, environmental remediation, medicine, and engineering. Beneath that classification, nanosubstances needs simply one of their quality scopes to be in the range of 1–100 nm to be classified as nanoparticles, even if its other dimensions are exterior to that range. Choi [1] proposed the idea of nanofluids. Numerous researchers have gained interest in the study of nanofluids due to their significant technical processes and applications. The coefficient of heat transfer between the medium and the soil is considered to be important. Nanofluids are commonly used in automotive and industrial cooling, tumor care, tracking, cooling towers, new fuel types, micro electron cooling, imaging, the heating/cooling of home appliances, and many more. Buongiorno [2] presented a two-phase model depicting the salient features of Brownian diffusion and thermophoresis. He determined that the properties of thermophoretic and Brownian motion are enough to establish a viable model. Khan et al. [3] discussed the characteristics of Brownian motion and thermophoresis in Jeffery nanofluid flow by a stretched surface. The article reveals the impacts of drag force and heat transfer rate for numerous physical variables. Nanofluid heat transfer in a porous duct in the presence of Lorentz forces by using the lattice Boltzmann method was examined by Li et al. [4]. Several researchers made serious attempts on nanofluids subjected to different configurations [5,6,7,8,9,10,11,12]. Alamri et al. [13] described the impacts of a second-order slip on plane Poiseuille nanofluid flow under the influence of Stefan blowing in a channel. Kumar [14] explored the behavior of heat transmit flow of non-Newtonian nanofluid through a stretching sheet. Farahangmehr et al. [15] scrutinized numerically the continuous boundary layer of an incompressible viscous nanofluid with heat and mass transfer over the horizontal surface. Tlili et al. [16] discussed boundary layer flow in the presence of convection and mass transfer features of the Buongiorno nanofluid model over a wedge with Navier slip and Biot number. Khan et al. [17] inspected the mass and heat transfer mechanism in Carreau nanofluid flow across a wedge. The unsteady slip flow of Carreau nanofluid over a moving/ static wedge was numerically investigated by Masood et al. [18]. Khan et al. [19] examined the significance of activation energy and nonlinear thermal radiation on modified second grade fluid flow in the presence of nanoparticles. Rehman et al. [20] scrutinized the impacts of heat transfer and heterogeneous–homogeneous reactions in modified second-grade nanofluid over a stretching sheet. Khan et al. [21] investigated the multi-dimensional axisymmetric flow and convection of the modified second-grade fluid through the nonlinear radially expanding oscillatory surface. Hassan et al. [22] investigated the mass and heat transport aspects in viscoelastic nanofluid flows under the presence of velocity slip conditions. Afterwards, the methods of combining second-grade fluid and a power law model resulted in the generalized second-grade fluid as described by Man and Sun [23].
Bioconvection caused by the density differences of motile microorganisms is effectively assorted in the field of environmental systems, biofuels, and industry. The presence of these microorganisms increases the primary density of the fluid and creates a density gradient by swimming, which leads to bioconvection. This interesting observation inevitably leads to an unstable, lower density surface. Microorganisms are typically classified according to the impellent force of various types due to gyrotactic, oxytactic, and gravitactic microorganisms. There are some distinct and similar characteristics of nanoparticles and motile microorganisms. Immunology microsystems (i.e., biomaterials, tissue engineering, protein engineering, synthetic biology, and drug delivery systems) such as enzyme biosensors are typically involved in bioconvection applications. Kuznetsov [24] therefore suggested that the growth of microorganisms in the biomicrosystems plays a major role in the processing and aggregation of mass transport. Khan et al. [25] analyzed the influence of bio-convection in a rheology of magnetized nanofluid. Waqas et al. [26] investigated the heat and mass transfer phenomenon of Williamson’s nanofluid flow in the existence of gyrotactic motile microorganisms. Waqas et al. [27] considered viscous nanofluid flow across a flexible stretching disk with a mixture of microorganisms. Waqas et al. [28] also visualized the flow of modified second-grade nanofluid in the presence of motile microorganisms with heat and mass transfer rates across the stretching boundary. Waqas et al. [29] examined the magnetohydrodynamics (MHD) flow of rate-type nanofluid in the presence of gyrotactic microorganisms with activation energy. Recently, Waqas et al. [30] securitized the novel biofuel significance in viscoelastic nanofluid under the bioconvection process.
This research paper provides a more compelling approach to the nanoliquid heat and mass transfer phenomena through the bioconvection of self-motile microorganisms, which helps to avoid the agglomeration of nanoparticles. New results are established with inclusion of Wu’s slip (second-order velocity slip), activation energy, and the thermal radiation. The manuscript is organized as follows. Section 1 describes the introduction, Section 2 contains a mathematical model, Section 3 yields the numerical solution, Section 4 present the graphical analysis, and the conclusions of the graphical and tabular results are exhibited in Section 5. Observation of the results shows that the buoyancy ratio parameter R b , bioconvection Rayleigh number N r , and Wu’s slip parameter inhibit the fluid flow and cause a reduction in the velocity. However, the parameters R b and N r augmented the thermal distribution.

2. Mathematical Model

Let us consider the two-dimensional, steady-state modified second-grade nanofluid with activation energy and gyrotactic motile microorganisms. Moreover, the flow is incompressible across a stretching surface. Additionally, non-linear thermal radiation and Wu’s slip (second-order velocity slip) are taken into consideration. Let us assume that the surface stretches as U = c x ˜ α , where α > 0 is the power law exponent parameter, and c > 0 is symbolized for stretching strength. Furthermore, V = V u , v is the flow velocity, the temperature of the fluid is T ˜ , and the volume fraction for nanoparticles and the motility of the microorganisms are symbolized as C ˜ , N ˜ , respectively. The Cauchy stress tensor τ ˜ ˙ for modified second-grade fluid is taken as [31]
τ ˜ ˙ = P I + u m 2 A ˜ 1 + β 1 A ˜ 2 + β 2 A ˜ 1 .
Here, I is the identity tensor, P is the pressure, μ is viscosity coefficient, β 1 & β 2 are the conventional stress coefficients, m is the material parameter, = 1 2 t r A ˜ 1 2 is the velocity gradient for the symmetric part of second invariant, and A ˜ 1 & A ˜ 2 are termed as kinematical tensors and defined as
A ˜ 1 = L + L T ,   A ˜ 2 = d A ˜ 1 d t + A ˜ 1 L + L T A ˜ 1
where L = g r a d V .
Man and Sun [23] studied the Cauchy stress tensor and described two models as terminated in modified forms of second-grade fluid and power law fluid model, which are:
Model (a):
φ 1 = P I + μ f m 2 A ˜ 1 + β 1 A ˜ 2 + β 2 A ˜ 1 2
Model (b):
φ 1 = P I + m 2 [ μ f A ˜ 1 + β 1 A ˜ 2 + β 2 A ˜ 1 2 ] .
Model (a) is about second-grade fluid.
Model (b) is representing the power-law index model.
The boundary layer formulation for mass conservation, momentum, energy, concentrations of nanoparticles, and motile microorganisms’ distributions are [32,33,34]:
u ˜ x ˜ + v ˜ y ˜ = 0 ,
u ˜ u ˜ x ˜ + v ˜ u ˜ y ˜ = ν 1 + m 2 u ˜ y ˜ 2 u ˜ y ˜ m + α 1 ρ f v ˜ 3 u ˜ y ˜ 3 + u ˜ 3 u ˜ x ˜ y ˜ 2 + u ˜ x ˜ 2 u ˜ y ˜ 2 u ˜ y ˜ 2 u ˜ x ˜ y ˜ + 1 ρ f 1 C ˜ g 1 * ρ f β 1 * T ˜ T ˜ ρ p ρ f g 1 * C ˜ C ˜ ρ p ρ f γ 1 * g 1 * N ˜ N ˜ ,
u ˜ T ˜ x ˜ + v ˜ T ˜ y ˜ = α f + 16 σ s T ˜ 3 3 k ρ c f 2 T ˜ y ˜ 2 + ρ c p ρ c f D B C ˜ y ˜ T ˜ y ˜ + D T T ˜ T ˜ y ˜ 2 ,
u ˜ C ˜ x ˜ + v ˜ C ˜ y ˜ = D B 2 C ˜ y ˜ 2 + D T T 2 T ˜ y ˜ 2 K r 2 C ˜ C ˜ T ˜ T ˜ n exp E a κ T ˜ ,
u ˜ N ˜ x ˜ + v ˜ N ˜ y ˜ + b 1 W c C ˜ w C ˜ y ˜ N ˜ C ˜ y ˜ = D m 2 N ˜ y ˜ 2 .
The boundary conditions entertained the following relations:
u = u w + U ¯ s l i p ,   v = 0 ,   k T ˜ y = h f T ˜ f T ˜ ,   D B C ˜ y + D T T ˜ T ˜ y = 0 ,   N ˜ N ˜ w y = 0 ,
u U ˜ = 0 ,   u y 0 ,   T ˜ T ˜ ,   C ˜ C ˜ ,   N ˜ N ˜     a t     y .
The temperature of the surface is assumed as T ˜ w , while the volume fraction for nanoparticles and the density of the microorganisms at the surface are denoted as C ˜ w , N ˜ w respectively. T ˜ , C ˜ , N ˜ are the free-stream temperature, concentration of nanoparticles, and the density of the microorganisms. Some researcher used the U ˜ s l i p as [35,36]:
U ˜ s l i p = 2 3 3 ϖ l 2 ϖ 3 2 1 p 2 K n σ *   y u 1 4 p 4 + 2 K n 2 1 p 2 σ 2 y y u ,
U ˜ s l i p = σ   y u + ε 1   y y u ,
where K n notifies the Knudsen number, σ and ε 1 stand for constant numbers, and β expresses the free path associated with the molecular mean motion. In the above equations, α 1 is the material-related parameter, m represents the power law index, ρ m represents the motile microorganism particles density, ρ p represents the nanoparticles density, γ 1 * represents the average volume of the microorganism, g 1 represents the gravity, T ˜ notifies the temperature of the nanofluid, C ˜ reflects the nanoparticles’ volume resistance, α f stands for the friction coefficient, ρ c f represents the heat capability of fluids, ρ c p represents the impressive heat capacity of the nanomaterial, D B symbolizes the Brownian diffusivity, σ * represents the Stefan Boltzmann constant, k * represents the mean absorption coefficient, D T stands for the thermophoresis diffusion coefficient, E a reflects the activation energy coefficient, K r 2 reflects the reaction rate constant, b 1 represents the chemotaxis constant, W c describes the superlative cell swimming rapidity, and D m is the microorganism’s diffusion coefficient. Equations (5)–(9) are non-dimensionalized by the following variables [22]:
ψ = x ˜ U R e 1 2 + m f ( ζ ) ,   ζ = y ˜ x ˜ Re 1 2 + m ,   θ ( ζ ) = T ˜ T ˜ T ˜ w T ˜ ,   ϕ ( ζ ) = C ˜ C ˜ C ˜ w C ˜ ,   χ ( ζ ) = N ˜ N ˜ N ˜ w N ˜
Equation (5) is satisfied. The re-established dimensionless form of Equations (6)–(9) respectively yield:
1 + α + 2 m α f f + m + 2 m + 1 f f m f 2 α * 3 α 1 f 2 2 f f + 1 + α + 2 m α f i v f   + Λ 1 θ N r ϕ R b χ = 0 ,
1 + ε θ θ + 4 3 R d 1 + θ w 1 θ 3 θ + ε + 1 + θ w 1 θ 2 4 R d θ w 1 θ 2 + Pr 1 + α + 2 m α m + 2 f θ + Pr N b θ ϕ + N t θ 2 = 0 ,
ϕ + Pr L e 1 + α + 2 m α m + 2 f ϕ + N t N b θ Pr L e σ 1 + δ θ n exp E 1 + δ θ ϕ = 0 ,
χ + L b 1 + α + 2 m α m + 2 f χ P e ϕ χ + δ 1 + χ ϕ = 0 .
Similarly, the boundary conditions (10) and (11) are transformed as follows:
f ( ζ )   = 0 ,   f ( ζ )   =   1 + Γ f ζ   + β f η ,   θ ( 0 ) = B i θ ζ 1 , N b θ ( ζ ) + N t ϕ ( ζ ) = 0 ,   χ ( ζ )   = 1 a t ζ = 0 ,
f   0 , θ 0 , ϕ 0 , χ 0 ,   as   ζ ,
where α * = α 1 Re 2 2 + m ρ f x ˜ 2 , is the generalized second-grade fluid parameter, Λ 1 = g β 1 C ˜ T ˜ f T ˜ x ˜ 3 c 2 r 2 α 1 is the mixed convection parameter, N r = ρ p ρ f C ˜ w C ˜ x ˜ 3 T ˜ w T ˜ β c r 2 α 1 1 C ˜ ρ f is the buoyancy ratio parameter, R b = γ ρ m ρ f N ˜ w N ˜ x ˜ 3 ρ f 1 C ˜ c r 2 α 1 β T ˜ w T ˜ is the bioconvection Rayleigh number, Pr = U x ˜ Re 2 2 + m k 1 / ρ c p is the Prandtl number, R d = 16 σ * T 3 3 k * K is the radiation parameter, N b = ρ c p C ˜ w D B C ˜ D B ρ c f ν is the Brownian motion constraint, N t = ρ c p D T T ˜ f D T T ˜ ρ c f T ˜ ν is the thermophoresis parameter, the temperature ratio parameter is θ w = T w T , L e = ν D B is the Lewis number, δ 1 = N N w N is the motile microorganism difference parameter, E = E a κ T is the activation energy, L b = ν D m is the bioconvection Lewis number, P e = b 1 W c D m is the Peclet number, Γ = C 1 Re 1 2 + m x ˜ is the first-order velocity slip, the second-order velocity slip is represented as β = C 2 Re 2 2 + m x ˜ 2 , and the Biot number is symbolized as B i = h f k 1 x ˜ Re 1 2 + m .
In order to express the consequences of wall shear stress, we have following mathematical relations:
C f = τ w 0.5 ρ U 2 , τ w = μ u ˜ y ˜ n u ˜ y ˜ + α 1 u ˜ 2 u ˜ x ˜ y ˜ + v ˜ 2 u ˜ y ˜ 2 + 2 u ˜ x ˜ u ˜ y ˜ y ˜ = 0 ,
where τ w is the wall shear stress. By using the following dimensionless variable, Equation (21) yields:
1 2 C f x Re x 1 / m + 2 = 1 m + 2 1 7 + 2 m α α 1 f 0 + f 0 m + 1 , .
Moreover, the local Sherwood number, motile density number, and local Nusselt number are stated as follows:
N u x = x ˜ T ˜ w T ˜ k + 16 σ * T ˜ 3 3 k * T ˜ y ˜ y ˜ = 0 ,   S h = x j w D B C w C ,   j w = D B C y y = 0 , N h = x ˜ j n D B N w N ,   j n = D B N y y = 0 , .
The above quantities in dimensionless forms are:
Re 1 m + 2 N u x = 1 + 4 3 R d θ 0 , Re 1 m + 2 S h x = ϕ 0 , Re 1 m + 2 N n x = χ 0 .

3. Numerical Procedure

The system of the non-linear coupled ordinary differential Equations (15)–(18) with the boundary conditions (19) and (20) has been computed by bvp4c [37,38] in MATLAB with a tolerance level of 10 5 . In the bvp4c routine, a three-stage Lobatto-IIIA method is employed in the development of this collocation technique. The bvp4c function of computational software MATLAB only resolves the first-order ordinary differential equations. To approximate the solution of the distorted non-linear ODEs (ordinary differential equations), appropriate initial guesses satisfying the pertinent boundary conditions are required. For further approximation, the initial guess is modified with the help of shooting technique. Let us consider
f = l 1 ,   d f d ζ = l 2 ,   d 2 f d ζ 2 = l 3 ,   d 3 f d ζ 3 = l 4 ,   d 4 f d ζ 4 = l 4 ,   θ = l 5 ,   d θ d ζ = l 6 , d 2 θ d ζ 2 = l 6 ,   ϕ = l 7 ,   d ϕ d ζ = l 8 ,   d 2 ϕ d ζ 2 = l 8 ,   χ = l 9 ,   d χ d ζ = l 10 ,   d 2 χ d ζ 2 = l 10 .
Equations (15)–(18) respectively are presented below:
l 3 = α * 3 α 1 l 3 2 l 4 l 2 1 + α + 2 m α l 1 l 3 m + 2 m + 1 l 4 l 3 m l l 2 2 Λ 1 l 5 N r l 7 R b l 9 α * l 1 1 + α + 2 m α
l 6 = 1 1 + ε l 5 + 4 3 R d 1 + θ w 1 l 5 3 Pr 1 + α + 2 m α m + 2 l 1 l 6 + Pr N b l 6 l 8 + N t l 6 2 + ε + 1 + θ w 1 l 5 2 4 R d θ w 1 l 6 2
l 8 = N t N b l 6 Pr L e 1 + α + 2 m α m + 2 l 1 l 8 + L e Pr σ 1 + δ l 5 n exp E 1 + δ l 5 l 7
l 10 = P e l 8 l 11 + δ 1 + l 8 l 10 L b 1 + α + 2 m α m + 2 l 1 l 10 .
The associated boundary conditions are:
l 1 ζ = 0 ,   l 2 ζ = 1 + Γ l 3 ζ + β l 4 ,   l 6 ζ = B i l 5 ζ 1 , N b l 6 ζ + N t l 8 ζ = 0 , l 9 ζ = 1 ,   a s   ζ = 0 ,
l 2 ζ 0 ,   l 5 ζ 0 ,   l 7 ζ 0 ,   l 9 ζ 0 ,   a s   ζ .

4. Graphical Analysis

This section is set to explain the emerging consequences of the involved parameters on the physical quantities of interest. Similar to the traditional way, the involved parameters have been assigned some fixed values such as α = 0.2 , Λ 1 = 0.3 , N r = 0.2 , R b = 0.2 , R d = 0.4 , Pr = 0.7 , N t = 0.3 , N b = 0.2 , δ = 0.1 , P e = 0.2 , and L b = 2.0 , Γ = 1.0 and β = 1.0 . Furthermore, the complete analysis has been performed in view of the shear thinning case m = 0.5 , second-grade fluid m = 0.0 , and shear thickening case m = 0.5 . The variation in the impact of the buoyancy ratio parameter over velocity distribution is drawn in Figure 1a. It illustrates that the velocity component f ζ decreases on the rising variation of the buoyancy ratio parameter R b with the power law index m = 0.5 , 0   &   0.5 . The physically buoyancy ratio parameter R b is the phenomenon of upthrust of fluid to objects, which exerts pressure on it. The discussion of the buoyancy ratio parameter rose to another next level when it was considered for the temperature profile. The behavior of temperature distribution θ ζ in relation to the variation of the buoyancy ratio parameter can be traced from Figure 1b. It can be noticed that the curve of temperature profile θ ζ rises when the buoyancy ratio parameter is increased. Regarding the concentration profile, the buoyancy ratio parameter has a unique trend, as graphed by Figure 1c. It can be viewed that the concentration profile gets an uplifting effect with the buoyancy ratio parameter. When the buoyancy ratio parameter R b is varied by higher values, the concentration distribution ϕ ζ shows a rising inclination in the graph, which depicts that the fluid is concentrating more and more in the boundary layer region. The motility profile is the characteristic of moving microorganisms mixed in the fluid. Figure 1d is the description of this whole effect. It can be noticed that when the buoyancy ratio parameter is enhanced, then the motility of the microorganism also increases; i.e., the motile particles become more active. The bioconvection Rayleigh number N r is a dimensionless parameter that is used to measure the fluid instability due to the density differences as described in study [39]. The effect of N r on velocity distribution is considered to be very interesting, as it is the effect of instability that is caused by the dense layer of microorganisms at the top of the fluid. This top dense layer breaks, and ultimately, the microorganisms fell down. Then, they move up to attain their motive (gravitative, phototactic, chemotactic, etc.). Figure 2a is the exhibition of the result between the Rayleigh number and velocity distribution. It displays that when the bioconvection Rayleigh number is uplifted, it causes a reduction in the velocity distribution. The investigation of N r on the temperature profile for the power law index m = 0.5 , 0   &   0.5 is illustrated by Figure 2b. The conduct of N r on thermal distribution disclosed that the temperature profile rises when we give a higher value to the Rayleigh number. The concentration profile is one of the main characteristics of the Rayleigh number, as it causes the instability of the fluid. Figure 2c depicts that N r tends to grow the concentration of nanoparticles in fluids. Actually, this is one of the main interests of our work to control the sedimentation of added/mixed nanoparticles. Figure 2d illustrates the behavior of a motile microorganism versus the bioconvection Rayleigh number. It can be noticed that the density function of the motile microorganism is increased with higher Rayleigh number values. Generally, the motile microorganism characteristics of the mobile particles in a fluid are used for an increased heat transfer phenomenon. The behavior of stretching parameter α has been examined in Figure 3a. Regarding increments in the stretching parameter α , the velocity f of the fluid particles diminishes. Here, the power law index is applied on discussed m = 0.5 , 0 , 0.5 respectively. We can say that physically, the stretching parameter α provides more strength to velocity distribution f . The impact of the stretching parameter α on the concentration function ϕ with the same fixed values of the power law index parameter m is plotted in Figure 3b. The plot depicts that when m is taken to be 0.5 , the concentration profile ϕ increases as we raise the stretching parameter α , but a contrasting effect is observed when we take m = 0 , 0.5 , so that the volumetric nanoparticle concentration distribution ϕ shows a retarding effect for increasing the value of stretching parameter α . The phenomenon of the motile microorganism of the fluid for the flow of a modified second-grade nanofluid is measured on assorted values of stretching parameter α for three individual values of the power law index m = 0.5 , 0 , 0.5 . The complete sketch is figured out in Figure 3c. When the power law index is taken to be m = 0.5 , then the motility of the liquid intensifies for rising values of stretching parameter α , but when m = 0 , 0.5 , the trend for the concentration of motile microorganisms χ goes down as we uplift the stretching parameter α . Hence, we conclude that the profile of χ increases for shear thickening and decreases for modified second-grade and shear thinning fluids, as described in Figure 3d. The effect of the first-order velocity slip parameter with the power law index on velocity is sketched via Figure 4a. The first-order velocity slip parameter is also a medium characteristic affecting the fluid of the velocity distribution. The sketched graph shows that for all the three chosen cases of power law index m = 0.5 , 0 , 0.5 , when we uplift the first-order velocity slip, it causes a reduction in the velocity distribution. Figure 4b illustrates the behavior of the second-order velocity slip β parameter on velocity distribution for an addition with power law index m = 0.5 , 0 , 0.5 . Actually, this slip parameter is a medium characteristic that affects fluid flow. The velocity distribution retards when we vary β with higher values. Generally, β gives a slope of second order for the fluid flow with different characteristics, which reduces the velocity distribution. Hence, enhancement in the value of β causes decay in the velocity contour. The mixed convection parameter Λ 1 is the phenomenon of the combined interaction of the pressure forces and buoyant forces altogether. For power law index m , the effect of Λ 1 becomes more influential for velocity distribution. The impact of Λ 1 on the velocity profile in the presence of three different priorities for the power law index, i.e., ( m = 0.5 , 0 , 0.5 ) , is depicted via Figure 4c. It can be visualized that when the mixed convection is intensified, it boosts up the velocity profile of the fluid in the presence of all the values of m . Figure 4d elucidates the behavior of Lewis number L e on the non-dimensional temperature function θ . It is inspected that the curve of the temperature distribution rises with the growing Lewis number L e values for three cases of the power law index m = 0.5 , 0 , 0.5 . The impact of the Lewis number L e on the concentration profile ϕ ζ is exhibited in Figure 5b. It is found that the concentration profile ϕ shows retardation with the rise of Lewis number L e for various power law index values: m = 0.5 , 0 , 0.5 . The variation of dimensionless temperature distribution versus Prandtl number is delineated in Figure 6a. It is observed that the θ ζ shows a diminishing pattern for larger Prandtl numbers, Pr . The thermal boundry layer thickness is contracted as the fluid charactaricstics change from shear thinning to shear thickninng. This is explained by the fact that the thermal diffusivity decreases, and therefore, the thermal distribution declines. The effect of temperature ratio parameter θ w on the temperature profile in the presence of the power law index is drawn graphically through Figure 6b. It is depicted that the strengthening of the temperature ratio parameter tends to improve the temperature of the flowing fluid. When it is gradually rises, the fluid also become hot. The radiation parameter is the phenomenon of the emission of heat waves. The effect of radiation parameter R d on the temperature in the presence of three different values of power law index ( m = 0.5 , 0 , 0.5 ) is sketched by Figure 6c. The sketched information reveals that the power law index varied from 0.5 to 0.5 and the larger value of the radiation parameter causes an enhancement of the temperature function. Generally, it can be said that when the radiation parameter increases, it directly increases the temperature of the fluid. The effect of Biot number on the temperature field can be predicted in the presence of the power law index. The predicted result between the Biot number and temperature profile is sketched by Figure 7a. The figure represents that when the Biot number is raised, it causes the temperature profile of the liquid to rise as well. Generally, we can say that the increment in the value of the Biot number enhances the temperature of the fluid. The Biot number is also marked as a key influence on the variation of the concentration profile. Figure 7b illustrates the effect of the Biot number on the concentration profile in the presence of the power law index. The value of the power law index is varied from m = 0.5 to 0.5 , and the effect is noticed when the Biot number is uplifted; thus, it increases the concentration profile. Figure 8a visualizes the perception of the Prandtl number on the concentration of nanoparticles. As we vary the Prandtl number, the concentration field reduces from shear thinning to shear thickening. This is caused by the fact that the rise in Prandtl number stands for a reduction in thermal diffusivity, which creates a decline in the concentration region. Figure 8b illustrates the impact of activation energy on concentration distribution in the presence of the power law index. The sketched information indicates that the activation energy has a unique effect on the concentration field. It is seen that when the amount of activation energy is increased, it directly enhances the concentration of nanoparticles in the fluid flow. Figure 9a depicts the influence of Peclet number P e against the density of the motile microorganism profile. It is observed that an increment in the value of the Peclet number P e causes decay in motility distribution. In addition, the motile microorganism profile declines as we enhance the value of the bioconvection Lewis number L b as interpreted through Figure 9b. In Table 1, the present numerical results are verified by comparing the exact and numerical solutions as available in the existing literature for special cases by Masood et al. [33]. The numerical values of the local skin friction and local Nusselt number are tabulated in the absence of nanoparticles. Both solutions are noticed to be in strong agreement, and this verifies our numerical technique. Table 2 describes the effects of the local skin friction coefficient f 0 for the shear thinning and shear thickening fluid against the different prominent parameters such as stretching parameter α , mixed convection parameter Λ 1 , buoyancy ratio parameter R b , bioconvection Rayleigh number N r , and Wu’s slip (second-order velocity slip) parameters Γ and β , respectively. The increasing behavior of f 0 is observed for m = 0.5 as compared to m = 0.5 , 0 when it decreases. As the value of the slip parameters enhances, the value for the skin friction coefficient declines. Table 3, Table 4 and Table 5 show the local Nusselt number θ 0 , local Sherwood number ϕ 0 , and the density of the local motile microorganism χ 0 respectively against the dimensionless parameters such as R d , N t , N b , P e , L b , Γ , β . The uplifted behavior of θ 0 , ϕ 0 , and χ 0 is seen for m = 0.5 as compared to m = 0.5 , 0 .

5. Conclusions

In this article, we have developed the numerical investigation of the non-linear radiation and the activation energy effects on the bioconvection of generalized second-grade nanofluid flow across a stretching surface with convective condition, zero nanoparticles mass flux condition, and Wu’s slip (second-order slip) on the boundary. One of the most prominent features of generalized second-grade fluid is that it depicts the effects of shear thinning and shear thickening as well as stress effects. For m > 0 , the fluid is shear thickening, while for m < 0 , the fluid is shear thinning, and when m = 0 , it becomes a second-grade fluid. The tables and graphs depict the sensitivity of heat transfer and flow characteristics of nanofluids. Some major remarks are listed below:
The first-order and second-order slip parameters reduced the velocity profile for both shear thinning and shear thickening, but the opposite behavior is observed for larger values of the mixed convection parameter.
The temperature distribution exhibits an improvement when the values of the temperature ratio parameter and the Lewis number are hosted.
The increment in stretching parameter discloses interesting effects. It is noticed that for shear thickening m > 0 , the associated thermal boundary layer become thicker, and for the shear thinning cases m < 0 , this boundary layer becomes thinner. However, the momentum boundary layer is reduced for both the cases.
The nanoparticles concentration profile declines as the Lewis number and Brownian motion parameter values increase.
As the buoyancy ratio parameter, the bioconvection Rayleigh number, and thermophoresis parameter values intensify, the temperature function is boosted up as well.
Mounting Prandtl number values reduce the thermal distribution, but the radiation parameter enhances the temperature distribution.
The volumetric nanoparticles concentration profile retards as the Prandtl number rises, but it can be more effectively boosted up in the presence of activation energy.
The gyrotactic microorganism concentration profile retards by enlarging the Peclet number and bioconvection Lewis number.

Author Contributions

Y.L.; H.W.; M.I.; U.F.; F.M. and I.T. modeled the problem, numerically computed results, discussed the results physically, computed the tabulated results, wrote the manuscript and proof read it. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

This project was funded by the Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah, Saudi Arabia under grant no. (KEP-16-130-40). Therefore, the authors acknowledge with thanks DSR technical and financial support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Choi, S.U.; Eastman, J.A. Enhancing thermal conductivity of fluids with nanoparticles. Argonne Natl. Lab. Il USA 1995, 1–8. [Google Scholar]
  2. Buongiorno, J. Convective transport in nanofluids. J. Heat Transf. 2006, 128, 240–250. [Google Scholar] [CrossRef]
  3. Khan, M.I.; Hayat, T.; Waqas, M.; Alsaedi, A.; Khan, M.I. Effectiveness of radiative heat flux in MHD flow of Jeffrey-nanofluid subject to Brownian and thermophoresis diffusions. J. Hydrodyn. 2019, 31, 421–427. [Google Scholar] [CrossRef]
  4. Li, Z.; Sheikholeslami, M.; Mittal, A.S.; Shafee, A.; Haq, R.U. Nanofluid heat transfer in a porous duct in the presence of Lorentz forces using the lattice Boltzmann method. Eur. Phys. J. Plus 2019, 134, 30. [Google Scholar] [CrossRef]
  5. Aziz, A.; Muhammad, T.; Alsaedi, A.; Hayat, T. An optimal study for 3D rotating flow of Oldroyd-B nanofluid with convectively heated surface. J. Braz. Soc. Mech. Sci. Eng. 2019, 41, 236. [Google Scholar] [CrossRef]
  6. Li, Z.; Sheikholeslami, M.; Shafee, A.; Ramzan, M.; Kandasamy, R.; Al-Mdallal, Q.M. Influence of adding nanoparticles on solidification in a heat storage system considering radiation effect. J. Mol. Liquids. 2019, 273, 589–605. [Google Scholar] [CrossRef]
  7. Safaei, M.R.; Karimipour, A.; Abdollahi, A.; Nguyen, T.K. The investigation of thermal radiation and free convection heat transfer mechanisms of nanofluid inside a shallow cavity by lattice Boltzmann method. Phys. A Stat. Mech. Its Appl. 2018, 509, 515–535. [Google Scholar] [CrossRef]
  8. Sivasankaran, S.; Alsabery, A.I.; Hashim, I. Internal heat generation effect on transient natural convection in a nanofluid-saturated local thermal non-equilibrium porous inclined cavity. Phys. A Stat. Mech. Its Appl. 2018, 509, 275–293. [Google Scholar] [CrossRef]
  9. Sheikholeslami, M.; Keramati, H.; Shafee, A.; Li, Z.; Alawad, O.A.; Tlili, I. Nanofluid MHD forced convection heat transfer around the elliptic obstacle inside a permeable lid drive 3D enclosure considering lattice Boltzmann method. Phys. A: Stat. Mech. Its Appl. 2019, 523, 87–104. [Google Scholar] [CrossRef]
  10. Farshad, S.A.; Sheikholeslami, M. Simulation of nanoparticles second law treatment inside a solar collector considering turbulent flow. Physics A 2019, 525, 1–12. [Google Scholar] [CrossRef]
  11. Khan, S.U.; Waqas, H.; Shehzad, S.A.; Imran, M. Theoretical analysis of tangent hyperbolic nanoparticles with combined electrical MHD, activation energy and Wu’s slip features: A mathematical model. Phys. Scr. 2019, 94, 125211. [Google Scholar] [CrossRef]
  12. Li, F.; Sheikholeslami, M.; Dara, R.N.; Jafaryar, M.; Shafee, A.; Nguyen-Thoi, T.; Li, Z. Numerical study for nanofluid behavior inside a storage finned enclosure involving melting process. J. Mol. Liquids. 2020, 297, 111939. [Google Scholar] [CrossRef]
  13. Alamri, S.Z.; Ellahi, R.; Shehzad, N.; Zeeshan, A. Convective radiative plane Poiseuille flow of nanofluid through porous medium with slip: An application of Stefan blowing. J. Mol. Liq. 2019, 273, 292–304. [Google Scholar] [CrossRef]
  14. Kumar, K.G. Exploration of flow and heat transfer of non-Newtonian nanofluid over a stretching sheet by considering slip factor. Int. J. Numer. Methods Heat Fluid Flow 2019. [Google Scholar] [CrossRef]
  15. Farhangmehr, V.; Moghadasi, H.; Asiaei, S. A nanofluid MHD flow with heat and mass transfers over a sheet by nonlinear boundary conditions: Heat and mass transfer’s enhancement. J. Cent. South Univ. 2019, 26, 1205–1217. [Google Scholar] [CrossRef]
  16. Tlili, I.; Hamadneh, N.N.; Khan, W.A. Thermodynamic Analysis of MHD Heat and Mass Transfer of Nanofluids Past a Static Wedge with Navier Slip and Convective Boundary Conditions. Arab. J. Sci. Eng. 2019, 44, 1255–1267. [Google Scholar] [CrossRef]
  17. Khan, M. Effects of multiple slip on flow of magneto-Carreau fluid along wedge with chemically reactive species. Neural Comput. Appl. 2018, 30, 2191–2203. [Google Scholar] [CrossRef]
  18. Khan, M.; Azam, M.; Alshomrani, A.S. Unsteady slip flow of Carreau nanofluid over a wedge with nonlinear radiation and new mass flux condition. Results Phys. 2017, 2261–2270. [Google Scholar] [CrossRef]
  19. Khan, S.U.; Tlili, I.; Waqas, H.; Imran, M. Effects of nonlinear thermal radiation and activation energy on modified second-grade nanofluid with Cattaneo–Christov expressions. J. Therm. Anal. Calorim. 2020, 139, 1–12. [Google Scholar] [CrossRef]
  20. Rahman, M.; Manzur, M.; Khan, M. Mixed convection heat transfer to modified second grade fluid in the presence of thermal radiation. J. Mol. Liq. 2016, 223, 217–223. [Google Scholar] [CrossRef]
  21. Khan, M.; Rahman, M.; Manzur, M. Axisymmetric flow and heat transfer to modified second grade fluid over a radially stretching sheet. Results Phys. 2017, 7, 878–889. [Google Scholar] [CrossRef]
  22. Waqas, H.; Imran, M.; Khan, S.U.; Shehzad, S.A.; Meraj, M.A. Slip flow of Maxwell viscoelasticity-based micropolar nanoparticles with porous medium: a numerical study. Appl. Math. Mech. 2019, 40, 1255–1268. [Google Scholar] [CrossRef]
  23. Man, C.S.; Sun, Q.X. On the significance of normal stress effects in the flow of glaciers. J. Glacial. 1987, 33, 268–273. [Google Scholar] [CrossRef] [Green Version]
  24. Kuznetsov, A.V. Thermo-bioconvection in a suspension of oxytactic bacteria. Int. Commun. Heat Mass Transf. 2005, 32, 991–999. [Google Scholar] [CrossRef]
  25. Khan, S.U.; Waqas, H.; Bhatti, M.M.; Imran, M. Bio-convection in the Rheology of Magnetized Couple Stress Nanofluid Featuring Activation Energy and Wu’s Slip. J. Non Equilib. Thermodyn. 2020, 28, 81–95. [Google Scholar] [CrossRef]
  26. Waqas, H.; Khan, S.U.; Imran, M.; Bhatti, M.M. thermally developed Falkner–Skan bio-convection flow of a magnetized nanofluid in the presence of a motile gyrotactic microorganism: Buongiorno’s nanofluid model. Phys. Scr. 2019, 94, 115304. [Google Scholar] [CrossRef]
  27. Waqas, H.; Shehzad, S.A.; Khan, S.U.; Imran, M. Novel Numerical Computations on Flow of Nanoparticles in Porous Rotating Disk with Multiple Slip Effects and Microorganisms. J. Nanofluids 2019, 8, 1423–1432. [Google Scholar] [CrossRef]
  28. Waqas, H.; Khan, S.U.; Hassan, M.; Bhatti, M.M.; Imran, M. Analysis on the bio-convection flow of modified second-grade nanofluid containing gyrotactic microorganisms and nanoparticles. J. Mol. Liq. 2019, 291, 111231. [Google Scholar] [CrossRef]
  29. Waqas, H.; Khan, S.U.; Shehzad, S.A.; Imran, M. Radiative flow of Maxwell nanofluid containing gyrotactic microorganism and energy activation with convective Nield conditions. Heat Transf. Asian Res. 2019, 48, 1663–1687. [Google Scholar] [CrossRef] [Green Version]
  30. Waqas, H.; Khan, S.U.; Tlili, I.; Muhammad, A.; Mustafa, S.S. Significance of Bioconvective and thermally dissipation flow of viscoelastic nanoparticles with activation energy features: Novel biofuel significance. Symmetry 2020, 12, 214. [Google Scholar] [CrossRef] [Green Version]
  31. Khan, M.; Manzur, M.; Rahman, M. On axisymmetric flow and heat transfer of Cross fluid over a radially stretching sheet. Results Phys. 2017, 7, 3767–3772. [Google Scholar] [CrossRef]
  32. Aksoy, Y.; Pakdemirli, M.; Khalique, C.M. Boundary layer equations and stretching sheet solutions for the modified second grade fluid. Int. J. Eng. Sci. 2007, 45, 829–841. [Google Scholar] [CrossRef]
  33. Khan, M.; Rahman, M.U. Flow and heat transfer to modified second grade fluid over a non-linear stretching sheet. Aip Adv. 2015, 19, 087157. [Google Scholar] [CrossRef]
  34. Rahman, M.; Khan, M.; Manzur, M. Boundary layer flow and heat transfer of a modified second grade nanofluid with new mass flux condition. Results Phys. 2018, 10, 594–600. [Google Scholar] [CrossRef]
  35. Fang, T.; Yao, S.; Zhang, J.; Aziz, A. Viscous flow over a shrinking sheet with a second order slip flow model. Commun. Nonlinear Sci. Numer. Simul. 2010, 15, 1831–1842. [Google Scholar] [CrossRef]
  36. Nandeppanavar, M.M.; Vajravelu, K.; Abel, M.S.; Siddalingappa, M.N. Second order slip flow and heat transfer over a stretching sheet with non-linear Navier boundary condition. Int. J. Therm. Sci. 2012, 58, 143–150. [Google Scholar] [CrossRef]
  37. Khan, M.; Irfan, M.; Khan, W.A. Numerical assessment of solar energy aspects on 3D magneto-Carreau nanofluid: A revised proposed relation. Int. J. Hydrogen Energy 2017, 22054–22065. [Google Scholar] [CrossRef]
  38. Hayat, T.; Ali, S.; Farooq, M.A.; Alsaedi, A. On comparison of series and numerical solutions for flow of Eyring-Powell fluid with Newtonian heating and internal heat generation/absorption. PLoS ONE 2015, 10. [Google Scholar] [CrossRef] [Green Version]
  39. Avramenko, A.A.; Kuznetsov, A.V. Bio-thermal convection caused by combined effects of swimming of oxytactic bacteria and inclined temperature gradient in a shallow fluid layer. Int. J. Numer. Methods Heat Fluid Flow 2010, 20, 157–173. [Google Scholar] [CrossRef]
Figure 1. (ad): Variations of f , θ , ϕ , χ for various values of R b when N r = 0.2 ,   α = 0.2 ,   Γ = 1.0 ,   β = 1.0 , Λ 1 = 0.2 ,   L e = 2.0 ,   Pr = 1.0 ,   θ w = 1.5 ,   R d = 0.4 ,   B i = 2.0 ,   P r = 0.7 ,   E = 0.2 ,   P e = 0.1 ,   L b = 2.0 .
Figure 1. (ad): Variations of f , θ , ϕ , χ for various values of R b when N r = 0.2 ,   α = 0.2 ,   Γ = 1.0 ,   β = 1.0 , Λ 1 = 0.2 ,   L e = 2.0 ,   Pr = 1.0 ,   θ w = 1.5 ,   R d = 0.4 ,   B i = 2.0 ,   P r = 0.7 ,   E = 0.2 ,   P e = 0.1 ,   L b = 2.0 .
Symmetry 12 00393 g001
Figure 2. (ad): Variations of f , θ , ϕ , χ for various values of N r when R b = 0.2 ,   α = 0.2 ,   Γ = 1.0 ,   β = 1.0 , Λ 1 = 0.2 ,   L e = 2.0 ,   Pr = 1.0 ,   θ w = 1.5 ,   R d = 0.4 ,   B i = 2.0 ,   P r = 0.7 ,   E = 0.2 ,   P e = 0.1 ,   L b = 2.0 .
Figure 2. (ad): Variations of f , θ , ϕ , χ for various values of N r when R b = 0.2 ,   α = 0.2 ,   Γ = 1.0 ,   β = 1.0 , Λ 1 = 0.2 ,   L e = 2.0 ,   Pr = 1.0 ,   θ w = 1.5 ,   R d = 0.4 ,   B i = 2.0 ,   P r = 0.7 ,   E = 0.2 ,   P e = 0.1 ,   L b = 2.0 .
Symmetry 12 00393 g002
Figure 3. (ad): Variations of f , θ , ϕ , χ for various values of α when R b = 0.2 ,   N r = 0.2 ,   Γ = 1.0 ,   β = 1.0 , Λ 1 = 0.2 ,   L e = 2.0 ,   Pr = 1.0 ,   θ w = 1.5 ,   R d = 0.4 ,   B i = 2.0 ,   P r = 0.7 ,   E = 0.2 ,   P e = 0.1 ,   L b = 2.0 .
Figure 3. (ad): Variations of f , θ , ϕ , χ for various values of α when R b = 0.2 ,   N r = 0.2 ,   Γ = 1.0 ,   β = 1.0 , Λ 1 = 0.2 ,   L e = 2.0 ,   Pr = 1.0 ,   θ w = 1.5 ,   R d = 0.4 ,   B i = 2.0 ,   P r = 0.7 ,   E = 0.2 ,   P e = 0.1 ,   L b = 2.0 .
Symmetry 12 00393 g003aSymmetry 12 00393 g003b
Figure 4. (ac): Variations of f for various parameters Γ , β , Λ 1 when R b = 0.2 ,   N r = 0.2 ,   α = 0.2 ,   L e = 2.0 , Pr = 1.0 ,   θ w = 1.5 ,   R d = 0.4 ,   B i = 2.0 ,   P r = 0.7 ,   E = 0.2 ,   P e = 0.1 ,   L b = 2.0 .
Figure 4. (ac): Variations of f for various parameters Γ , β , Λ 1 when R b = 0.2 ,   N r = 0.2 ,   α = 0.2 ,   L e = 2.0 , Pr = 1.0 ,   θ w = 1.5 ,   R d = 0.4 ,   B i = 2.0 ,   P r = 0.7 ,   E = 0.2 ,   P e = 0.1 ,   L b = 2.0 .
Symmetry 12 00393 g004
Figure 5. (a,b): Variations of θ , ϕ for various values of L e when R b = 0.2 ,   N r = 0.2 ,   α = 0.2 ,   Γ = 1.0 ,   β = 1.0 , Λ 1 = 0.2 ,   Pr = 1.0 ,   θ w = 1.5 ,   R d = 0.4 ,   B i = 2.0 ,   P r = 0.7 ,   E = 0.2 ,   P e = 0.1 ,   L b = 2.0 .
Figure 5. (a,b): Variations of θ , ϕ for various values of L e when R b = 0.2 ,   N r = 0.2 ,   α = 0.2 ,   Γ = 1.0 ,   β = 1.0 , Λ 1 = 0.2 ,   Pr = 1.0 ,   θ w = 1.5 ,   R d = 0.4 ,   B i = 2.0 ,   P r = 0.7 ,   E = 0.2 ,   P e = 0.1 ,   L b = 2.0 .
Symmetry 12 00393 g005
Figure 6. (ac): Variations of θ , for various parametrs Pr , θ w , R d , when R b = 0.2 ,   N r = 0.2 ,   α = 0.2 ,   Γ = 1.0 , β = 1.0 ,   Λ 1 = 0.2 ,   L e = 2.0 ,   B i = 2.0 ,   P r = 0.7 ,   E = 0.2 ,   P e = 0.1 ,   L b = 2.0 .
Figure 6. (ac): Variations of θ , for various parametrs Pr , θ w , R d , when R b = 0.2 ,   N r = 0.2 ,   α = 0.2 ,   Γ = 1.0 , β = 1.0 ,   Λ 1 = 0.2 ,   L e = 2.0 ,   B i = 2.0 ,   P r = 0.7 ,   E = 0.2 ,   P e = 0.1 ,   L b = 2.0 .
Symmetry 12 00393 g006aSymmetry 12 00393 g006b
Figure 7. (a,b): Variations of θ , ϕ , for various values of B i when R b = 0.2 ,   N r = 0.2 ,   α = 0.2 ,   Γ = 1.0 , β = 1.0 ,   Λ 1 = 0.2 ,   L e = 2.0 ,   Pr = 1.0 ,   θ w = 1.5 ,   R d = 0.4 ,   P r = 0.7 ,   E = 0.2 ,   P e = 0.1 ,   L b = 2.0 .
Figure 7. (a,b): Variations of θ , ϕ , for various values of B i when R b = 0.2 ,   N r = 0.2 ,   α = 0.2 ,   Γ = 1.0 , β = 1.0 ,   Λ 1 = 0.2 ,   L e = 2.0 ,   Pr = 1.0 ,   θ w = 1.5 ,   R d = 0.4 ,   P r = 0.7 ,   E = 0.2 ,   P e = 0.1 ,   L b = 2.0 .
Symmetry 12 00393 g007
Figure 8. (a,b): Variations of ϕ for various values of parameters Pr , E when R b = 0.2 ,   N r = 0.2 ,   α = 0.2 ,   Γ = 1.0 , β = 1.0 ,   Λ 1 = 0.2 ,   L e = 2.0 ,   Pr = 1.0 ,   θ w = 1.5 ,   R d = 0.4 ,   B i = 2 , 0 ,   P e = 0.1 ,   L b = 2.0 .
Figure 8. (a,b): Variations of ϕ for various values of parameters Pr , E when R b = 0.2 ,   N r = 0.2 ,   α = 0.2 ,   Γ = 1.0 , β = 1.0 ,   Λ 1 = 0.2 ,   L e = 2.0 ,   Pr = 1.0 ,   θ w = 1.5 ,   R d = 0.4 ,   B i = 2 , 0 ,   P e = 0.1 ,   L b = 2.0 .
Symmetry 12 00393 g008
Figure 9. (a,b): Variations of χ for various values of parameters P e , L b when R b = 0.2 ,   N r = 0.2 ,   α = 0.2 ,   Γ = 1.0 , β = 1.0 ,   Λ 1 = 0.2 ,   L e = 2.0 ,   Pr = 1.0 ,   θ w = 1.5 ,   R d = 0.4 ,   B i = 2 , 0 ,   P r = 0.7 ,   E = 0.2 .
Figure 9. (a,b): Variations of χ for various values of parameters P e , L b when R b = 0.2 ,   N r = 0.2 ,   α = 0.2 ,   Γ = 1.0 , β = 1.0 ,   Λ 1 = 0.2 ,   L e = 2.0 ,   Pr = 1.0 ,   θ w = 1.5 ,   R d = 0.4 ,   B i = 2 , 0 ,   P r = 0.7 ,   E = 0.2 .
Symmetry 12 00393 g009
Table 1. Comparison for the results of f 0 and θ 0 in the case of second-grade fluid m = 0 , α = 1 , λ = N r = N c = , θ w = R d = E = P e = L b = 0 .
Table 1. Comparison for the results of f 0 and θ 0 in the case of second-grade fluid m = 0 , α = 1 , λ = N r = N c = , θ w = R d = E = P e = L b = 0 .
α * Pr Masood et al. [39]Present Results
Exact SolutionNumerical SolutionExact SolutionNumerical Solution
f 0 θ 0 f 0 θ 0 f 0 θ 0 f 0 θ 0
0.5100.816496582.34787450.8164511602.34787040.8164511612.34787010.816496452.3478744
1.0 0.707106782.37156830.707161772.37155440.707161702.37155420.707106962.3715684
1.5 0.632455532.38770340.632576702.38767360.632576722.38767300.63245572.3877036
2.0 0.577350272.3995950.577557262.39954500.57755722.39954520.57735012.3995950
2.0 0.95141934 0.9514135 0.9514137 0.9514193
5.0 1.6081636 1.6081591 1.6081599 1.6081636
7.0 1.9354025 1.9353982 1.9353985 1.9354025
Table 2. Variations in f 0 against α , Λ 1 , R b , N r , Γ , and β .
Table 2. Variations in f 0 against α , Λ 1 , R b , N r , Γ , and β .
α Λ 1 R b N r Γ β f 0
m = 0.5 m = 0 m = 0.5
0.10.20.20.21.0−1.00.34150.35590.3691
0.40.32770.34240.3546
0.80.31580.33000.3407
0.50.20.20.21.0−1.00.33880.35470.3684
0.40.33540.35420.3680
0.60.33390.35590.3674
0.50.20.10.21.0−1.00.34420.35710.3699
0.50.34400.35690.3697
1.00.34380.35670.3695
0.50.20.20.11.0−1.00.34170.35600.3692
0.50.34150.35580.3690
1.00.34130.35560.3688
0.50.20.20.22.0−1.00.24870.25770.2643
3.00.19650.20260.2067
4.00.16240.16700.1699
0.50.20.20.21.0−2.00.27320.28290.2941
−3.00.23370.23920.2474
−4.00.20710.20980.2155
Table 3. Variations in θ 0 against Pr , ε , θ w , N t , N b , L e , Λ 1 , R b , and N r .
Table 3. Variations in θ 0 against Pr , ε , θ w , N t , N b , L e , Λ 1 , R b , and N r .
Pr ε θ w N t N b L e Λ 1 R d R b N r θ 0
m = 0.5 m = 0 m = 0.5
110.50.30.220.10.50.20.20.27870.29260.3024
20.36830.39190.4077
30.43360.46360.4824
100.10.50.30.220.10.50.20.20.38330.35630.3489
0.40.36870.34240.3355
0.80.35140.32600.3193
1010.10.30.220.10.50.20.20.35330.32760.3209
0.40.33320.30830.3019
0.80.31330.28920.2832
1010.50.10.220.10.50.20.20.55770.59860.6238
0.40.48950.51960.5394
0.50.42130.43810.4506
1010.50.30.120.10.50.20.20.52280.55840.5810
0.50.52230.55800.5807
1.00.52210.55780.5805
1010.50.30.210.20.50.20.20.55320.59150.6151
1.50.54040.57730.6005
1.80.53240.56870.5916
1010.50.30.220.20.50.20.20.52360.55890.5812
0.30.52640.56070.5818
0.40.52770.56130.5816
1010.50.30.220.20.10.20.20.56680.60460.6286
0.40.51000.54500.5672
0.50.46800.50030.5211
1010.50.30.220.20.50.10.20.51660.55260.5761
0.50.52430.55990.5823
1.00.53070.56620.5879
1010.50.30.220.20.50.20.10.52210.55790.5806
0.50.52260.55820.5809
1.00.52310.55860.5811
Table 4. Variations in ϕ 0 against Pr , E , σ , N t , N b , L e , Λ 1 , R d , R b , and N r .
Table 4. Variations in ϕ 0 against Pr , E , σ , N t , N b , L e , Λ 1 , R d , R b , and N r .
Pr E σ N t N b L e Λ 1 R d R b N r ϕ 0
m = 0.5 m = 0 m = 0.5
110.50.30.220.10.50.20.20.41810.43880.4536
20.55250.58780.6115
30.65040.69480.7236
100.10.50.30.220.10.50.20.20.56090.52070.5100
0.50.56330.53200.5110
1.00.56500.52300.5120
1010.10.30.220.10.50.20.20.55280.51580.5062
0.50.53710.50450.4968
1.00.53050.49920.4920
1010.50.10.220.10.50.20.20.27890.29930.3119
0.41.22381.29901.3485
0.52.10652.19072.2531
1010.50.30.120.10.50.20.21.56851.67521.7429
0.50.31340.33480.3484
1.00.15670.16740.1742
1010.50.30.21.00.20.50.20.20.89970.88720.9226
1.40.81050.86600.9007
1.70.79860.85310.8875
1010.50.30.220.20.50.20.20.78540.83830.8718
0.30.78960.84100.8724
0.40.79160.84200.8727
1010.50.30.220.20.10.20.20.85010.90690.9428
0.40.76500.81740.8509
0.50.70200.75050.7817
1010.50.30.220.20.50.10.20.77490.82890.8642
0.50.78640.83980.8735
1.00.79600.84930.8818
1010.50.30.220.20.50.20.10.78320.83680.8710
0.50.78390.83740.8713
1.00.78460.83790.8717
Table 5. Variations in χ 0 against P e , Γ , β , Λ 1 , L b , R b , and N r .
Table 5. Variations in χ 0 against P e , Γ , β , Λ 1 , L b , R b , and N r .
P e Γ β Λ 1 L b R b N r χ 0
m = 0.5 m = 0 m = 0.5
0.310.50.120.20.20.55020.64170.6981
0.50.62700.72350.7831
0.70.70390.86350.9632
0.51.20.50.120.20.20.43210.50000.5367
1.60.40600.46310.4915
2.00.38790.43790.4613
0.510.50.220.20.20.43890.51050.5534
0.41600.47780.5131
0.38980.43630.4620
0.510.50.220.20.20.47810.56350.6155
0.30.49030.57290.6215
0.40.49740.57840.6247
0.510.50.21.00.20.20.38460.44030.4741
1.40.42970.50130.5452
1.80.47370.56020.6133
0.510.50.220.10.20.45680.54470.6003
0.50.47870.56490.6173
1.00.49660.58190.6322
0.510.50.220.20.10.47260.55940.6127
0.50.47400.56040.6135
1.00.47530.56150.6142

Share and Cite

MDPI and ACS Style

Li, Y.; Waqas, H.; Imran, M.; Farooq, U.; Mallawi, F.; Tlili, I. A Numerical Exploration of Modified Second-Grade Nanofluid with Motile Microorganisms, Thermal Radiation, and Wu’s Slip. Symmetry 2020, 12, 393. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12030393

AMA Style

Li Y, Waqas H, Imran M, Farooq U, Mallawi F, Tlili I. A Numerical Exploration of Modified Second-Grade Nanofluid with Motile Microorganisms, Thermal Radiation, and Wu’s Slip. Symmetry. 2020; 12(3):393. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12030393

Chicago/Turabian Style

Li, Yurong, Hassan Waqas, Muhammad Imran, Umar Farooq, Fouad Mallawi, and Iskander Tlili. 2020. "A Numerical Exploration of Modified Second-Grade Nanofluid with Motile Microorganisms, Thermal Radiation, and Wu’s Slip" Symmetry 12, no. 3: 393. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12030393

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop