Skip to main content

ORIGINAL RESEARCH article

Front. Phys., 12 April 2023
Sec. Fluid Dynamics

Radiative simulation of non-Newtonian MHD fluid over a boundary-driven multi-physical curved mechanism: Keller–Box evidence

Kehinde M. Sanni&#x;Kehinde M. Sanni1Saleem Asghar&#x;Saleem Asghar2Isra Al-Shbeil&#x;Isra Al-Shbeil3Adriana C&#x;ta
&#x;Adriana Cătaş4*
  • 1Department of Mathematics, Chak Shahzad, COMSATS University, Islamabad, Pakistan
  • 2Department of Mathematics, COMSATS University Islamabad, Islamabad, Pakistan
  • 3Department of Mathematics, Faculty of Science, The University of Jordan, Amman, Jordan
  • 4Department of Mathematics and Computer Science, University of Oradea, Oradea, Romania

This study is numerically driven to ascertain the flow of two-dimensional heat transfer of an incompressible electrically conducting non-Newtonian fluid over a continuous power-law stretching curved surface. The flow model considers rheological fluid viscosity using curvilinear (r −, s −) coordinates. The energy equation for the curved mechanism is examined in two streams: the prescribed surface temperature and the prescribed heat flux. Surface frictional heating is influenced by thermal radiation and viscous dissipation. Similarity transformations are executed to reduce partial differential equations into ordinary differential equations. The Keller–Box shooting method with the Jacobi iterative techniques is numerically computed for the degenerated nonlinear system of the boundary value problem. The associated boundary-layer thickness and flow fields- velocity and temperature are analyzed against characterizing parameters. Significant results are obtained and discussed with graphical plots showing that fluid velocity can be controlled by virtue of fluid parameters and stretching power index. These results are useful in polymer dynamics involving the melting and manufacturing of stretchable sheets.

1 Introduction

The heat transfer analysis for non-Newtonian fluid flow driven by a nonlinear boundary surface has received great attention for decades. Most fluids found in science, engineering, and manufacturing industry, such as biological fluids, foodstuff, cosmetics and toiletries, particulate slurries, multi-phase fluid mixtures, synthetic lubricants, paints, and numerous polymeric fluids, display characteristics that do not obey classical Newton’s law of viscosity. In other words, the complex behavior of these fluids is evidenced by normal stresses, viscosity, and time-dependent elastic effect, which can be traced to being a non-Newtonian fluid. Without a unique general model to exhibit all the properties of various non-Newtonian fluids, many rheology viscous models have been proposed to capture the dynamic phenomena for the completeness of flow studies in fluid mechanics. The simplest non-Newtonian fluid model is documented as Ostwald–de Waele fluid (known as power-law fluid), which explains the occurrence of a real non-Newtonian fluid. It classifies non-Newtonian fluids as shear thinning (n < 1) and shear thickening (n > 1) with the limitation of not being able to predict the stress relaxation time and elasticity performance of viscoelastic fluids. The second-, third-, and fourth-grade fluids are subclasses of differential-type materials that can predict normal stress, shear thinning, and shear thickening but do not capture time-dependent relaxation. A linear viscoelastic (Maxwell fluid) model is characterized by relaxation time and explains the properties of polymeric fluids but fails in cognizance of the shear thinning and thickening properties of the fluids. Other non-Newtonian fluids, including Jeffrey fluid, Williamson fluid, Casson fluid, Reiner–Philipoff fluid, Ellis fluid, Eyring viscous fluid, Power–Eyring fluid, Prandtl–Erying fluid, Sisco fluid, Carreau fluid, Nanofluid, and Meter, are characterized by shearing stress, implicit and explicit expression of strain rate, and composite function of strain rate [18]. The study of non-Newtonian fluids and heat transfer has numerous applications in engineering productions and the polymer industry. It enhances the quality of the finished product through heating and cooling processes. Some recent studies in this direction include a stagnation stretching sheet for a flow of Maxwell fluid in the presence of radiative heat flux analyzed by [9]. Jawadi et al. [10] investigated the asymptotic power-law fluid flow with a numerical technique for the unsteady scenario. The Rayleigh–Stokes problem with the energy balance of a Maxwell fluid is studied by [11]. Hayat et al. [12] analyzed the flow of MHD Cross fluid with heat transfer over a stretching sheet through numerical computations. Ijaz et al. [13] considered the numerical simulation of MHD Cross fluid flow over a point-stagnation stretched sheet with heat transfer. Manzur et al. [14] examined the flow of Cross fluid influenced by thermal radiation with mixed convective boundary conditions and buoyancy effects. Khan et al. [15] presented an axisymmetric flow of Cross fluid over a radially varying stretching plate in the presence of heat transfer phenomena. For more studies, readers can visit [1622]; [23]. Recently, the non-Newtonian viscous model given by Sutterby [24] has attracted the keen interest of researchers to capture some limitations witnessed in non-Newtonian fluids. For example, the power-law (Ostwald–de Waele) model, which cannot describe zero and high viscous shear rates (i.e., η0 and η, respectively), only applies to the intermediate shear rate. The Ellis model covers the zero shear rate η0 and cannot account for the infinite shear rate η. The Power–Erying and Seely models need to be amended to fit viscosity data on the entire range of shear rates. Furthermore, the Meter and Ellis models may not have reproduced Newtonian fluid in the limit case without meeting certain conditions. However, the Sutterby fluid model returns to the Newtonian fluid when m = 0, reduces to Erying fluid if η = 0, and Erying–Power fluid as η ≠ 0 with significant flexibility of curve fitting parameters. It is noted that the infinite shear rate is not seen as a limitation due to rare practical occurrences in an industrial situation. Gangadhar et al. [25] considered the numerical Keller–Box simulation of thermo-diffusion effects on the MHD flow of Casson fluid past a non-flatness stretching surface. The Sutterby fluid flow over a helical channel is documented by [26]. Hayat et al. [27] examined the heat transfer of a Sutterby fluid flow using a rotating disc with chemically reactive species. A numerical solution of the point stagnation flow of the Sutterby fluid in the presence of the Cattaneo–Christov heat flux model is discussed by [28]. Other relevant studies can be found [2932]. In the last decade, attention has been given to boundary-driven flow by a curved stretching sheet. Starting from the fundamental studies in [33], [34], where viscous fluid over curved stretching sheets is investigated, these studies showed that the pressure distribution inside the boundary-layer region could not be neglected compared to the flat sheet. For more related and recent studies, readers should visit [3539], Hayat et al. [40], and [41,42]. In the aforementioned literature, to the best of our knowledge, no consideration has been given to the heat transfer analysis of the Sutterby fluid over a curved stretching sheet. The prime motive of this work is to fill this gap for the completeness of fluid flow and examine the non-Newtonian fluid using curvilinear geometry. However, the Sutterby fluid viscosity is a generalized model that yields shear thinning fluid for m < 0, shear thickening fluid for m > 0, and a Newtonian fluid when m = 0. Therefore, the heat transfer of a magnetohydrodynamic Sutterby fluid past a continuous curved stretched moving surface is investigated in this study. The features of the heat flow are parameterized to prescribe surface temperature, PST, heat flux, and PHF. The flow situations are induced by radiation, viscous dissipation, and frictional heating from the surface mechanism. Evidence of the Keller–Box method enhanced with the Jacobi iterative techniques is numerically implemented. Error analysis is provided to ascertain the convergence of the scheme (see Supplementary Appendix SA1). The velocity field results show that fluid flow can be controlled by virtue of fluid parameters and stretching power index. These observations have profound applications for the flow of the polymer solution and polymer dynamics involving melting and manufacturing of stretchable mechanisms.

2 Basic model

The fundamental equation of the viscosity model is considered as follows [24]:

S=η+η0ηsinh1αΔαΔmA1,(1)

where m, α, Δ, η0, η, and A1 represent the constant index, characteristic time, dynamic viscosity, zero and infinite viscous limit, and first Rivlin–Erickson tensor, respectively. It is worth stating that practically, η does not approach infinity in the form (η) but at an extremely high shear rate. In neglecting this, Eq. 1 reduces to

S=η0sinh1αΔαΔmA1,(2)

such that Δ=trace(A1)2/2 is the symmetry rate of deformation taking in curvilinear geometry (r-,s-). Using Navier–Stokes equation given by

ρDv̄Dt=̄P+̄.S+F̄,(3)

where D/Dt defines the time derivative, ρ is the density of the Sutterby fluid, P represents the pressure, F̄ gives the body force, and ̄=êrr+λêss for which êr and ês are orthogonal radial and axial unit vectors. Using Eqs 2, 3, the steady flow momentum equations can be expressed as

vrv+λsvλRu2=1ρrP+2Srr+Ssrrη02sinh1αΔαΔm,+sinh1αΔαΔmrSrr+λSrrR+λsSsrλSssR,+F̄(4)
vru+λsu+λRuv=λρsP+Srs+2Sssλsη02sinh1αΔαΔm,+sinh1αΔαΔmλrSrs+λSrsr+R+λ2sSss+λSsrr+R,+F̄(5)

where η0 is the viscosity of the Sutterby fluid, and stress components are given as Srr = rv, Srs=Ssr=ru+λsvλRu, Sss=λru+λRv, and λ=Rr+R.

3 Problem description

Consider the steady flow of an electrically conducting Sutterby fluid past a power-law boundary-driven curved surface of radius R. The applied magnetic field B̄(r)=RB0r+Rêr acts upwardly perpendicular to the direction of the fluid flow. The current density J=σ(V×B̄), and neglecting the electric field effect (E ≈ 0), the Lorentz force F̄=J×B̄ can be expressed as

F̄=σRB0r+Ru,0,0,(6)

such that σ is the electrical conductivity of the fluid and B0 is the strength of the applied magnetic field. The velocity vector V = (v, u, 0) gives corresponding velocities v and u in the r − and s − directions. The physical flow geometry is shown in Figure 1.Eqs 4, 5, 6 and the continuity equation yield the following:

Rsur+R+vr+R+rv=0,(7)
vrv+Rsvr+Ru2r+R=rPρ,N14Srrrrv+2Srsrru+Rrsvr+RRsvr+R2,N14SssRrsur+RRsur+R2+rvr+Rvr+R2,+N12Srsrur+Rur+R2,N24Rrvrsvr+R+2RSrssSrsr+R+4RSsssSssr+R,+η0N3ρrr+RSrrr+R+RsSrsr+RSssr+R,(8)
vru+Rsur+R+uvr+R=RsPρr+R,N24Srrrrv+2Srsrru+Rrsvr+RRsvr+R2,N24SssRrsur+RRsur+R2+rvr+Rvr+R2,+N22Srsrur+Rur+R2σR2B02uρr+R2,N24Rrvrsvr+R+2RSrssSrsr+R+4RSsssSrsr+R,+η0N3ρrr+RSsrr+R+RsSssr+R+Ssrr+R,(9)

where N1=mη0α2Srr6ρ, N2=mη0α2Srs6ρ, and N3=(2Srr2+Srs2+2Sss2)(1mα26).

FIGURE 1
www.frontiersin.org

FIGURE 1. Geometry of the physical problem.

Impose boundary conditions relevant to this problem as

uw|r=0=bsn,vw|r=0=0,(10)
u|r=0,ru|r=0,(11)

such that b (t−1l1−n) is a constant and uw and vw are velocities at the surface of the curved sheet. By preserving the flow properties with appropriate scaling, the non-dimensional variables are defined as

ū=uU,v̄=vlUδ,s̄=sl,r̄=rδ,R̄=Rδ,P̄=Pρuw2,η1=η0Ulmα26δ2.(12)

After using Eq. 12, Eqs 79 reduce to the boundary-layer equations:

Rr+Ru2=rP,(13)
vru+Rr+Rsu+uvr+R=Rr+RsP,+N43rrurur+R+ur+R2,+η0ρrru+rur+Rur+R2,σR2B02uρr+R2,(14)

where N4=η1(ruur+R)2 and the continuity equation is identically satisfied.

The similarity variables are introduced in the following form:

η=rδ,ξ=Rδ,uw=bsnfη,P=b2s2P*,(15)
δ=sνbsn+1,H0=σB02b2η0,Res12=bsn+1ν,(16)
v=bRsnr+Rνbsn+1n+12dfdη+ηn12dfdη,(17)

where f(η) represents the stream function, whereas the primes denote the derivative with respect to η. Using Eqs 1517 in Eqs 13, 14, one can get the following:

1η+ξdfdη2=dP*ηdη,(18)
2nξη+ξP*+ηξη+ξn12dP*ηdηd3fdη31η+ξd2fdη2+1η+ξ2dfdη,=ξη+ξn+12fd2fdη2ndfdη2H02ξ2η+ξ2dfdη,mβ0d2fdη21η+ξdfdη23df3dη31η+ξdf2dη2+1η+ξ2dfdη,+ξη+ξ2n+12fdfdη+ηn12dfdη2,(19)

such that H0 and β0=mb4α26 represent the Hartmann number and Sutterby fluid parameter, respectively. Equation 19 in limiting scenario as ξ and P* = 0 yields a generalized Sutterby fluid model over a flat surface:

mβ0d2fdη22d3fdη3=d3fdη3+n+12fd2fdη2ndfdη2H02dfdη.(20)

This equation confirms the present physical models. Numerical solutions have been provided for the Newtonian fluid case, β0 = 0, and the non-Newtonian case, β0 > 0, over a flat stretching sheet [1,5].

Removing the pressure term, after incorporating Eq. 18, Eq. 19 gives

d4fdη4+2η+ξd3fdη3+1η+ξ3dfdη+ξη+ξn+12fd3fdη33n12dfdηd2fdη2,1η+ξ2d2fdη2+ξη+ξ2n+12fd2fdη23n12dfdη2,N53d2fdη2d4fdη43η+ξdfdηd4fdη42η+ξd2fdη2d3fdη3+2η+ξ2dfdηd3fdη3,N53d3fdη32+2η+ξ2d2fdη22+3η+ξ4dfdη2,+4N5η+ξ3dfdηd2fdη2ξη+ξ3n+12fdfdηH02ξ2η+ξ2d2fdη21η+ξdfdη,=0,(21)

subject to

fη|η=0=0,dfdη|η=0=1,dfdη|η==0,d2fdη2|η==0,(22)

where N5=mβ0(d2fdη21η+ξdfdη).

Validating the present model, the second tier over curved stretching mechanism occurs when β0 = 0, Ho = 0, and n ≥ 1 give

d4fdη4+2η+ξd3fdη31η+ξ2d2fdη2+1η+ξ3dfdη,+ξη+ξn+12fd3fdη33n12dfdηd2fdη2,+ξη+ξ2n+12fd2fdη23n12dfdη2,=ξη+ξ3n+12fdfdη.(23)

It is worth mentioning that numerical solutions and homotopy analysis solutions have been presented so far in response to Eq. 23, see ref. [33]

4 Heat transfer analysis

Energy equation incorporating radiation, viscous dissipation, and surface frictional heating can be expressed as

ρCpDTDtK0T=η0sinh1αΔαΔm2Srr2+2Srs2+2Sss2,rQ0σB02R2r+R2u2,(24)

such that T = Tw represents the surface temperature; TT, the ambient temperature; K0, the thermal conductivity of the Sutterby fluid; Cp, the specific heat of the fluid at constant pressure; and Q0, the heat flux. The last term in Eq. 24 contributes surface frictional heating effect. Taking the curvilinear coordinate expansion, this equation can be expressed as

vrT+Rur+RsT=K0ρCprrT+RrTr+R+R2ssTr+R2,+η0Srr2+Srs2+Sss2ρCp1mα262Srr2+Srs2+2Sss2,rQ0ρCp+σB02R2ρCpr+R2u2,(25)

with the condition specify as

PST:T=Tw|r=0=0,T|rT,(26)
PHF:K0Tr|r=0=Bsl2,T|rT.(27)

Imposing the Rosseland approximation for heat flux,

Q0=4σ*3k*T4r,(28)

but the linear radiative flux in Eq. 28 about T gives

Q0r=r4σ*3k*r4T3T3T4.(29)

Using Eq. 29 in Eq. 25, one can get

vrT+Rur+RsT=K0ρCp1+16σ*T23k*K0rrT+K0ρCprTr+R+R2ssTr+R2,+η0ρCpruur+R21mα26ruur+R2,+σB02R2ρCpr+R2u2.(30)

Define similarity transitions of the form

PST:θη=TTTwT,(31)
PHF:K0TT=Dsl2νbsn+1θη.(32)

Applying Eqs 15-17 and Eqs 31, 32 in Eq. 25, one can get

1+Rdd2θdη2+1η+ξdθdη+Prξη+ξn+12fdθdη+ηn12dθdηdfdη,+Brd2fdη21η+ξdfdη21mβ0d2fdη21η+ξdfdη2,=τξ2η+ξ2dfdη2,(33)

where Pr=η0CpK0 is the Prandtl number, Ec=uw2Cp(TwT) the Eckert number, Rd=16σ*T33k*K0 the radiation parameter, τ=σuB02U2nρCp(TwT) the frictional heating parameter, and Br = PrEc the Brinkman number. It should be noted that, for PHF analysis, θ(η) becomes ϑ(η) in Eq. 33.

The corresponding boundary conditions become

PST:θ0|η=0=1,θ|η==0,(34)
PHF:dϑ0dη|η=0=1,ϑ|η==0.(35)

In view of engineering and industrial interest, the physical quantities, such as skin friction coefficient Cf, Nusselt number Nu, and local Nusselt Nu*, are calculated as follows:

Cf=τrs|r=012ρuw2,Nu=sQ0K0TTw,Nu*=sQ0Dsl2,(36)

using τrs=η0A1[sinh1(αΔ)αΔ]m(ruur+R)|r=0, Q0 = −K0rT|r=0, and uw = bsn,and these yield

Res12Cf=d2f0dη21ξ1β0d2f0dη21ξ2,(37)
Res12Nu=dθ0dη,Res12Nu*=ϑ0.(38)

Eq. 38 gives the surface rate of the heat transfer for PST and PHF, respectively.

5 Numerical computations

The numerical technique to solve the higher order nonlinear coupled differential Eqs 21, 33 subjected to Eqs 22, 34, 35 is considered in this section. Initialize the boundary value problem (BVP) into the initial value problem (IVP) using the variables

f,dfdη,d2fdη2,d3fdη3,θ,dθdη,ϑ,dϑdηT=y,y1,y2,y3,y4,y5T.(39)

Linearizing into the system of first-order equations can be expressed as

yy1y2y3y4y5=y1y2y32λy3+λ2y2λ3y1λξm+12yy33m12y1y2+ξλ3m+12yy2λ2ξm+12yy23m12y12+mβy2λy13y2y3+3λy1y32λy2y3+2λ2y1y3+2λ3y32+mβy2λy13λ4y22+3y324λ3y1y2H02y2ξ2y2λy1y511+Rdλy6+λξPrm+12yy5+ηm12y1y5+τξ2λ2y12Br1+Rdy2λy121βy2λy1,(40)

subject to the boundary conditions

y10y20y30y40y50y60=01c1c21c3,(41)

where λ=1η+ξ and c1, c2, and c3 are missing initial conditions. It is worth noting that θ and ϑ are interchanged to obtain the solution of the PHF case.

Implementing the standard Keller–Box shooting method and the Jacobi iterative technique [18], let

w=w1,w1=w2,w2=w3,,wn1=Hξ,w1,w2,w3,,wn1,(42)
win+1wi1nδh=w1i1/2n,
w1in+1w1i1nδh=w2i1/2n,
,
wn1in+1wn1i1nδh=H2ξi1/2n,wi1/2n,w1i1/2n,,wn1i1/2n.(43)

Equation 43 is explicitly implemented into Eq. 40 subject to conditions Eq. 41 as follows:

yin+1yi1nδh+y1in+y1i1n2=0,(44)
y0n+1=0,(45)
y1in+1y1i1nδh+y2in+y2i1n2=0,(46)
y10n+1=1,(47)
y2in+1y2i1nδh+y3in+y3i1n2=0,(48)
y20n+1=c1,(49)
y3in+1δh=y3i1nδhy3i+y3i1η+ξN5y1i+y1i1y2i+y2i1η+ξ3n+1yi+yi1y3i+y3i1ξ8η+ξy1i+y1i12η+ξ3+3n1y1i+y1i1y2i+y2i1ξ8η+ξ+y2i+y2i12η+ξ2n+1yi+yi1y2i+y2i1ξ8η+ξ2+3n1y1i+y1i12ξ8η+ξ2+n+1yi+yi1y1i+y1i1ξ8η+ξ3+H02ξ2y2i+y2i12η+ξ2+3N5y2i+y2i1y3in+1y3i1n4δh+3N5y3i+y3i1243N5y1i+y1i1y3in+1y3i1n4δhη+ξ+N5y2i+y2i122η+ξ2N5y2i+y2i1y3i+y3i12η+ξ+3N5y1i+y1i124η+ξ4+N5y1i+y1i1y3i+y3i12η+ξ2H02ξ2y1i+y1i12η+ξ3,(50)
y30n+1=c2,(51)
y4in+1y4i1nδh+y51in+y51i1n2=0,(52)
y40n+1=1,(53)
y5in+1y5i1nδh=y5i+y5i121+Rdη+ξτξ2y1i+y1i12η+ξ2n+1yi+yi1y5i+y5i1Prξ8η+ξηn1y1i+y1i1y5i+y5i1Prξ8η+ξBrN52m2β2+BrN54m3β3,(54)
y50n+1=c3,(55)

where c1, c2, and c3 are initial guesses obtained through MATLAB fsolve. The accuracy of the initial guesses is checked against endpoint ηmax. Otherwise, the values are updated, and iterations continue. A step size, δh = 10–2, is used until the solution is obtained, which satisfies the endpoint initially at infinity and is subject to the stopping criterion given as

|wn+12wn2|<ϵ,(56)

such that ϵ remains the fixed tolerance value of 10–8 to ascertain the numerical results. The algorithm is implemented through the ode45 built-in MATLAB command.

6 Results and graphical analysis

The behaviors of the field quantities, velocity, and temperature against various flow parameters are presented in this section. Fluid parameters characterized in this problem include the radius of curvature, ξ, magnetic parameter, Ho, Sutterby fluid parameters m and β0, and stretching power index n. The behavior of temperature distribution is examined against Prandtl number Pr, Brinkman number Br, magnetic parameter Ho, radiation parameter Rd, and surface frictional heating parameter τ. Flow stream function is plotted in Figures 26 against the radius of curvature, Lorentz force, Sutterby fluid parameters, and stretching power index, respectively. In these graphs, the enhancement of any flow parameter decreases the flow trajectories. Physically, this can be seen as a means of controlling the flow field and boundary-layer region. Figures 7, 8 show that the velocity decreases for increasing ξ and Ho and consequently reduces the boundary-layer thickness. Thus, the flow field and flow trajectories can be significantly controlled via the size of these parameters. The impacts of Sutterby fluid parameters, β0 and m, on the velocity profile are shown in Figures 9, 10, respectively. The graphs infer that velocity profiles decrease for optimizing the fluid parameters. Physically, this observation shows that increasing fluid parameters enhances the viscous resistance (shear-thickening) property of the fluid and invariably controls the flow velocity and associated boundary-layer thickness.

FIGURE 2
www.frontiersin.org

FIGURE 2. Effects of ξ on f(η).

FIGURE 3
www.frontiersin.org

FIGURE 3. Effects of Ho on f(η).

FIGURE 4
www.frontiersin.org

FIGURE 4. Effects of β0 on f(η).

FIGURE 5
www.frontiersin.org

FIGURE 5. Effects of m on f(η).

FIGURE 6
www.frontiersin.org

FIGURE 6. Effects of n on f(η).

FIGURE 7
www.frontiersin.org

FIGURE 7. Effects of ξ on f′(η).

FIGURE 8
www.frontiersin.org

FIGURE 8. Effects of Ho on f′(η).

FIGURE 9
www.frontiersin.org

FIGURE 9. Effects of β0 on f′(η).

FIGURE 10
www.frontiersin.org

FIGURE 10. Effects of m on f′(η).

Figure 11 explains the effect of stretching power n on the velocity distribution. The graph shows a decrease in velocity for large n. Thus, it confirms that nonlinear stretching velocity for the boundary-driven flow could serve as a means of controlling the flow field. This observation is significant and can play a vital role in flow past stretchable materials. The effect of thermal kinetics parameters Pr, Br, Ho, Rn, τ, and ξ on both PST, θ(η), and PHF, ϑ(η) scenarios are presented in Figures 1217. Figure 12 shows the impact of the Prandtl number Pr on the temperature distribution, revealing that the temperature region with the thermal boundary-layer thickness decay significantly. In other words, it is caused by the increased thermal conductivity of the Sutterby fluid. Furthermore, this observation is prominently seen in PHF, meaning that heat flow can be enhanced by virtue of the prescribed heat flux than surface temperature. Figure 13 elucidates the effect of the Brinkman number Br on the temperature distribution. In this graph, the temperature region is expanded by increasing Br in the PST and PHF cases. However, it infers that large Br weakens the fluid’s viscous bond and generates additional heat from the fluid. Furthermore, it is noted that the heat occurrence is more intense in the PHF than in the PST region. An increase in the temperature distribution and the thermal boundary-layer thickness for higher magnetic number Ho is witnessed in Figure 14. The reason is that quantum heat is produced as the flow velocity is disrupted by the opposing large Lorentz force and then enlarges the temperature region. Figure 15 presents an optimizing effect of increasing the radiation parameter Rd on the temperature profile. This observation agrees with the usual and existing fact that radiation serves as an additional source of heat generation. Thus, the increase is more rapid for the PHF scenario compared to the PST region. Increasing the surface frictional heating parameter, τ optimizes the temperature profile in both heat situations, as shown in Figure 16. This impact is caused by surface induction and adding more heat from the solid surface to the fluid. However, a further effect is displayed in the PHF than in the PST profile. The impact of the dimensionless radius of curvature on the temperature profile is given in Figure 17. The temperature distribution with associated thermal boundary-layer thickness is enhanced to increase the curvature (reducing the radius of curvature). This observation concludes that heat travels faster over a curved surface compared to a flat surface and is more profound in the PHF than in the PST cases. It shows an important result and a substantial contribution to this study.

FIGURE 11
www.frontiersin.org

FIGURE 11. Effects of n on f′(η).

FIGURE 12
www.frontiersin.org

FIGURE 12. Effects of Pr on temperature.

FIGURE 13
www.frontiersin.org

FIGURE 13. Effects of Br on temperature.

FIGURE 14
www.frontiersin.org

FIGURE 14. Effects of Ho on temperature.

FIGURE 15
www.frontiersin.org

FIGURE 15. Effects of Rd on temperature.

FIGURE 16
www.frontiersin.org

FIGURE 16. Effects of τ on temperature.

FIGURE 17
www.frontiersin.org

FIGURE 17. Effects of ξ on temperature.

The behavior of pressure is significant to the flow field over the curved structure; therefore, pressure variations need to be discussed. Figures 1821 depict the pressure variations due to the Lorentz force, stretching the power index, Sutterby fluid parameter, and radius of curvature. Clearly, Figure 18 elucidates a diminished pressure profile for a large Lorentz force. This observation confirms the earlier behavior of flow trajectory and velocity field (see Figures 3, 8). The gradual decrease in pressure profiles is shown in Figures 1821. These observations support the graphs in Figures 8, 9, 11 and further accentuate flow behaviors over the curved stretching mechanism. In other words, fluid flow can be regulated and controlled by means of the stretching power index. Likewise, Figure 20 elucidates a slight decrease in the flow pressure for a large fluid parameter, which shows a powerful influence of the fluid parameter for regulating flow velocity. It further substantiates the observations made in Figures 4, 9. Figure 21 is plotted to authenticate the present model; see Eq. 21. In validation of the present investigation, we set the pressure to zero, P(η) = 0 and the large dimensionless radius of curvature, ξ (surface becomes flat). A valuable known result for fluid flow over the flat stretching surface is achieved [23–25]. Table 1 shows an excellent agreement of the computational values for surface drags compared to the published studies. Table 2 is provided in the remark of the surface rate of heat transfer for the PST and PHF cases. It also presents the impact of the magnetic field on skin friction coefficient and rate of heat transfer for PST and PHF cases.

FIGURE 18
www.frontiersin.org

FIGURE 18. Effects of Ho on P(η).

FIGURE 19
www.frontiersin.org

FIGURE 19. Effects of n on P(η).

FIGURE 20
www.frontiersin.org

FIGURE 20. Effects of β0 on P(η).

FIGURE 21
www.frontiersin.org

FIGURE 21. Effects of ξ on P(η).

TABLE 1
www.frontiersin.org

TABLE 1. Comparison of numerical results ξ with published results.

TABLE 2
www.frontiersin.org

TABLE 2. Numerical values of the surface rate of heat transfer for various ξ, β0, Ho, Pr, Br, Rd, and τ when m =0.2.

7 Concluding remark

The magnetohydrodynamic flow of the non-Newtonian Sutterby fluid driven by nonlinear power-law velocity over a continuous stretching curved mechanism has been addressed. The analysis is presented in a radially varying magnetic field acting perpendicular to the flow direction. Viscous dissipation, surface frictional heating, and radiation serve as the additional source of heat for prescribed surface temperature, PST, and prescribed heat flux, PHF. Governing models are theoretically formulated under desirable physical phenomena to ensure realistic flow behaviors. Computation is performed using the standard Keller–Box method and the Jacobi iterative scheme. The following results were plausibly recorded and summarized.

1. The velocity field decreases to increase Sutterby fluid parameters and stretch the nonlinear power index, whereas the opposite effect occurs to enhance the magnetic field parameter and dimensionless radius of curvature.

2. The temperature distribution increases for large Brinkman number, radiation, and surface frictional heating parameters.

3. Increasing magnetic field generates a large Lorenz force, which enhances the temperature profile.

4. Additional heat is transferred over the curved surface compared to the flat surface, and it is more significant in the PHF than in the PST regions.

5. Heat transfer through the prescribed heat flux is more important to enhance heat flow due to the additional heat generated from the solid surface to the fluid.

6. The results are validated by comparing the surface drags with the published work, authenticating our numerical method and concluding by presenting numerical values for the heat transfer rate in the PST and PHF cases.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.

Author contributions

All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.

Funding

The research was funded by the University of Oradea, Romania.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphy.2023.1126003/full#supplementary-material

Abbreviations

η0, zero viscous limit; η, infinite viscous limit; m, fluid behavioral index; △, dynamic viscosity; ρ, fluid density; R, radius of curvature; σ, fluid electrical conductivity; B, magnetic field; B0, applied magnetic field strength; (v, u), velocity in radial and axial directions; Ho, Hartmann number; β0, Sutterby fluid parameter; ξ, dimensionless curvature parameter; K0, fluid thermal conductivity; Cp, fluid specific heat capacity; θ(η), dimensionless temperature for PST; ϑ(η), dimensionless temperature for PHF.

References

1. Jalil M, Asghar S, Mushtaq M. Analytical solutions of the boundary layer flow of power-law fluid over a power-law stretching surface. Commun Nonlinear Sci Numer Simulat (2013) 18:1143–50. doi:10.1016/j.cnsns.2012.09.030

CrossRef Full Text | Google Scholar

2. Patil VS, Patil NS, Timol MG. A remark on similarity analysis of boundary layer equations of a class of non-Newtonian fluids. Int J Nonlinear Mech (2015) 71:127–31. doi:10.1016/j.ijnonlinmec.2014.10.022

CrossRef Full Text | Google Scholar

3. Anuar NS, Bachok N, Mustafa T, Md Arifin N, Haliza R. Analytical and stability analysis of MHD flow past a nonlinearly deforming vertical surface in Carbon Nanotubes. A Eng J (2020) 59(1):497–507. doi:10.1016/j.aej.2020.01.024

CrossRef Full Text | Google Scholar

4. Kumaran G, Sivaraj R, Reddy AS, Kumar R, Prasad VR. Hydromagnetic forced convective flow of Carreau nanofluid over a wedge plate/stagnation of the plate. Eur Phys J Spec Top (2019) 228:2647–59. doi:10.1140/epjst/e2019-900069-2

CrossRef Full Text | Google Scholar

5. Jalil M, Asghar S. Flow of power-law fluid over a stretching surface: A lie group analysis. Int J Nonlinear Mech (2013) 48:65–71. doi:10.1016/j.ijnonlinmec.2012.07.004

CrossRef Full Text | Google Scholar

6. Wahid NS, Md Arifin N, Turkyilmazoglu M, Hafidz MEH, Abd Rahmin NA. MHD hybrid Cu-Al&lt;sub&gt;2&lt;/sub&gt;O&lt;sub&gt;3&lt;/sub&gt;/water nanofluid flow with thermal radiation and partial slip past a permeable stretching surface: Analytical solution. J Nano Res (2020) 64:75–91. doi:10.4028/www.scientific.net/jnanor.64.75

CrossRef Full Text | Google Scholar

7. Hamad AA. Analytical solution of natural convection flow of a nanofluid over a linearly stretching sheet in the presence of magnetic field. Int Commun Heat Mass (2011) 38:487–92. doi:10.1016/j.icheatmasstransfer.2010.12.042

CrossRef Full Text | Google Scholar

8. Gorla RS, Sidawi I. Free convection on a vertical stretching surface with suction and blowing. Appl Sci Res (1994) 52:247–57. doi:10.1007/bf00853952

CrossRef Full Text | Google Scholar

9. Mustaq A, Mustafa M, Hayat T, Alsaedi A Effect of thermal Radiation on the stagnation-point flow of upper-convected Maxwell fluid over a stretching sheet. J Aerosp Eng, 27(4), 1–6. (2014). doi:10.1061/(ASCE)AS.1943-5525.0000361

CrossRef Full Text | Google Scholar

10. Jawadi A, Boutyour H, Cadou JM. Asymptotic numerical method for steady flow of power-law fluids. J Non-newton Fluid Mech (2013) 202:22–31. doi:10.1016/j.jnnfm.2013.09.005

CrossRef Full Text | Google Scholar

11. Zierep J, Fetecau C. Energetic balance for the Rayleigh-Stokes problem of a Maxwell fluid. Int J Eng Sci (2007) 45:617–27. doi:10.1016/j.ijengsci.2007.04.015

CrossRef Full Text | Google Scholar

12. Hayat T, Ijaz MK, Tamoor M, Waqas M, Alsaedi A. Numerical Simulation of heat transfer in MHD Stagnation point flow of Cross fluid model towards a stretched surface. Results Phys (2017) 7:1824–7. doi:10.1016/j.rinp.2017.05.022

CrossRef Full Text | Google Scholar

13. Ijaz MK, Waqas W, Hayat T, Alsaedi A. Magneto-hydrodynamical numerical simulation of heat transfer in MHD stagnation point flow of Cross fluid model towards a stretched surface. Phys Chem Liq (2017) 56(5):584–95. doi:10.1080/00319104.2017.1367791

CrossRef Full Text | Google Scholar

14. Manzur M, Khan M, Rahman M. Mixed convection heat transfer to Cross fluid with thermal radiation: Effects of buoyancy assisting and opposing flows. Int J Mech Sci (2018) 138-139:515–23. doi:10.1016/j.ijmecsci.2018.02.010

CrossRef Full Text | Google Scholar

15. 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–72. doi:10.1016/j.rinp.2017.08.039

CrossRef Full Text | Google Scholar

16. Zheng LC, Wang KN, Gao YT. Unsteady flow and heat transfer of a generalized Maxwell fluid due to a hyperbolic sine accelerating plate. Comp Math Appl (2011) 61:2209–12. doi:10.1016/j.camwa.2010.09.017

CrossRef Full Text | Google Scholar

17. Rout BC, Mishra SR. Thermal energy transport on MHD nanofluid flow over a stretching surface: A comparative study. Eng Sci Tech Int J (2018) 21(1):60–9. doi:10.1016/j.jestch.2018.02.007

CrossRef Full Text | Google Scholar

18. Yasir N, Shoaib Arif M. Keller-Box shooting method and its application to nanofluid flow over convectively heated sheet with stability and convergence. Num Heat Transfer, B: Fundamentals (2019) 76:152–80. doi:10.1080/10407790.2019.1644924

CrossRef Full Text | Google Scholar

19. Mohamed R, Eid Kasseb ML, Mohammed T, Sheikholeslami M. Numerical treatment for Carreau nanofluid flow over a porous nonlinear stretching surface. Results Phys (2018) 8:1185–93. doi:10.1016/j.rinp.2018.01.070

CrossRef Full Text | Google Scholar

20. Prasannakumara BC, Gireesha BJ, Krishnamurthy MR, Kumar KG. MHD flow and nonlinear radiative heat transfer of Sisko nanofluid over a nonlinear stretching sheet. Inf Med Unlocked (2017) 9:123–32. doi:10.1016/j.imu.2017.07.006

CrossRef Full Text | Google Scholar

21. Feroz AS, Usman M, Ul Haq R, Wang W. Melting heat transfer analysis of Sisko fluid over a moving surface with nonlinear thermal radiation via collocation method. Heat Mass Transfer (2018) 126:1034–42. doi:10.1016/j.ijheatmasstransfer.2018.05.099

CrossRef Full Text | Google Scholar

22. Lund LA, Omar Z, llyas K. Analysis of dual solution for MHD flow of Williamson fluid with Slippage. Heliyon (2019) 5(3):e01345–20. doi:10.1016/j.heliyon.2019.e01345

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Khalil R, Qaiser A, Malik MY, Ali U. Numerical communication for MHD thermally stratified dual convection flow of Casson fluid yields by stretching cylinder. Chin J Phys (2017) 55(4):1605–14. doi:10.1016/j.cjph.2017.05.002

CrossRef Full Text | Google Scholar

24. Sutterby JL. Laminar converging flow of dilute polymer solutions in conical sections: Part I. Viscosity data, new viscosity model, tube flow solution. A Ch E J (1965) 12(1):63–8. doi:10.1002/aic.690120114

CrossRef Full Text | Google Scholar

25. Gangadhar K, Kannan T, Narasimharao NSIV, Kumar BR, Sakthivel G. Thermo diffusion effects on MHD Casson fluid flow over non-flatness stretching surface: Keller Box method. Int J Amb Energ (2021) 42(4):374–82. doi:10.1080/01430750.2018.1531263

CrossRef Full Text | Google Scholar

26. Batra RL, Eissa M. Helical flow of a sutterby model fluid. Polym-plast Technol Eng (1994) 33(4):489–501. doi:10.1080/03602559408010743

CrossRef Full Text | Google Scholar

27. Hayat T, Sidra A, Ijaz MK, Alsaedi A. Irreversibility aspect to flow of sutterby fluid subject to nonlinear heat flux and Ohmic heating. App Nanosci (2019) 9:1215–26. doi:10.1007/s13204-019-01015-3

CrossRef Full Text | Google Scholar

28. Azhar E, Iqbal Z, Ijaz S, Maraj EN. Numerical approach for stagnation point flow of sutterby fluid impinging to Cattaneo-Christov heat flux model. Pramana J Phys (2018) 91:61. doi:10.1007/s12043-018-1640-z

CrossRef Full Text | Google Scholar

29. Hayat T, Ahmad S, Ijaz MK, Alsaedi A. Modeling chemically reactive flow of sutterby nanofluid by a rotating disk in presence of heat generation/absorption. Commun Theo Phys (2018) 69:569. doi:10.1088/0253-6102/69/5/569

CrossRef Full Text | Google Scholar

30. Sher NA. Biomathematical study of sutterby fluid model for blood flow in stenosed arteries. Int J Biomaths (2015) 8(6):1550075. doi:10.1142/s1793524515500758

CrossRef Full Text | Google Scholar

31. Ahmed S, Farooq M, Javed M, Aisha A. Double stratification effects in chemically reactive squeezed sutterby fluid flow with thermal radiation and mixed convection. Results Phys (2018) 8:1250–9. doi:10.1016/j.rinp.2018.01.043

CrossRef Full Text | Google Scholar

32. Abbas Z, Shabbir MS, Ali N. Numerical Study of magnetohydrodynamic pulsatile flow of sutterby fluid through an inclined overlapping arterial stenosis in the presence of periodic body acceleration. Results Phys (2018) 9:753–62. doi:10.1016/j.rinp.2018.03.020

CrossRef Full Text | Google Scholar

33. Sajid M, Ali N, Javed T, Abbas Z. Stretching a curved surface in a viscous fluid. Chin Phy Lett (2010) 27:024703. doi:10.1088/0256-307x/27/2/024703

CrossRef Full Text | Google Scholar

34. Sanni KM, Asghar S, Jalil M, Okechi NF. Flow of viscous fluid along a nonlinearly stretching curved surface. Results Phys (2017) 7:1–4. doi:10.1016/j.rinp.2016.11.058

CrossRef Full Text | Google Scholar

35. Sanni KM, Hussain Q, Asghar S. Heat transfer analysis for non-linear boundary driven flow over a curved stretching sheet with a variable magnetic field. Front Phys (2020) 8:113. doi:10.3389/fphy.2020.00113

CrossRef Full Text | Google Scholar

36. Naveed M, Abbas Z, Sajid M. Hydromagnetic flow over an unsteady curved stretching surfce. Eng Sci Tech Int J (2016) 19:841–5. doi:10.1016/j.jestch.2015.11.009

CrossRef Full Text | Google Scholar

37. Maria I, Hayat T, Alsaedi A, Hobiny A. Homogeneous-heterogeneous reaction in MHD flow due to an unsteady curved stretching surface. J Mol Liq (2016) 221:245–53. doi:10.1016/j.molliq.2016.05.060

CrossRef Full Text | Google Scholar

38. Abbas Z, Naveed M, Sajid M. Hydromagnetic Slip flow of nanofluid over a Curved Surface with heat generation and thermal radiation. J Mol Liq (2016) 215:756–62. doi:10.1016/j.molliq.2016.01.012

CrossRef Full Text | Google Scholar

39. Sanni KM, Asghar S, Saima R, Wu YC. Nonlinear radiative treatment of hydromagnetic non-Newtonian fluid flow induced by nonlinear convective boundary driven curved sheet with dissipations and chemical reaction effects. Front Phys (2021) 1–15. doi:10.3389/fphy.2021.670930

CrossRef Full Text | Google Scholar

40. Hayat T, Sajjad R, Ellahi R, Alsaedi A, Muhammad T. Homogeneous-heterogeneous reaction in MHD flow of micropolar fluid by a curved stretching surface. J Mol Liq (2017) 240:209–20. doi:10.1016/j.molliq.2017.05.054

CrossRef Full Text | Google Scholar

41. Sanni KM, Hussain Q, Asghar S. Flow of magnetohydrodynamic viscous fluid by curved configuration with non-linear boundary driven velocity. J Taibah Uni Sci (2021) 15(1):589–98. doi:10.1080/16583655.2021.1991076

CrossRef Full Text | Google Scholar

42. Sanni M, Hussain Q, Asghar S. Thermal analysis of a hydromagnetic viscoelastic fluid over a continuous curved stretching surface in the presence of radiative heat flux. Arab J Sci Eng (2020) 46:1–8. doi:10.1007/s13369-020-04671-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: nonlinear, MHD, curved sheet, radiation, viscous dissipation, power-law, Keller-Box, error analysis

Citation: Sanni KM, Asghar S, Al-Shbeil I and Cătaş A (2023) Radiative simulation of non-Newtonian MHD fluid over a boundary-driven multi-physical curved mechanism: Keller–Box evidence. Front. Phys. 11:1126003. doi: 10.3389/fphy.2023.1126003

Received: 16 December 2022; Accepted: 15 March 2023;
Published: 12 April 2023.

Edited by:

Stefano Leonardi, The University of Texas at Dallas, United States

Reviewed by:

Mustafa Turkyilmazoglu, Hacettepe University, Türkiye
B. Rushi Kumar, VIT University, India

Copyright © 2023 Sanni, Asghar, Al-Shbeil and Cătaş. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Adriana Cătaş, acatas@gmail.com

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.