Thermal analysis for [Formula: see text]-sodium alginate magnetized Jeffrey's nanofluid flow past a stretching sheet embedded in a porous medium.

The magnetohydrodynamics (MHD) viscous Jeffrey heat transport flow past a permeable extending sheet is analyzed. The Alumina ([Formula: see text]) is chosen as nanoparticles immersed in sodium alginate ([Formula: see text]) as the based fluid. The effect of heat generation, Ohmic heating and viscous dissipation are also being investigated adopting Tiwari and Das model. The adequate similarity transformation is used to convert the governing equations to non-linear of higher-order ordinary differential equations (ODEs). The numerical solution of the transformed ODEs is accomplished using a finite-difference technique. The results are described in graphs according to selected parameters' values provided. The flow velocity reductions when the porosity parameter is augmented. The thermal distribution is affected by the presence of [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text] and [Formula: see text]. Deborah number and the volume fraction of nanoparticles affect the skin friction coefficient in opposite ways. A higher volume percentage of nanoparticles and a higher Deborah number are both shown to boost the rate of heat transfer. These findings suggest that the concentration of nanoparticles can be used to manipulate heat transport and nanofluid motions.

Introduction

Before nanotechnology was introduced, the reduced thermal conductivity of several frequently used liquids, such as water, ethylene glycol, kerosene oil, engine oil, has been a challenge for research workers and engineers. The inclusion of a suitable quantity of nanosized particles was proposed by Choi and Eastman[1] to solve the problem. Choi’s investigations paved the way for advancing and developing techniques to improve conventional fluids’ thermal properties, known as base fluids. Nanofluids keep acquired a lot of consideration in recent years because of their applications[2] in various fields. This application includes oil recovery, transportation, microfluidics, medical, microelectronics, fuel cells, and manufacturing. Some parameters that can contribute to improving the thermal conductivity exist, such as the Brownian motion effects (see Jang and Choi[3]), as well as nanoparticle size and temperature (see Ref.[4]). Elbashbeshy et al.[5] have looked at external forces’ contribution and coolant (nanofluid) on rolling cylinders’ mechanical properties. Kameswaran et al.[6] have undertaken work on magnetohydrodynamic (MHD) nanofluid flows across a stretching and shrinking surface with chemical reaction and viscous expansion. Ganga et al.[7] studied a nano liquid MHD’s behaviour due to stretching surface with heat generation/absorption effects. Seth and Mishra[8] have inspected MHD influence on fluid embedded with nanoparticles moving subjected to non-linear expandable surface.

Researchers are currently very interested in everything concerning non-Newtonian liquids. Indeed, these fluids are of considerable importance in various organic and natural, mechanical, and design processes, such as glass design, fibre sheet formation, subsistence items, wire drawing, paper creation, gemstone development, etc. Among the different types of non-Newtonian liquids, Jeffrey fluid model is presented mainly because of a less complicated linear model in which time derivatives are utilized rather than convective derivatives. Xu and Liao[9] envisioned an analytical solution of time-dependent MHD non-Newtonian flow resulting from an impulse-stretching wall surface. Hung[10] carried out an analysis of entropy production using non-Newtonian fluid within microchannels and viscous dissipation. Kamali and Binesh[11] provided a computational study on enhancing the heat transmission using non-Newtonian nanofluids of carbon nanotubes (CNTs) type. They found that the Nusselt number of CNT nanofluid was more remarkable than that of the base fluid. Rundora and Makinde[12] reviewed how aspiration/injection affects a 3rd-degree fluid with a variable viscosity with an irregular reaction within a channel containing a porous matrix. Sheikholeslami et al.[13] worked on the MHD flux of nanofluids together with heat transfer and the thermal radiation effect using the Buongiorno model. Eldabe and Abou-zeid[14] employed the homotopy perturbation technique to solve the magneto-nanofluid flow problem within a porous matrix. Shahsavani et al.[15] used the experimental results to calculate the heat transfer and pressure fall of nano liquid flow in a cylinder. By the lattice Boltzmann technique, Krishna and Reddy[16] simulated the magnetohydrodynamic forced convective flow due to stumpy spongy, porous media.

Additionally, boundary layer flows with the stretched surface have numerous applications, including paper production, wire drawing, cooling of a continuous strip, and polymers’ extrusions. Fluid dynamics flow due to an expandable medium was coined by Crane[17]. Bhatnagar et al.[18] examined non-Newtonian flow produced by a free-flowing stretched sheet. Raptis and Chemical reaction and electromagnetic force were both studied by Perdikis[19] to see how they affected the flow of viscous fluids across non-linearly stretched surfaces. Thermal and mass transfer parameters of MHD viscous flow generated by porous stretching surfaces were studied by Turkyilmazoglu[20]. Cattaneo–Christov thermal flux impact on viscoelastic MHD fluid flows and heat transfer across a vertical stretching surface were of interest to Li et al.[21]. After that, many studies appear by considering various aspects of the problem[22-29] to list only a few.

Magnetohydrodynamics refers to the fluid mechanics of the electrically conductive media having an interaction with magnetic fields. Its application in geothermal energy extraction, nuclear reactors, MHD pumps, blood flow measurements, etc. Attia and Kotb[30] examined the steady-state flow of an incompressible fluid, electrically conductive, flowing between two parallel plates and subjected to an external magnetic force. The effects of thermo-capillary and buoyancy on the movement of an electrical insulator liquid in a rectangular duct when internal heat production is taken into account were studied by Hossain et al.[31]. Malekzadeh et al.[32] investigate the effects of varying viscosity and radiation on MHD fluid velocity over a moving vertical surface. For the first time, Ghasemi et al.[33] have quantitatively investigated the effects of sun radiation on MHD flow. An inclined stretched porous sheet with unstable mixed convection of nanofluid flow is studied by Jain et al.[34]. An electrically conducting nanofluid travelling towards an extensible surface was proposed by Nayak et al.[35]. Abel[36] scrutinized the effect of viscous dissipation on fluid past an extending surface and found that fluid temperature amplifies in the case of enrichment in the heat dissipation parameter. Kishan and Deepa[37] contemplated micropolar fluid flow along an expandable sheet having pores embedded with heat dissipation, stagnation point as well as viscous dissipation effects. Alim et al.[38] pondered Newtonian-based fluid embedded with results like Ohmic and viscous dissipation over a vertical expandable sheet surface. Ferdows et al.[39] were interested in viscous dissipation effects associated with Hall current on liquid flow through an expanding medium.

A Casson fluid is a shear thinning liquid with an infinite viscosity at zero rate of shear and a zero viscosity at an infinite rate of shear[40]. Casson fluid is commonly found in honey, jelly, soup, tomato sauce, concentrated fruit liquids, and other foods. It is also the best rheological model for blood and chocolate[41]. Casson fluid also has yield stress and is very important in polymer processing industries and biomechanics[42]. The Jeffrey fluid model can represent the stress relaxation property of non-Newtonian fluids, whereas the traditional viscous fluid model cannot. The Jeffrey fluid model well describes a class of non-Newtonian fluids with the distinctive memory time scale, also known as the relaxation time[43]. Relaxation time is the period of time when a system relaxes in response to changes in external variables[44,45]. It refers to the time it takes for a polymer coil to relax from a distorted condition to its equilibrium configuration. It is an important parameter in determining the properties of a viscoelastic fluid. Non-Newtonian fluids include drilling muds, apple sauce, foams, soaps, sugar solution pastes, clay coating, ketchup, lubricant, some oils, colloidal and suspension solutions. Non-Newtonian fluids are classified into three types: differential, integral, and rate.

The goal of the analysis is to assess the boundary-layer –sodium alginate magnetized viscoelastic Jeffrey’s nanofluid flow due to a continuous surface. Moreover, this type of work for the non-Newtonian fluid model is even more narrow down. This work is primarily concerned with investigating the heat transfer flow of Jeffrey nanomaterial across a flexible surface. For the following reasons, the present investigation has a discernible novelty in it. (i) the incorporation of the mono nanoparticles with sodium alginate-based fluid, (ii) the embodiment of magnetohydrodynamics effect, (iii) the analysis of the heat transfer flow in the light of the Ohmic heating, viscous heat dissipation, porous medium and heat source and (iv) the utilization of Tiwari and Das[46] model with the thermophysical attributes of nanofluids. The modified ODEs were solved digitally using MATLAB software and the finite difference method. The stimulus of dimensionless factors on velocity and heat fields was studied in detail with graphs. In addition, the impact of various dimensionless factors on drag coefficient and local Nusselt number was investigated.

Mathematical modelling

We deliberate an –sodium alginate Jeffrey nanofluid laminarly and electrically conducting flow through a stretchable sheet. The nanofluid experiencing a magnetic field employed in the -axis direction, and flow is bounded in , as presented in Fig. 1. The stretchable sheet is along the -axis path, and the flow is known as linearly stretching with velocity such that is a positive constant. The plate wall temperature is known as and assumed to have a quadratic type expansion at , i.e., to administer the nanofluid flow. The induced magnetic field impact is ignored, which is justifiable in the small magnetic Reynolds number.

Figure 1

Physical interpretation of flow geometry.

Subject to these assumptions, the boundary layer equations administering the nanofluid flow and the heat fields can be presented in the dimensional form given below

The corresponding endpoint conditions are:

The thermophysical characteristics of nanofluids are described as follows (see Tables 1 and 2)[46,47]:

Table 1

Thermophysical characteristics of Jeffrey nanofluid.

Properties	Nanofluid
Dynamic viscosity (μ)	μ_nf = μ_f(1-ϕ)^-2.5
Density (ρ)	ρ_nf = (1-ϕ)ρ_f + ϕρ_s
Heat capacity (ρC_p)	(ρC_p)_nf = (1-ϕ)(ρC_p)_f + ϕ(ρC_p)_s
Electrical Conductivity (σ)	σ_nfσ_f=1+3(σ_s-σ_f)ϕ(σ_s+2σ_f)-(σ_s-σ_f)ϕ
Thermal Conductivity (κ)	κ_nfκ_f = (κ_s+2κ_f)-2ϕ(κ_f-κ_s)(κ_s+2κ_f)+ϕ(κ_f-κ_s)

Table 2

Values of thermophysical features[48].

Thermophysical	ρ(kg/m^3)	c_p(J/kgK)	k(W/mK)	σ(S/m)	Pr
Sodium alginate (SA)	989	4175	0.6376	2.6×10^-4	6.5
Alumina (Al_2O_3)	3970	765.0	40.000	2.7×10^-1	–

Initiate the following similarity variables, which transform Eqs. (1)–(3) into the ODEsand stream function is expressed as follows

So we have

Employing Eqs. (5) and (7), Eqs. (1)–(3) are turned intosubjected to BCs

Various dimensionless parameters occurring in Eqs. (8) and (9) are delineated as

The expressions regarding substantial surface drag and heat transfer , are manifested byexpressions regarding shear as well as heat flux at the wall are manifested by and respectively and moreover Eq. (12) can be transformed likewisewhereas the expression regarding Reynolds number is .

Solution methodology

The dimensionless ODEs (8)–(9) along with Eq. (10) can be solved using a reliable numerical technique named the Keller box method[49,50] with the aid of MATLAB software. The flow chart (see Fig. 2) mechanisms of this scheme are displayed below to achieve numerical outcomes:

Figure 2

Flow chart illustrating the Keller box method.

Numerical procedure

We introduce dependent variables , , and such that

So that Eqs. (8) and (9) can be written asanddimensionless BCs are enumerated by

The domain of the system can be discretized with the following nodes:

whereas the term indicates step-size. Equations (14)–(16) with the central utilization difference can be approximated at the midpoint , likewise

Equations (18)–(23) can be further linearized with the help of a well-known scheme called Newton’s method by introducing the substitutions mentioned underneath:

Putting these expressions in (18)–(23) and dropping terms having higher powers in terms of to get the system of equation mentioned belowwhere

After the linearizing process, the tridiagonal block matrix mentioned below is achieved.wherewhere the elements defined in Eq. (33) are

Now we factorize A aswherewherein A is a tridiagonal block structure of order with each block having a size of , and is a unit block of order 6. Equation (33) can be tackled using the LU factorization approach to give numerical results of . Because the physical domain of the issue is unbounded, but the computing domain must be limited, we apply far-field boundary constraints so . In order to compute the numerical results, the step size and error tolerance level can be adjusted at and . A healthy comparison with previously published research might be used to assess the credibility of forthcoming outcomes of Hayat et al.[51] and Ishak et al.[52] in the case of by retaining (Table 3). Table 3 further compares the bvp4c MATLAB code and KBM solutions for various embedding parameter values for . A comparison analysis reveals a high level of agreement. In the case of makes the problem identical to that of Ahmad et al.[53].

Table 3

Nusselt number comparison by keeping .

Pr	Hayat et al.[51]	Ishak et al.[52]	bvp4c code	Present study
0.5	0.72325	0.72323	0.72320	0.72327
5	3.16243	3.16242	3.16235	3.16249
10	4.64681	4.64682	4.64679	4.64686
15	5.78369	5.78367	5.78368	5.78373

Results and discussions

This research’s primary goal is to explore nanofluids’ flow and heat transfer characteristics due to a stretching sheet. The values of surface drag and heat transfer in terms of –sodium alginate nanofluid are computed and displayed in Table 4. Magnification in and magnetic term depreciates the surface drag effect but elevates for the case of positive variation in Deborah number . Amplification in , , and Prandtl number amplifies heat transfer phenomenon while depreciates for the case of magnification in , and Eckert number . In the present section, the physical significance of sundry parameters such as Deborah number , magnetic parameter , porosity parameter , nanoparticle volume fraction , Prandtl number and Eckert number versus velocity, temperature and heat transfer rate is examined through Figs. 3, 4, 5, 6, 7 and 8. The field is noticeably enhanced for the rising and this outcome is in the augmentation of momentum boundary layer thickness. As depends upon the retardation time , higher offers rise to the fluid flow, because of which the velocity field is elevated. The effect of is monitored from Fig. 3, and it is discovered that the existence of diminishes the thickness of the boundary layer as well as the velocity field. This phenomenon happens because an increase in the permeability parameter increases the resistance to fluid motion, resulting in a reduction in velocity. This occurrence reflects in slower fluidity over the surface for Jeffery nanofluid in comparison with that of base fluid may be due to added particles making the nanofluid combo flow slower. Figure 4 displays effect on the thermal field . Here, the thermal field is a decreasing function of because causes a decrease in the molecular movements, which eventually decreases the fluid’s temperature profile. From Fig. 4, it is pretty clear that magnification in amplifies the temperature profile. More significant values of mean a rise in Lorentz force that is a resistive force. It generates a kind of friction on the flow field; more heat energy is caused due to this friction. Hence thermal profile enhanced. Figure 5 shows that the nanoparticle concentration increases the fluid temperature. Physically, increasing brings about an enhancement in the thermal conductivity phenomenon, which elevates the thermal profile. Figure 6 reflects the performance of Prandtl number on the temperature field for both viscous nanofluid () and Jeffrey nanofluid (), respectively. It can be noted from the figure that the fluid temperature reduces with escalating . This incidence is because of a decrement in amplifying thermal conductivities. Hence heat can diffuse away from the heated surface more quickly than larger values of . Consequently, in the case of the lesser , the heat transfer rate is reduced. It is observed that the introduction of reduces the thermal profile. The fluctuation of the local skin friction coefficient with β for various values of M is seen in Fig. 7. It has been found that when the value of β increases, the skin friction coefficient increases as well. This prevalence is because an increase causes more fluid particles to move into the boundary layer, particularly in the region of the sheet surface. As a result, the velocity boundary layer thickness decreases, resulting in greater skin friction coefficient values. On the other hand, the skin friction coefficient falls as M increases. The magnetic field has the effect of slowing down the velocity. As a result, the magnetic field may be employed to regulate flow properties.

Table 4

Numerical outcomes for and for diverse values of sundry dimensionless parameters.

ϕ	β	M	Pr	Ec	C_fRex1/2	Nu_XReX-1/2
0	0.5	0.5	6.5	0.3	0.209005	3.226895
0.01					0.198873	3.328853
0.02					0.199823	3.417300
0.01	0.0	0.2	6.5	0.3	0.339209	1.473574
	0.5				0.433021	1.596740
	1				0.507677	1.693113
0.01	0.2	0.0	6.5	0.3	0.412718	4.933286
		0.5			0.359445	4.363074
		1			0.308883	3.811261
0.01	0.5	0.5	5	0.3	0.197353	0.993617
			10		0.197353	3.672235
			15		0.197353	4.195064
0.01	2	0.5	6.5	0.5	0.575514	1.554988
				1	0.575514	1.118344
				1.5	0.575514	0.681700

Figure 3

Sway of on .

Figure 4

Sway of on .

Figure 5

Sway of on .

Figure 6

Sway of on .

Figure 7

Sway of on the skin friction coefficient.

Figure 8

Sway of on the skin friction coefficient.

It is noticed from Fig. 8 that the local skin friction coefficient decreases with increase the magnitude of volume fraction of nanoparticles ф. Figure 9 depicts the heat transfer rate with for diverse values of . It shows that the heat transport rate at the wall enhances with Deborah number and declines with Eckert number parameter . It can be noticed from Fig. 10 that heat transfer rate amplify with a rise in . Also, it is identified from Fig. 10 that the heat transfer rate reduces for rising values of . Physically, thermal conductivity is reduced for a more excellent value of Pr. Consequently, their heat conduction ability lessens. Heat transport rate abates near the wall.

Figure 9

Sway of on .

Figure 10

Sway of on .

Concluding remarks

The numerical solutions of MHD Jeffery nanofluid flow and heat transfer past over a stretching sheet with Joule heating and viscous dissipation effects are investigated. The calculated solutions are validated by comparing them to previously published results. The following are some of the findings that have been explored and summarised:

Flow field declines for rising the values of porosity parameter.

Thermal profile climbs with a boost in the concentration of nanoparticles and Lorentz force.

The interaction of Deborah number is proven to be counterproductive in embellishing temperature distribution though salubrious in accomplishing higher velocity inside flow field.

The thermal profile is abbreviating function of the Prandtl number.

Nanoparticles volume fraction supports the heat transfer rate but tends to reduce the skin frictions of the –sodium alginate nanofluid. Independent study on heat transfer rate shows that the parameters like Deborah number favour the Nusselt numbers while the Prandtl number and Eckert number tend to oppose it.

–Sodium alginate nanofluid possesses a small drag coefficient and heat transfer rate compared to the base fluid.

By and large, –Sodium alginate Jeffrey nanofluid combinations carry an upper hand in the main aspects of heat transfer efficiency compared to the sodium alginate base fluid.