Original Research ARTICLE
Computational and Physical Examination About the Aspects of Fluid Flow Between Two Coaxially Rotated Disks by Capitalizing Non-fourier Heat Flux Theory: Finite Difference Approach
- 1Department of Mathematics, Air University, Islamabad, Pakistan
- 2College of Computer and Information Sciences, Majmaah University, Al Majma'ah, Saudi Arabia
- 3Department of Mathematics, College of Arts and Sciences, Prince Sattam Bin Abdulaziz University, Wadi Aldawaser, Saudi Arabia
- 4Department of Mathematics, City University of Science and Information Technology, Peshawar, Pakistan
- 5Department of Mathematics, College of Sciences, King Khalid University, Abha, Saudi Arabia
This pagination is executed to exemplify flow features exhibited by viscous fluid between two coaxially rotated disks. Thermal analysis is performed by using Cattaneo-Christov heat flux theory. Porosity aspects are also taken into account. Mathematically structured non-linear PDEs are transmuted into non-linear ODEs by employing Karman transformations. Afterward, solution is heeded by applying implicit finite difference scheme renowned as Keller box method. Interpretation of flow controlling parameters on axial, tangential, and radial components of velocity, thermal distribution is exhibited. Assurance of computed data is done by managing comparison for skin friction coefficients at walls of disks. From the attained outcomes, it is addressed that the magnitude of axial and radial velocities diminishes at lower disk contrary to upper disk for intensifying magnitude of Reynolds number. Increment in tangential component of velocity is also demonstrated for uplifts values of Reynolds number. It is also concluded that thermal field decrements for increasing of Pr and thermal relaxation parameter. It is worthy to mention that shear drag coefficient at wall of lower disk decreases conversely to the wall shear coefficient magnitude at wall of upper disk.
Rotational fluid flow generated by coaxial disks is one of the classical problems of fluid mechanics. In recent years, it has become a popular research area and has persuaded researchers due to magnificent theoretical and practical significance in engineering and applied sciences. Some important practical fields in which rotatory flow is capitalized are rotor–stator system, gas turbine engineering, air rotational cleaners, medical equipment, chemical engineering, and thermal power–generating systems. In view of its capitalization in various processes, researcher fraternity is examining such type of flows in current days. Inaugurated work on flow induced due to rotating disk is performed by Karman . He introduced transformations and provided a mathematical framework for construction of ordinary differential systems of rotational flows from Navier Stokes theory. Cochran  also used these transformations to scrutinize rotating disk flow by using numerical integration scheme. Batchelor  validated that Karman transformation can be evenly used for fluid flow between two coaxial rotating disks. Rotating flow by two coaxial disks is primarily examined by Stewartson . Chapple and Stokes  elucidated the flow features of fluid between two coaxially rotated disks. Mellor et al.  bestowed comprehensive treatment of fluid flow restricted between two coaxial infinite disks, one rotating, and other stationary. Thermal aspects of fluid between rotational disks were discussed by Arora and Stokes . Interpretation of flow phenomenon between porous stationary disk and solid rotating disk was manipulated by Kumar et al. . Xun et al.  considered rotating disk of variable thickness and adumbrated the flow features of Power law fluid. Hall effects on an unsteady MHD (magneto hydrodynamics) flow of viscous incompressible electrically conducting fluid between two rotating disks with non-coincident parallel axes embedded in a porous viscous medium were accorded by Das et al. . Asgher et al.  conducted Lie group analysis on the thermal features of fluid manifested by rotating disks. Elmaboud et al.  discussed peristaltic flow induced by sinusoidal wave propagating with constant speed on the walls of two-dimensional infinite rotating channel by heeding semi-analytical solutions.
In the most recent couple of decades, researcher fraternity has shown fantastic energy in exploring the heat propagation by means of a wave mechanism rather than essentially by diffusion. Late studies affirm that this is not just a low-temperature phenomenon but heat transfer mechanism also occurs at high temperature through diffusion. Just about 200 years prior, thermal features in various circumstances and especially in flowing fluid environment were interpreted by Fourier law of heat conduction . However, this law is inadequate in comprehending complete description about the heat exchange procedure among multiple connected surfaces in various conditions because of its disablement to fulfill the principle of causality. Later on, in 1948, Cattaneo  modified Fourier law by viewing the inadequateness generated by Fourier law of heat conduction and explored that this law explains the thermal attribute at low temperature because it generates parabolic heat equation in which initial disturbance are felt throughout the domain. After getting thorough analysis about Fourier law and viewing vector field aspect of heat flux, he included thermal relaxation time term to control generated thermal inertia, which is known as Maxwell-Cattaneo law. Afterward, Christov  proposed that objective time derivative instead of material time derivative is used for exact fulfillment of causality principle. He changed the time derivative in Maxwell-Cattaneo model by Oldryod upper convective derivative, which has successfully preserved the material invariant formulation and famously known as Cattaneo-Christov heat flux law. Cattaneo-Christov heat flux model has bounteous applications in engineering and modern industrial procedures like in skin burns and nanofluids, cooling of electronic devices, food technology, nuclear reactor cooling, power generation, heat exchangers, heat propagation in tissues, and so many. The uniqueness and stability of the solution for governing temperature equations by Cattaneo-Christov model in some initial and boundary value problems were proven by Straughan . Additionally, steadiness of structure of Cattaneo-Christov heat flux model with uniqueness was revealed by Ciarletta and Straughan . Tibullo and Zampoli  explicated the behavior of Cattaneo-Christov heat flux model in incompressible fluid flows. Aqsa  and Haddad  heeded numerical solution for thermal convection of an incompressible viscous fluid by obliging Cattaneo-Christov heat flux model. Mekheimer and Elmaboud  interpreted the aspects of temperature-dependent viscosity and thermal conductivity on peristaltic flow of a Newtonian fluid in a vertical asymmetric channel. Mekneimer  addressed heat transfer features of peristaltic couple stress fluid in asymmetric channel generated by wave with different phase and amplitudes. All of the abovementioned thought-provoking investigations have generated prodigious interest of researchers toward the analysis of flow in the presence of thermal aspects [23–26].
Transport procedures through porous space are commonly encountered in various chemical, mechanical, geophysical, electrochemical, and metallurgical routines. The theory about macroscopic movement of fluid in porous medium comprises differential equation that expresses linear relation between velocity and pressure gradient. Initially, Henry Darcy  (in 1856) presented a law to explicate the dynamic phenomenon in porous medium by working on the flow of sandy water through pebbles. Several technological processes depend on porous media theory, such as hydrology, oil exploration, solar collectors, porous insulations, packed beds, chromatography, heterogeneous catalysis, control of shear stresses at the seabed bottom, and oscillatory flow through seabed ripples. Darcy theory has promising applications in the field of biomedicine and the development of biological clogging and flow through tissues . Granular material  where significant amount of pore structures exist has application in manufacturing, paper, ceramic products, and textiles. Taseer et al.  addressed the flow behavior of Maxwell nanofluid in porous medium by implementing zero mass flux condition. They capitalized on Darcy-Forchheimer law to depict the flow pattern. They found that porosity parameter mounts the magnitude of temperature and concentration of particles. Seddeek  probed convective heat transfer in fluid immersed in porous medium. Analytical results for Darcy flow was described by Jha and Kaurangini . Aziz et al.  computed the traveling wave solution for the time-dependent viscoelastic fluid by way of a porous flat plate.
Magnetohydrodynamics is the study of the interaction between magnetic field and conductive fluid. The essence about magnetization is that the external magnetic field controls the turbulence in flow field. In addition, the magnetized flows differ from ordinary fluids because the generated current in the bulk fluid produces volumetric Lorentz force that extensively modifies the features. In recent years, magnetization and its impact on flow features have attained pervasive focus due to its extraordinary industrial applications, such as magnetized materials processes, manufacturing of glass, and MHD controlled electric generators. So the analysis of application of magnetic field has experienced great development and diversity. Andersson  performed exclusive study by manipulating electromagnet hydrodynamic waves mathematically. The stretched flow of two-dimensional Newtonian fluid under the effects of applied magnetic field was contemplated by Andersson . Liu  extended the work of Andersson  and described the heat and mass transfer of MHD viscous fluid flow over stretching surface. He computed exact solution of the problem by following the procedure of Andersson . The impact of normally impinging magnetic field on boundary layer flow of Newtonian fluid over permeable stretching sheet was analyzed by Kumaran et al. . Yirga and Tesfay  developed the numerical simulations for MHD viscous fluid flow over non-linear stretching sheet. The fluid flow equations were solved via Keller-Box method, and variations in physical quantities were presented regarding different parametric conditions. Recently, Yasin et al.  simulated the problem of two-dimensional MHD viscous nanofluid flow over porous stretched sheet. The formulated equations were solved by implementing well-known shooting technique. Mabood et al.  developed the approximate analytic solution of MHD boundary layer fluid flow over exponentially stretching surface. Some of the literature regarding the mentioned aspects is accessed through the references [40–43].
Present disquisition is addressed to excogitate thermophysical features exhibited in viscous fluid flow between two coaxially rotating disks embedded in permeable medium by obliging Cattaneo-Christov heat flux law. According to author's knowledge and available literature survey, it is found that very concise work is done so far in this direction. Tremendous engineering and practical application generated by rotating disk flows make present analysis highly potential. The authors hope that this manuscript will serve as a reference study for future researches. The article is strategized in such a way that the literature assessment is presented in section Introduction, whereas the mathematical structuring is provided in section Mathematical Model. The explanation about the solution methodology is debated in section Numerical Procedure. Comprehensive analysis and interpretation of flow controlling parameters are disclosed in section Results and Discussion. Last, the outcome norms are listed in section Conclusions.
Consider a steady, incompressible flow of viscous fluid between two coaxially rotated disks. The lower disk is placed at z = 0, whereas the distance between the disks is h units. Lower and upper disks possess angular velocities Ω1 and Ω2, respectively, and a1 and a2 are corresponding stretching rates (Figure 1). Porous medium between disks is considered, and Cattaneo-Christov heat flux model is obliged to analyze thermal features of fluid flow model.
We have used cylindrical coordinates (r, θ, z) with velocity components to the velocity profile and temperature equations as follows:
with boundary conditions:
where Equations (3–5) are referred to Hayat et al. , also pressure is expressed as and are the temperatures of upper and lower disks, and flux of heat satisfies:
where γ is the thermal relaxation parameter (It is defined as the parameter that controls the speed of heat waves produced within the system and makes them move with finite speed to follow the principle of causality), and k is the thermal conductivity. Now, we omit q from the Equations (5, 7) and obtain:
Equations (2–5) and Equation (8) are transformed into ordinary differential equations by obliging Von Karman transformations :
Mass conservation is identically satisfied, and Equations (2–4, 6, 8) take the following form:
where Re denotes Reynolds number, Pr is the Prandtl number, A1 and A2 are scaled stretching parameters, λ is the thermal relaxation, τ and β are rotational number and porosity parameter.
To make a simpler form of Equation (10), we removed ϵ.
The pressure parameter ϵ can be found by using Equations (10 and 14) as:
Equation (17) vanishes due to the given initial conditions in Equation (14).
The radial and tangential components of shear stress at lower disk are τzr and τzθ
where τw is the total shear stress, which is defined as:
Cf1 and Cf2 are the skin friction coefficients at lower and upper disks defined as:
where is the local Reynolds number.
Manipulation of accurate solution is necessary for physical interpretation of current work. Initially, equations are modeled by using Karman approximation and afterward, we have attained an intricate system of ordinary differential Equations (10–13) along with boundary conditions in Equation (14). We have applied Keller-box scheme referred to [40, 44–46], that is, the implicit finite difference scheme. For the implementation of this technique, first, we have to transform it into a system of first-order equations and define new variables u(y, ζ), v(y, ζ),
and Equations (11–13, 16) are reduced to
Similarly, the boundary conditions are converted into the following forms
Average and center difference gradients at the point of net derivatives are demarcated in Figure 2 and mathematically as below
Applying the Newton iteration fk+1 = fk+δfk, for all dependent variables involved in linearized non-linear algebraic equations and substituting these expressions in non-linear equations and neglecting quadratic and higher order terms in δ, a linear tridiagonal system is presented as follows:
with boundary conditions are:
Now, we consist the tridiagonal block matrices of given linearized equations in the form:
In Equation (28), the elements are defined as:
These computations are repeated until some convergence criteria are satisfied.
Results and Discussion
Current segment is dedicated to elucidate the numerical and graphical impact of velocity parameters involved and temperature profile and coefficient of skin friction. Variation in axial velocity against Reynolds number is adorned in Figure 3. It is found that axial velocity at lower disk decays with increment in Re. The reason behind this fact is that Re has direct relation with inertial forces. Therefore, with increase of Re, inertial effects dominate and cause the velocity of lower disk to decelerate. Figure 4 expresses the impact of Re on radial component of velocity. It is explored that with the increase of Re, velocity of lower disk diminishes and upper disk uplifts. It is seen from the picture that half of its portion from 0.0 to 0.6 represents velocity pattern of lower disk and from 0.6 to 1.0, it discloses radial velocity for upper plate. Impact of scaled stretching parameter on axial component of velocity at lower disk is exemplified in Figure 5. It is evidenced that with the increment of A1, f(ξ) mounts at lower disk because stretching rate is decreasing continuously. Figure 6 is portrayed to manifest the variation in radial against A1 velocity for lower and upper disks. It is found that radial component of velocity increases for the lower disk as compared to upper disk. This is due to the stretching rate of the lower disk that is continuously increasing and upper disk decrements. It is exhibited in Figure 7 that f′(ξ) decays with A2 for lower disk and mounts in the case of upper disk. The justification behind this impact is that stretching rate of the upper disk is more than that of the lower disk. The behavior of axial velocity with stretching scaled parameter A1of upper disk is depicted in Figure 8. For larger values of A2, the axial velocity of fluid decrements near the upper disk and upsurging behavior is noticed at the lower disk. By increasing A2, velocity magnitude along radial direction in the vicinity of the upper disk increases, so velocity along axial direction as an outcome depreciates. The impact of Re on tangential component of velocity g(ξ) is disclosed in Figure 9. It is fond that with increase of Re, g(ξ) suppresses. By increasing Re, inertial forces increase so more velocity is induced by the inertial forces. Behavior of tangential velocity g(ξ) against A2 is stretched in Figure 10. It is found that g(ξ) decrements and f′(ξ) uplifts for A2. With increase of A2, stretching rate of the upper disk increases, and as an outcome, axial velocity increases and tangential velocity decreases. The variation in g(ξ) with rotational parameter is revealed in Figure 11. Positive attribute in g(ξ) is observed against τ. With the increase of rotational parameter τ, centrifugal force is induced, which as an outcome uplifts the tangential component of velocity. Curves investigating the aspects of β on tangential component of velocity are adorned in Figure 12. It is justified by the fact that momentum equation is a dimensionless parameter, so with the increase of β, momentum profile is tangential and its direction diminishes. Positive impact on thermal distribution against Re is observed in Figure 13. With increase of Re, viscous forces decrement and velocity of fluid particle increases. Thus, the temperature is defined as average motion of fluid molecules so thermal field molecules due to uplifts of movement of particles. The impact of θ(ξ) against β is anticipated and shown in Figure 14. Increase in thermal magnitude is observed if β is increased. Since increase in β raises the rotation of disks, by increasing the rotation of disk, more rotational motion in fluid is generated and as a consequence kinetic energy of fluid molecule increases. The increase in kinetic motion raises the temperature profile. Variation in thermal profile by varying A1 in the range of (0.0 ≤ A1 ≤ 1.5) is revealed in Figure 15. It is observed that temperature of fluid boosts against the values of A1. As we increase the A1, the stretching rate increases, the fluid particles between disks exceeds, and hence temperature boosts up. The impact of Prandtl number Pr on θ(ξ) is exhibited in Figure 16. It is found that thermal distribution decreases with Pr. This is due to the fact that Pr is the ratio of viscous diffusion to thermal diffusion. Thus, by increasing Pr, thermal diffusion decreases so temperature decreases. Figure 17 is adorned to study the impact of thermal relaxation parameter on thermal distribution. Declined attribute in temperature against thermal relaxation parameter λ is depicted. It is because of the fact that with the increase of λ, fluid particles will take more time to transfer heat to its neighboring particles, thus the temperature decreases.
Table 1 numerically discloses the influence of porosity parameter β, stretching parameters A1 and A2, and Reynolds number Re. The skin friction coefficient increases for greater value of Re and stretching parameters A1 and A2, whereas it decreases for increasing values of porosity parameter β and rotating parameter τ at the upper and lower disks. Table 2 gives assurance of present work by constructing comparison with previously published literature for skin friction coefficient along radial and tangential components. Here, τ ≥ 0 shows the rotation of both disks in the same direction, τ ≤ 0 represents the direction of rotation of both disks in opposite direction, and τ = 0 means upper disk is fixed.
Current exertion is devoted to analyze the impact of Cattaneo-Christov heat flux theory on fluid flow between the two parallel rotating disks. Equations are modeled in the form of partial differential equations and then transformed into ordinary differential expressions. These ODE (ordinary differential equations) are tackled by Keller-box scheme. The key findings are summarized as follows:
• At the lower disk, the radial and axial velocity profile increases for maximum value of A1 while the same effects at upper disk for greater A2.
• For rotational and stretching parameters, the tangential velocity profile increases at disk with variation of parameters.
• Thermal effects are reduced for both thermal relaxation and Prandtl number.
• The skin friction coefficient at both disks is less for greater value of rotational parameter.
Data Availability Statement
All datasets generated for this study are included in the article/supplementary material.
All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.
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.
The authors extend their appreciation to the Deanship of Scientific Research at Majmaah University for funding this work under Project Number (RGP-2019-28).
9. Xun S, Zhao J, Zheng L, Chen X, Zhang X. Flow and heat transfer of Ostwald-de Waele fluid over a variable thickness rotating disk with index decreasing. Int J Heat Mass Transf . (2016) 103:1214–24. doi: 10.1016/j.ijheatmasstransfer.2016.08.066
10. Das S, Maji SL, Guria M, Jana RN. Hall effects on unsteady MHD flow between two disks with non-coincident parallel axes of rotation. Int J Appl Mech Eng. (2010) 15:5–18. doi: 10.1007/s40010-012-0049-1
11. Asgher S, Jalil M, Hussan M, Turkyimazoglu M. Lie group analysis of flow and heat transfer over a stretching rotating disk. Int J Heat Mass Transf. (2014) 69:140–6. doi: 10.1016/j.ijheatmasstransfer.2013.09.061
19. Aqsa, Malik M, Imtiaz A, Awais M. Rheology of Burgers' model with Cattaneo-Christov heat flux in the presence of heat source/sink and magnetic field. Scientia Iranica. (2019) 26:323–30. doi: 10.24200/sci.2018.20165
21. Mekheimer KS, Elmaboud YA. Simultaneous effects of variable viscosity and thermal conductivity on peristaltic flow in a vertical asymmetric channel. Cand J Physc. (2014) 92:1541–55. doi: 10.1139/cjp-2013-0465
24. Hayat T, Khan MI, Farooq M, Alsaedi A, Waqas M, Yasmeen T. Impact of Cattaneo-Christov heat flux model in flow of variable thermal conductivity fluid over a variable thicked surface. Int J Heat Mass Transf. (2016) 66:702–10. doi: 10.1016/j.ijheatmasstransfer.2016.04.016
29. Chen D, Mioshi H, Akai T, Yazawa Y. Colorless transparent fluorescence material: sintered porous glass containing rare-earth and transition-metal ions. Appl Phys Lett. (2005) 86:231–8. doi: 10.1063/1.1946897
30. Muhammad T, Alsaedi A, Shahzad SA, Hayat T. A revised model for Darcy-Forchheimer flow of Maxwell nanofluid subject to convective boundary condition. Chin J Phys. (2017) 55:963–76. doi: 10.1016/j.cjph.2017.03.006
31. Seddeek MA. Influence of viscous dissipation and thermophoresis on Darcy-Forchheimer mixed convection in a fluid saturated porous media. J Colloids Interface Sci. (2006) 293:137–42. doi: 10.1016/j.jcis.2005.06.039
33. Aziz T, Mabood FM, Shahzad A. Travelling wave solutions for the unsteady flow of a third grade fluid induced due to impulsive motion of flat porous plate embedded in a porous medium. J Mech. (2017) 30:527–35. doi: 10.1017/jmech.2014.17
38. Yasin MHM, Ishak A, Pop I. MHD heat and mass transfer flow over a permeable stretching/shrinking sheet with radiation effect. J Magnet Mag Mater. (2016) 407:235–40. doi: 10.1016/j.jmmm.2016.01.087
40. Khalid A, Khan I, Khan A, Shafie S. Unsteady MHD free convection flow of Casson fluid past over an oscillating vertical plate embedded in a porous medium. Eng Sci Tech Int J. (2015) 18:307–19. doi: 10.1016/j.jestch.2014.12.006
41. Abro KA, Khan I, Nisar KS. Novel technique of Atangana and Baleanu for heat dissipation in transmission line of electrical circuit. Chaos Solit Fractals. (2017) 129:40–5. doi: 10.1016/j.chaos.2019.08.001
42. Saqib M, Ali F, Khan I, Sheikh NA, Aftab S, Jan SAA, et al. Exact solutions for free convection flow of generalized Jeffrey fluid: a Caputo-Fabrizio fractional model. Alex Eng J. (2017) 57:1849–58. doi: 10.1016/j.aej.2017.03.017
Keywords: Cattaneo Christov heat flux model, permeable medium, fluid flow with coaxially rotated disks, implicit finite difference scheme, coaxially rotated disks, viscous fluid
Citation: Bilal S, Tassaddiq A, Majeed AH, Nisar KS, Ali F and Malik MY (2020) Computational and Physical Examination About the Aspects of Fluid Flow Between Two Coaxially Rotated Disks by Capitalizing Non-fourier Heat Flux Theory: Finite Difference Approach. Front. Phys. 7:209. doi: 10.3389/fphy.2019.00209
Received: 21 September 2019; Accepted: 21 November 2019;
Published: 09 January 2020.
Edited by:Muhammad Mubashir Bhatti, Shanghai University, China
Reviewed by:Anwar Shahid, Nanjing University of Aeronautics and Astronautics, China
Ilyas Khan, Ton Duc Thang University, Vietnam
Kh S. Mekheimer, Al-Azhar University, Egypt
Copyright © 2020 Bilal, Tassaddiq, Majeed, Nisar, Ali and Malik. 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: Sardar Bilal, email@example.com