Computer Modeling in Engineering & Sciences

DOI: 10.32604/cmes.2022.017539

ARTICLE

Flow and Melting Thermal Transfer Enhancement Analysis of Alumina, Titanium Oxide-Based Maxwell Nanofluid Flow Inside Double Rotating Disks with Finite-Element Simulation

Liangliang Chen1, Madeeha Tahir2,*, Sumeira Yasmin3, Taseer Muhammad4, Muhammad Imran5,* and Fenghua Liu1

1Chongqing Water Resources and Electric Engineering Collage, Chongqing, 402160, China
2Department of Mathematics, Government College Women University, Faisalabad, 38000, Pakistan
3Department of Mathematics, Government College University Faisalabad, Layyah Campus, Layyah, 31200, Pakistan
4Department of Mathematics, College of Sciences, King Khalid University, Abha, 61413, Saudi Arabia
5Department of Mathematics, Government College University, Faisalabad, 38000, Pakistan
*Corresponding Authors: Madeeha Tahir. Email: dr.madeeha.2017@gmail.com; Muhammad Imran Email: drmimranranchaudhry@gmail.com
Received: 19 May 2021; Accepted: 17 September 2021

1  Introduction

Advanced technologies realize the critical value of a unique type of energy transport fluids known as nanoliquid, because of the growing requirements of heat. The most significant aspect of the heat transfer system is the heat efficiency of the base fluid. Since non-metallic materials have poorer thermal efficiency than metallic substances. High thermal efficiency nanoparticles dispersed in regular fluids, which are typically made up of metals and oxides, greatly improve the heat proficiency of the host liquid. Therefore, metals are more useful to improving the thermal transfer rate. Choi and Eastman in (1995) coined the term “nanofluid” to describe the regular fluid that contained very tiny (1–100 nm) nanomaterials. Hayat et al. [1] reviewed the analysis of activation in Ree-Eyring nanofluid flow inside double disks. Qayyum et al. [2] scrutinized the entropy production inspirations on Williamson nanoliquid flow insides double rotating disks. Muhammad et al. [3] observed the slip impacts with activation energy across a three-dimensional sheet. Rafiq et al. [4] examined the numerical computations effects of nanofluid containing six different particles namely Ag, Cu, CuO, Fe3O4,TiO2 and Al2O3. Hassan et al. [5] debated Cu−Ag/water hybrid nanofluid flow through a cone. Sheikholeslami et al. [6] illustrate the impact of electro-hydrodynamic flow of Fe3O4−ethylene glycol nanofluid over an enclosure. The effects of magnetic fields on the hybrid nanomaterials are scrutinized by Shah et al. [7]. Arefmanesh et al. [8] discussed the influence of mixed convection in nanofluid through wavy wall cavity. Ashraf et al. [9] introduced the thermal transfer improvement in blood-base nanofluid across wavy tube. The theoretical analysis of three-dimensional Newtonian nanofluids with permeability porous medium is also examined by Ullah et al. [10]. Lu et al. [11] examined the squeezing flow of nanoliquid by adopting MDP.

The development of hybrid nanofluids, which are essentially an aqueous mixture of two or more forms of nanostructures in mixture or composite shape, is the next advancement in nanofluids technology. Hybrid nanofluids are being developed to solve the drawbacks of single suspension and to take benefit of nanoparticle synergy. Hybrid nanomaterials demonstrated that nanofluids have enhanced energy transfer and thermal conductivity, resulting in cost savings in commercial applications. There is little experimental, theoretical, and computational research on hybrid nanofluids. Li et al. [12] simulated the fractional study of hybrid nanoliquid flow across a spinning disk. Gul et al. [13] highlights the conical gap between cone and disk under hybrid nanofluid flow by adopting Homotopy analysis method. Waqas et al. [14] introduced the hybrid nanofluid flow over radiative disk. Shafee et al. [15] introduced the Entropy generation effects along NEPCM charging mechanics utilizing hybrid nanomaterials. Armaghani et al. [16] discussed the role of magnetized hybrid nanofluid flow through L-shape cavity. Shoaib et al. [17] explored the joule heating effects on magneto hybrid nanofluid over radiative rotating disk. The significance of magnetic dipole in hybrid nanofluid across stretchable surface is illustrated by Gul et al. [18]. Several evaluations have been discussed and presented in the literature, see for illustration [19,20].

Magnetohydrodynamics (MHD) is a class of physics that describes the magnetic properties of electrically interacting liquids. It usually influences thermal transport and manifests itself as Joule heating and Lorentz strength. MHD includes things like refrigerator cooling, saltwater, plasma, tumor therapy, X-ray radiation, and electrolytes. Many excellent works describing various aspects of MHD have been released. Reddy et al. [21] explored the MHD Cu−Ag/water mono nanoliquid flow. The significant features of flow and thermal transfer in Cu/water nanofluid inside double plates are scrutinized by Wakif et al. [22]. The effects of entropy generation on magnetized hybrid nanoliquid flowing configured by a sheet are investigated by Mumraiz et al. [23].

The concept of porous media is used in various fields of technology and applied science, including mechanics, filtration, petroleum engineering, construction engineering, hydrogeology, and biophysics. MHD CuO−water nanofluid flow within the porous cavity was explored by Sheikholeslami [24]. Waqas et al. [25] analyzed the viscous nanofluid flowing across a porous stretchable revolving disk embedded in porous medium. The features of second order slip on nanofluid contain nano-sized nanomaterials past curved channels under porous permeability are explored by Riaz et al. [26].

The thermal radiation behavior in thermal transfer has applications in a variety of thermal engineering fields, including hybrid power solar systems, rocket propulsion, nuclear reactors, rockets, aircraft, and communication technology. It is worth noting that linear radiation is adequate to produce thermal equipment and is widely used in several technological processes. Hayat et al. [27] analyzed the impacts of thermal radiation on Carbon/water nanofluid. The behavior of thermal radiation in hydrothermal nature of nanofluid is scrutinized by Sheikholeslami et al. [28]. Li et al. [29] numerically investigate the radiative modified second grade nanofluid aspects. Naqvi et al. [30] studied the heat radiation effects on Casson nanofluid in the occurrence of magnetic field.

Melting heat has piqued the interest of scientists and researchers owing to its various uses in novel industrial processes. In recent times, scientists have concentrated their efforts on creating more long-term, reliable, and affordable energy storing technologies. Unintended thermal efficiency, solar energy and power are all examples of such technologies. Roberts [31] the first to depict the nature of melting thermal transfer in a hot air stream of ice surface. The behavior of melting impact with bioconvection in generalized second grade nanofluid configured by stretching surface is illustrated by Waqas et al. [32]. Hayat et al. [33] introduced the melting thermal transport in stagnation point flow of nanofluid with zero mass flux condition across a stretching surface. Ullah et al. [34] discussed the numerical solution of SWCNTs−MWCNTs/Engineoil hybrid nanofluid flow with melting mechanism. The behavior of melting heat on magneto nanofluid flow past stretching sheet is analyzed by Sharma et al. [35].

As reviewed above, the non-Newtonian nature of nanofluid has vital applications. The key determination of this article is to scrutinize the MHD flow and melting heat transfer for a nanofluid model through a two parallel rotating disks. Novelty of the communication is firstly to scrutinize the axisymmetric flow of Maxwell-water based nanofluid with thermal radiation. Secondly the numerical solution of current problem is solved numerically by applied Finite element method. The numerical solutions are obtained and the role of physical factors on the subjective dimensionless profiles is scrutinized and reflected by graphs.

2  Mathematical Modeling

2.1 Mathematical and Physical Flow Description

Considerable the effect of steady incompressible, axisymmetric electromagnetic flow with heat transformation by using different particles based on Maxwell nanofluid between both revolving stretching disks at the distance h from upper to lower disks. The upper and lower disks have the distinct angular velocities denoted by (ω1&ω2) and both disks have distinct stretching rates (b1&b2) along the radial direction and axial direction correspondingly. The magnetic field strength B0 alongz−axis is presented and the cylindrical coordinate system (r−,θ−,z−) is considered in Fig. 1.

Figure 1: Physical configuration of flow issue

2.2 Dimensional Governing Boundary Flow Expressions

The governing flow equations are given bellow [36,37]:

∂u∂r+ur+∂w∂z=0,(1)

u∂u∂r−v2r+w∂u∂z+λ1(v2∂2u∂z2+u2∂2u∂r2+2uv∂2u∂r∂z)=−1ρnf∂p∂r+μnfρnf[∂2u∂r2+∂2u∂z2+1r∂u∂r−ur2]−νnfk∗u+σnfρnfβ02(−u−λ1v∂u∂z),(2)

u∂v∂r+uvr+w∂v∂z+λ1(v2∂2u∂z2+u2∂2u∂r2+2uv∂2u∂r∂z)=−1ρnf∂p∂r+μnfρnf[−vr2+∂2v∂r2+1r∂v∂r+∂2v∂z2]−νnfk∗v+σnfρnfβ02(−v−λ1u∂v∂z),(3)

u∂w∂r+w∂w∂z=−1ρnf∂p∂z+μnfρnf[1r∂w∂r+∂2w∂r2+∂2w∂z2],(4)

(u∂T∂r+w∂T∂z)=knf(ρcp)nf[∂2T∂r2+1r∂T∂r+∂2T∂z2]−1(ρcp)nf∂qr∂z,(5)

With,

u=rb1,v=rω1,w=0,T=T0,atz=0,u=rb2,v=rω2,w=0,T=T1,atz=h.(6a)

The melting condition is addressed as

knf(∂T∂z)z=0=ρnf[L0+(cp)s(T0−T1)]w.(6b)

Here the electrical conductivity represented by σnf, the thermal diffusivity simplified as αnf, the density for nanofluid is ρnf, T0 be the melting temperature, the nanofluid heat capacitance symbolized as (ρcp)nf, the dynamic viscosity denoted by μnf, the thermal conductivity for nanofluid demonstrated by knf the kinematic viscosity of base nanofluid characterized by νf.

where nf elucidated the thermophysical properties of nanofluid, solid particles symbolized as s, base fluid can be written as f and the solid volume fraction denoted by ϕ for fluid nanoparticles. “For the said flow the equation of continuity is satisfied for detail see Appendix”.

Now by employing the Rosseland evaluation of radiative heat flux qr is considered as

qr=−4σ∗3K∗∂T4∂z(7)

Here consider that the heat transfers in the flow field are such that the expression T4 may be designed as thermal dependent linearly functions. This is capable by expending T4 in Taylor sequence about the ambient heat T∞ as bellow described.

T4=T∞4+4T∞3(T−T∞)+6T∞2(T−T∞)2+⋅⋅⋅⋅⋅⋅(8)

In above Eq. (8) eliminate the terms of multi order and outer the first degree in (T−T∞) then we get:

T4=~4T∞3T−3T∞4,(9)

Now putting Eq. (9) into Eq. (7) then we have,

qr=−16T∞3σ∗3K∗∂T∂z.(10)

2.3 Similarity Transformations

The stream function ψ can be demarcated as follows [38]:

u=∂ψ∂y,v=−∂ψ∂x.(11)

The transformations are as follows:

u=rω1f′(ζ),w=−2ω1hf(ζ),v=rω1g(ζ),ζ=zh,θ(ζ)=T−T1T0−T1,p=ρfω1νf(p(ζ)+12r2h2+ε).(12)

2.4 Dimension-Less Equations

After applying the similarity variables, we obtain the non-dimensional expressions are as follows:

f′′′−A12Re(−2ff′′+(f′)2−g2)−ReβA1(f2f′′′−2ff′f′′)−A2KRef′−A1A2M(f′+βff′′)−A1A2ε=0,(13)

g′′+A1Re(2fg′−2f′g)−A2KReg−A1A2M(g+βfg′)=0,(14)

P′=4A2Reff′+2A2A1f′′,(15)

(1+Rd)θ′′+2RePrA3A4fθ′=0,(16)

With

ζ=0,Prρnfρff(0)+Maknfkfθ′(0),f′=B1,g=1,θ′=1,ζ=1,f′=B2,g=B3,θ′=0.(17)

2.5 Non-Dimensional Parameters

The dimension-less sundry parameters are given bellow:

[Re=ω1h2νf,Rd=16σ∗T∞33kfK∗,M=σβ022ω1ρf,Pr=νfαf,β=λ1a,B3=ω1ω2,(B1=b1ω1,B2=b2ω1),Ma=(cp)f(T1−T1)L∗+(cp)S(T0−T1)andA1=(1−ϕ)2.5(1−ϕ+ϕ(ρsρf)),A2=(1−ϕ)+ϕ(ρsρf),A3=(1−ϕ)+ϕ(ρCp)s(ρCp)f,A4=knfkf.](18)

Here the Eq. (13) differentiated w.r.t. to ζ and eliminate ε then we get shorter expression is as follows:

fiv+A1Re(ff′′′+gg′)−A2KRef′′+ReβA1(f2f′′′+2f′2f′′+2ff′′2)−A1A2M(f′′+βf′f′′+βff′′′)=0.(19)

Through Eqs. (13) and (17) the parameter ε which shows the pressure can be expressed as:

ε=A2A1f′′′(0)−A2KRef′′(0)−A22Re[(f′(0))2−2f(0)f′′(0)−g2(0)]−ReβA2[(f(0))2f′′′(0)−2f(0)f′(0)f′′(0)]−M(f′(0)+βf′(0)f′′(0)).(20)

The above Eq. (15) integrated with respect to ζ and determined the pressure equation, by taking the limits starting (0toζ).

p+2A2Ref2+2A2A1(f′−f′(0))=0,(21)

2.6 Physical Quantities

The Physical industrial inters in which the coefficients of skin frictions (C1&C2) and the coefficients of Nusselt numbers (Nux1&Nux2) for both discs are described.

C1=τrzρf(rω1)2,C2=τθzρf(rω2)2,(22)

Nux1=hqwkf(T0−T1),Nux2=hqwkf(T0−T1).(23)

Here we take stress (τzrandτzθ) elucidated the shear stresses for the lower disc along radial as well as tangential direction and the heat flux qw correspondingly.

The skin frictions coefficients in non-dimensional forms are

C1=τw|z=0ρf(rω1)2=1Rer(1−ϕ)2.5[(f′′(0))2+(g′(0))2]12,(24)

C2=τw|z=hρf(rω1)2=1Rer(1−ϕ)2.5[(f′′(1))2+(g′(1))2]12,(25)

And the coefficients of Nusselt numbers in non-dimensional forms are

Nu1=−(A4+Rd)θ′(0),Nu2=−(A4+Rd)θ′(1),(26)

Here the total shear stress τw addressed as:

τw=τzr2+τzθ2(27)

3  Numerical Process: Finite Element Method (FEM)

Here, the nonlinear dimensionless ODE's (14)–(16) and (19) with specific boundary conditions (17) are tackled by utilizing the finite element method (FEM) or finite element analysis (FEA). These ODEs are highly nonlinear, so analytical techniques are not useful to solve these expressions. Therefore, the (FEM) is used to compute the solution of highly nonlinear system. The FEM is more effective than other technique.

3.1 Finite Element Technique (FET)

We get the exact solution through the finite element analysis the coupled nonlinearly ODE's structure (14)–(16) and (19) with specified boundary conditions (17) firstly we let

f′=s,f′′=s1,f′′′=s2,fiv=s3,(28)

s3+A1Re(fs2+gg′)−A2KRes2+ReβA1(f2s2+2s2s1+2fs12)−A1A2M(s1+βss1+βfs2)=0,(29)

g′′+A1Re(2fg′−2f′g)−A2KReg−A1A2M(g+βfg′)=0,(30)

P′−4A2Refs−2A2A1s1=0,(31)

(1+Rd)θ′′+2RePrA3A4fθ′=0.(32)

With

ζ=0,Prρnfρff(0)+Maknfkfθ′(0),s=B1,g=1,θ′=1,ζ=1,s=B2,g=B3,θ′=0.(33)

3.2 Variational Formulation

The variational method associated with the expressions (28)–(32) by typical linear component (ζe,ζe+1) is explain as

∫ζeζe+1ω1(∂f∂ζ−s)dζ=0,(34)

∫ζeζe+1ω2[s′′−A12Re[s2−2fs′−g2]−A2KRes′+−ReβA1[f2s′′−2fss′]−A1A2M(s+βfs′)]dζ=0,(35)

∫ζeζe+1ω3[g′′+A1Re[2fg′−2hg]−A2KReg−A1A2M(g+βgg′)]dζ=0,(36)

∫ζeζe+1ω4[(1+Rd)θ′′+2A3A4PrRefθ′]dζ=0,(37)

Here (ω1,ω2,ω3andω4) described the arbitrary test functions and may be observed as the values in (f,s,gandθ) etc.

3.3 Finite-Element Description

Substitute finite element estimation of the form declared below; the finite element method may be attained from the abovementioned expressions:

f=∑j=13fjψj,s=∑j=13sjψj,g=∑j=13gjψj,θ=∑j=13θjψj,(38)

By

ω1=ω2=ω3=ω4=ψi(i=1,2,3).

For typical components(ζe,ζe+1), the above-represented form function ψi can be defined as

ψ1e=(ζe+1+ζe−2ζ)(ζe+1−ζ)(ζe+1−ζ)2,ψ2e=4(ζ−ζe)(ζe+1−ζ)(ζe+1−ζ)2,ψ3e=(ζe+1+ζe−2ζ)(ζ−ζe)(ζe+1−ζ)2,ζe≤ζ≤ζe+1,(39)

The structure of finite-element method for nonlinear governing equations is mentioned as:

[[K11][K12][K13][K14][K21][K22][K23][K24][K31][K32][K33][K34][K41][K42][K43][K44]][{f}{s}{g}{θ}]=[{r1}{r2}{r3}{r4}],

Here [Knm]and[rn](m,n=1,2,3,4) are established as

Kij11=∫ζeζe+1ψi∂ψj∂ζdζ,Kij12=−∫ζeζe+1ψiψjdζ,Kij13=0,Kij14=0,Kij21=0,Kij23=A12Re∫ζeζe+1ψi(ψj)2dζ,Kij24=0,

Kij22=−∫ζeζe+1∂ψi∂ζ∂ψj∂ζdζ−βRe[∫ζeζe+1(ψi)2∂ψi∂ζ∂ψj∂ζdζ+∫ζeζe+1ψiψj∂ψj∂ζdζ]+A12Ref¯1∫ζeζe+1ψiψ1∂ψj∂ζdζ+A12Ref¯2∫ζeζe+1ψiψ2∂ψj∂ζdζ−KReA2∫ζeζe+1ψiψjdζ−A1A2(M)∫ζeζe+1ψiψjdζ−A12Re∫ζeζe+1ψi(ψj)2dζ+A1A2(M)βf¯1∫ζeζe+1ψiψ1∂ψj∂ζdζ+A1A2Mβf¯2∫ζeζe+1ψiψ2∂ψj∂ζdζ,

Kij31=0,Kij32=−A1Reh¯1∫ζeζe+1ψiψ1dζ+A1Reh¯2∫ζeζe+1ψiψ2dζ,Kij34=0,

Kij33=−∫ζeζe+1∂ψi∂ζ∂ψj∂ζdζ+A1Ref¯1∫ζeζe+1ψiψ1∂ψj∂ζdζ+A1Ref¯2∫ζeζe+1ψiψ2∂ψj∂ζdζ−A1A2

Kij41=0,Kij42=0,Kij43=(1+Rd)∫ζeζe+1∂ψi∂ζ∂ψj∂ζdζ+2PrReA3A4f¯1∫ζeζe+1ψiψ1∂ψj∂ζdζ+34A3A4f¯2∫ζeζe+1ψiψ2∂ψj∂ζdζ,Kij44=0,

ri2=0,ri2=−(ψidψi∂ζ)ζeζe+1,ri3=−(ψidψi∂ζ)ζeζe+1,ri4=−(ψidψi∂ζ)ζeζe+1.

3.4 Validation of Results

We used finite element method for tackling the nonlinear set of ODEs. Because fewer nodes are required in finite element, less memory is required to perform the full program, resulting in improved accuracy and, as a result, a shorter calculation time than FEM. Physical features like density, heat capacity, thermal conductivity and electrical conductivity of some particles and base liquid are given in Table 1. Characteristics for some physical properties of nanofluids are given in Table 2. For validation of code we compare current outcomes with the numerical outcomes found by Lance et al. [40], Turkyilmazoglu [41] and Kumar et al. [42] for the rotation parameter B3 by fixed B1=B2=M=ϕ=0 and Re=1.0 as summarized in Table 3. We observed good agreement between our results and published literature. Here minimal percentage errors are also display in Table 3.

4  Results and Discussion

Fig. 2 indicates the radial velocity profile f′ for varying variations in magnetic parameter M for Al2O3/water and TiO2/water nanofluid flow between double radiative nanofluids. It is clearly illustrated that radial velocity field f′ is reduced via magnetic parameter M. Magnetic parameter develops a resistance force known Lorentz force that is a delaying force in flow of fluid. The radial velocity component f′ of Maxwell nanoliquid with change in Reynolds number Re is depicted in Fig. 3. From this picture it is learnt that the Reynolds number Re depresses the radial component of velocity f′. The physical reason for this is because when the lower disk rotates; the inertial impacts grow, slowing down the flow rate. Fig. 4 reveals the vital role of the radial velocity of fluid f′ for the various variations of the nanoparticles volume fraction ϕ. The increment in the velocity component f′ is analyzed for higher nanoparticles volume friction ϕ. The converse relationship between the volume fraction and the dynamic viscosity of nanofluid causes the increase in velocity as a physical phenomenon. As a result, the viscosity of ordinary fluid is reduced in response to an increase in nanoparticle volume fraction, and fluid flow is increased [43]. The variation in radial velocity f′ of water-based Maxwell single nanofluid with stretching parameter B2 for upper disk is clarified in Fig. 5. It is mentioned that radial velocity f′~of nanoliquid is depresses with stretching parameter for upper disk B2. Fig. 6 illustrate the effect of fluid parameter β against radial flow f′ of Al2O3/water and 44 nanofluid for distinguish values of fluid parameter β. Here we observed that larger fluid parameter β, decays the velocity of fluid. The radial velocity f′ for Al2O3/water and TiO2/water nanofluid against melting parameter Ma is elucidated in Fig. 7. Improve in melting parameter Ma reduces the radial velocity field f′ of single nanoliquid inside parallel two disks. The significant decrement in radial velocity profile f′ depicts in Fig. 8, due to growing variations of porosity parameter K. Fig. 9 is plotted to show the behavior of stretching parameter for upper disk B2 against tangential velocity g of Al2O3/water and TiO2/water nanofluid flow. It is illustrated that with an increment in variation of stretching parameter B2, the tangential flow g of nanofluid increases. Fig. 10 signifies the significant features of melting parameter Ma vs. tangential velocity g of Al2O3/water and TiO2/water nanofluid flow over two parallel disks. It is examined that tangential velocity of fluid g is decreases via larger variations in the melting parameter Ma.

Figure 2: Importance of M via f′

Figure 3: Importance of Re via f′

Figure 4: Importance of ϕ via f′

Figure 5: Importance of B2 via f′

Figure 6: Importance of β via f′

Figure 7: Importance of Ma via f′

Figure 8: Importance of K via f′

Figure 9: Importance of B2 via g

Figure 10: Importance of Ma via g

The effect of thermal radiation parameter Rd against thermal field of nanomaterials θ is disclosed in Fig. 11. It is analyzed that heat transfer θ is increased with greater thermal radiation parameter Rd. Physically the absorption coefficient reduces for greater thermal radiation parameter and therefore the thermal field enhances. Fig. 12 examines the estimation in energy field θ with different amounts of melting parameter Ma for Al2O3/water and TiO2/water nanofluid flow over two parallel disks. We noted that temperature distribution θ is decays via larger melting parameter Ma. Fig. 13 is designed to scrutinize the behavior of Prandtl number Pr against heat profile θ. It is mentioned that by growing the Prandtl number Pr, the solutal field of species is declined. Since the Prandtl number indicates how quickly thermal diffusion occurs in compared to momentum diffusion, increasing the Prandtl number causes the thermal diffusivity to decrease and the temperature to drop. Fig. 14 shows that temperature distribution θ is decreased with solid volume fraction ϕ of Al2O3/water and TiO2/water nanofluid flow.

Figure 11: Importance of Rd via θ

Figure 12: Importance of Ma via θ

Figure 13: Importance of Pr via θ

Figure 14: Importance of ϕ via θ

5  Conclusions

The current study novelty is recent progress in melting heat transfer transportation of Maxwell Al2O3−TiO2/water nanofluid flow through double parallel rotating disks. The conclusive remarks of the present work are as follows:

•   An increase in radial velocity of single Maxwell water-based nanofluid is analyzed vs. an improvement in nanoparticles volume friction.

•   Tlhe radial velocity of nanofluid is reducing function of larger magnitudes of the melting parameter.

•   The fluid flow declines with fluid parameter.

•   The significant value of the tangential velocity is noticed via stretching parameter for upper disk.

•   Temperature distribution raises for thermal radiation parameter.

•   The thermal field of the fluid reduces as Prandtl number [44] increases.

•   Larger solid volume friction decreases the temperature distribution.

Acknowledgement: This work was sponsored in part by National Natural Science Foundation of China (No. 51869031); and Science and Technology Research Program of Chongqing Municipal Education Commission (Grant No. KJQN201903801) and Huzhou Key Laboratory of Green Building Technology.

Funding Statement: This work is financially supported by the Government College University, Faisalabad and Higher Education Commission, Pakistan.

Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the current analysis.

Appendix

∂u∂r+ur+∂w∂z=0,(41)

We have

∂u∂r=rω1f′,∂u∂z=−2ω1f′,andur=ω1f′

Putting all required values in Eq. (A1) we, have

=rω1f′+ω1f′−2ω1f′=0

Hence equation of continuity trivially satisfied.

This work is licensed under a Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.