Abstract
Viscoelastic fluids, such as polymers, paints, and DNA suspensions, are almost everywhere and very useful in the industry. This article aims to study the significance of ramped temperature in the dynamics of viscoelastic fluid. Magnetohydrodynamic (MHD) effect is considered in the presence of Lorentz force. The flow is considered in a porous medium using generalized Darcy’s law. Heat transfers through convection, and the fluid near the plate takes heat in a ramped nature. Instead of the classical fluid model which has certain limitations, a generalized model is considered with fractional derivatives of the Atangana–Baleanu type. The well-known technique of Laplace transform was adopted to obtain the solutions which are displayed in various plots with detailed discussion analysis. From the graphical analysis, it is worth noting that the Atangana–Baleanu fractional model shows a good memory effect on the dynamics of the viscoelastic fluid as compared to its classical form.
1 Introduction
Viscoelastic fluids (VFs) form a subclass of non-Newtonian fluids such as the blend of a solvent and some polymer. Other examples include paints, DNA suspensions, some biological fluids, and other products from the chemical industry. One of the well-known examples from the list of VFs is an Oldroyd-B fluid (OBF). Due to the non-linear relationship between stress and rate of strain, OBF has the same nature as that of a non-Newtonian fluid. The idea of OBF was introduced in 1950, in the pioneering work [1]. The OBF has many applications in fluid dynamics. This model is the generalization of the Maxwell model. Many researchers analyzed OBFs based on their real-world applications: the authors in Ref. [2] studied unsteady OBF over a vertical plate in a porous medium. The authors in Ref. [3] analyzed the influence of generalized OBFs using numerical simulations. The authors in Ref. [4] developed analytical and numerical solutions of a two-dimensional multi-term time‐fractional OBF model. The authors in Ref. [5] considered generalized OBF flow between two infinite parallel rigid plates and numerically simulated the results. In the literature, the flow of Oldroyd-B is less investigated using ramped wall velocity and temperature. These fluids have many applications in different practical situations and real-world problems [6]. Some researchers tried to investigate ramped wall velocity and temperature in their studies. The authors in Ref. [7] examined natural convection in polyethylene glycol-based molybdenum disulfide nanofluid under ramped wall condition and other additional effects. The authors in Ref. [8] discussed ramped wall temperature and ramped wall concentration with the Hall effect in an unsteady MHD flow of a second-grade fluid through a porous medium. The authors in Ref. [9] developed a general solution of MHD free convection flows over a vertical plate with ramped wall temperature. The authors in Ref. [10] examined entropy generation in MHD flow over a vertical plate embedded in a porous medium. The authors in Ref. [11] considered ramped wall temperature conditions and studied irreversibility analysis in time-dependent flow past a plate (vertical) with arbitrary wall shear stress. The authors in Ref. [12] studied MHD flow over an inclined plate under a conjugate effect and ramped wall temperature condition in a porous medium.
Recently, many fractional derivative operators have been used for different scientific reasons in the field of fluid dynamics and in many other dynamical systems. Some of these applications are found in fractional-order neurons for parameter estimation, fractional viscoelasticity models, fractional single-phase-lag models of heat conduction, and the eigenproblem of molecular alignment, as discussed in Ref. [13] and in many other references. Sene [14, 15] used fractional derivatives and developed two different fluid models. He obtained exact analytical solutions and plotted results graphically with a detailed discussion analysis. In Ref. [16], the authors studied fractional derivatives together with their various applications in reservoir engineering problems. In Ref. [17], the authors investigated radial basis functions using a fractional derivative approach. In Ref. [18],the authors applied fractional derivatives in the field of physics and modern sciences. In Ref. [19],the authors discussed the theory of fractional derivatives and its application in mathematics, physics, chemistry, and engineering. In Ref. [20],the authors introduced new fractional derivatives (AB fractional derivatives) with non-local and no-singular kernel. They presented some useful and important properties of the new derivative with application in a fractional heat transfer model. After that, this new idea was applied to several other problems. In Ref. [21], the authors applied the Atangana–Baleanu fractional derivative in the Caputo sense to the convective flow of CMC-based CNT nanofluid in a vertical microchannel. In Ref. [22],the authors studied the magnetic field effect on the convection flow of Newtonian viscous fluid past a moving plate such that its temperature is constant at the boundary and its concentration depends on time. The problem is modeled using the Caputo–Fabrizio time fractional derivative. In Ref. [23], the authors used the Caputo time-fractional derivative and examined the natural convection flow through a vertical cylinder. In Ref. [24], Sheikh et al. provided a comparative analysis of Caputo–Fibrizio and Atangana–Baleanu derivatives for a generalized Casson fluid model with heat generation and chemical reaction.
In this work, the idea of fractional derivatives (Atangana–Baleanu fractional derivative in the Caputo sense) is applied to develop a fractional model of OBF. The OBF is considered as electrically conducting and passing through a porous medium. The flow is considered over a vertical plate with ramped heating, which sufficiently influences thermal analysis and the flow itself due to the involvement of convection term in momentum equation (25). This article is arranged in the following sections. Section 2 includes the mathematical modeling of the problem. The governing equations are derived from the constitutive equations of OBF. Initial and boundary conditions are defined. The next section (Section 3) includes solution methodology with a basic definition of the AB fractional derivatives. The solution of energy equation is given in sub-Section 3.1, and the solution of momentum equation is given in sub-Section 3.2. The results are computed and plotted in Section 4 with a detailed discussion analysis. At the end, conclusion remarks are added in Section 5.
2 Mathematical Modeling
This report focuses on the dynamics of incompressible unsteady OBF over a vertical plate. The flow is considered in a porous medium with the magnetic field applied in a perpendicular direction. The equations of continuity and momentum are given as follows:
Here, is the velocity vector and is the Cauchy stress tensor T.
The velocity field (for one-dimensional and uni-directional flow) is given as follows:
The corresponding Cauchy stress tensor is given as follows:
Here, represent the indeterminate stress tensor and extra stress tensor, respectively. In addition, extra stress tensor is given by the following relation:where show dynamic viscosity, retardation time, and relaxation time, respectively. Furthermore, material time derivative and Rivlin–Ericksen tensor can be written as follows:
In the given study, we have considered Oldroyd-B fluid; therefore, modified Darcy’s law can be written in the following form:where represents the porous medium and represents permeability of the porous medium.
From Maxwell equations,Here, is the total electric field. By Ohm’s law (generalized form),
The cross product with the magnetic field givesHere, is the applied magnetic field and is the induced magnetic field by polarization (perturbation produced by fluid motion).
Eq. 12 further reduces to
In view of Eqs 3–14, equations take the following form:
Similarly, the temperature equation in the presence of thermal radiation and heat generation is given as follows:
Here, are heat capacitance, radiation, thermal conductivity, and heat generation constant, respectively. Eq. 16 is the same as Eq. 8 in Ref. [24].
The radiative heat flux is given as follows:
Using the Taylor series expansion to linearized given in Eq. 17, we get
Substituting the result obtained in Eq. 18 into Eq. 17 yields
Substituting Eq. 19 in Eq. 16 gives
The following constraints on the physical model are imposed. Initially, when time , the fluid and plate were at rest at initial temperature . At time , both the temperature and velocity are higher or lower than and when , respectively. This condition is also known as ramped condition for temperature and velocity. But when , the temperature and velocity remain unchanged. In other words, for small values of time, the temperature and velocity are referred to as ramped velocity and temperature, while for greater values of time, the velocity and temperature are known as isothermal velocity and temperature.
Under the aforementioned assumption, the following initial and boundary conditions are defined:
The following dimensionless variables are used for dimensional analysis:
Introducing the aforementioned equation in Eq. 15, 20, the following non-dimensional system is obtained:with transformed initial and boundary conditionsand the non-dimensional parameters and constants are as follows:
3 Research Methodology: Atangana–Baleanu Fractional Derivative
To transform the OBF classical model into the Atangana–Baleanu (AB) fractional model of order , it is required to first define the Atangana–Baleanu fractional derivative as follows:
Definition: The Atangana–Baleanu (AB) time fractional derivative in the Caputo sense, with fractional order , is defined as [20, 21]where is the normalization function, and
Here, is the generalized form of the Mittag–Leffler function, which is defined as follows:
The Laplace transform (LT) of the AB fractional derivative in the Caputo sense is given as follows:
3.1 Solution of Energy Equation
The non-dimensional classical heat Eq. 24 is transformed into the fractional form by incorporating the AB fractional derivative of order as given in the following:
Note that Eq. 30 is similar to Eq. 12 in Ref. [24]. The Laplace transform is applied to Eq. 30, taking into account the initial condition from Eq. 25, the following transformed equation is obtained:
Eq. 31 can be written in a more suitable form as follows:
The solution of Eq. 32 giveswhere and
Now using the boundary conditions from Eq. 26, Eq. 33 gives the following solution:
After the Laplace inversion, Eq. 34 giveswhere
The symbol shows the convolution product.
3.2 Solution of Momentum Equation
The fractional form of Eq. 23, after using the Laplace transform method, reduces to the following form:equivalently
The aforementioned equation can be written in more appropriate forms as follows:where ; ; and .
The solution of the homogeneous part of Eq. 38 isand the corresponding non-homogeneous part is
Eqs 39, 40 give the total solution asWith the help of the boundary condition given in Eq. 26, and after finding constants, the final solution is
The inverse Laplace transform of Eq. 42 is with
4 Results and Discussion
This section provides the graphical analysis of OBF with MHD and porous medium effects with ramped wall velocity and temperature. The problem is formulated in terms of Atangana–Baleanu fractional derivatives. To find the analytical solutions, the Laplace transform technique is applied. All the flow parameters are discussed by graphical analysis. The solutions obtained for the Atangana–Baleanu fractional OBF model have been discussed, and the influence of all parameters is shown. The influence of Grashof number Gr, magnetic number M, Prandtl number Pr, radiation number Nr, and Atangana–Baleanu β is studied in various plots.
The influence of fractional parameters on velocity profile is highlighted in Figure 1. It is found that an increase in fractional parameters results in a decrease in velocity for both cases of ramped and isothermal temperature. This figure clearly shows that as , the fractional model reduces to the classical OBF model.
FIGURE 1
Influence of the AB fractional operator on velocity for ramped and isothermal wall temeprature.
The impact of on velocity profiles is shown in Figure 2. From this figure, it is clear that increasing results in an increase in the buoyancy force, which enables to speed up the fluid motion, thus resulting in an increase in the velocity of OBF. Figure 3 shows the influence of on AB fractional velocity profiles. From this figure, it can be seen that increasing results in a decrease in the OBF velocity. Note that the Prandtl number variation is shown only to see its effect on velocity, and its chosen values do not correspond to a specific non-Newtonian liquid. The effect of magnetic parameter on OBF fractional velocity is highlighted in Figure 4. The influence of on velocity in both cases of ramped and constant temperature increases, resulting in a decrease in the velocity. However, the variation in velocity profiles in case of constant temperature is more visible. It is true physically as increasing leads to an increase in the Lorentz forces, and as a result, the retardation force increases and hence the fluid velocity decreases.
FIGURE 2
Influence of Gr on velocity for constant and ramped temperature.
FIGURE 3
Influence of Pr on velocity for constant and ramped temperature.
FIGURE 4
Influence of M on velocity for constant and ramped temperature.
The impact of porosity parameter is shown in Figure 5. It is found that velocity increases with larger values of porosity parameter. The impact of on the fluid velocity is same as expected. With increasing porosity, the retardation impact decreases to a more permeable surface, and hence the velocity increases. Figure 6 shows the influence of relaxation time on the velocity profile for OBF AB fractional velocity. It is found that by increasing , the magnitude of velocity profile increases. The influence of retardation time is highlighted in Figure 7. From this figure, it is found that by increasing , the magnitude of the velocity profile decreases.
FIGURE 5
Influence of K on velocity for constant and ramped temperature.
FIGURE 6
Influence of on velocity for constant and ramped temperature.
FIGURE 7
Influence of λ1 on velocity for constant and ramped temperature.
The influence of the AB fractional operator on temperature distributions is highlighted in Figure 8. From this figure, it is clear that with the increase in the fractional parameter in the case of the AB derivative, the temperature of the OBF decreases. It is worth noting that for , the temperature goes to the classical form, which provides a comparison between the classical and fractional forms. Figure 9 shows the influence of Pr on temperature distribution. From this graph, it is clear that with increasing Pr, the temperature of the fluid decreases. The influence of Nr on temperature is highlighted in Figure 10. From this figure, it is found that with the increase in the values of Nr, the temperature of the fluid increases. Figure 11 shows the influence of H1 on the temperature profile. It is found that an increase in H1 results in an increase in OBF temperature.
FIGURE 8
Impact of the AB fractional operator on temperature.
FIGURE 9
Impact of Pr on temperature (AB fractional model).
FIGURE 10
Impact of Nr on temperature (AB fractional model).
FIGURE 11
Impact of H1 on temperature (AB fractional model).
5 Conclusion
This study aimed to investigate the dynamics of OBF flowing over a vertical plate with ramped heating and time-dependent velocity. The effect of MHD is considered under the Lorentz force. Additional effects of porosity, thermal radiation, and heat generation are also considered. The fractional model of OBF was first developed using Atangana–Baleanu fractional derivatives and then the well-known technique of Laplace transform was adopted to obtain the solutions, which are displayed in various plots with detailed discussion analysis. The significance of rising ramped temperature on the dynamics of the unsteady viscoelastic fluid subject to Lorentz force is concluded with the following main points.
• Atangana–Baleanu fractional parameter reduces OBF velocity in both isothermal wall temperature and ramped heating; however, in ramped heating, the velocity profiles are closer than in the isothermal heating case.
• Increasing , the OBF velocity increases in both isothermal wall temperature and ramped heating.
• Increasing , the OBF velocity decreases in both isothermal wall temperature and ramped heating.
• The greater values of increase the OBF temperature in both isothermal wall temperature and ramped heating.
• The greater values of decrease the temperature of the fluid in both isothermal wall temperature and ramped heating.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.
Author contributions
IK formulated, solved the problem, and wrote the manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Nomenclature
Extra stress tensor
Specific heat (pressure constant)
Skin friction coefficient
Cauchy stress tensor
Material time derivative
Rivlin–Ericksen tensor
Mean spectral absorption coefficient
Permeability
Radiation number
Prandtl number
Radiative energy flux
Transform parameter
Temperature
Temperature away from the boundary
Temperature at the boundary ()
Velocities along the x-axis
Uniform velocity
Co-ordinates (space and time)
Acceleration due to gravity
- H0
Heat generation coefficient in the dimensional form
- H1
Heat generation coefficient in the non-dimensional form
Grashof number (
Greek letters
Fractional parameter (
Volumetric coefficient of thermal expansion
Retardation time
Relaxation time ()
Thermal conductivity
Dynamic viscosity
Kinematic viscosity
Density of the fluid
Stefan–Boltzmann constant
Non-dimensional velocity
Non-dimensional temperature
Non-dimensional space variable (-)
References
1.
OldroydJG. On the Formulation of Rheological Equations of State. In: Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences (1950) 200(1063):523–41.
2.
ZhaoJZhengLZhangXLiuFChenX. Unsteady Natural Convection Heat Transfer Past a Vertical Flat Plate Embedded in a Porous Medium Saturated with Fractional Oldroyd-B Fluid. J Heat Transfer (2017) 139(1). 10.1115/1.4034546
3.
LiuYLiuFFengLXinB. Novel Numerical Analysis for Simulating the Generalized 2D Multi-Term Time Fractional Oldroyd-B Fluid Model. arXiv preprint arXiv:1903.07816 (2019). 10.48550/arXiv.1903.07816
4.
ZhangJLiuFAnhVV. Analytical and Numerical Solutions of a Two‐Dimensional Multi‐Term Time‐Fractional Oldroyd‐B Model. Numer Methods Partial Differential Eq (2019) 35(3):875–93. 10.1002/num.22327
5.
FengLLiuFTurnerIZhuangP. Numerical Methods and Analysis for Simulating the Flow of a Generalized Oldroyd-B Fluid between Two Infinite Parallel Rigid Plates. Int J Heat Mass Transfer (2017) 115:1309–20. 10.1016/j.ijheatmasstransfer.2017.08.105
6.
FengLLiuFTurnerIZhengL. Novel Numerical Analysis of Multi-Term Time Fractional Viscoelastic Non-Newtonian Fluid Models for Simulating Unsteady MHD Couette Flow of a Generalized Oldroyd-B Fluid. Fractional Calculus Appl Anal (2018) 21(4):1073–103. 10.1515/fca-2018-0058
7.
AliFArifMKhanISheikhNSaqibM. Natural Convection in Polyethylene Glycol Based Molybdenum Disulfide Nanofluid With Thermal Radiation, Chemical Reaction and Ramped Wall Temperature. Int J Heat Technology (2018):619–631. 10.18280/ijht.360227
8.
VeeraKrishnaMChamkhaAJ. Hall Effects on Unsteady MHD Flow of Second Grade Fluid Through Porous Medium with Ramped wall Temperature and Ramped Surface Concentration. Phys Fluids (2018) 30(5):053101. 10.1063/1.5025542
9.
ShahNAZafarAAAkhtarS. General Solution for MHD-Free Convection Flow Over a Vertical Plate with Ramped Wall Temperature and Chemical Reaction. Arab J Math (2018) 7(1):49–60. 10.1007/s40065-017-0187-z
10.
KhanAKhanIAliFShafieS. A Note on Entropy Generation in MHD Flow Over a Vertical Plate Embedded in a Porous Medium with Arbitrary Shear Stress and Ramped Temperature. J Porous Media (2016) 19(2):175. 10.1615/jpormedia.v19.i2.50
11.
KhanAul KarimFKhanIAliFKhanD. Irreversibility Analysis in Unsteady Flow Over a Vertical Plate with Arbitrary Wall Shear Stress and Ramped Wall Temperature. Results Phys (2018) 8:1283–90. 10.1016/j.rinp.2017.12.032
12.
KhanAKhanIAliFShafieS. Effects of wall Shear Stress on MHD Conjugate Flow over an Inclined Plate in a Porous Medium with Ramped wall Temperature. Math Probl Eng (2014) 2014:15. 10.1155/2014/861708
13.
YavuzMSeneNYıldızM. Analysis of the Influences of Parameters in the Fractional Second-Grade Fluid Dynamics. Mathematics (2022) 10:1125. 10.3390/math10071125
14.
SeneN. Fractional Model and Exact Solutions of Convection Flow of an Incompressible Viscous Fluid under the Newtonian Heating and Mass Diffusion. J Mathematics (2022) 2022:8785197. 10.1155/2022/8785197
15.
SeneN. Analytical Solutions of a Class of Fluids Models with the Caputo Fractional Derivative. Fractal Fract (2022) 6:35. 10.3390/fractalfract6010035
16.
ObembeADAl-YousefHYHossainMEAbu-KhamsinSA. Fractional Derivatives and Their Applications in Reservoir Engineering Problems: a Review. J Pet Sci Eng (2017) 157:312–27. 10.1016/j.petrol.2017.07.035
17.
LiuQZhuangPLiuFLaiJAnhVChenS. An Investigation of Radial Basis Functions for Fractional Derivatives and Their Applications. Comput Mech (2020) 65(2):475–86. 10.1007/s00466-019-01779-z
18.
HilferR. Applications of Fractional Calculus in Physics. Singapore: World Scientific (2000).
19.
YangXJ. General Fractional Derivatives: Theory, Methods and Applications. Chapman and Hall: CRC Press (2019): 380.
20.
AtanganaABaleanuD. New Fractional Derivatives with Nonlocal and Non-singular Kernel: Theory and Application to Heat Transfer Model. Therm Sci (2016) 20(2):763–9. 10.2298/tsci160111018a
21.
SaqibMKhanIShafieS. Application of Atangana-Baleanu Fractional Derivative to MHD Channel Flow of CMC-Based-CNT's Nanofluid through a Porous Medium. Chaos, Solitons and Fractals (2018) 116:79–85. 10.1016/j.chaos.2018.09.007
22.
ShahNAElnaqeebTAnimasaunIL. Insight into the Natural Convection Flow through a Vertical Cylinder Using Caputo Time-Fractional Derivatives. Int J Appl Comput Math (2018) 4:80. 10.1007/s40819-018-0512-z
23.
ShahNAAhmadHHajizadehAZebMAhmadSMahsudYet alEffect of Magnetic Field on Double Convection Flow of Viscous Fluid over a Moving Vertical Plate with Constant Temperature and General Concentration by Using New Trend of Fractional Derivative. Open J Math Sci (2018) 2018:253–65. 10.30538/oms2018.0033
24.
SheikhNAAliFSaqibMKhanIJanSAAAlshomraniASet alComparison and Analysis of the Atangana-Baleanu and Caputo-Fabrizio Fractional Derivatives for Generalized Casson Fluid Model with Heat Generation and Chemical Reaction. Results Phys (2017) 7:789–800. 10.1016/j.rinp.2017.01.025
25.
AnimasaunILShahNAWakifAMahantheshBSivarajRKorikoOK. Ratio of Momentum Diffusivity to Thermal Diffusivity: Introduction, Meta-Analysis, and Scrutinization. New York: Chapman and Hall/CRC.
Summary
Keywords
viscoelastic fluid, ramped heating, MHD, fractional model, heat transfer, porous medium
Citation
Khan I (2022) Significance of Ramped Temperature in the Dynamics of Unsteady Viscoelastic Fluid Subjected to Lorentz Force. Front. Phys. 10:924910. doi: 10.3389/fphy.2022.924910
Received
20 April 2022
Accepted
09 May 2022
Published
07 July 2022
Volume
10 - 2022
Edited by
Animasaun I. L., Federal University of Technology, Nigeria
Reviewed by
Ndolane Sene, Cheikh Anta Diop University, Senegal
Muhammad Junaid, Nanjing Normal University, China
Updates
Copyright
© 2022 Khan.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ilyas Khan, i.said@mu.edu.sa
This article was submitted to Interdisciplinary Physics, a section of the journal Frontiers in Physics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.