Biological Fluid Flow Project 07

Total Page:16

File Type:pdf, Size:1020Kb

Biological Fluid Flow Project 07 Development of a Numerical Model to Describe Steady Streaming in Unsteady Axisymmetric Biological Flows in a Geometry with a Slowly Varying Cross-Section Laboratory of Computational Engineering Dr. Mark Sawley Laboratory of Fluid Mechanics Prof. Peter Monkewitz - Radboud Nelissen Development of a Numerical Model to Describe Steady Streaming in Unsteady Axisymmetric Biological Flows in a Geometry with a Slowly Varying Cross-Section Laboratory of Computational Engineering Dr. Mark Sawley Laboratory of Fluid Mechanics Prof. Peter Monkewitz - Radboud Nelissen Anthony Claude Abstract Cerebrospinal fluid moves from the cranial cavities all along the cerebrospinal canal in a pulsatile way. Its flow has been studied here through various geometries. From the basic case of a tube to a model of an idealized anatomical case, velocity profile and flow behaviour were plotted, analyzed and compared. This bi-dimensional study takes into account variation of geometry of the cerebrospinal canal and pulsatile character of the flow. A data base relating variation of internal and external wall for an adult has been realized with Visual Human Server1. This research allows determining Reynolds and Womersley numbers average for a given region of the spinal canal. Quantitative solutions of unsteady analysis through a cylindrical tube and qualitative results concerning steady streaming in geometries with slowly variations of section have been found . Finally comparison were carried out between computational calculations made in 1998 by A. Sarkar and G. Jayaraman2 and Fluent calculations for a similar issue. 1http://visiblehuman.epfl.ch/vhve.php 2Correction to flow rate - pressure drop relation in coronary angioplasty: steady straeming effect Contents I Basic Fluid Mechanics Applied to Biological Flows 4 II Womersley Number’s Influence on Biological Flows 5 1 Introduction 5 2 Pulsatile flow through a tube 6 2.1 Hypotheses . 6 2.2 Velocity’s expression . 6 2.3 Results’ Interpretation . 9 2.4 RMS-velocity’s expression . 11 3 Pulsatile flow through an annular tube 12 3.1 Velocity’s expression . 12 3.2 Results’ Interpretation . 13 III Numerical Flow Modelling 14 4 Introduction 14 5 Steady state flow - analytical vs. numerical results 15 5.1 Analytical solution . 15 5.2 Numerical solution . 18 5.3 Results’ interpretation . 20 6 CSF flow through an axisymmetric annular model 21 6.1 Generating the mesh . 22 6.2 Results’ interpretation . 23 7 CSF flow through 2-D models with varying cross-sections 33 7.1 Visual Human Server - Database . 33 7.2 Theoritecal part . 35 7.3 Axisymmetric geometry with a linear variation of section . 37 7.4 Axisymmetric geometry with a half cosinusoidal variation of section . 39 7.5 Axisymmetric geometry with a cosinusoidal variation of geometry . 40 7.6 Sarkar and Jayaraman case . 43 8 Conclusion 46 9 Bibliography 47 Introduction Biofliudics has gained in importance in recent years, forcing engineers to use fluid me- chanic and apply it to biological situation. Cardiovascular and more generally every human flow circulation problems are very important nowadays. Because of the multiple diseases caused by flow comportment, it is fundamental to understand its characteristics and the rules by which they operate.. Why does an artery get stenosed? How could we optimized drug dispersion in the cerebrospinal canal? What is the influence of a cardiac valve on the blood flow? By answering these questions, we could offer better treatments and solutions to patients. At the EPFL, some experimental researches are conducted on the spread of drugs ad- ministered directly in the spinal canal (with a catheter). Furthermore, this project will also focus on cerebrospinal fluid flow. Its aim is to evaluate the potential of Fluent solver to calculate accurate solution for oscillatory flow. With the development of numerical simulation possibilities and the increase of preci- sion and liability of computer programs, research progresses faster. But technology progress is not enough, we still need an engineer to check, analyse and interpret results. That is why good knowledge in fluid mechanics is needed before solving a problem numerically. To what extent is it possible to consider the results obtained as an exact solution? Is it possible to control the precision, and if so, in what way? This project follows this path. Unfortunately the spinal canal geometry is extremely complex and varies between peo- ple. Considering that and the calculation time needed, it is not appropriate to model precisely the cerebrospinal canal. In order to study flow behaviour of these geometries with varying cross-section, an axisymmetric analytical model will be used. 4 Part I Basic Fluid Mechanics Applied to Biological Flows As every incompressible fluid, biological one obeys the different conservation principles. Traditionally we consider the conservation of the mass, the momentum and the energy: dρ + ρ∇ · v = 0 (1) dt Considering the incompressibility condition (ρ = cte) we obtain: ∇ · v = 0 (2) The conservation of the momentum using the Gauss theorem (local form) is ∂ (ρv) + ∇ · (ρvv) = ∇ · σ (3) ∂t Biological fluid flows are unsteady. Thus, they are driven by the gravitational force and the pressure gradient force (cardiac cycle governs this pulsatile flow). Some shear forces are also noticed(caused by the fluid’s viscosity).Many biological flow issue can be discribed analytically. For instance, recirculation caused by a stenosis in an artery. This phenomenon is important in hemodynamic because there is a risk of canal’s obstruction. Considering the cerebrospinal canal, we encounter the same problem. The difference is that the variation of geometry is due to the change of internal and external radius or to the nerve’s insertions and not to deposits.Most of the biological incompressible problems are solved using the Navier-Stokes equations with cylindrical coordinate. Navier-Stokes equation (in cylindrical coordinate) are: 2 ∂vr ∂vr vθ ∂vr ∂vr vθ • ρ ∂t + vr ∂r + r ∂θ + vz ∂z − r 2 2 2 ∂p ∂ vr 1 ∂ vr ∂ vr 1 ∂vr 2 ∂vθ vr = − ∂r + η ∂r2 + r2 ∂θ2 + ∂z2 + r ∂r − r2 ∂θ − r2 + ρfr ∂vθ ∂vθ vθ ∂vθ ∂vθ vθvr • ρ ∂t + vr ∂r + r ∂θ + vz ∂z + r 2 2 2 1 ∂p ∂ vθ 1 ∂ vθ ∂ vθ 1 ∂vθ 2 ∂vr vθ = − r ∂θ + η ∂r2 + r2 ∂θ2 + ∂z2 + r ∂r + r2 ∂θ − r2 + ρfθ ∂vz ∂vz vθ ∂vz ∂vz • ρ ∂t + vr ∂r + r ∂θ + vz ∂z 2 2 2 ∂p ∂ vz 1 ∂ vz ∂ vz 1 ∂vz = − ∂z + η ∂r2 + r2 ∂θ2 + ∂z2 + r ∂r + ρfz 5 Part II Womersley Number’s Influence on Biological Flows 1 Introduction Before undertaking more complex analysis of the cerebro-spinal canal, a study of an oscillatory fluid flow through a simple geometry is recommanded. The first part of this project is dedicated to an analysis of a flow moving by a pressure gradient (representing the cardiac cycle) in a cylindrical tube. It is crucial to understand the basic physics of the issue before developing a numeri- cal model. Therefore, the following development can be viewed as a theoretical approach of unsteady flows, particularly, of biological flows. The first step will be to solve Navier-Stokes equation in order to find an expression for the velocity. Once this is done, the profile shape can be plotted and observed. Finally, a discussion of the physical interpretation and of the parameters’ influence will take place. The aim here is to define the influence of the Womersley number on a fluid flow and its effect on the velocity profile. Cerebrospinal fluid - functions and properties The main function of the cerebrospinal fluid is to protect the brain. In fact, by bathing in this fluid, it does not feel effect of the pressure fluctuations. Obviously there are some limits and this protection is efficient as far as we are in a standard situation. Normally, CSF volume and blood volume vary inversely, maintaining the intracranial pressure within normal limits. CSF has also several other important functions: it regulates chemical environment of the central nervous system and it is active in the transportation of biological substances like proteins. It is principally composed by plasma with few proteins and cells and it is considered as a Newtonian flow g m2 Density ρ [ m3 ] Viscosity µ [ s ] CSF 993 0.0067 Table 1: CSF’s properties 2 PULSATILE FLOW THROUGH A TUBE 6 2 Pulsatile flow through a tube As a first approximation, the cerebrospinal canal will be modeled as a cylindrical tube. By doing so, the equations are simpler but further from reality. The main objective here is to study the influence of the Womersley number on a fluid flow and on its velocity profile. Figure 1: One-dimensional fluid flow in a tube 2.1 Hypotheses • Incompressible flow, no external forces ∂P iωt • Pressure grandient is periodic (sinus): ∂z = Ge , with G constant • One-dimensional fluid flow: (0, 0,W ) avec W = W (R, t) ∂ • Axisymmetric flow: ∂θ = 0 D • Boundary condition: W ( 2 , t) = 0 2.2 Velocity’s expression After some simplifications, Navier-Stokes equations become: ∂W ∂P 1 ∂ ∂W ρ = − + µ (R ) (4) ∂t ∂z R ∂R ∂R ∂P ∂P = = 0 (5) ∂R ∂θ ∂P = Geiωt (6) ∂z To introduce the Womersley number and to make the calculation easier, it is prefer- able to make this problem adimensional: τ t = ω R = R0r W = W0w 2 PULSATILE FLOW THROUGH A TUBE 7 q ω Womersley number α: α = R0 ν W0R0 Reynolds number Re: Re = ν ∂W 1 ∂P 1 ∂ ∂W = − + ν (R ) (7) ∂t ρ ∂z R ∂R ∂R ∂w 2 ∂p ν 1 ∂ ∂w ωW0 = −W0 + rW0 (8) ∂τ ∂z R0r R0 ∂r ∂r ∂w ∂p ∂2w 1 ∂w α2 = −ReR + + (9) ∂τ 0 ∂z ∂r2 r ∂r As mentioned aboved, the pressure gradient oscillates. It’s expression is: ∂p = Geiωt = Geiτ (10) ∂z The consequence of this is that the velocity will be oscillatory too: w(r, t) =w ˆ(r)eiωt =w ˆ(r)eiτ (11) To simplify the notation, wˆ(r, t) will be noted w hereafter.
Recommended publications
  • Wall Thermal Inertia Effects of Pulsatile Flow in a Ribbed Tube: a Numerical Approach Amin Kardgar*,A, Davood Domiri Ganjib
    Journal of Heat and Mass Transfer Research 7 (2020) 85-94 Journal of Heat and Mass Transfer Research Journal homepage: http://jhmtr.semnan.ac.ir Semnan University Wall Thermal Inertia Effects of Pulsatile Flow in a Ribbed Tube: A Numerical Approach Amin Kardgar*,a, Davood Domiri Ganjib aDepartment of Mechanical Engineering, University of Mazandaran, Babolsar, Iran. bDepartment of Mechanical Engineering, Babol University of Technology, Babol, Iran. PAPER INFO ABSTRACT Paper history: In the present paper, the heat transfer of pulsatile flow in ribbed tube was investigated Received: 2020-02-24 numerically by considering the effect of thermal inertia of solid wall thickness. To this purpose, the CVFV technique with collocated grids arrangement was adopted to discretize Revised: 2020-07-01 momentum and energy equations. In order to avoid checker-board of pressure field in Accepted: 2020-10-16 numerical simulation, Chow and Rhie interpolation scheme was employed. The well- established SIMPLE (Semi-Implicit Method for Pressure Linked Equations) method was Keywords: utilized to handle velocity and pressure coupling in the momentum equation. Stone’s Strongly Pulsatile flow; Implicit Procedure (SIP) was used to solve the set of individual linear algebraic equations. Conjugate heat transfer; Womersley number, Reynolds number, velocity amplitude, and wall thickness ratio are four Womersley number; essential parameters that influence heat transfer and Nusselt number in pulsatile flow in a Wall thickness; ribbed tube. It was deduced by varying Womersley number Nu does not change. Nu enhances ribbed tube. almost 19% by augmentation of wall thickness ratio from 0.125 to 1. It was shown by increasing velocity amplitude from 0.1 to 0.8; Nu reduces almost 4.7%.
    [Show full text]
  • Introduction to Biofluid Mechanics Portonovo S
    To protect the rights of the author(s) and publisher we inform you that this PDF is an uncorrected proof for internal business use only by the author(s), editor(s), reviewer(s), Elsevier and typesetter TNQ Books and Journals Pvt Ltd. It is not allowed to publish this proof online or in print. This proof copy is the copyright property of the publisher and is confidential until formal publication. CHAPTER 16 Introduction to Biofluid Mechanics Portonovo S. Ayyaswamy OUTLINE 16.1 Introduction e2 Exercises e71 16.2 The Circulatory System in the Acknowledgments e72 Human Body e2 Literature Cited e72 16.3 Modeling of Flow in Blood Vessels e18 Supplemental Reading e73 16.4 Introduction to the Fluid Mechanics of Plants e65 CHAPTER OBJECTIVES • To properly introduce the subject of biofluid • To review the parametric impact of the prop- mechanics including the necessary language erties of rigid, flexible, branched, and curved fl • To describe the components of the human cir- tubes on blood ow culation system and document their nominal • To provide an overview of fluid transport in characteristics plants • To present analytical results of relevant models of steady and pulsatile blood flow Fluid Mechanics http://dx.doi.org/10.1016/B978-0-12-405935-1.00016-2 e1 Copyright © 2016 Elsevier Inc. All rights reserved. To protect the rights of the author(s) and publisher we inform you that this PDF is an uncorrected proof for internal business use only by the author(s), editor(s), reviewer(s), Elsevier and typesetter TNQ Books and Journals Pvt Ltd. It is not allowed to publish this proof online or in print.
    [Show full text]
  • CFD Modelling of Air Flow and Fine Powder Deposition in the Respiratory Tract
    CFD modelling of air flow and fine powder deposition in the respiratory tract Yun Hwan Kim A thesis in fulfilment of the requirements for the degree of Master of Philosophy School of Materials Science and Engineering Faculty of Science April 2017 THE UNIVERSITY OF NEW SOUTH WALES Thesis/Dissertation Sheet Surname or Family name: Kim First name: Yun Hwan Other name/s: Abbreviation for degree as given in the University calendar: M.Phil. School: School of Materials Science and Engineering Faculty: Science Title: CFD modelling of air flow and fine powder deposition in the respiratory tract Abstract This project was to investigate and observe characteristics of micro particles suspended in the ambient air or pharmaceutical aerosols with respect to the mechanisms of deposition in human airways under different inspiratory conditions. Such determination includes pattern observations of inspiratory flow-field of the air, particle trajectory during inspiratory conditions and particle deposition. Computational fluid dynamic (CFD) was employed to simulate above problems, aiming to observe flow-field of the inspiratory air and characteristic of flow turbulence in the respiratory tract as well as particle behaviour in the respiratory tract regarding to the particle deposition. In order to do so, three different airway models were employed for the simulations: two realistic airway models introduced by Kitaoka and Weibel airways model. The motion of micro-sized particles between 1~20 µm were simulated under the steady state two inlet- inspiratory conditions – inhalation condition (60 L/min) and breathing condition (18 L/min); to evaluate deposition efficiency. Inertial impaction was dominantly caused high density deposition of particles in upper tracheobronchial region, particularly in regions where daughter airways bifurcate.
    [Show full text]
  • Mass Transfer and Separation of Species in Oscillating Flows
    MASS TRANSFER AND SEPARATION OF SPECIES IN OSCILLATING FLOWS By AARON M. THOMAS A DISSERTATION PRESENTED TO THE GRADUATE SCHOOL OF THE UNIVERSITY OF FLORIDA IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY UNIVERSITY OF FLORIDA 2001 Copyright 2001 by AARON M. THOMAS To Susan, Yesenia, and Dakota. ACKNOWLEDGMENTS guidance on this work I am most grateful to my advisor, Dr. Narayanan, for his and his counsel on a variety of topics that arose during our numerous conversations. some of I also thank Dr. Kurzweg and Dr. Jaeger for their help on understanding the aspects of this research. for a large part of my I am fortunate that Dr. Weaver has taken time to mentor me education. My parents and sisters have been very supportive of me in all of my work and decisions, and I know I can never repay everything that they have done for me. have made The friends that I have made both inside and outside this department my experience in Gainesville most memorable. Most importantly, my wife continually makes me realize that nothing has as much meaning without her by my side. My wife and my children properly balance my life and make me understand what is truly important. from under I gratefully acknowledge the support for this project received NASA grant number NAG3 2415 and for my financial support through an NSF Fellowship under grant number 9616053. IV TABLE OF CONTENTS page ACKNOWLEDGMENTS i v LIST OF TABLES viii LIST OF FIGURES « LIST OF SYMBOLS xiii 1.3 ABSTRACT xvii 1 INTRODUCTION OF THE PHYSICS 1 1.1 Physics of the
    [Show full text]
  • A Computer Simulation of the Non-Newtonian Blood Flow at the Aortic Bifurcation
    Do21 9290/93 %S.W+.oO J. Biomechmics Vol. 26, NO. I, pp. 37-49. 1993. c#; 1992 Pergamon Press Ltd Pnnted m Great Britain A COMPUTER SIMULATION OF THE NON-NEWTONIAN BLOOD FLOW AT THE AORTIC BIFURCATION ZHENG Lou* and WEN-JEI YANG? *Transportation Research Institute, University of Michigan, Ann Arbor, MI 48109, U.S.A. and tMechanicil Engineering and Applied Mechanics, University of Michigan, Ann Arbor, MI 48109, U.S.A. Abstract-A two-dimensional numerical model was developed to determine the effect of the non-New- tonian behavior of blood on a pulsatile flow at the aortic bifurcation. The blood rheology was described by a weak-form Casson equation. The successive-over-relaxation (SOR) method was used to solve both the vorticity and Poisson equations numerically. It was disclosed that the non-Newtonian property of blood did not drastically change the flow patterns, but caused an appreciable increase in the shear stresses and a slightly higher resistance to both flow separations and the phase shifts between flow layers. NOMENCLATURE INTRODUCTION AR area ratio The curvatures, junctions and bifurcations of large D diameter and medium arteries are severely affected by athero- second invariant of the strain rate tensor Dll sclerosis. It has been suggested that fluid dynamics H half-width of a channel J, .i and J’ coordinate transformation coefficients may play a role in the genesis and progression of the I, reference length disease. The blood flows in these regions are disturbed Lb min minimum branch length and the arterial walls are exposed to either high or low parameter in the Papanastasiou’s weak form shear stresses.
    [Show full text]
  • Numerical, Analytical, and Experimental Studies Of
    Florida International University FIU Digital Commons FIU Electronic Theses and Dissertations University Graduate School 11-14-2017 Numerical, Analytical, and Experimental Studies of Reciprocating Mechanism Driven Heat Loops for High Heat Flux Cooling Olubunmi Tolulope Popoola Florida International University, [email protected] DOI: 10.25148/etd.FIDC004065 Follow this and additional works at: https://digitalcommons.fiu.edu/etd Part of the Heat Transfer, Combustion Commons, and the Other Mechanical Engineering Commons Recommended Citation Popoola, Olubunmi Tolulope, "Numerical, Analytical, and Experimental Studies of Reciprocating Mechanism Driven Heat Loops for High Heat Flux Cooling" (2017). FIU Electronic Theses and Dissertations. 3505. https://digitalcommons.fiu.edu/etd/3505 This work is brought to you for free and open access by the University Graduate School at FIU Digital Commons. It has been accepted for inclusion in FIU Electronic Theses and Dissertations by an authorized administrator of FIU Digital Commons. For more information, please contact [email protected]. FLORIDA INTERNATIONAL UNIVERSITY Miami, Florida NUMERICAL, ANALYTICAL, AND EXPERIMENTAL STUDIES OF RECIPROCATING MECHANISM DRIVEN HEAT LOOPS FOR HIGH HEAT FLUX COOLING A dissertation submitted in partial fulfilment of the requirements for the degree of DOCTOR OF PHILOSOPHY in MECHANICAL ENGINEERING by Olubunmi Tolulope Popoola 2017 To: Dean John L. Volakis College of Engineering and Computing, This dissertation, written by Olubunmi Tolulope Popoola, Numerical, Analytical, and
    [Show full text]
  • Vortexformationinthecardiovasc.Pdf
    Vortex Formation in the Cardiovascular System Arash Kheradvar • Gianni Pedrizzetti Vortex Formation in the Cardiovascular System Authors Arash Kheradvar, M.D., Ph.D. Gianni Pedrizzetti, Ph.D. University of California, Irvine Università degli Studi di Trieste 2410 Engineering Hall Dipto. Ingegneria Civile e The Edwards Lifesciences Center for Architettura Advanced Cardiovascular Technology P. Europa 1 Irvine 34127 Trieste California 92697 Italy USA ISBN 978-1-4471-2287-6 e-ISBN 978-1-4471-2288-3 DOI 10.1007/978-1-4471-2288-3 Springer London Dordrecht Heidelberg New York British Library Cataloguing in Publication Data A catalogue record for this book is available from the British Library Library of Congress Control Number: 2011944219 © Springer-Verlag London Limited 2012 Apart from any fair dealing for the purposes of research or private study, or criticism or review, as permitted under the Copyright, Designs and Patents Act 1988, this publication may only be reproduced, stored or transmitted, in any form or by any means, with the prior permission in writing of the publishers, or in the case of reprographic reproduction in accordance with the terms of licenses issued by the Copyright Licensing Agency. Enquiries concerning reproduction outside those terms should be sent to the publishers. The use of registered names, trademarks, etc., in this publication does not imply, even in the absence of a specifi c statement, that such names are exempt from the relevant laws and regulations and therefore free for general use. Product liability: The publisher can give no guarantee for information about drug dosage and application thereof contained in this book.
    [Show full text]
  • Linear and Nonlinear Dynamics of Pulsatile Channel Flow Benoît Pier, Peter J
    Linear and nonlinear dynamics of pulsatile channel flow Benoît Pier, Peter J. Schmid To cite this version: Benoît Pier, Peter J. Schmid. Linear and nonlinear dynamics of pulsatile channel flow. Journal of Fluid Mechanics, Cambridge University Press (CUP), 2017, 815, pp.435-480. 10.1017/jfm.2017.58. hal-01385093v3 HAL Id: hal-01385093 https://hal.archives-ouvertes.fr/hal-01385093v3 Submitted on 6 Mar 2017 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. J. Fluid Mech. (2017), vol. 815, pp. 435–480. 435 doi: 10.1017/jfm.2017.58 Linear and nonlinear dynamics of pulsatile channel flow Benoˆıt P I E R1 and Peter J. S C H M I D2 † 1Laboratoire de m´ecanique des fluides et d’acoustique, CNRS, Ecole´ centrale de Lyon, Universit´ede Lyon 1, INSA Lyon, 36 avenue Guy-de-Collongue, 69134 Ecully,´ France 2Department of Mathematics, Imperial College, South Kensington Campus, London SW7 2AZ, UK (Received 14 July 2016; revised 19 January 2017; accepted 19 January 2017) The dynamics of small-amplitude perturbations as well as the regime of fully developed nonlinear propagating waves is investigated for pulsatile channel flows.
    [Show full text]
  • Basic Field Equations
    APPEND IX Basic Field Equations A.I Introduction In this appendix we list the basic equations of mechani cs. We present these equations in order to exhibit the mathematical expressions of the basic laws of physics (conservation of mass, momentum, and energy) and to show the notations used in this book. How easily the material presented here will be assimilated depends on the background of the reader. For some, a quick glance will be sufficient. For others, some hard study and further reading will be necessary. The author' s book, A First Course in Continuum M echanics (Fung, 1977), is one of the shortest and most con venient references for this purpose. Other larger treati ses are listed at the end of this appendix. The mechanical properties of materials are expressed in constitutive equations. The constitutive equations of biological materials are discussed in the author's book, Biomechanics: Mechanical Properties ofLiving Tissues (Fung, 1981). They are often quoted below without discussion. A.2 Conservation of Mass and Momentum Let us consider the motion of a solid or a fluid. Let XI> X2 'X 3 or X, y, z be rectangular Cartesian coordinates . Let the velocity components along the X-, y-, z-axis directions be denoted by VI'V2' V3 or u, v, IV, respectively. Let p denote pressure, and uij or (ux x ' uYY' Uzz ' Ux }.' uyx ' u zz) be the stress components. Throughout this book the index notation will be used. Unless stated otherwise, all indices will range over 1,2,3. The summation convention will be used : repetition of an index means summation over I, 2, 3.
    [Show full text]
  • Effect of Prandtl Number on Internal Convective Heat Transfer in Laminar Single-Phase Pulsating Flows
    Effect of Prandtl number on internal convective heat transfer in laminar single-phase pulsating flows Sushant Kumar 1, Balkrishna Mehta 2, Sameer Khandekar 2a 1Birla Institute of Technology & Science Pilani – Goa Campus, Goa, India 2Indian Institute of Technology Kanpur, Kanpur, India (aCorresponding author, email: [email protected]) Abstract. In this paper, heat transfer during simultaneously developing single- phase laminar flow inside a pipe has been investigated numerically. A two di- mensional axisymmetric model of a circular mini-channel is used, with a con- stant wall temperature boundary condition. In the present simulations, flow Reynolds number ( Re ) and amplitude ratio ( A) are kept constant to 200 and 0.2, respectively. The variable quantities are imposed pulsating frequency ( f) or non- dimensional frequency (Womersley number, Wo ) and Prandtl number ( Pr ). Study is conducted for three different imposed flow frequencies of 5 Hz, 15 Hz and 25 Hz respectively while Prandtl number values are varied from 2.43 to 3404. These flow parameters can be grouped together to form different time scales which are important to explain the heat transfer behavior in laminar sin- gle-phase pulsating flow. It is observed that for the range of frequency (f) and Prandtl number ( Pr ) used in present work, change in heat transfer is marginal which is of no practical significance. Enhancement of heat transfer coefficient due to such periodic pulsatile internal flows, over and above the non-pulsatile regular flow conditions, is questionable, and at best, rather limited. 1 Introduction The occurrence of single-phase pulsating flow can be encountered in many engi- neering systems, such as bio-fluidic devices, cooling systems for miniaturized elec- tronic circuits, reciprocating pumps, etc.
    [Show full text]
  • MHD Pulsatile Flow of Engine Oil Based Carbon Nanotubes Between Two Con- Centric Cylinders
    Accepted Manuscript MHD Pulsatile Flow of Engine oil based Carbon Nanotubes between Two Con- centric Cylinders Rizwan Ul Haq, Faisal Shahzad, Qasem M. Al-Mdallal PII: S2211-3797(16)30327-8 DOI: http://dx.doi.org/10.1016/j.rinp.2016.11.057 Reference: RINP 457 To appear in: Results in Physics Received Date: 7 October 2016 Revised Date: 23 November 2016 Accepted Date: 25 November 2016 Please cite this article as: Ul Haq, R., Shahzad, F., Al-Mdallal, Q.M., MHD Pulsatile Flow of Engine oil based Carbon Nanotubes between Two Concentric Cylinders, Results in Physics (2016), doi: http://dx.doi.org/10.1016/ j.rinp.2016.11.057 This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain. MHD Pulsatile Flow of Engine oil based Carbon Nanotubes between Two Concentric Cylinders Rizwan Ul Haq1,∗, Faisal Shahzad 2 1Department of Electrical Engineering, Bahria University, Islamabad Campus, Islamabad, 44000, Pakistan 2Department of Mathematics, Capital University of Science and Technology, Islamabad 44000, Pakistan”” Abstract Present analysis is dedicated to examine the thermal performance of engine oil in the presence of both single and multiple wall carbon nanotubes (SWCNTs and MWCNTs) between two concentric cylinders. Flow is driven with oscillatory pressure gradient and magneto- hydrodynamics (MHD) effects are also introduced to control the random motion of the nanoparticles.
    [Show full text]
  • Effect of Grid Topology and Resolution on Computation of Steady and Unsteady Internal Flows
    The Pennsylvania State University The Graduate School College of Engineering EFFECT OF GRID TOPOLOGY AND RESOLUTION ON COMPUTATION OF STEADY AND UNSTEADY INTERNAL FLOWS A Thesis in Mechanical Engineering by Steven P. McHale c 2009 Steven P. McHale Submitted in Partial Fulllment of the Requirements for the Degree of Master of Science May 2009 The Thesis of Steven P. McHale has been reviewed and approved* by the follow- ing: Eric G. Paterson Associate Professor of Mechanical Engineering Thesis Advisor Laura L. Pauley, P.E. Arthur L. Glenn Professor of Engineering Education Professor of Mechanical Engineering Karen A. Thole Professor of Mechanical Engineering Department Head of Mechanical and Nuclear Engineering *Signatures are on le in the Graduate School Abstract Internal ow has been studied under steady and unsteady conditions. The work was motivated by previous research in canine olfaction which showed that the structure of the dog's nose consists of a highly complex pipe network in which the unsteady air ow transports odorant molecules to the sensory region. Due to geometric complexity of the canine turbinates, a body-tted hex- dominant mesher was used, thus generating the question: What is the impact of grid topology and resolution on solution accuracy? The governing equations are dened and reduced to an analytical solution for a circular pipe with both steady and unsteady components. The numerical methods and solution algorithms are described. Simple test cases are performed using an open source computational uid dynamics (CFD) code, OpenFOAM, to understand the underlying physics and computational challenges inherent in internal pipe ows. The characteristic steady, Poiseuille ow was found and a grid study was performed over a range of Womersley numbers (Wo) for the oscillatory, unsteady calculations of the straight pipe.
    [Show full text]