Abstract
An attempt has been made to describe the effects of geothermal viscosity with viscous dissipation on the three dimensional time dependent boundary layer flow of magnetic nanofluids due to a stretchable rotating plate in the presence of a porous medium. The modelled governing time dependent equations are transformed a from boundary value problem to an initial value problem, and thereafter solved by a fourth order Runge-Kutta method in MATLAB with a shooting technique for the initial guess. The influences of mixed temperature, depth dependent viscosity, and the rotation strength parameter on the flow field and temperature field generated on the plate surface are investigated. The derived results show direct impact in the problems of heat transfer in high speed computer disks (Herrero et al. [1]) and turbine rotor systems (Owen and Rogers [2]).
1 Introduction
In the last few years, deep interest has been shown in investigating and understanding the boundary layer flow and convective heat transfer in magnetic nanofluid flow, which is a topic of major contemporary interest in both science and engineering. Magnetic nanofluids, popularly known as Ferrofluids, are fluids in which nano-sized particles of Fe3O4, γ-Fe2O3 or CoFe2O4, having an average size of range about 3-15 nm, are dispersed in a carrier fluid (e.g. water, kerosene, ethylene glycol, toluene, biofluids, or lubricants). These particles are coated with a stabilizing dispersing agent called surfactant (antimony, oleic acid, glycerol, fatty acids, etc.) to prevent agglomeration even when a strong magnetic field gradient is applied to the fluid. For the past few decades, ferrofluids have been widely research for their numerous heat exchanger applications in engineering, industry, the physical and medical sciences, etc. The pioneer study of ordinary viscous fluid flow over a rotating disk was described by Karman [3]. After his work, many scientists focused on the study of boundary layer flow due to a rotating disk (Cochran [4], Benton [5], El-Mistikawy and Attia [6], and El-Mistikawy et al. [7]). Turkyilmazoglu [8] focused on a class of exact solutions to the steady Navier-Stokes equations for incompressible Newtonian viscous fluid flow motion due to a rotating porous disk. By considering the field dependent viscosity, Ram and Sharma [9] described the behaviour of ferrofluid flow over a rotating disk. Soundalgekar [10] studied the effect ofviscous dissipation on unsteady free convection flow past an infinite, vertical porous plate with uniform suction.
The effects of viscosity variation due to the geophysical position and temperature, taken one at a time, have been carried out by many researchers in MHD and few in FHD. However, the effects of mixed depth and temperature dependent viscosity have not been popularly analyzed in the rotating disk problem yet. Geothermal (depth and temperature dependent) viscosity has a number of fruitful applications in geophysics and the geosciences. Important applications of flow induced by a stretching boundary are found in extrusion processes in the plastic and metal industries (Altan et al. [11]). Recently, Turkyilmazoglu [12] examined the steady magnetohydrodynamic(MHD) boundary layer flow over a radially stretchable rotating disk in the presence of a transverse magnetic field together with viscous dissipation and Joule heating. Ram and Kumar [13] obtained the heat flow pattern of the boundary layer flow of a ferrofluid over a stretchable rotating disk. Nawaz et al. [14] investigated the laminar boundary layer flow of nanofluid induced by a radially stretching sheet. Turkyilmazoglu [15] examined the behaviour of stagnation point flow over a stretchable rotating disk subjected to a uniform vertical magnetic field. Viscous dissipation effects on nonlinear MHD flow in a porous medium over a stretching porous surface have been studied by Devi and Ganga [16]. Chen [17] studied the effects of viscous dissipation on radiative heat transfer in MHD flow past a stretching surface. Siddheshwar et al. [18] studied temperature dependent flow behaviour and heat transfer over an exponential stretching sheet in a Newtonian liquid. Ram and Kumar [19] described the temperature dependent viscosity on steady revolving axi-symmetric laminar boundary layer flow of an incompressible ferrofluid over a stationary disk. Maleque [20] investigated the unsteady flow of an electrically conducting fluid over a rotating disk by considering the viscosity of the fluids dependent on depth and temperature together. Al-Hadhrami et al. [21] proposed a model for viscous dissipation in a porous medium which is probably adequate for most of the practical purposes. Attia [22] studied the effect of a porous medium and temperature dependent viscosity on unsteady flow and heat transfer for a viscous laminar incompressible fluid due to an impulsively-started rotating, infinite disc. Mukhopadhyay [23] presented similarity solutions for unsteady flow and heat transfer over a stretching sheet. Turkyilmazoglu [24] obtained an exact numerical solution of steady fluid flow and heat transfer induced between two stretchable co-axial disks rotating in same or opposite direction. The case of an unsteady viscous fluid flow past a horizontal stretching sheet through a porous medium in the presence of a viscous dissipation effect and a heat source was investigated by Sharma [25]. The modelled rotating disk problem of our manuscript has a number of real world engineering and industrial applications. Disk-shaped bodies are often encountered in many real world engineering applications, and heat transfer problems of free convection boundary layer flow over a rotating disk, which occur in rotating heat exchangers, rotating disk reactors for bio-fuel production, and gas or marine turbines, are extensively used by the energy, chemical, and automobile industries [26]. Many application areas, such as rotating machinery, viscometry, computer storage devices, and crystal growth processes, require the study of rotating flows [27]. In this problem, the effects of geothermal viscosity on the time dependent boundary layer flow of an electrically non-conducting, magnetic nanofluid on a stretchable rotating plate in the presence of a porous medium and viscous dissipation have been investigated. Here, the viscosity of the fluid is considered as depth and temperature dependent. The plate is maintained at a temperature Tw and subjected to a magnetic field H (Hr, Ηφ, 0). After transforming the governing equations into non-dimensionalized ordinary differential equations (ODEs) using the Von-Karman transformations, the resultant modelled set of equations is solved in Matlab to obtain the flow pattern. To study the effect of various entities of physical interest such as viscosity variation parameters, rotation parameter, porosity etc., on a hydrocarbon based magnetic nanofluid, EFH1, having the properties: fluid density ρ = 1.21x103kg/ m3, dynamic viscosity μ = 0.006kg/(ms), Prandtl number Pr = 101.2 at an average temperature of 250C, a magnetic particle concentration of 7.9%, is used. The nano particles contain a mixture of magnetite (Fe3O4 about 80%) with maghemite (Fe2O3 about 20%).
2 Mathematical Formulation of the Problem
Cylindrical polar coordinates (r,φ,z) are used to represent the problem mathematically, with an electrically nonconducting infinite rotating plate fixed at z = 0. The plate, subjected to a magnetic field H, is considered to revolve about z-axis with a constant angular velocity Ω at a uniform temperature Tw (at the plate) and temperature T∞. (far away from it).
The viscosity of the magnetic nanofluid, considered to be both depth and temperature dependent [28], is given as:
where µ∞ is the uniform viscosity of the fluid and α ≥0 is a constant [29].
The governing equations of unsteady FHD boundary layer flow along with the boundary conditions in component form are given as (Ram and Kumar [19]):
Equation of Continuity:
Equations of Momentum:
Energy Equation:
where ρ is the fluid density, µ is the magnetic permeability in free space, k is the thermal conductivity, and Cp is the specific heat at constant pressure. The boundary conditions for the flow are:
2.1 Solution of the Problem
The magnetic scalar potential due to the magnetic dipole m is given by
The corresponding magnetic field H, considering negligible variation along z-axis, has the components,
The magnetic field H is assumed to be sufficiently strong to produce the state of maximum magnetization of the magnetic nanofluid and the magnetization M is approximated, as considered by Neuringer [30] in his work, by
where Tc is he Curie temperature and K is the pyromagnetic coefficient.
The outward flow of the fluid particles due to the centrifugal force is balanced by the radial pressure force. So, the boundary layer approximation for equation (3) is
We introduce the following similarity transformations to non-dimensionalize the governing equations:
with ∆T = Tw − T∞ and δ(t) is a scalar factor.
Using eq. (11) in the Equation of Continuity (2), we get
Using the above transformations (8)–(12) in the system of governing equations (3)–(6), we have the transformed ODEs as:
Energy Equation:
where ε = α∆T and ε1 = αδ are the viscosity variation parameter and modified viscosity variation parameter, respectively,
where L represents the length scale of steady flow. Introducing (16) in equations 13)–(15, we get the following set of modelled equations:
The boundary conditions (7) reduce to
where
where Cfr and Cfϕ are the radial and tangential skin frictions coefficients, respectively, Nu is the Nusselt number, and Re is the local rotational Reynolds number.
3 Numerical Results and Discussion
Using the finite difference scheme in the MATLAB tool ODE45, the effects of various parameters on the flow behaviour near the plate for the considered flow equations (2)–(6) are discussed. From the present simulation (Table 1), it is observed that the skin friction coefficients increase with increasing viscosity variation parameters ε, ε1, and modified rotational parameter R. The rate of heat transfer is represented by the value of −θη(0). The present computation shows that the rate of heat transfer increases with increasing modified rotation parameter R and Prandtl number Pr, while it decreases with increasing viscosity variation parameters.
ε | ε1 | R | Pr | U″(0) | V′(0) | θ′(0) |
---|---|---|---|---|---|---|
0.0 | 0.0 | 0.1 | 101.2 | 0.349533145015 | 0.477964239598 | 0.228767975032 |
0.1 | 0.1 | 0.1 | 101.2 | 0.375095331450 | 0.511642395981 | 0.512479202287 |
0.5 | 0.5 | 0.1 | 101.2 | 0.465096593315 | 0.623501942396 | 1.798767975032 |
1.0 | 1.0 | 0.1 | 101.2 | 0.549465965933 | 0.725019423959 | 3.769876797503 |
0.5 | 0.5 | 0.1 | 101.2 | 0.465096593315 | 0.623501942396 | 1.798767975031 |
0.5 | 0.5 | 0.2 | 101.2 | 0.639596593315 | 0.870896250194 | 1.908984279977 |
0.5 | 0.5 | 0.3 | 101.2 | 0.781797996593 | 1.059601942396 | 2.253298984279 |
0.5 | 0.5 | 0.4 | 101.2 | 0.901797996593 | 1.217205960194 | 2.470789842799 |
0.5 | 0.5 | 0.1 | 78.4 | 0.465046965933 | 0.623219423959 | 1.283167975031 |
0.5 | 0.5 | 0.1 | 79.3 | 0.464696593315 | 0.623219423959 | 1.316797503213 |
0.5 | 0.5 | 0.1 | 101.2 | 0.465096593315 | 0.623501942396 | 1.798767975030 |
0.5 | 0.5 | 0.1 | 176.4 | 0.466024596593 | 0.624059670194 | 3.342973423798 |
Figures (2)–(5) exhibit the trends of the velocity and temperature profiles vs. the non-dimensional parameter η for viscosity variation parameters (ε & ε1) with modified rotational parameter R = 0.1, stretching parameter R1 = 1.0, Reynolds number Re = 40.0, FHD interaction parameter B = 2.0, permeability parameter β = 1.0, Prandtl number Pr = 101.2, and Eckert number Ec = 0.25. The parameter of viscosity variation ε depends on the temperature of the heated surface. The positive values of ε taken in the problem are for a heated surface, while for ε = ε1 = 0 the problem reduces to the case of uniform viscosity flow. An increase in the viscosity variations causes a decrease in the velocity profiles, while the dual nature of the temperature profile is noted along the critical point (0.62, 0.4393). The temperature profile changes its nature at the critical point from increasing to decreasing with increased viscosity variation. The fact is due to the increase in the skin friction which is caused by variation in the viscosity.
Figures 6–9 depict the influence of the modified rotational parameter on the boundary layer flow profiles. The trend shows that the effect of the rotation parameter is to decrease the velocity at all layers of fluid within the boundary layer. From Figures 6–8, it is observed that the velocity increases with decreasing values of the rotation parameter. Figure 9 shows that temperature profile decreases with increasing rotation of the plate. Also, the rate of reaching of the state of free steam velocity is faster in the velocity and temperature profiles as the value of modified rotation parameter increases.
To explicitly illustrate the influences of the Prandtl number Pr on the temperature profiles with a specified set of values of various physical parameters, such as ε = ε1 = 0.5; R = 0.1; R1 = 1.0; Re = 40.0; B = 2.0; β=1.0 & Ec = 0.25, Equations (17)–(19) have been solved and results are numerically presented in Figure 10. It is observed that the Prandtl number shows negligible effects on the velocity profile, however, an interesting behaviour is noticed in the temperature profile. A critical point (0.26, 0.9655) is found before which temperature the profile increases while after this point the temperature profile has a reverse trend. The numerical results reveal that an increase in Prandtl number results in an increase in the rate of heat transfer and a decrease in the thermal boundary layer thickness.
4 Conclusion
When we move from uniform viscosity to variable viscosity (temperature and depth dependent), the velocity profiles thins and approaches quickly to its steady state. The temperature profile increases with increasing viscosity variation parameter, causing thickening of the thermal boundary layer.
The temperature profile has a dual nature for viscosity variation and Prandtl number along the critical points (0.62, 0.4393) and (0.26, 0.9655) respectively.
In the present model, fast cooling of the device can be achieved by increasing the modified rotation parameter (R) and Prandtl number (Pr).
The present study considers the range of Prandtl numbers which have a direct application in controlling heat losses or in keeping an instrument cool and avoiding the damage caused due to the heat generation by the motion of its blades/shafts such as in the cases of thermal-power generating systems, high speed rotating machinery, and aerodynamic extrusion of plastic sheets.
References
[1] Herrero J., Humphrey J.A.C., Giralt F., Comparative analysis of coupled flow and heat transfer between co-rotating discs in rotating and fixed cylindrical enclosures, Heat Transfer in Gas Turbines, HTD- 300, 1994, 111-121.Search in Google Scholar
[2] Owen J.M., Rogers R.H., Flow and heat transfer in rotating disc system, Rotor-stator Systems, Research Studies, Taunton, UK and Wiley, New York, 1989, 1.Search in Google Scholar
[3] Karman T.V., Uber laminare und turbulente reibung, Z. Angew. Math. Mech., 1921, 1, 233-255.10.1002/zamm.19210010401Search in Google Scholar
[4] Cochran W.G., The flow due to a rotating disk, Proceedings of the Cambridge Philosophical Society, 1934, 30, 365-375.10.1017/S0305004100012561Search in Google Scholar
[5] Benton E.R., On the flow due to a rotating disk, Journal of Fluid Mechanics, 1966, 24 (4), 781-800.10.1017/S0022112066001009Search in Google Scholar
[6] El-Mistikawy T.M.A., Attia H.A., The rotating disk flow in the presence of strong magnetic field, Proceedings of the Third International Congress of Fluid Mechanics, Cairo, Egypt, 1990, 3, 1211-1222.Search in Google Scholar
[7] El-Mistikawy T.M.A., Attia H.A., Megahed A.A., The rotating disk flow in the presence of weak magnetic field, Proceedings of the Fourth Conference on Theoretical and Applied Mechanics, Cairo, Egypt, 1991, 69-82.Search in Google Scholar
[8] Turkyilmazoglu M., Effects of Partial Slip on the Analytic Heat and Mass Transfer for the Incompressible Viscous Fluid of a Porous Rotating Disk Flow, ASME Journal of Heat Transfer, 2011, 133, 122602 (5 Pages).10.1115/1.4004558Search in Google Scholar
[9] Ram P., Sharma K., On the revolving ferrofluid flow due to a rotating disk, International Journal of Nonlinear Science, 2012, 13 (3), 317-324.Search in Google Scholar
[10] Soundalgekar V.M., Viscous dissipation effect on unsteady free convection flow past an infinite vertical porous plate with variable suction, Int. J. Heat and Mass Transfer, 1974, 17, 85-92.10.1016/0017-9310(74)90041-6Search in Google Scholar
[11] Altan T., Oh S., Gegel H.L., Metal Forming, Fundamentals and Applications, American Society for Metals, Novelty, OH, 1983.Search in Google Scholar
[12] Turkyilmazoglu M., MHD fluid flow and heat transfer due to a stretching rotating disk, Int. J. Therm. Sci., 2012, 51, 195-201.10.1016/j.ijthermalsci.2011.08.016Search in Google Scholar
[13] Ram P., Kumar V., FHD flow with heat transfer over as tretchable rotating disk, Multidiscipline Modeling in Materials & Structures, 2013, 9 (4), 524-537.10.1108/MMMS-03-2013-0013Search in Google Scholar
[14] Nawaz M., Hayat T., Axisymmetric Stagnation-Point Flow of Nanofluid over a Stretching Surface, Adv. Appl. Math. Mech., 2014, 6 (2), 220-232.10.4208/aamm.2013.m93Search in Google Scholar
[15] Turkyilmazoglu M., Three dimensional MHD stagnation flow due to a stretchable rotating disk, International Journal of Heat and Mass Transfer, 2012, 55 (23), 6959-6965.10.1016/j.ijheatmasstransfer.2012.05.089Search in Google Scholar
[16] Devi S.P.A., Ganga B., Effects of Viscous and Joules dissipation on MHD flow, heat and mass transfer past a stretching porous surface embedded in a porous medium, Nonlinear Analysis: Modeling and Control, 2009, 14 (3), 303-314.10.15388/NA.2009.14.3.14497Search in Google Scholar
[17] Chen C.H., Combined effects of Joule heating and viscous dissipation on magneto hydrodynamic flow past a permeable, stretching surface with free convection and radiative heat transfer, ASME Journal of Heat Transfer, 2010, 132 (6), 064503 (5 pages).10.1115/1.4000946Search in Google Scholar
[18] Siddheshwar P.G., Sekhar G.N., Chethan A.S., Flow and Heat Transfer in a Newtonian Liquid with Temperature Dependent Properties over an Exponential Stretching Sheet, Journal of Applied Fluid Mechanics, 2014, 7 (2), 367-374.10.36884/jafm.7.02.20304Search in Google Scholar
[19] Ram P., Kumar V., Effect of temperature dependent viscosity on revolving axi-symmetric ferrofluid flow with heat transfer, Applied Mathematics and Mechanics, 2012, 33 (11), 1441-1452.10.1007/s10483-012-1634-8Search in Google Scholar
[20] Maleque K.A., Effects of combined temperature- and depth-dependent viscosity and Hall current on an unsteady MHD lami-nar convective flow due to a rotating disk, Chemical Engineering Communications, 2010, 197, 506-521.10.1080/00986440903288492Search in Google Scholar
[21] Al-Hadhrami A.K., Ellott L., Ingham D.B., A new model for viscous dissipation in porous media across a range of permeability values, Transp. Por. Media, 2003, 53, 117-122.10.1023/A:1023557332542Search in Google Scholar
[22] Attia H.A., Unsteady flow and heat transfer of viscous incompressible fluid with temperature-dependent viscosity due to a rotating disc in a porous medium, J. Phys. A: Math. Gen., 2006, 39, 979.10.1088/0305-4470/39/4/017Search in Google Scholar
[23] Mukhopadhyay S., Effect of thermal radiation on unsteady mixed convection flow and heat transfer over a porous stretching surface in porous medium, Int. J. Heat Mass Trans., 2009, 52 (13-14), 3261-3265.10.1016/j.ijheatmasstransfer.2008.12.029Search in Google Scholar
[24] Turkyilmazoglu M., Flow and heat simultaneously induced by two stretchable rotating disks, Physics of Fluids (1994-present), 2016, 28 (4), 043601.10.1063/1.4945651Search in Google Scholar
[25] Sharma R., Effect of viscous dissipation and heat source on unsteady boundary layer flow and heat transfer past a stretching surface embedded in a porous medium using element free Galerkin method, Applied Mathematics and Computation, 2012, 219 (3), 976-987.10.1016/j.amc.2012.07.002Search in Google Scholar
[26] Osalusi E., Jonathan S., Robert H., Thermal-diffusion and diffusion-thermo effects on combined heat and mass transfer of a steady MHD convective and slip flow due to a rotating disk with viscous dissipation and Ohmic heating, International Communications in Heat and Mass Transfer, 2008, 35 (8), 908-915.10.1016/j.icheatmasstransfer.2008.04.011Search in Google Scholar
[27] Brady J.F., Louis D., On rotating disk flow, Journal of Fluid Mechanics, 1987, 175, 363-394.10.1017/S0022112087000430Search in Google Scholar
[28] Lai F.C., Kulacki F.A., The effect of variable viscosity on convective heat and mass transfers along a vertical surface in saturated porous media, Int. J. Heat Mass Transfer, 1990, 33, 1028.10.1016/0017-9310(90)90084-8Search in Google Scholar
[29] Ling J.X., Dybbs A., Forced convection flow over a flat plate submerged in a porous medium with variable viscosity case, ASME paper 87-WA=TH-23, 1987.Search in Google Scholar
[30] Neuringer J.L., Some viscous flows of a saturated ferrofluid under the combined influence of thermal and magnetic field gradients, Int. J. Non-linear Mechanics, 1966, 1, 123-137.10.1016/0020-7462(66)90025-4Search in Google Scholar
[31] Schlichting H., Boundary Layer Theory, McGraw Hill, New York, 1968.Search in Google Scholar
© 2016 P. Ram et al.
This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 3.0 License.