DOCSLIB.ORG

Novel Aspects of Soret and Dufour in Entropy Generation Minimization for Williamson Uid Ow

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Novel Aspects of Soret and Dufour in Entropy Generation Minimization for Williamson Uid Ow Scientia Iranica B (2020) 27(5), 2451{2464 Sharif University of Technology Scientia Iranica Transactions B: Mechanical Engineering http://scientiairanica.sharif.edu Novel aspects of Soret and Dufour in entropy generation minimization for Williamson uid ow T. Hayata,b, F. Masooda, S. Qayyuma,*, M. Ijaz Khana, and A. Alsaedib a. Department of Mathematics, Quaid-I-Azam University 45320, Islamabad 44000, Pakistan. b. Nonlinear Analysis and Applied Mathematics (NAAM) Research Group, Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia. Received 27 December 2018; received in revised form 4 February 2019; accepted 27 April 2019 KEYWORDS Abstract. Soret and Dufour e ects on MHD ow of Williamson uid between two rotating disks are examined in this study. Impacts of strati cation, viscous dissipation, Williamson uid; and activation energy are also considered. Bejan number and entropy generation for Rotating disks; the strati ed ow are discussed. The governing Partial Di erential Equations (PDEs) Activation energy; are converted into ODE using von Karmantransformations. The convergent solution Viscous dissipation; of complicated ODE is found using a homotopic procedure. The results of physical Entropy generation; quantities are discussed through plots and numerical values. It is noted that axial and Soret and Dufour radial velocities are higher in the case of greater Weissenberg number. Temperature and e ects; concentration pro les are the decreasing functions of thermal and solutal strati cation Strati cation. parameters, respectively. Entropy and Bejan number show the opposite trends for higher Weissenberg number and Brinkman number. © 2020 Sharif University of Technology. All rights reserved. 1. Introduction was numerically analyzed by Zehra et al. [3]. Khan et al. [4] and Qayyum et al. [5] studied entropy gen- Recent researchers have particularly focused on the eration during heat transfer of Williamson nano uid ows of non-Newtonian uids because these liquids due to stretching sheet and rotating disks, respectively. have penetrated industrial and technological processes. Further studies on non-Newtonian uids can be seen In this regard, they have suggested a number of models in [6{10]. for such liquids due to their diverse properties. In Non-Newtonian uid ow with rotating disks has general, subclasses of di erential-type liquids, known many applications in engineering, e.g., air cleaning as second, third, and fourth grades, are extensively machine, medical equipment and aerodynamical engi- analyzed. However, Williamson uid model has re- neering processes, etc. Karman [11] rstly examined ceived insigni cant attention. Williamson [1] studied the steady ow using an in nite rotating disk. The the pseudoplastic material experimentally. Nadeem et ow generated by rotating disks with radiative heat al. [2] discussed Williamson uid ow that results ux and variable thickness was examined by Hayat from the stretching velocity of a surface. This ow, et al. [12]. Khan et al. [13] presented MHD Eyring- which is characterized by pressure-dependent viscosity, Powell uid ow with a rotating disk. Doh and Muthtamilselvan discussed Micropolar uid ow by a *. Corresponding author. rotating disk with MHD and thermophoretic particle E-mail address: [email protected] (S. Qayyum) deposition e ects [14]. Khan et al. examined magnetic eld and double di usion in the couple stress uid ow doi: 10.24200/sci.2019.52553.2772 by a rotating disk [15]. Griths et al. addressed the 2452 T. Hayat et al./Scientia Iranica, Transactions B: Mechanical Engineering 27 (2020) 2451{2464 stability analysis of ow using a rotating disk [16]. 2. Formulation Qayyum et al. [17] discussed disorder in the system based on Williamson uid motion. Three-dimensional steady ow of Williamson uid be- It is a well-known fact that mass transfer exists in tween two rotating stretchable disks was considered. view of species concentration di erence in a mixture. This study aims to scrutinize FGM in ow with e ects The species characterized by di erent concentrations of Dufour/Soret, strati cation, and viscous dissipation. transport themselves in a mixture from a region with Chemical reaction with activation energy was also higher concentration to a region with lower concen- investigated. Flow takes place due to the stretching tration. Moreover, activation energy is the minimum of disks. The lower (z = 0) and upper (z = h) quantity of energy that must be possessed by reactants disks are characterized by respective angular velocities before any speci ed chemical reaction occurs. This of 1 and 2. Flow is caused by the stretching of lower disk. Lower and upper disks correspond to process, which is followed by a chemical reaction with ^ ^ ^ ^ activation energy in mass transfer, is usually used temperatures (T1; T2) and concentration (C1; C2). A in food processing, geothermal reservoirs, chemical uniform magnetic eld of strength (B0) is exerted in engineering, etc. Bestman investigated convection of the Z-direction (see Figure 1). a binary mixture owing in a porous space with acti- An extra stress tensor () of Williamson uid is given below: vation energy [18]. Makinde et al. discussed unsteady radiative ow in the presence of chemical reaction [19]. (~ ~ ) = ~ + 0 1 A; (1) Awad et al. described unsteady rotating ow with 1 1 _ activation energy and chemical reaction [20]. Sha que et al. studied the rotating ow containing chemically where A denotes the rst Rilvin-Erickson tensor, ~0 is reactive species and activation energy [21]. the zero shear rate viscosity, ~1 is the in nite shear In recent years, to nd optimal designs for en- rate viscosity, and > 0 is a time constant. Herein, _ is de ned below: gineering system, the second law of thermodynamics r has been applied to the analysis of Entropy Generation 1 _ = (trA2): (2) Minimization (EGM). EGM helps determine the grow- 2 ing rate of irreversibility during a process. Reversibility Through Eq. (1), one obtains the following: of heat and mass transfer and irreversibility of viscous dissipation can be determined by entropy generation. = (~0 + (~0 ~1) _ ) A; (3) During the convection process, irreversibility takes place inside the cavity. It is quite necessary to diminish which further yields: the irreversibility process to ensure energy conserva- @u^ = 2((~ ~ ) _ + ~ ) ; tion. To study the e ects of entropy generation inside rr 0 1 0 @r a thermal system, the second law of thermodynamics is considered. Firstly, Bejan introduced entropy gen- eration [22] and explained that entropy was generated due to viscous e ects and thermal conductivity. Amani and Nobari investigated entropy generation in a curved pipe at a constant wall temperature [23]. Hayat et al. studied peristaltic rotating ow of nanoparticles with entropy generation [24]. Shit et al. found that nano uid ow resulted from exponential stretching sheet and entropy generation [25]. Some of the recent works about entropy generation can be seen in [26{30]. Here, this study analyzes the e ects of entropy generation on MHD Williamson uid ow between two rotating disks. Further strati cation, viscous dissipa- tion, and activation energy e ects are considered. To the best of the authors' knowledge, such an attempt at studying Williamson uid has not been considered yet. Homotopy technique [31{42] is used to de- velop convergent solutions. Behaviors of temperature, velocity, entropy generation, Nusselt number, Bejan number, and skin friction are discussed via graphs and tabulated values. At the end, concluding remarks are given. Figure 1. Schematic diagram of the problem. T. Hayat et al./Scientia Iranica, Transactions B: Mechanical Engineering 27 (2020) 2451{2464 2453 !n 1 @v^ u^ T^ E = 2((~0 ~1) _ + ~0) + ; k2(C^ C^ ) exp ; (11) r @ r r 2 ^ ^ T2 T @w^ = 2 ((~ ~ ) _ + ~ ) ; zz 0 1 0 @z :L =[(~0 ~1) _ + ~0] 1 @u^ @v^ v^ @u^ 2 @w^ 2 v^ @v^ @v^ 2 = =((~ ~ ) _ + ~ ) + ; 2 + 2 2 + r r 0 1 0 r @ @r r @r @z r @r @r @w^ @u^ v^2 u^2 @u^ 2 @v^ 2 rz = zr = ((~0 ~1) _ + ~0) + ; + + 2 + + ; (12) @r @z r2 r2 @z @z 1 @w^ @v^ = = ((~ ~ ) _ + ~ ) + : (4) u^ = a1r; v^ = 1r; w^ = 0; z z 0 1 0 r @ @z ^ ^ ^ ^ ^ ^ Eq. (5) is shown in Box I. Mathematical statements T = T1 = m1r + T0; C = C1 = m3r + C0; of the problem under consideration satisfy the follow- ing [5]: at z = 0; u^ = a2r; v^ = 2r; p = 0; @u^ u^ @w^ + + = 0; (6) C^ = C^ = m r + C^ ; T^ = T^ = m r + T^ @r r @z 2 4 0 2 2 0 @u^ @u^ v^2 at z = h; (13) w +u ^ @z @r r where, in (x; y; z) directions, the velocity components @p @ @ are (u; v; w^), respectively, p the hydrostatic pressure, = + rr + zr + rr B2u;^ (7) @r @r @z r 0 Cp speci c heat, density of uid, k thermal conduc- tivity, Dm the e ective di usivity rate of mass, Cs the @v^ @v^ v^u^ u^ + w + susceptibility of concentration, KT thermal-di usion @r @z r ratio, kr the reaction rate, electrical conductivity, E activation energy where = 8:61 105eV/K the @r @z r 2 Boltzmann constant, n the tted rate constant, and = + + 2 B0 v;^ (8) @r @z r Tm mean temperature of uid. Suitable transformations for this analysis are @w @w @p @ @ u^ + w = + rz + zz + rz ; (9) given as follows: @r @z @z @r @z r ^ ^ @T^ @T^ k @2T^ @2T^ 1 @T^ u^ = r 1 f(); v^ = r 1 g^(); w^ = h 1 h(); u^ +w ^ = + + @r @z C @z2 @r2 r @r p C^ C^ p^ '^() = 2 ; P^() = ; ^ ^ 2 2 D K @2C^ @2C^ 1 @C^ C1 C0 h 1 +:L + m T + + ; (10) 2 2 CpCs @r @z r @r T^ T^ z ^() = 2 ; = : (14) T^ T^ h @C^ @C^ @2C^ @2C^ 1 @C^ 1 0 u^ +w ^ = D + + @r @z @r2 @z2 r @r Eq.
Load more

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]