# Exploration of Temperature-Dependent Thermal Conductivity and Diffusion Coefficient for Thermal and Mass Transportation in Sutterby Nanofluid Model over a Stretching Cylinder

Authors: Rabeeah Raza; Muhammad Sohail; Thabet Abdeljawad; Rahila Naz; Phatiphat Thounthong
Journal: Complexity
 (2021)

---

## Abstract

This declaration ponders the impacts of Joule warm, separation, and warming radiation for the progression of MHD Sutterby nanofluid past over an all-inclusive chamber. The wonder of warmth and mass conduction is demonstrated under warm conductivity relying upon temperature and dispersion coefficients individually. Besides, the conventional Fourier and Fick laws have been applied in the outflows of warm and mass transport. The control model comprising of a progression of coupled incomplete differential conditions is changed over into a standard arrangement of nonlinear coupled differential conditions by reasonable likeness changes. The subsequent arrangement of articulations is systematically treated through an ideal homotopic method. The impacts of various dimensionless stream boundaries on the speed, temperature, and focus fields are delineated through diagrams. The range of some parameters involved is assumed for the convergent solution as0<Re<10, 0<Pr<6.5, 0<Ec<40, 0<Rd<1.5, 0<S1<0.5, 0<S2<0.5, 0<Le<0.5, 0<Nt<2.5, and 0<Nb<2.0. The patterns of skin friction coefficient, local Nusselt, and Sherwood numbers are examined via bar charts. The principle consequence of the proposed study is that the decay of the speed for the Sutterby liquid boundary, the deterioration of the variable warm conductivity, the temperature, and the radiation increase the framework temperature. The delineation boundaries show the opposite conduct for the temperature and fixation outskirts layers.


---

## Body

## 1. Introduction

In the momentum period, non-Newtonian liquids are progressively significant for analysts as a result of their assorted applications in the mechanical field. For example, the assembling of polymers, gas turbines, power generators, glass fabrics, paper creation, wire drawing, and the sky is the limit from there. Sutterby liquid is one of the significant non-Newtonian fluids under the perceptions of the analysts. Among these, the Sutterby model addresses the high polymer fluid arrangements. It is one of the liquid models which expand the huge conduct of non-Newtonian liquids like pseudoplastic and dilatant liquids and it has the attributes to catch both the shear diminishing and shear thickening nature of the stream. Bilal et al. [1] researched Sutterby’s MHD liquid stream in the Darcy medium. They performed OHAM to appraise the outcomes and found that speed demonstrated the contrary conduct for Darcy’s opposition. Stale Sutterby liquid stream disregarding a contracting sheet highlighting the Cattaneo–Christov hypothesis going with heat retention/age was dissected by Rehman et al. [2]. Nawaz [3] contemplated the warm execution of the mixture Sutterby nanofluid stream and applied the FEM strategy to process results. He noticed that surface warmth transition is higher on account of crossbreed nanofluid than that of the base liquid. Rana and Nawaz [4] delineated the warmth improvement model in Sutterby nanofluid stream. They likewise considered Koo–Kleinstreuer and Li (KKL) Cattaneo–Christov heat motion models and found that temperature is altogether affected by the Joule warming impact. Sutterby nanofluid electromagnetohydrodynamic limit layer Darcy Forchheimer stream getting over an extending surface with a slip limit was clarified by Bilal et al. [5]. They applied the shooting technique modernized with Cash and Crap coefficients. They remarked that, in correlation with a traditional case, the drag on a superficial level for the Sutterby case shows a piddling pattern.Warmth move is a significant wonder as it is an approach to move energy to obtain the necessary errand through various cycles. Warmth transport is needed in chilling of electrical thing, inner purifier with radiator, heat exchanger, chilling of reactor, and so forth. Transport of warmth and mass exchange assume a vital part in designing and modern fields. Trials show that steady warm conductivity and mass dispersion exist for extremely uncommon cases and results got with consistent properties do not approve for all circumstances. That is the reason there is a need to present a variable warm model. Variable thermophysical properties in the liquid stream are essential to expanding the productivity of the framework. Hayat et al. [6] examined blended, convected Casson liquid stream with temperature-subordinate warm conductivity. They contemplated heat move within the sight of gooey dispersal and find that warm limit upgrades for developing estimations of little boundary use to fluctuate warm conductivity. Temperature-subordinate warm conductivity in the stagnation guide stream of Jeffrey liquid to nonstraight extended surface having variable divider thickness is examined by Hayat et al. [7]. They expressed that temperature shows a decrease in the rising warm unwinding boundary. Khan et al. [8] endeavor to delineate the stagnation direct progression of the 2nd grade to an impermeable extending chamber. They consider temperature-subordinate warm conductivity and warm definition and infer that Fourier’s law model has a higher temperature than the Cattaneo–Christov heat motion model. Hammad et al. [9] examined radiation impact for the gooey liquid stream with temperature-subordinate thickness and warm conductivity disregarding a penetrable plate and pronounce that convective warmth move rots the temperature and speed of the stream.Nanofluid is standing out enough to be noticed by scientists around the globe to expand the productivity of energy. Nanoparticles can assimilate straightforwardly from the wellspring of the episode radiation and improve the fluids radiative properties. Likewise, warm radiations alongside nanoparticles are primarily centered around the examination. They help in using sun-oriented energy as radiation heat moves. The sun is the less expensive source of the energy. Ijaz et al. [10] investigated the entropy and system of energy actuation and twofold compound response for the progression of nonlinear radiative Sisko liquid ignoring a radially stretchable rotatory circle. The huge ascribes of this examination are thick dispersal, Joule warming, and blended convection. Stagnation point blended convective MHD stream of cross nanofluid toward the surface which is extending straightly highlighting the idea of initiation energy shown by Khan et al. [11]. Parts of thermophoresis, Brownian movement, and nonlinear radiation are additionally introduced. They advise that for higher assessment of actuation energy variable types of focus rises. Stagnation point upper convective Maxwell nanofluid streaming over directly extending sheet alongside the credits of sun-oriented radiation and gooey dispersal was broken down by Khan et al. [12]. The shooting strategy is utilized to register results and they found that temperature is emphatically affected for radiative warmth move. Waqas et al. [13] noticed the effects of the attractive dipole and nonlinear warm radiation on liquid whose rheology was described via Carreau nanofluid. They applied the shooting technique to get an answer and expressed that the bigger ferrohydrodynamic connection variable tumbles down the speed.Stratification has a critical part in the advancement of businesses. It applies in the dismissal framework, condensers, and nuclear power. Khan et al. [14] analyzed the impacts of variable consistency and warm separation on Williamson nanofluid over a nondirect extending sheet and reasoned that speed and temperature rot for bigger warm delineation boundaries. Axisymmetric Jeffrey liquid stagnation point stream over a slanted porous stretchable chamber was concentrated by Ijaz and Ayub [15]. The ramifications which are thought about are the Cattaneo–Christov hypothesis, warm separation, and homogeneous-heterogeneous. They uncovered that temperature is diminishing the capacity of warm definition. Rehman et al. [16] outlined the blended convection digression exaggerated stream by extending round and hollow surface which is drenched in twofold defining media with no-slip conditions. They guaranteed that warm and solutal delineation is the reason for a decrease in temperature and convergence of the stream. Nanofluid extended stream with motile microorganisms with the ascribes of a slanted attractive field, definition, and gooey dissemination being inspected by Alsaedi et al. [17]. For additional investigations and examinations, perusers allude to the investigations announced in [18, 19].Existing writing demonstrates that there is no examination so far to see the highlights of temperature-subordinate warm conductivity and mass dispersion coefficient for the progression of magnetohydrodynamic Sutterby nanofluid past over a stretchable chamber where warm radiation as sunlight based energy is assuming a significant part in the stream model with respect to the fact that its application is concerned. This correspondence fills this gap.


## 2. Mathematical Modelling under Boundary Layer Theory (MMUBLT)

The transportation of warmth and mass for the steady and incompressible two-dimensional Sutterby nanofluid stream on a stretch chamber involving variable thermophysical qualities is considered in the current examination in Figure1. It additionally considers the ramifications of straight warm radiation, layering warm and focus, Joule warming, thermophoresis, Brownian movement, and slanted attractive power. The complete description of the modelled equations is explained in Appendix. Sutterby’s liquid model [1, 4, 20, 21] is characterized as(1)S=μ2sinh−1bγ˙bγ˙mA1.Figure 1
Geometry of the model under consideration.(2)∂u∂r+ur+∂w∂z=0,(3)u∂w∂r+w∂w∂z=νf2r∂w∂r+νf2∂2w∂r2−νfmb24∂w∂r2∂2w∂r2−δfB02ρfsin2αw,(4)ρcpfu∂T∂r+w∂T∂z=1r∂∂rrKT∂T∂r+16σ∗T∞33k∗∂2T∂r2+δfB02sin2αw2+ρcppDB∂T∂r∂C∂r+DtT∞∂T∂r2,(5)u∂C∂r+w∂C∂z=1r∂∂rrDT∂C∂r+DtT∞∂2T∂r2+1r∂T∂r.Here,γ˙=1/2Π, and Π=TrA12. Under the suspicion of limit layer, the demonstrated incomplete differential condition framework relating to the protection laws of mass, movement, warmth, and species transport is given asThe marvel considered with various actual perspectives considered is in Figure1. Related limit conditions are (6)u=0,w=Uw=U0zl,T=Tw=T0+a1zl,C=Cw=C0+b1zlatr=R,w⟶0,T⟶T∞=T0+a2zl,C⟶C∞=C0+b2zlasr⟶∞.The numerical models for warm conductivity and mass dissemination coefficient are referenced as [22, 23](7)KT=k∞1+εT−T∞Tw−T0,DT=D∞1+ε1T−T∞Tw−T0.The similarity transformations mentioned below [1, 20] are used to convert equations (3)–(6) into coupled ordinary differential equations:(8)η=U0νflr2−R22R,u=−νfU0lRrfη,w=U0zlf′η,θη=T−T∞Tw−T0,φη=C−C∞Cw−C0.Exploiting the above-declared renovations, equations (3)–(6) are as the following formation:(9)1+2γηf‴+2γ+1f″−2f′2−2βγRe1+2γηf″3−2βRe1+2γη2f″2f‴−Msin2αf′=0,(10)1+2γηθ″+Pr∞fθ′+2γθ′+Nb∞1+2γηθ′φ′+Nt1+2γηθ′2+Rd∞1+2γηθ″+γθ′+Pr∞EcMsin2αf′2−Pr∞S1+θf′=0,,(11)Nb∞1+2γηφ″+2Nb∞γφ′+Nb∞Lefφ′−S2+φf′+Nt1+2γηθ″+2γθ′=0.Transformed conditions for equations (9)–(11) are(12)f=0,f′=1,θ=1−S1,φ=1−S2atη=0,f′⟶0,θ⟶0,φ⟶0asη⟶∞,where(13)Pr∞=Pr1+εθ,Nb∞=Nb1+ε1θ,Rd∞=Rd1+εθ.Incorporating equation (13) in equations (10)–(11), we have(14)1+2γηθ″+Pr1+εθfθ′+2γθ′+Nb1+ε1θ1+2γηθ′φ′+Nt1+2γηθ′2+Rd1+εθ1+2γηθ″+γθ′+Pr1+εθEcMsin2αf′2−Pr1+εθS1+θf′=0,Nb1+2γηφ″+2Nbγφ′+NbLefφ′−S2+φf′+Nt1+ε1θ1+2γηθ″+2γθ′=0.The dimensionless parameters which are shown in the above equations are defined as follows:(15)γ=1RνflU0,β=mb2U024l2,M=δfB02lρfU0,Re=U0Z2νfl,Pr=νfρcpfk∞,Nt=ρcppD∞Tw−T0νfρcpfT∞,Nb=ρcppDBCw−C0νfρcpf,Rd=16σ∗T∞33K∗k∞,Le=νfDB,Ec=Uw2cpTw−T0,S1=a1a2,S2=b1b2.Coefficient of Skin erosion [1] and neighborhood warmth and mass transitions boundaries in nondimensional structures are deciphered as (16)cf=Re1/2f″0−β4f″04f′0+4γ2f20−4f0f′0+λf′′20,Nuz=−Re1/21+43Rd+εθθ′0,Shz=−Re1/21+εθφ′0,where λ=z2U0/νfl is a dimensionless parameter.


## 3. Solution via OHAM

This plan is broadly used to deal with nonlinear coupled issues that emerge in numerical material science. This technique is utilized by numerous specialists in their investigations. The itemized plans of this technique are treated in [1, 20–27]. The underlying assessments and their comparing direct administrators for speed, temperature, and dimensionless fixation fields are(17)f0η=1−e−η,θ0η=1−S1e−η,φ0η=1−S1e−η,Ł1=D3−Df,Ł2=D2−1θ,Ł3=D2−1φ.


## 4. Graphical Analysis and Physical Interpretation and Justification

### 4.1. Influence of Involved Parameters on Fluid Velocity

The emerging nonlinear arrangement of conditions has been handled by means of OHAM. The union of the applied calculation by processing the blunders is introduced in Table1. In the current segment effect of various boundaries on stream, profiles have been examined as charts and their representation is incorporated. Figures 2–5 are plotted to examine the conduct of speed concerning diverse dimensionless boundaries. In Figures 2 and 3 one can see the combined impact of Reynold number Re, Sutterby fluid parameter β , and curvature parameter γ on both dimensionless components of velocities fη and f′η, respectively. The fluid becomes denser for higher values of Re which as a result slows downs the fluid motion due to the fact that velocity shows a decline curve. Also, growing value of fluid parameter shows the same reversing behaviour on fluid motion because physically increasing β means enhancement in relaxation time which dominates the viscous effect and produces resistance to fluid flow. So, the impact of curvature parameter is concerned with higher curvature showing the fast fluid flow in both velocity profiles. Figures 4 and 5 are presented to show the impact of magnetic parameter M along with its angle of inclination α on both velocity profiles fη and f′η, respectively. The growth in magnetic parameter results in arising the opposing Lorentz force to the fluid which downshifts the fluid motion. The same is the case with α as increasing the angle enhances the hold of the magnetic force on the flow. On the whole, velocity of the flow shows downturns when any parameter becomes strong and creates resistance to the flow.Table 1
Convergence examination of established series solutions.
mE^kfE^kθE^kφ20.01126140.00469210.0052450540.003102990.0003172580.0017800660.00114390.00008523810.0010251380.0003604550.00004478330.000739084100.00008628070.00002790750.00059564120.00002061960.00001955690.000510457146.83992×10−60.00001488870.000453768162.71667×10−70.00001168630.000412904189.56869×10−79.26413×10−60.000381734202.94672×10−77.39892×10−60.00035696Figure 2
Influence ofRe,β,andγonfηFigure 3
Influence ofRe,β,andγonf′η.Figure 4
Influence ofMandαonfη.Figure 5
Influence ofMandαonf′η.


### 4.2. Comportment of Contributing Parameters on Temperature Field

Figures6–14 are designed for the impact of several flow parameters on dimensionless temperature θη. Figure 6 shows the view of curvature impact on temperature θη. As velocity is an upsurging function of curvature and temperature is also defined through average kinetic energy. Higher curvature implies higher velocity which implies higher kinetic energy results in the rising temperature profile. The declining outturns of a small parameter ε on temperature are shown in Figure 7. For ε=0, the case reduces to normal temperature-independent thermal conductivity case and shows the maximum temperature on it which implies that variable thermal conductivity cools down the flow with higher temperature. Figure 8 shows the impact of the small parameter ε1 (used to vary the diffusion coefficient) on temperature θη. For higher values of ε1, the graph shows a rise in temperature. Figure 9 is portraying the declining behaviour of temperature curves for growing values of Prandtl number Pr since the thickness of the thermal boundary layer reduced by enhancing Pr. As Prandtl number is inversely related to thermal diffusivity, enhancing it results in cooling the flow. The combined impact of the Brownian motion parameter Nb and thermophoresis parameter Nt is depicted in Figure 10. Speeding the Brownian motion leads to the faster random motion of nanoparticles in a flow which shows an extension in thermal boundary layer thickness and enhances the temperature of flow more rapidly. A similar pattern is observed for growing values Nt. As in process of thermophoresis, more heated particles near the surface travel away from heated regions toward the cold region and raise temperature there and collective temperature of the whole system rises. Figure 11 illustrates the impression of radiation parameter Rd on temperature θη. Radiation parameter tells about the relative contribution of heat transfer through conduction to thermal radiation. So, a clear upsurge in temperature curves is observed for escalating values of Rd. Enhancing thermal radiation means adding energy through radiation in the flow which eventually boosts the temperature distribution of the fluid. Figure 12 shows the impression of the Eckert number Ec on dimensionless temperature θη. The growth in Ec has a positive influence on heating the system due to the fact that it strengthens the kinetic energy of the flow. The graphical view for the behaviour of temperature when magnetic parameter M and its inclination angle α are taken into account is illustrated in Figure 13. The resistive force which arises due to the presence of a magnetic field is the cause of boosting the temperature distribution of flow. Figure 14 is a pictorial illustration of the stratification parameter S1 due to temperature difference. The thermal stratification process causes layer formation due to temperature variation, the enlarging stratification parameter falls down the temperature gradient between the wall and reference point, and thus the temperature of the system drops.Figure 6
Influence ofγonθη.Figure 7
Influence ofεonθη.Figure 8
Influence ofε1onθη.Figure 9
Influence ofPronθη.Figure 10
Influence ofNbandNtonθη.Figure 11
Influence ofRdonθη.Figure 12
Influence ofEconθη.Figure 13
Influence ofMandαonθη.Figure 14
Influence ofS1onθη.


### 4.3. Bearing of Involved Parameters on Fluid Concentration Field

Figures15–18 are mapped to describe the etiquette of dimensionless concentration ϕη by the variety of different flow parameters. Figure 15 is an image of thermophoresis parameter Nt influence on concentration. Boosting thermophoresis process boosts the traveling of particles from higher to lower temperature difference which as a result maximizes the concentration of nanoparticles in the flow regime. Figure 16 is plotted to understand the impact of Lewis's number Le on concentration ϕη. As Le gives the relation between the thermal boundary layer and the concentration boundary layer, physically, escalation in Le shrinks the thickness of the boundary layer of concentration. Figure 17 is planned to study the behaviour of concentration ϕη on the basis of a small meter ε1 (used to vary diffusion coefficient). An increasing behaviour is observed for elevating values of ε1, Figure 18 is sketched to illustrate the impact of the stratification parameter S2 due to concentration difference. The higher values of S2 are a reason for the decline in concentration difference between wall and reference point which as a result lessens the concentration boundary layer thickness.Figure 15
Influence ofNtonφη.Figure 16
Influence ofLeonφη.Figure 17
Influence ofε1onφη.Figure 18
Influence ofS2onφη.


### 4.4. Impacts of Involved Parameters oncfRe−1/2,NuzRe−1/2,andShzRe−1/2

Figures19–21 are bar chart representation of coefficient of skin friction, local Nusselt number, and local Sherwood number based on several parameters. Figure 19 shows the illustration of the effect of the Sutterby fluid parameter β on skin friction coefficient cfRe−1/2. Enhancing β has a reverse reaction on flow velocity as it enhances the viscous nature of the fluid due to which opposing force arises which is the cause of boosting the skin friction coefficient. Figure 20 gives a pictorial view for the impressions of radiation parameter Rd on the heat flux in the form of local Nusselt number NuzRe−1/2. Higher radiation parameter lowers the local Nusselt number as it also has a reverse effect on the momentum boundary layer. Physical temperature gradient at the wall surface becomes small which lowers natural convection ability. Figure 21 shows the footprints of stratification concentration parameter S2 on local Sherwood number ShzRe−1/2and the bar chart shows a reverse behaviour for the escalation in S2. In Table 2, a limiting case of dimensionless stress is computed and a comparison is established. It is found that, due to the factor 2 in the momentum equation, our results are differing by the multiple of 2. If we ignore 2, our results exactly match with the precious finding.Figure 19
Influence ofβ on skin friction coefficient cfRe−1/2.Figure 20
Influence ofRdon local Nusselt numberNuzRe−1/2.Figure 21
Influence ofS2on local Sherwood number ShzRe−1/2.Table 2
Comparative study for skin friction coefficient.
γf″0 [28]f″0 present0–1.000−2.4159799501377160.25–1.094378−2.4727793625851120.5–1.188715−2.5408757184912120.75–1.281833−2.59795493209095161–1.459308−2.599693623769820


## 4.1. Influence of Involved Parameters on Fluid Velocity

The emerging nonlinear arrangement of conditions has been handled by means of OHAM. The union of the applied calculation by processing the blunders is introduced in Table1. In the current segment effect of various boundaries on stream, profiles have been examined as charts and their representation is incorporated. Figures 2–5 are plotted to examine the conduct of speed concerning diverse dimensionless boundaries. In Figures 2 and 3 one can see the combined impact of Reynold number Re, Sutterby fluid parameter β , and curvature parameter γ on both dimensionless components of velocities fη and f′η, respectively. The fluid becomes denser for higher values of Re which as a result slows downs the fluid motion due to the fact that velocity shows a decline curve. Also, growing value of fluid parameter shows the same reversing behaviour on fluid motion because physically increasing β means enhancement in relaxation time which dominates the viscous effect and produces resistance to fluid flow. So, the impact of curvature parameter is concerned with higher curvature showing the fast fluid flow in both velocity profiles. Figures 4 and 5 are presented to show the impact of magnetic parameter M along with its angle of inclination α on both velocity profiles fη and f′η, respectively. The growth in magnetic parameter results in arising the opposing Lorentz force to the fluid which downshifts the fluid motion. The same is the case with α as increasing the angle enhances the hold of the magnetic force on the flow. On the whole, velocity of the flow shows downturns when any parameter becomes strong and creates resistance to the flow.Table 1
Convergence examination of established series solutions.
mE^kfE^kθE^kφ20.01126140.00469210.0052450540.003102990.0003172580.0017800660.00114390.00008523810.0010251380.0003604550.00004478330.000739084100.00008628070.00002790750.00059564120.00002061960.00001955690.000510457146.83992×10−60.00001488870.000453768162.71667×10−70.00001168630.000412904189.56869×10−79.26413×10−60.000381734202.94672×10−77.39892×10−60.00035696Figure 2
Influence ofRe,β,andγonfηFigure 3
Influence ofRe,β,andγonf′η.Figure 4
Influence ofMandαonfη.Figure 5
Influence ofMandαonf′η.


## 4.2. Comportment of Contributing Parameters on Temperature Field

Figures6–14 are designed for the impact of several flow parameters on dimensionless temperature θη. Figure 6 shows the view of curvature impact on temperature θη. As velocity is an upsurging function of curvature and temperature is also defined through average kinetic energy. Higher curvature implies higher velocity which implies higher kinetic energy results in the rising temperature profile. The declining outturns of a small parameter ε on temperature are shown in Figure 7. For ε=0, the case reduces to normal temperature-independent thermal conductivity case and shows the maximum temperature on it which implies that variable thermal conductivity cools down the flow with higher temperature. Figure 8 shows the impact of the small parameter ε1 (used to vary the diffusion coefficient) on temperature θη. For higher values of ε1, the graph shows a rise in temperature. Figure 9 is portraying the declining behaviour of temperature curves for growing values of Prandtl number Pr since the thickness of the thermal boundary layer reduced by enhancing Pr. As Prandtl number is inversely related to thermal diffusivity, enhancing it results in cooling the flow. The combined impact of the Brownian motion parameter Nb and thermophoresis parameter Nt is depicted in Figure 10. Speeding the Brownian motion leads to the faster random motion of nanoparticles in a flow which shows an extension in thermal boundary layer thickness and enhances the temperature of flow more rapidly. A similar pattern is observed for growing values Nt. As in process of thermophoresis, more heated particles near the surface travel away from heated regions toward the cold region and raise temperature there and collective temperature of the whole system rises. Figure 11 illustrates the impression of radiation parameter Rd on temperature θη. Radiation parameter tells about the relative contribution of heat transfer through conduction to thermal radiation. So, a clear upsurge in temperature curves is observed for escalating values of Rd. Enhancing thermal radiation means adding energy through radiation in the flow which eventually boosts the temperature distribution of the fluid. Figure 12 shows the impression of the Eckert number Ec on dimensionless temperature θη. The growth in Ec has a positive influence on heating the system due to the fact that it strengthens the kinetic energy of the flow. The graphical view for the behaviour of temperature when magnetic parameter M and its inclination angle α are taken into account is illustrated in Figure 13. The resistive force which arises due to the presence of a magnetic field is the cause of boosting the temperature distribution of flow. Figure 14 is a pictorial illustration of the stratification parameter S1 due to temperature difference. The thermal stratification process causes layer formation due to temperature variation, the enlarging stratification parameter falls down the temperature gradient between the wall and reference point, and thus the temperature of the system drops.Figure 6
Influence ofγonθη.Figure 7
Influence ofεonθη.Figure 8
Influence ofε1onθη.Figure 9
Influence ofPronθη.Figure 10
Influence ofNbandNtonθη.Figure 11
Influence ofRdonθη.Figure 12
Influence ofEconθη.Figure 13
Influence ofMandαonθη.Figure 14
Influence ofS1onθη.


## 4.3. Bearing of Involved Parameters on Fluid Concentration Field

Figures15–18 are mapped to describe the etiquette of dimensionless concentration ϕη by the variety of different flow parameters. Figure 15 is an image of thermophoresis parameter Nt influence on concentration. Boosting thermophoresis process boosts the traveling of particles from higher to lower temperature difference which as a result maximizes the concentration of nanoparticles in the flow regime. Figure 16 is plotted to understand the impact of Lewis's number Le on concentration ϕη. As Le gives the relation between the thermal boundary layer and the concentration boundary layer, physically, escalation in Le shrinks the thickness of the boundary layer of concentration. Figure 17 is planned to study the behaviour of concentration ϕη on the basis of a small meter ε1 (used to vary diffusion coefficient). An increasing behaviour is observed for elevating values of ε1, Figure 18 is sketched to illustrate the impact of the stratification parameter S2 due to concentration difference. The higher values of S2 are a reason for the decline in concentration difference between wall and reference point which as a result lessens the concentration boundary layer thickness.Figure 15
Influence ofNtonφη.Figure 16
Influence ofLeonφη.Figure 17
Influence ofε1onφη.Figure 18
Influence ofS2onφη.


## 4.4. Impacts of Involved Parameters oncfRe−1/2,NuzRe−1/2,andShzRe−1/2

Figures19–21 are bar chart representation of coefficient of skin friction, local Nusselt number, and local Sherwood number based on several parameters. Figure 19 shows the illustration of the effect of the Sutterby fluid parameter β on skin friction coefficient cfRe−1/2. Enhancing β has a reverse reaction on flow velocity as it enhances the viscous nature of the fluid due to which opposing force arises which is the cause of boosting the skin friction coefficient. Figure 20 gives a pictorial view for the impressions of radiation parameter Rd on the heat flux in the form of local Nusselt number NuzRe−1/2. Higher radiation parameter lowers the local Nusselt number as it also has a reverse effect on the momentum boundary layer. Physical temperature gradient at the wall surface becomes small which lowers natural convection ability. Figure 21 shows the footprints of stratification concentration parameter S2 on local Sherwood number ShzRe−1/2and the bar chart shows a reverse behaviour for the escalation in S2. In Table 2, a limiting case of dimensionless stress is computed and a comparison is established. It is found that, due to the factor 2 in the momentum equation, our results are differing by the multiple of 2. If we ignore 2, our results exactly match with the precious finding.Figure 19
Influence ofβ on skin friction coefficient cfRe−1/2.Figure 20
Influence ofRdon local Nusselt numberNuzRe−1/2.Figure 21
Influence ofS2on local Sherwood number ShzRe−1/2.Table 2
Comparative study for skin friction coefficient.
γf″0 [28]f″0 present0–1.000−2.4159799501377160.25–1.094378−2.4727793625851120.5–1.188715−2.5408757184912120.75–1.281833−2.59795493209095161–1.459308−2.599693623769820


## 5. Conclusion and Key Points of Current Research

In the course of the study, a two-dimensional stationary radiation flux of Sutterby nanofluid with an oblique magnetic field and with variable thermophysical properties was taken under observation. To solve the system of transformed ODEs, the OHAM scheme was applied. The main findings of this study are as follows:Larger curvatureγ boosts both velocity and temperature profiles.Small parameterε reduces the temperature and ε1 shows opposite behaviour.Radiation parameterRd has a positive impact on the temperature profile.Velocity field slows down for higher Reynolds numberRe, magnetic parameter M, and Sutterby fluid parameter β.Larger values of stratification parametersS1 and S2 reduce the thermal and concentration boundary layer thickness in flow.Skin friction increases with an increase in Sutterby fluid parameter and radiation decays the local Nusselt number.Augmentation in Eckert number and radiation parameter enhance the thermal transport.



---

*Source: 6252864-2021-03-02.xml*