
The current research aims to investigate thermodynamic responses to thermal media based on a modified mathematical model in the field of thermoelasticity. In this context, it was considered to present a new model with a fractional time derivative that includes Caputo-Fabrizio and Atangana-Baleanu fractional differential operators within the framework of the two-phase delay model. The proposed mathematical model is employed to examine the problem of an unbounded material with a spherical hole experiencing a reduced moving heat flow on its inner surface. The problem is solved analytically within the modified space utilizing the Laplace transform as the solution mechanism. An arithmetic inversion of the Laplace transform was performed and presented visually and tabularly for the studied distributions. In the tables, specific comparisons are introduced to evaluate the influences of different fractional operators and thermal properties on the response of all the fields examined.
Citation: Osama Moaaz, Ahmed E. Abouelregal. Multi-fractional-differential operators for a thermo-elastic magnetic response in an unbounded solid with a spherical hole via the DPL model[J]. AIMS Mathematics, 2023, 8(3): 5588-5615. doi: 10.3934/math.2023282
[1] | Mohammed Alrehili . Managing heat transfer effectiveness in a Darcy medium with a vertically non-linear stretching surface through the flow of an electrically conductive non-Newtonian nanofluid. AIMS Mathematics, 2024, 9(4): 9195-9210. doi: 10.3934/math.2024448 |
[2] | Nadeem Abbas, Wasfi Shatanawi, Taqi A. M. Shatnawi . Innovation of prescribe conditions for radiative Casson micropolar hybrid nanofluid flow with inclined MHD over a stretching sheet/cylinder. AIMS Mathematics, 2025, 10(2): 3561-3580. doi: 10.3934/math.2025164 |
[3] | Asifa, Poom Kumam, Talha Anwar, Zahir Shah, Wiboonsak Watthayu . Analysis and modeling of fractional electro-osmotic ramped flow of chemically reactive and heat absorptive/generative Walters'B fluid with ramped heat and mass transfer rates. AIMS Mathematics, 2021, 6(6): 5942-5976. doi: 10.3934/math.2021352 |
[4] | W. Abbas, Ahmed M. Megahed, Osama M. Morsy, M. A. Ibrahim, Ahmed A. M. Said . Dissipative Williamson fluid flow with double diffusive Cattaneo-Christov model due to a slippery stretching sheet embedded in a porous medium. AIMS Mathematics, 2022, 7(12): 20781-20796. doi: 10.3934/math.20221139 |
[5] | Khalil Ur Rehman, Wasfi Shatanawi, Zeeshan Asghar, Haitham M. S. Bahaidarah . Neural networking analysis for MHD mixed convection Casson flow past a multiple surfaces: A numerical solution. AIMS Mathematics, 2023, 8(7): 15805-15823. doi: 10.3934/math.2023807 |
[6] | S. R. Mishra, Subhajit Panda, Mansoor Alshehri, Nehad Ali Shah, Jae Dong Chung . Sensitivity analysis on optimizing heat transfer rate in hybrid nanofluid flow over a permeable surface for the power law heat flux model: Response surface methodology with ANOVA test. AIMS Mathematics, 2024, 9(5): 12700-12725. doi: 10.3934/math.2024621 |
[7] | Khalil Ur Rehman, Nosheen Fatima, Wasfi Shatanawi, Nabeela Kousar . Mathematical solutions for coupled nonlinear equations based on bioconvection in MHD Casson nanofluid flow. AIMS Mathematics, 2025, 10(1): 598-633. doi: 10.3934/math.2025027 |
[8] | Fahad Alsharari, Mohamed M. Mousa . New application of MOL-PACT for simulating buoyancy convection of a copper-water nanofluid in a square enclosure containing an insulated obstacle. AIMS Mathematics, 2022, 7(11): 20292-20312. doi: 10.3934/math.20221111 |
[9] | C. S. K. Raju, S.V. Siva Rama Raju, S. Mamatha Upadhya, N. Ameer Ahammad, Nehad Ali Shah, Thongchai Botmart . A numerical study of swirling axisymmetric flow characteristics in a cylinder with suspended PEG based magnetite and oxides nanoparticles. AIMS Mathematics, 2023, 8(2): 4575-4595. doi: 10.3934/math.2023226 |
[10] | Subhan Ullah, Hassan Ali Ghazwani, Dolat Khan, Zareen A. Khan . Heat transfer augmentation of Jeffery–Hamel hybrid nanofluid in a stretching convergent/divergent channel through porous medium. AIMS Mathematics, 2025, 10(1): 388-402. doi: 10.3934/math.2025018 |
The current research aims to investigate thermodynamic responses to thermal media based on a modified mathematical model in the field of thermoelasticity. In this context, it was considered to present a new model with a fractional time derivative that includes Caputo-Fabrizio and Atangana-Baleanu fractional differential operators within the framework of the two-phase delay model. The proposed mathematical model is employed to examine the problem of an unbounded material with a spherical hole experiencing a reduced moving heat flow on its inner surface. The problem is solved analytically within the modified space utilizing the Laplace transform as the solution mechanism. An arithmetic inversion of the Laplace transform was performed and presented visually and tabularly for the studied distributions. In the tables, specific comparisons are introduced to evaluate the influences of different fractional operators and thermal properties on the response of all the fields examined.
a1 | Amplitude of the upper wall, m; |
Re | Reynolds number; |
b1 | Amplitude of the lower wall, m; |
Pr | Prandtl number; |
a | Unit-less amplitude of the upper wall; |
Ec | Eckert number; |
b | Unit-less amplitude of the lower wall; |
Br | Brinkman number; |
d1+d2 | Channel width, m; |
Sij | Cauchy stress tensor; |
d | Unit-less half channel width; |
RPh | Reiner-Philippoff; |
(⌢X,⌢Y) | Cartesian coordinates, m; |
⌢t | Time, s; |
u | Velocity vector; |
t | Unit-less time; |
(⌢U,⌢V) | Velocity components, m/s; |
T | Temperature of fluid, K; |
(x,y) | Unit-less Cartesian coordinates; |
Nb | Unit-less Brownian parameter; |
(u,v) | Unit-less velocity components. |
k1 | Dimensionless non-uniformity parameter; |
δ | Wave number, S/m; |
⌢H1 | Geometry of the upper wall, m; |
ζ | Phase difference; |
⌢H2 | Geometry of the lower wall, m; |
λ | Wavelength, m; |
c | Peristaltic wave speed, m/s; |
ψ | Stream function; |
h1 | Unit-less geometry of the upper wall; |
Φ | Viscous dissipation; |
h2 | Unit-less geometry of the lower wall; |
λ1 | RPh fluid parameter; |
ˉC | Specific heat, j/kg.K; |
ρ | Density, kg/m3; |
C | Concentration, mol/m3; |
θ | Unit-less temperature; |
T0 | Reference temperature at upper wall, K; |
μ | Dynamic viscosity, Ns/m2; |
T1 | Reference temperature at lower wall, K; |
ν | Kinematics viscosity; |
C0 | Reference concentration at upper wall; |
⌢τ | Shear stress, N/m2; |
C1 | Reference concentration at lower wall; |
τ | Dimensionless shear stress; |
⌢P | Fluid's pressure, N/m2; |
⌢τ0 | Reference shear stress, N/m2; |
p | Dimensionless fluid's pressure; |
ϕ | Dimensionless concentration; |
⌢Q | Instantaneous volumetric flow rate; |
α | non-uniformity angle, radian; |
ˉQ | Time average flow rate; |
γ | Shear stress parameter; |
η | Dimensionless mean flow rate; |
μ0 | zero shear viscosity, Ns/m2; |
DB | Brownian diffusion coefficient. |
DT | Thermophoretic diffusion coefficient; |
f | Base fluid; |
K | Thermal conductivity, W/mK; |
p | Nanoparticles; |
Nt | Unit-less thermophoretic parameter. |
Demand of getting insight of the peristaltic flow is increasing day by day due to its applications in many fields of science and engineering especially in the fields of biology and medical engineering. Spectrum of research on the topic of peristaltic flow got widen by involving non-Newtonian fluids since most of the fluid involved in different biological systems and other field of engineering are non-Newtonian in nature. Peristalsis is a phenomenon in which the sinusoidal waves arise due to the continuous relaxation and diminution of smooth muscular tissues. Some of the natural and engineering systems involving the application of peristalsis motion are locomotion of earthworm, artificial heart pump, transportation of urine starting from the kidneys to the urinary bladder, vasomotion of tiny blood arteries, motion of the food along the digestive tube and blood transport vessels. The Rotho-peristaltic pump is also one of the applications of peristalsis in chemical industry.
Before Thurston [1], blood was treated as a Newtonian fluid. However, viscoelasticity is a basic rheological property of blood which makes the human blood non-Newtonian. Recently interest in problems of non-Newtonian fluid has grown and many mathematical models for describing the rheological behavior of blood have been extensively developed. Many mathematical structures for elaborating the rheological characteristics of blood have been briefly investigated in references [2,3,4,5]. Latham [6] and Shapiro et al. [7] primary reported the fluid dynamical simulation of peristalsis for 2D Newtonian fluid. Tripathi [8] addressed the peristaltic transport of high viscosity across a finite length of channel under the influence of porosity and heat transmission. Hayat et al. [9] explored the significance of hall current on peristaltic flow of Maxwell's fluid passing via porous medium. The latest research on the peristaltic subject can also be seen through references [10,11,12,13,14].The Reiner-Philippoff fluid model is one of the most challenging non-Newtonian fluid models due to implicit relationship of shear stress and strain. Due to this mathematical difficulty, limited literature is available on Reiner-Philippoff fluid model. In 1994, Na [15] was the first to examine the boundary layer flow of Reiner-Philippoff fluid due to free stream velocity. The boundary layer flow of Reiner-Philippoff fluid over a stretching sheet was studied by Yam et al. [16] Ahmad [17] investigated the boundary layer flow of Reiner-Philippoff based nanofluid past a non-linear stretching sheet. Impact of Brownian and thermophoresis diffusions on thin film flow of Reiner-Philippoff fluid was studied by Asad Ullah [18]. Sajid et al. [19] investigates the impact of triple diffusive convective flowing on Reiner-Philippoff (RPh) fluid past a stretching surface with Avramenko-Blinov-Shevchuk boundary restrictions. Also Reiner-Philippoff (RPh) fluid model is discussed by considering blood as a base fluid by Sajid et al. [20]. Tahir and Ahmad [21] initiated the study of peristaltic transportation of Reiner-Philippoff fluid.
Heat transfer plays a significant role in equipment efficiency and material selection. The importance of peristaltic fluid motion with heat transfer is noticeable in biomedical sciences, metabolic heat generation, blood transport, perfusion hypothermia, thermoregulation and in bio-heat exchange processes. The working capability of every biological and mechanical process decays time to time due to the production of heat by such processes. To improve the efficiency of such mechanism it is compulsory to eliminate the heat regularly from the system. Nanofluids are innovative compounds that can control the heat transfer. Nanofluids are made by nanoprecipitation method in ordinary fluids. The ordinary fluids like water, polymer solutions, engine oil and some lubricants have less thermal conductivity values as compared to solids. Including a specific quantity of solid particles to these fluids, can easily enhance the heat conductance of these fluids. These particles are usually of sizes ranging 1-100 nanometer and known as nanoparticles. The consequential mixture of fluid and nanoparticles is called nanofluid. This idea of constructing nanofluids was primarily introduced by Choi [22,23] in 1995. Xuan [24] talked about the improvements of heat transfers for nanofluids. Buongiorno [25] proposed a model that incorporates thermophoresis and Brownian motion phenomena to evaluate the transport of nanofluids. Nanofluids are also helpful and applicable in the understanding of a variety of phenomenon such as the minimization/enhancement of the thermal magnitude of the system, microchannel clogging, miniaturization of the system and minimal clogging. Recently Rasool et al. [26] investigates and provides a numerical analysis of electro-magneto-hydrodynamic (EMHD) nanofluid flows past a Riga pattern embedded horizontally in a Darcy-Forchheimer porous medium. The more literature regarding heat transfer through nanofluids can be seen through references [27,28,29,30,31,32,33]. Many academics are interested in the combined impacts of heat and mass transfer because of its vast spectrum of uses in industry, science, and engineering. These procedures play a critical role in enhancing the product quality. They're used in a variety of pharmaceutical and chemical operations, as well as magnetic drug targeting and electronic devices cooling. Thermal insulation, thermal energy storage, fossil fuel extraction, the paper industry, food processing, oxygenation, porous solid drying, hemodialysis, subterranean energy transmission and geothermal energy recovery are some additional physical examples. Hina [34] in 2012 explored how the compliant walls affected peristaltic flow, heat transfer, mass transport and chemical reactions. Peristaltic movement study of Eyring-Prandtl fluid in a tube together with heat and mass transfer analysis was investigated by Iftikhar et al. [35]. Researcher's contributions to peristaltic flows are fast expanding due to the relevance of heat and mass transport in several fields of engineering and science (see citations [36,37,38,39]).
As discussed earlier that negligible research was available on non-Newtonian Reiner-Philippoff fluid model. And particularly no literature was available on peristaltic motion of Reiner-Philippoff fluid. In previous work [40] heat transfer analysis of Renner-Philippoff nanofluid was discussed by using conventional nanofluid model. Considering the importance of non-Newtonian fluid in peristaltic flow and role of nanoparticles as an efficient agent of heat transfer control, the aim of this article is to explore heat and mass transfer analysis of non-Newtonian Reiner-Philippoff fluid in peristaltic motion in the presence of nanoparticles by using Buongiorno's nanofluid model. The impact of thermophoretic diffusion parameter, Brownian diffusion parameter and Reiner-Philippoff fluid parameter on different quantities of interest is discussed in detail. I hope this study will be step forward in the field of peristaltic flow, heat and mass transfer of non-Newtonian fluids and also provides new avenues of research for scholars to the way of RPh fluid academics.
Consider Reiner-Philippoff (RPh) fluid flowing through a non-uniform (tapered) two-dimensional asymmetric channel to be incompressible. The fluid movement is caused by the sinusoidal waves having wavelength λ and constant speed c along the channel walls. The walls of the non-uniform channel are separated by width d1(⌢X)+d2(⌢X). The coordinate system for the given problem is Cartesian system (⌢X,⌢Y) introduced in such a way that ⌢Y−axis is chosen alongside the width and the ⌢X−axis is along the length of the channel. Geometrical deformation of the walls for discussed problem can be put mathematically as [41]:
⌢H1(⌢X,⌢t)=d1(⌢X,⌢t)+a1cos(2πλ(⌢X−c⌢t)),upperwall | (1) |
⌢H2(⌢X,⌢t)=−d2(⌢X,⌢t)−b1cos(2πλ(⌢X−c⌢t)+ζ),lowerwall, | (2) |
where, di(⌢X,⌢t)=di+(⌢X−c⌢t)tan(α),(i = 1, 2), is the half width of channel, α is the non-uniformity angle of the channel, ⌢H1(⌢X,⌢t) and ⌢H2(⌢X,⌢t) symbolized to show the geometry of the said problem for topper and lower walls respectively corresponding to the coordinate ⌢Y>0 and ⌢Y<0. a1 and b1 are the amplitudes of the topper and lower waves respectively, ⌢t is for the time and ζ(0≤ζ≤π) represents phase difference between the waves, for ζ=0 the channel is reduced to symmetric in which waves are totally out of phase and for ζ=π waves are totally in phase (see Figure 1). Superscript ⌢. is used to indicate the dimensional quantities. Further, d1,a1,d2 and b1 satisfies the following [42]:
a21+b21+2a21b21cosζ≤(d1+d2)2. | (3) |
Two-dimensional velocity fields are taken with the velocity components as ⌢U(⌢X,⌢Y,⌢t) and ⌢V(⌢X,⌢Y,⌢t) along the ⌢X−axis and ⌢Y−axis respectively.
The governing equations for an incompressible RPh fluid flow can be expressed as the following set of equations [43]:
∇.u=0, | (4) |
ρDuD⌢t=∇.Sij, | (5) |
(ρˉC)fDTD⌢t=Kf∇2T+Φ+(ρˉC)p(DB(∇C.∇T)+DT∇C.∇TTm), | (6) |
DCD⌢t=DTTm∇2T+DB∇2C, | (7) |
where, ρf signifies the base fluid density and ρp is the density of particles, u=u(⌢U(⌢X,⌢Y,⌢t),⌢V(⌢X,⌢Y,⌢t),0) indicate velocity vectors, ⌢t is for time, ˉCf is specific heat of base fluid and ˉCp is particles specific heat, Kf indicates the thermal conductivity, T is temperature, C is concentration, DB and DT are stands for Brownian and thermophoretic diffusion coefficients respectively. Sij and Φ are the Cauchy stress tensor and viscous dissipation terms and are given as [44]:
Sij=−⌢PI+⌢τij, | (8) |
Φ=(⌢τ⌢Y⌢Y−⌢τ⌢X⌢X)∂⌢V∂⌢Y+⌢τ⌢X⌢Y(∂⌢U∂⌢Y+∂⌢V∂⌢X). | (9) |
Where for the Reiner-Philippoff (RPh) fluid model, the stress strain relation is specified as [18]:
⌢τij=[μ∞+μ0−μ∞1+(12⌢τ20)(∑2l=1∑2m=1⌢τlm⌢τml)]⌢eij, | (10) |
with,
∑2l=1∑2m=1⌢τlm⌢τml=∑2l=1[⌢τl1⌢τ1l+⌢τl2⌢τ2l] |
= ⌢τ11⌢τ11+⌢τ21⌢τ12+⌢τ12⌢τ21+⌢τ22⌢τ22 |
= ⌢τ112+2⌢τ122+⌢τ222. |
Using ⌢τ12=⌢τ21.
The RPh fluid model (Eq (10)) contain three positive parameters μ∞,μ0 and ⌢τ0. Here μ∞ represents the higher Newtonian limiting viscosity, μ0 is the viscosity at zero shear stress rate and ⌢τ0 indicates the reference shear stress.
In two dimensional Cartesian coordinates system, the continuity Eq (4), momentum Eq (5), the thermal energy Eq (6) and the concentration Eq (7) in laboratory frame are expressed as:
∂⌢U∂⌢X+∂⌢V∂⌢Y=0, | (11) |
(ρf)(∂∂⌢t+⌢U∂∂⌢X+⌢V∂∂⌢Y)⌢U=−∂⌢P∂⌢X+∂⌢τ⌢X⌢X∂⌢X+∂⌢τ⌢X⌢Y∂⌢Y, | (12) |
(ρf)(∂∂⌢t+⌢U∂∂⌢X+⌢V∂∂⌢Y)⌢V=−∂⌢P∂⌢Y+∂⌢τ⌢Y⌢X∂⌢X+∂⌢τ⌢Y⌢Y∂⌢Y, | (13) |
ρfˉCf(∂T∂⌢t+⌢U∂T∂⌢X+⌢V∂T∂⌢Y)=Kf(∂2T∂⌢X2+∂2T∂⌢Y2)+[(⌢τ⌢Y⌢Y−⌢τ⌢X⌢X)∂⌢V∂⌢Y+⌢τ⌢X⌢Y(∂⌢U∂⌢Y+∂⌢V∂⌢X)] |
+ΓρfˉCf[DB(∂T∂⌢Y∂C∂⌢Y+∂C∂⌢X∂T∂⌢X)+DTTm(∂T∂⌢X+∂T∂⌢Y)2], | (14) |
(∂∂⌢t+⌢U∂∂⌢X+⌢V∂∂⌢Y)C=DTTm(∂2T∂⌢X2+∂2T∂⌢Y2)+DB(∂2C∂⌢X2+∂2C∂⌢Y2), | (15) |
where, Γ=ρpˉCpρfˉCf and Eq (10) in component form can be written as:
⌢τ⌢X⌢X=[μ∞+μ0−μ∞1+(12⌢τ20)(⌢τ⌢X⌢X2+2⌢τ⌢X⌢Y2+⌢τ⌢Y⌢Y2)]⌢e⌢X⌢X, | (16) |
⌢τ⌢X⌢Y=[μ∞+μ0−μ∞1+(12⌢τ20)(⌢τ⌢X⌢X2+2⌢τ⌢X⌢Y2+⌢τ⌢Y⌢Y2)]⌢e⌢X⌢Y, | (17) |
⌢τ⌢Y⌢Y=[μ∞+μ0−μ∞1+(12⌢τ20)(⌢τ⌢X⌢X2+2⌢τ⌢X⌢Y2+⌢τ⌢Y⌢Y2)]⌢e⌢Y⌢Y, | (18) |
where,
⌢e⌢X⌢X=2.∂⌢U∂⌢X,⌢e⌢Y⌢Y=2.∂⌢V∂⌢Y and ⌢e⌢X⌢Y=(∂⌢V∂⌢X+∂⌢U∂⌢Y). | (19) |
The subsequent transformations are used to change the investigated problem from stationary reference frame to moving frame [45]:
{⌢x=⌢X−c⌢t,⌢y=⌢Y,̑v(⌢x,⌢y)=⌢V(⌢X,⌢Y,⌢t)⌢p(⌢x,⌢y)=⌢P(⌢X,⌢Y,⌢t)⌢u(⌢x,⌢y)=⌢U(⌢X,⌢Y,⌢t)−c. | (20) |
And the following dimensionless quantities are used to convert the above dimensional system in dimensionless form:
{x=⌢xλ,y=⌢yd1,t=c⌢tλ,u=⌢uc,v=⌢vcδ,δ=d1λ,τij=d1⌢τijcμf,p=d21⌢pcμfλ,h1=⌢H1d1,h2=⌢H2d1,a=a1d1,b=b1d1,d=d2d1,ν=μfρf,Pr=μfCfKf,Re=cd1ρfμf,Ec=c2CfΔT,Br=Pr.Ec,θ=T−TmT1−T0,ϕ=C−CmC1−C0,Nb=ΓDB(C1−C0)ν,Nt=ΓDT(T1−T0)νTm. | (21) |
Here h1 and h2 represents the dimension free forms of both walls, ν is kinematics viscosity, d is the channel width ratio, ϕ indicates the concentration of nanoparticles, Nb is the parameter of Brownian motion, Nt is indicate the thermophoresis diffusion parameter, θ correspond to dimension free temperature, Pr, Re, Ec and Br are the Prandtl number, Reynolds number, Eckert number and Brinkman number respectively.
The dimensionless form of the governing equations in the wave frame in accounts of Eq (20) and the quantities (21) is as:
∂u∂x+.∂v∂y=0, | (22) |
δRe{[∂∂t+(u+1)∂∂x+v∂∂y](u+1)}=−∂p∂x+δ(∂τxx∂x)+∂τxy∂y, | (23) |
δ3Re{[∂∂t+(u+1)∂∂x+v∂∂y]v}=−∂p∂y+δ2(∂τyx∂x)+δ∂τyy∂y, | (24) |
δRePr{[∂θ∂t+(u+1)∂θ∂x+v∂θ∂y]} |
=(δ2∂2θ∂x2+∂2θ∂y2)+δBr(τyy−τxx)∂v∂y+Brτxy(∂u∂y+δ2∂v∂x) |
+PrNb(∂θ∂y∂ϕ∂y+δ2∂ϕ∂x∂θ∂x)+PrNb(δ2∂θ∂x+∂θ∂y)2, | (25) |
δRe(ΔC.Γ)[∂ϕ∂t+(u+1)∂ϕ∂x+v∂ϕ∂y]=Nt(δ2∂2θ∂x2+∂2θ∂y2)+Nb(δ2∂2ϕ∂x2+∂2ϕ∂y2), | (26) |
and corresponding dimensionless form of Eqs (16)-(18) are as under:
τxx=δ[1+λ1−11+12τ20(τ2xx+2τ2xy+τ2yy)]∂u∂x, | (27) |
τyy=δ[1+λ1−11+12τ20(τ2xx+2τ2xy+τ2yy)]∂v∂y, | (28) |
τxy=[1+λ1−11+12τ20(τ2xx+2τ2xy+τ2yy)]∂u∂y. | (29) |
In above equations λ1=μ0μ∞. Now by using the estimations of lengthy wavelength (δ<<1), small Reynolds number and by using the description of following stream functions:
u=∂ψ∂y and v=−∂ψ∂x. | (30) |
Equation (22) satisfies identically and Eqs (23)-(29) simplifies to the following form:
∂p∂x=∂∂yτxy, | (31) |
∂p∂y=0, | (32) |
∂2θ∂y2+Brτxy(∂2ψ∂y2)+PrNb(∂θ∂y∂ϕ∂y)+PrNt(∂θ∂y)2=0, | (33) |
NbNt(∂2ϕ∂y2)+∂2θ∂y2=0, | (34) |
τxx=0 and τyy=0, | (35) |
and
τxy=(1+λ1−11+(τ2xyτ20))∂2ψ∂y2. | (36) |
From Eqs (31) and (32) we have:
∂2∂y2τxy=0. | (37) |
Which depicts that τxy is a linear function of y and can be written as τxy=A1y+B1, where A1 and B1 are arbitrary parameters. And Eq (36) can be written as:
∂2ψ∂y2=τxy[γ+τ2xyτ2xy+γλ1]. | (38) |
This is explicit relationship between deformation rate and shear stress for RPh nanofluid. Where, λ1=μ0μ∞ and γ=τ20 are RPh fluid parameter and shear stress parameter respectively. For λ1=1 this relation converges to Newtonian fluid.
The dimension freeform of the walls h1(x) and h2(x) are is given as:
{h1(x)=1+k1x+acos[2πx],h2(x)=−d−k1x−bcos[2πx+ζ], | (39) |
where, k1=λtanαd1 is the non-uniformity parameter of the channel.
The instantaneous volumetric flow rate in the laboratory frame of reference is given by:
⌢Q(⌢X,⌢t)=∫⌢H1⌢H2⌢U(⌢X,⌢Y,⌢t)d⌢Y, | (40) |
where, ⌢H1 and ⌢H2 are the functions of ⌢X and ⌢t. Now using Eq (19), Eq (40) becomes:
⌢Q(⌢X,⌢t)=∫⌢h1⌢h2⌢u(⌢x,⌢y)d⌢y+c(⌢h1−⌢h2)=⌢q+c(⌢h1−⌢h2), | (41) |
where, ⌢h1(⌢x) and ⌢h2(⌢x) represents the upper, lower channel walls for wave frame and ⌢q(⌢x,⌢y) is the instantaneous volumetric flow rate in the moving frame and is given by:
⌢q(⌢x,⌢y)=∫⌢h1⌢h2⌢u(⌢x,⌢y)d⌢y. | (42) |
The time average flow rate (over one period T) is given by:
ˉQ=1T∫T0⌢Qd⌢t. | (43) |
Using Eq (41) in Eq (43) we have:
ˉQ=1T∫T0⌢qd⌢t+c(⌢h1−⌢h2) |
=⌢q+cd1+cd2+2c⌢xλtanα |
+ca1cos(2πλ⌢x)+cb1cos(2πλ⌢x+ζ). | (44) |
We note that h1(x) and h2(x) given in Eq (39) are the corresponding dimension free forms of the non-uniform peristaltic channel. So the corresponding dimension free form of Eq (44) is given as:
η=F+1+d+2k1x+acos(2πx)+bcos(2πx+ζ), | (45) |
where, η=ˉQcd1 and F=⌢qcd1 are the dimension free mean flow rates in the fixed and moving frame respectively and the expression of F is given as:
F=∫h1h2u(x,y)dy=∫h1h2∂ψ∂ydy=ψ(h1)−ψ(h2). | (46) |
The associated dimensionless boundary conditions are:
{ψ=F2,ψ'=−1,θ=−0.5,ϕ=−0.5 at upper wall h1(x),ψ=−F2,ψ'=−1,θ=0.5,ϕ=0.5 at lower wall h2(x). | (47) |
Equation (31) depicts that the expression of the pressure function is free from 'y' variable and Eq (36) demonstrates that τxy is a linear function of 'y'. Thus the solution of Eqs (32), (33) and (37) corresponding to the boundary conditions given in Eq (46) is retrieved numerically by the help of built in NDSolve function of MATHEMATICA software. Outcomes of the study are discussed thoroughly by plotting the graphs in rest of this section. We will talk about the effects of RPh nanofluid properties on physical quantities of interest like the axial velocity u, streamline functions ψ, concentration of nanoparticles ϕ, temperature profiles θ, heat transfer rate −θ'(h) and mass transportation rate −ϕ'(h) at the upper boundary and pressure gradient dpdx.
Graphical visualization of velocity profiles against dissimilar values of RPh fluid parameter λ1, shear stress parameter γ, for both pseudoplastic and dilatants fluids and arbitrary parameter A1 is shown in Figures 2-5. Reiner-Philippoff (RPh) fluid behaves like dilatants fluid for λ1<1, as pseudoplastic for λ1 > 1 and as Newtonian fluid for λ1=1. The influence of RPh fluid parameter λ1 on the fluid's velocity profiles is drawing in Figure 2 the blue line in the figure against the value λ1=1 is showing the velocity of Newtonian fluid. Further it is shown that the maximum velocity for pseudoplastic (shear thinning) fluid, corresponding to the value λ1>1, is larger than the velocity of dilatants fluid corresponding to the value of λ1<1. As RPh fluid proceeds as shear thinning (pseudoplastic) fluid for the value λ1>1, therefore its viscosity decreases as compared to shear thickening (dilatants) fluid and that's why the fluid's deformation rate increases and consequently velocity profile improves. RPh fluid model (see Eq (35)) behaves like Newtonian fluid for extreme shear rates. The same case is analyzed in Figures 3 and 4. From these two figures we can see that for γ→0 and for γ→∞. in both cases the fluid's velocity is same for both pseudoplastic (λ1=2.0) and dilatants (λ1=0.25) fluid. The other survey from these two graphs is that velocity is increasing function of shear stress parameter γ for pseudoplastic fluid and decreasing function for dilatants fluid. As in shear thinning case, fluid decreases its viscosity against applied shear stress and in case of shear thickening fluid, fluid gains viscosity that's why in case of pseudoplastic fluid velocity increases and in case of dilatants fluid velocity decreases against shear stress parameter. Figure 5 depicts that by escalating the value of arbitrary parameter A1 the fluid's velocity increase. Initially this increment is notable but after that there is no significant effect of arbitrary parameter A1 on velocity graph.
This section covers physical assessment of parameters like RPh fluid parameter λ1 shear stress parameter γ, Brownian motion parameter Nb, thermophoresis diffusion parameter Nt, Brinkman number Br and Prandtl number Pr on temperature profiles. As a result, Figures 6-12 are required to look for such an analysis briefly. Figure 6 demonstrates that by ascending the value of RPh fluid parameter λ1 the fluid's temperature increases. As at λ1=1 the fluid behaves likes the Newtonian fluid. And so the study shows that in same circumstances the temperature for dilatants fluid is less than the Newtonian's fluid temperature and Newtonian's fluid temperature is smaller than the pseudoplastic fluid's temperature. In Figures 7 and 8 the influences of shear stress parameter γ are deliberated on temperature profiles for pseudoplastic and dilatants fluids respectively. These graphs reveal that the temperature is an increasing function of shear stress parameter γ for pseudoplastic fluids and it behaves totally opposite for dilatants fluids for same values of shear stress parameter γ as decreasing function. The relevance of Brinkman's number on fluid hotness can be seen through Figure 9. This figure proves that fluid becomes hotter by enlarging the value of Brinkman number. Because the Brinkman number occurs as a result of viscous dissipation therefore, extra heat is created with greater Brinkman number, and therefore the temperature rises and hence temperature profile get enlarged. Also see reference [45]. Figure 10 represents that for greater Brownian motion parameter the temperature of the liquid increases. The consequence of Brownian scale tends to increase the relative mobility of nanoparticles compounds from the walls to the substance in response to which an increment in the temperature profile occurs. Figure 11 confirms a similar pattern of temperature distributions in relation to the thermophoresis parameter Nt. Physically greater Nt implies the occurrence of temperature differential phenomenon in which high-temperature nanoparticles transported to cooler region of the channel. As a result temperature is intensified. Figure 12 signifies that by increasing the amount of Prandtl number Pr the temperature of the nanofluid also increases. Because Prandtl number Pr has direct connection with specific heat of the liquid so increased specific heat raise the temperature of fluid.
The influence of sundry factors on the concentration of nanoparticles graphs are scrutinized through Figures 13 to 19. These graphs pointed out that the concentration of nanoparticles at the center of the channel is highly effected by the changing the values of implanted parameters. Figure 13 is plotted to show the effect of RPh fluid parameter λ1 on the concentration of nanoparticles. A decreasing tendency in nanoparticles mass transport rate is observed for enhancing values of λ1. From this we can conclude that concentration of nanoparticles for dilatants fluids is larger as compared to pseudoplastic fluids. Shear stress parameter γ have opposite concern on concentration for pseudoplastic and dilatants fluid which can be seen through Figures 14 and 15. For pseudoplastic fluids the concentration of nanoparticles is decreasing graph of shear stress parameter while for dilatants fluids it depicts increasing behavior. Figure 16 illustrates that with an increment in Brinkman number Br results with a decrease in nanoparticles mass transfer. Figures 17 and 18 reveal the antagonistic effects of Brownian motion parameter Nb and thermophoretic diffusions parameter Nt on concentration graphs. Declining impact is marked by ϕ for increasing thermophoretic diffusion parameter because with the strength of thermophoretic forces the action of diffusion is turned on which depressed the concentration profile of nanofluid. Contrary to this a rising effect in concentration is achieved by enhancing Brownian motion as nanoparticles attained higher densities for stronger Brownian motion which in turn develops concentration. Figure 19 depicts the Prandtl number Pr effects on concentration profile. From figure we can see that by growing values of Prandtl number Pr yields a depressed behavior of concentration profiles near the upper wall but this observation is not hold near lower wall where the significance of Pr is directly proportional to concentration profile.
In Figures 20-23 we will discussed about the effects of and on the pressure gradient profiles within channel. The common observation from these graphs is that the pressure gradient gains its maximum value at the narrowest section of the channel. Figure 20 demonstrates that by intensifying the value of the pressure gradient of RPh nanofluid increases. The rising behavior of pressure gradient field is larger in the contracted as compared to the broader section of the channel. In this figure the blue thick line is the representation of Newtonian fluid's pressure gradient. The lines upper the blue lines against the values correspond to the pseudoplastic fluid and the lines below the blue line against the values characterized the dilatants fluid. In Figure 21 the significance of on the pressure gradient is evaluated for shear thinning fluid. It is investigated from this sketch that pressure gradient of RPh fluid raises by increasing the value of. In Figure 22 the same significance of on pressure gradient is studied for the dilatants fluid. The influence of on pressure gradient for dilatants fluid is quite dissimilar as compared to the pseudoplastic fluid. In the case of dilatants, by increasing the value of the amount of pressure gradient reduces. Figure 23 examined that for, the height of every crest of the pressure gradient is same at narrow section of the channel. However, as we increase the value of parameter, that is, when we increase the width of channel, the height of each succeeding crest of the pressure gradient decreases than the preceding crest. This is because of by adding the value of the width of the channel increases and so the pressure gradient decreases in the wider part. This whole phenomenon can be seen in Figure 23.
This part consists of the physical study of parameters like RPh fluid parameter λ1, phase difference of the waves ζ, shear stress parameter γ and non-uniformity parameter k1 on stream line contours. From Figure 24, it is experimented that by improving the value of λ1 the size of the bolus increases. Bolus is the result of fluid's trapped volume by streamlines with the motion of the fluid, this generates an intriguing phenomenon known as trapping, the construction of an inner flowing bolus of nanofluid which moved together with the waves with wave speed at a fixed mean flow rate. Figure 25 indicates that phase difference of the channel has opposite effect on fluid trapped volume than the λ1 as by increasing the angle of phase difference, the volume of the bolus decreases. Figures 26 and 27 are drawn to understand the impact shear stress parameter γ on streamline functions for shear thinning fluids and shear thickening fluids respectively. From Figure 26 we can see that bolus for pseudoplastic fluids have increasing behavior with greater shear stress parameter. As this fluid is shear thinning fluid so by increasing shear stress parameter the fluid's viscosity decreases and fluid moves more freely, and more fluid's volume get trapped. Contrary to this fluid the size of bolus decreases in case of dilatants fluids for same values of shear stress parameter as we can observe from Figure 27. Figure 28 signifies that the amount of trapping volume of fluid get increases by increasing the non-uniformity of the channel and consequently the dimension of the bolus get increased.
Table 1 presents the numerical assessment of heat and mass transfer rates at the upper border of the channel for various physical parameters. First four lines of the table exhibits that heat transfer (−θ'(h1)) and mass transfer (−ϕ'(h1)) rates are increasing function of RPh fluid parameter λ1, which concludes that RPh fluid parameter improves heat and mass transfer process between the solid boundaries and the base fluid. Secondly, the table 1 predicts that in case of shear thinning fluids, heat (−θ'(h1)) and mass transfer (−ϕ'(h1)) rates raise for greater shear stress parameter and on the other hand for shear thickening fluid these transfer rates decline. Next four lines represent that heat transfer (−θ'(h1)) and mass transfer (−ϕ'(h1)) rates are increases by enlarging Nt. Moreover the ratio of mass transfer (−ϕ'(h1)) is much significant as compared to heat transfer (−θ'(h1)) rate. The table's following four lines anticipated that the Brownian motion parameter has opposing impacts on heat (−θ'(h1)) and mass transfer (−ϕ'(h1)) rates. As with increased Brownian motion parameter, the mass transfer rate (−ϕ'(h1)) drops drastically while the heat transfer rate (−θ'(h1)) increases steadily. The very next four lines of the table demonstrate that the Brinkman's number has a direct influence on heat (−θ'(h1)) and mass transfer (−ϕ'(h1)) rates. As by enlarging the Brinkman number, both the heat (−θ'(h1)) and mass transport (−ϕ'(h1)) rates improved. The last four lines of the table reveals that raising the Prandtl number improved the heat (−θ'(h1)) and mass transfer (−ϕ'(h1)) rates.
λ1 | γ | Nt | Nb | Br | Pr | −θ'(h1) | −ϕ'(h1) |
0.5 | 1.0 | 0.8 | 0.8 | 1.5 | 1.0 | 3.94339 | 3.25946 |
1.0 | 4.24148 | 3.55755 | |||||
1.5 | 4.52987 | 3.84593 | |||||
2.0 | 4.81133 | 4.12740 | |||||
1.5 | 0.5 | 4.39119 | 3.70726 | ||||
1.0 | 4.52987 | 3.84593 | |||||
1.5 | 4.64350 | 3.95957 | |||||
2.0 | 4.73910 | 4.05516 | |||||
0.5 | 0.5 | 4.09858 | 3.41464 | ||||
1.0 | 4.94339 | 3.25946 | |||||
1.5 | 4.81863 | 3.13470 | |||||
2.0 | 3.71590 | 3.03197 | |||||
1.5 | 1.0 | 0.4 | 4.30704 | 1.64057 | |||
0.8 | 4.52987 | 3.84593 | |||||
1.2 | 4.75514 | 6.27779 | |||||
1.6 | 4.98124 | 8.93658 | |||||
0.8 | 0.4 | 4.30703 | 7.58817 | ||||
0.8 | 4.52987 | 3.84593 | |||||
1.2 | 4.75514 | 2.60015 | |||||
1.6 | 4.98124 | 1.97764 | |||||
0.8 | 0.5 | 1.66990 | 1.28304 | ||||
1.0 | 3.28843 | 2.56540 | |||||
1.5 | 4.52987 | 3.84593 | |||||
2.0 | 5.81130 | 5.12737 | |||||
1.5 | 0.5 | 4.08841 | 3.40448 | ||||
1.0 | 4.52987 | 3.84593 | |||||
1.5 | 4.98124 | 4.29731 | |||||
2.0 | 5.43078 | 7.74984 |
Peristaltic flow of non-Newtonian based nanofluid in a tapered asymmetric channel is examined in this research. The governing equations incorporates dispersion model for the nanofluid and RPh fluid model for non-Newtonian fluid. Simplified form of governing equation based on long wavelength and small Reynolds number approximations are solved numerically. The consequential effects of different physical quantities on the velocity of the nanofluid, pressure gradient, temperature and concentration profiles are portrayed graphically. The most important outcomes of the study are given below:
• Velocity is increasing function of shear stress parameter for pseudoplastic fluid and decreasing function of shear stress parameter for dilatants fluid.
• It is noted that by increasing the value of RPh fluid parameter λ1 the temperature of the nanofluid increases.
• It is worth noting that both thermophoretic and Brownian diffusions causes temperature graphs to climb.
• A decreasing tendency in nanoparticles mass transfer is observed for enhancing values of RPh fluid parameter. From this we can conclude that concentration of nanoparticles for dilatants fluid is greater as compared to pseudoplastic fluid.
• Thermophoretic diffusion and Brownian motion factors have opposing influence on concentration of nanoparticles profiles. Concentration graph decreases by increasing thermophoretic diffusion parameter and increases by enhancing Brownian motion parameter.
• Further it is noticed that the addition of Brinkman number results in an enhancement in the nanofluid temperature.
• The pressure gradient of the nanofluid is higher in the occluded section of the channel as compared to the broader section of the channel.
• By increasing the non-uniformity of the channel, the pressure gradient of the nanofluid decreases as the width of the channel increases.
• The addition in the value of RPh fluid parameter λ1 the trapped volume of the nanofluid increases.
• It is also important to observe that the RPh fluid parameter enables heat and mass transfer rates to increase.
• Brownian motion parameter has opposing impacts on heat and mass transfer rates. As with increased Brownian motion parameter, the mass transfer rate drops drastically while the heat transfer rate increases steadily.
• It is analyzed from the study represented that heat transfer and mass transfer rates are increased by enlarging the thermophoretic diffusion parameter.
I hope this study will be step forward in the field of peristaltic flow, heat and mass transfer of non-Newtonian fluids and also provides new avenues of research for scholars to the way of RPh fluid academics.
This research work was funded by institutional fund projects under no. (IFP-A-2022-2-5-24). Therefore, authors gratefully acknowledge technical and financial support from the ministry of education and University of Hafr Al Batin, Saudi Arabia.
The authors declare no conflict of interest.
[1] |
M. Alquran, F. Yousef, F. Alquran, T. A. Sulaiman, A. Yusuf, Dual-wave solutions for the quadratic-cubic conformable-Caputo time-fractional Klein-Fock-Gordon equation, Math. Comput. Simulat., 185 (2021), 62–76. https://doi.org/10.1016/j.matcom.2020.12.014 doi: 10.1016/j.matcom.2020.12.014
![]() |
[2] | A. Hussanan, M. Z. Ismail, Samiulhaq, I. Khan, S. Sharidan, Radiation effect on unsteady MHD free convection flow in a porous medium with Newtonian heating, Int. J. Appl. Math. Stat., 42 (2013), 474–480. |
[3] |
A. Khan, K. Ali Abro, A. Tassaddiq, I. Khan, Atangana–Baleanu and Caputo Fabrizio analysis of fractional derivatives for heat and mass transfer of second grade fluids over a vertical plate: a comparative study, Entropy, 19 (2017), 279. https://doi.org/10.3390/e19080279 doi: 10.3390/e19080279
![]() |
[4] | A. A. Shaikh, S. Qureshi, Comparative analysis of Riemann Liouville, Caputo-Fabrizio, and Atangana-Baleanu integrals, J. Appl. Math. Comput. Mech., 21 (2022), 91–101. |
[5] | M. Caputo, M. Fabrizio, A new definition of fractional derivative without singular kernel, Progr. Fract. Differ. Appl., 1 (2015), 1–13. |
[6] |
A. Atangana, D. Baleanu, New fractional derivatives with nonlocal and non-singular kernel: theory and application to heat transfer model, Therm. Sci., 20 (2016), 763–769. https://doi.org/10.48550/arXiv.1602.03408 doi: 10.48550/arXiv.1602.03408
![]() |
[7] | A. Atangana, D. Baleanu, Caputo-Fabrizio derivative applied to groundwater flow within confined aquifer, J. Eng. Mech., 143 (2016), D4016005. |
[8] |
A. Atangana, I. Koca, Chaos in a simple nonlinear system with Atangana-Baleanu derivatives with fractional order, Chaos Soliton. Fract., 89 (2016), 447–454. https://doi.org/10.1016/j.chaos.2016.02.012 doi: 10.1016/j.chaos.2016.02.012
![]() |
[9] |
O. J. J. Algahtani, Comparing the Atangana-Baleanu and Caputo-Fabrizio derivative with fractional order: Allen Cahn model, Chaos Soliton. Fract., 89 (2016), 552–559. https://doi.org/10.1016/j.chaos.2016.03.026 doi: 10.1016/j.chaos.2016.03.026
![]() |
[10] |
H. Shatha, Atangana–Baleanu fractional framework of reproducing kernel technique in solving fractional population dynamics system, Chaos Soliton. Fract., 133 (2020), 109624. https://doi.org/10.1016/j.chaos.2020.109624 doi: 10.1016/j.chaos.2020.109624
![]() |
[11] |
A. Atangana, J. F. Gómez-Aguilar, Numerical approximation of Riemann–Liouville definition of fractional derivative: from Riemann–Liouville to Atangana–Baleanu, Numer. Meth. Part. Differ. Equ., 34 (2018), 1502–1523. https://doi.org/10.1002/num.22195 doi: 10.1002/num.22195
![]() |
[12] |
T. Abdeljawad, M. A. Hajji, Q. M. Al-Mdallal, F. Jarad, Analysis of some generalized ABC-fractional logistic models, Alex. Eng. J., 59 (2020), 2141–2148. https://doi.org/10.1016/j.aej.2020.01.030 doi: 10.1016/j.aej.2020.01.030
![]() |
[13] |
H. Abboubakar, P. Kumar, N. A. Rangaig, S. Kumar, A malaria model with Caputo–Fabrizio and Atangana–Baleanu derivatives, Int. J. Model. Simul. Sci., 12 (2021), 2150013. https://doi.org/10.1142/S1793962321500136 doi: 10.1142/S1793962321500136
![]() |
[14] |
T. Sitthiwirattham, R. Gul, K. Shah, I. Mahariq, J. Soontharanon, K. J. Ansari, Study of implicit-impulsive differential equations involving Caputo-Fabrizio fractional derivative, AIMS Math., 7 (2022), 4017–4037. https://doi.org/10.3934/math.2022222 doi: 10.3934/math.2022222
![]() |
[15] |
D. Baleanu, S. S. Sajjadi, A. Jajarmi, Z. Defterli, On a nonlinear dynomical system with both chaotic and nonchaotic behaviors: a new fractional analysis and control, Adv. Differ. Equ., 2021 (2021), 234. https://doi.org/10.1186/s13662-021-03393-x doi: 10.1186/s13662-021-03393-x
![]() |
[16] |
D. Baleanu, S. S. Sajjadi, J. H. Asad, A. Jajarmi, E. Estiri, Hyperchaotic behaviors, optimal control and synchronization of a nonautonomous cardiac conduction system, Adv. Differ. Equ., 2021 (2021), 175. https://doi.org/10.1186/s13662-021-03320-0 doi: 10.1186/s13662-021-03320-0
![]() |
[17] |
D. Baleanu, S. Zibaei, M. Namjoo, A. Jajarmi, A nonstandard finite difference scheme for the modeling and nonidentical synchronization of a noval fractional chaotic system, Adv. Differ. Equ., 2021 (2021), 308. https://doi.org/10.1186/s13662-021-03454-1 doi: 10.1186/s13662-021-03454-1
![]() |
[18] |
R. B. Hetnarski, J. Ignaczak, Generalized thermoelasticity, J. Therm. Stress., 22 (1999), 451–476. https://doi.org/10.1080/014957399280832 doi: 10.1080/014957399280832
![]() |
[19] |
H. W. Lord, Y. Shulman, A generalized dynamical theory of thermoelasticity, J. Mech. Phys. Solids, 15 (1967), 229–309. https://doi.org/10.1016/0022-5096(67)90024-5 doi: 10.1016/0022-5096(67)90024-5
![]() |
[20] |
A. E. Green, K. A. Lindsay, Thermoelasticity, J. Elasticity, 2 (1972), 1–7. https://doi.org/10.1007/BF00045689 doi: 10.1007/BF00045689
![]() |
[21] |
S. Chen, F. Liu, V. Anh, A novel implicit finite difference method for the one-dimensional fractional percolation equation, Numer. Algor., 56 (2011), 517–535. https://doi.org/10.1007/s11075-010-9402-0 doi: 10.1007/s11075-010-9402-0
![]() |
[22] |
R. Metzler, J. Klafter, The random walk's guide to anomalous diffusion: a fractional dynamics approach, Phys. Rep., 339 (2000), 1–77, 2000. https://doi.org/10.1016/S0370-1573(00)00070-3 doi: 10.1016/S0370-1573(00)00070-3
![]() |
[23] |
I. Podlubny, A. Chechkin, T. Skovranek, Y. Chen, B. M. Vinagre Jara, Matrix approach to discrete fractional calculus. Ⅱ. partial fractional differential equations, J. Comput. Phys., 228 (2009), 3137–3153. https://doi.org/10.1016/j.jcp.2009.01.014 doi: 10.1016/j.jcp.2009.01.014
![]() |
[24] |
H. Jafari, A. Golbabai, S. Seifi, K. Sayevand, Homotopy analysis method for solving multi-term linear and nonlinear diffusion-wave equations of fractional order, Comput. Math. Appl., 59 (2010), 1337–1344. https://doi.org/10.1016/j.camwa.2009.06.020 doi: 10.1016/j.camwa.2009.06.020
![]() |
[25] |
S Momani, Z. Odibat, Homotopy perturbation method for nonlinear partial differential equations of fractional order, Phys. Lett. A, 365 (2007), 345–350. https://doi.org/10.1016/j.physleta.2007.01.046 doi: 10.1016/j.physleta.2007.01.046
![]() |
[26] |
J. S. Duan, M. Li, Y. Wang, Y. L. An, Approximate solution of fractional differential equation by quadratic splines, Fractal Fract., 6 (2022), 369. https://doi.org/10.3390/fractalfract6070369 doi: 10.3390/fractalfract6070369
![]() |
[27] | S. K. Lydia, M. M. Jancirani, A. A. Anitha, Numerical solution of nonlinear fractional differential equations using kharrat-toma iterative method, Nat. Volatiles Essent. Oils, 8 (2021), 9878–9890. |
[28] |
N. A. Zabidi, Z. A. Majid, A. Kilicman, Z. B. Ibrahim, Numerical solution of fractional differential equations with Caputo derivative by using numerical fractional predict–correct technique, Adv. Cont. Discr. Mod., 2022 (2022), 26. https://doi.org/10.1186/s13662-022-03697-6 doi: 10.1186/s13662-022-03697-6
![]() |
[29] |
H. Wang, F. Wu, D. Lei. A novel numerical approach for solving fractional order differential equations using hybrid functions, AIMS Math., 6 (2021), 5596–5611. https://doi.org/10.3934/math.2021331 doi: 10.3934/math.2021331
![]() |
[30] |
Z. F. Bonab, M. Javidi, Higher order methods for fractional differential equation based on fractional backward differentiation formula of order three, Math. Comput. Simul., 172 (2020), 71–89. https://doi.org/10.1016/j.matcom.2019.12.019 doi: 10.1016/j.matcom.2019.12.019
![]() |
[31] |
A. E. Green, P. M. Naghdi, A Re-examination of the basic postulates of thermomechanics, P. Roy. Soc. A Math. Phy., 432 (1991), 171–194. https://doi.org/10.1098/rspa.1991.0012 doi: 10.1098/rspa.1991.0012
![]() |
[32] |
A. E. Green, P. M. Naghdi, Thermoelasticity without energy dissipation, J. Elasticity, 31 (1993), 189–208. https://doi.org/10.1007/BF00044969 doi: 10.1007/BF00044969
![]() |
[33] |
. E. Green, P. M. Naghdi, On undamped heat waves in an elastic solid, J. Therm. Stresses, 15 (1992), 253–264. https://doi.org/10.1080/01495739208946136 doi: 10.1080/01495739208946136
![]() |
[34] |
D. Y. Tzou, A unified approach for heat conduction from macro- to micro-scales, J. Heat Transfer., 117 (1995), 8–16. https://doi.org/10.1115/1.2822329 doi: 10.1115/1.2822329
![]() |
[35] |
D. Y. Tzou, The generalized lagging response in small-scale and high-rate heating, Int. J. Heat Mass Transf., 38 (1995), 3231–3240. https://doi.org/10.1016/0017-9310(95)00052-B doi: 10.1016/0017-9310(95)00052-B
![]() |
[36] | D. Y. Tzou, Macro-to microscale heat transfer: the lagging behavior, New York: Taylor & Francis, 1997. |
[37] |
S. K. Roychoudhuri, On a thermoelastic three-phase-lag model, J. Therm. Stresses, 30 (2007), 231–238. https://doi.org/10.1080/01495730601130919 doi: 10.1080/01495730601130919
![]() |
[38] | A. E. Abouelregal, On Green and Naghdi thermoelasticity model without energy dissipation with higher order time differential and phase-lags, J. Appl. Comput. Mech., 6 (2020), 445–456. |
[39] |
A. E. Abouelregal, Two-temperature thermoelastic model without energy dissipation including higher order time-derivatives and two phase-lags, Mater. Res. Express, 6 (2019), 116535. https://doi.org/10.1088/2053-1591/ab447f doi: 10.1088/2053-1591/ab447f
![]() |
[40] | A. E. Abouelregal, Generalized mathematical novel model of thermoelastic diffusion with four phase lags and higher-order time derivative, Eur. Phys. J. Plus, 135 (2020), 263. |
[41] |
A. E. Abouelregal, A novel generalized thermoelasticity with higher-order time-derivatives and three-phase lags, Multidiscip. Model. Ma., 16 (2019), 689–711. https://doi.org/10.1108/MMMS-07-2019-0138 doi: 10.1108/MMMS-07-2019-0138
![]() |
[42] |
A. E. Abouelregal, Ö. Civalek, H. F. Oztop, Higher-order time-differential heat transfer model with three-phase lag including memory-dependent derivatives, Int. Commun. Heat Mass, 128 (2021), 105649. https://doi.org/10.1016/j.icheatmasstransfer.2021.105649 doi: 10.1016/j.icheatmasstransfer.2021.105649
![]() |
[43] |
Y. Z. Povstenko, Fractional heat conduction equation and associated thermal stress, J. Therm. Stresses, 28 (2004), 83–102. https://doi.org/10.1080/014957390523741 doi: 10.1080/014957390523741
![]() |
[44] |
Y. Z. Povstenko, Fractional radial heat conduction in an infinite medium with a cylindrical cavity and associated thermal stresses, Mech. Res. Commun., 37 (2010), 436–440. https://doi.org/10.1016/j.mechrescom.2010.04.006 doi: 10.1016/j.mechrescom.2010.04.006
![]() |
[45] |
Y. Povstenko, Non-axisymmetric solutions to time-fractional diffusion-wave equation in an infinite cylinder, Fract. Calc. Appl. Anal., 14 (2011), 418–435. https://doi.org/10.2478/s13540-011-0026-4 doi: 10.2478/s13540-011-0026-4
![]() |
[46] |
Y. Povstenko, T. Kyrylych, Fractional thermoelasticity problem for an infinite solid with a penny-shaped crack under prescribed heat flux across its surfaces, Phil. Trans. R. Soc. A, 378 (2020), 20190289. https://doi.org/10.1098/rsta.2019.0289 doi: 10.1098/rsta.2019.0289
![]() |
[47] |
Y. Povstenko, T. Kyrylych, B. Woźna-Szcześniak, R. Kawa, A. Yatsko, An external circular crack in an infinite solid under axisymmetric heat flux loading in the framework of fractional thermoelasticity, Entropy, 24 (2022), 70. https://doi.org/10.3390/e24010070 doi: 10.3390/e24010070
![]() |
[48] |
Y. Qiao, X. Wang, H. Qi, H. Xu, Numerical simulation and parameters estimation of the time fractional dual-phase-lag heat conduction in femtosecond laser heating, Int. Commun. Heat Mass, 125 (2021), 105355. https://doi.org/10.1016/j.icheatmasstransfer.2021.105355 doi: 10.1016/j.icheatmasstransfer.2021.105355
![]() |
[49] |
Y. J. Yu, L. J. Zhao, Fractional thermoelasticity revisited with new definitions of fractional derivative, Eur. J. Mech. A Solid., 84 (2020), 104043. https://doi.org/10.1016/j.euromechsol.2020.104043 doi: 10.1016/j.euromechsol.2020.104043
![]() |
[50] |
Y. J. Yu, Z. C. Deng, Fractional order theory of Cattaneo-type thermoelasticity using new fractional derivatives, Appl. Math. Model., 87 (2020), 731–751. https://doi.org/10.1016/j.apm.2020.06.023 doi: 10.1016/j.apm.2020.06.023
![]() |
[51] |
Z. Xue, J. Liu, X. Tian, Y. Yu, Thermal shock fracture associated with a unified fractional heat conduction, Eur. J. Mech. A Solids, 85 (2021), 104129. https://doi.org/10.1016/j.euromechsol.2020.104129 doi: 10.1016/j.euromechsol.2020.104129
![]() |
[52] |
Y. Yu, Z. C. Deng, New insights on microscale transient thermoelastic responses for metals with electron-lattice coupling mechanism, Eur. J. Mech. A Solid., 80 (2020), 103887. https://doi.org/10.1016/j.euromechsol.2019.103887 doi: 10.1016/j.euromechsol.2019.103887
![]() |
[53] |
Y. Yu, Z. C. Deng, Fractional order thermoelasticity for piezoelectric materials, Fractals, 29 (2021), 2150082. https://doi.org/10.1142/S0218348X21500821 doi: 10.1142/S0218348X21500821
![]() |
[54] | C. Li, F. Zeng, Numerical methods for fractional calculus, Boca Raton: CRC Press, 2019. |
[55] |
A. Atangana, On the new fractional derivative and application to nonlinear Fisher's reaction–diffusion equation, Appl. Math. Comput., 273 (2016), 948–956. https://doi.org/10.1016/j.amc.2015.10.021 doi: 10.1016/j.amc.2015.10.021
![]() |
[56] |
A. E. Abouelregal, M. Alesemi, Vibrational analysis of viscous thin beams stressed by laser mechanical load using a heat transfer model with a fractional Atangana-Baleanu operator, Case Stud. Therm. Eng., 34 (2022), 102028. https://doi.org/10.1016/j.csite.2022.102028 doi: 10.1016/j.csite.2022.102028
![]() |
[57] | N. Sarkar, A. Lahiri, Eigenvalue approach to two-temperature magneto-thermoelasticity, Vietnam J. Math. Math., 40 (2012), 13–30. |
[58] |
A. Sur, Nonlocal memory-dependent heat conduction in a magneto-thermoelastic problem, Wave. Random Complex, 32 (2020), 251–271. https://doi.org/10.1080/17455030.2020.1770369 doi: 10.1080/17455030.2020.1770369
![]() |
[59] |
A. E. Aboueregal, H. M. Sedighi, The effect of variable properties and rotation in a visco-thermoelastic orthotropic annular cylinder under the Moore–Gibson–Thompson heat conduction model, P. I. Mech. Eng. L J. Mat., 235 (2021), 1004–1020. https://doi.org/10.1177/14644207209858 doi: 10.1177/14644207209858
![]() |
[60] |
A. E. Abouelregal, R. Alanazi, H. M. Sedighid, Thermal plane waves in unbounded nonlocal medium exposed to a moving heat source with a non-singular kernel and higher order time derivatives, Eng. Anal. Bound. Elem., 140 (2022), 464–475. https://doi.org/10.1016/j.enganabound.2022.04.032 doi: 10.1016/j.enganabound.2022.04.032
![]() |
[61] |
G. Honig, U. Hirdes, A method for the numerical inversion of Laplace transform, J. Comput. Appl. Math., 10 (1984), 113–132. https://doi.org/10.1016/0377-0427(84)90075-X doi: 10.1016/0377-0427(84)90075-X
![]() |
[62] |
H. Dubner, J. Abate, Numerical inversion of Laplace transforms by relating them to the finite Fourier cosine transform, J. Assoc. Comp. Mach., 15 (1968), 115–123. https://doi.org/10.1145/321439.321446 doi: 10.1145/321439.321446
![]() |
[63] |
F. R. De Hoog, J. H. Knight, A. N. Stokes, An improved method for numerical inversion of Laplace transforms, SIAM J. Sci. Statist. Comput., 3 (1982), 357–366. https://doi.org/10.1137/0903022 doi: 10.1137/0903022
![]() |
[64] |
A. E. Abouelregal, H. M. Sedighi, Magneto-thermoelastic behaviour of a finite viscoelastic rotating rod by incorporating Eringen's theory and heat equation including Caputo–Fabrizio fractional derivative, Eng. Comput., 2022. https://doi.org/10.1007/s00366-022-01645-2 doi: 10.1007/s00366-022-01645-2
![]() |
[65] |
T. He, L. Cao, S. Li, Dynamic response of a piezoelectric rod with thermal relaxation, J. Sound Vib., 306 (2007), 897–907. https://doi.org/10.1016/j.jsv.2007.06.018 doi: 10.1016/j.jsv.2007.06.018
![]() |
[66] | K. Cole, J. Beck, A. Haji-Sheikh, B. Litkouhi, Heat conduction using green's functions, 2nd, New York: Taylor & Francis, 2010. |
[67] |
Z. B. Hou, R. Komanduri, General solutions for stationary/moving plane heat source problems in manufacturing and tribology, Int. J. Heat Mass, 43 (2000), 1679–1698. https://doi.org/10.1016/S0017-9310(99)00271-9 doi: 10.1016/S0017-9310(99)00271-9
![]() |
[68] | R. Viskanta, T. L. Bergman, Heat transfer in materials processing, New York: McGraw-Hill, 1998. |
[69] |
G. Araya, G. Gutierrez, Analytical solution for a transient, three-dimensional temperature distribution due to a moving laser beam, Int. J. Heat Mass, 49 (2006), 4124–4131. https://doi.org/10.1016/j.ijheatmasstransfer.2006.03.026 doi: 10.1016/j.ijheatmasstransfer.2006.03.026
![]() |
1. | Bahram Jalili, Amirhossein Rezaeian, Payam Jalili, Fathollah Ommi, Davood Domiri Ganji, Numerical modeling of magnetic field impact on the thermal behavior of a microchannel heat sink, 2023, 45, 2214157X, 102944, 10.1016/j.csite.2023.102944 | |
2. | Iskandar Waini, Nurul Amira Zainal, Najiyah Safwa Khashi'ie, Khairum Bin Hamzah, Abdul Rahman Mohd Kasim, Anuar Ishak, Ioan Pop, Mixed convection of MHD Reiner–Philippoff fluid flow past a vertical shrinking plate with radiative heat transfer, 2023, 05779073, 10.1016/j.cjph.2023.04.001 | |
3. | J. Iqbal, F. M. Abbasi, I. Ali, Heat transfer analysis for magnetohydrodynamic peristalsis of Reiner–Philippoff fluid: Application of an artificial neural network, 2024, 36, 1070-6631, 10.1063/5.0207600 | |
4. | Rajib Kumar Mandal, Hiranmoy Maiti, Samir Kumar Nandy, Bioconvective MHD flow of Williamson nanofluid past an expandable Riga wedge in the presence of activation energy, mass suction and velocity slip, 2023, 1040-7782, 1, 10.1080/10407782.2023.2263155 | |
5. | M. Tahir, A. Ahmad, Y. Khan, R. Nawaz, S. Ahmad, Role of minimization entropy generation and heat transfer on the nanofluid peristaltic flow via Reiner-Philippoff model , 2024, 1040-7790, 1, 10.1080/10407790.2024.2354493 | |
6. | M.P. Mkhatshwa, M. Khumalo, Double diffusion and Hall effects on MHD sinusoidal natural convection flow of silver water-based nanofluid from a porous vertical plate, 2023, 7, 26668181, 100516, 10.1016/j.padiff.2023.100516 | |
7. | Shuguang Li, Sania Naseer, Nazia Shahmir, Muhammad Ramzan, Shalan Alkarni, Seifedine Kadry, Influence of magnetic dipole and thermophoretic particle deposition on the flow of Reiner–Philippoff fluid with Thompson and Troian slip, 2024, 0044-2267, 10.1002/zamm.202400098 | |
8. | M. Alkinidri, Jamshaid Iqbal, F. M. Abbasi, Analysis of thermal transport and wall stresses for MHD peristaltic motion of Reiner‐Philippoff fluid with viscous dissipation and mixed convection effects, 2025, 105, 0044-2267, 10.1002/zamm.202400329 |
λ1 | γ | Nt | Nb | Br | Pr | −θ'(h1) | −ϕ'(h1) |
0.5 | 1.0 | 0.8 | 0.8 | 1.5 | 1.0 | 3.94339 | 3.25946 |
1.0 | 4.24148 | 3.55755 | |||||
1.5 | 4.52987 | 3.84593 | |||||
2.0 | 4.81133 | 4.12740 | |||||
1.5 | 0.5 | 4.39119 | 3.70726 | ||||
1.0 | 4.52987 | 3.84593 | |||||
1.5 | 4.64350 | 3.95957 | |||||
2.0 | 4.73910 | 4.05516 | |||||
0.5 | 0.5 | 4.09858 | 3.41464 | ||||
1.0 | 4.94339 | 3.25946 | |||||
1.5 | 4.81863 | 3.13470 | |||||
2.0 | 3.71590 | 3.03197 | |||||
1.5 | 1.0 | 0.4 | 4.30704 | 1.64057 | |||
0.8 | 4.52987 | 3.84593 | |||||
1.2 | 4.75514 | 6.27779 | |||||
1.6 | 4.98124 | 8.93658 | |||||
0.8 | 0.4 | 4.30703 | 7.58817 | ||||
0.8 | 4.52987 | 3.84593 | |||||
1.2 | 4.75514 | 2.60015 | |||||
1.6 | 4.98124 | 1.97764 | |||||
0.8 | 0.5 | 1.66990 | 1.28304 | ||||
1.0 | 3.28843 | 2.56540 | |||||
1.5 | 4.52987 | 3.84593 | |||||
2.0 | 5.81130 | 5.12737 | |||||
1.5 | 0.5 | 4.08841 | 3.40448 | ||||
1.0 | 4.52987 | 3.84593 | |||||
1.5 | 4.98124 | 4.29731 | |||||
2.0 | 5.43078 | 7.74984 |
λ1 | γ | Nt | Nb | Br | Pr | −θ'(h1) | −ϕ'(h1) |
0.5 | 1.0 | 0.8 | 0.8 | 1.5 | 1.0 | 3.94339 | 3.25946 |
1.0 | 4.24148 | 3.55755 | |||||
1.5 | 4.52987 | 3.84593 | |||||
2.0 | 4.81133 | 4.12740 | |||||
1.5 | 0.5 | 4.39119 | 3.70726 | ||||
1.0 | 4.52987 | 3.84593 | |||||
1.5 | 4.64350 | 3.95957 | |||||
2.0 | 4.73910 | 4.05516 | |||||
0.5 | 0.5 | 4.09858 | 3.41464 | ||||
1.0 | 4.94339 | 3.25946 | |||||
1.5 | 4.81863 | 3.13470 | |||||
2.0 | 3.71590 | 3.03197 | |||||
1.5 | 1.0 | 0.4 | 4.30704 | 1.64057 | |||
0.8 | 4.52987 | 3.84593 | |||||
1.2 | 4.75514 | 6.27779 | |||||
1.6 | 4.98124 | 8.93658 | |||||
0.8 | 0.4 | 4.30703 | 7.58817 | ||||
0.8 | 4.52987 | 3.84593 | |||||
1.2 | 4.75514 | 2.60015 | |||||
1.6 | 4.98124 | 1.97764 | |||||
0.8 | 0.5 | 1.66990 | 1.28304 | ||||
1.0 | 3.28843 | 2.56540 | |||||
1.5 | 4.52987 | 3.84593 | |||||
2.0 | 5.81130 | 5.12737 | |||||
1.5 | 0.5 | 4.08841 | 3.40448 | ||||
1.0 | 4.52987 | 3.84593 | |||||
1.5 | 4.98124 | 4.29731 | |||||
2.0 | 5.43078 | 7.74984 |