[277-286]-Kunt Atallk.Fm

Total Page:16

File Type:pdf, Size:1020Kb

[277-286]-Kunt Atallk.Fm Korea-Australia Rheology Journal, Vol.24, No.4, pp.277-286 (2012) www.springer.com/13367 DOI: 10.1007/s13367-012-0034-x Isothermal and non-isothermal viscoelastic flow of PTT fluid in lid-driven polar cavity Hatice Mercan and Kunt Atal k* o Mechanical Engineering Department, Bog aziçi University, 34342 Bebek, Istanbul, Turkey (Received January 24, 2012; final revision received October 16, 2012; accepted October 31, 2012) The isothermal and non-isothermal viscoelastic flow of Phan-Thien-Tanner (PTT) fluids is considered in lid- driven polar cavity geometry, using a numerical solution method with parameter continuation technique. Thermoelastic effects, in terms of elastic/elongational effects and viscous dissipation, are demonstrated by the changes in vortical structure, temperature/stress distributions and heat transfer characteristics in the curved cavity. Central vortex/maximum temperature location shifts are observed under elastic and elon- gational (strain hardening and strain softening/shear thinning) effects for isothermal and non-isothermal conditions. The growth in size and strength of a secondary vortex is denoted in the downstream stationary corner of the cavity for the viscoelastic fluid under strain hardening, which also introduces an increase in stress gradients. Viscous heating is observed with elongational effects near the central vortex in the cavity. Stress components and their gradients decrease under viscous dissipation. The changes in temperature field and heat transfer properties in the cavity are revealed. Keywords: Lid-driven polar cavity, viscoelastic, thermoelastic, non-isothermal, viscous dissipation, Phan- Thien-Tanner fluid 1. Introduction tally and numerically the increase of temperature due to viscous dissipation in Newtonian and non-Newtonian In this study we investigate non-isothermal effects on power-law (shear thinning) fluids for the flow in a scraped two-dimensional viscoelastic flow in a polar cavity, under surface heat exchanger. They noticed an increase in vis- low inertia. The polar cavity geometry can also be con- cous heating as the rotor speed increases. However, they sidered as a typical cross section of a single screw reported that the effect of viscous dissipation was not extruder, and is suitable to illustrate the curved geometry observed for the power-law fluid. They attributed this effect on an internal flow field. This geometry (with cavity behavior to the shear thinning character of the fluid which angle of 1 radian and gap ratio of 1/2) has been described considerably reduces the apparent viscosity and therefore and used in a study by Fuchs and Tillmark (1985), where the viscous dissipation in the high shear rate zones near experimental and numerical simulation results have been the blade. obtained for Newtonian flow at moderate Reynolds num- Peters and Baaijens (1997) considered the modeling of bers. Wu and co-workers (2007) solved numerically the non-isothermal viscoelastic flows using different consti- polar cavity flow of a Newtonian fluid up to Reynolds tutive equations and analyzed the partition between dis- number Re = 1000 to validate their Lagrange interpolation sipated and elastically stored energy in terms of entropy method. Also, Kim (1989), Lei and co-workers (2000), and energy elasticity. Yeçs ilata (2002) formulated the prob- and Darbandi and Vikilipour (2008) used this geometry as lem for the non-isothermal flow of an Oldroyd-B fluid a test and example case for their numerical schemes to between two rotating parallel plates and investigated how develop solvers for Navier-Stokes equations. to predict the material properties of viscoelastic fluids The effects of viscous dissipation on Newtonian and under viscous dissipation effects. non-Newtonian flow fields have been studied by many Purely elastic and thermoelastic instabilities have been researchers, for curved geometries. Al-Mubaiyedh and co- revealed in flows with curved streamlines like Taylor- workers (2002) investigated the influence of viscous heat- Couette, parallel plates, cone-and-plate and Dean flow ing on the stability of Taylor-Couette flow for Newtonian geometries (Muller, 2008). In addition to purely elastic fluids with thermally sensitive viscosity. Based on the lin- instabilities, viscous heating leads to temperature gradi- ear stability analysis, they showed that viscous heating ents in the flow field which in turn produce gradients in leads to significant destabilization in Taylor-Couette flow. fluid viscosity and elasticity. These gradients are observed Yataghane and co-workers (2009) analyzed experimen- to lead to new thermoelastic instability modes. Thomas and co-workers (2004) investigated the effect of base flow temperature and inertia on the stability of Taylor-Couette *Corresponding author: [email protected] flow and noted that inertia effects may become important © 2012 The Korean Society of Rheology and Springer 277 Hatice Mercan and Kunt Atali k for sufficiently large temperature ranges. Coelho and co-workers (2003) presented a theoretical study of viscoelastic thermal entry flow (Graetz problem) in channels and tubes for the PTT model. For the imposed temperature and fluid heating case, they observed that vis- cous heating due to dissipation tends to increase in regions where the velocity gradients are steeper, especially near the wall. They also mentioned that as elongational/elastic effects increase the overall heat transfer tends to increase in fully developed flow regime. Pinho and Coelho (2006) investigated theoretically a simplified viscoelastic PTT flow in annuli under viscous dissipation and found that for imposed wall temperature case, shear thinning and elas- ticity lead to lower heat transfer when viscous dissipation is weak and to higher heat transfer when viscous dissi- pation is strong. In this study we investigate numerically the isothermal and non-isothermal viscoelastic flow of a PTT fluid in a lid-driven polar cavity. Elastic, elongational (strain hard- ening and strain softening/shear thinning) and viscous dis- Fig. 1. Cavity geometry. sipation effects are considered in the analysis. As a general model the Williams-Lendel-Ferry equation is used application, (Karwe and Jaluria, 1990). to describe the temperature dependence of the material The mass and momentum conservation equations in parameters such as the viscosity and the relaxation time. dimensionless vector form for unsteady, non-isothermal, Thermoelastic effects on the flow and thermal fields are incompressible viscoelastic flow can be written as, shown and compared with respect to the Newtonian case. ∇⋅ V = 0 (1) The results are presented and discussed in terms of the flow structure, the temperature and stress distributions, as DV well as the heat transfer properties in the cavity. Re-------- = –∇p+∇⋅ []∇21()–ωr D + ⋅ T (2) The isothermal and non-isothermal modeling of vis- Dt coelastic flow for PTT fluids is given in Section 2 in terms where V is the velocity vector, D is the rate of deformation of non-dimensional parameters and equations. The details tensor, T is the viscoelastic extra-stress tensor, p is the of the numerical solution method together with numerical pressure, T is the temperature, D/Dt denotes the material tests and comparisons to available data in the literature are derivative. The Reynolds number is defined as presented in Section 3. The results are shown and dis- Re = ρV0ri ⁄µ0 , where ρ is the fluid density, V0 is the char- cussed in Section 4, before concluding with final remarks acteristic lid velocity and µ0 is the total viscosity calcu- in Section 5. lated at reference temperature. The pressure and viscoelastic stress are scaled with µ0V0 ⁄ri . The non-dimen- ˜ 2. Mathematical Modeling sional temperature is defined as TT= ()–Ti ⁄()To –Ti where T˜ is the dimensional temperature. The retardation param- The two-dimensional flow of a polymer solution is con- eter ωr is defined as the ratio of zero shear rate polymer sidered in a polar cavity with cavity angle θp = π ⁄2 and viscosity µp at reference temperature, to µ0 . The vis- inner radius ri = 1 , outer radius ro = 1.75 (Fig. 1). The coelastic stress (T) is defined using linear PTT constitutive outer and side walls are fixed, while the inner lid moves relation. In dimensionless form the non-isothermal PTT in the transverse θ direction with a predefined velocity model reads, distribution. The fluid is initially at rest and no-slip bound- ∇ ε ()⎛⎞() ω () ary conditions are specified at the walls. T+2We T ⎝⎠T+ ------------- tr T T = r T D (3) ωr()T The outer wall is kept at a constant high temperature To , while the moving inner lid is kept at a constant low tem- where tr denotes the trace. In the above constitutive rela- perature Ti (To >Ti ). The stationary side walls are insu- tion, the temperature dependency of the connector force lated. Since this geometry is a better approximation in related to polymeric fluid structure in kinetic theory, is two-dimensional space for single screw extruders com- neglected following the final formulation in Peters and pared to the rectangular cavity geometry, the thermal Baaijens (1997). The Weissenberg number is defined as boundary conditions are determined according to this We = λVo ⁄ri , where λ is the relaxation time of the fluid 278 Korea-Australia Rheology J., Vol. 24, No. 4 (2012) Isothermal and non-isothermal viscoelastic flow of PTT fluid in lid-driven polar cavity and ε is the elongational parameter of the PTT model. This lowing the assumption that the temperatures are close to parameter characterizes the strain hardening and strain glass transition temperature, for strong thermorheological softening/shear ∇thinning feature of the polymer additive. coupling (Peters and Baaijens, 1997). In the above for- The operator () denotes the upper convected derivative mulation (Eq.6), the temperature dependence of the steady as follows, state recoverable compliance (the ratio of polymer relax- ∇ ation time to polymer viscosity) which is found to be weak DT T T = -------- –L ⋅TTL– ⋅ (4) in temperature ranges close to glass transition temperature Dt and often neglected in the literature (Ngai and Plazek, where L is the velocity gradient tensor.
Recommended publications
  • Soret and Dufour Effects on MHD Boundary Layer Flow of Non
    Pramana – J. Phys. (2020) 94:108 © Indian Academy of Sciences https://doi.org/10.1007/s12043-020-01984-z Soret and Dufour effects on MHD boundary layer flow of non-Newtonian Carreau fluid with mixed convective heat and mass transfer over a moving vertical plate ANIL KUMAR GAUTAM1, AJEET KUMAR VERMA1, KRISHNENDU BHATTACHARYYA1 ,∗ and ASTICK BANERJEE2 1Department of Mathematics, Institute of Science, Banaras Hindu University, Varanasi 221 005, India 2Department of Mathematics, Sidho-Kanho-Birsha University, Purulia 723 104, India ∗Corresponding author. E-mail: [email protected], [email protected] MS received 1 August 2019; revised 24 April 2020; accepted 3 May 2020 Abstract. In this analysis, the mixed convection boundary layer MHD flow of non-Newtonian Carreau fluid subjected to Soret and Dufour effects over a moving vertical plate is studied. The governing flow equations are converted into a set of non-linear ordinary differential equations using suitable transformations. For numerical computations, bvp4c in MATLAB package is used to solve the resulting equations. Impacts of various involved parameters, such as Weissenberg number, power-law index, magnetic parameter, thermal buoyancy parameter, solutal buoyancy parameter, thermal radiation, Dufour number, Soret number and reaction rate parameter, on velocity, temperature and concentration are shown through figures. Also, the local skin-friction coefficient, local Nusselt number and local Sherwood number are calculated and shown graphically and in tabular form for different parameters. Some important facts are revealed during the investigation. The temperature and concentration show decreasing trends with increasing values of power-law index, whereas velocity shows reverse trend and these trends are more prominent for larger values of Weissenberg number.
    [Show full text]
  • Thermal Radiations and Mass Transfer Analysis of the Three-Dimensional Magnetite Carreau Fluid Flow Past a Horizontal Surface of Paraboloid of Revolution
    processes Article Thermal Radiations and Mass Transfer Analysis of the Three-Dimensional Magnetite Carreau Fluid Flow Past a Horizontal Surface of Paraboloid of Revolution T. Abdeljawad 1,2,3 , Asad Ullah 4 , Hussam Alrabaiah 5,6, Ikramullah 7, Muhammad Ayaz 8, Waris Khan 9 , Ilyas Khan 10,∗ and Hidayat Ullah Khan 11 1 Department of Mathematics and General Sciences, Prince Sultan University, Riyadh 11 586, Saudi Arabia; [email protected] 2 Department of Medical Research, China Medical University, Taichung 40402, Taiwan 3 Department of Computer Science and Information Engineering, Asia University, Taichung, Taiwan 4 Institute of Numerical Sciences, Kohat University of Science & Technology, Kohat 26000, Khyber Pakhtunkhwa, Pakistan; [email protected] 5 College of Engineering, Al Ain University, Al Ain 64141, UAE; [email protected] 6 Department of Mathematics, Tafila Technical University, Tafila 66110, Jordan 7 Department of Physics, Kohat University of Science & Technology, Kohat 26000, Khyber Pakhtunkhwa, Pakistan; [email protected] 8 Department of Mathematics, Abdul Wali Khan University, Mardan 23200, Khyber Pakhtunkhwa, Pakistan; [email protected] 9 Department of Mathematics, Islamia College University Peshawar, Peshawar 25000, Khyber Pakhtunkhwa, Pakistan; [email protected] 10 Faculty of Mathematics and Statistics, Ton Duc Thang University, Ho Chi Minh City 72915, Vietnam 11 Department of Economics, Abbottabad University of Science and Technology (AUST), Abbottabad Havelian 22500, Khyber Pakhtunkhwa, Pakistan; [email protected] * Correspondence: [email protected] Received: 18 April 2020; Accepted: 22 May 2020; Published: 1 June 2020 Abstract: The dynamics of the 3-dimensional flow of magnetized Carreau fluid past a paraboloid surface of revolution is studied through thermal radiation and mass transfer analysis.
    [Show full text]
  • Research Article a Combined Convection Carreau–Yasuda Nanofluid Model Over a Convective Heated Surface Near a Stagnation Point: a Numerical Study
    Hindawi Mathematical Problems in Engineering Volume 2021, Article ID 6665743, 14 pages https://doi.org/10.1155/2021/6665743 Research Article A Combined Convection Carreau–Yasuda Nanofluid Model over a Convective Heated Surface near a Stagnation Point: A Numerical Study Azad Hussain,1 Aysha Rehman ,1 Sohail Nadeem,2 M. Y. Malik ,3 Alibek Issakhov,4,5 Lubna Sarwar,1 and Shafiq Hussain6 1Department of Mathematics, University of Gujrat, Gujrat 50700, Pakistan 2Department of Mathematics, Quaid-I-Azam University, Islamabad 44000, Pakistan 3Department of Mathematics, College of Sciences, King Khalid University, Abha 61413, Saudi Arabia 4Department of Mathematical and Computer Modeling, Al-Farabi Kazakh National University, Almaty, Kazakhstan 5Department of Mathematical and Computer Modeling, Kazakh British-Technical University, Almaty, Kazakhstan 6Department of Computer Science, University of Sahiwal, Sahiwal, Pakistan Correspondence should be addressed to Aysha Rehman; [email protected] Received 18 November 2020; Revised 1 March 2021; Accepted 22 March 2021; Published 5 April 2021 Academic Editor: Adrian Neagu Copyright © 2021 Azad Hussain et al. 'is is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. 'e focus of this manuscript is on two-dimensional mixed convection non-Newtonian nanofluid flow near stagnation point over a stretched surface with convectively heated boundary conditions.
    [Show full text]
  • On Dimensionless Numbers
    chemical engineering research and design 8 6 (2008) 835–868 Contents lists available at ScienceDirect Chemical Engineering Research and Design journal homepage: www.elsevier.com/locate/cherd Review On dimensionless numbers M.C. Ruzicka ∗ Department of Multiphase Reactors, Institute of Chemical Process Fundamentals, Czech Academy of Sciences, Rozvojova 135, 16502 Prague, Czech Republic This contribution is dedicated to Kamil Admiral´ Wichterle, a professor of chemical engineering, who admitted to feel a bit lost in the jungle of the dimensionless numbers, in our seminar at “Za Plıhalovic´ ohradou” abstract The goal is to provide a little review on dimensionless numbers, commonly encountered in chemical engineering. Both their sources are considered: dimensional analysis and scaling of governing equations with boundary con- ditions. The numbers produced by scaling of equation are presented for transport of momentum, heat and mass. Momentum transport is considered in both single-phase and multi-phase flows. The numbers obtained are assigned the physical meaning, and their mutual relations are highlighted. Certain drawbacks of building correlations based on dimensionless numbers are pointed out. © 2008 The Institution of Chemical Engineers. Published by Elsevier B.V. All rights reserved. Keywords: Dimensionless numbers; Dimensional analysis; Scaling of equations; Scaling of boundary conditions; Single-phase flow; Multi-phase flow; Correlations Contents 1. Introduction .................................................................................................................
    [Show full text]
  • Free Convective Heat Transfer Over Truncated Cone with Convective Heating: Numerical Analysis
    International Journal of Advanced Scientific Research and Management, Volume 4 Issue 3, Mar 2019 www.ijasrm.com ISSN 2455-6378 Free Convective Heat Transfer over Truncated Cone with Convective Heating: Numerical Analysis G. Raju Department of Mathematics, IIIT Ongole Campus, RGUKT-516330, Andhra Pradesh, India. Abstract cone embedded in a porous medium. Elbashbeshy et This paper is focused on the study of heat transfer al. [6] analyzed radiation effects non-Newtonian characteristics of MHD non-Newtonian Tangent fluid flow past a vertical truncated cone. Subba Rao Hyperbolic boundary layer flow over a truncated et al. [7] discussed the Newtonian heating on cone with Biot number effect. The governing partial hydrodynamic fluid flow over a circular cylinder. differential equations of the flow field are converted Reddy and Pradeepa [8] studied viscous disspation to a system of non-linear coupled non-similarity and sorret effects on natural convection flow over a ordinary differential equations. Employing Shooting truncated cone in the presence of biot number. technique followed by finite centered difference Several authors [9-11] investigated heat transfer of technique (Keller Box Method), the system is solved the non-Newtonian fluid of various physical numerically. The numerical code is validated with problems. previous studies. The results of flow and heat analysis on velocity and temperature profiles are given diagrammatically. The numerical results point The present work, we examine theoretically and out that, increasing MHD parameter (M) serves to computationally the steady-state transport decelerate the flow, but increased temperature values phenomena in magneto-hydrodynamic non- and increase in the Biot number (Bi) is observed to Newtonian fluid flow from a truncated cone with enhance velocity and temperature.
    [Show full text]
  • The British Society of Rheology
    THE BRITISH SOCIETY OF RHEOLOGY To cite: R J Poole (2012), The Deborah and Weissenberg numbers. The British Society of Rheology, Rheology Bulletin . 53(2) pp 32-39 1 The Deborah and Weissenberg numbers Rob Poole School of Engineering University of Liverpool, United Kingdom Since their introduction to the literature nearly fifty years ago the Deborah and Weissenberg numbers have proved invaluable for rheologists in quantifying viscoelastic effects. Despite their different origins, and the fact that they quantify different effects, they are frequently used as synonyms. This situation is perhaps partly explainable given that in a wide range of flows of complex fluids characterised by a relaxation time ( λ) the use of simple single characteristic length ( L) and velocity scales ( U) does indeed result in identical definitions! In this short article we return to the original papers to highlight the different effects that the Deborah and Weissenberg numbers measure. The Deborah number As every good rheologist knows the Deborah number owes its name to the Prophetess Deborah who, in the Book of Judges, proclaimed “The mountains flowed before the lord” A. The definition is due to Reiner and we can do no better than to quote his original paper B “Deborah knew two things. First, that the mountains flow, as everything flows. But, secondly, that they flowed before the Lord, and not before man, for the simple reason that man in his short lifetime cannot see them flowing, while the time of observation of God is infinite . We may therefore well define a nondimensional number the Deborah number D = time of relaxation/time of observation.” AAs Huilgol notes “Reiner quoted the Prophetess Deborah as saying “The mountains flowed before the Lord” while the King James version says “The mountains melted from before the Lord”, and the New English Bible has “Mountains shook in fear before the Lord,” (See Judges 5:5) [R R Huigol (1975) “On the concept of the Deborah number”.
    [Show full text]
  • Conjugate Heat Transfer in the Unbounded Flow of a Viscoelastic Fluid Past a Sphere
    Conjugate heat transfer in the unbounded flow of a viscoelastic fluid past a sphere F. Pimenta, M.A. Alves* CEFT, Departamento de Engenharia Química, Faculdade de Engenharia da Universidade do Porto, Rua Dr. Roberto Frias, 4200-465 Porto, Portugal Abstract This work addresses the conjugate heat transfer of a simplified PTT fluid flowing past an unbounded sphere in the Stokes regime (Re = 0.01). The problem is numerically solved with the finite-volume method assuming axisymmetry, absence of natural convection and constant physical properties. The sphere generates heat at a constant and uniform rate, and the analysis is conducted for a range of Deborah (0 ≤ De ≤ 100), Prandtl (100 ≤ Pr ≤ 105) and Brinkman (0 ≤ Br ≤ 100) numbers, in the presence or absence of thermal contact resistance at the solid-fluid interface and for different conductivity ratios (0.1 ≤ κ ≤ 10). The drag coefficient shows a monotonic decrease with De, whereas the normalized stresses on the sphere surface and in the wake first increase and then decrease with De. A negative wake was observed for the two solvent viscosity ratios tested (β = 0.1 and 0.5), being more intense for the more elastic fluid. In the absence of viscous dissipation, the average Nusselt number starts to decrease with De after an initial increase. Heat transfer enhancement relative to an equivalent Newtonian fluid was observed for the whole range of conditions tested. The dimensionless temperature of the sphere decreases and becomes more homogeneous when its thermal conductivity increases in relation to the conductivity of the fluid, although small changes are observed in the Nusselt number.
    [Show full text]
  • Impact of Double-Diffusive Convection and Motile Gyrotactic
    Open Physics 2020; 18: 74–88 Review Article Tanveer Sajid*, Muhammad Sagheer, Shafqat Hussain, and Faisal Shahzad Impact of double-diffusive convection and motile gyrotactic microorganisms on magnetohydrodynamics bioconvection tangent hyperbolic nanofluid https://doi.org/10.1515/phys-2020-0009 to be explored widely because of its enormous applications in received July 07, 2019; accepted February 10, 2020 the field of pharmaceutical industry, purification of cultures, fl Abstract: The double-diffusive tangent hyperbolic nanofluid micro uidic devices, mass transport enhancement and containing motile gyrotactic microorganisms and magneto- mixing, microbial enhanced oil recovery and enzyme hydrodynamics past a stretching sheet is examined. By biosensors. Bioconvection systems could be categorized ff adopting the scaling group of transformation, the governing based on the directional motion of di erent species of equations of motion are transformed into a system of microorganisms. In particular, gyrotactic microorganisms are nonlinear ordinary differential equations. The Keller box the ones whose swimming direction is dependent on a [ ] scheme, a finite difference method, has been employed for the balance between gravitational and viscous torques 4,5 . [ ] solution of the nonlinear ordinary differential equations. The Oyelakin et al. 6 pondered the impact of bioconvection and fl behaviour of the working fluid against various parameters of motile gyrotactic microorganisms on the Casson nano uid physical nature has been analyzed through graphs and tables. past a stretching sheet and observed that the microorganism fi The behaviour of different physical quantities of interest pro le decreases as a result of an increment in the Peclet [ ] ff such as heat transfer rate, density of the motile gyrotactic number.
    [Show full text]
  • Numerical Study of Particle Suspensions in Newtonian and Non-Newtonianfluids
    Numerical study of particle suspensions in Newtonian and non-Newtonianfluids by Dhiya Abdulhussain Alghalibi December 2019 Technical Reports Royal Institute of Technology Department of Mechanics SE-100 44 Stockholm, Sweden Akademisk avhandling som med tillst˚andav Kungliga Tekniska H¨ogskolan i Stockholm framl¨agges till offentlig granskning f¨or avl¨aggande av teknologie doctorsexamen Fredag den 6 December 2019 kl 10:15 i Angdomen˚ (Rumsnr: 5209), Osquars backe 31, KTHB, v˚aningsplan 2, KTH Campus, Stockholm. TRITA-MEK 2019;55 ISSN 0384-467X ISRN KTH/MEK/TR-19/55-SE ISBN 978-91-7873-385-9 Cover: Instantaneous snapshot of the particle arrangement for a laminar shear thinning Couetteflow, ˙γ= 0.1, Rep = 0.5 andΦ = 40%. Particles are colored according to their velocity in the streamwise direction, where the red and blue colors denote the negative and positive directions, respectively. Particle diameters are shown at their actual size. c Dhiya Abdulhussain Alghalibi 2019 � Universitetsservice US–AB, Stockholm 2019 Dedicated to my lovely family ... Numerical study of particle suspensions in Newtonian and non-Newtonianfluids Dhiya Abdulhussain Alghalibi Linn´eFLOW Centre, KTH Mechanics, Royal Institute of Technology SE-100 44 Stockholm, Sweden. Abstract Solid or deformable particles suspended in a viscousfluid are of scientific and technological interest in a broad range of applications. Pyroclasticflows from volcanoes, sedimentationflows in river bed, food industries, oil-well drilling, as well as bloodflow in the human body and the motion of suspended micro- organisms in water (like plankton) are among the possible examples. Often, in these particulateflows, the carrierfluid might exhibit an inelastic or a visco-elastic non-Newtonian behavior.
    [Show full text]
  • Effects of Viscoelasticity in the High Reynolds Number Cylinder Wake
    J. Fluid Mech. (2012), vol. 693, pp. 297–318. c Cambridge University Press 2012 297 doi:10.1017/jfm.2011.531 Effects of viscoelasticity in the high Reynolds number cylinder wake David Richter, Gianluca Iaccarino and Eric S. G. Shaqfeh† Department of Mechanical Engineering, Stanford University, CA 94305, USA (Received 26 May 2011; revised 17 October 2011; accepted 29 November 2011; first published online 16 January 2012) At Re 3900, Newtonian flow past a circular cylinder exhibits a wake and detached shearD layers which have transitioned to turbulence. It is the goal of the present study to investigate the effects which viscoelasticity has on this state and to identify the mechanisms responsible for wake stabilization. It is found through numerical simulations (employing the FENE-P rheological model) that viscoelasticity greatly reduces the amount of turbulence in the wake, reverting it back to a state which qualitatively appears similar to the Newtonian mode B instability which occurs at lower Re. By focusing on the separated shear layers, it is found that viscoelasticity suppresses the formation of the Kelvin–Helmholtz instability which dominates for Newtonian flows, consistent with previous studies of viscoelastic free shear layers. Through this shear layer stabilization, the viscoelastic far wake is then subject to the same instability mechanisms which dominate for Newtonian flows, but at far lower Reynolds numbers. Key words: instability, viscoelasticity, wakes 1. Introduction In Newtonian flow past a circular cylinder, transition to turbulence occurs in several distinct stages as the Reynolds number is increased, and this path of transition is outlined nicely by Williamson (1996b).
    [Show full text]
  • Mathematical Study of Laminar Boundary Layer Flow and Heat Transfer of Tangenthyperbolic Fluid Pasta Vertical Porous Plate with Biot Number Effects
    Journal of Applied Fluid Mechanics, Vol. 9, No. 3, pp. 1297-1307, 2016. Available online at www.jafmonline.net, ISSN 1735-3572, EISSN 1735-3645. DOI: 10.18869/acadpub.jafm.68.228.24217 Mathematical Study of Laminar Boundary Layer Flow and Heat Transfer of Tangenthyperbolic Fluid Pasta Vertical Porous Plate with Biot Number Effects V. Ramachandra Prasad1, S. Abdul Gaffar2 and O. Anwar Bég3 1 Department of Mathematics, Madanapalle Institute of Technology and Science, Madanapalle-517325, India 2Department of Mathematics, Salalah College of Technology, Salalah, Oman 3 Gort Engovation Research (Aerospace), 15 Southmere Avenue, Great Horton, Bradford, BD7 3NU, West Yorkshire, UK. †Corresponding Author Email: [email protected] (Received December 12, 2014; accepted July 29, 2015) ABSTRACT In this article, we investigate the nonlinear steady boundary layer flow and heat transfer of an incompressible Tangent Hyperbolicnon-Newtonian fluid from a vertical porous plate. The transformed conservation equations are solved numerically subject to physically appropriate boundary conditions using a second-order accurate implicit finite-difference Keller Box technique. The numerical code is validated with previous studies. The influence of a number of emerging non-dimensional parameters, namely the Weissenberg number (We), the power law index (n), Prandtl number (Pr), Biot number (), and dimensionless local suction parameter()on velocity and temperature evolution in the boundary layer regime are examined in detail. Furthermore the effects of these parameters on surface heat transfer rate and local skin friction are also investigated. Validation with earlier Newtonian studies is presented and excellent correlation achieved. It is found that velocity, Skin friction and Nusselt number (heat transfer rate) are reduced with increasing Weissenberg number (We), whereas, temperature is enhanced.
    [Show full text]
  • Numerical Study of Entropy Generation and Melting Heat Transfer on MHD Generalised Non-Newtonian fluid (GNF): Application to Optimal Energy
    Pramana – J. Phys. (2018) 90:64 © Indian Academy of Sciences https://doi.org/10.1007/s12043-018-1557-6 Numerical study of entropy generation and melting heat transfer on MHD generalised non-Newtonian fluid (GNF): Application to optimal energy Z IQBAL, ZAFFAR MEHMOOD and BILAL AHMAD∗ Department of Mathematics, HITEC University, Taxila, Pakistan ∗Corresponding author. E-mail: [email protected] MS received 6 September 2017; revised 17 November 2017; accepted 23 November 2017; published online 16 April 2018 Abstract. This paper concerns an application to optimal energy by incorporating thermal equilibrium on MHD- generalised non-Newtonian fluid model with melting heat effect. Highly nonlinear system of partial differential equations is simplified to a nonlinear system using boundary layer approach and similarity transformations. Numerical solutions of velocity and temperature profile are obtained by using shooting method. The contribution of entropy generation is appraised on thermal and fluid velocities. Physical features of relevant parameters have been discussed by plotting graphs and tables. Some noteworthy findings are: Prandtl number, power law index and Weissenberg number contribute in lowering mass boundary layer thickness and entropy effect and enlarging thermal boundary layer thickness. However, an increasing mass boundary layer effect is only due to melting heat parameter. Moreover, thermal boundary layers have same trend for all parameters, i.e., temperature enhances with increase in values of significant parameters. Similarly, Hartman and Weissenberg numbers enhance Bejan number. Keywords. Tangent hyperbolic fluid; numerical solutions; melting heat transfer; entropy generation; optimal energy. PACS Nos 44.25.+f; 47.10.ad; 47.50.−d 1. Introduction science and engineering.
    [Show full text]