Basics of Fluid Mechanics Genick Bar-Meir, Ph. D. 1107 16thAve S. E. Minneapolis, MN 55414-2411 email:barmeir@gmail.com Copyright ©@2008, 2007, and 2006 by Genick Bar-Meir See the file copying.fdl or copyright.tex for copying conditions. Version (0.1.6 August 11, 2008)  'We are like dwarfs sitting on the shoulders of giants" from The Metalogicon by John in 1159  CONTENTS Nomenclature xi GNU Free Documentation License . . . . . . . . . . . . . . . . . . . . . xv 1. APPLICABILITY AND DEFINITIONS . . . . . . . . . . . . . . xvi 2. VERBATIM COPYING . . . . . . . . . . . . . . . . . . . . . . . xvii 3. COPYING IN QUANTITY . . . . . . . . . . . . . . . . . . . . . . xviii 4. MODIFICATIONS . . . . . . . . . . . . . . . . . . . . . . . . . . . xviii 5. COMBINING DOCUMENTS . . . . . . . . . . . . . . . . . . . . xx 6. COLLECTIONS OF DOCUMENTS . . . . . . . . . . . . . . . . xxi 7. AGGREGATION WITH INDEPENDENT WORKS . . . . . . . xxi 8. TRANSLATION . . . . . . . . . . . . . . . . . . . . . . . . . . . xxi 9. TERMINATION . . . . . . . . . . . . . . . . . . . . . . . . . . . xxi 10. FUTURE REVISIONS OF THIS LICENSE . . . . . . . . . . . . . xxii ADDENDUM: How to use this License for your documents . . . . . . xxii How to contribute to this book . . . . . . . . . . . . . . . . . . . . . . . xxiii C redits . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . xxiii Steven from artofproblemsolving.com . . . . . . . . . . . . . . . . xxiii Dan Olsen . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . xxiv Richard Hackbarth . . . . . . . . . . . . . . . . . . . . . . . . . . . xxiv Tousher Yang . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . xxiv Your name here . . . . . . . . . . . . . . . . . . . . . . . . . . . . xxiv Typo corrections and other "minor" contributions............xxiv Version 0.1.8 August 6, 2008 . . . . . . . . . . . . . . . . . . . . . . . . xxxiii pages 189 size 2.6M . . . . . . . . . . . . . . . . . . . . . . . . . . . xxxiii Version 0.1 April 22, 2008 . . . . . . . . . . . . . . . . . . . . . . . . . . . xxxiii pages 151 size 1.3M . . . . . . . . . . . . . . . . . . . . . . . . . . . xxxiii Properties ..... ................................ xxxix  ii CONTENTS Open Channel Flow . . . . . . . . . . . . . . . . . . . . . . . . . . . xxxix 1 Introduction 1 1 .1 What is Fluid Mechanics? . . . . . . . . . . . . . . . . . . . . . . . . 1 1.2 Brief History. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3 1.3 Kindsof Fluids . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 1.4 Shear Stress . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6 1.5 Viscosity. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9 1.5.1 General . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9 1.5.2 Non-Newtonian Fluids. . . . . . . . . . . . . . . . . . . . . . 10 1.5.3 Kinematic Viscosity . . . . . . . . . . . . . . . . . . . . . . . 11 1.5.4 Estimation of The Viscosity . . . . . . . . . . . . . . . . . . . 12 1.5.5 Bulk Modulus . . . . . . . . . . . . . . . . . . . . . . . . . . . 19 1.6 Surface Tension . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22 1.6.1 Wetting of Surfaces . . . . . . . . . . . . . . . . . . . . . . . 24 2 Review of Thermodynamics 33 2.1 Basic Definitions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33 3 Review of Mechanics 41 3.1 Center of Mass . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41 3.1.1 Center of the Mass . . . . . . . . . . . . . . . . . . . . . . . . 41 3.1.2 Center of Area . . . . . . . . . . . . . . . . . . . . . . . . . . 42 3.2 Moment of Inertia . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43 3.2.1 Moment of Inertia for Mass . . . . . . . . . . . . . . . . . . . 43 3.2.2 Moment of Inertia for Area . . . . . . . . . . . . . . . . . . . . 44 3.2.3 Examples of Moment of Inertia . . . . . . . . . . . . . . . . . 46 3.2.4 Product of Inertia . . . . . . . . . . . . . . . . . . . . . . . . . 48 3.2.5 Principal Axes of Inertia . . . . . . . . . . . . . . . . . . . . . 50 3.3 Newton's Laws of Motion . . . . . . . . . . . . . . . . . . . . . . . . 50 3.4 Angular Momentum and Torque . . . . . . . . . . . . . . . . . . . . . 51 3.4.1 Tables of geometries . . . . . . . . . . . . . . . . . . . . . . 52 4 Fluids Statics 55 4.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55 4.2 The Hydrostatic Equation . . . . . . . . . . . . . . . . . . . . . . . . 55 4.3 Pressure and Density in a Gravitational Field . . . . . . . . . . . . . 57 4.3.1 Constant Density in Gravitational Field . . . . . . . . . . . . . 57 4.3.2 Pressure Measurement . . . . . . . . . . . . . . . . . . . . . 59 4.3.3 Varying Density in a Gravity Field . . . . . . . . . . . . . . . . 61 4.3.4 The Pressure Effects Because Temperature Variations . . . . 65 4.3.5 Gravity Variations Effects on Pressure and Density . . . . . . 69 4.3.6 Liquid Phase . . . . . . . . . . . . . . . . . . . . . . . . . . . 71 4.4 Fluid in a Accelerated System . . . . . . . . . . . . . . . . . . . . . . 72 4.4.1 Fluid in a Linearly Accelerated System . . . . . . . . . . . . . 72  CONTENTS III 4.4.2 Angular Acceleration Systems: Constant Density . 4.5 Fluid Forces on Surfaces . . . . . . . . . . . . . . . . . . 4.5.1 Fluid Forces on Straight Surfaces . . . . . . . . . . 4.5.2 Force on Curved Surfaces . . . . . . . . . . . . . . 4.6 Buoyancy and Stability . . . . . . . . . . . . . . . . . . . . 4.6.1 Stability . . . . . . . . . . . . . . . . . . . . . . . . 4.6.2 Surface Tension . . . . . . . . . . . . . . . . . . . 4.7 Rayleigh-Taylor Instability . . . . . . . . . . . . . . . . . . 5 Multi-Phase Flow 5.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . 5.2 History . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5.3 What to Expect From This Chapter . . . . . . . . . . . . . 5.4 Kind of Multi-Phase Flow . . . . . . . . . . . . . . . . . . 5.5 Classification of Liquid-Liquid Flow Regimes . . . . . . . 5.5.1 Co-Current Flow . . . . . . . . . . . . . . . . . . . 5.6 Multi-Phase Flow Variables Definitions . . . . . . . . . . . 5.6.1 Multi-Phase Averaged Variables Definitions . . . . 5.7 Homogeneous Models . . . . . . . . . . . . . . . . . . . . 5.7.1 Pressure Loss Components . . . . . . . . . . . . . 5.7.2 Lockhart Martinelli Model . . . . . . . . . . . . . . 5.8 Solid-Liquid Flow . . . . . . . . . . . . . . . . . . . . . . . 5.8.1 Solid Particles with Heavier Density ps > PL . . . 5.8.2 Solid With Lighter Density ps < p and With Gravity 5.9 Counter-Current Flow . . . . . . . . . . . . . . . . . . . . 5.9.1 Horizontal Counter-Current Flow . . . . . . . . . . 5.9.2 Flooding and Reversal Flow . . . . . . . . . . . . . 5.10 Multi-Phase Conclusion . . . . . . . . . . . . . . . . . . . Index Subjects Index . . . . . . . . . . . . . . . . . . . . . . . . . . . Authors Index . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74 75 75 85 92 98 107 108 113 113 113 114 115 116 117 122 122 125 126 128 129 130 132 133 136 136 143 145 .145 .147  iv CONTENTS  LIST OF FIGURES 1.1 Diagram to explain part of relationships of fluid mechanics branches. 2 1.2 Density as a function of the size of sample . . . . . . . . . . . . . . 6 1.3 Schematics to describe the shear stress in fluid mechanics. . . . . . 6 1.4 The deformation of fluid due to shear stress as progression of time.. 7 1.5 the difference of power fluids . . . . . . . . . . . . . . . . . . . . . 9 1.6 Nitrogen (left) and Argon (right) viscosity as a function of the tem- perature and pressure after Lemmon and Jacobsen. .......... 10 1.7 The shear stress as a function of the shear rate . . . . . . . . . . . . 11 1.8 Air viscosity as a function of the temperature. . . . . . . . . . . . . 12 1.9 Water viscosity as a function temperature.. . . . . . . . . . . . . . 12 1.10 Liquid metals viscosity as a function of the temperature. . . . . . . 13 1.11 Reduced viscosity as function of the reduced temperature. . . . . . 17 1.12 Reduced viscosity as function of the reduced temperature. . . . . . 18 1.13 Surface Tension control volume analysis . . . . . . . . . . . . . . . 21 1.14 Forces in Contact angle.. . . . . . . . . . . . . . . . . . . . . . . . 24 1.15 Description of wetting and non-wetting fluids. . . . . . . . . . . . . 24 1.16 Description of liquid surface.. . . . . . . . . . . . . . . . . . . . . . 26 1.17 The raising height as a function of the radii . . . . . . . . . . . . . . 29 1.18 The raising height as a function of the radius. . . . . . . . . . . . . 29 3.1 Description of how the center of mass is calculated. . . . . . . . 42 3.2 Thin body center of mass/area schematic.. . . . . . . . . . . . . . 42 3.3 The schematic that explains the summation of moment of inertia. . . 44 3.4 The schematic to explain the summation of moment of inertia.. . . . 45 3.5 Cylinder with the element for calculation moment of inertia. . . . . . 45 3.6 Description of rectangular in x-y plane for calculation of moment of inertia . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4 6 v  vi LIST OF FIGURES 3.7 A square element for the calculations of inertia of two-dimensional to three-dimensional deviations. . . . . . . . . . . . . . . . . . . . . 47 3.8 The ratio of the moment of inertia of two-dimensional to three-dimensional. 47 3.9 Description of parabola for calculation of moment of inertia and cen- ter of area... . . .................................48 3.10 Product of inertia for triangle. . . . . . . . . . . . . . . . . . . . . . . 49 4.1 Description of a fluid element in accelerated system under body forces. 55 4.2 Pressure lines a static fluid with a constant density.. . . . . . . . . . 58 4.3 A schematic to explain the measure of the atmospheric pressure. . . 58 4.4 Schematic of gas measurement utilizing the "U" tube. . . . . . . . . 59 4.5 Schematic of sensitive measurement device. . . . . . . . . . . . . . 60 4.6 Hydrostatic pressure when there is compressibility in the liquid phase. 64 4.7 Two adjoin layers for stability analysis. . . . . . . . . . . . . . . . . . 67 4.8 The varying gravity effects on density and pressure. . . . . . . . . . 69 4.9 The effective gravity is for accelerated cart. . . . . . . . . . . . . . . 73 4.10 A cart slide on inclined plane . . . . . . . . . . . . . . . . . . . . . . 73 4.11 Forces diagram of cart sliding on inclined plane . . . . . . . . . . . . 74 4.12 Schematic to explain the angular angle. . . . . . . . . . . . . . . . . 74 4.13 Rectangular area under pressure. . . . . . . . . . . . . . . . . . . . 75 4.14 Schematic of submerged area to explain the center forces and mo- m ents. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77 4.15 The general forces acting on submerged area. . . . . . . . . . . . . 78 4.16 The general forces acting on non symmetrical straight area. . . . . . 79 4.17 The general forces acting on non symmetrical straight area. . . . . . 80 4.18 The effects of multi layers density on static forces. . . . . . . . . . . 83 4.19 The forces on curved area... . ........................85 4.20 Schematic of Net Force on floating body.. . . . . . . . . . . . . . . . 86 4.21 Dam is a part of a circular shape.. . . . . . . . . . . . . . . . . . . . 87 4.22 Area above the dam arc subtract triangle. . . . . . . . . . . . . . . . 87 4.23 Area above the dam arc calculation for the center. . . . . . . . . . . 88 4.24 Moment on arc element around Point "0." . . . . . . . . . . . . . . . 89 4.25 Polynomial shape dam description for the moment around point "0" and force calculations. . . . . . . . . . . . . . . . . . . . . . . . . . . 90 4.26 The difference between the slop and the direction angle.. . . . . . . 90 4.27 Schematic of Immersed Cylinder. . . . . . . . . . . . . . . . . . . . . 92 4.28 The floating forces on Immersed Cylinder. . . . . . . . . . . . . . . . 93 4.29 Schematic of a thin wall floating body. . . . . . . . . . . . . . . . . . 94 4.30 Schematic of floating bodies. . . . . . . . . . . . . . . . . . . . . . . 98 4.31 Schematic of floating cubic. . . . . . . . . . . . . . . . . . . . . . . . 98 4.32 Stability analysis of floating body. . . . . . . . . . . . . . . . . . . . . 99 4.33 Cubic body dimensions for stability analysis.. . . . . . . . . . . . . . 100 4.34 Stability of cubic body infinity long. . . . . . . . . . . . . . . . . . . . 101 4.35 The maximum height reverse as a function of density ratio. . . . . . 102  LIST OF FIGURES VII 4.36 The effects of liquid movement on the GM. . . . . . . . . . . . . . . 103 4.37 Measurement of GM of floating body.. . . . . . . . . . . . . . . . . . 105 4.38 Calculations of GM for abrupt shape body. . . . . . . . . . . . . . . 106 4.39 A heavy needle is floating on a liquid. . . . . . . . . . . . . . . . . . 107 4.40 Description of depression to explain the Rayleigh-Taylor instability. . 108 4.41 Description of depression to explain the instability. . . . . . . . . . . 110 4.42 The cross section of the interface. The purple color represents the maximum heavy liquid raising area. The yellow color represents the maximum lighter liquid that are "going down." . . . . . . . . . . . . . 111 5.1 lines a static fluid with a constant density. . . . . . . . . . . . . . . . 115 5.2 Stratified flow in horizontal tubes when the liquids flow is very slow. . 117 5.3 Kind of Stratified flow in horizontal tubes. . . . . . . . . . . . . . . . 118 5.4 Plug flow in horizontal tubes when the liquids flow is faster slow. . . 118 5.5 Modified Mandhane map for flow regime in horizontal tubes.. . . . . 119 5.6 lines a static fluid with a constant density. . . . . . . . . . . . . . . . 120 5.7 A dimensional vertical flow map under very low gravity against the gravity.. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 121 5.8 The terminal velocity that left the solid particles.. . . . . . . . . . . . 131 5.9 The flow patterns in solid-liquid flow. . . . . . . . . . . . . . . . . . . 132 5.10 Counter-current flow in a can (the left figure) has only one hole thus pulse flow and a flow with two holes (right picture). . . . . . . . . . . 134 5.11 Counter-flow in vertical tubes map.. . . . . . . . . . . . . . . . . . . 134 5.12 Pictures of Counter-current flow in liquid-gas and solid-gas config- urations. The container is made of two compartments. The upper compartment is filled with the heavy phase (liquid, water solution, or small wood particles) by rotating the container. Even though the solid-gas ratio is smaller, it can be noticed that the solid-gas is faster than the liquid-gas flow. . . . . . . . . . . . . . . . . . . . . 135 5.13 Flood in vertical pipe. . . . . . . . . . . . . . . . . . . . . . . . . . . 135 5.14 A flow map to explain the horizontal counter-current flow. . . . . . . 136 5.15 A diagram to explain the flood in a two dimension geometry.. . . . . 137 5.16 General forces diagram to calculated the in a two dimension geometry.142  viii LIST OF FIGURES  LIST OF TABLES 1 Books Under Potto Project . . . . . . . . . . . . . . . . . . . . . . . xxix 1 continue .. . ................................ xxx 1.1 Sutherland's equation coefficients . . . . . . . . . . . . . . . . . . 13 1.2 Viscosity of selected gases . . . . . . . . . . . . . . . . . . . . . 14 1.3 Viscosity of selected liquids . . . . . . . . . . . . . . . . . . . . . 14 1.4 Properties at the critical stage . . . . . . . . . . . . . . . . . . . . . 15 1.5 Bulk modulus for selected materials . . . . . . . . . . . . . . . . . 20 1.6 The contact angle for air/water with selected materials.. . . . . . . 25 1.7 The surface tension for selected materials . . . . . . . . . . . . . . 31 1.7 continue.. ..................................32 2.1 Properties of Various Ideal Gases [300K] . . . . . . . . . . . . . . 38 3.1 Moments of Inertia for various plane surfaces about their center of gravity (full shapes) . . . . . . . . . . . . . . . . . . . . . . . . . . . 53 3.2 Moment of inertia for various plane surfaces about their center of gravity . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54 ix  x LIST OF TABLES  NOMENCLATURE R Universal gas constant, see equation (2.26), page 38 £ Units length., see equation (2.1), page 33 viscosity at input temperature T, see equation (1.17), page 12 po reference viscosity at reference temperature, To, see equation (1.17), page 12 Martinelli parameter, see equation (5.43), page 129 A The area of surface, see equation (4.117), page 85 a The acceleration of object or system, see equation (4.0), page 55 Bf Body force, see equation (2.9), page 35 C, Specific pressure heat, see equation (2.23), page 37 C, Specific volume heat, see equation (2.22), page 37 Eu Internal energy, see equation (2.3), page 34 Eu Internal Energy per unit mass, see equation (2.6), page 34 E System energy at state i, see equation (2.2), page 34 G The gravitation constant, see equation (4.62), page 70 9G general Body force, see equation (4.0), page 55 H Enthalpy, see equation (2.18), page 36 h Specific enthalpy, see equation (2.18), page 36 xi  xii LIST OF TABLES k the ratio of the specific heats, see equation (2.24), page 37 L Angular momentum, see equation (3.38), page 51 Patmos Atmospheric Pressure, see equation (4.85), page 77 q Energy per unit mass, see equation (2.6), page 34 Q12 The energy transfered to the system between state 1 and state 2, see equa- tion (2.2), page 34 R Specific gas constant, see equation (2.27), page 38 S Entropy of the system, see equation (2.13), page 36 Suth Suth is Sutherland's constant and it is presented in the Table 1.1, see equa- tion (1.17), page 12 TT Torque, see equation (3.40), page 51 T2o reference temperature in degrees Kelvin, see equation (1.17), page 12 Tin input temperature in degrees Kelvin, see equation (1.17), page 12 U velocity , see equation (2.4), page 34 w Work per unit mass, see equation (2.6), page 34 W12 The work done by the system between state 1 and state 2, see equa- tion (2.2), page 34 z the coordinate in z direction, see equation (4.14), page 58  The Book Change Log Version 0.1.8 Aug 6, 2008 (2.6 M 189 pages) - Add the chapter on introduction to muli-phase flow - Again additional improvement to the index (thanks to Irene). - Add the Rayleigh-Taylor instability. - Improve the doChap scrip to break up the book to chapters. Version 0.1.6 Jun 30, 2008 (1.3 M 151 pages) - Fix the English in the introduction chapter, (thanks to Tousher). - Improve the Index (thanks to Irene). - Remove the multiphase chapter (it is not for public consumption yet). Version 0.1.5a Jun 11, 2008 (1.4 M 155 pages) - Add the constant table list for the introduction chapter. - Fix minor issues (English) in the introduction chapter. xiii  xiv LIST OF TABLES Version 0.1.5 Jun 5, 2008 (1.4 M 149 pages) - Add the introduction, viscosity and other properties of fluid. - Fix very minor issues (English) in the static chapter. Version 0.1.1 May 8, 2008 (1.1 M111 pages) - Major English corrections for the three chapters. - Add the product of inertia to mechanics chapter. - Minor corrections for all three chapters. Version 0.1a April 23, 2008 Version 0.1a April 23, 2008 - The Thermodynamics chapter was released. - The mechanics chapter was released. - The static chapter was released (the most extensive and detailed chapter).  Notice of Copyright For This Document: This document published Modified FDL. The change of the license is to prevent from situations where the author has to buy his own book. The Potto Project Li- cense isn't long apply to this document and associated docoments. GNU Free Documentation License The modification is that under section 3 "copying in quantity" should be add in the end. "If you print more than 200 copies, you are required to furnish the author with two (2) copies of the printed book." Version 1.2, November 2002 Copyright ©2000,2001,2002 Free Software Foundation, Inc. 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA Everyone is permitted to copy and distribute verbatim copies of this license document, but changing it is not allowed. Preamble The purpose of this License is to make a manual, textbook, or other functional and useful document "free" in the sense of freedom: to assure everyone the effective freedom to copy and redistribute it, with or without modifying it, either xv  xvi LIST OF TABLES commercially or noncommercially. Secondarily, this License preserves for the au- thor and publisher a way to get credit for their work, while not being considered responsible for modifications made by others. This License is a kind of "copyleft", which means that derivative works of the document must themselves be free in the same sense. It complements the GNU General Public License, which is a copyleft license designed for free software. We have designed this License in order to use it for manuals for free software, because free software needs free documentation: a free program should come with manuals providing the same freedoms that the software does. But this License is not limited to software manuals; it can be used for any textual work, regardless of subject matter or whether it is published as a printed book. We recommend this License principally for works whose purpose is instruction or ref- erence. 1. APPLICABILITY AND DEFINITIONS This License applies to any manual or other work, in any medium, that contains a notice placed by the copyright holder saying it can be distributed under the terms of this License. Such a notice grants a world-wide, royalty-free license, unlimited in duration, to use that work under the conditions stated herein. The "Document", below, refers to any such manual or work. Any member of the public is a licensee, and is addressed as "you". You accept the license if you copy, modify or distribute the work in a way requiring permission under copyright law. A "Modified Version" of the Document means any work containing the Document or a portion of it, either copied verbatim, or with modifications and/or translated into another language. A "Secondary Section" is a named appendix or a front-matter section of the Document that deals exclusively with the relationship of the publishers or authors of the Document to the Document's overall subject (or to related matters) and contains nothing that could fall directly within that overall subject. (Thus, if the Document is in part a textbook of mathematics, a Secondary Section may not explain any mathematics.) The relationship could be a matter of historical connec- tion with the subject or with related matters, or of legal, commercial, philosophical, ethical or political position regarding them. The "Invariant Sections" are certain Secondary Sections whose titles are designated, as being those of Invariant Sections, in the notice that says that the Document is released under this License. If a section does not fit the above definition of Secondary then it is not allowed to be designated as Invariant. The Document may contain zero Invariant Sections. If the Document does not identify any Invariant Sections then there are none. The "Cover Texts" are certain short passages of text that are listed, as Front-Cover Texts or Back-Cover Texts, in the notice that says that the Document is released under this License. A Front-Cover Text may be at most 5 words, and a Back-Cover Text may be at most 25 words.  GNU FREE DOCUMENTATION LICENSE xvi i A "Transparent" copy of the Document means a machine-readable copy, represented in a format whose specification is available to the general public, that is suitable for revising the document straightforwardly with generic text editors or (for images composed of pixels) generic paint programs or (for drawings) some widely available drawing editor, and that is suitable for input to text formatters or for automatic translation to a variety of formats suitable for input to text formatters. A copy made in an otherwise Transparent file format whose markup, or absence of markup, has been arranged to thwart or discourage subsequent modification by readers is not Transparent. An image format is not Transparent if used for any substantial amount of text. A copy that is not "Transparent" is called "Opaque". Examples of suitable formats for Transparent copies include plain ASCII without markup, Texinfo input format, LaTeX input format, SGML or XML using a publicly available DTD, and standard-conforming simple HTML, PostScript or PDF designed for human modification. Examples of transparent image formats include PNG, XCF and JPG. Opaque formats include proprietary formats that can be read and edited only by proprietary word processors, SGML or XML for which the DTD and/or processing tools are not generally available, and the machine-generated HTML, PostScript or PDF produced by some word processors for output purposes only. The "Title Page" means, for a printed book, the title page itself, plus such following pages as are needed to hold, legibly, the material this License re- quires to appear in the title page. For works in formats which do not have any title page as such, "Title Page" means the text near the most prominent appearance of the work's title, preceding the beginning of the body of the text. A section "Entitled XYZ" means a named subunit of the Document whose title either is precisely XYZ or contains XYZ in parentheses following text that translates XYZ in another language. (Here XYZ stands for a specific section name mentioned below, such as "Acknowledgements", "Dedications", "En- dorsements", or "History".) To "Preserve the Title" of such a section when you modify the Document means that it remains a section "Entitled XYZ" according to this definition. The Document may include Warranty Disclaimers next to the notice which states that this License applies to the Document. These Warranty Disclaimers are considered to be included by reference in this License, but only as regards dis- claiming warranties: any other implication that these Warranty Disclaimers may have is void and has no effect on the meaning of this License. 2. VERBATIM COPYING You may copy and distribute the Document in any medium, either com- mercially or noncom mercially, provided that this License, the copyright notices, and the license notice saying this License applies to the Document are reproduced in all copies, and that you add no other conditions whatsoever to those of this Li- cense. You may not use technical measures to obstruct or control the reading or further copying of the copies you make or distribute. However, you may accept  xviii LIST OF TABLES compensation in exchange for copies. If you distribute a large enough number of copies you must also follow the conditions in section 3. You may also lend copies, under the same conditions stated above, and you may publicly display copies. 3. COPYING IN QUANTITY If you publish printed copies (or copies in media that commonly have printed covers) of the Document, numbering more than 100, and the Document's license notice requires Cover Texts, you must enclose the copies in covers that carry, clearly and legibly, all these Cover Texts: Front-Cover Texts on the front cover, and Back-Cover Texts on the back cover. Both covers must also clearly and legibly identify you as the publisher of these copies. The front cover must present the full title with all words of the title equally prominent and visible. You may add other material on the covers in addition. Copying with changes limited to the covers, as long as they preserve the title of the Document and satisfy these conditions, can be treated as verbatim copying in other respects. If the required texts for either cover are too voluminous to fit legibly, you should put the first ones listed (as many as fit reasonably) on the actual cover, and continue the rest onto adjacent pages. If you publish or distribute Opaque copies of the Document number- ing more than 100, you must either include a machine-readable Transparent copy along with each Opaque copy, or state in or with each Opaque copy a computer- network location from which the general network-using public has access to down- load using public-standard network protocols a complete Transparent copy of the Document, free of added material. If you use the latter option, you must take rea- sonably prudent steps, when you begin distribution of Opaque copies in quantity, to ensure that this Transparent copy will remain thus accessible at the stated location until at least one year after the last time you distribute an Opaque copy (directly or through your agents or retailers) of that edition to the public. It is requested, but not required, that you contact the authors of the Doc- ument well before redistributing any large number of copies, to give them a chance to provide you with an updated version of the Document. 4. MODIFICATIONS You may copy and distribute a Modified Version of the Document under the conditions of sections 2 and 3 above, provided that you release the Modified Version under precisely this License, with the Modified Version filling the role of the Document, thus licensing distribution and modification of the Modified Version to whoever possesses a copy of it. In addition, you must do these things in the Modified Version: A. Use in the Title Page (and on the covers, if any) a title distinct from that of the Document, and from those of previous versions (which should, if there were any, be listed in the History section of the Document). You may use the  GNU FREE DOCUMENTATION LICENSE xix same title as a previous version if the original publisher of that version gives permission. B. List on the Title Page, as authors, one or more persons or entities responsible for authorship of the modifications in the Modified Version, together with at least five of the principal authors of the Document (all of its principal authors, if it has fewer than five), unless they release you from this requirement. C. State on the Title page the name of the publisher of the Modified Version, as the publisher. D. Preserve all the copyright notices of the Document. E. Add an appropriate copyright notice for your modifications adjacent to the other copyright notices. F. Include, immediately after the copyright notices, a license notice giving the public permission to use the Modified Version under the terms of this License, in the form shown in the Addendum below. G. Preserve in that license notice the full lists of Invariant Sections and required Cover Texts given in the Document's license notice. H. Include an unaltered copy of this License. I. Preserve the section Entitled "History", Preserve its Title, and add to it an item stating at least the title, year, new authors, and publisher of the Modified Version as given on the Title Page. If there is no section Entitled "History" in the Document, create one stating the title, year, authors, and publisher of the Document as given on its Title Page, then add an item describing the Modified Version as stated in the previous sentence. J. Preserve the network location, if any, given in the Document for public access to a Transparent copy of the Document, and likewise the network locations given in the Document for previous versions it was based on. These may be placed in the "History" section. You may omit a network location for a work that was published at least four years before the Document itself, or if the original publisher of the version it refers to gives permission. K. For any section Entitled "Acknowledgements" or "Dedications", Preserve the Title of the section, and preserve in the section all the substance and tone of each of the contributor acknowledgements and/or dedications given therein. L. Preserve all the Invariant Sections of the Document, unaltered in their text and in their titles. Section numbers or the equivalent are not considered part of the section titles. M. Delete any section Entitled "Endorsements". Such a section may not be in- cluded in the Modified Version.  xx LIST OF TABLES N. Do not retitle any existing section to be Entitled "Endorsements" or to conflict in title with any Invariant Section. 0. Preserve any Warranty Disclaimers. If the Modified Version includes new front-matter sections or appendices that qualify as Secondary Sections and contain no material copied from the Docu- ment, you may at your option designate some or all of these sections as invariant. To do this, add their titles to the list of Invariant Sections in the Modified Version's license notice. These titles must be distinct from any other section titles. You may add a section Entitled "Endorsements", provided it contains nothing but endorsements of your Modified Version by various parties-for example, statements of peer review or that the text has been approved by an organization as the authoritative definition of a standard. You may add a passage of up to five words as a Front-Cover Text, and a passage of up to 25 words as a Back-Cover Text, to the end of the list of Cover Texts in the Modified Version. Only one passage of Front-Cover Text and one of Back- Cover Text may be added by (or through arrangements made by) any one entity. If the Document already includes a cover text for the same cover, previously added by you or by arrangement made by the same entity you are acting on behalf of, you may not add another; but you may replace the old one, on explicit permission from the previous publisher that added the old one. The author(s) and publisher(s) of the Document do not by this License give permission to use their names for publicity for or to assert or imply endorse- ment of any Modified Version. 5. COMBINING DOCUMENTS You may combine the Document with other documents released under this License, under the terms defined in section 4 above for modified versions, provided that you include in the combination all of the Invariant Sections of all of the original documents, unmodified, and list them all as Invariant Sections of your combined work in its license notice, and that you preserve all their Warranty Disclaimers. The combined work need only contain one copy of this License, and multiple identical Invariant Sections may be replaced with a single copy. If there are multiple Invariant Sections with the same name but different contents, make the title of each such section unique by adding at the end of it, in parentheses, the name of the original author or publisher of that section if known, or else a unique number. Make the same adjustment to the section titles in the list of Invariant Sections in the license notice of the combined work. In the combination, you must combine any sections Entitled "History" in the various original documents, forming one section Entitled "History"; likewise combine any sections Entitled "Acknowledgements", and any sections Entitled "Dedications". You must delete all sections Entitled "Endorsements".  GNU FREE DOCUMENTATION LICENSE xxi 6. COLLECTIONS OF DOCUMENTS You may make a collection consisting of the Document and other doc- uments released under this License, and replace the individual copies of this Li- cense in the various documents with a single copy that is included in the collection, provided that you follow the rules of this License for verbatim copying of each of the documents in all other respects. You may extract a single document from such a collection, and distribute it individually under this License, provided you insert a copy of this License into the extracted document, and follow this License in all other respects regarding verbatim copying of that document. 7. AGGREGATION WITH INDEPENDENT WORKS A compilation of the Document or its derivatives with other separate and independent documents or works, in or on a volume of a storage or distribution medium, is called an "aggregate" if the copyright resulting from the compilation is not used to limit the legal rights of the compilation's users beyond what the individual works permit. When the Document is included in an aggregate, this Li- cense does not apply to the other works in the aggregate which are not themselves derivative works of the Document. If the Cover Text requirement of section 3 is applicable to these copies of the Document, then if the Document is less than one half of the entire aggregate, the Document's Cover Texts may be placed on covers that bracket the Document within the aggregate, or the electronic equivalent of covers if the Document is in electronic form. Otherwise they must appear on printed covers that bracket the whole aggregate. 8. TRANSLATION Translation is considered a kind of modification, so you may distribute translations of the Document under the terms of section 4. Replacing Invariant Sections with translations requires special permission from their copyright holders, but you may include translations of some or all Invariant Sections in addition to the original versions of these Invariant Sections. You may include a translation of this License, and all the license notices in the Document, and any Warranty Dis- claimers, provided that you also include the original English version of this License and the original versions of those notices and disclaimers. In case of a disagree- ment between the translation and the original version of this License or a notice or disclaimer, the original version will prevail. If a section in the Document is Entitled "Acknowledgements", "Dedica- tions", or "History", the requirement (section 4) to Preserve its Title (section 1) will typically require changing the actual title. 9. TERMINATION  xxii LIST OF TABLES You may not copy, modify, sublicense, or distribute the Document except as expressly provided for under this License. Any other attempt to copy, modify, sublicense or distribute the Document is void, and will automatically terminate your rights under this License. However, parties who have received copies, or rights, from you under this License will not have their licenses terminated so long as such parties remain in full compliance. 10. FUTURE REVISIONS OF THIS LICENSE The Free Software Foundation may publish new, revised versions of the GNU Free Documentation License from time to time. Such new versions will be similar in spirit to the present version, but may differ in detail to address new prob- lems or concerns. See http://www.gnu.org/copyleft/. Each version of the License is given a distinguishing version number. If the Document specifies that a particular numbered version of this License "or any later version" applies to it, you have the option of following the terms and conditions either of that specified version or of any later version that has been published (not as a draft) by the Free Software Foundation. If the Document does not specify a version number of this License, you may choose any version ever published (not as a draft) by the Free Software Foundation. ADDENDUM: How to use this License for your documents To use this License in a document you have written, include a copy of the License in the document and put the following copyright and license notices just after the title page: Copyright ©YEAR YOUR NAME. Permission is granted to copy, dis- tribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or any later version published by the Free Software Foundation; with no Invariant Sections, no Front- Cover Texts, and no Back-Cover Texts. A copy of the license is included in the section entitled "GNU Free Documentation License". If you have Invariant Sections, Front-Cover Texts and Back-Cover Texts, replace the "with...Texts." line with this: with the Invariant Sections being LIST THEIR TITLES, with the Front- Cover Texts being LIST, and with the Back-Cover Texts being LIST. If you have Invariant Sections without Cover Texts, or some other com- bination of the three, merge those two alternatives to suit the situation. If your document contains nontrivial examples of program code, we rec- ommend releasing these examples in parallel under your choice of free software license, such as the GNU General Public License, to permit their use in free soft- ware.  CONTRIBUTOR LIST How to contribute to this book As a copylefted work, this book is open to revisions and expansions by any inter- ested parties. The only "catch" is that credit must be given where credit is due. This is a copyrighted work: it is not in the public domain! If you wish to cite portions of this book in a work of your own, you must follow the same guidelines as for any other GDL copyrighted work. Credits All entries have been arranged in alphabetical order of surname (hopefully. Major contributions are listed by individual name with some detail on the nature of the contribution(s), date, contact info, etc. Minor contributions (typo corrections, etc.) are listed by name only for reasons of brevity. Please understand that when I classify a contribution as "minor," it is in no way inferior to the effort or value of a "major" contribution, just smaller in the sense of less text changed. Any and all contributions are gratefully accepted. I am indebted to all those who have given freely of their own knowledge, time, and resources to make this a better book! - Date(s) of contribution(s): 1999 to present - Nature of contribution: Original author. - Contact at: barmeir at gmail. com Steven from artofproblemsolving.com - Date(s) of contribution(s): June 2005 xxiii  xxiv LIST OF TABLES - Nature of contribution: LaTeX formatting, help on building the useful equa- tion and important equation macros. Dan Olsen - Date(s) of contribution(s): April 2008 - Nature of contribution: Some discussions about chapter on mechanics and correction of English. Richard Hackbarth - Date(s) of contribution(s): April 2008 - Nature of contribution: Some discussions about chapter on mechanics and correction of English. Tousher Yang - Date(s) of contribution(s): April 2008 - Nature of contribution: Major review of statics and thermo chapters. Your name here - Date(s) of contribution(s): Month and year of contribution - Nature of contribution: Insert text here, describing how you contributed to the book. - Contact at: myemail@provider .net Typo corrections and other "minor" contributions - R. Gupta, January 2008, help with the original img macro and other ( LaTeX issues).  About This Author Genick Bar-Meir holds a Ph.D. in Mechanical Engineering from University of Min- nesota and a Master in Fluid Mechanics from Tel Aviv University. Dr. Bar-Meir was the last student of the late Dr. R.G.E. Eckert. Much of his time has been spend doing research in the field of heat and mass transfer (related to renewal energy issues) and this includes fluid mechanics related to manufacturing processes and design. Currently, he spends time writing books (there are already two very popu- lar books) and softwares for the POTTO project (see Potto Prologue). The author enjoys to encourage his students to understand the material beyond the basic re- quirements of exams. In his early part of his professional life, Bar-Meir was mainly interested in elegant models whether they have or not a practical applicability. Now, this author's views had changed and the virtue of the practical part of any model becomes the essential part of his ideas, books and software. He developed models for Mass Transfer in high concentration that be- came a building blocks for many other models. These models are based on an- alytical solution to a family of equations1. As the change in the view occurred, Bar-Meir developed models that explained several manufacturing processes such the rapid evacuation of gas from containers, the critical piston velocity in a par- tially filled chamber (related to hydraulic jump), application of supply and demand to rapid change power system and etc. All the models have practical applicability. These models have been extended by several research groups (needless to say with large research grants). For example, the Spanish Comision Interministerial provides grants TAP97-0489 and PB98-0007, and the CICYT and the European Commission provides 1 FD97-2333 grants for minor aspects of that models. More- over, the author's models were used in numerical works, in GM, British industry, Spain, and Canada. ' Where the mathematicians were able only to prove that the solution exists. xxv  xxvi LIST OF TABLES In the area of compressible flow, it was commonly believed and taught that there is only weak and strong shock and it is continue by Prandtl-Meyer func- tion. Bar-Meir discover the analytical solution for oblique shock and show that there is three solutions, and there is a gap between the oblique shock and Prandtl- Meyer. He also build analytical solution to several moving shock cases. He de- scribed and categorized the filling and evacuating of chamber by compressible fluid in which he also found analytical solutions to cases where the working fluid was ideal gas. The common explanation to Prandtl-Meyer function shows that flow can turn in a sharp corner. Engineers have constructed design that based on this conclusion. Bar-Meir demonstrated that common Prandtl-Meyer explanation vio- lates the conservation of mass and therefor the turn must be around a finite radius. The author's explanations on missing diameter and other issues in fanno flow and ""naughty professor's question"" are used in the industry. In his book "Basics of Fluid Mechanics", Bar-Meir demonstrated that flu- ids must have wavy surface when the materials flow together. All the previous models for the flooding phenomenon did not have a physical explanation to the dryness. He built a model to explain the flooding problem (two phase flow) based on the physics. He also build many categories for two flow regimes. The author lives with his wife and three children. A past project of his was building a four stories house, practically from scratch. While he writes his programs and does other computer chores, he often feels clueless about computers and programing. While he is known to look like he knows about many things, the author just know to learn quickly. The author spent years working on the sea (ships) as a engine sea officer but now the author prefers to remain on solid ground.  Prologue For The POTTO Project This books series was born out of frustrations in two respects. The first issue is the enormous price of college textbooks. It is unacceptable that the price of the college books will be over $150 per book (over 10 hours of work for an average student in The United States). The second issue that prompted the writing of this book is the fact that we as the public have to deal with a corrupted judicial system. As individuals we have to obey the law, particularly the copyright law with the "infinite2" time with the copyright holders. However, when applied to "small" individuals who are not able to hire a large legal firm, judges simply manufacture facts to make the little guy lose and pay for the defense of his work. On one hand, the corrupted court system defends the "big" guys and on the other hand, punishes the small "entrepreneur" who tries to defend his or her work. It has become very clear to the author and founder of the POTTO Project that this situation must be stopped. Hence, the creation of the POTTO Project. As R. Kook, one of this author's sages, said instead of whining about arrogance and incorrectness, one should increase wisdom. This project is to increase wisdom and humility. The POTTO Project has far greater goals than simply correcting an abu- sive Judicial system or simply exposing abusive judges. It is apparent that writing textbooks especially for college students as a cooperation, like an open source, is a new idea3. Writing a book in the technical field is not the same as writing a novel. The writing of a technical book is really a collection of information and prac- tice. There is always someone who can add to the book. The study of technical 2After the last decision of the Supreme Court in the case of Eldred v. Ashcroff (see http: //cyber. law. harvard. edu/openlaw/eldredvashcroft for more information) copyrights practically remain in- definitely with the holder (not the creator). 3|n some sense one can view the encyclopedia Wikipedia as an open content project (see http: //en. wikipedia. org/wiki/Main_Page). The wikipedia is an excellent collection of articles which are written by various individuals. xxvii  xxviii LIST OF TABLES material isn't only done by having to memorize the material, but also by coming to understand and be able to solve related problems. The author has not found any technique that is more useful for this purpose than practicing the solving of prob- lems and exercises. One can be successful when one solves as many problems as possible. To reach this possibility the collective book idea was created/adapted. While one can be as creative as possible, there are always others who can see new aspects of or add to the material. The collective material is much richer than any single person can create by himself. The following example explains this point: The army ant is a kind of carnivorous ant that lives and hunts in the tropics, hunting animals that are even up to a hundred kilograms in weight. The secret of the ants' power lies in their collective intelligence. While a single ant is not intelligent enough to attack and hunt large prey, the collective power of their networking creates an extremely powerful intelligence to carry out this attack ( see for information http://www.ex.ac.uk/ bugclub/raiders.html)4. When an insect which is blind can be so powerful by networking, So can we in creating textbooks by this powerful tool. Why would someone volunteer to be an author or organizer of such a book? This is the first question the undersigned was asked. The answer varies from individual to individual. It is hoped that because of the open nature of these books, they will become the most popular books and the most read books in their respected field. In a way, the popularity of the books should be one of the incentives for potential contributors. The desire to be an author of a well-known book (at least in his/her profession) will convince some to put forth the effort. For some authors, the reason is the pure fun of writing and organizing educational material. Experience has shown that in explaining to others any given subject, one also begins to better understand the material. Thus, contributing to these books will help one to understand the material better. For others, the writing of or contributing to this kind of books will serve as a social function. The social function can have at least two components. One component is to come to know and socialize with many in the profession. For others the social part is as simple as a desire to reduce the price of college textbooks, especially for family members or relatives and those students lacking funds. For some contributors/authors, in the course of their teaching they have found that the textbook they were using contains sections that can be improved or that are not as good as their own notes. In these cases, they now have an opportunity to put their notes to use for others. Whatever the reasons, the undersigned believes that personal intentions are appropriate and are the author's/organizer's private affair. If a contributor of a section in such a book can be easily identified, then that contributor will be the copyright holder of that specific section (even within question/answer sections). The book's contributor's names could be written by their sections. It is not just for experts to contribute, but also students who hap- pened to be doing their homework. The student's contributions can be done by 4see also in Franks, Nigel R.; "Army Ants: A Collective Intelligence," American Scientist, 77:139, 1989  CREDITS xxix adding a question and perhaps the solution. Thus, this method is expected to accelerate the creation of these high quality books. These books are written in a similar manner to the open source software process. Someone has to write the skeleton and hopefully others will add "flesh and skin." In this process, chapters or sections can be added after the skeleton has been written. It is also hoped that others will contribute to the question and answer sections in the book. But more than that, other books contain data5 which can be typeset in LATEX. These data (tables, graphs and etc.) can be redone by anyone who has the time to do it. Thus, the contributions to books can be done by many who are not experts. Additionally, contributions can be made from any part of the world by those who wish to translate the book. It is hoped that the books will be error-free. Nevertheless, some errors are possible and expected. Even if not complete, better discussions or better ex- planations are all welcome to these books. These books are intended to be "con- tinuous" in the sense that there will be someone who will maintain and improve the books with time (the organizer(s)). These books should be considered more as a project than to fit the tradi- tional definition of "plain" books. Thus, the traditional role of author will be replaced by an organizer who will be the one to compile the book. The organizer of the book in some instances will be the main author of the work, while in other cases only the gate keeper. This may merely be the person who decides what will go into the book and what will not (gate keeper). Unlike a regular book, these works will have a version number because they are alive and continuously evolving. The undersigned of this document intends to be the organizer-author- coordinator of the projects in the following areas: Table -1: Books under development in Potto project. Availability Project Progress Remarks Version for Name Public Download Compressible beta 0.4.8.2 / Flow Die Casting alpha 0.0.3 / Dynamics NSY 0.0.0 X Fluid Mechanics alpha 0.1.8 / Heat Transfer NSY Based 0.0.0 X on Eckert Mechanics NSY 0.0.0 X 5 Data are not copyrighted.  xxx LIST OF TABLES Table -1: Books under development in Potto project. (continue) Availability Project Progress Remarks Version for Name Public Download Open Channel NSY 0.0.0 X Flow Statics early first 0.0.1 X alpha chapter Strength of Mate- NSY 0.0.0 X rial Thermodynamics early 0.0.1 X alpha Two/Multi phases NSY Tel- 0.0.0 X flow Aviv'notes NSY = Not Started Yet The meaning of the progress is as: - The Alpha Stage is when some of the chapters are already in a rough draft; - in Beta Stage is when all or almost all of the chapters have been written and are at least in a draft stage; - in Gamma Stage is when all the chapters are written and some of the chap- ters are in a mature form; and - the Advanced Stage is when all of the basic material is written and all that is left are aspects that are active, advanced topics, and special cases. The mature stage of a chapter is when all or nearly all the sections are in a mature stage and have a mature bibliography as well as numerous examples for every section. The mature stage of a section is when all of the topics in the section are written, and all of the examples and data (tables, figures, etc.) are already pre- sented. While some terms are defined in a relatively clear fashion, other definitions give merely a hint on the status. But such a thing is hard to define and should be enough for this stage. The idea that a book can be created as a project has mushroomed from the open source software concept, but it has roots in the way science progresses. However, traditionally books have been improved by the same author(s), a process in which books have a new version every a few years. There are book(s) that have continued after their author passed away, i.e., the Boundary Layer Theory originated6 by Hermann Schlichting but continues to this day. However, projects 6Originally authored by Dr. Schlichting, who passed way some years ago. A new version is created every several years.  CREDITS xxxi such as the Linux Documentation project demonstrated that books can be written as the cooperative effort of many individuals, many of whom volunteered to help. Writing a textbook is comprised of many aspects, which include the ac- tual writing of the text, writing examples, creating diagrams and figures, and writing the LATEX macros7 which will put the text into an attractive format. These chores can be done independently from each other and by more than one individual. Again, because of the open nature of this project, pieces of material and data can be used by different books. 7One can only expect that open source and readable format will be used for this project. But more than that, only LATEX, and perhaps troff, have the ability to produce the quality that one expects for these writings. The text processes, especially LATEX, are the only ones which have a cross platform ability to produce macros and a uniform feel and quality. Word processors, such as OpenOffice, Abiword, and Microsoft Word software, are not appropriate for these projects. Further, any text that is produced by Microsoft and kept in "Microsoft" format are against the spirit of this project In that they force spending money on Microsoft software.  xxxii LIST OF TABLES  Prologue For This Book Version 0.1.8 August 6, 2008 pages 189 size 2.6M When this author was an undergraduate student, he spend time to study the wave phenomenon at the interface of open channel flow. This issue is related to renewal energy of extracting energy from brine solution (think about the Dead Sea, so much energy). The common explanation to the wave existence was that there is always a disturbance which causes instability. This author was bothered by this explanation. Now, in this version, it was proven that this wavy interface is created due to the need to satisfy the continuous velocity and shear stress at the interface and not a disturbance. Potto project books are characterized by high quality which marked by presentation of the new developments and clear explanations. This explanation (on the wavy interface) demonstrates this characteristic of Potto project books. The introduction to multi-phase is another example to this quality. While it is a hard work to discover and develop and bring this information to the students, it is very satisfying for the author. The number of downloads of this book results from this quality. Even in this early development stage, number of downloads per month is about 5000 copies. Version 0.1 April 22, 2008 pages 151 size 1.3M The topic of fluid mechanics is common to several disciplines: mechanical engi- neering, aerospace engineering, chemical engineering, and civil engineering. In fact, it is also related to disciplines like industrial engineering, and electrical en- xxxiii  xxxiv LIST OF TABLES gineering. While the emphasis is somewhat different in this book, the common material is presented and hopefully can be used by all. One can only admire the wonderful advances done by the previous geniuses who work in this field. In this book it is hoped to insert, what and when a certain model is suitable than other models. One of the difference in this book is the insertion of the introduction to multiphase flow. Clearly, multiphase is an advance topic. However, some minimal familiarity can be helpful for many engineers who have to deal with non pure single phase fluid. This book is the third book in the series of POTTO project books. POTTO project books are open content textbooks so everyone are welcome to joint in. The topic of fluid mechanics was chosen just to fill the introduction chapter to com- pressible flow. During the writing it became apparent that it should be a book in its own right. In writing the chapter on fluid statics, there was a realization that it is the best chapter written on this topic. It is hoped that the other chapters will be as good this one. This book is written in the spirit of my adviser and mentor E.R.G. Eckert. Eckert, aside from his research activity, wrote the book that brought a revolution in the education of the heat transfer. Up to Egret's book, the study of heat transfer was without any dimensional analysis. He wrote his book because he realized that the dimensional analysis utilized by him and his adviser (for the post doc), Ernst Schmidt, and their colleagues, must be taught in engineering classes. His book met strong criticism in which some called to "burn" his book. Today, however, there is no known place in world that does not teach according to Eckert's doctrine. It is assumed that the same kind of individual(s) who criticized Eckert's work will criti- cize this work. Indeed, the previous book, on compressible flow, met its opposition. For example, anonymous Wikipedia user name EMBaero claimed that the material in the book is plagiarizing, he just doesn't know from where and what. Maybe that was the reason that he felt that is okay to plagiarize the book on Wikipedia. These criticisms will not change the future or the success of the ideas in this work. As a wise person says "don't tell me that it is wrong, show me what is wrong"; this is the only reply. With all the above, it must be emphasized that this book is not expected to revolutionize the field but change some of the way things are taught. The book is organized into several chapters which, as a traditional text- book, deals with a basic introduction to the fluid properties and concepts (under construction). The second chapter deals with Thermodynamics. The third book chapter is a review of mechanics. The next topic is statics. When the Static Chap- ter was written, this author did not realize that so many new ideas will be inserted into this topic. As traditional texts in this field, ideal flow will be presented with the issues of added mass and added forces (under construction). The classic issue of turbulence (and stability) will be presented. An introduction to multi-phase flow, not a traditional topic, will be presented next (again under construction). The next two chapters will deals with open channel flow and gas dynamics. At this stage, dimensional analysis will be present (again under construction).  How This Book Was Written This book started because I needed an introduction to the compressible flow book. After a while it seems that is easier to write a whole book than the two original planned chapters. In writing this book, it was assumed that introductory book on fluid mechanics should not contained many new ideas but should be modern in the material presentation. There are numerous books on fluid mechanics but none of which is open content. The approach adapted in this book is practical, and more hands-on approach. This statement really meant that the book is intent to be used by students to solve their exams and also used by practitioners when they search for solutions for practical problems. So, issue of proofs so and so are here only either to explain a point or have a solution of exams. Otherwise, this book avoids this kind of issues. The structure of Hansen, Streeter and Wylie, and Shames books were adapted and used as a scaffolding for this book. This author was influenced by Streeter and Wylie book which was his undergrad textbooks. The chapters are not written in order. The first 4 chapters were written first because they were sup- posed to be modified and used as fluid mechanics introduction in "Fundamentals of Compressible Flow." Later, multi-phase flow chapter was written. The presentation of some of the chapters is slightly different from other books because the usability of the computers. The book does not provide the old style graphical solution methods yet provides the graphical explanation of things. Of course, this book was written on Linux (Micro$oftLess book). This book was written using the vim editor for editing (sorry never was able to be com- fortable with emacs). The graphics were done by TGIF, the best graphic program that this author experienced so far. The figures were done by gle. The spell check- ing was done by ispell, and hope to find a way to use gaspell, a program that currently cannot be used on new Linux systems. The figure in cover page was created by Genick Bar-Meir, and is copyleft by him. xxxv  xxxvi LIST OF TABLES  Preface "In the beginning, the POTTO project was without form, and void; and emptiness was upon the face of the bits and files. And the Fingers of the Author moved upon the face of the keyboard. And the Author said, Let there be words, and there were words. "$. This book, Basics of Fluid Mechanics, describes the fundamentals of fluid mechanics phenomena for engineers and others. This book is designed to replace all introductory textbook(s) or instructor's notes for the fluid mechanics in undergraduate classes for engineering/science students but also for technical peo- ples. It is hoped that the book could be used as a reference book for people who have at least some basics knowledge of science areas such as calculus, physics, etc. The structure of this book is such that many of the chapters could be usable independently. For example, if you need information about, say, statics' equations, you can read just chapter (4). I hope this makes the book easier to use as a reference manual. However, this manuscript is first and foremost a textbook, and secondly a reference manual only as a lucky coincidence. I have tried to describe why the theories are the way they are, rather than just listing "seven easy steps" for each task. This means that a lot of information is presented which is not necessary for everyone. These explanations have been marked as such and can be skipped.9 Reading everything will, naturally, increase your understanding of the many aspects of fluid mechanics. This book is written and maintained on a volunteer basis. Like all volun- teer work, there is a limit on how much effort I was able to put into the book and its organization. Moreover, due to the fact that English is my third language and time limitations, the explanations are not as good as if I had a few years to perfect them. Nevertheless, I believe professionals working in many engineering fields will benefit from this information. This book contains many worked examples, which can be very useful for many. 8To the power and glory of the mighty God. This book is only to explain his power. 9At the present, the book is not well organized. You have to remember that this book is a work in progress. xxxvii  xxxvii i LIST OF TABLES I have left some issues which have unsatisfactory explanations in the book, marked with a Mata mark. I hope to improve or to add to these areas in the near future. Furthermore, I hope that many others will participate of this project and will contribute to this book (even small contributions such as providing examples or editing mistakes are needed). I have tried to make this text of the highest quality possible and am in- terested in your comments and ideas on how to make it better. Incorrect language, errors, ideas for new areas to cover, rewritten sections, more fundamental material, more mathematics (or less mathematics); I am interested in it all. I am particularly interested in the best arrangement of the book. If you want to be involved in the editing, graphic design, or proofreading, please drop me a line. You may contact me via Email at "barmeir@gmail.com". Naturally, this book contains material that never was published before (sorry cannot avoid it). This material never went through a close content review. While close content peer review and publication in a professional publication is excellent idea in theory. In practice, this process leaves a large room to blockage of novel ideas and plagiarism. If you would like be "peer reviews" or critic to my new ideas please send me your comment(s). Even reaction/comments from individuals like David Marshall10. Several people have helped me with this book, directly or indirectly. I would like to especially thank to my adviser, Dr. E. R. G. Eckert, whose work was the inspiration for this book. I also would like to thank to Jannie McRotien (Open Channel Flow chapter) and Tousher Yang for their advices, ideas, and assistance. The symbol META was added to provide typographical conventions to blurb as needed. This is mostly for the author's purposes and also for your amuse- ment. There are also notes in the margin, but those are solely for the author's pur- poses, ignore them please. They will be removed gradually as the version number advances. I encourage anyone with a penchant for writing, editing, graphic ability, LATEX knowledge, and material knowledge and a desire to provide open content textbooks and to improve them to join me in this project. If you have Internet e-mail access, you can contact me at "barmeir@gmail.com". 10Dr. Marshall wrote to this author that the author should review other people work before he write any thing new (well, literature review is always good, isn't it?). Over ten individuals wrote me about this letter. I am asking from everyone to assume that his reaction was innocent one. While his comment looks like unpleasant reaction, it brought or cause the expansion of the explanation for the oblique shock. However, other email that imply that someone will take care of this author aren't appreciated.  To Do List and Road Map This book isn't complete and probably never will be completed. There will always new problems to add or to polish the explanations or include more new materials. Also issues that associated with the book like the software has to be improved. It is hoped the changes in TEX and LATEX related to this book in future will be min- imal and minor. It is hoped that the style file will be converged to the final form rapidly. Nevertheless, there are specific issues which are on the "table" and they are described herein. At this stage, many chapters are missing. Specific missing parts from every chapters are discussed below. These omissions, mistakes, approach prob- lems are sometime appears in the book under the Meta simple like this Meta sample this part. Meta End You are always welcome to add a new material: problem, question, illustration or photo of experiment. Material can be further illuminate. Additional material can be provided to give a different angle on the issue at hand. Properties The chapter isn't in development stage yet. Open Channel Flow The chapter isn't in the development stage yet. Some parts were taken from Fun- damentals of Die Casting Design book and are in a process of improvement. xxxix  xl LIST OF TABLES  CHAPTER 1 Introduction 1.1 What is Fluid Mechanics? Fluid mechanics deals with the study of all fluids under static and dynamic situa- tions. Fluid mechanics is a branch of continuous mechanics which deals with a relationship between forces, motions, and statical conditions in continuous mate- rial. This study area deals with many and diversified problems such as surface tension, fluid statics, flow in enclose bodies, or flow round bodies (solid or other- wise), flow stability, etc. In fact, almost any action a person is doing involves some kind of a fluid mechanics problem. Furthermore, the boundary between the solid mechanics and fluid mechanics is some kind of gray shed and not a sharp distinc- tion (see Figure (1.1) for the complex relationships between the different branches which only part of it should be drawn in the same time.). For example, glass ap- pears as a solid material, but a closer look reveals that the glass is a liquid with a large viscosity. A proof of the glass "liquidity" is the change of the glass thickness in high windows in European Churches after hundred years. The bottom part of the glass is thicker than the top part. Materials like sand (some call it quick sand) and grains should be treated as liquids. It is known that these materials have the ability to drown people. Even material such as aluminum just below the mushy zone also behaves as a liquid similarly to butter. After it was established that the boundaries of fluid mechanics aren't sharp, the discussion in this book is limited to simple and (mostly) Newtonian (sometimes power fluids) fluids which will be defined later. The fluid mechanics study involve many fields that have no clear boundary between them. Researchers distinguish between orderly flow and chaotic flow as the laminar flow and the turbulent flow. The fluid mechanics can also be distinguish between a single phase flow and multiphase flow (flow made more than one phase or single distinguishable material). The last boundary (as all the boundaries in fluid 1  2 CHAPTER 1. INTRODUCTION W"- Turbulent Flow Laminar Flow Fig. -1.1: Diagram to explain part of relationships of fluid mechanics branches. mechanics) isn't sharp because fluid can go through a phase change (condensa- tion or evaporation) in the middle or during the flow and switch from a single phase flow to a multi phase flow. Moreover, flow with two phases (or materials) can be treated as a single phase (for example, air with dust particle). After it was made clear that the boundaries of fluid mechanics aren't sharp, the study must make arbitrary boundaries between fields. Then the dimensional analysis will be used explain why in certain cases one distinguish area/principle is more relevant than the other and some effects can be neglected. Or, when a gen- eral model is need because more parameters are effecting the situation. It is this author's personal experience that the knowledge and ability to know in what area the situation lay is one of the main problems. For example, engineers in software company (EKK Inc, http: //ekkinc. com/HTML) analyzed a flow of a complete still  1.2. BRIEF HISTORY 3 liquid assuming a complex turbulent flow model. Such absurd analysis are com- mon among engineers who do not know which model can be applied. Thus, one of the main goals of this book is to explain what model should be applied. Before dealing with the boundaries, the simplified private cases must be explained. There are two main approaches of presenting an introduction of fluid me- chanics materials. The first approach introduces the fluid kinematic and then the basic governing equations, to be followed by stability, turbulence, boundary layer and internal and external flow . The second approach deals with the Integral Anal- ysis to be followed with Differential Analysis, and continue with Empirical Analysis. These two approaches pose a dilemma to anyone who writes an introductory book for the fluid mechanics. These two approaches have justifications and positive points. Reviewing many books on fluid mechanics made it clear, there isn't a clear winner. This book attempts to find a hybrid approach in which the kinematic is presented first (aside to standard initial four chapters) follow by Integral analysis and continued by Differential analysis. The ideal flow (frictionless flow) should be expanded compared to the regular treatment. This book is unique in providing chapter on multiphase flow. Naturally, chapters on open channel flow (as a sub class of the multiphase flow) and compressible flow (with the latest developments) are provided. 1.2 Brief History The need to have some understanding of fluid mechanics started with the need to obtain water supply. For example, people realized that wells have to be dug and crude pumping devices need to be constructed. Later, a large population created a need to solve waste (sewage) and some basic understanding was created. At some point, people realized that water can be used to move things and provide power. When cities increased to a larger size, aqueducts were constructed. These aqueducts reached their greatest size and grandeur in those of the City of Rome and China. Yet, almost all knowledge of the ancients can be summarized as appli- cation of instincts, with the exception Archimedes (250 B.C.) on the principles of buoyancy. For example, larger tunnels built for a larger water supply, etc. There were no calculations even with the great need for water supply and transportation. The first progress in fluid mechanics was made by Leonardo Da Vinci (1452-1519) who built the first chambered canal lock near Milan. He also made several at- tempts to study the flight (birds) and developed some concepts on the origin of the forces. After his initial work, the knowledge of fluid mechanics (hydraulic) in- creasingly gained speed by the contributions of Galileo, Torricelli, Euler, Newton, Bernoulli family, and D'Alembert. At that stage theory and experiments had some discrepancy. This fact was acknowledged by D'Alembert who stated that, "The theory of fluids must necessarily be based upon experiment." For example the concept of ideal liquid that leads to motion with no resistance, conflicts with the reality.  4 CHAPTER 1. INTRODUCTION This discrepancy between theory and practice is called the "D'Alembert paradox" and serves to demonstrate the limitations of theory alone in solving fluid problems. As in thermodynamics, two different of school of thoughts were cre- ated: the first believed that the solution will come from theoretical aspect alone, and the second believed that solution is the pure practical (experimental) aspect of fluid mechanics. On the theoretical side, considerable contribution were made by Euler, La Grange, Helmhoitz, Kirchhoff, Rayleigh, Rankine, and Kelvin. On the "ex- perimental" side, mainly in pipes and open channels area, were Brahms, Bossut, Chezy, Dubuat, Fabre, Coulomb, Dupuit, d'Aubisson, Hagen, and Poisseuille. In the middle of the nineteen century, first Navier in the molecular level and later Stokes from continuous point of view succeeded in creating governing equa- tions for real fluid motion. Thus, creating a matching between the two school of thoughts: experimental and theoretical. But, as in thermodynamics, people cannot relinquish control. As results it created today "strange" names: Hydrodynamics, Hydraulics, Gas Dynamics, and Aeronautics. The Navier-Stokes equations, which describes the flow (or even Euler equations), were considered unsolvable during the mid nineteen century because of the high complexity. This problem led to two consequences. Theoreticians tried to simplify the equations and arrive at approximated solutions representing specific cases. Examples of such work are Hermann von Helmholtz's concept of vortexes (1858), Lanchester's concept of circulatory flow (1894), and the Kutta-Joukowski circulation theory of lift (1906). The experimentalists, at the same time proposed many correlations to many fluid mechanics problems, for example, resistance by Darcy, Weisbach, Fanning, Ganguillet, and Manning. The obvious happened without theoretical guidance, the empirical formulas generated by fitting curves to experimental data (even sometime merely presenting the results in tabular form) resulting in formulas that the relationship between the physics and properties made very little sense. At the end of the twenty century, the demand for vigorous scientific knowl- edge that can be applied to various liquids as opposed to formula for every fluid was created by the expansion of many industries. This demand coupled with new sev- eral novel concepts like the theoretical and experimental researches of Reynolds, the development of dimensional analysis by Rayleigh, and Froude's idea of the use of models change the science of the fluid mechanics. Perhaps the most rad- ical concept that effects the fluid mechanics is of Prandtl's idea of boundary layer which is a combination of the modeling and dimensional analysis that leads to modern fluid mechanics. Therefore, many call Prandtl as the father of modern fluid mechanics. This concept leads to mathematical basis for many approxima- tions. Thus, Prandtl and his students Blasius, von Karman, Meyer, and Blasius and several other individuals as Nikuradse, Rose, Taylor, Bhuckingham, Stanton, and many others, transformed the fluid mechanics to modern science that we have known today. While the understanding of the fundamentals did not change much, af- ter World War Two, the way how it was calculated changed. The introduction of  1.3. KINDS OF FLUIDS 5 the computers during the 60s and much more powerful personal computer has changed the field. There are many open source programs that can analyze many fluid mechanics situations. Today many problems can be analyzed by using the numerical tools and provide reasonable results. These programs in many cases can capture all the appropriate parameters and adequately provide a reasonable description of the physics. However, there are many other cases that numerical analysis cannot provide any meaningful result (trends). For example, no weather prediction program can produce good engineering quality results (where the snow will fall within 50 kilometers accuracy. Building a car with this accuracy is a dis- aster). In the best scenario, these programs are as good as the input provided. Thus, assuming turbulent flow for still flow simply provides erroneous results (see for example, EKK, Inc). 1.3 Kinds of Fluids Some differentiate fluid from solid by the reaction to shear stress. It is a known fact said that the fluid continuously and permanently deformed under shear stress while solid exhibits a finite deformation which does not change with time. It is also said that liquid cannot return to their original state after the deformation. This dif- ferentiation leads to three groups of materials: solids and liquids. This test creates a new material group that shows dual behaviors; under certain limits; it behaves like solid and under others it behaves like liquid (see Figure 1.1). The study of this kind of material called rheology and it will (almost) not be discussed in this book. It is evident from this discussion that when a liquid is at rest, no shear stress is applied. The fluid is mainly divided into two categories: liquids and gases. The main difference between the liquids and gases state is that gas will occupy the whole volume while liquids has an almost fix volume. This difference can be, for most practical purposes considered, sharp even though in reality this difference isn't sharp. The difference between a gas phase to a liquid phase above the crit- ical point are practically minor. But below the critical point, the change of water pressure by 1000% only change the volume by less than 1 percent. For example, a change in the volume by more 5% will required tens of thousands percent change of the pressure. So, if the change of pressure is significantly less than that, then the change of volume is at best 5%. Hence, the pressure will not affect the volume. In gaseous phase, any change in pressure directly affects the volume. The gas fills the volume and liquid cannot. Gas has no free interface/surface (since it does fill the entire volume). There are several quantities that have to be addressed in this discussion. The first is force which was reviewed in physics. The unit used to measure is [N]. It must be remember that force is a vector, e.g it has a direction. The second quantity discussed here is the area. This quantity was discussed in physics class but here it has an addtional meaning, and it is referred to the direction of the area. The direction of area is perpendicular to the area. The area is measured in [m2]. Area  6 CHAPTER 1. INTRODUCTION of three-dimensional object has no single direction. Thus, these kinds of areas should be addressed infinitesimally and locally. The traditional quantity, which is force per area has a new meaning. This is a result of division of a vector by a vector and it is referred to as tensor. In this book, the emphasis is on the physics, so at this stage the tensor will have to be broken into its components. Later, the discussion on the mathematic meaning will be presented (later version). For the discussion here, the pressure has three components, one in the area direction and two perpendicular to the area. The pressure component in the area direction is called pressure (great way to confuse, isn't it?). The other two components are referred as the shear stresses. The units used for the pressure components is [N/rM2]. The density is a property which requires that liquid to be continuous. The density can be changed and it is a function of time and space (location) but must have a continues property. It doesn't mean that a sharp and abrupt change in the density cannot occur. It referred to density that is independent of the sampling size. Figure 1.2 shows the density as a function of the sam- 1OgE ple size. After certain sample size, the density remains constant. Thus, the density is defined as .F p = rn Am AV-E AV Fig. -1.2: Density as a function of the size of sample. (1.1) It must be noted that E is chosen so that the continuous assumption is not broken, that is, it did not reach/reduced to the size where the atoms or molecular statistical calculations are significant (see Figure 1.2 for point where the green lines con- verge to constant density). When this assumption is broken, then, the principles of statistical mechanics must be utilized. 1.4 Shear Stress The shear stress is part of the pressure ten- sor. However, here it will be treated as a separate issue. In solid mechanics, the shear stress is considered as the ratio of the force acting on area in the direction of the forces perpendicular to area. Different from solid, fluid cannot pull directly but through a solid surface. Consider liquid that under- goes a shear stress between a short dis- tance of two plates as shown in Figure (1.3). The upper plate velocity generally will Uo0 F h Fig. -1.3: Schematics to describe the shear stress in fluid mechanics. U= f(A,F,h) (1.2)  1.4. SHEAR STRESS 7 Where A is the area, the F denotes the force, h is the distance between the plates. From solid mechanics study, it was shown that when the force per area increases, the velocity of the plate increases also. Experiments show that the increase of height will increase the velocity up to a certain range. Consider moving the plate with a zero lubricant (h ~ 0) (results in large force) or a large amount of lubricant (smaller force). In this discussion, the aim is to develop differential equation, thus the small distance analysis is applicable. For cases where the dependency is linear, the following can be written hF U oA (1.3) Equations (1.3) can be rearranged to be U F h a A (1.4) Shear stress was defined as F (1.5) From equations (1.4) and (1.5) it follows that ratio of the velocity to height is pro- portional to shear stress. Hence, applying the coefficient to obtain a new equality as U Txy h (1.6) Where p is called the absolute viscosity or dynamic viscosity which will be dis- cussed later in this chapter in great length. In steady state, the distance the-to <"t2 TO (1.14) Shear--ate___- _ dU = 0 if Tyl 0 (2.11) The integration symbol with the circle represent integral of cycle (therefor circle) in with system return to the same condition. If there is no lost, it is referred as a reversible process and the inequality change to equality. 1> 0 (2.12) The last integral can go though several states. These states are independent of the path the system goes through. Hence, the integral is independent of the path. This observation leads to the definition of entropy and designated as S and the derivative of entropy is (6Q ds T (2.13) rev  36 CHAPTER 2. REVIEW OF THERMODYNAMICS Performing integration between two states results in S2 -1 j2()rev= dS (2.14) One of the conclusions that can be drawn from this analysis is for reversible and adiabatic process dS = 0. Thus, the process in which it is reversible and adi- abatic, the entropy remains constant and referred to as isentropic process. It can be noted that there is a possibility that a process can be irreversible and the right amount of heat transfer to have zero change entropy change. Thus, the reverse conclusion that zero change of entropy leads to reversible process, isn't correct. For reversible process equation (2.12) can be written as SQ = TdS (2.15) and the work that the system is doing on the surroundings is SW = PdV (2.16) Substituting equations (2.15) (2.16) into (2.10) results in TdS = d Eu + PdV (2.17) Even though the derivation of the above equations were done assuming that there is no change of kinetic or potential energy, it still remail valid for all situa- tions. Furthermore, it can be shown that it is valid for reversible and irreversible processes. Enthalpy It is a common practice to define a new property, which is the combination of already defined properties, the enthalpy of the system. H=Eu+PV (2.18) The specific enthalpy is enthalpy per unit mass and denoted as, h. Or in a differential form as dH = dEu + dPV + PdV (2.19) Combining equations (2.18) the (2.17) yields (2.20)  2.1. BASIC DEFINITIONS 37 For isentropic process, equation (2.17) is reduced to dH = VdP. The equation (2.17) in mass unit is Tds = du+ Pdv = dh - dP(2.21) p when the density enters through the relationship of p = 1/v. Specific Heats The change of internal energy and enthalpy requires new definitions. The first change of the internal energy and it is defined as the following (2.22) And since the change of the enthalpy involve some kind of work is defined as (2.23) The ratio between the specific pressure heat and the specific volume heat is called the ratio of the specific heat and it is denoted as, k. k --(2 .2 4 ) For solid, the ratio of the specific heats is almost 1 and therefore the differ- ence between them is almost zero. Commonly the difference for solid is ignored and both are assumed to be the same and therefore referred as C. This approx- imation less strong for liquid but not by that much and in most cases it applied to the calculations. The ratio the specific heat of gases is larger than one. Equation of state Equation of state is a relation between state variables. Normally the relation- ship of temperature, pressure, and specific volume define the equation of state for gases. The simplest equation of state referred to as ideal gas. and it is defined as P = pRT (2.25) Application of Avogadro's law, that "all gases at the same pressures and tempera- tures have the same number of molecules per unit of volume," allows the calculation of a "universal gas constant." This constant to match the standard units results in R = 8.3145 kmo K (2.26)  38 CHAPTER 2. REVIEW OF THERMODYNAMICS Thus, the specific gas can be calculate as R R M (2.27) The specific constants for select gas at 300K is provided in table 2.1. Table -2.1: Properties of Various Ideal Gases [300K] Gas Chemical Molecular RkR l C ] [Ceki] k Formula Weight [KgKj (KgKj [KgKj Air - 28.970 0.28700 1.0035 0.7165 1.400 Argon Ar 39.948 0.20813 0.5203 0.3122 1.400 Butane C4H10 58.124 0.14304 1.7164 1.5734 1.091 Carbon CO2 44.01 0.18892 0.8418 0.6529 1.289 Dioxide Carbon CO 28.01 0.29683 1.0413 0.7445 1.400 Monoxide Ethane C2H6 30.07 0.27650 1.7662 1.4897 1.186 Ethylene C2H4 28.054 0.29637 1.5482 1.2518 1.237 Helium He 4.003 2.07703 5.1926 3.1156 1.667 Hydrogen H2 2.016 4.12418 14.2091 10.0849 1.409 Methane CH4 16.04 0.51835 2.2537 1.7354 1.299 Neon Ne 20.183 0.41195 1.0299 0.6179 1.667 Nitrogen N2 28.013 0.29680 1.0416 0.7448 1.400 Octane 08H18 114.230 0.07279 1.7113 1.6385 1.044 Oxygen 02 31.999 0.25983 0.9216 0.6618 1.393 Propane C3H8 44.097 0.18855 1.6794 1.4909 1.327 Steam H20 18.015 0.48152 1.8723 1.4108 1.327 From equation (2.25) of state for perfect gas it follows d(Pv) = RdT For perfect gas dh = dEu + d(Pv) = dEu + d(RT) = f(T) (only) (2.28) (2.29)  2.1. BASIC DEFINITIONS 39 From the definition of enthalpy it follows that d(Pv) = dh - dEu (2.30) Utilizing equation (2.28) and subsisting into equation (2.30) and dividing by dT yields Cp - C = R (2.31 ) This relationship is valid only for ideal/perfect gases. The ratio of the specific heats can be expressed in several forms as (2.32) U--(2.33) The specific heat ratio, k value ranges from unity to about 1.667. These values depend on the molecular degrees of freedom (more explanation can be obtained in Van Wylen "F. of Classical thermodynamics." The values of several gases can be approximated as ideal gas and are provided in Table (2.1). The entropy for ideal gas can be simplified as the following 82 -s81 i( )(2.34) Using the identities developed so far one can find that lf dT f2 R dP T2 P2 82 - 81 ]CPT P = C1ln - R In (2.35) Or using specific heat ratio equation (2.35) transformed into s2_-__ __k ln-- In- (2.36) R k-1 T1 P1 For isentropic process, As = 0, the following is obtained IlnlyT2 (n )2 k(2.37) Tit Pl There are several famous identities that results from equation (2.37) as (2.38)  40 CHAPTER 2. REVIEW OF THERMODYNAMICS The ideal gas model is a simplified version of the real behavior of real gas. The real gas has a correction factor to account for the deviations from the ideal gas model. This correction factor referred as the compressibility factor and defined as Z RT(2.39)  CHAPTER 3 Review of Mechanics This author would like to express his gratitute to Dan Olsen (former Minneapolis city Engineer) and his friend Richard Hackbarth. This chapter provides a review of important definitions and concepts from Mechanics (statics and dynamics). These concepts and definitions will be used in this book and a review is needed. 3.1 Center of Mass The center of mass is divided into two sections, first, center of the mass and two, center of area (two-dimensional body with equal distribution mass). 3.1.1 Center of the Mass In many engineering problems, the center of mass is required to make the calcula- tions. This concept is derived from the fact that a body has a center of mass/gravity which interacts with other bodies and that this force acts on the center (equivalent force). It turns out that this concept is very useful in calculating rotations, moment of inertia, etc. The center of mass doesn't depend on the coordinate system and on the way it is calculated. The physical meaning of the center of mass is that if a straight line force acts on the body in away through the center of gravity, the body will not rotate. In other words, if a body will be held by one point it will be enough to hold the body in the direction of the center of mass. Note, if the body isn't be held through the center of mass, then a moment in addtional to force is required (to pre- vent the body for rotating). It is convenient to use the Cartesian system to explain this concept. Suppose that the body has a distribution of the mass (density, rho) as a function of the location. The density "normally" defined as mass per volume. Here, the the line density is referred to density mass per unit length in the x direc- tion. 41  42 CHAPTER 3. REVIEW OF MECHANICS In x coordinate, the center will be de- fined as m Iv dm x p(x)dV (3.1) Here, the dV element has finite dimen- dv sions in y-z plane and infinitesimal di- mension in x direction see Figure 3.1. Also, the mass, m is the total mass of the object. It can be noticed that center Fig. -3.1: Description of how the center of of mass in the x-direction isn't affectedF masscistcalculated. by the distribution in the y nor by z di- rections. In same fashion the center of mass can be defined in the other directions as following (3.2) where x is the direction of either, x, y or z. The density, p(xi) is the line density as function of xi. Thus, even for solid and uniform density the line density is a function of the geometry. 3.1.2 Center of Area In the previous case, the body was a three dimensional shape. There are cases where the body can be approxi- mated as a two-dimensional shape be- cause the body is with a thin with uni- form density. Consider a uniform thin body with constant thickness shown in Figure 3.2 which has density, p. Thus, equation (3.1) can be transfered into It y X dm 1 xp V Fig. -3.2: Thin body center of mass/area schematic. (3.3) The density, p and the thickness, t, are constant and can be canceled. Thus equation (3.3) can be transfered into (3.4)  3.2. MOMENT OF INERTIA 43 when the integral now over only the area as oppose over the volume. Finding the centeroid location should be done in the most convenient coordi- nate system since the location is coordinate independent. 3.2 Moment of Inertia As it was divided for the body center of mass, the moment of inertia is divided into moment of inertia of mass and area. 3.2.1 Moment of Inertia for Mass The moment of inertia turns out to be an essential part for the calculations of rotating bodies. Furthermore, it turns out that the moment of inertia has much wider applicability. Moment of inertia of mass is defined as Irrm = pr2dm (3.5) If the density is constant then equation (3.5) can be transformed into Irr= p r2dV (3.6) The moment of inertia is independent of the coordinate system used for the cal- culation, but dependent on the location of axis of rotation relative to the body. Some people define the radius of gyration as an equivalent concepts for the cen- ter of mass concept and which means if all the mass were to locate in the one point/distance and to obtain the same of moment of inertia. rk = (3.7) The body has a different moment of inertia for every coordinate/axis and they are Ixx = fy rx2dm = f(yy2 + z2) dm Iy =fy ry2dm = fy(x2 + z2) dm (3.8) Izz = fy rz2dm = f (x2 + y2) dm  44 CHAPTER 3. REVIEW OF MECHANICS 3.2.2 Moment of Inertia for Area 3.2.2.1 General Discussion For body with thickness, t and uniform density the following can be written moment of inertia for area Ixxm = LIr2dm pt jr2dA (3.9) The moment of inertia about axis is x can be defined as I. = r2dA = '"xm / Itl (3.10) where r is distance of dA from the Any point distance can be calculated from axis x as axis x and t is the thickness. x = y +z (3.11) Thus, equation (3.10) can be written as I (y2 + z2) dA (3.12) X' In the same fashion for other two coor- dinates as IYY (x2 + z2) dA Fig. -3.3: The schematic that explains the summation of moment of inertia. (3.13) Izz (x2 + y2) dA (3.14) 3.2.2.2 The Parallel Axis Theorem The moment of inertial can be calculated for any axis. The knowledge about one axis can help calculating the moment of inertia for a parallel axis. Let Ix the moment of inertia about axis xx which is at the center of mass/area. The moment of inertia for axis x' is Ixax _ ' r,2dAJfy,2+,2dA S[y+y)2+ (z+z)2] dA (3.15)  3.2. MOMENT OF INERTIA 45 equation (3.15) can be expended as I j(y2 + z2) dA+2j (yAy + z Az) dA+ A((y)2 + (Az)2) dA (3.16) The first term in equation (3.16) on the right hand side is the about axis x and the second them is zero. The second therm integral of center about center thus is zero. The third term is ar be written as moment of inertia is zero because it new term and can r2 A A constant ((AY)2 + (Az)2) dA ( 2 r2A (3.17) Hence, the relationship between the moment of inertia at xx and parallel axis x x is Igg = I.. + r2 A The moment of inertia of several areas is the sum of moment inertia of each area see Figure 3.4 and therefore, (3.18) Ixx = [Ixxi (2 i=1 If the same areas are similar thus n i 1 3.19) Fig. -3.4: The schematic to explain the sum- mation of moment of inertia. 3.20) Equation (3.20) is very useful in the calculation of the moment of inertia utilizing the moment of inertia of known bodies. For example, the moment of in- ertial of half a circle is half of whole cir- cle for axis a the center of circle. The moment of inertia can then move the center of area. of the Fig. -3.5: Cylinder with the element for calcu- lation moment of inertia.  46 CHAPTER 3. REVIEW OF MECHANICS 3.2.3 Examples of Moment of Inertia Example 3.1: Calculate the moment of inertia for the mass of the cylinder about center axis which height of h and radius, ro, as shown in Figure 3.5. The material is with an uniform density and homogeneous. SOLUTION The element can be calculated using cylindrical coordinate. Here element is a shell of thickness dr which shown in Figure 3.5 as the convenient Ir = p r2dm pj r2 h2wrrdr 4 2 r 2 The radius of gyration is mr02 m Example 3.2: Calculate the moment of inertia of the rectangular shape shown in Figure 3.6. SOLUTION ro The moment of inertia is calculated uti- lizing equation (3.12) as following 0 Ixx = y2 A +z2 dA dA a A 2 bdz 0 Fig. -3.6: Description of rectangular in x-y plane for calculation of moment of inertia. a3 b 3 Example 3.3: As measure to study the assumption of zero thickness consider a simple shape to see the effects of this assumption. Calculate the moment of inertia about the center of mass of a square shape with a thickness, t compare the results to a square shape with zero thickness. SOLUTION The moment of inertia of transverse slice about y' (see Figure mech:fig:squareEll)  3.2. MOMENT OF INERTIA 47 is A ba3 dlxxm = p dy 12 (3.21) The transformation into f x to center axis, x' can b lowing b a3 dI; ;i , = pdy 12 + rom local axis e done as fol- r2 A r2 A 11 (3.22) Fig. -3.7: A square element for the calcula- tions of inertia of two-dimensional to ined by integraionofqa rtion 2) to The totalr write as moment of inertia can be obta /t/2 ba3 I jm=p+z 2 b adz. -t/2 12 Comparison with the thin body results in I. pt ba3 Ixxm t2ba+ba3 a b t2 + a3 b pt 12 (3.23) 1 1+ t (3.24) It can be noticed right away that equation (3.24) indicates that ratio ap- proaches one when thickness ratio is approaches zero, Ixzm(t 0) 1. Ad- ditionally it can be noticed that the ratio a2/t2 is the only contributor to the er- ror1. The results are present in Figure 3.8. I can be noticed that the error is significant very fast even for small val- ues of t/a while the with of the box, b has no effect on the error. Fig. -3.8: The ratio of the moment of in- ertia of two-dimensional to three- dimensional. Example 3.4: Calculate the center of area and moment of inertia of the parabola, y =ax2, shown in Figure 3.9. 'This ratio is a dimensionless number that commonly has no special name. This author suggests to call this ratio as the B number.  48 CHAPTER 3. REVIEW OF MECHANICS SOLUTION For y = b the value of x = b/a. First the area inside the parabola calculated XX as dA/2 c A f bc a2(3a-1) b A= 3 ) a=2h-. The center of area can be calculated uti- lizing equation (3.4). The center of ev- Fig. -3.9: Description of parabola for calcula- ery element is at, (a2 + 2j 2 the el- tion of moment of inertia and center ement area is used before and therefore of area. 1 fb/t(2 (b - a 2 a 3 a b C= AL a + 2 ) (b-a dt=15a-5 (3.25) The moment of inertia of the area about the center can be found using in equation (3.25) can be done in two steps first calculate the moment of inertia in this coor- dinate system and then move the coordinate system to center. Utilizing equation (3.12) and doing the integration from 0 to maximum y provides dA Ib 2 2 b 72 Utilizing equation (3.18) I x/,/xA ( x x) I , AAX2_ 4b7/2 3a-1 (b (3ab 2 7 a 3 af) 15a - 5 or after working the details results in vb (20b3 - 14 b2) 35 a 3.2.4 Product of Inertia In addition to the moment of inertia, the product of inertia is commonly used. Here only the product of the area is defined and discussed. The product of inertia de- fined as I 2A z dA(3.26)  3.2. MOMENT OF INERTIA 49 For example, the product of inertia for x and y axises is I jy = x ydA (3.27) JA Product of inertia can be positive or negative value as oppose the moment of inertia. The calculation of the product of inertia isn't different much for the calcula- tion of the moment of inertia. The units of the product of inertia are the same as for moment of inertia. Transfer of Axis Theorem Same as for moment of inertia there is also similar theorem. Igjy/ = x' y'dA j(x + Ax) (y + Ay)dA (3.28) expanding equation (3.28) results in 0 0 syfxdA Az ydA Ixy JA JA AxAyA I jydA+ xAydA+JAxydA+AxAYdA (3.29) JAY JA JA JA The final form is Ig/Y/ = Ixy + AOxAy A (3.30) There are several relationships should be mentioned I y = Iyx (3.31) Symmetrical area has zero product of inertia. Example 3.5: y Calculate the product of inertia of straight edge triangle. SOLUTION The equation of the line is a y= x+a The product of inertia at the center is zero. The total Fig. -3.10:tiafor P triangle.r product of inertia is Ax Ay A a b (aba2b2 Ig/yg=0+ 3 3 2 18  50 CHAPTER 3. REVIEW OF MECHANICS 3.2.5 Principal Axes of Inertia The inertia matrix or inertia tensor is Ixx -Ixy -Ixz -Iyx Iyy -Iyz (3.32) -Izx -Izy Izz In linear algebra it was shown that for some angle equation (3.32) can be transform into Ixx 0 0 0 Ii'-o 0 (3.33) 0 0 I z'z' System which creates equation (3.33) referred as principle system. 3.3 Newton's Laws of Motion These laws can be summarized in two statements one, for every action by body A on Body B there is opposite reaction by body B on body A. Two, which can expressed in mathematical form as FD(mU)(3.34 Dt It can be noted that D replaces the traditional d since the addtional meaning which be added. Yet, it can be treated as the regular derivative. This law apply to any body and any body can "broken" into many small bodies which connected to each other. These small "bodies" when became small enough equation (3.34) can be transformed to a continuous form as D (pU)dv ( F =Jy Dt (3.35) The external forces are equal to internal forces the forces between the "small" bodies are cancel each other. Yet this examination provides a tool to study what happened in the fluid during operation of the forces. Since the derivative with respect to time is independent of the volume, the derivative can be taken out of the integral and the alternative form can be written as 5 F IpD UdV (3.36)  3.4. ANGULAR MOMENTUM AND TORQUE 51 The velocity, U is a derivative of the location with respect to time, thus, D2 SF =D2 JprdV (3.37) where r is the location of the particles from the origin. The external forces are typically divided into two categories: body forces and surface forces. The body forces are forces that act from a distance like magnetic field or gravity. The surface forces are forces that act on the surface of the body (pressure, stresses). The same as in the dynamic class, the system accelera- tion called the internal forces. The acceleration is divided into three categories: Centrifugal, wx(r x w), Angular, r x cZ, Coriolis, 2 (Ur x w). The radial velocity is denoted as Ur. 3.4 Angular Momentum and Torque The angular momentum of body, dm, is defined as L = r x Udm (3.38) The angular momentum of the entire system is calculated by integration (summa- tion) of all the particles in the system as L = r x U dm (3.39) The change with time of angular momentum is called torque, in analogous to the momentum change of time which is the force. D L D TT (r x Udm) _D( (3.40) where TT is the torque. The torque of entire system is T DL Df(rxUdm) (3.41) TS~mDt Dt i It can be noticed (well, it can be proved utilizing vector mechanics) that Dr D DDr D2r (3.42) TT Dt(r x U) Dt(r x Dt Dt2 To understand these equations a bit better, consider a particle moving in x-y plane. A force is acting on the particle in the same plane (x-y) plane. The velocity can be written as U = u+v) and the location from the origin can be written as r = xi+yj. The force can be written, in the same fashion, as F = Fxi + F j. Utilizing equation (3.38) provides L =rx U= x y 0 =(x v- yu)fc (3.43) U v 0  52 CHAPTER 3. REVIEW OF MECHANICS Utilizing equation (3.40) to calculate the torque as i j k TT =rxF (x y 0 =(F-yF)I(3.44) F F F 0 Since the torque is a derivative with respect to the time of the angular momentum it is also can be written as D xF, - yF = Dt [(xv - yu) dm] (3.45) The torque is a vector and the various components can be represented as Tx =i Dt m rx Udm (3.46) In the same way the component in y and z can be obtained. 3.4.1 Tables of geometries Th following tables present several moment of inertias of commonly used geome- tries.  3.4. ANGULAR MOMENTUM AND TORQUE 53 Table -3.1: Moments of Inertia for various plane surfaces about their center of gravity (full shapes) Shape Picture zc, yc A Ixz Name description XX b b CiRectangle V ab/2ba Elis -b/2 221 a XX b a a b ab ' Triangle T/ 3 36 a XX b b wb2nU Cicl b/2 2 4 64 a XX b b b wab6 Ab2 -b Z 2 2 4 64 1-5 6a-2 x V7 O3_14 b2) Parabola a  54 CHAPTER 3. REVIEW OF MECHANICS Table -3.2: Moment of inertia for various plane surfaces about their center of gravity Shape Picture Name description zc' yc A Ixx r Quadrant X 4---r(n of Circle4 3 r Ellipsoidal xx b 4 b a b Elisia 34b 37r 4 °bis Quadrant b a Half of .3 4ab Elliptic-3w 37 a Circular XX 0 2a r2 4 (-sin2) Sector 0 2r j-sn~ r xx 12r sin a 3 a I/X Circular Sector 4 (a+sin2I)  CHAPTER 4 Fluids Statics 4.1 Introduction The simplest situation that can occur in the study of fluid is when the fluid is at rest or queasy rest. This topic was introduced to most students in previous study of rigid body. However, here this topic will be more vigorously examined. Further- more, the student will be exposed to stability analysis probably for the first time. Later, the methods discussed here will be expanded to more complicated dynam- ics situations. 4.2 The Hydrostatic Equation _ , + dy) dxdz A fluid element with dimensions of dx, (adz dxdY dy, and dz is motionless in the accel- erated system, with acceleration, a as "±> (x+Pdx) dydz shown in Figure 4.1. The system is in a body force field, 9G(x,gy, z). The combination of an acceleration and the body force results in effective body force which is gG- a eff.1) Fig. -4.1: Description of a fluid element in ac- celerated system under body forces. Equation (4.1) can be reduced and simplified for the case of no acceleration, a = 0. In these derivations, several assumptions must be made. The first assump- tion is that the change in the pressure is a continuous function. There is no require- ment that the pressure has to be a monotonous function e.g. that pressure can 55  56 CHAPTER 4. FLUIDS STATICS increase and later decrease. The changes of the second derivative pressure are not significant compared to the first derivative (DP/&rn x di >>82P/&r2). where n is the steepest direction of the pressure derivative and di is the infinitesimal length. This mathematical statement simply requires that the pressure can deviate in such a way that the average on infinitesimal area can be found and expressed as only one direction. The net pressure force on the faces in the x direction results in dF = dydx i (4.2) In the same fashion, the calculations of the three directions result in the total net pressure force as 7OP - OP - F = ( i + j+ N (4.3) surface The term in the parentheses in equation (4.3) referred to in the literature as the pressure gradient. This mathematical operation has a geometrical interpreta- tion. If the pressure, P, was a two-dimensional height (that is only a function of x and y) then the gradient is the steepest ascent of the height (to the valley). The second point is that the gradient is a vector (that is, it has a direction). Even though, the pressure is treated, now, as a scalar function (there no reference to the shear stress in part of the pressure) the gradient is a vector. For example, the dot product of the following is i * gradP =(4.4) &x In general, if the coordinates were to "rotate/transform" to a new system which has a different orientation the dot product results in i gradP =(4.5) where i, is the unit vector in the n direction and O/On is a derivative in that direc- tion. As before, the effective gravity force is utilized in case where the gravity is the only body force and in an accelerated system. The body (element) is in rest and therefore the net force is zero 5F= F +(F (4.6) total surface body Hence, the utilizing the above derivations one can obtain -gradPdx dy dz + p genf dx dy dz = 0 (4.7)  4.3. PRESSURE AND DENSITY IN A GRAVITATIONAL FIELD 57 or grad e p af (4.8) Some refer to equation (4.8) as the Fluid Static Equation. This equation can be integrated and therefore solved. However, there are several physical implica- tions to this equation which should be discussed and are presented here. First, a discussion on a simple condition and will continue in more challenging situations. 4.3 Pressure and Density in a Gravitational Field In this section, a discussion on the pressure and the density in various conditions is presented. 4.3.1 Constant Density in Gravitational Field The simplest case is when the density, p, pressure, P, and temperature, T (in a way no function of the location) are constant. Traditionally, the z coordinate is used as the (negative) direction of the gravity1. The effective body force is ge =-g k (4.9) Utilizing equation (4.9) and substituting it into equation simple partial differential equations. These equations are (4.8) results into three OP Ox DP ay0 (4.10) and Prz Equations (4.10) can be integrated to yield (4.11) P(x, y) = constant (4.12) and constant in equation (4.12) can be absorbed by the (4.11) and therefore integration of equation P(x, y, z) = -pgz + constant (4.13) The integration constant is determined from the initial conditions or another 'This situation were the tradition is appropriated, it will be used. There are fields where x or y are designed to the direction of the gravity and opposite direction. For this reason sometime there will be a deviation from the above statement.  58 CHAPTER 4. FLUIDS STATICS Fig. -4.2: Pressure lines a static fluid with a constant density. point. For example, if at point comes zo the pressure is Po then the equation (4.13) be- P(z) - Po ---pg(z - z0) It is evident from equation (4.13) that the pressure depends only on z and/or the constant pressure lines are in the plane of x and y. Figure 4.2 describes the constant pressure lines in the con- tainer under the gravity body force. The pressure lines are continuous even in area where there is a discontinuous fluid. The reason that a solid boundary doesn't break the continuity of the pres- sure lines is because there is always a path to some of the planes. (4.14) pgh It is convenient to reverse the di- rection of z to get rid of the negative sign and to define h as the dependent of the fluid that is h ---(z - zo) so equation (4.14) becomes Fig. -4.3: A schematic to explain the measure of the atmospheric pressure. (4.15) In the literature, the right hand side of the equation (4.15) is defined as piezo- metric pressure.  4.3. PRESSURE AND DENSITY IN A GRAVITATIONAL FIELD 59 4.3.2 Pressure Measurement 4.3.2.1 Measuring the Atmospheric Pressure One of the application of this concept is the idea of measuring the atmospheric pressure. Consider a situation described in Figure 4.3. The liquid is filling the tube and is brought into a steady state. The pressure above the liquid on the right side is the vapor pressure. Using liquid with a very low vapor pressure like mercury, will result in a device that can measure the pressure without additional information (the temperature). Example 4.1: Calculate the atmospheric pressure at 20 C. The high of the Mercury is 0.76 [m] and the gravity acceleration is 9.82[m/sec]. Assume that the mercury vapor pressure is 0.0001 79264[kPa]. The description of the height is given in Figure 4.3. The mercury density is 13545.85[kg/m3]. SOLUTION The pressure is uniform or constant plane perpendicular to the gravity. Hence, knowing any point on this plane provides the pressure anywhere on the plane. The atmospheric pressure at point a is the same as the pressure on the right hand side of the tube. Equation (4.15) can be utilized and it can be noticed that pressure at point a is Pa = pgh+Pvapor (4.16) The density of the mercury is given along with the gravity and therefore, Pa = 13545.85 x 9.82 x 0.76 ~ 101095.39 [Pa] ~ 1.01 [Bar] The vapor pressure is about 1 x 10-4 percent of the total results. The main reason the mercury is used because its of large density and the fact that it is in a liquid phase in most of the mea- surement range. The third reason is the low vapor (partial) pressure of the mercury. The Gas 2 partial pressure of mercury is in the range of The pressure, P the 0.000001793[Bar] which is insignificant compared to the total measurement as can be observed from the above example. 1 4.3.2.2 Pressure Measurement The idea describes the atmospheric mea- Fig. -4.4: Schematic of gas measure- surement can be extended to measure the ment utilizing the "U" tube.  60 CHAPTER 4. FLUIDS STATICS pressure gas chambers. Consider a cham- ber filled with gas needed to be measured (see Figure 4.4). One technique is to attached "U" tube to the chamber and measure the pressure. This way, the gas is prevented from escaping and its pressure can be measured with a minimal inter- ference to the gas (some gas enters to the tube). The gas density is significantly lower than the liquid density and therefore can be neglected. The pressure at point "1" is P1 =Patmos+pgh (4.17) Since the atmospheric pressure was measured previously (the technique was shown in the previous section) the pressure of the chamber can be measured. 4.3.2.3 Magnified Pressure Measurement For situations where the pres- A, P, h4 A1 P sure difference is very small, en- gineers invented more sensitive measuring device. This device is build around the fact that the height is a function of the den-Ah2 sities difference. In the pre- vious technique, the density of one side was neglected (the gas side) compared to other side (liquid). This technique utilizes Fig. -4.5: Schematic of sensitive measurement de- the opposite range. The two vice. sides of the densities are very close to each other so the height became large. Figure 4.5 shows a typical and simple schematic of such an instrument. If the pressure differences between P1 and P2 is small this instrument can "magnified" height, h1 and provide "better" accuracy reading. This device is based on the fol- lowing mathematical explanation. In steady state, the pressure balance (only differences) is P1+gp1(hi+h2) =P2 +gh2p2 (4.18) It can be noticed that the "missing height" is canceled between the two sides. It can be noticed that h1 can be positive or negative or zero and it depends on the ratio that two containers filled with the light density liquid. Additionally, it can be observed that h1 is relatively small because A1 >> A2. The densities of the liquid are chosen so that they are close to each other but not equal. The densities of the liquids are chosen to be much heavier than the measured gas density. Thus, in writing equation (4.18) the gas density was neglected. The pressure difference can be expressed as PI - P2 =g [p2 h2 - pi(h + h2)] (4.19)  4.3. PRESSURE AND DENSITY IN A GRAVITATIONAL FIELD 61 If the light liquid volume in the two containers is known, it provides the relationship between h1 and h2. For example, if the volumes in two containers is equal then h2 A2 (.0 -h1 A1 =h2 A2 -0 h _21= A(4) A1 Liquid volumes do not necessarily have to be equal. Additional parameter, the volume ratio, will be introduced when the volumes ratio isn't equal. The calculations as results of this additional parameter does not cause significant complications. Here, this ratio equals to one and it simplify the equation (4.20). But this ratio easily can be inserted into the derivations. With the equation for height (4.20) equation (4.18) becomes P1 - gP2 gh2 P2 -p1 (1 - (4.21) or the height is h2 (4.22) 9 [(P2 - P1) + p1 A For the small value of the area ratio, A2/A1 «1, then equation (4.22) becomes h P2P) (4.23) g (p2 -p1) Some refer to the density difference shown in equation (4.23) as "magnification factor" since it replace the regular density, P2. 4.3.3 Varying Density in a Gravity Field There are several cases that will be discussed here which are categorized as gases, liquids and other. In the gas phase, the equation of state is simply the ideal gas model or the ideal gas with the compressibility factor (sometime referred to as real gas). The equation of state for liquid can be approximated or replaced by utilizing the bulk modulus. These relationships will be used to find the functionality between pressure, density and location. 4.3.3.1 Gas Phase under Hydrostatic Pressure Ideal Gas under Hydrostatic Pressure The gas density vary gradually with the pressure. As first approximation, the ideal gas model can be employed to describe the density. Thus equation (4.11) becomes -P RT (4.24) (9z RT  62 CHAPTER 4. FLUIDS STATICS Separating the variables and changing the partial derivatives to full derivative (just a notation for this case) results in dP_ gdz - - - (4.25) P R T Equation (4.25) can be integrated from point "0" to any point to yield P g InP - (z-z°) (4.26) P0 R T It is convenient to rearrange equation (4.26) to the following P g(z-zo) RTe (4.27) P0 Here the pressure ratio is related to the height exponentially. Equation (4.27) can be expanded to show the difference to standard assumption of constant pressure as -h P0 P0 P _(z - zo)g (z-zo)2g P0 1 RT + 6RT +(4.28) Po R T 6RT Or in a simplified form where the transformation of h = (z - zo) to be correction factor ++.--- (4.29) P0 P Equation (4.29) is useful in mathematical derivations but should be ignored for practical use2. Real Gas under Hydrostatic Pressure The mathematical derivations for ideal gas can be reused as a foundation for the real gas model (P = ZpRT). For a large range of P/Pc and T/Tc, the value of the compressibility factor, Z, can be assumed constant and therefore can be swallowed into equations (4.27) and (4.28). The compressibility is defined in equation (2.39). The modified equation is (g(z-zo)) e ZRT )(4.30) Po 2These derivations are left for a mathematical mind person. These deviations have a limited practical purpose. However, they are presented here for students who need to answer questions on this issue.  4.3. PRESSURE AND DENSITY IN A GRAVITATIONAL FIELD 63 Or in a series form which is (z-zo)g + (z-z)2 + (4.31) P0 ZRT 6ZRT Without going through the mathematics, the first approximation should be noticed that the compressibility factor, Z enter the equation as h/Z and not just h. Another point that is worth discussing is the relationship of Z to other gas properties. In general, the relationship is very complicated and in some ranges Z cannot be assumed constant. In these cases, a numerical integration must be carried out. 4.3.3.2 Liquid Phase Under Hydrostatic Pressure The bulk modulus was defined in equation (1.28). The simplest approach is to assume that the bulk modulus is constant (or has some representative average). For these cases, there are two differential equations that needed to be solved. Fortunately, here, only one hydrostatic equation depends on density equation. So, the differential equation for density should be solved first. The governing differential density equation (see equation (1.28)) is p=BT (4.32) The variables for equation (4.32) should be separated and then the integration can be carried out as JdP fj Bdp (4.33) Po Po P The integration of equation (4.33) yields P - Po = BT In-- (4.34) Po Equation (4.34) can be represented in a more convenient form as (4.35) Equation (4.35) is the counterpart for the equation of state of ideal gas for the liquid phase. Utilizing equation (4.35) in equation (4.11) transformed into OP_-gpo e BT (4.36) Oz Equation (4.36) can be integrated to yield P-Po BT e BT = z + Constant (4.37) g po  64 CHAPTER 4. FLUIDS STATICS It can be noted that BT has units of pressure and therefore the ratio in front of the exponent in equation (4.37) has units of length. The integration constant, with units of length, can be evaluated at any specific point. If at z = 0 the pressure is P0 and the density is po then the constant is Constant =BT (4.38) 9 po This constant, BT/g po, is a typical length of the problem. Additional discussion will be presented in the dimensionless issues chapter (currently under construction). The solution becomes -T - B - 1 = z (4.39) 9 Po Or in a dimensionless form (4.40) The solution is presented in equation (4.39) and is plotted in Figure 4.6. The solution is a reverse function (that is not P = f(z) but z = f (P)) it is a 4.0 monotonous function which is easy to3- solve for any numerical value (that is 30 only one z corresponds to any Pres- 2.5 sure). Sometimes, the solution is pre- 1.5 sented as- 5 - 1 15 20 25 30 - P BjIn(gpoz+1 + 1 (4.41) Mrh1, 208T Po Po BT An approximation of equation Fig. -4.6: Hydrostatic pressure when there is (4.40) is presented for historical rea- compressibility in the liquid phase. sons and in order to compare the constant density assumption. The exponent can be expanded as /piezometric corrections pressure BT P -Po)2 +BTP -Po)p (P -Po) + 2 B + 6 B +---. =z gpo (4.42) 2 T T B  4.3. PRESSURE AND DENSITY IN A GRAVITATIONAL FIELD 65 It can be noticed that equation (4.42) is reduced to the standard equation when the normalized pressure ratio, P/BT is small (<< 1). Additionally, it can be observed that the correction is on the left hand side and not as the "traditional" correction on the piezometric pressure side. 4.3.4 The Pressure Effects Because Temperature Variations 4.3.4.1 The Basic Analysis There are situations when the main change of the density results from other effects. For example, when the temperature field is not uniform, the density is affected and thus the pressure is a location function (for example, in the atmosphere the temperature is assumed to be a linear with the height under certain conditions.). A bit more complicate case is when the gas is a function of the pressure and another parameter. Air can be a function of the temperature field and the pressure. For the atmosphere, it is commonly assumed that the temperature is a linear function of the height. Here, a simple case is examined for which the temperature is a linear function of the height as dT -C.(4.43) dh where h here referred to height or distance. Hence, the temperature-distance function can be written as T = Constant - C, h (4.44) where the Constant is the integration constant which can be obtained by utilizing the initial condition. For h = 0, the temperature is To and using it leads to (4.45) Combining equation (4.45) with (4.11) results in 8P gP Oh -(T-C h)(4.46) 3h R (To - C h) Separating the variables in equation (4.46) and changing the formal 3 to the infor- mal d to obtain dP gdh P R(To-Cxh)  66 CHAPTER 4. FLUIDS STATICS Defining a new variable3 as = (To - C, h) for which o = To - C, ho and d/d= -C, d/dh. Using these definitions results in dP _ g d( P RCX (4.48) After integration of equation (4.47) and reusing (the reverse definitions) the vari- ables transformed the result into n = n To - Ch (449) Po R C, To Or in a more convenient form as P' 0 (4.50) It can be noticed that equation (4.50) is a monotonous function which decreases with height because the term in the brackets is less than one. This situation is roughly representing the pressure in the atmosphere and results in a temperature decrease. It can be observed that C, has a "double role" which can change the pressure ratio. Equation (4.50) can be approximated by two approaches/ideas. The first approximation for a small distance, h, and the second approximation for a small temperature gradient. It can be recalled that the following expansions are hPO correction factor PR. C 2g h (RgCa,-g2)h2 P= hm 1 -Cxh = 1- h-... g ) (4.51) Po h->O To To R 2 To 2R2 Equation (4.51) shows that the first two terms are the standard terms (negative sign is as expected i.e. negative direction). The correction factor occurs only at the third term which is important for larger heights. It is worth to point out that the above statement has a qualitative meaning when additional parameter is added. However, this kind of analysis will be presented in the dimensional analysis chap- ter4. The second approximation for small C, is R rnC_ 9 gh h2C gh~ lim(1 -Ch e 77 - g 2C R-o - ... (4.52) P cx->0 To 2To02 R 3A colleague asked this author to insert this explanation for his students. If you feel that it is too simple, please, just ignore it. 4These concepts are very essential in all the thermo-fluid science. I am grateful to my adviser E.R.G. Eckert who was the pioneer of the dimensional analysis in heat transfer and was kind to show me some of his ideas.  4.3. PRESSURE AND DENSITY IN A GRAVITATIONAL FIELD 67 Equation (4.52) shows that the correction factor (lapse coefficient), Cx, influences at only large values of height. It has to be noted that these equations (4.51) and (4.52) are not properly represented without the characteristic height. It has to be inserted to make the physical significance clearer. Equation (4.50) represents only the pressure ratio. For engineering pur- poses, it is sometimes important to obtain the density ratio. This relationship can be obtained from combining equations (4.50) and (4.45). The simplest assumption to combine these equations is by assuming the ideal gas model, equation (2.25), to yield P T p PTOR_ Ch _RxC h p PoT To 1 (4.53) *Advance material can be skiped* 4.3.4.2 The Stability Analysis It is interesting to study whether this solution (4.50) is stable and if h + dh so under what conditions. Suppose that for some reason, a small slab h of material moves from a layer at height, h, to layer at height h + dh (see Figure 4.7) What could hap- Fig. -4.7: Two adjoin layers for stability analysis. pen? There are two main possibili- ties one: the slab could return to the original layer or two: stay at the new layer (or even move further, higher heights). The first case is referred to as the stable condition and the second case referred to as the unstable condition. The whole system falls apart and does not stay if the analysis predicts unstable conditions. A weak wind or other disturbances can make the unstable system to move to a new condition. This question is determined by the net forces acting on the slab. Whether these forces are toward the original layer or not. The two forces that act on the slab are the gravity force and the surroundings pressure (buoyant forces). Clearly, the slab is in equilibrium with its surroundings before the movement (not necessarily stable). Under equilibrium, the body forces that acting on the slab are equal to zero. That is, the surroundings "pressure" forces (buoyancy forces) are equal to gravity forces. The buoyancy forces are proportional to the ratio of the density of the slab to surrounding layer density. Thus, the stability question is whether the slab density from layer h, p'(h) undergoing a free expansion is higher or lower than the density of the layer h + dh. If p'(h) > p(h + dh) then the situation is stable. The term p'(h) is slab from layer h that had undergone the free expansion. The reason that the free expansion is chosen to explain the process that the slab undergoes when it moves from layer h to layer h + dh is because it is the  68 CHAPTER 4. FLUIDS STATICS simplest. In reality, the free expansion is not far way from the actual process. The two processes that occurred here are thermal and the change of pressure (at the speed of sound). The thermal process is in the range of [cm/sec] while the speed of sound is about 300 [m/sec]. That is, the pressure process is about thousands times faster then the thermal process. The second issue that occurs during the "expansion" is the shock (in the reverse case [h + dh] -> h). However, this shock is insignificant (check book on Fundamentals of Compressible Flow Mechanics by this author on the French problem). The slab density at layer h + dh can be obtained using equation (4.53) as following p(h) PoT Cdh c)(1+CTdh -h P - 1+ C(4.54) The pressure and temperature change when the slab moves from layer at h to layer h + dh. The process, under the above discussion and simplifications, can be assumed to be adiabatic (that is, no significant heat transfer occurs in the short period of time). The little slab undergoes isentropic expansion as following for which (see equation (2.25)) p'(h+dh) (P'(h+dh)/k(4.55) p(h) P(h) When the symbol ' denotes the slab that moves from layer h to layer h + dh. The pressure ratio is given by equation (4.50) but can be approximated by equation (4.51) and thus p'(h + dh) gdh 1/k p(h) T(h) R)(4.56 Again using the ideal gas model for equation (4.57) transformed into p'(h+dh) pgdh 1k (4.57) p(h) P ) Expanding equation (4.57) in taylor series results in -pdh) / pdh (g2p2 k - g2 p2) dh2 1-( Pk 2P2k2 ... (4.58) The density at layer h+dh can be obtained from (4.54) and then it is expanded in taylor series as p(h + dh) ( C xdh a(R) ( Cx dh g p Cx p(h) )TO/ 1+ 7, ~1- dh+--- (4.59)  4.3. PRESSURE AND DENSITY IN A GRAVITATIONAL FIELD 69 The comparison of the right hand terms of equations (4.59) and (4.58) pro- vides the conditions to determine the stability. From a mathematical point of view, to keep the inequality for a small dh only the first term need to be compared as gp> gp _ C Pk P T (4.60) After rearrangement of the inequality (4.60) and using the ideal gas identity, it transformed to C> (k-1)gp T kP C < k - g (4.61) The analysis shows that the maximum amount depends on the gravity and gas properties. It should be noted that this value should be changed a bit since the k should be replaced by polytropic expansion n. When lapse rate C, is equal to the right hand side of the inequality, it is said that situation is neutral. However, one has to bear in mind that this analysis only provides a range and isn't exact. Thus, around this value additional analysis is needed 5. One of the common question this author has been asked is about the forces of continuation. What is the source of the force(s) that make this situation when unstable continue to be unstable? Supposed that the situation became unstable and the layers have been exchanged, would the situation become stable now? One has to remember that temperature gradient forces continuous heat transfer which the source temperature change after the movement to the new layer. Thus, the unstable situation is continuously unstable. 4.3.5 Gravity Variations Effects on Pressure and Density Until now the study focus on the change of density and pressure of the fluid. Equation (4.11) has two terms on the right hand side, the density, p and the body force, g. The body force was assumed until now to be constant. This assumption must be deviated when the distance from the body source is significantly change. At first glance, the body force is independent of the fluid. The source of the gravity force in gas is another body, while the gravity force source in liquid can be the liquid itself. Thus, the discussion Fig. -4.8: The varying gravity effects on density and pressure. 5The same issue of the floating ice. See example for the floating ice in cup.  70 CHAPTER 4. FLUIDS STATICS is separated into two different issues. The issues of magnetohydrodynamics are too advance for undergraduate student and thereforewill not be introduced here. 4.3.5.1 Ideal Gas in Varying Gravity In physics, it was explained that the gravity is a function of the distance from the center of the plant/body. Assuming that the pressure is affected by this gravity/body force. The gravity force is reversely proportional to r2. The gravity force can be assumed that for infinity, r - oc the pressure is about zero. Again, equation (4.11) can be used (semi one directional situation) when r is used as direction and thus DP G =_-p (4.62) where here G denotes the general gravity constant. The regular method of sepa- ration is employed to obtain fPdP - r r(4.63) JPbPRT ,b 2 where the subscript b denotes the conditions at the body surface. The integration of equation (4.63) results in in - = - - - - -(4.64) Pb RT rb r Or in a simplified form as _G r-r P RT '' (4.65) pb b Equation (4.65) demonstrates that the pressure is reduced with the distance. It can be noticed that for r - rb the pressure is approaching P - Pb. This equation confirms that the density in outer space is zero p(oc) = 0. As before, equation (4.65) can be expanded in taylor series as standard correction factor p - P G (r-rb)' (2GRT+G2rb) (rrb)2 Pb Pb R T 2rb (R T) 2 Notice that G isn't our beloved and familiar g and also that Grb/RT is a dimen- sionless number (later in dimensionless chapter about it and its meaning). 4.3.5.2 Real Gas in Varying Gravity The regular assumption of constant compressibility, Z, is employed. It has to re- member when this assumption isn't accurate enough, numerical integration is a  4.3. PRESSURE AND DENSITY IN A GRAVITATIONAL FIELD 71 possible solution. Thus, equation (4.63) is transformed into dP _ G dr f -dP C - (4.67) IPb PZRT Lb r With the same process as before for ideal gas case, one can obtain _G r-r Pb- ZRT rrb(4.68) pb b Equation (4.65) demonstrates that the pressure is reduced with the distance. It can be observed that for r - rb the pressure is approaching P - Pb. This equation confirms that the density in outer space is zero p(oo) = 0. As before taylor series for equation (4.65) is standard correctionfactor p P G G(r-rb)' (2GZRT+G2rb)(r-rb)2 - P-Z=T12-+ ... (4.69) pb Pb ZRT 2rb(ZRT)2 It can be noted that compressibility factor can act as increase or decrease of the ideal gas model depending on whether it is above one or below one. 4.3.5.3 Liquid in Under Varying Gravity For comparison reason consider the deepest location in the ocean which is about 11,000 [m]. If the liquid "equation of state" (4.35) is used with the hydrostatic fluid equation results in -rPoe (4.70) which the solution of equation (4.70) is e BT = Constant - BT g p0 (4.71) Since this author is not aware to which practical situation this solution should be applied, it is left for the reader to apply according to problem, if applicable. 4.3.6 Liquid Phase While for most practical purposes, the Cartesian coordinates provides sufficient treatment to the problem, there are situations where the spherical coordinates must be considered and used.  72 CHAPTER 4. FLUIDS STATICS Derivations of the fluid static in spherical coordinates are 1 d r2 dP r2 dr p dr + 47r Gp = 0 (4.72) Or in a vector form as Ve ( VP + 47Gp = 0 (4.73) 4.4 Fluid in a Accelerated System Up to this stage, body forces were considered as one-dimensional. The general, the linear acceleration have three components as opposed to the previous case of only one. However, the previous derivations can be easily extended. Equation (4.8) can be transformed into a different coordinate system where the main coordinate is in the direction of the effective gravity. Thus, the previous method can be used and there is no need to solve new three (or two) different equations. As before, the constant pressure plane is perpendicular to the direction of the effective gravity. Generally the acceleration is divided into two categories: linear and angular and they will be discussed in this order. 4.4.1 Fluid in a Linearly Accelerated System For example, in a two dimensional system, for the effective gravity geff = a i + g k (4.74) where the magnitude of the effective gravity is geff= =g2 + a2 (4.75) and the angle/direction can be obtained from tan/3 = (4.76) g Perhaps the best way to explain the linear acceleration is by examples. Con- sider the following example to illustrate the situation. Example 4.2: What is the angle of the liquid surface for a container in an accelerated system of a = 5[m/sec]?  4.4. FLUID IN A ACCELERATED SYSTEM 73 SOLUTION This question is one of the traditional question the fluid static and is straight for- Fig. -4.9: The effective gravity is for accelerated cart. ward. This requires to find the angle of the effective body force. This angle/direction can be found using the following tan-1 tan-a g 5 27.01° 9.81 The magnitude of the effective acceleration is ge ff l=/52 + 9.812 =11.015[m/sec2] Example 4.3: Cart partially is filled with liquid and is slid- ing on an inclined plane as shown in Figure 4.10. Calculate the shape of the surface. If there is a resistance what will be the an- gle? What happen when the slope angle is straight (the cart is dropping straight down)? Fig. -4.10: A cart slide on inclined plane SOLUTION (a) The angle can be found when the acceleration of the cart is found. If there is no resistance, the acceleration in the cart direction is determined from a = g sin #3 (4.77)  74 CHAPTER 4. FLUIDS STATICS The effective body force is acting perpendicular to the slope. Thus, the liquid surface is parallel to the surface of the inclination surface. (b) In case of resistance force (either of fric- tion due to the air or resistance in the wheels) reduces the acceleration of the cart. In that case the effective body moves closer to the gravity forces. The net body force depends on the mass anFnet of the liquid and the net acceleration is - m geff a = g - (4.78) g m The angle of the surface, a 4 gh Patmos p1 + Patmos 2 The full solution is (wR 4ghPatmo p+Patos 2 -2 g h w R pi -Pat mos w R) T1+2 g h Patmos w pi tmax (2 g w+4 g h) R pi ps T1 In this analysis the air temperature in the container immediately after inser- tion in the liquid has different value from the final temperature. It is reasonable as the first approximation to assume that the process is adiabatic and isentropic. Thus, the temperature in the cavity immediately after the insertion is  96 CHAPTER 4. FLUIDS STATICS Ti Pi Tf Pf The final temperature and pressure were calculated previously. The equation of state is mair RTi Pi = Vi The new unknown must provide additional equation which is Vi-= w2 hi Thickness Below The Maximum For the half thickness t =t2 the general solution for any given thickness below maximum is presented. The thickness is known, but the liquid displacement is still unknown. The pressure at the interface (after long time) is 7Tw2 hPHWsRT1 plghin +Patmos = (hin+h1 w which can be simplified to h Patmos p1 g hin + Patmos = a+ h The second equation is Archimedes' equation, which is 7rw2(h - h - h1) (w2 + 2wh)tpg) +w2 h Patmos g Example 4.11: A body is pushed into the liquid to a distance, ho and left at rest. Calculate accel- eration and time for a body to reach the surface. The body's density is a p1 , where a is ratio between the body density to the liquid density and (0 < a < 1). Is the body volume important? SOLUTION The net force is liquid body weight weight F= Vgp -Vgap =Vgpl(1-a) But on the other side the internal force is m F ma Vapi a  4.6. BUOYANCY AND STABILITY 97 Thus, the acceleration is a = g If the object is left at rest (no movement) thus time will be (h = 1/2 a t2) 2 ha g(1-a) If the object is very light (a - 0) then 2ha 2gh a2 3 2gha2 5 2gha2 ** g 2g 8g 16g From the above equation, it can be observed that only the density ratio is important. This idea can lead to experiment in "large gravity" because the acceleration can be magnified and it is much more than the reverse of free falling. Example 4.12: In some situations, itis desired to find equivalent of force of a certain shape to be replaced by another force of a "standard" shape. Consider the force that acts on a half sphere. Find equivalent cylinder that has the same diameter that SOLUTION The force act on the half sphere can be found by integrating the forces around the sphere. The element force is dAx h dA dF = (pL - ps) g r cos # cos 8cosO cos 0# r2 d8 d# The total force is then Fj(pj - ps) g cos2O#cos2r 83dod The result of the integration the force on sphere is 7r2 (pL - PS) r3 F f The force on equivalent cylinder is Fe =r2 (pL - ps) h  98 CHAPTER 4. FLUIDS STATICS These forces have to be equivalent and thus Thus the height is h 7r r 4 4.6.1 Stability Figure 4.30 shows a body made of hollow balloon and a heavy sphere Empty buoyancy connected by a thin and light rod. center This arrangement has mass cen- troid close to the middle of the sphere. The buoyant center is be- low the middle of the balloon. If this gravity arrangement is inserted into liquid center and it will be floating, the balloon Full will be on the top and sphere on a b c the bottom. Tilting the body with a small angle from its resting po- sition creates a shift in the forces Fig. -4.30: Schematic of floating bodies. direction (examine Figure 4.30b). These forces create a moment which wants to return the body to the resting (orig- inal) position. When the body is at the position shown in Figure 4.30c ,the body is unstable and any tilt from the original position creates moment that will further continue to move the body from its original position. This analysis doesn't violate the second law of thermodynamics. Moving bodies from an unstable position is in essence like a potential. A wooden cubic (made of pine, for ex- ample) is inserted into water. Part of the block floats above water line. The cubic mass (gravity) centroid is in the middle of G the cubic. However the buoyant center is the middle of the volume under the water (see Figure 4.31). This situation is similar B to Figure 4.30c. However, any experiment of this cubic wood shows that it is stable lo- cally. Small amount of tilting of the cubic Fig. -4.31: Schematic of floating cubic. results in returning to the original position. When tilting a larger amount than 7r/4 , it results in a flipping into the next stable  4.6. BUOYANCY AND STABILITY 99 position. The cubic is stable in six positions (every cubic has six faces). In fact, in any of these six positions, the body is in situation like in 4.30c. The reason for this local stability of the cubic is that other positions are less stable. If one draws the stability (later about this criterion) as a function of the rotation angle will show a sinusoidal function with four picks in a whole rotation. So, the body stability must be based on the difference between the body's local positions rather than the "absolute" stability. That is, the body is / "stable" in some points more than oth- M oAF ers in their vicinity. These points are raised from the buoyant force analysis. \F F GM dA When the body is tilted at a small angle, B #, the immersed part of the body center changes to a new location, B' as shown in Figure 4.32. The center of the mass (gravity) is still in the old location since the body did not change. The stability Fig. -4.32: Stability analysis of floating body. of the body is divided into three cate- gories. If the new immerse volume creates a new center in such way that the couple forces (gravity and buoyancy) try to return the body, the original state is re- ferred as the stable body and vice versa. The third state is when the couple forces do have zero moment, it is referred to as the neutral stable. The body, shown in Figure 4.32, when given a tilted position, move to a new buoyant center, B'. This deviation of the buoyant center from the old buoyant center location, B, should to be calculated. This analysis is based on the difference of the displaced liquid. The right green area (volume) in Figure 4.32 is displaced by the same area (really the volume) on left since the weight of the body didn't change10 so the total immersed section is constant. For small angle, #Q, the moment is calculated as the integration of the small force shown in the Figure 4.32 as AF. The displacement of the buoyant center can be calculated by examining the moment these forces are creating. The body weight creates opposite moment to balance the moment of the displaced liquid volume. BB' W = M (4.132) Where M is the moment created by the displaced areas (volumes), BB' is the distance between points B and point B', and, W referred to the weight of the body. It can be noticed that the distance BB' is an approximation for small angles (ne- glecting the vertical component.). So the perpendicular distance, BB', should be BB' = W (4.133) 10lt is correct to state: area only when the body is symmetrical. However, when the body is not symmetrical, the analysis is still correct because the volume and not the area is used.  100 CHAPTER 4. FLUIDS STATICS The moment M can be calculated as 6F M = fgpi x/3dA x = g pQ/x2dA (4.134) A A dV The integral in the right side of equation (4.134) is referred to as the area moment of inertia and was discussed in Chapter 3. The distance, BB' can be written from equation (4.134) as BB'= =g PI (4.135) Ps body The point where the gravity force direction is intersecting with the center line of the cross section is referred as metacentric point, M. The location of the meta- centric point can be obtained from the geometry as BM =(4.136) sin/# And combining equations (4.135) with (4.136) yields BM --=P1 1(4.137) f Ps Sill body Ps body For small angle (#Q~ 0) lim s ~N1(4.138) It is remarkable that the results is inde- pendent of the angle. Looking at Figure 4.32 ,the geometrical quantities can be related as h BMh1-.- GM- -BG (4.139) Ps body' a To understand these principles con- sider the following examples. Fig. -4.33: Cubic body dimensions for stabil- ity analysis. Example 4.13: A solid block of wood of uniform density, ps = a pi where (0 h1= P"h Pi Thus, the distance BG is (see Figure 4.31) hi BGU= hps7 h 1ps (4.140) 2 pi 2 2 p Stability of Square Block 3.0 2.5 - VPl 12.0h GM= 12- - 1- |psai 2 pl 50 . a 0.3 Simplifying the above equation provides 0.50 GM 1(a 2 1 0.0 - h - (1 - a) 0 (4.14 1) 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 where a is the density ratio. Notice April16, 2008 that GM/h isn't a function of the depth, L. This equation leads to the condition where the maximum height above which Fig. -4.34: Stability of cubic body infinity long. the body is not stable any more as h 6 (1 - a)a (4.142) One of the interesting point that the above analysis was that there is a point above which the ratio of the height to the body width is not stable any more. In cylindrical shape equivalent equation to equation (4.142) can be expressed. For cylinder (circle) the moment of inertia is Ixx = 7br4/64. The distance BG is the same as for the square shape (cubic) (see above (4.140)). Thus, the equation is GM g b2 h 64a hJ- -(1 -a)  102 CHAPTER 4. FLUIDS STATICS And the condition for maximum height for stability is b b> s32 (1- a) a h-> This kind of analysis can be carried for different shapes and the results are shown for these two shapes in Figure 4.35. Stability of Solid Blocks It can be noticed that the square 3.0c-square body is more stable than the circular 2.5 circle shape body. 2.0 Principle Main Axises .5 Any body has infinite number of 10 different axises around which moment of inertia can be calculated. For each of12 41 these axises there is a different moment April 16, 2008 0 of inertia. With the exception of the cir- cular shape, every geometrical shape Fig. -4.35: The maximum height reverse as a has an axis in which the moment of in- function of density ratio. ertia is without the product of inertia. This axis is where the main rotation of the body will occur. Some analysis of floating bodies are done by breaking the rotation of arbitrary axis to rotate around the two main axises. For stability analysis, it is enough to find if the body is stable around the smallest moment of inertia. For example, a square shape body has larger moment of inertia around diagonal. The difference between the previous calculation and the moment of inertia around the diagonal is I diagonal axis "normal"axis 2a(\a)4 AI~~G 5 a 3.~a Alxx= 6 12 0.07 a4 Which show that if the body is stable at main axises, it must be stable at the "diagonal" axis. Thus, this problem is reduced to find the stability for principle axis. Unstable Bodies What happen when one increases the height ratio above the maximum height ratio? The body will flip into the side and turn to the next stable point (angle). This is not a hypothetical question, but rather practical. This happens when a ship is overloaded with containers above the maximum height. In commercial ships, the fuel is stored at the bottom of the ship and thus the mass center (point G) is changing during the voyage. So, the ship that was stable (positive GM) leaving the initial port might became unstable (negative GM) before reaching the destination port. Example 4.14: In order to make the ship construction to be a hydrodynamic body is to make the  4.6. BUOYANCY AND STABILITY 103 body as narrow as possible. What is the minimum width of the two opposite sides triangular attached to each other see Figure ??. What is the minimum ratio width of the triangle to have a stable body up to a/h = 1. SOLUTION under construction 4.6.1.1 Stability of Body with Shifting Mass Centroid Ships and other floating bodies carry liquid or have a load which changes the mass location during tilting of the float- 'M ing body. For example, a ship that car- ries wheat grains where the cargo is Gc not properly secured to the ship. The G G' movement of the load (grains, furni- ture, and/or liquid) does not occur in the B ' same speed as the body itself or the dis- placed outside liquid. Sometimes, the slow reaction of the load, for stability Fig. -4.36: The effects of liquid movement on analysis, is enough to be ignored. Ex- the GM. act analysis requires taking into account these shifting mass speeds. However, here, the extreme case where the load re- acts in the same speed as the tilting of the ship/floating body is examined. For practical purposes, it is used as a limit for the stability analysis. There are situ- ations where the real case approaches to this extreme. These situations involve liquid with a low viscosity (like water, alcohol) and ship with low natural frequency (later on the frequency of the ships). Moreover, in this analysis, the dynamics are ignored and only the statics is examined (see Figure 4.36). A body is loaded with liquid "B" and is floating in a liquid "A" as shown in Figure 4.36. When the body is given a tilting position the body displaces the liquid on the outside. At the same time, the liquid inside is changing its mass centroid. The moment created by the inside displaced liquid is Min = g9p01I--B (4.1 43) Note that IxB isn't the same as the moment of inertia of the outside liquid inter- face.  104 CHAPTER 4. FLUIDS STATICS The change in the mass centroid of the liquid "A" then is G1G1= /f0x x B __Ix xB 1G1 VixB V- (4.144) VB Inside liquid weight Equation (4.144) shows that GG' is only a function of the geometry. This quantity, G1G1, is similar for all liquid tanks on the floating body. The total change of the vessel is then calculated similarly to center area calculations. gmtotai C b= gd-#mfG1G1 (4.145) For more than one tank, it can be written as GG' W GiGiV gVbi2p7 (4.146) Wtotai Wtotali= Vbi A new point can be defined as Gc. This point is the intersection of the center line with the vertical line form G'. GG' G Gc_ =(4.147) sin 1# The distance that was used before GM is replaced by the criterion for stability by Gc M and is expressed as Ge M = PA A - BG - 1 (4.148) psVbody mtotal Vb If there are more than one tank partially filled with liquid, the general formula is G M=pA IxxA BG- 1aib (4.149) Ps Vooiy m ita 1 Vbi One way to reduce the effect of the moving mass center due to liquid is done by substituting a single tank with several tanks. The moment of inertial of the combine two tanks is smaller than the moment of inertial of a single tank. Increasing the number of tanks reduces the moment of inertia. The engineer could design the tanks in such a way that the moment of inertia is operationally changed. This control of the stability, GM, can be achieved by having some tanks spanning across the entire body with tanks spanning on parts of the body. Movement of the liquid (mostly the fuel and water) provides way to control the stability, GM, of the ship. 4.6.1.2 Metacentric Height, GM, Measurement  4.6. BUOYANCY AND STABILITY 105 The metacentric height can be measured by finding the change in the angle when a weight is moved on the floating body.h Moving the weight, T a dis- tance, d then the moment created is Mweight = T d (4.150) This moment is balanced by Fig. -4.37: Measurement of GM of floating body. Mrighting = WtotaiGMnew 0 (4.151) Where, Wtotai, is the total weight of the floating body including measuring weight. The angle, 0, is measured as the difference in the orientation of the floating body. The metacentric height is GMnew Td (4.152) Wtotai 0 If the change in theGM can be neglected, equation (4.152) provides the solution. The calculation of GM can be improved by taking into account the effect of the measuring weight. The change in height of G is g mtotai Gnew = mship Gactuai + # T h (4.153) Combining equation (4.153) with equation (4.152) results in G3mactual = Gonewmtai - hT (4.154) mship mship The weight of the ship is obtained from looking at the ship depth. 4.6.1.3 Stability of Submerged Bodies The analysis of submerged bodied is different from the stability when the body lay between two fluid layers with different density. When the body is submerged in a single fluid layer, then none of the changes of buoyant centroid occurs. Thus, the mass centroid must be below the buoyant centroid in order to have stable condition. However, all fluids have density varied in some degree. In cases where the density changes significantly, it must be taken into account. For an example of such a case is an object floating in a solar pond where the upper layer is made of water with lower salinity than the bottom layer(change up to 20% of the density). When the floating object is immersed in the two layers, the stability analysis must  106 CHAPTER 4. FLUIDS STATICS take into account the changes of the displaced liquids of the two liquid layers. The calculations for such cases are a bit more complicated but based on the similar principles. Generally, this density change helps to increase the stability of the floating bodies. This analysis is out of the scape of this book (for now). 4.6.1.4 Stability of None Systematical or "Strange" Bodies While most floating bodies are symmet- rical or semi-symmetrical, there are sit- uations where the body has a "strange" AF and/or un-symmetrical body. Consider M M13 the first strange body that has an abrupt step change as shown in Figure 4.38. AF The body weight doesn't change during GM the rotation that the green area on the B left and the green area on right are the same (see Figure 4.38). There are two b situations that can occur. After the tilt- ing, the upper part of the body is above the liquid or part of the body is sub- merged under the water. The mathe- Fig. -4.38: Calculations of GM for abrupt matical condition for the border is when shape body. b = 3 a. For the case of b < 3 a the calculation of moment of inertia are similar to the previous case. The moment created by change in the displaced liquid (area) act in the same fashion as the before. The center of the moment is needed be found. This point is the intersection of the liquid line with the brown middle line. The moment of inertia should be calculated around this axis. For the case where b < 3 a x some part is under the liquid. The amount of area under the liquid section depends on the tilting angle. These calculations are done as if none of the body under the liquid. This point is intersection point liquid with lower body and it is needed to be calculated. The moment of inertia is calculated around this point (note the body is "ended" at end of the upper body). However, the moment to return the body is larger than actually was calculated and the bodies tend to be more stable (also for other reasons). 4.6.1.5 Neutral frequency of Floating Bodies This case is similar to pendulum (or mass attached to spring). The governing equation for the pendulum is - g = 0 (4.155) Where here f is length of the rode (or the line/wire) connecting the mass with the rotation point. Thus, the frequency of pendulum is 2 \ which measured in Hz.  4.6. BUOYANCY AND STABILITY 107 The period of the cycle is 27 /g. Similar situation exists in the case of floating bodies. The basic differential equation is used to balance and is rotation rotating moment I/ - V ps GM/3 =0 (4.156) In the same fashion the frequency of the floating body is 1 Vps G 217 VpGM (4.157) and the period time is 27 'rbody(4.158) VpsGM In general, the larger GM the more stable the floating body is. Increase in GM increases the frequency of the floating body. If the floating body is used for transport humans and/or other creatures or sensitive cargo it requires to reduce the GM so that the traveling will be smother. 4.6.2 Surface Tension The surface tension is one of the mathematically complex topic and related to many phenomena like boiling, coating, etc. In this section, only simplified topics like constant value will be discussed. Example 4.15: In interaction of the molecules shown in Figure ? describe the existence of surface tension. Explain why this description is erroneous? SOLUTION The upper layer of the molecules have unbalanced force towards the liquid phase. Newton's law states when there is unbalanced force, the body should be acceler- ate. However, in this case, the liquid is not in motion. Thus, the common explana- tion is wrong. Fig. -4.39: A heavy needle is floating on a liquid.  108 CHAPTER 4. FLUIDS STATICS Example 4.16: Needle is made of steel and is heavier than water and many other liquids. However, the surface tension between the needle and the liquid hold the needle above the liquid. After certain diameter, the needle cannot be held by the liquid. Calculate the maximum diameter needle that can be inserted into liquid without drowning. SOLUTION Under Construction 4.7 Rayleigh-Taylor Instability RayleighTaylor instability (or RT instability) is named after Lord Rayleigh and G. I. Taylor. There are situations where a heavy liquid layer is placed over a lighter fluid layer. This situation has engineering implications in several industries. For example in die casting, liquid metal is injected in a cavity filled with air. In poor designs or other situations, some air is not evacuated and stay in small cavity on the edges of the shape to be casted. Thus, it can create a situation where the liquid metal is above the air but cannot penetrate into the cavity because of instability. This instability deals with a dense, heavy fluid that is being placed above a lighter fluid in a gravity field perpendicular to interface. Example for such systems are dense water over oil (liquid-liquid), or water over air(gas-liquid). The original Rayleigh's paper deals with the dynamics and density variations. For example, density variations according to the bulk modulus (see section 4.3.3.2) are always stable but unstable of the density is in the reversed order. Supposed that a liquid density is arbitrary function of the height. This distor- tion can be as a result of heavy fluid above the lighter liquid. This analysis asks the question what happen when a small amount of liquid from above layer enter into the lower layer? Whether this liquid continue and will grow or will it return to its original conditions? The surface tension is the opposite mechanism that will returns the liquid to its original place. This analysis is referred to the case of infinite or very large surface. The simplified case is the two different uniform densities. For example a heavy fluid density, PL, above lower fluid with lower density, pG. For perfectly straight interface, the heavy fluid will stay above the lighter fluid. If the surface will disturbed, some of heavy liquid moves down. This disturbance can grow or returned to its original situation. This condition is determined by com- peting forces, the surface density, and the buoyancy forces. The fluid above the depression is in equilibrium with the sounding pressure since the material is ex- tending to infinity. Thus, the force that acting to get the above fluid down is the buoyancy force of the fluid in the depression. The depression is returned to its original position if the surface forces L h are large enough. In that case, this sit- uation is considered to be stable. On 17 I the other hand, if the surface forces Fig. -4.40: Description of depression to explain the Rayleigh-Taylor instability.  4.7. RAYLEIGH-TAYLOR INSTABILITY 109 (surface tension) are not sufficient, the situation is unstable and the heavy liq- uid enters into the liquid fluid zone and vice versa. As usual there is the neutral stable when the forces are equal. Any con- tinues function can be expanded in serious of cosines. Thus, example of a cosine function will be examined. The conditions that required from this function will be required from all the other functions. The disturbance is of the following h = -hmax COS 2x (4.159) L ( where hmax is the maximum depression and L is the characteristic length of the depression. The depression has different radius as a function of distance from the center of the depression, x. The weakest point is at x = 0 because symmetrical reasons the surface tension does not act against the gravity as shown in Figure (4.40). Thus, if the center point of the depression can "hold" the intrusive fluid then the whole system is stable. The radius of any equation is expressed by equation (1.47). The first deriva- tive of cos around zero is sin which is approaching zero or equal to zero. Thus, equation (1.47) can be approximated as 1 d~h 1 - d(4.160) R dx2 For equation (??) the radius is 1 _ 47 2hmax (4.161) R L2 - According to equation (1.38) the pressure difference or the pressure jump is due to the surface tension at this point must be PH-PL 4hmax 07r2 (4.162) The pressure difference due to the gravity at the edge of the disturbance is then PH - PL g (PH - PL) hmax (4.163) Comparing equations (4.162) and (4.163) show that if the relationship is 4o-72 L2 > g(pH-pL) (4.164) It should be noted that hmax is irrelevant for this analysis as it is canceled. The point where the situation is neutral stable 47r20-4.6 Lc = 4p2 (4.165) 9(EPH -- pL)  110 CHAPTER 4. FLUIDS STATICS An alternative approach to analyze this instability is suggested here. Con- sider the situation described in Figure 4.41. If all the heavy liquid "attempts" to move straight down, the lighter liquid will "prevent" it. The lighter liquid needs to move up at the same time but in a different place. The heavier liquid needs to move in one side and the lighter liquid in another location. In this process the heavier liq- uid "enter" the lighter liquid in one point and creates a depression as shown in Figure 4.41. To analyze it, considered two con- trol volumes bounded by the blue lines in Figure 4.41. The first control vol- 2r ume is made of a cylinder with a ra- dius r and the second is the depres- y sion below it. The "extra" lines of the depression should be ignored, they are not part of the control volume. The hori- zontal forces around the control volume are canceling each other. At the top, the force is atmospheric pressure times the area. At the cylinder bottom, the force Fig. -4.41: Description of depression to ex- is pgh x A. This acts against the grav- plain the instability. ity force which make the cylinder to be in equilibrium with its surroundings if the pressure at bottom is indeed pgh. For the depression, the force at the top is the same force at the bottom of the cylinder. At the bottom, the force is the integral around the depression. It can be approximated as a flat cylinder that has depth of r 7/4 (read the explanation in the example 4.12) This value is exact if the shape is a perfect half sphere. In reality, the error is not significant. Additionally when the depression occurs, the liquid level is reduced a bit and the lighter liquid is filling the missing portion. Thus, the force at the bottom is Fbottom ~w7 r2 + h) (pL - Pc) g + Patmos] (4.166) The net force is then Fbottom ~ 7r 2 (>)4(PL - Pc) g (4.167) The force that hold this column is the surface tension. As shown in Figure 4.41, the total force is then FF = 27,r6 cosO (4.168) The forces balance on the depression is then 2 x r cosO6 2(w4)(pL - pc) 9( (4-169)  4.7. RAYLEIGH-TAYLOR INSTABILITY 111 The radius is obtained by r~ 27rcosO (4.170) (PL - Pc) 9 The maximum surface tension is when the angle, 8 = 7r/2. At that case, the radius is is ~(4.171) PL- Pc) g 40 Fig. -4.42: The cross section of the interface. The purple color represents the maximum heavy liquid raising area. The yellow color represents the maximum lighter liquid that are "going down." The maximum possible radius of the depression depends on the geometry of the container. For the cylindrical geometry, the maximum depression radius is about half for the container radius (see Figure 4.42). This radius is limited because the lighter liquid has to enter at the same time into the heavier liquid zone. Since the "exchange" volumes of these two process are the same, the specific radius is limited. Thus, it can be written that the minimum radius is 27x6 rmin(tube 4.1 72) The actual radius will be much larger. The heavier liquid can stay on top of the lighter liquid without being turned upside down when the radius is smaller than the equation 4.172. This analysis introduces a new dimensional number that will be discussed in a greater length in the Dimensionless chapter. In equation (4.172) the angle was assumed to be 90 degrees. However, this angle is never obtained. The actual value of this angle is about r/4 to 7r/3 and in only extreme cases the angle exceed this value (considering dynamics). In Figure 4.42, it was shown that the depression and the raised area are the same. The actual area of the depression  112 CHAPTER 4. FLUIDS STATICS is only a fraction of the interfacial cross section and is a function. For examplethe depression is larger for square area. These two scenarios should be inserting into equation 4.168 by introducing experimental coefficient. Example 4.17: Estimate the minimum radius to insert liquid aluminum into represent tube at tem- perature of 600[K]. Assume that the surface tension is 400[mN/m]. The density of the aluminum is 2400kg/n3. SOLUTION The depression radius is assume to be significantly smaller and thus equation (4.171) can be used. The density of air is negligible as can be seen from the temperature compare to the aluminum density. a 87r 0.4 2400 x 9.81 The minimum radius is rr~ 0.02[m] which demonstrates the assumption of h >> r was appropriate.  CHAPTER 5 Multi-Phase Flow 5.1 Introduction Traditionally, the topic of multi-phase flow is ignored in an introductory class on fluid mechanics. For many engineers, this class will be the only opportunity to be exposed to this topic. The knowledge in this topic without any doubts, is required for many engineering problems. Calculations of many kinds of flow deals with more than one phase or material flow1. The author believes that the trends and effects of multiphase flow could and should be introduced and considered by engineers. In the past, books on multiphase flow were written more as a literature review or heavy on the mathematics. It is recognized that multiphase flow is still evolving. In fact, there is not a consensus to the exact map of many flow regimes. This book attempts to describe these issues as a fundamentals of physical aspects and less as a literature review. This chapter provides information that is more or less in consensus2. Additionally, the nature of multiphase flow requires solving many equations. Thus, in many books the representations is by writing the whole set governing equations. Here, it is believed that the interactions/calculations requires a full year class and hence, only the trends and simple calculations are described. 5.2 History The study of multi-phase flow started for practical purposes after World War II. Initially the models were using simple assumptions. For simple modelsthere are 1An example, there was a Ph.D. working for the government who analyzed filing cavity with liquid metal (aluminum), who did not consider the flow as two-phase flow and ignoring the air. As result, his analysis is in the twilight zone not in the real world. 20r when the scientific principles simply dictate. 113  114 CHAPTER 5. MULTI-PHASE FLOW two possibilities (1) the fluids/materials are flowing in well homogeneous mixed (where the main problem to find the viscosity), (2) the fluids/materials are flowing separately where the actual total loss pressure can be correlated based on the separate pressure loss of each of the material. If the pressure loss was linear then the total loss will be the summation of the two pressure losses (of the lighter liquid (gas) and the heavy liquid). Under this assumption the total is not linear and experimental correlation was made. The flow patterns or regimes were not considered. This was suggested by Lockhart and Martinelli who use a model where the flow of the two fluids are independent of each other. They postulate that there is a relationship between the pressure loss of a single phase and combine phases pressure loss as a function of the pressure loss of the other phase. It turned out this idea provides a good crude results in some cases. Researchers that followed Lockhart and Martinelli looked for a different map for different combination of phases. When it became apparent that specific models were needed for different situations, researchers started to look for different flow regimes and provided different models. Also the researchers looked at the situa- tion when the different regimes are applicable. Which leads to the concept of flow regime maps. Taitle and Duckler suggested a map based on five non-dimensional groups which are considered as the most useful today. However, Taitle and Duck- er's map is not universal and it is only applied to certain liquid-gas conditions. For example, Taitle-Duckler's map is not applicable for microgravity. 5.3 What to Expect From This Chapter As oppose to the tradition of the other chapters in this book and all other Potto project books, a description of what to expect in this chapter is provided. It is an attempt to explain and convince all the readers that the multi-phase flow must be included in introductory class on fluid mechanics3. Hence, this chapter will explain the core concepts of the multiphase flow and their relationship, and importance to real world. This chapter will provide: a category of combination of phases, the concept of flow regimes, multi-phase flow parameters definitions, flow parameters effects on the flow regimes, partial discussion on speed of sound of different regimes, double choking phenomenon (hopefully), and calculation of pressure drop of simple ho- mogeneous model. This chapter will introduce these concepts so that the engineer not only be able to understand a conversation on multi-phase but also, and more importantly, will know and understand the trends. However, this chapter will not provide a discussion of transient problems, phase change or transfer processes during flow, and actual calculation of pressure of the different regimes. 3This author feels that he is in an unique position to influence many in the field of fluid mechanics. This fact is due to the shear number of the downloaded Potto books. The number of the downloads of the book on Fundamental of compressible flow has exceed more than 100,000 in about two and half years. It also provides an opportunity to bring the latest advances in the fields since this author does not need to "sell" the book to a publisher or convince a "committee."  5.4. KIND OF MULTI-PHASE FLOW 5.4 Kind of Multi-Phase Flow Gas 115 Liquid Solid Liquid - Liquid - Liqu'd Gas Solid Fig. -5. 1: lines a static fluid with a constant density. All the flows are a form of multiphase flow. The discussion in the previous chapters is only as approximation when multiphase can be "reduced" into a single phase flow. For example, consider air flow that was discussed and presented earlier as a single phase flow. Air is not a pure material but a mixture of many gases. In fact, many proprieties of air are calculated as if the air is made of well mixed gases of Nitrogen and Oxygen. The results of the calculations of a mixture do not change much if it is assumed that the air flow as stratified flow 4 of many concentration layers (thus, many layers (infinite) of different materials). Practically for many cases, the homogeneous assumption is enough and suitable. However, this assumption will not be appropriate when the air is stratified because of large body forces, or a large acceleration. Adopting this assumption might lead to a larger error. Hence, there are situations when air flow has to be considered as multiphase flow and this effect has to be taken into account. In our calculation, it is assumed that air is made of only gases. The creation of clean room is a proof that air contains small particles. In almost all situations, the cleanness of the air or the fact that air is a mixture is ignored. The engineering accuracy is enough to totally ignore it. Yet, there are situations where cleanness of 4Different concentration of oxygen as a function of the height. While the difference of the concentra- tion between the top to button is insignificant, nonetheless it exists.  116 CHAPTER 5. MULTI-PHASE FLOW the air can affect the flow. For example, the cleanness of air can reduce the speed of sound. In the past, the breaks in long trains were activated by reduction of the compressed line (a patent no. 360070 issued to George Westinghouse, Jr., March 29, 1887). In a four (4) miles long train, the breaks would started to work after about 20 seconds in the last wagon. Thus, a 10% change of the speed of sound due to dust particles in air could reduce the stopping time by 2 seconds (50 meter difference in stopping) and can cause an accident. One way to categorize the multiphase is by the materials flows, For example, the flow of oil and water in one pipe is a multiphase flow. This flow is used by engineers to reduce the cost of moving crude oil through a long pipes system. The "average" viscosity is meaningless since in many cases the water follows around the oil. The water flow is the source of the friction. However, it is more common to categorize the flow by the distinct phases that flow in the tube. Since there are three phases, they can be solid-liquid, solid-gas, liquid-gas and solid-liquid- gas flow. This notion eliminates many other flow categories that can and should be included in multiphase flow. This category should include any distinction of phase/material. There are many more categories, for example, sand and grain (which are "solids") flow with rocks and is referred to solid-solid flow. The category of liquid-gas should be really viewed as the extreme case of liquid-liquid where the density ratio is extremely large. The same can be said for gas-gas flow. For the gas, the density is a strong function of the temperature and pressure. Open Channel flow is, although important, is only an extreme case of liquid-gas flow and is a sub category of the multiphase flow. The multiphase is an important part of many processes. The multiphase can be found in nature, living bodies (bio-fluids), and industries. Gas-solid can be found in sand storms, and avalanches. The body inhales solid particle with breathing air. Many industries are involved with this flow category such as dust collection, fluidized bed, solid propellant rocket, paint spray, spray casting, plasma and river flow with live creatures (small organisms to large fish) flow of ice berg, mud flow etc. The liquid-solid, in nature can be blood flow, and river flow. This flow also appears in any industrial process that are involved in solidification (for example die casting) and in moving solid particles. Liquid-liquid flow is probably the most common flow in the nature. Flow of air is actually the flow of several light liquids (gases). Many natural phenomenon are multiphase flow, for an example, rain. Many industrial process also include liquid-liquid such as painting, hydraulic with two or more kind of liquids. 5.5 Classification of Liquid-Liquid Flow Regimes The general discussion on liquid-liquid will be provided and the gas-liquid flow will be discussed as a special case. Generally, there are two possibilities for two differ- ent materials to flow (it is also correct for solid-liquid and any other combination). The materials can flow in the same direction and it is referred as co-current flow. When the materials flow in the opposite direction, it is referred as counter-current.  5.5. CLASSIFICATION OF LIQUID-LIQUID FLOW REGIMES 117 In general, the co-current is the more common. Additionally, the counter-current flow must have special configurations of long length of flow. Generally, the counter- current flow has a limited length window of possibility in a vertical flow in conduits with the exception of magnetohydrodynamics. The flow regimes are referred to the arrangement of the fluids. The main difference between the liquid-liquid flow to gas-liquid flow is that gas density is extremely lighter than the liquid density. For example, water and air flow as oppose to water and oil flow. The other characteristic that is different between the gas flow and the liquid flow is the variation of the density. For example, a reduction of the pressure by half will double the gas volumetric flow rate while the change in the liquid is negligible. Thus, the flow of gas-liquid can have several flow regimes in one situation while the flow of liquid-liquid will (probably) have only one flow regime. 5.5.1 Co-Current Flow In Co-Current flow, two liquids can have three main categories: vertical, horizontal, and what ever between them. The vertical configuration has two cases, up or down. It is common to differentiate between the vertical (and near vertical) and horizontal (and near horizontal). There is no exact meaning to the word "near vertical" or "near horizontal" and there is no consensus on the limiting angles (not to mention to have limits as a function with any parameter that determine the limiting angle). The flow in inclined angle (that not covered by the word "near") exhibits flow regimes not much different from the other two. Yet, the limits between the flow regimes are considerably different. This issue of incline flow will not be covered in this chapter. 5.5.1.1 Horizontal Flow The typical regimes for horizontal flow are stratified flow (open channel flow, and non open channel flow), dis-L- persed bubble flow, plug flow, and an- HeavyLiquid nular flow. For low velocity (low flow rate) of the two liquids, the heavy liq- uid flows on the bottom and lighter liq- Fig. -5.2: Stratified flow in horizontal tubes uid flows on the top5 as depicted in when the liquids flow is very slow. Figure 5.2. This kind of flow regime is referred to as horizontal flow. When the flow rate of the lighter liquid is almost zero, the flow is referred to as open channel flow. This definition (open channel flow) continues for small amount of lighter liq- uid as long as the heavier flow can be calculated as open channel flow (ignoring the lighter liquid). The geometries (even the boundaries) of open channel flow are very diverse. Open channel flow appears in many nature (river) as well in industrial 5With the exception of the extremely smaller diameter where Rayleigh-Taylor instability is an impor- tant issue.  118 CHAPTER 5. MULTI-PHASE FLOW process such as the die casting process where liquid metal is injected into a cylin- der (tube) shape. The channel flow will be discussed in a greater detail in Open Channel Flow chapter. As the lighter liquid (or the gas phase) flow rate increases (superficial veloc- ity), the friction between the phases increase. The superficial velocity is referred to as the velocity that any phase will have if the other phase was not exist. This friction is one of the cause for the instability which manifested itself as waves and changing the surface from straight line to a different configuration (see Figure 5.3). The wave shape is created to keep the gas and the liquid velocity equal and at the same time to have shear stress to be balance by surface tension. The configura- tion of the cross section not only depend on the surface tension, and other physical properties of the fluids but also on the material of the conduit. As the lighter liquid velocity increases two things can happen (1) wave size increase and (2) the shape of cross section continue to de- LightLiquid LightLiquid form. Some referred to this regime as wavy HeavyLiqui eavyLiqui stratified flow but this definition is not accepted by all as a category by itself. In fact, all the two phase flow are categorized by wavy flow which will proven later. There are two paths that can Fig. -5.3: Kind of Stratified flow in occur on the heavier liquid flow rate. If the heav- horizontal tubes. ier flow rate is small, then the wave cannot reach to the crown and the shape is deformed to the point that all the heavier liquid is around the periphery. This kind of flow regime is referred to as annular flow. If the heavier liquid flow rate is larger6 than the distance, for the wave to reach the conduit crown is smaller. At some point, when the lighter liquid flow increases, the heavier liquid wave reaches to the crown of the pipe. At this stage, the flow pattern is referred to as slug flow or plug flow. Plug flow is characterized by regions of lighter liquid filled with drops of the heavier liquid with Plug (or Slug) of the heavier liquid (with bubble of the lighter liquid). These plugs are separated by large "chunks" that almost fill the entire tube. The plugs are flowing in a succession (see Figure 5.4). The pressure drop of this kind of regime is significantly larger than the stratified flow. The slug flow cannot be assumed to be as homogeneous flow nor it can exhibit some average viscosity. The "average" viscosity depends on the flow and thus making it as insignificant way to do the calculations. Further increase of the lighter liquid flow rate move the flow regime into annular flow. Thus, the possibility to go through slug flow regime depends on if there is enough liquid flow rate. Choking occurs in compressible flow when the flow rate is above a cer- tain point. All liquids are compress- yiu - ible to some degree. For liquid which the density is a strong and primary function of the pressure, choking oc- Fig. -5.4: Plug flow in horizontal tubes when the liquids flow is faster slow. 6The liquid level is higher.  5.5. CLASSIFICATION OF LIQUID-LIQUID FLOW REGIMES 119 curs relatively closer/sooner. Thus, the flow that starts as a stratified flow will turned into a slug flow or stratified wavy7 flow after a certain distance depends on the heavy flow rate (if this category is accepted). After a certain distance, the flow become annular or the flow will choke. The choking can occur before the annular flow regime is obtained depending on the velocity and compressibility of the lighter liquid. Hence, as in compressible flow, liquid-liquid flow has a maximum combined of the flow rate (both phases). This maximum is known as double choking phenomenon. The reverse way is referred to the process where the starting point is high flow rate and the flow rate is decreasing. As in many fluid mechanics and magnetic fields, the return path is not move the exact same way. There is even a possibility to return on different flow regime. For example, flow that had slug flow in its path can be returned as stratified wavy flow. This phenomenon is refer to as hysteresis. Flow that is under small angle from the horizontal will be similar to the hor- izontal flow. However, there is no consensus how far is the "near" means. Quali- tatively, the "near" angle depends on the length of the pipe. The angle decreases with the length of the pipe. Besides the length, other parameters can affect the "near." Dispersed Bubble Elongated Bubble Slug Flow -' Annular Q Flow Stratified -. Flow J Wavy Stratified Open Channel Gas Superficial Velocity Flow Fig. -5.5: Modified Mandhane map for flow regime in horizontal tubes. The results of the above discussion are depicted in Figure 5.5. As many things in multiphase, this map is only characteristics of the "normal" conditions, e.g. in normal gravitation, weak to strong surface tension effects (air/water in "normal" 7Well, all the flow is wavy, thus it is arbitrary definition.  120 CHAPTER 5. MULTI-PHASE FLOW gravity), etc. 5.5.1.2 Vertical Flow Bubble Slug or Churn Annular Dispersed Flow Plug Flow Flow Flow Flow Fig. -5.6: lines a static fluid with a constant density. The vertical flow has two possibilities, with the gravity or against it. In en- gineering application, the vertical flow against the gravity is more common used. There is a difference between flowing with the gravity and flowing against the grav- ity. The buoyancy is acting in two different directions for these two flow regimes. For the flow against gravity, the lighter liquid has a buoyancy that acts as an "extra force" to move it faster and this effect is opposite for the heavier liquid. The oppo- site is for the flow with gravity. Thus, there are different flow regimes for these two situations. The main reason that causes the difference is that the heavier liquid is more dominated by gravity (body forces) while the lighter liquid is dominated by the pressure driving forces. Flow Against Gravity For vertical flow against gravity, the flow cannot start as a stratified flow. The heavier liquid has to occupy almost the entire cross section before it can flow be- cause of the gravity forces. Thus, the flow starts as a bubble flow. The increase of the lighter liquid flow rate will increase the number of bubbles until some bubbles start to collide. When many bubbles collide, they create a large bubble and the flow is referred to as slug flow or plug flow (see Figure 5.6). Notice, the different mechanism in creating the plug flow in horizontal flow compared to the vertical flow. Further increase of lighter liquid flow rate will increase the slug size as more bubbles collide to create "super slug"; the flow regime is referred as elongated bubble flow. The flow is less stable as more turbulent flow and several "super  5.5. CLASSIFICATION OF LIQUID-LIQUID FLOW REGIMES 121 slug" or churn flow appears in more chaotic way, see Figure 5.6. After additional increase of "super slug" , all these "elongated slug" unite to become an annular flow. Again, it can be noted the difference in the mechanism that create annular flow for vertical and horizontal flow. Any further increase transforms the outer liquid layer into bubbles in the inner liquid. Flow of near vertical against the gravity in two- phase does not deviate from vertical. The choking can occur at any point depends on the fluids and temperature and pressure. 5.5.1.3 Vertical Flow Under Micro Gravity The above discussion mostly explained the flow in a vertical configuration when Dipeed Dispesed the surface tension can be neglected. In cases where the surface tension is very Pulsing important. For example, out in space be- tween gas and liquid (large density dif- Pulsing& ference) the situation is different. The . Bubbling flow starts as dispersed bubble (some Trickling call it as "gas continuous") because the Flow Spray or Mist Flow gas phase occupies most of column. The liquid flows through a trickle or chan- Gas Flow Rage neled flow that only partially wets part of the tube. The interaction between the Fig. -5.7: A dimensional vertical flow map phases is minimal and can be considered under very low gravity against the as the "open channel flow" of the vertical gravity. configuration. As the gas flow increases, the liquid becomes more turbulent and some parts enter into the gas phase as drops. When the flow rate of the gas increases further, all the gas phase change into tiny drops of liquid and this kind of regime referred to as mist flow. At a higher rate of liquid flow and a low flow rate of gas, the regime liquid fills the entire void and the gas is in small bubble and this flow referred to as bubbly flow. In the medium range of the flow rate of gas and liquid, there is pulse flow in which liquid is mov- ing in frequent pulses. The common map is based on dimensionless parameters. Here, it is presented in a dimension form to explain the trends (see Figure 5.7). In the literature, Figure 5.7 presented in dimensionless coordinates. The abscissa is a function of combination of Froude ,Reynolds, and Weber numbers. The ordinate is a combination of flow rate ratio and density ratio. Flow With The Gravity As opposed to the flow against gravity, this flow can starts with stratified flow. A good example for this flow regime is a water fall. The initial part for this flow is more significant. Since the heavy liquid can be supplied from the "wrong" point/side, the initial part has a larger section compared to the flow against the grav- ity flow. After the flow has settled, the flow continues in a stratified configuration. The transitions between the flow regimes is similar to stratified flow. However, the  122 CHAPTER 5. MULTI-PHASE FLOW points where these transitions occur are different from the horizontal flow. While this author is not aware of an actual model, it must be possible to construct a model that connects this configuration with the stratified flow where the transitions will be dependent on the angle of inclinations. 5.6 Multi-Phase Flow Variables Definitions Since the gas-liquid system is a specific case of the liquid-liquid system, both will be united in this discussion. However, for the convenience of the terms "gas and liquid" will be used to signify the lighter and heavier liquid, respectively. The liquid- liquid (also gas-liquid) flow is an extremely complex three-dimensional transient problem since the flow conditions in a pipe may vary along its length, over its cross section, and with time. To simplify the descriptions of the problem and yet to retain the important features of the flow, some variables are defined so that the flow can be described as a one-dimensional flow. This method is the most common and im- portant to analyze two-phase flow pressure drop and other parameters. Perhaps, the only serious missing point in this discussion is the change of the flow along the distance of the tube. 5.6.1 Multi-Phase Averaged Variables Definitions The total mass flow rate through the tube is the sum of the mass flow rates of the two phases m=mrnG-| L (5.1) It is common to define the mass velocity instead of the regular velocity because the "regular" velocity changes along the length of the pipe. The gas mass velocity is GG = (5.2) Where A is the entire area of the tube. It has to be noted that this mass velocity does not exist in reality. The liquid mass velocity is GL=L (5.3) The mass flow of the tube is then G = (5.4) It has to be emphasized that this mass velocity is the actual velocity. The volumetric flow rate is not constant (since the density is not constant) along the flow rate and it is defined as QG = = UsG (5.5) PG  5.6. MULTI-PHASE FLOW VARIABLES DEFINITIONS 123 and for the liquid QL = G(5.6) PL For liquid with very high bulk modulus (almost constant density), the volumetric flow rate can be considered as constant. The total volumetric volume vary along the tube length and is Q=QL+QG (5.7) Ratio of the gas flow rate to the total flow rate is called the 'quality' or the "dryness fraction" and is given by X= =_C0(5.8) In a similar fashion, the value of (1 - X) is referred to as the "wetness fraction." The last two factions remain constant along the tube length as long the gas and liquid masses remain constant. The ratio of the gas flow cross sectional area to the total cross sectional area is referred as the void fraction and defined as AG a A (5.9) This fraction is vary along tube length since the gas density is not constant along the tube length. The liquid fraction or liquid holdup is L H1-aL A (5.10) It must be noted that Liquid holdup, LH is not constant for the same reasons the void fraction is not constant. The actual velocities depend on the other phase since the actual cross sec- tion the phase flows is dependent on the other phase. Thus, a superficial velocity is commonly defined in which if only one phase is using the entire tube. The gas superficial velocity is therefore defined as GG X m UsG A QG (5.11) P pGA The liquid superficial velocity is GL (1- X) m s PL PLA QL (5.12) Since UsL = QL and similarly for the gas then Um = UsG + UsL (5.13)  124 CHAPTER 5. MULTI-PHASE FLOW Where Um is the averaged velocity. It can be noticed that Um is not constant along the tube. The average superficial velocity of the gas and liquid are different. Thus, the ratio of these velocities is referred to as the slip velocity and is defined as the following SLP UL (5.14) Slip ratio is usually greater than unity. Also, it can be noted that the slip velocity is not constant along the tube. For the same velocity of phases (SLP = 1), the mixture density is defined as pm =apG + (1 -a)pL (5.15) This density represents the density taken at the "frozen" cross section (assume the volume is the cross section times infinitesimal thickness of dx). The average density of the material flowing in the tube can be evaluated by looking at the definition of density. The density of any material is defined as p = rn/V and thus, for the flowing material it is p = Q(5.16) Where Q is the volumetric flow rate. Substituting equations (5.1) and (5.7) into equation (5.16) results in mG iL X n+(1 -X) i Paverage G + QL Xm+ (1 -X~i h(1-X~ii PG PL QC QL (5.17) X)+X = Equation (5.17) can be simplified by canceling the ?i and noticing the (1- 1 to become (5.18) The average specific volume of the flow is then 1 X+ (1l-X)x ±( Vaverage P P+ P = X VG + paverage PG PL X) VL (5.19)  5.7. HOMOGENEOUS MODELS 125 The relationship between X and a is AG X - - pc UcAa pcUca a5.20 mG + rL PLULA(1-a)+pcUGAa pLUL(1-a)+pcUa AL If the slip is one SLP = 1, thus equation (5.20) becomes X =p pcLa(5.21) PLt(1 -a) +pc a 5.7 Homogeneous Models Before discussing the homogeneous models, it is worthwhile to appreciate the complexity of the flow. For the construction of fluid basic equations, it was as- sumed that the flow is continuous. Now, this assumption has to be broken, and the flow is continuous only in many chunks (small segments). Furthermore, these seg- ments are not defined but results of the conditions imposed on the flow. In fact, the different flow regimes are examples of typical configuration of segments of contin- uous flow. Initially, it was assumed that the different flow regimes can be neglected at least for the pressure loss (not correct for the heat transfer). The single phase was studied earlier in this book and there is a considerable amount of information about it. Thus, the simplest is to used it for approximation. The average velocity (see also equation (5.13)) is Um = QL +Qc = UsL+ UsG = Um (5.22) It can be noted that the continuity equation is satisfied as 5 =pm Um A (5.23) Example 5.1: Under what conditions equation (5.23) is correct? SOLUTION Under construction The governing momentum equation can be approximated as TdUm dP m d = -A - S &- Apm gsinO (5.24) dxdx  126 CHAPTER 5. MULTI-PHASE FLOW or modifying equation (5.24) as |dP _S midUm . - =- - + pmg sing dx -A A ddmI The energy equation can be approximated as dy dw d Ur2 =_ x fx2hx hm+ + g x sinO dx dx dx 2 ~ (5.25) (5.26) 5.7.1 Pressure Loss Components In a tube flowing upward in incline angle 0, the pressure loss is affected by fric- tion loss, acceleration, and body force(gravitation). These losses are non-linear and depend on each other. For example, the gravitation pressure loss reduce the pressure and thus the density must change and hence, acceleration must occur. However, for small distances (dx) and some situations, this dependency can be neglected. In that case, from equation (5.25), the total pressure loss can be written as (5.27) Every part of the total pressure loss will be discussed in the following section. 5.7.1.1 Friction Pressure Loss The frictional pressure loss for a conduit can be calculated as dP S dx f A Where S is the perimeter of the fluid. For calculating the frictional the pipe is (5.28) pressure loss in dP dx, The wall shear stress can be estimated by 4Tw D (5.29) Pm Urn2 T f2 (5.30)  5.7. HOMOGENEOUS MODELS 127 The friction factor is measured for a single phase flow where the average velocity is directly related to the wall shear stress. There is not available experimental data for the relationship of the averaged velocity of the two (or more) phases and wall shear stress. In fact, this friction factor was not measured for the "averaged" viscosity of the two phase flow. Yet, since there isn't anything better, the experimental data that was developed and measured for single flow is used. The friction factor is obtained by using the correlation f = C PmUmD (5.31) Where C and n are constants which depend on the flow regimes (turbulent or laminar flow). For laminar flow C = 16 and n = 1. For turbulent flow C = 0.079 and n = 0.25. There are several suggestions for the average viscosity. For example, Duckler suggest the following Q+G QG /LQL (5.32) QG -| L QG -| L Duckler linear formula does not provide always good approximation and Cichilli suggest similar to equation (5.33) average viscosity as /1average= x +1-X) (5.33) G + L Or simply make the average viscosity depends on the mass fraction as pm = X pG + (1 - X) pL (5.34) Using this formula, the friction loss can be estimated. 5.7.1.2 Acceleration Pressure Loss The acceleration pressure loss can be estimated by =dP _ ndm (5.35) dx a dx The acceleration pressure loss (can be positive or negative) results from change of density and the change of cross section. Equation (5.35) can be written as dP d __d(5.36) dx a dx ApmJ Or in an explicit way equation (5.36) becomes pressure loss due to pressure loss due to dP _2density change area change = [r2(5.37) dx a l d 1 1 dA A dx pm) + pm A2 dx  128 CHAPTER 5. MULTI-PHASE FLOW There are several special cases. The first case where the cross section is constant, dA/ dx= 0. In second case is where the mass flow rates of gas and liquid is constant in which the derivative of X is zero, dX/ dxo= 0. The third special case is for constant density of one phase only, dpL/ dxo= 0. For the last point, the private case is where densities are constant for both phases. 5.7.1.3 Gravity Pressure Loss Gravity was discussed in Chapter 4 and is dP =d gpmsinO (5.38) The density change during the flow can be represented as a function of density. The density in equation (5.38) is the density without the "movement" (the "static" density). 5.7.1.4 Total Pressure Loss The total pressure between two points, (a and b) can be calculated with integration as rb dP APab = dx (5.39) Iad and therefore f riction acceleration gravity APab-- APab-f + APaba +APabg (5.40) 5.7.2 Lockhart Martinelli Model The second method is by assumption that every phase flow separately One such popular model by Lockhart and Martinelli8. Lockhart and Martinelli built model based on the assumption that the separated pressure loss are independent from each other. Lockhart Martinelli parameters are defined as the ratio of the pressure loss of two phases and pressure of a single phase. Thus, there are two parameters as shown below. dP dP C dx TP dx(5.41) e ct p SGt 8This method was considered a military secret, private communication with Y., Taitle  5.8. SOLID-LIQUID FLOW 129 Where the TP denotes the two phases and SG denotes the pressure loss for the single gas phase. Equivalent definition for the liquid side is dP dP dP dP(5.42) L dx TP dx SL f Where the SL denotes the pressure loss for the single liquid phase. The ratio of the pressure loss for a single liquid phase and the pressure loss for a single gas phase is dP dP dx SL dxSG (543) where u is Martinelli parameter. It is assumed that the pressure loss for both phases are equal. dP _ dP dx dx(5.44) SG SL The pressure loss for the liquid phase is dP 2 f L UL 2 pl dP fUD p(5.45) dxL DL For the gas phase, the pressure loss is dP 2 fG UG2p( dxD(5.46) dG DG Simplified model is when there is no interaction between the two phases. To insert the Diagram. 5.8 Solid-Liquid Flow Solid-liquid system is simpler to analyze than the liquid-liquid system. In solid- liquid, the effect of the surface tension are very minimal and can be ignored. Thus, in this discussion, it is assumed that the surface tension is insignificant compared to the gravity forces. The word "solid" is not really mean solid but a combination of many solid particles. Different combination of solid particle creates different "liquid." Therefor,there will be a discussion about different particle size and different geometry (round, cubic, etc). The uniformity is categorizing the particle sizes, distribution, and geometry. For example, analysis of small coal particles in water is different from large coal particles in water. The density of the solid can be above or below the liquid. Consider the case where the solid is heavier than the liquid phase. It is also assumed that the "liquids" density does not change significantly and it is far from the choking point. In that case there are four possibilities for vertical flow:  130 CHAPTER 5. MULTI-PHASE FLOW 1. The flow with the gravity and lighter density solid particles. 2. The flow with the gravity and heavier density solid particles. 3. The flow against the gravity and lighter density solid particles. 4. The flow against the gravity and heavier density solid particles. All these possibilities are different. However, there are two sets of similar characteristics, possibility, 1 and 4 and the second set is 2 and 3. The first set is similar because the solid particles are moving faster than the liquid velocity and vice versa for the second set (slower than the liquid). The discussion here is about the last case (4) because very little is known about the other cases. 5.8.1 Solid Particles with Heavier Density ps > PL Solid-liquid flow has several combination flow regimes. When the liquid velocity is very small, the liquid cannot carry the solid parti- cles because there is not enough resistance to lift up the solid particles. A particle in a middle of the vertical liquid flow experience several forces. The force balance of spherical particle in field viscous fluid (creeping flow) is gravity and buoyancy forces drag forces rD3g(ps-pL) _ CDOOWD2pLUL2 (5.47) 6 8 Where CDo is the drag coefficient and is a function of Reynolds number, Re, and D is the equivalent radius of the particles. The Reynolds number defined as Re = UL D (5.48) Inserting equating (5.48) into equation (5.47) become CDo(UL ) f(Re) U2D 3PL)(5.49) 3 pL Equation (5.49) relates the liquid velocity that needed to maintain the particle "floating" to the liquid and particles properties. The drag coefficient, CDo is com- plicated function of the Reynolds number. However, it can be approximated for several regimes. The first regime is for Re < 1 where Stokes' Law can be approxi- mated as 24 ~Re (5.50)  5.8. SOLID-LIQUID FLOW 131 In transitional region 1 < Re < 1000 24 (_1 3 RDo-e 1 ) e/ (5.51) For larger Reynolds numbers, the Newton's Law region, CDo, is nearly constant as CD o = 0.44 (5.52) In most cases of solid-liquid system, the Reynolds number is in the second range9. For the first region, the velocity is small to lift the particle unless the density differ- ence is very small (that very small force can lift the particles). In very large range (especially for gas) the choking might be approached. Thus, in many cases the middle region is applicable. So far the discussion was about single particle. When there are more than one particle in the cross section, then the actual velocity that every particle expe- rience depends on the void fraction. The simplest assumption that the change of the cross section of the fluid create a parameter that multiply the single particle as CD oa CD.o f(a) (5.53) When the subscript a is indicating the void, the function f(a) is not a linear func- tion. In the literature there are many functions for various conditions. Minimum velocity is the velocity when the particle is "floating". If the velocity is larger, the particle will drift with the liquid. When the velocity is lower, the particle will sink into the liquid. When the velocity of liquid is higher than the minimum velocity many particles will be floating. It has to remember that not all the particle are uniform in size or shape. Consequently, the minimum velocity is a range of velocity rather than a sharp transition point. As the solid particles are not pushed by a pump but moved by the forces the fluid applies to them. Thus, the only velocity that can be applied is the fluid velocity. Yet, the solid parti- Trasiton cles can be supplied at different rate. Packed Thus, the discussion will be focus on partialy Fully the fluid velocity. For small gas/liquid solid Fuze velocity, the particles are what some particles fluidized call fixed fluidized bed. Increasing the flow fluid velocity beyond a minimum will move the particles and it is referred to as mix fluidized bed. Additional in- APtube crease of the fluid velocity will move 91t be wonderful if flow was in the last range? TI Ari- v t} u 199996Yr solid particles.  132 CHAPTER 5. MULTI-PHASE FLOW all the particles and this is referred to as fully fluidized bed. For the case of liquid, further increase will create a slug flow. This slug flow is when slug shape (domes) are almost empty of the solid particle. For the case of gas, additional increase create "tunnels" of empty almost from solid particles. Additional increase in the fluid velocity causes large turbulence and the ordinary domes are replaced by churn type flow or large bubbles that are almost empty of the solid particles. Further increase of the fluid flow increases the empty spots to the whole flow. In that case, the sparse solid particles are dispersed all over. This regimes is referred to as Pneumatic conveying (see Figure 5.9). Fixed Mixed Slug or Turbulent Fast Pneumatic Bed Bed Plug Flow Regimes Fluidization Conveying Fig. -5.9: The flow patterns in solid-liquid flow. One of the main difference between the liquid and gas flow in this category is the speed of sound. In the gas phase, the speed of sound is reduced dramatically with increase of the solid particles concentration (further reading Fundamentals of Compressible Flow" chapter on Fanno Flow by this author is recommended). Thus, the velocity of gas is limited when reaching the Mach somewhere between 1/ k and 1 since the gas will be choked (neglecting the double choking phenomenon). Hence, the length of conduit is very limited. The speed of sound of the liquid does not change much. Hence, this limitation does not (effectively) exist for most cases of solid-liquid flow. 5.8.2 Solid With Lighter Density ps < p and With Gravity This situation is minimal and very few cases exist. However, it must be pointed out that even in solid-gas, the fluid density can be higher than the solid (especially with  5.9. COUNTER-CURRENT FLOW 133 micro gravity). There was very little investigations and known about the solid-liquid flowing down (with the gravity). Furthermore, there is very little knowledge about the solid-liquid when the solid density is smaller than the liquid density. There is no known flow map for this kind of flow that this author is aware of. Nevertheless, several conclusions and/or expectations can be drawn. The issue of minimum terminal velocity is not exist and therefor there is no fixed or mixed fluidized bed. The flow is fully fluidized for any liquid flow rate. The flow can have slug flow but more likely will be in fast Fluidization regime. The forces that act on the spherical particle are the buoyancy force and drag force. The buoyancy is accelerating the particle and drag force are reducing the speed as ,rwD3g(ps - pL) _ CD 7rD2pL (Us - UL)2 (5.54) 6 8 From equation 5.54, it can observed that increase of the liquid velocity will increase the solid particle velocity at the same amount. Thus, for large velocity of the fluid it can be observed that UL/Us -- 1. However, for a small fluid velocity the velocity ratio is very large, UL/Us -- 0. The affective body force "seems" by the particles can be in some cases larger than the gravity. The flow regimes will be similar but the transition will be in different points. The solid-liquid horizontal flow has some similarity to horizontal gas-liquid flow. Initially the solid particles will be carried by the liquid to the top. When the liquid velocity increase and became turbulent, some of the particles enter into the liquid core. Further increase of the liquid velocity appear as somewhat similar to slug flow. However, this author have not seen any evidence that show the annular flow does not appear in solid-liquid flow. 5.9 Counter-Current Flow This discussion will be only on liquid-liquid systems (which also includes liquid- gas systems). This kind of flow is probably the most common to be realized by the masses. For example, opening a can of milk or juice. Typically if only one hole is opened on the top of the can, the liquid will flow in pulse regime. Most people know that two holes are needed to empty the can easily and continuously. Otherwise, the flow will be in a pulse regime. In most cases, the possibility to have counter-current flow is limited to having short length of tubes. In only certain configurations of the infinite long pipes the counter-current flow can exist. In that case, the pressure difference and gravity (body forces) dominates the flow. The inertia components of the flow, for long tubes, cannot compensate for the pressure gradient. In short tube, the pressure difference in one phase can be positive while the pressure difference in the other phase can be negative. The pressure difference in the interface must be finite. Hence, the counter-current flow can have opposite pressure gradient for short conduit. But in most cases, the heavy phase (liquid) is pushed by the gravity and lighter phase (gas) is driven by the pressure difference.  134 CHAPTER 5. MULTI-PHASE FLOW Fig. -5. 10: Counter-current flow in a can (the left figure) has only one hole thus pulse flow and a flow with two holes (right picture). The counter-current flow occurs, for example, when cavity is filled or emptied with a liquid. The two phase regimes "occurs" mainly in entrance to Annular the cavity. For example, Figure 5.10 Flow Open depicts emptying of can filled with liq- ° """'' uid. The air is "attempting" to enter Pulse the cavity to fill the vacuum created Inpossible Flow thus forcing pulse flow. If there are oppi two holes, in some cases, liquid flows Flow through one hole and the air through f(D/L,physicalproperties) the second hole and the flow will be continuous. It also can be noticed that if there is one hole (orifice) and a long Fig. -5.11: Counter-flow in vertical tubes map. and narrow tube, the liquid will stay in the cavity (neglecting other phenom- ena such as dripping flow.). There are three flow regimes10 that have been observed. The first flow pat- tern is pulse flow regime. In this flow regime, the phases flow turns into different direction (see Figure 5.12). The name pulse flow is used to signify that the flow is flowing in pulses that occurs in a certain frequency. This is opposed to counter- current solid-gas flow when almost no pulse was observed. Initially, due to the gravity, the heavy liquid is leaving the can. Then the pressure in the can is reduced compared to the outside and some lighter liquid (gas)entered into the can. Then, the pressure in the can increase, and some heavy liquid will starts to flow. This process continue until almost the liquid is evacuated (some liquid stay due the sur- face tension). In many situations, the volume flow rate of the two phase is almost equal. The duration the cycle depends on several factors. The cycle duration can be replaced by frequency. The analysis of the frequency is much more complex issue and will not be dealt here. 10Caution! this statement should be considered as "so far found". There must be other flow regimes that were not observed or defined. For example, elongated pulse flow was observed but measured. This field hasn't been well explored. There are more things to be examined and to be studied.  5.9. COUNTER-CURRENT FLOW 135 Fig. -5.12: Pictures of Counter-current flow in liquid-gas and solid-gas configurations. The container is made of two compartments. The upper compartment is filled with the heavy phase (liquid, water solution, or small wood particles) by rotating the container. Even though the solid-gas ratio is smaller, it can be noticed that the solid-gas is faster than the liquid-gas flow. Annular Flow in Counter-current flow The other flow regime is annular flow in which the heavier phase is on the periph- Water ery of the conduit (In the literature, there Flow are someone who claims that heavy liquid will be inside). The analysis is provided, but somehow it contradicts with the experimen- tal evidence. Probably, one or more of the assumptions that the analysis based is er- Steam roneous). In very small diameters of tubes Flow the counter-current flow is not possible be- cause of the surface tension (see section 4.7). The ratio of the diameter to the length Fig. -5.13: Flood in vertical pipe. with some combinations of the physical properties (surface tension etc) determines the point where the counter flow can start. At this point, the pulsing flow will start and larger diameter will increase the flow and turn the flow into annular flow. Ad- ditional increase of the diameter will change the flow regime into extended open channel flow. Extended open channel flow retains the characteristic of open chan- nel that the lighter liquid (almost) does not effect the heavier liquid flow. Example of such flow in the nature is water falls in which water flows down and air (wind) flows up. The driving force is the second parameter which effects the flow existence. When the driving (body) force is very small, no counter-current flow is possible. Consider the can in zero gravity field, no counter-current flow possible. However, if the can was on the sun (ignoring the heat transfer issue), the flow regime in the can moves from pulse to annular flow. Further increase of the body force will move  136 CHAPTER 5. MULTI-PHASE FLOW the flow to be in the extended "open channel flow." In the vertical co-current flow there are two possibilities, flow with gravity or against it. As opposed to the co-current flow, the counter-current flow has no possibility for these two cases. The heavy liquid will flow with the body forces (gravity). Thus it should be considered as non existent flow. 5.9.1 Horizontal Counter-Current Flow Up to this point, the discussion was focused on Sigph the vertical tubes. In horizontal tubes, there oDlpping is an additional flow regime which is stratified'l Horizontal flow is different from vertical flow from the stability issues. A heavier liquid layer stratied can flow above a lighter liquid. This situation is unstable for large diameter but as in static (see section (4.7) page 108) it can be considered sta- f(D/L,physicalproperties) ble for small diameters. A flow in a very nar- row tube with heavy fluid above the lighter fluid Fig. -5.14: A flow map to explain should be considered as a separate issue. the horizontal counter- When the flow rate of both fluids is very current flow. small, the flow will be stratified counter-current flow. The flow will change to pulse flow when the heavy liquid flow rate increases. Further increase of the flow will result in a single phase flow regime. Thus, closing the window of this kind of flow. Thus, this increase terminates the two phase flow possibility. The flow map of the horizontal flow is different from the vertical flow and is shown in Figure 5.14. A flow in an angle of inclination is closer to vertical flow unless the angle of inclination is very small. The stratified counter flow has a lower pressure loss (for the liquid side). The change to pulse flow increases the pressure loss dramatically. 5.9.2 Flooding and Reversal Flow The limits of one kind the counter-current flow regimes, that is stratified flow are discussed here. This problem appears in nuclear engineering (or boiler engineer- ing) where there is a need to make sure that liquid (water) inserted into the pipe reaching the heating zone. When there is no water (in liquid phase), the fire could melt or damage the boiler. In some situations, the fire can be too large or/and the water supply failed below a critical value the water turn into steam. The steam will flow in the opposite direction. To analyze this situation consider a two dimensional conduit with a liquid inserted in the left side as depicted in Figure 5.13. The liquid velocity at very low gas velocity is constant but not uniform. Further increase of the gas velocity will reduce the average liquid velocity. Additional increase of the gas velocity will bring it to a point where the liquid will flow in a reverse direction and/or disappear (dried out).  5.9. COUNTER-CURRENT FLOW 137 A simplified model for this situation is for a two dimensional configuration where the liquid is flowing down and the gas is flowing up as shown in Figure 5.15. It is assumed that both fluids are flowing in a laminar regime and steady state. Additionally, it is assumed that the entrance ef- fects can be neglected. The liquid flow rate, QL, is unknown. However, the pressure difference in the (x direction) is known and equal to zero. The boundary conditions for the liquid is that velocity at the wall is zero and the velocity at the inter- face is the same for both phases UG = UL or rlG= L. As it will be shown later, both con- ditions cannot coexist. The model can be im- proved by considering turbulence, mass trans- y D h r - L Liquid Flow Gas Flow Fig. -5.15: A diagram to explain the flood in a two dimension geometry. fer, wavy interface, etc'1. This model is presented to exhibits the trends and the special features of counter-current flow. Assuming the pressure difference in the flow direction for the gas is constant and uniform. It is assumed that the last as- sumption does not contribute or change significantly the results. The underline rational for this assumption is that gas density does not change significantly for short pipes (for more information look for the book "Fundamentals of Compress- ible Flow" in Potto book series in the Fanno flow chapter.). The liquid film thickness is unknown and can be expressed as a function of the above boundary conditions. Thus, the liquid flow rate is a function of the boundary conditions. On the liquid side, the gravitational force has to be balanced by the shear forces as dTx_ dx The integration of equation (5.55) results in x= PL g x + C1 (5.55) (5.56) The integration constant, C1, can be found from the boundary condition where rxy(x = h) =Ti. Hence, i = pL g h+ C1 The integration constant is then Ci = i - PL g h which leads to xy =pLg(x-h)+Ti (5.57) (5.58) Substituting the newtonian fluid relationship into equation (5.58) to obtained dU L pid * =p g (x -h) +Tri (5.59) 11The circular configuration is under construction and will be appeared as a separated article mo- mentarily.  138 CHAPTER 5. MULTI-PHASE FLOW or in a simplified form as dUypL g (x - h) Tr dxp + E ot toLyeL Equation (5.60) can be integrate to yield (5.60) The liquid velocity at the wall, [U(x = 0) can be found to be hx) + +C2 pL (5.61) = 0], is zero and the integration coefficient C2 = 0 (5.62) The liquid velocity profile is then (5.63) The velocity at the liquid-gas interface is Uy(x = h) = i h _ pL g h2 PL 21pL (5.64) The velocity can vanish (zero) inside the obtained from film in another point which can be O= PL9(x2 Si 2 The solution for equation (5.65) is hx +TZx L (5.65) xu o = 0 2h 2 r IlL gPL (5.66) The maximum x value is limited by the liquid film thickness, h. The minimum shear stress that start to create reversible velocity is obtained when x = h which is 0 _ L g(h2 hh + h / L (5.67) h g PL 2 If the shear stress is below this critical shear stress Tio then no part of the liquid will have a reversed velocity. The notation of Tio denotes the special value at which a starting shear stress value is obtained to have reversed flow. The point where the liquid flow rate is zero is important and it is referred to as initial flashing point.  5.9. COUNTER-CURRENT FLOW 139 The flow rate can be calculated by integrating the velocity across the entire liquid thickness of the film. Q w h Uydx= 0o .h uL 9 x2 Jo . L 2 hx + Ti' dx L (5.68) Where w is the thickness of the conduit (see Figure 5.15). Integration equation (5.68) results in Q _ h2 (37 - 2ghpL) w 6 Lp (5.69) It is interesting to find the point where the liquid mass flow rate is zero. This point can be obtained when equation (5.69) is equated to zero. There are three solutions for equation (5.69). The first two solutions are identical in which the film height is h = 0 and the liquid flow rate is zero. But, also, the flow rate is zero when 3 i = 2 g h pL. This request is identical to the demand in which - --(5.70) This critical shear stress, for a given film thickness, reduces the flow rate to zero or effectively "drying" the liquid (which is different then equation (5.67)). For this shear stress, the critical upward interface velocity is f fa/r (5.71) The wall shear stress is the last thing that will be done on the shear stress is liquid side. The wall __ (5.72) TL Iwau dU IL dx IL + 3 3 =0 IL x= 0 Simplifying equation (5.72)12 becomes (notice the change of for the direction) the sign accounting TL @ alo =ghpL TL ~wa11 3 (5.73) 12AIso noticing that equation (5.70) has to be equal g h pL to support the weight of the liquid.  140 CHAPTER 5. MULTI-PHASE FLOW Again, the gas is assumed to be in a laminar flow as well. The shear stress on gas side is balanced by the pressure gradient in the y direction. The momentum balance on element in the gas side is dTxG _GdP (5.74) dx dy The pressure gradient is a function of the gas compressibility. For simplicity, it is assumed that pressure gradient is linear. This assumption means or implies that the gas is incompressible flow. If the gas was compressible with an ideal gas equation of state then the pressure gradient is logarithmic. Here, for simplicity reasons, the linear equation is used. In reality the logarithmic equation should be used (a discussion can be found in "Fundamentals of Compressible Flow" a Potto project book). Thus, equation (5.74) can be rewritten as X - X(5.75) dx Ay L Where Ay = L is the entire length of the flow and AP is the pressure difference of the entire length. Utilizing the Newtonian relationship, the differential equation is d2UG /XP dPdx2 0L (5.76) Equation (5.76) can be integrated twice to yield /xP UG= P 2 +C1 x+C2 (5.77) PG L This velocity profile must satisfy zero velocity at the right wall. The velocity at the interface is the same as the liquid phase velocity or the shear stress are equal. Mathematically these boundary conditions are UG(x =D) =0 (5.78) and UG (x = h) = UL(x = h) (a) or (5.79) TG (x =h) = TL(x = h) (b) Applying B.C. (5.78) into equation (5.77) results in UG = 0 = D2 + C1 D + C2 (5.80) iMG L C2 = - D2 +C1 D G L  5.9. COUNTER-CURRENT FLOW 141 Which leads to U(2 - D2) +C1 (x -D) (5.81) G L At the other boundary condition, equation (5.79)(a), becomes pLgh2 IP (h2 - D2) +C1 (h - D) (5.82) 6/pL pG0L The last integration constant, C1 can be evaluated as pLg h2 _/P (h+D) (5.83) 6 pL (h - D) pG L With the integration constants evaluated, the gas velocity profile is AP 2- pLgh2 (x-D) _P (h+ D) (x-D) UG 0 2 ( - D2) + (5.84 PG L 6 Lp(h - D) pG L The velocity in Equation (5.84) is equal to the velocity equation (5.63) when (x = h). However, in that case, it is easy to show that the gas shear stress is not equal to the liquid shear stress at the interface (when the velocities are assumed to be the equal). The difference in shear stresses at the interface due to this assumption, of the equal velocities, cause this assumption to be not physical. The second choice is to use the equal shear stresses at the interface, condi- tion (5.79)(b). This condition requires that dUG dUL(5.85 G x L (5.85) The expressions for the derivatives are gas side liquid side 2h/P 2ghpL L +GC1- ghp=(5.86) As result, the integration constant is C1 =2IPL 2hAP (5.87) 1 3 pG MG) The gas velocity profile is then AP 2 2ghpL (hDP UGP2 - D2)+2 (x - D) (5.88) PGL (3pG1G)L  142 CHAPTER 5. MULTI-PHASE FLOW The gas velocity at the interface is then UG ,axn _ = P (h2 D2) + (2ghPL 2hAP (h G L ( D) (5.89) This gas interface velocity is different than the velocity of the liquid side. The velocity at interface can have a "slip" in very low density and for short distances. The shear stress at the interface must be equal, if no special effects occurs. Since there no possibility to have both the shear stress and velocity on both sides of the interface, different thing(s) must happen. It was assumed that the interface is straight but is impossible. Then if the interface becomes wavy, the two conditions can co-exist. The wall shear stress is TG Qwall G d or in a simplified form as Ji(P2x (2hPL 2 hAP G L x= (5.90) =D 2AP (D - h) 2ghpL TG Qwall L + 3 (5.91) The Required Pressure Difference The pressure difference to create the flooding (drying) has to take into account the fact that the surface is wavy. T However, as first estimate the waviness of the surface can be L7- J1pgLh neglected. The estimation of 1LTwG the pressure difference under the assumption of equal shear stress can be applied. In the DAP same fashion the pressure dif- ference under the assumption the equal velocity can be calcu- Fig. -5.16: General forces diagram to calculated the in lated. The actual pressure dif- a two dimension geometry. ference can be between these two assumptions but not must be between them. This model and its assumptions are too simplistic and the actual pressure difference is larger. However, this expla- nation is to show magnitudes and trends and hence it provided here. To calculate the required pressure that cause the liquid to dry, the total bal- ance is needed. The control volume include the gas and liquid volumes. Figure 5.16 describes the general forces that acts on the control volume. There are two  5.10. MULTI-PHASE CONCLUSION 143 forces that act against the gravity and two forces with the gravity. The gravity force on the gas can be neglected in most cases. The gravity force on the liquid is the liquid volume times the liquid volume as Volme/w FgL = pg h L (5.92) The total momentum balance is (see Figure 5.16) A/w A/w force due to pressure FgL + L TwG= L TwL+ D/XP (5.93) Substituting the different terms into (5.93) result in pgLh++L (2AP (D-h) +2gp) Lg3 +DXP (5.94) Simplifying equation (5.94) results in 4pgLh=(2hD)AP (5.95) 3 or AP = 4 pg Lh(.6 3 (2h-D) (5.96) This analysis shows far more reaching conclusion that initial anticipation ex- pected. The interface between the two liquid flowing together is wavy. Unless the derivations or assumptions are wrong, this analysis equation (5.96) indicates that when D > 2 h is a special case (extend open channel flow). 5.10 Multi-Phase Conclusion For the first time multi-phase is included in a standard introductory textbook on fluid mechanics. There are several points that should be noticed in this chapter. There are many flow regimes in multi-phase flow that "regular" fluid cannot be used to solve it such as flooding. In that case, the appropriate model for the flow regime should be employed. The homogeneous models or combined models like Lockhart-Martinelli can be employed in some cases. In other case where more accurate measurement are needed a specific model is required. Perhaps as a side conclusion but important, the assumption of straight line is not appropriate when two liquid with different viscosity are flowing.  144 CHAPTER 5. MULTI-PHASE FLOW  SUBJECTS INDEX Subjects Index 145 A absolute viscosity, 7, 8, 12 Arc shape, 88 Archimedes, 3 Atmospheric pressure, 59 B Bingham's model, 11 Body force, 55, 56, 58, 67, 69 effective, 57 bulk modulus, 19, 21 buoyancy, 3, 92, 93 buoyant force, 67, 99 C Co-current flow, 117 Compressibility factor, 61, 71 Correction factor, 66 Counter-current Pulse flow, 134 Counter-current flow, 117, 133 Annular flow, 135 Extended Open channel flow, 135 Cut-out shapes, 86 D D'Alembert paradox, 4 Density ratio, 67, 101 dilettante, 11 Double choking phenomenon, 119 F Fixed fluidized bed, 131 Flow regime map, 114 Flow regimes in one pipe, 118 Free expansion, 67-69 Fully fluidized bed, 132 G Gas-gas flow, 116 H horizontal counter-current flow, 136 Horizontal flow, 117 Hydrostatic pressure, 55, 88 Ideal gas, 61 K kinematic viscosity, 11 L Lapse rate, 69 Linear acceleration, 72 Liquid phase, 63 Liquid-Liquid Regimes, 116 Lockhart martinelli model, 128 M "Magnification factor", 61 Mass velocity, 122 Metacentric point, 100 Minimum velocity solid-liquid flow, 130 Mixed fluidized bed, 131 Multi-phase flow, 113 Multiphase flow against the gravity, 120 N Neutral moment Zero moment, 99 Neutral stable, 69, 99, 109 Newtonian fluids, 1, 8 0 Open channel flow, 117 p Pendulum action, 106 Piezometric pressure, 58 Pneumatic conveying, 132 Polynomial function, 89 Pressure center, 81 pseudoplastic, 11 Pulse flow, 134  146 CHAPTER 5. MULTI-PHASE FLOW purely viscous fluids, 11 Vertical counter-current flow, 134 Vertical flow, 117 R Rayleigh-Taylor Instability, 108 W Rayleigh-Taylor instability, 117 Watson's method, 19 Real gas, 61 Westinghouse patent, 116 Return path for flow regimes, 119 S Scalar function, 56, 85 shear stress, 6 Solid-fluid flow Gas dynamics aspects, 132 Solid-fluid flow, 129 Solid-liquid flow, 129 Solid-solid flow, 116 Spherical coordinates, 71 Stability analysis, 67 stability analysis, 69, 92 cubic, 99 Stability in counter-current flow, 136 Stable condition, 67, 105 stratified flow, 117 Superficial velocity, 118 Sutherland's equation, 12 T Terminal velocity, 130 thixotropic, 11 Total moment, 75 Triangle shape, 88, 95 Two-Phase Gas superficial velocity, 123 Liquid holdup, 123 Quality of dryness, 123 Reversal flow, 136 Slip velocity, 124 Void Fraction, 123 Wetness fraction, 123 U Unstable condition, 67 V Vapor pressure, 59  AUTHORS INDEX Authors Index B Bhuckingham, 4 Blasiu, 4 Blasius, 4 C Cichilli, 127 D Darcy, 4 Duckler, 114, 127 F Fanning, 4 Froude, 4 G Ganguillet, 4 H Helmholtz, Hermann von, 4 K Kutta-Joukowski, 4 L Lockhart, 114 M Manning, 4 Martinelli, 114 Meye, 4 N Nikuradse, 4 P Prandtl, 4 R Rayleigh, 4,108 Rose, 4 147 S Stanton, 4 T Taitle, 114 Taylor, G.I., 108 V von Karma, 4 W Weisbach, 4 Westinghouse, 116