Mathematical Problems in Engineering
Volume 2021 (2021), Article ID 5520780, 22 pages
https://doi.org/10.1155/2021/5520780

Research Article

A Numerical Study of MHD Carreau Nanofluid Flow with Gyrotactic Microorganisms over a Plate, Wedge, and Stagnation Point

R. M. Akram Muntazir, M. Mushtaq, and K. Jabeen

Department of Mathematics, University of the Engineering and Technology, Lahore 54890, Pakistan

Correspondence should be addressed to K. Jabeen; kanwaljabeen44@yahoo.com

Received 10 February 2021; Revised 10 March 2021; Accepted 26 March 2021; Published 13 April 2021

Academic Editor: Ali Ahmadian

Copyright © 2021 R. M. Akram Muntazir et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This article addresses the numerical exploration of steady and 2D flow of MHD Carreau nanofluid filled with motile microorganisms over three different geometries, i.e., plate, wedge, and stagnation point of a flat plate. The influence of magnetic field, viscous dissipation, thermophoresis, and Brownian motion is considered for both cases, i.e., shear thinning and shear thickening. A set of relevant similarity transformations are utilized to obtain dimensionless form of governing coupled nonlinear partial differential equations (PDEs). The transformed system of ordinary differential equations (ODEs) is then numerically solved by bvp4c via MATLAB based on shooting technique and Runge–Kutta–Fehlberg (RKF) scheme via MAPLE. Also, a numerical analysis has been made for skin friction factor, heat, and mass transfer rates. Results elucidate that all the profiles except velocity show decreasing behavior for higher values of magnetic field parameter. Among all three flow geometries for both shear thinning and shear thickening cases, the flow over a plate has lesser skin friction factor. The nanoparticle concentration and density of motile microorganism decrease in both the shear thinning and shear thickening cases, for increasing values of Brownian motion (Nb), but reverse trend is observed for rising values of thermophoresis parameter (Nt). Furthermore, it is observed that, as we increase the values of suction/injection parameter (S), the velocity of fluid increases but decreases the fluid temperature, concentration of mass and density of motile organisms over a plate, wedge, and stagnation point of a flat plate. Also, we observed that shear thinning nanofluid has higher rate of heat, mass, and motile microorganisms mass transfers than shear thickening fluid. Both shear thinning and thickening nanofluid have a low rate of heat/mass and gyrotactic microorganisms mass transfer over plate among wedge and stagnation point flow.

1. Introduction

The heat transfer phenomenon has gained extensive attraction for researchers due to its vital range of applications and implementations in industrial and manufacturing processes particularly in thermal insulations, nuclear waste storage, crude oil extraction, geothermal systems, enhanced oil recovery, heat exchangers, ground water pollution, etc. In prior work, boundary layer (BL) theory was related the flows driven by stream-wise pressure gradient. In fluid dynamics, a classical problem of fluid flow past a stationary plate is the Blasius problem in which the free stream is parallel to the plate and its velocity is constant. Falkner and Skan [1] had proved that their problem confessed similarity solution same like Blasius problem. Initially, they employed the Prandtl theory to transform the governing system and considered 2D flows over a wedge. They introduced a suitable similarity method in which the PDEs were converted to a nonlinear -order ODEs and then solved by numeric techniques.

In 1937, Hartee [2] discussed thoroughly the same problem and presented some useful results. Later, Lin et al. [3] had explored the heat transfer characteristics in force convectional flow for variable Prandtl number along a static wedge. Na [4] by using a set of transformations converted boundary value problem (BVP) to an initial value problems and applied forward integration scheme. Rajagopal et al. [5] explored Falkner–Skan flow for 2nd-grade fluid over a wedge positioned in symmetric way with respect to the flow direction. In addition, Kuo [6] and many more researchers presented numerical and analytical precise solutions of BVP and discussed heat/mass transfer behaviors of the flow along a wedge. Ishak et al. [7] explored 2D steady MHD flow past a wedge and plate along a variable wall temperature.

The fluids are considered as non-Newtonian fluids, when “Newton law of viscosity” is disobeyed. In non-Newtonian fluids, the shear stress and shear rate are unsystematic (highly nonlinear); in contrast, the Newtonian fluids are linear. The custard, blood, toothpaste, shampoo, honey, and paint all are non-Newtonian fluids. The classification of non-Newtonian fluids comprises viscoelastic fluid, and it covers lubricants, whipped cream, etc. The rheopectic and thixotropic (including yogurt, gelatin gels, pectin gels, peanut butter) non-Newtonian fluids can also be categorized into shear thickening (dilatant) like oobleck (suspension of cornstarch and water) and shear thinning (pseudoplastic) such as nail polish, sand in water, paper pulp in water, and micropolar fluids, but from research point of view, non-Newtonian micropolar fluids yet encapsulated a small attention.

The fundamental equations for the non-newtonian fluid flow are much complicated and highly nonlinear and imply the high-order terms. The momentum equations are incompetent to explicate the flow ability of non-Newtonian fluids. This fact is on account of complex rheological effects. This complexity can be overwhelmed as different rheological models are presented. Few of the models are Casson fluid, Maxwell fluid, Walter’s B fluid, Carreau fluid, cross fluid, the power law fluid, etc. The constitutive relation of power-law is uncomplicated but it has the deficiency that it cannot accurately predict viscosity for extremely high or low rates of shear. To overcome the limitation, another model was considered, i.e., Carreau rheological model (a combination of Newtonian and power law properties) proposed by Carreau [8]. Due to multiple characteristics, Carreau fluid has its own importance in literature. Carreau fluid is better to depict several kinds of physical issues because this fluid model has the ability to reveal the rheology of multiple specific fluids such as fluids with brief-chain suspension particles, detergents, fluid crystals, and blood in humans and animals. This flow model has gained notable attention due to its significance in extrusion of a polymer, tumors treatment, bitumen for road construction, cosmetics, etc. In the Carreaue model, the essential conditions retain at high and low shear rates unlike power law fluid. Accordingly, this fluid has received a widespread adoption in research these days. For free surface flow, the fundamental relation of Carreau model is suitable enough. The reason behind the fact is that, in Carreau model, when the shear rate approaches zero, viscosity remains finite. This fluid possesses shear thinning shear thickening and Newtonian fluid properties. At high shear rate, the non-Newtonian fluids exhibit thickening while at low shear rate and exhibit the thinning properties which makes sense that its viscosity is highly dependent on shear rate. Different researchers in different times discussed Carreau fluid. In [9], Raju explained MHD Carreau fluid with cross diffusion over a wedge and concluded that velocity field grows for increasing wedge angle parameter. Sandeep [10] had a detailed work for the consequences of nonlinear thermal radiation and an uncertain heat source, sink in the presence of homogenous and heterogeneous phenomenon. Khan et al. [11] used Buongiorno model to predict the time-dependent Carreau nanofluid along a wedge and found that Weissenberg number enhances the temperature of fluid. Waqas et al. [12] have explored two different characteristics of MHD Carreau nanofluid. Khan et al. [13] have studied and analyzed numerically the effects of heat generation/absorption and melting phenomenon of Carreau nanofluid over a wedge and noted that nanoparticle concentration field decreases by enhancing melting parameter. Khan et al. [14] have studied the MHD Carreau nanofluid transport properties for solar energy application over a 3D extending sheet. Tanmay [15] represented his work for the stability criteria of the non-Newtonian fluid taking into account the whole analysis of arbitrary Reynolds number. In [16–18], authors addressed 2D and 3D flow of MHD Carreau nanofluid with chemical reaction over an exponential surface and presented detailed discussion upon physical outcomes of various parameters such as Weissenberg number, thermal and concentration Biot numbers, Brownian movement thermophoresis, Schmidt and Prandtl numbers, thermal radiation and temperature ratio parameter for primary and secondary velocities, concentration and temperature, and for both n = 1.5 and n = 0.5 flows.

Bhatti [19] explored the heat transfer on the electromagnetohydrodynamic (EMHD) Carreau fluid flow over a rectangular plates. For physical modeling, he considered a Darcy–Brinkman–Forchheimer medium and solved the formulated flow system by numerical and semi-analytical schemes. In his work, the author observed the Forchheimer and MF suppress the flow. Also, he concluded that the We and n retard the velocity of fluid but enhance the temperature profile. He declared that DTM is more powerful technique [20, 21]. Due to vast number of applications in the industrial, technological, and engineering developments, the Carreau model was further worked out by numerous researchers/authors in the last decades [22–26].

A little while back, nanotechnology played an indispensable role of heat dissipation processes in numerous fields. In this developed era, the ultimatum of energy increasing was unwaveringly. Nanotechnology made its indigenous efforts for this ultimate demand and its most credible achievement leads us to nanoparticles. The word “nano” means “billionth part of a meter” and is denoted by a factor 0.000000001 m. The nanofluids contain nanosized particles. A very small amount, almost less than 100 nm of guest nanoparticles dispersed uniformly and suspended stably in the host liquid, comes up with a dramatic improvement in the thermal possessions of the host fluid and finally gives nanofluids which are most stable and hold high possible thermal properties with smallest volume fraction of the nanoparticles and expose rapid heat dissipation than conventional fluids. The commonly known base (conventional) fluids are oil, water, and some other lubricants, polymer solutions, biofluids, ethylene glycol, etc. This sudden rise in thermal conductive property after the emulsion of nanoparticles into the conventional fluids is owing to the fact of Brownian motion that develop the micromixing, volume fraction of particles, interfacial layer. The term was brought into being by Choi [27] at Argonne National Laboratory of USA. Nanofluids have innovative possessions that make them capable of being potentially convenient. The more relative surface area of nanoparticles, unlike the ordinary particles, not only enhances the credibility of heat dissipation but also improves the strength and stability of suspension. The desirability of nanofluids is having high specified surface area, high dispersion stability, adaptable properties which differ in the concentration of nanoparticles, and exhibiting lower pumping power as compared to the conventional fluid. The particular applications of nanofluid are in engine cooling, cooling of transformer oil, cooling of electronics, improving efficiency of diesel generator, in grinding, fuel cells, domestic refrigerator, chillers, space and defence, cooling of microchips, etc. Initially, the usage of nanofluids in heat transmit applications was too much common, but with the enhancement in worldwide development, nanofluids find their way also in biomedical engineering such as in diagnosis and therapy. Magnetic nanofluid system is a consequential track for targeted drug transportation, hyperthermia (a physiological state), and differential diagnosis. Nanofluid found its extreme way as a powerful antiseptic factor to overthrow antibacterial resistance. Therefore, this transformation of conventional fluids to the nanofluids has stimulated chemists, physicists, and engineers around the globe and is a stunning out-turn of such continuance. New areas of research in which nanofluids are developed include bioconvection, inorganic-based nanomaterials, tissue engineering, regenerative medicine, stem cells, etc. Falkner–Skan nanofluid flows along a wedge with various properties, such as heat source/sink, suction/injection, magnetic field, slip/no slip, thermophoresis, and chemical reaction were studied by many authors [28–31].

In viscous fluid flow, microorganisms are considered as the research domain to researchers, scientists, chemists, and engineers to an extensive range in the fields like agronomy, industrial point of view, or scientifically establishing a power in research area, as microorganisms play an indispensable role in plenty ways like to get rid of toxins, in digestion of food, etc. Moreover, magnetic field is essential in some of the applications of fluid flow. Furthermore, microorganism particles can be used in various commercial and industrial products, for example, biofuel, biodiesel, biofertilizers, ethanol, and bio-microsystems.

An erratic and crooked pattern of countless swimming microoraganisms is the cause for bioconvection. It is the phenomenon in which microorganism shows movement in upward and downward directions. Just because of the denseness behavior, bioconvection is basically correlated with the minuscule convection of these microorganisms. Generally, the motion of such unicellular microbes is categorized as taxis based on its incentive and stimulant motion. Moreover, these taxis are further compartmentalized into gravitaxis, gyrotaxis, and phototaxis. Such distinguished and eminent applications stimulate many of the analysts to put their benefactions in this domain. The conceptualization of nanofluid bioconvection was first presumably given by “Kuznetsov” [32, 33], and his center of attraction is to work on nanofluids containing gyrotactic microorganisms, corroborating that the locomotion of the fluid on extensive scale due to self-sustained motile microorganisms intensifies the merging and prohibits the agglomeration of nanoparticles in nanofluid. This phenomenon arises as the infinitesimal microorganisms, that are more dense with the base fluid, and the base fluid here is water float upwardly due to density gradient establishing the fact that very small microorganism of insignificant settling velocity cast no distinguishable outcome on bioconvection, also the particles with extraordinary settlement have no further possessions on bioconvection, regardless this, the particles with much appropriate and optimal size and its gravitational settlement must emulate with Brownian diffusion, so very these particles can proficiently slows down bioconvection. In short, the phenomenon of bioconvection come about to happen due to unreliable stratification, that sequentially leads to the evolution of indiscernible radiative-motion. Bioconvection involves wide range of applications in biological and bio-microsystems such as enzymes, biofuels, biosensors, oil recovery, fertilizer, and various fields of life where it shows significant and notable applications. Beside this, bioconvection flows play a key role in fuel cells, biological polymer synthesis, and environmental systems [34]. Other important advantages of bioconvection are to enhance the mass transfer process and being used to solve the mixing hurdles in microsystems.

Bhatti [35] solved the flow system by perturbation method to investigate the bioconvection flow of blood nanofluid in an annulus with gyrotactic microorganism and have seen that, for enhancing the Peclet number, microorganism density decays. Raju [36] studied a numerical study to explore the Casson nanofluid transport properties filled by gyrotactic microorganisms [37] and declared that increasing wedge angle parameter causes reducing the density of microorganisms. Ullah [38] scrutinized his model for protrusion system regarding with the bio-convention flow, using the same idea as he proposed, which may be used for intensification of extrusion system. Wang [39] did great work in medicinal field and in rheology carrying the square channel for simulation. Kotha [40] employed magnetohydrodynamics and heat and mass transfer sensations assuming the gyrotactic microorganisms in a water based nanofluid suspension in a vertical manner of a plane that comes out in form of partial differential fundamental equations, that are transformed using similarity possessions and turn into non-linear ODEs. Few works associated with bioconvection Carreau nanofluid over a wedge/plate/stagnation of the plate can be seen in [41–50].

The present survey explores the impacts of shear-thinning/shear-thickening characteristics of Carreau nanofluid filled with gyrotactic microorganisms over three different geometries, i.e., (i) flat plate, (ii) over a wedge, and (iii) at stagnation point in the presence of magnetic field, suction/injection, thermophoresis, and Brownian motion. This analysis casts a noteworthy attention as the present work involves nanofluid in view of the fact that the nanofluids are of significant worth due to their high thermal conductivity. The substructure of the whole work is to visualize the impacts on Falkner–Skan flow assuming a wedge, stagnation point, and plate of the plate with shear thinning and shear thickening cases. Appropriate similarity transformations are applied to transmute the fundamental PDEs into ODEs and used RKF and bvp4c techniques that leads to the convergence. The fluid velocity, concentration, temperature, stream function, rate of heat, and mass transfer are scrutinized with rising parameters, and their evaluations and outcomes are properly discussed using graphical approach. Evaluated graphical scheme evidently announces that the perspective shear thinning/thickening fluid exhibit an opposed behavior in both fluid velocity and temperature for inflated values of Weissenberg for the cases of shear thinning and shear thickening, respectively. The flow of Carreau nanofluid with magnetic field and suction/injection can be beneficial for improving solar energy efficiency. The scientific results investigated in this work may provide improvement in energy consumption, solar systems, and heat extrusion systems.

This paper is divided into three sections: (i) formulation of flow model, (ii) graphical and tabular representation, and (iii) results and discussion.

2. Formulation of Flow Model

The schematic configuration of three different geometries can be seen in Figure 1. The flow of 2D, viscous, incompressible, and time independent non-Newtonian Carreau nanofluid having gyrotactic microorganisms for both shear thinning and thickening cases is considered in this study. The mass/heat transfer characteristics under the impacts of viscous dissipation, suction/injection, thermal radiation, Brownian, and thermophoresis motion are studied thoroughly.

Figure 1: Schematic configuration of the flow problem.

The Carreau nanofluid particles get influenced by magnetic force of strength , Brownian motion, and thermophoresis mechanism. Total wedge angle is where is the parameter of wedge angle or measure of pressure gradient. In other words, is the Hartree pressure gradient. Here, represents the fluid flow over a plate, wedge, and stagnation point. Further suppose that is the free stream velocity of the potential flow where a and are constants. The surface temperature, density of motile microorganisms, and nanoparticle concentration at free stream and boundary are and , respectively. Also, it is assumed are fixed and greater than (Figure 1).

Based on aforementioned situation, the problem is formulated upon the following assumptions:

(i)The MHD Carreau nanofluid flow filled with gyrotactic microorganisms is steady, laminar, and incompressible(ii)Flow has no slip behavior(iii)No body force is considered in the momentum equation(iv)On surface, suction/injection can happen(v)The famous Buongiorno nanofluid model is considered to investigate the Nanofluid features, i.e., Brownian and thermophoresis movement(vi)The modeled equations for the abovementioned flow problem are numerically handled by imposing RKF and bvp4c

The fluid transport system, containing the gyrotactic microorganisms, nanoparticles concentration, momentum, and energy equations, is expressed as follows:

The constitutive problem proceeds with the following BCs:

2.1. Similarity Transformations

We introduce the similarity transformations [36, 37]:

Using PDEs ((1)–(5)) and associated BCs (7), the transformed system of ODEs is as follows:

Also, the transformed boundary conditions are

where is Weissenberg number and . Notice that is for mass injection, is for suction, is the wedge angle parameter, and is Prandtl , thermal radiation, is viscous dissipation, is Brownian movement, thermophoretic movement parameter, is Schmidt number is bioconvection Peclet number, and is difference of microorganisms concentration parameter.

The radioactive heating flux in the temperature equation is calculated and simplified by the Rosseland approximation, i.e., , where constant and  = mean absorption coefficient. The temperature variation within the flow is supposed to be small enough so can be taken as a linear function of temperature neglecting the higher order term while expanding in Taylor-series about is obtained as [20, 21]

Neglecting the higher order terms of the above said statement beyond the first degree in the series, we get

2.1.1. Physical Quantities

The significant physical parameters of engineering concerns are dimensionless skin friction coefficient , Nusselt numbers , and Sherwood number that physically indicate shear stress, rate of heat, and mass transfer, respectively, and are defined as follows:

where is the shear stress at the sheet along x-axis:

and , , and are surface heat and mass fluxes defined as

Upon fixing the local Reynold’s number and using the above defined relations, we have

3. Results and Discussion

In the presence of the above literature and the present work, the ingenuity of the conventional problem is to address the impacts of active parameters on fluid properties over a wedge where thermal jump and viscous dissipation exist and all the exertions have been made to get solution for this problem via similarity function. Furthermore, the considered physical problem gives highly nonlinear and coupled temperature dependent flow equations, while these converted ODEs subject to BCs are then computed by two different methods, bvp4c (using MATLAB) and dsolve technique (on MAPLE). Our basic assumptions are a wedge, plate, and stagnation point with shear thinning and shear thickening cases. The fixed parameters are and we elucidate their effects on fluid velocity , temperature , concentration , and gyrotactic motile , whereas the effects of the above-mentioned parameters are thoroughly investigated and are given in Tables 1–4 and their graphical outlook. To monitor the reliability of present work, we expose our results with the prior researches, Rajagopal [5], Kuo [6], and Masood [24] et al. After scrutinizing, our results are in good agreement with their results. Tables 1 and 2 are created for proving the authenticity of the present work with the already published work. It would be worth mentioning that all physical behavior of this flow problem coincides very well with already published literature either in graphical or in tabular representation.

Table 1: Calculation of of Carreau nanofluid for different values of when .

[5]
[6]
[24]
MATLAB
MAPLE

0
—
0.469600
0.469600
0.469600
0.469600

0.3
0.474755
0.775524
0.474755
0.774755
0.774755

0.6
0.995836
0.995757
0.995836
0.995836
0.995836

1.2
1.335722
1.333833
1.335722
1.335721
1.335721

Table 2: Calculation of of Carreau nanofluid for different values of when .

[6]
[37]
Present
[6]
[37]
Present

0.72
0.41809
0.41886
0.418091
0.418091
0.5298
0.5296
0.529608
0.529608

1.0
0.46960
0.46960
0.469600
0.469600
0.60541
0.6052
0.605197
0.605197

10
1.02974
1.02974
1.029747
1.029747
1.4561
1.4588
1.455750
1.455750

30
1.4873
1.4873
1.487319
1.487319
2.1582
2.1580
2.157737
2.157737

100
2.2229
2.2229
2.222906
2.222906
3.2869
3.2873
3.286250
3.286250

1000
4.7901
4.7901
4.790062
4.790062
7.2225
7.2334
7.221173
7.221173

Table 3: Calculation of , , , and for different values of .

Stagnation

0.779657
0.779656
1.098775
1.098775
1.544226
1.544226

0.742201
0.742201
1.036636
1.036636
1.448262
1.448262

0.712195
0.712198
0.991033
0.9910334
1.381483
1.381483

0.689058
.689059
0.957061
.957061
1.332688
1.332688

Shear thickening

0.863839
0.863837
1.260130
1.260130
1.822569
1.822569

0.946347
0.946347
1.399239
1.399239
2.042439
2.042439

1.018201
1.018201
1.515076
1.515076
2.220557
2.220557

1.080080
1.080080
1.612923
1.612925
2.369274
2.369275

Shear thinning

0.359687
0.359688
0.401976
0.4019761
0.473756
0.4737565

0.370326
0.370327
0.417153
0.4171527
0.495277
0.4952776

0.379118
0.379120
0.428921
0.4289213
0.511341
0.511341

0.386060
0.386061
0.437996
0.4379966
0.523556
0.5235560

Shear thickening

0.337408
0.337408
0.366525
0.366525
0.419309
0.419309

0.317775
0.317775
0.338210
0.338210
0.378466
0.378466

0.302065
0.302066
0.316093
0.316095
0.346995
0.346995

0.289385
0.289386
0.298345
0.298346
0.321817
0.321818

Shear thinning

0.592494
0.592493
0.744520
0.7445205
0.970189
0.9701885

0.596715
0.596715
0.747650
.7476495
0.971433
0.9714334

0.599205
0.599205
0.748616
0.7486163
0.970419
0.970419

0.600708
0.600709
0.748788
0.748788
0.968957
0.9689572

Shear thickeing

0.581245
0.581245
0.731875
0.731875
0.957420
0.957420

0.573951
0.573952
0.726529
0.726529
0.95531
0.955311

0.569376
0.569376
0.724228
0.724227
0.956262
0.956262

0.566316
0.566316
0.723249
0.723249
0.958158
0.958157

Shear thinning

0.703413
0.703413
0.856238
0.8562379
1.087959
1.087959

0.714318
0.714319
0.869533
0.8695330
1.104199
1.104199

0.722319
0.722317
0.878400
0.8784001
1.114421
1.114421

0.728184
0.728186
0.884688
0.8846888
1.121540
1.121540

Shear thickening

0.678455
0.678455
.820446
.820446
1.038309
1.038309

0.659986
.659987
0.797830
0.797830
1.010170
1.010170

0.646851
.646851
.782502
.782503
0.991659
0.991659

0.637043
0.637043
0.771270
.771270
.978246
.978246

Table 4: Calculation of and for different values of Nt and Pe.

Stagnation

Shear thinning

0.381564
0.381566
0.427964
0.427965
0.506142
0.506144

0.359687
0.359688
0.401976
.401977
0.473756
0.473757

0.338546
0.338547
0.376876
0.376877
0.442485
0.442486

0.318183
0.318183
0.352709
.352710
0.412382
0.412383

Shear thickening

0.358586
0.358587
0.391480
0.391482
0.450192
0.450194

0.337408
0.337408
0.366525
0.366526
0.419309
0.419311

0.316979
0.316980
0.342489
0.342490
0.38959
0.389600

0.297334
0.297336
0.319408
0.319410
0.361101
0.361102

Shear thinning

0.703413
.703413
0.856238
.856238
1.087959
1.087959

0.800922
0.800926
0.980401
.980401
1.251472
1.251472

0.900387
0.900388
1.106908
1.106908
1.417923
1.417923

1.001706
1.001706
1.235637
1.235637
1.587159
1.587160

Shear thickening

0.678455
0.678455
0.820446
0.820449
1.038309
1.038309

0.774200
.774197
0.942647
0.942646
1.199928
1.199928

0.871902
0.871898
1.067214
1.067214
1.364543
1.364543

0.971456
.971459
1.194023
1.194023
1.531992
1.531992

Table 3 is created to show the numerically calculated results of , , , and to see the impacts of Weissenberg number (We) for both shear thinning and shear thickening cases with n = 0.75 and n = 1.75, respectively, by MATLAB bvp4c and MAPLE dsolve techniques. Both methods are in good agreement, and this shows the validity of present results. From this table, we observed that the skin friction factor is the decreasing function of We in case of shear thinning fluid but in shear thickening case skin friction factor accelerates with rising values of We. But we noticed an opposite behavior in , , and . We also noticed that, among all the flow geometries (plate, stagnation point, and wedge), the plate experienced less skin friction when compared with wedge and stagnation point.

Table 4 is constructed to show the numerical results of for varying values of thermophoresis parameter (Nt) and for selected values of Peclet number (Pe) in both cases when and . The density of motile microorganism increases in both the shear thinning and shear thickening cases, for increasing values of Pe, but reverse trend is observed for rising values of Nt while calculating thermal gradient.

Figures 2–11 illustrate the fluid transport properties. The graphs are presented by dotted, dash, and solid lines orderly managed for plate, wedge, and stagnation point, manifesting the nanofluid properties that whether these properties are increasing or decreasing, or remain constant throughout the phenomenon.

(a)

(b)

(c)

(d)

Figure 2: We effects on and .

(a)

(b)

(c)

(d)

Figure 3: We effects on and .

(a)

(b)

(c)

(d)

Figure 4: Effects of M on and .

(a)

(b)

(c)

(d)

Figure 5: Effects of M on and .

(a)

(b)

(c)

(d)

Figure 6: Effects of and .

(a)

(b)

(c)

(d)

Figure 7: Effects of S and .

(a)

(b)

(c)

(d)

Figure 8: Effects of R for and on .

(a)

(b)

(c)

(d)

Figure 9: Effects of Sc and Pe on and .

(a)

(b)

(c)

(d)

(e)

(f)

Figure 10: Effects of Nb , and .

(a)

(b)

(c)

(d)

(e)

(f)

Figure 11: Effects of Nt on , , and .

Figures 2 and 3 illustrate the impacts on fluid velocity, temperature, concentration, and the gyrotactic motile microorganisms distribution signified by for both shear thinning and shear thickening cases with n = 0.75 and n = 1.75, respectively. It was noticed that for varying values of We the velocity of the fluid declines in the case when but it increases for increasing values of We when . Actually, We is the relation between relaxation time and the process in which time grows the fluid’s viscosity. Increasing the values of We enhances the relaxation time. So, there is a decrement in the fluid’s velocity due to dominating effects of viscosity which generates resistance in the fluid particles. Moreover, the plots of the variation in the thermal, solutal, and motile microorganisms distributions for various values of We can also be seen in Figures 2 and 3 for both types of fluids for three different geometries . These figures exhibit that the thermal, concentration, and motile microorganism layers thickness increases with increasing We in shear thinning case. Scientifically, ratio of elastic and viscous forces is known as We, so it elevates the thickness of fluid so temperature increases, while both profiles show decreasing trend for shear thickening fluids. It is necessary to mention here that We = 0 recovered Newtonian fluid.

The influence of magnetic parameter on velocity, temperature, concentration, and gyrotactic motile microorganisms contours is visualized graphically through Figures 4 and 5. In the BL flow of various conducting fluids, magnetic field utilization can play a vital role in controlling the momentum and heat transfer. From Figures 4 and 5, we can see the impacts of M on for distinct values of M in both shear thinning and thickening cases. It is noteworthy to mention here that the fluid velocity rises for rising values of M in both flow cases but increasing M causes decelerating the profiles for both cases. Physically, if less suppression of drag force is involved, it reduces the thermal, concentration, and gyrotactic motile microorganism profiles but elevates the velocity distribution with higher values of M.

Figures 6 and 7 depict the behavior of suction parameter S for different flow geometries on boundary layer profiles for shear thinning and shear thickening cases of MHD Carreau nanofluid. It is seen that all the profiles are decreasing except the velocity distribution that is the only increasing function of S for and when the fluid is injected through a surface, but escalating values of S depreciating the boundary layer thickness of over the plate, wedge, and stagnation point. Further, it is noticed that the stagnation point of a flat plate has a low temperature, concentration, and gyrotactic motile microorganism when compared with the plate and wedge.

Figure 8 depicts the effects of distinct values of Eckert and radiation parameter on temperature profile , where in the figure the temperature propagates. The Eckert number pops up as coefficient of internal energy changes; therefore, the upsurge of temperature is obvious. The reason behind this upsurge temperature is that when heat is dispersed because of viscous force, the heat starts flowing towards the fluid, ensuring the case of increase in Eckert number, and follows to rise the temperature of the fluid. This divulges that temperature of fluid is going to be high when change in internal energy is negligible compared with the temperature of fluid where this change is consequent. So, for both discussed cases, temperature increases. Also, Figure 8 represents the curves of temperature field for several values of radiation parameter . Due to rise in radiation parameter, this will uplift the temperature of the fluid. As more heat is transmitted to the fluid when the phenomenon of radiation occurs, in shear thinning and shear thickening of the fluid, the temperature will rise when parameter rises.

Figure 9 depicts the graphical behavior of Schmidt (Sc), and bioconvectional Peclet (Pe) numbers scrutinize that for positive values of Sc and . As Sc increases, the concentration profile decreases. The above behavior of concentration profile is according to the existing physical situations. On explanatory note, the rate of mass transfer gets smaller when diffusion coefficient increases; this retards the concentration profile for both cases. Similarly, higher molecular motion leads to sharpening process of mass transportation which results in the deceleration in the concentration profile. The motile microorganism profile also shows a declined graph. This declined scratch for and is due to the fact that motile density dwindles due to the observed relation of number with microorganism diffusivity.

Figure 10 highlights the effects of Brownian motion on temperature, concentration, and density of motile microorganism profiles of the fluid for shear thinning and shear thickening. Nanoparticles with incremented values of cast a stronger impact on heat transfer. The temperature of nanofluid escalates as the parameter increases. However, causes declining of the concentration BL because the nanoparticles diffuse away from the surface into the fluid which causes retardation in nanoparticle concentration and hence results in deceleration. Furthermore, it has been observed that, with increasing magnitude of parameter, this uplifts the rate at which nanoparticles in the conventional fluids start to move themselves in irregular directions with disparate velocities. Also, Figure 10 demonstrates the graphical illustration of density of motile microorganisms of the flow due to uplifting Brownian motion parameter . We see that a strength in parameter results in stifling the dispersion of nanoparticles quite away from the respective region that result in declining the -curves.

Figure 11 shows the effects on due to gradual elevation of thermophoresis parameter . Gradually, for different increasing values of , the thermal, concentration, and density of motile microorganism BLs tend to be enhanced in both flow cases. These nanoparticles show a powerful impact of heat transport properties of Carreau fluid. The thermophoretic force develops as a consequence of temperature gradient which in turn results in increasing the flow rate of the fluid away from plate, wedge, and stagnation point. The outcome of this rapid flow of the fluid is that as the parameter increases, the thermal boundary layer increases. The concentration increases as increases. Basically, the nanoparticles due to this thermophoresis force are activated and start moving from hot to cold surface that comes out with the consequence that increases the mass related thickness. So, rise in parameter upsurges the concentration curves of fluid flow. Figure 11 also highlights the graphical impacts of varying on profile that results in increasing the density of motile microorganism. All these figures indicate the effects for both shear thinning and thickening cases for plate, wedge, and stagnation point orderly.

4. Conclusion

In this paper, a numerical investigation on MHD Carreau nanofluid flow comprising gyrotactic microorganisms over a wedge/plate/stagnation point of the plate is performed under the influence of suction/injection, viscous dissipation, Brownian, and thermophoresis motion. The results are carried out by utilizing the bvp4c and RKF techniques. Impacts of salient flow parameters on velocity, temperature, nanoparticle concentration, and density of motile microorganism to control the fluid transport properties are discussed. Some important findings drawn from the present work are summarized as follows:

(i)The influence of Weissenberg number for various flow geometries upon Carreau nanofluid has an opposite behavior for shear thinning/shear thickening fluids on velocity, temperature, concentration, and motile microorganism distributions; i.e., as We increases, velocity profile increases in shear thinning case but decreases in shear thickening case. But, for rising values of We, the thermal, concentration, and motile microorganism profiles depreciate for shear thinning case but quite reverse behavior is seen for shear thickening fluid.(ii)A rise in suction/injection parameter S leads to enhancing the fluid velocity but decreases the temperature, nanoparticle concentration, and density of motile microorganism of the fluid. Suction or injection can be used to control boundary layer flows as it helps to bring stability in the flow.(iii)In both the shear thinning and thickening cases, the temperature increases with increasing Ec, R, Nb, and Nt.(iv)The nanoparticle concentration and density of motile microorganism decrease in both the shear thinning and shear thickening cases, for increasing values of Nb, but reverse trend is observed for rising values of Nt.(v)The shear thinning nanofluid has slightly higher rate of heat, mass, and motile microorganisms mass transfers than shear thickening fluid.(vi)All the profiles except velocity field show decreasing behavior for higher values of magnetic field parameter.(vii)Both profiles, i.e., nanoparticle concentration and density of motile microorganism profile, show decreasing trend for higher values of Sc and Pe in both flow situations and .(viii)Among all three flow geometries, i.e., plate, wedge, and stagnation point of a flat plate, for both shear thinning and shear thickening cases, the flow over a plate has lesser skin friction factor.(ix)Both shear thinning and thickening nanofluid have a low rate of heat/mass and gyrotactic microorganisms mass transfer over plate among wedge and stagnation point flow.(x)The bioconvection with nanoparticle interaction can be used in thermal transpiration for engineering and industrial processes.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

V. M. Falkner and S. W. Skan, “Some approximate solutions of the boundary-layer equations,” Philosophical Magazine, vol. 12, pp. 865–896, 1931.

D. R. Hartree, “On an equation occurring in Falkner and Skanʼs approximate treatment of the equations of the boundary layer,” Mathematical Proceedings of the Cambridge Philosophical Society, vol. 33, no. 2, pp. 223–239, 1937.

H.-T. Lin and L.-K. Lin, “Similarity solutions for laminar forced convection heat transfer from wedges to fluids of any Prandtl number,” International Journal of Heat and Mass Transfer, vol. 30, no. 6, pp. 1111–1118, 1987.

T. Y. Na, Computational Methods in Engineering Boundary Value Problems, Academic Press, New York, NY, USA, 1979.

K. R. Rajagopal, A. S. Gupta, and T. Y. Na, “A note on the Falkner-Skan flows of a non-Newtonian fluid,” International Journal of Non-Linear Mechanics, vol. 18, no. 4, pp. 313–320, 1983.

B.-L. Kuo, “Application of the differential transformation method to the solutions of Falkner-Skan wedge flow,” Acta Mechanica, vol. 164, no. 3-4, pp. 161–174, 2003.

A. Ishak, R. Nazar, and I. Pop, “Moving wedge and flat plate in amicropolar fluid,” International Journal of Engineering Science, vol. 44, no. 18-19, pp. 1225–1236, 2006.

P. J. Carreau, “Rheological equations from molecular network theories,” The Society of Rheology, vol. 116, no. 1, Article ID 99127, 1972.

C. S. K. Raju and N. Sandeep, “Falkner-Skan flow of a magnetic-Carreau fluid past a wedge in the presence of cross diffusion effects,” The European Physical Journal Plus, vol. 131, p. 267, 2016.

C. S. K. Raju and N. Sandeep, Unsteady Three Dimensional Flow of Carreau-Casson Fluid Past Stretching Surface, VIT university, Vellore, India, 2016.

M. Khan, M. Azam, and A. Munir, “On unsteady Falkner-Skan flow of MHD Carreau nanofluid past a static/moving wedge with convective surface condition,” Journal of Molecular Liquids, vol. 230, p. 4858, 2017.

M. Waqas, M. I. Khan, T. Hayat, and A. Alsaedi, “Numerical simulation for magneto Carreau nanofluid model with thermal radiation: a revised model,” Computer Methods in Applied Mechanics and Engineering, vol. 324, pp. 640–653, 2017.

M. Khan, M. Azam, and A. S. Alshomrani, “Effects of melting and heat generation/absorption on unsteady Falkner-Skan flow of Carreau nanofluid over a wedge,” International Journal of Heat and Mass Transfer, vol. 110, pp. 437–446, 2017.

M. Khan, M. Azam, and A. S. Alshomrani, “Unsteady slip flow of Carreau nanofluid over a wedge with nonlinear radiation and new mass flux condition,” Results in Physics, vol. 7, pp. 2261–2270, 2017.

V. S. Tanmay, “Stability of plane Couette flow of Carreau fluid past a deformable solid at arbitrary Reynolds number,” Physics of Fluids, vol. 30, pp. 74–103, 2018.

M. Irfan, K. Rafiq, W. A. Khan, and M. Khan, “Numerical analysis of unsteady Carreau nanofluid flow with variable conductivity,” Applied Nanoscience, vol. 10, no. 8, pp. 3075–3084, 2020.

M. Irfan, M. Khan, W. A. Khan et al., “Physical aspects of shear thinning/thickening behavior in radiative flow of magnetite Carreau nanofluid with nanoparticle mass flux conditions,” Applied Nanoscience, vol. 10, pp. 3021–3033, 2020.

M. Irfan, “tudy of Brownian motion and thermophoretic diffusion on nonlinear mixed convection flow of carreau nanofluid subject to variables properties,” Surfaces and Interfaces, vol. 23, Article ID 100926, 2021.

M. M. Bhatti, L. Phali, and C. M. Khalique, “Heat transfer effects on electro-magnetohydrodynamic Carreau fluid flow between two micro-parallel plates with DarcyBrinkman Forchheimer medium,” Archive of Applied Mechanics, 2021.

K. Jabeen, M. Mushtaq, and R. M. Akram, “Analysis of the MHD boundary layer flow over a nonlinear stretching sheet in a porous medium using semianalytical approaches,” Mathematical Problems in Engineering, vol. 2020, no. 1, Article ID 3012854, 9 pages, 2020.

K. Jabeen, M. Mush, and R. M. Akram Muntazir, “Analysis of MHD fluids around a linearly stretching sheet in porous media with thermophoresis, radiation, and chemical reaction,” Mathematical Problems in Engineering, vol. 2020, no. 1, Article ID 9685482, 14 pages, 2020.

R. K. Kumar, C. Raju, B. Mahanthesh, B. Gireesha, and S. Varma, “Chemical reaction effects on nano carreau liquid flow past a cone and a wedge with Cattaneo-Christov heat flux model,” International Journal of Chemical Reactor Engineering, vol. 16, no. 4, 2017.

K. Ganesh Kumar, B. Gireesha, N. Rudraswamy, and S. Manjunatha, “Radiative heat transfers of Carreau fluid flow over a stretching sheet with fluid particle suspension and temperature jump,” Results in Physics, vol. 7, pp. 3976–3983, 2017.

M. Khan and H. Sardar, “On steady two-dimensional Carreau fluid flow over a wedge in the presence of infinite shear rate viscosity,” Results in Physics, vol. 8, pp. 516–523, 2018.

M. R. Eid, K. L. Mahny, T. Muhammad, and M. Sheikholeslami, “Numerical treatment for Carreau nanofluid flow over a porous nonlinear stretching surface,” Results in Physics, vol. 8, pp. 1185–1193, 2018.

Hashim, “Multiple nature analysis of Carreau nanomaterial flow due to shrinking geometry with heat transfer,” Applied Nanoscience, vol. 1, pp. 3305–3314, 2020.

S. U. S. Choi, “Enhancing thermal conductivity of fluids with nanoparticle,” in Development and applications of non-Newtonian flow, FED-vol. 231/MD-vol. 66, pp. 99–105, ASME, 1995.

S. M. S. Murshed, K. C. Leong, and C. Yang, “Enhanced thermal conductivity of TiO₂-water based nanofluids,” International Journal of Thermal Sciences, vol. 44, no. 4, pp. 367–373, 2005.

R. Ahmad and W. A. Khan, “Effect of viscous dissipation and internal heat generation/absorption on heat transfer flow over A moving wedge with convective boundary condition,” Heat Transfer-Asian Research, vol. 42, no. 7, pp. 589–602, 2013.

A. A. Afify and M. A. A. Bazid, “MHD Falkner-Skan flow and heat transfer characteristics of nanofluids over a wedge with heat source/sink effects,” Journal of Computational and Theoretical Nanoscience, vol. 11, pp. 1844–1852, 2014.

R. Ellahil, “A mathematical study of non-newtonian micropolar fluid in arterial blood flow through composite stenosis,” Applied Mathematics & Information Sciences, vol. 8, pp. 1567–1573, 2014.

A. V. Kuznetsov and A. V. Avramenko, “Effect of small particles on the stability of bioconvection in a suspension of gyrotactic microorganisms in a layer of finite depth,” International Communications in Heat and Mass Transfer, vol. 31, p. 110, 2004.

A. V. Kuznetsov, “Nanofluid bioconvection in water-based suspensions containing nanoparticles and oxytactic microorganisms: oscillatory instability,” Nanoscale Research Letters, vol. 6, p. 100, 2011.

K. Das, P. R. Duari, and P. K. Kundu, “Nanofluid bioconvection in presence of gyrotactic microorganisms and chemical reaction in a porous medium,” Journal of Mechanical Science and Technology, vol. 29, pp. 4841–4849, 2015.

M. M. Bhatti, A. Zeeshan, and R. Ellahi, “Simultaneous effects of coagulation and variable magnetic field on peristaltically induced motion of Jeffrey nanofluid containing gyrotactic microorganism,” Microvascular Research, vol. 110, p. 3242, 2017.

C. S. K. Raju, S. M. Ibrahim, S. Anuradha, and P. Priyadharshini, “Bio-convection on the nonlinear radiative flow of a Carreau fluid,” The European Physical Journal Plus, vol. 131, p. 409, 2016.

C. S. K. Raju, M. Mainul Hoque, and T. Sivasankar, “Radiative flow of Casson fluid over a moving wedge filled with gyrotactic microorganisms Advanced,” Powder Technology, vol. 28, pp. 575–583, 2017.

S. U. Khan and S. A. Shehzad, “Anaysis for time dependent flow of carreau nanofluid over an accelerating surface with gyrotactic microorganism,” Advances in Mechanical Engineering, vol. 12, pp. 1350–1358, 2019.

E. Wang, “Liquid of aloe extract with shear thinning carreau fluid,” Journal of Molecular Liquids, vol. 309, Article ID 113011, 2020.

G. Kotha, V. R. Kolipaula, M. V. S. Rao, and S. Penki, “Internal heat generation on bioconvection of an MHD nanofluid flow due to gyrotactic microorganisms,” European Physical Journal Plus, vol. 135, no. 7, p. 600, 2020.

R. Naz, S. Tariq, M. Sohail, and Z. Shah, “Investigation of entropy generation in stratified MHD Carreau nanofluid with gyrotactic microorganisms under Von Neumann similarity transformations,” The European Physical Journal Plus, vol. 135, p. 178, 2020.

A. Majeed, A. Zeeshan, N. Amin, N. Ijaz, and T. Saeed, “Thermal analysis of radiative bioconvection magnetohydrodynamic flow comprising gyrotactic microorganism with activation energy,” Journal of Thermal Analysis and Calorimetry, vol. 143, pp. 2545–2556, 2021.

Z. Abdelmalek and S. U. Khan, “A proposed unsteady bioconvection model for transient thin film flow of rate-type nanoparticles configured by rotating disk,” Journal of Thermal Analysis and Calorimetry, pp. 1–16, 2020.

M. Taseer, S. Z. Alamri, H. Waqas, D. Habib, and R. Ellahi, “Bioconvection flow of magnetized Carreau nanofluid under the influence of slip over a wedge with motile microorganisms,” Journal of Thermal Analysis and Calorimetry, vol. 143, pp. 945–957, 2021.

M. Aisha, “Alqahtani, numerical investigation of heat and mass transport in the flow over a magnetized wedge by incorporating the effects of cross-diffusion gradients: applications in multiple engineering systems,” Mathematical Problems in Engineering, vol. 2020, Article ID 2475831, 10 pages, 2020.

G. Kumaran and R. Sivaraj, “Hydromagnetic forced convective flow of Carreau nano fluid over a wedge/plate/stagnation of the plate,” The European Physical Journal Special Topics, vol. 228, pp. 2647–2659, 2019.

H. T. Basha, R. Sivaraj, A. Subramanyam Reddy, and A. J. Chamkha, “SWCNH/diamond-ethylene glycol nano fluid flow over a wedge, plate and stagnation point with induced magnetic field and nonlinear radiation-solar energy application,” The European Physical Journal Special Topics, vol. 228, pp. 2531–2551, 2019.

H. T. Basha, R. Sivaraj, A. Subramanyam Reddy, A. J. Chamkha, and H. M. Baskonus, “A numerical study of the ferromagnetic flow of Carreau nanofluid over a wedge, plate and stagnation point with a magnetic dipole,” AIMS Mathematics, vol. 5, no. 5, pp. 4197–4219, 2020.

H. T. Basha and R. Sivaraj, “Numerical simulation of blood nanofluid flow over three different geometries by means of gyrotactic microorganisms: applications to the flow in a circulatory system,” Journal of Mechanical Engineering Science, pp. 1–20, 2020.

R. M. Muntazir, M. Mushtaq, S. Shahzadi, and K. Jabeen, “Influence of chemically reacting ferromagnetic carreau nanofluid over a stretched sheet with magnetic dipole and viscous dissipation,” Mathematical Problems in Engineering, vol. 2021, Article ID 6652522, 2021.