Force analysis of underwater object with supercavitation evolution
|
|
- Ronald Haynes
- 5 years ago
- Views:
Transcription
1 Indian Journal of Geo-Marine Sciences Vol. 42(8), December 2013, pp Force analysis of underwater object with supercavitation evolution B C Khoo 1,2,3* & J G Zheng 1,3 1 Department of Mechanical Engineering, National University of Singapore, Singapore Singapore-MIT Alliance, National University of Singapore, Singapore Temasek Laboratories, National University of Singapore, Singapore * [ mpekbc@nus.edu.sg] Reeived 5 December 2012; revised 9 September 2013 Supercavitation generally occurs as a result of flow acceleration along underwater body surface and is numerically investigated in this study using a compressible Navier-Stokes equations solver. Here, the supercavitating flow is assumed to be the homogeneous mixture of pure liquid water and vapour which are in kinematic and thermodynamic equilibrium. Liquid phase and cavitation are modeled by Tait equation of state (EOS) and isentropic one-fluid formulation, respectively. Convective terms of the governing equations are numerically integrated using Godunov-type, cell-centered finite volume MUSCL scheme on unstructured triangular mesh, whereas time integration is handled with the second-order accurate Runge-Kutta approach. Our interest is focused on the force analysis of traveling object with the formation, growth, evolution and even collapse of supercavity enveloping the object. It is found that skin friction drag exerted on the object can be reduced significantly by the formation of supercavity where viscosity of vapour is much smaller than that of liquid water. It is also observed that form drag acting on the object is influenced by the supercavitation. Collapse of supercavity over the body due to external perturbation not only damages underwater object but also alters form drag on it. [Keywords: Supercavitating flow, Isentropic cavitation model, Underwater object, Navier-Stokes equations] Introduction The numerical simulation and prediction of cavitation/supercavitation are of great importance for efficient design of engineering devices such as hydrofoil ship, hydraulic turbine and propeller which may suffer from the damage due to cavitation bubble collapse. However, a lot of numerical and physical modeling difficulties are encountered in the computational fluid dynamics simulation and analysis. The discontinuities of fluid properties are present across cavitation bubble interface, which include large density ratio of liquid to vapour, drastic sound speed variation and completely different equations of state for liquid and vapour. On the other hand, the cavitation phenomenon is related to a number of complex physics such as significant mass transfer rate, viscous effect, interfacial dynamics and so on. All the issues mentioned above constitute challenges for any numerical methods which should be able to model the cavitation problem accurately and be robust as well as efficient. Various numerical methods have been developed to resolve the cavitation. Among them, multi-phase method is implemented by solving the separate system of governing equations for each phase. This type of approach is able to model fluid dynamics and thermodynamics of liquid and vapour, respectively. Although the multi-phase model is physically reasonable, it is usually quite complex, involving various transfer terms, and therefore not easy to implement. On the contrary, one-fluid model treats the cavitating flow as the homogeneous mixture of liquid and vapour and is used widely 1,2,3,4. In most cavitation simulations, the liquid flow is assumed to be incompressible as the flow speed is relatively low. This assumption is consistent with the experimental observation that the fluid density almost does not change when the flow remains in low or moderate Mach number regime. A consequence of the incompressible flow assumption is that the method cannot model the perturbation propagation through the flow at finite speed. However, under a number of situations such perturbations like shock wave and other pressure wave may have significant effect on the cavitation evolution and should be modeled. In this paper, a homogeneous cavitating flow model based on compressible Navier-Stokes equations is employed. The liquid is described as
2 958 INDIAN J. MAR. SCI., VOL. 42, NO. 8, DECEMBER 2013 being compressible and modeled by Tait equation of state (EOS). The cavitation is treated as the homogeneous mixture of isentropic liquid as well as vapour and its inception and evolution is modeled by isentropic cavitation model. The phase transition is assumed to occur instantaneously and other physics associated with the phase change are neglected. The present method allows us to model the shock wave propagation through liquid media and its interaction with cavitation bubble. This study is concentrated on the simulation of pressure wave induced supercavitation deformation and collapse and on their influence on the underwater object. Materials and Methods In the present study, the cavitating flow is assumed to be the homogeneous mixture of isentropic liquid and vapour which remain in kinematic and thermodynamic equilibrium, that is, liquid and vapour share the same velocity, pressure and temperature. Baseline partial differential equation system describing the cavitating flow is the compressible Navier-Stokes equations, U F G i Fv G v i S = + + S t x y y x y y v (1) where i takes value of 0 and 1 for planar and axisymmetric flows, respectively. Here, U, F, G and S are the vector of conserved variables, inviscid flux vectors in x and y direction, geometric source term vector associated with axisymmetric flow, respectively, and are given by, ρ ρu ρv ρv ρ ρ ρ ρ 2 U= u F= u + p G= uv S= uv 2 2 ρv ρuv ρv + p ρv (2) In (2), ρ is the averaged density, u is the averaged x velocity component, v is the averaged y velocity component and p is the averaged pressure. Energy equation is neglected in model system (1) as the flow is assumed to be isentropic and the pressure depends on the density only. The viscous flux vectors F v, G v and S v are not presented here to save space. The liquid phase is modeled by Tait EOS, 1/N p + B A ρ = ρ0, p psat. B (3) When the flow pressure drops below the saturated vapour pressure p sat, the cavitation occurs and is modeled by isentropic cavitation model 3, kρ + ρ =, p < p p + B A p + k p + B A p cav cav g l ρ 1/ N 1/ γ cav cav sat (4) where k=α 0 /(1-α 0 ) and α 0 is the known void fraction of the mixture density at p cav. Here, ρ and ρ are cav g cav l the associated vapour and liquid densities at the cavitation pressure p cav. For liquid water, values of parameters in (3) are set to be N=7.15, A=10 5 Pa, B = Pa, p o =1000 Kg/m 3, respectively. The sound speed is calculated through Wallis s formula 5, α ( 1 α ) ρ ρ ρ a = sv. av sw. aw 1 2 (5) where void fraction α is determined from the mixture density definition, sv ( 1 ) ρ = αρ + α ρ (6) sw In (6), p sw and p sv are the saturated liquid and vapour densities, respectively. Viscosity coefficient for the two-phase mixture in (1) is given by, v ( 1 ) µ = α µ + α µ (7) l where µ v and µ 1 are viscosity coefficients of vapour and liquid, respectively. The convective and viscous terms of model system (1) closed by (3) and (4) are discretized using secondorder MUSCL 6,7 scheme and an algorithm based on calculation of average of gradients 8, respectively, while the time-marching is dealt with using two-stage Runge-Kutta method. Results and Discussion Supercavitating flow around a hemispherical head cylinder To validate our method, a supercavitating flow around an underwater cylinder with hemispherical head is resolved. The computational domain is shown in Fig. 1 where the inflow boundary condition is imposed on the left hand side of the domain. The
3 KHOO & ZHENG: FORCE ANALYSIS OF UNDERWATER OBJECT 959 Fig. 1 The supercavitating flow around a hemispherical head cylinder. Upper left: the triangular mesh; upper right: close-up view of the mesh around the cylinder head; lower left: the intermediate stage of supercavitation evolution; lower right: the steady state supercavitation. In all the coordinate systems shown in the figures of this paper, unless stated otherwise, unit of length is meter non-slip boundary is implemented along the cylinder surface while the remaining sides are treated as open boundary. The freestream pressure and velocity are set to be P =10 5 Pa, U =30 m/s, respectively. The flow is assumed to be symmetric. The triangular mesh employed in the simulation and its close-up view around cylinder head are illustrated in the first row of Fig. 1. The unstructured mesh is better able to conform to the geometric features. The density contour maps showing the supercavitation evolution process is demonstrated in the second row of Fig. 1. In the lower left frame of Fig. 1, the cavitation bubble has appeared. The cavitation bubble continues growing and finally develops into a supercavity enveloping the whole cylinder as illustrated in the lower right frame of Fig. 1 where the flow has reached steady state. The density in the cavitation region is several orders of magnitude smaller than that of the liquid outside the supercavitation bubble. In all the simulations presented in the paper, Courant-Friedrichs-Lewy (CFL) number is set to be 0.8. The code is seen to accurately resolve the supercavitating flow and is stable and robust. Supercavitating flow past a high-speed underwater projectile The second validation case is on the computation of a high-speed conical-shaped projectile travelling at velocity of 970 m/s in water. A disk cavitator of radius 0.71 mm is mounted on the head of the projectile which has length of mm and base radius of 6.56 mm. The freestream pressure is Pa. The experimental 9 and numerical images for supercavity encompassing the entire object are compared visually in Fig. 2(a,b). The supercavitation
4 960 INDIAN J. MAR. SCI., VOL. 42, NO. 8, DECEMBER 2013 Fig. 2a,b Experimental (a) and numerical (b) images for a high-subsonic underwater projectile at speed of 970 m/s and comparison between experimental and calculated supercavity profiles (c). Fig. 2c profile and wake resolved with isentropic model are qualitatively comparable to their counterparts in experimental image, as shown in Fig. 2. In Fig. 2(c), the computed supercavity thickness concurs quantitatively with Hrubes supercavity size measurements 9, which demonstrates the validity of our method. Supercavitating flow interacting with a pressure wave In this case, the supercavitation bubble impacted by a pressure wave is numerically investigated. Initially, a blunt head cylinder with radius of 10 mm is immersed in a uniform water flow at pressure of 10 5 Pa and velocity of 100 m/s. The cylinder is at the angle of attack of 0. The boundary conditions are the same as those in the first case except that the freestream pressure and velocity at inlet are set to be P =10 5 Pa, U =100 m/s, respectively. Due to the symmetry of the flow, only the upper half of the flow field is calculated. The above initial conditions will give rise to a supercavitation bubble enveloping the entire cylinder. Fig. 3 depicts part of the density field where time t=0 corresponds to the time instant when the numerical solution is converged. The supercavity interface is smooth and there are no any pressure oscillations. The supercavitating flow is well resolved as illustrated in Fig. 3.
5 KHOO & ZHENG: FORCE ANALYSIS OF UNDERWATER OBJECT 961 In the steady state calculation, the pressure and velocity at inlet are fixed as constants. For unsteady simulation, after the supercavity is formed and the solution is converged, the freestream pressure is increased suddenly to a higher value and then returns to its initial value after a short time period of 66 µs. The pressure pulse will generate complex pressure wave which is to collide with the supercavity. The situation is schematized in Fig. 4. Four sets of density evolution maps corresponding to different pressure pulses are illustrated in Fig. 5. In Fig. 5 (a1-e1), the freestream pressure is abruptly increased by 100 times, from 10 5 Pa to 10 7 Pa. It is observed that the supercavity is deformed by the pressure wave. The perturbation travels through water quite fast and the locally deformed boundary gradually recovers to its original profile, see Fig. 5 (a1-e1). In Fig. 5 (a2-e2), the pressure is increased to P= Pa. The supercavitation is seen to undergo more severe deformations. The qualitative behaviors of supercavitation evolution associated with the two pressure pulses are similar. With increasing peak pressure to P= Pa, the supercavity locally collapses and is divided into two parts as shown in Fig. 5 (a3-e3). The right cavity fully collapses with time advancing while the left one expands downstream and finally develops into a new supercavity. The evolution process for P= Pa is similar to that with the peak pressure of Pa. The skin friction drag on the cylinder is reduced significantly by the supercavitation formation as the viscosity in vapour is much smaller than that in water. However, the form drag exerted on the body increases due to impact of generated pressure wave on cylinder head. Supercavitating flow around a 2D plate Next, we consider the supercavitating flow past two dimensional (2D) plate inclined relative to the flow Fig. 3 The density field map of a supercavitation around a blunt head cylinder. Here, P =10 5 Pa and U =100 m/s. Fig. 4 The schematic of simulation setup for the supercavitation interacting with a shock wave. Fig. 5 The density field evolution for the supercavitation bubble impacted by a pressure wave. Frames (a1-e1): P=10 7 Pa; frames (a2-e2): P= Pa; frames (a3-e3): P= Pa; frames (a4-e4): P= Pa.
6 962 INDIAN J. MAR. SCI., VOL. 42, NO. 8, DECEMBER 2013 Fig. 6 The supercavitating 2D plate at different angles of attack (the first row) and shock induced supercavity collapse (the second row) at angle of attack α=5. direction. The flow is not symmetric and therefore assumed to be planar. The plate is 150 mm long and 20 mm high. The freestream pressure and velocity are 10 5 Pa and 100 m/s, respectively. The density contour images associated with angles of attack of 5, 10 and 15 at steady state are presented in Fig. 6. The angles of attack are selected in such a way that supercavitating flow exhibits different patterns. It is obvious that cavity shape depends on the angle of attack. The upper cavity over the upper surface of the plate becomes thicker with increasing attack angle. The cavity below lower surface appears to be more affected by inclined angle. At α=10, the lower cavity is narrower than that at α=5. For the larger angle of 15, only a thin layer of cavitation is formed near plate head. It is known that the drag acting on an object may be influenced by cavity over it. The form drag coefficients for α=5 and α=10 are and 0.063, respectively. Here, the drag coefficient is 2 defined as C = F / (0.5P U A ) where F d is the d d ref form drag, P and U are the freestream density and velocity, respectively, and A ref is the reference area. For the two angles of attack of 5 and 10, the force on the left side of the plate is the major source of drag and hence the two coefficients are close to each other. However, at α=15, the pressure contribution to the drag on the lower surface outside the cavity has to be accounted for and therefore larger coefficient of results. The calculation indicates that in engineering practice, devices may need to operate within a certain range of angle of attack to avoid the possible situation like that at α=15, which may lead to larger drag and cause a degradation of devices efficiency. For the case with α=5, a Mach 1.1 incident shock is introduced after a steady state solution is obtained and the numerical results are presented in Fig. 6. As illustrated in Fig. 6(d), at t=2 τ, the cavities originally generated and anchored by plate corners are shrinking from their leading edges. Here, τ represents a time period of 0.1 ms. Simultaneously, a relatively small cavity appears near the upper left corner of plate and grows. The lower cavity appears to shrink more rapidly than the upper one as time advances, see Fig. 6(e) at t=4 τ. On the other hand, two water jets are formed along top and bottom of the plate and the lower jet penetrates deeper into cavity than the upper one. By t=12 τ, the cavities have completely collapsed as shown Fig. 6(f). After the cavities vanish, the form drag coefficient for the flow configuration with 5 angle of attack is around 0.13 which is about 2 times the value with the cavities.
7 KHOO & ZHENG: FORCE ANALYSIS OF UNDERWATER OBJECT 963 Conclusion The supercavitation interacting with the pressure wave is numerically simulated using a homogeneous cavitation fluid method. It is found that the supercavity may be locally distorted or destroyed by the imposed pressure wave. The skin friction drag can be reduced significantly by the supercavity and the form drag can also be affected by the cavitation evolution. References 1 Kunz R F, Boger D A & Stinebring D R, A preconditioned Navier-Stokes method for two-phase flows with application to cavitation prediction, Comput. Fluid, 29 (2000) Saurel R, Petitpas F & Berry R A, Simple and efficient relaxation methods for interfaces separating compressible fluids, cavitating flows and shocks in multiphase mixtures, J Comput Phys, 228 (2009) Liu T G, Khoo B C & Xie W F, Isentropic one-fluid modeling of unsteady cavitating flow, J Comput Phys, 201 (2004) Zheng J G, Khoo B C & Hu Z M, Simulation of Wave-Flow- Cavitation Interaction Using a Compressible Homogenous Flow Method, Commun Comput Phys, 14(2) (2013) Wallis G B, One-dimensional two-phase flow, (McGill-Hill, New York) Van Leer B, Upwind and High-Resolution Methods for Compressible Flow: From Donor Cell to Residual- Distribution Schemes, Commun Comput Phys, 1(2) (2006) Toro E F, Riemann solvers and numerical methods for fluid Dynamics: A practical introduction, (Springer-Verlag, Berlin Heidelberg) Blazek J, Computational fluid dynamics: principles and applications, (ELSEVIER) Hrubes J D, High-speed imaging of supercavitating underwater projectiles, Exp. Fluids, 30 (2001)
Application of a Laser Induced Fluorescence Model to the Numerical Simulation of Detonation Waves in Hydrogen-Oxygen-Diluent Mixtures
Supplemental material for paper published in the International J of Hydrogen Energy, Vol. 30, 6044-6060, 2014. http://dx.doi.org/10.1016/j.ijhydene.2014.01.182 Application of a Laser Induced Fluorescence
More informationEFFECT OF LIQUID PHASE COMPRESSIBILITY ON MODELING OF GAS-LIQUID TWO-PHASE FLOWS USING TWO-FLUID MODEL
EFFECT OF LIQUID PHASE COMPRESSIBILITY ON MODELING OF GAS-LIQUID TWO-PHASE FLOWS USING TWO-FLUID MODEL Vahid SHOKRI 1*,Kazem ESMAEILI 2 1,2 Department of Mechanical Engineering, Sari Branch, Islamic Azad
More informationCavitation modeling using compressible Navier Stokes and Korteweg equations
Computational Methods in Multiphase Flow VIII 425 Cavitation modeling using compressible Navier Stokes and Korteweg equations T. Ait Ali 1, S. Khelladi 1, L. Ramirez 2 & X. Nogueira 2 1 DynFluid Laboratory,
More informationNumerical investigation of cavitation-regimes in a converging-diverging nozzle
Numerical investigation of cavitation-regimes in a converging-diverging nozzle 1 Polina Gorkh, 1 Steffen J. Schmidt, and 1 Nikolaus A. Adams 1 Institute of Aerodynamics and Fluid Mechanics, Department
More informationNUMERICAL SIMULATION OF TRANSIENT SLURRY-CAVITATED MULTIPHASE FLOWS
ISTP-1, 005, PRAGUE 1 TH INTERNATIONAL SYMPOSIUM ON TRANSPO PHENOMENA NUMERICAL SIMULATION OF TRANSIENT SLURRY-CAVITATED MULTIPHASE FLOWS Yang-Yao Niu and Yee-Ming Lin Institute of Mechanical and Aerospace
More informationA Multi-Dimensional Limiter for Hybrid Grid
APCOM & ISCM 11-14 th December, 2013, Singapore A Multi-Dimensional Limiter for Hybrid Grid * H. W. Zheng ¹ 1 State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy
More informationComparison of Two Computational Methods for High-Velocity Cavitating Flows around Conical Projectiles : OTi-HULL Hydrocode and Two-Phase Flow Method
CAV200:sessionB3.003 - -Schmitt Page 29/03/0 Comparison of Two Computational Methods for High-Velocity Cavitating Flows around Conical Projectiles : OTi-HULL Hydrocode and Two-Phase Flow Method M. Schaffar
More informationA Numerical Study of Compressible Two-Phase Flows Shock and Expansion Tube Problems
A Numerical Study of Compressible Two-Phase Flows Shock and Expansion Tube Problems Dia Zeidan,a) and Eric Goncalves 2 School of Basic Sciences and Humanities, German Jordanian University, Amman, Jordan
More informationTHE numerical simulation of the creation and evolution
Proceedings of the World Congress on Engineering Vol III WCE, July 4-6,, London, U.K. Numerical Simulation of Compressible Two-phase Flows Using an Eulerian Type Reduced Model A. Ballil, Member, IAENG,
More informationInvestigation of an implicit solver for the simulation of bubble oscillations using Basilisk
Investigation of an implicit solver for the simulation of bubble oscillations using Basilisk D. Fuster, and S. Popinet Sorbonne Universités, UPMC Univ Paris 6, CNRS, UMR 79 Institut Jean Le Rond d Alembert,
More informationTransport equation cavitation models in an unstructured flow solver. Kilian Claramunt, Charles Hirsch
Transport equation cavitation models in an unstructured flow solver Kilian Claramunt, Charles Hirsch SHF Conference on hydraulic machines and cavitation / air in water pipes June 5-6, 2013, Grenoble, France
More informationNumerical investigation of shock wave dense particles cloud interaction
25 th ICDERS August 2 7, 215 Leeds, UK Numerical investigation of shock wave dense particles cloud interaction Pavel S. Utkin Moscow Institute of Physics and Technology Dolgoprudny, Moscow region, Russia
More information1. Introduction Some Basic Concepts
1. Introduction Some Basic Concepts 1.What is a fluid? A substance that will go on deforming in the presence of a deforming force, however small 2. What Properties Do Fluids Have? Density ( ) Pressure
More informationWave propagation methods for hyperbolic problems on mapped grids
Wave propagation methods for hyperbolic problems on mapped grids A France-Taiwan Orchid Project Progress Report 2008-2009 Keh-Ming Shyue Department of Mathematics National Taiwan University Taiwan ISCM
More informationExperimental and numerical study of the initial stages in the interaction process between a planar shock wave and a water column
Experimental and numerical study of the initial stages in the interaction process between a planar shock wave and a water column Dan Igra and Kazuyoshi Takayama Shock Wave Research Center, Institute of
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationIntroduction to Aerodynamics. Dr. Guven Aerospace Engineer (P.hD)
Introduction to Aerodynamics Dr. Guven Aerospace Engineer (P.hD) Aerodynamic Forces All aerodynamic forces are generated wither through pressure distribution or a shear stress distribution on a body. The
More informationMULTIGRID CALCULATIONS FOB. CASCADES. Antony Jameson and Feng Liu Princeton University, Princeton, NJ 08544
MULTIGRID CALCULATIONS FOB. CASCADES Antony Jameson and Feng Liu Princeton University, Princeton, NJ 0544 1. Introduction Development of numerical methods for internal flows such as the flow in gas turbines
More informationTarget Simulations. Roman Samulyak in collaboration with Y. Prykarpatskyy, T. Lu
Muon Collider/Neutrino Factory Collaboration Meeting May 26 28, CERN, Geneva U.S. Department of Energy Target Simulations Roman Samulyak in collaboration with Y. Prykarpatskyy, T. Lu Center for Data Intensive
More informationCavitation Induction by Projectile Impacting on a Water Jet
Cavitation Induction by Projectile Impacting on a Water Jet 1 E. Stavropoulos Vasilakis*; 1 P. Koukouvinis; 2 M. Farhat; 1 M. Gavaises; 1 School of Mathematics, Computer Science and Engineering, CITY,
More informationValidation 3. Laminar Flow Around a Circular Cylinder
Validation 3. Laminar Flow Around a Circular Cylinder 3.1 Introduction Steady and unsteady laminar flow behind a circular cylinder, representing flow around bluff bodies, has been subjected to numerous
More informationA drop forms when liquid is forced out of a small tube. The shape of the drop is determined by a balance of pressure, gravity, and surface tension
A drop forms when liquid is forced out of a small tube. The shape of the drop is determined by a balance of pressure, gravity, and surface tension forces. 2 Objectives 3 i i 2 1 INTRODUCTION Property:
More informationDetailed Outline, M E 320 Fluid Flow, Spring Semester 2015
Detailed Outline, M E 320 Fluid Flow, Spring Semester 2015 I. Introduction (Chapters 1 and 2) A. What is Fluid Mechanics? 1. What is a fluid? 2. What is mechanics? B. Classification of Fluid Flows 1. Viscous
More informationA Computational Study of High-Speed Droplet Impact
Copyright 2011 Tech Science Press FDMP, vol.7, no.4, pp.329-340, 2011 A Computational Study of High-Speed Droplet Impact T. Sanada 1, K. Ando 2 and T. Colonius 2 Abstract: When a droplet impacts a solid
More informationProjection Dynamics in Godunov-Type Schemes
JOURNAL OF COMPUTATIONAL PHYSICS 142, 412 427 (1998) ARTICLE NO. CP985923 Projection Dynamics in Godunov-Type Schemes Kun Xu and Jishan Hu Department of Mathematics, Hong Kong University of Science and
More informationTurbulent Boundary Layers & Turbulence Models. Lecture 09
Turbulent Boundary Layers & Turbulence Models Lecture 09 The turbulent boundary layer In turbulent flow, the boundary layer is defined as the thin region on the surface of a body in which viscous effects
More informationOn the Modeling and Simulation of a Laser-Induced Caviation Bubble
1 On the Modeling and Simulation of a Laser-Induced Caviation Bubble Ali Zein Supervisor: Prof. Dr. Gerald Warnecke Co-supervisor: Dr. Maren Hantke Institut für Analysis und Numerik Otto-von-Guericke-Universität
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More information[N175] Development of Combined CAA-CFD Algorithm for the Efficient Simulation of Aerodynamic Noise Generation and Propagation
The 32nd International Congress and Exposition on Noise Control Engineering Jeju International Convention Center, Seogwipo, Korea, August 25-28, 2003 [N175] Development of Combined CAA-CFD Algorithm for
More informationn v molecules will pass per unit time through the area from left to
3 iscosity and Heat Conduction in Gas Dynamics Equations of One-Dimensional Gas Flow The dissipative processes - viscosity (internal friction) and heat conduction - are connected with existence of molecular
More information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationPart 3. Stability and Transition
Part 3 Stability and Transition 281 Overview T. Cebeci 1 Recent interest in the reduction of drag of underwater vehicles and aircraft components has rekindled research in the area of stability and transition.
More informationDetailed Outline, M E 521: Foundations of Fluid Mechanics I
Detailed Outline, M E 521: Foundations of Fluid Mechanics I I. Introduction and Review A. Notation 1. Vectors 2. Second-order tensors 3. Volume vs. velocity 4. Del operator B. Chapter 1: Review of Basic
More informationLiquids and solids are essentially incompressible substances and the variation of their density with pressure is usually negligible.
Properties of Fluids Intensive properties are those that are independent of the mass of a system i.e. temperature, pressure and density. Extensive properties are those whose values depend on the size of
More informationEffect of Applied Magnetic Field on Shock Boundary Layer Interaction
University of Kentucky UKnowledge Mechanical Engineering Faculty Publications Mechanical Engineering 1-2012 Effect of Applied Magnetic Field on Shock Boundary Layer Interaction Ovais U. Khan University
More informationImprovements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme
Improvements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme Oren Peles and Eli Turkel Department of Applied Mathematics, Tel-Aviv University In memoriam of
More informationSimulation of unsteady muzzle flow of a small-caliber gun
Advances in Fluid Mechanics VI 165 Simulation of unsteady muzzle flow of a small-caliber gun Y. Dayan & D. Touati Department of Computational Mechanics & Ballistics, IMI, Ammunition Group, Israel Abstract
More informationConstraints on Equation of State for Cavitating Flows with Thermodynamic Effects
Constraints on Equation of State for Cavitating Flows with Thermodynamic Effects Eric Goncalvès da Silva, Régiane Fortes Patella To cite this version: Eric Goncalvès da Silva, Régiane Fortes Patella. Constraints
More informationBoundary-Layer Theory
Hermann Schlichting Klaus Gersten Boundary-Layer Theory With contributions from Egon Krause and Herbert Oertel Jr. Translated by Katherine Mayes 8th Revised and Enlarged Edition With 287 Figures and 22
More informationNumerical studies of real-gas effects on two-dimensional hypersonic shock-wave/boundary-layer interaction
Numerical studies of real-gas effects on two-dimensional hypersonic shock-wave/boundary-layer interaction Gregory H. Furumoto, a) Xiaolin Zhong, and John C. Skiba Department of Mechanical and Aerospace
More informationHybrid RANS/LES simulations of a cavitating flow in Venturi
Hybrid RANS/LES simulations of a cavitating flow in Venturi TSFP-8, 28-31 August 2013 - Poitiers - France Jean Decaix jean.decaix@hevs.ch University of Applied Sciences, Western Switzerland Eric Goncalves
More informationSimulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions
Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,
More informationComputational Fluid Dynamics Study Of Fluid Flow And Aerodynamic Forces On An Airfoil S.Kandwal 1, Dr. S. Singh 2
Computational Fluid Dynamics Study Of Fluid Flow And Aerodynamic Forces On An Airfoil S.Kandwal 1, Dr. S. Singh 2 1 M. Tech Scholar, 2 Associate Professor Department of Mechanical Engineering, Bipin Tripathi
More informationSimulating Interfacial Tension of a Falling. Drop in a Moving Mesh Framework
Simulating Interfacial Tension of a Falling Drop in a Moving Mesh Framework Anja R. Paschedag a,, Blair Perot b a TU Berlin, Institute of Chemical Engineering, 10623 Berlin, Germany b University of Massachusetts,
More informationTo study the motion of a perfect gas, the conservation equations of continuity
Chapter 1 Ideal Gas Flow The Navier-Stokes equations To study the motion of a perfect gas, the conservation equations of continuity ρ + (ρ v = 0, (1.1 momentum ρ D v Dt = p+ τ +ρ f m, (1.2 and energy ρ
More informationA Simple Model for Cavitation with Non-condensable Gases
A Simple Model for Cavitation with Non-condensable Gases Mathieu Bachmann, Siegfried Müller, Philippe Helluy, Hélène Mathis To cite this version: Mathieu Bachmann, Siegfried Müller, Philippe Helluy, Hélène
More informationCFD SIMULATIONS OF FLOW, HEAT AND MASS TRANSFER IN THIN-FILM EVAPORATOR
Distillation Absorption 2010 A.B. de Haan, H. Kooijman and A. Górak (Editors) All rights reserved by authors as per DA2010 copyright notice CFD SIMULATIONS OF FLOW, HEAT AND MASS TRANSFER IN THIN-FILM
More informationFundamentals of Fluid Mechanics
Sixth Edition Fundamentals of Fluid Mechanics International Student Version BRUCE R. MUNSON DONALD F. YOUNG Department of Aerospace Engineering and Engineering Mechanics THEODORE H. OKIISHI Department
More informationCENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer
CENG 501 Examination Problem: Estimation of Viscosity with a Falling - Cylinder Viscometer You are assigned to design a fallingcylinder viscometer to measure the viscosity of Newtonian liquids. A schematic
More informationPhysical Diffusion Cures the Carbuncle Phenomenon
Physical Diffusion Cures the Carbuncle Phenomenon J. M. Powers 1, J. Bruns 1, A. Jemcov 1 1 Department of Aerospace and Mechanical Engineering University of Notre Dame, USA Fifty-Third AIAA Aerospace Sciences
More informationNUMERICAL SIMULATION OF CAVITATING FLOW AROUND A 2D HYDROFOIL: A BAROTROPIC APPROACH
Cav03-OS--005 Fifth International Symposium on Cavitation (CAV003) Osaka, Japan, November -4, 003 NUMERICAL SIMULATION OF CAVITATING FLOW AROUND A D HYDROFOIL: A BAROTROPIC APPROACH B. Pouffary Laboratoire
More informationFluid Mechanics. Spring 2009
Instructor: Dr. Yang-Cheng Shih Department of Energy and Refrigerating Air-Conditioning Engineering National Taipei University of Technology Spring 2009 Chapter 1 Introduction 1-1 General Remarks 1-2 Scope
More informationNUMERICAL SIMULATION OF UNSTEADY CAVITATING FLOWS
NUMERICAL SIMULATION OF UNSTEADY CAVITATING FLOWS Charles C.S. Song and Qiao Qin St. Anthony Falls Laboratory, University of Minnesota Mississippi River at 3 rd Ave. SE, Minneapolis, MN 55414, USA ABSTRACT
More informationCavitation Control on Hydrofoils
Proceedings of the International Conference on Heat Transfer and Fluid Flow Prague, Czech Republic, August 11-12, 2014 Paper No. 181 Cavitation Control on Hydrofoils Mohammad Mortezazadeh, Ali Katal, Khodayar
More informationMass flow determination in flashing openings
Int. Jnl. of Multiphysics Volume 3 Number 4 009 40 Mass flow determination in flashing openings Geanette Polanco Universidad Simón Bolívar Arne Holdø Narvik University College George Munday Coventry University
More informationWilliam В. Brower, Jr. A PRIMER IN FLUID MECHANICS. Dynamics of Flows in One Space Dimension. CRC Press Boca Raton London New York Washington, D.C.
William В. Brower, Jr. A PRIMER IN FLUID MECHANICS Dynamics of Flows in One Space Dimension CRC Press Boca Raton London New York Washington, D.C. Table of Contents Chapter 1 Fluid Properties Kinetic Theory
More informationENGINEERING FLUID MECHANICS. CHAPTER 1 Properties of Fluids
CHAPTER 1 Properties of Fluids ENGINEERING FLUID MECHANICS 1.1 Introduction 1.2 Development of Fluid Mechanics 1.3 Units of Measurement (SI units) 1.4 Mass, Density, Specific Weight, Specific Volume, Specific
More informationHydraulic Modelling for Drilling Automation
Hydraulic Modelling for Drilling Automation CASA Day Harshit Bansal April 19, 2017 Where innovation starts Team: Supervisors at TU/e : W.H.A. Schilders, N. van de Wouw, B. Koren, L. Iapichino Collaborators:
More informationURANS Computations of Cavitating Flow around a 2-D Wedge by Compressible Two-Phase Flow Solver
URANS Computations of Cavitating Flow around a 2-D Wedge by Compressible Two-Phase Flow Solver *Yohan Choe 1), Hyeongjun Kim 1), Chongam Kim 2) 1), 2) Department of Aerospace Engineering, Seoul National
More informationLecture1: Characteristics of Hypersonic Atmosphere
Module 1: Hypersonic Atmosphere Lecture1: Characteristics of Hypersonic Atmosphere 1.1 Introduction Hypersonic flight has special traits, some of which are seen in every hypersonic flight. Presence of
More informationApplied Fluid Mechanics
Applied Fluid Mechanics 1. The Nature of Fluid and the Study of Fluid Mechanics 2. Viscosity of Fluid 3. Pressure Measurement 4. Forces Due to Static Fluid 5. Buoyancy and Stability 6. Flow of Fluid and
More information2.The lines that are tangent to the velocity vectors throughout the flow field are called steady flow lines. True or False A. True B.
CHAPTER 03 1. Write Newton's second law of motion. YOUR ANSWER: F = ma 2.The lines that are tangent to the velocity vectors throughout the flow field are called steady flow lines. True or False 3.Streamwise
More informationDirect numerical simulation of compressible multi-phase ow with a pressure-based method
Seventh International Conference on Computational Fluid Dynamics (ICCFD7), Big Island, Hawaii, July 9-13, 2012 ICCFD7-2902 Direct numerical simulation of compressible multi-phase ow with a pressure-based
More informationGODUNOV-TYPE SOLUTIONS FOR TWO-PHASE WATER HAMMER FLOWS
GODUNOV-TYPE SOLUTIONS FOR TWO-PHASE WATER HAMMER FLOWS ARTURO S. LEON Dept. of Civil and Envir. Engng., Univ. of Illinois at Urbana-Champaign, 2519 Hydrosystems Lab., MC-250. 205 North Mathews Av., Urbana,
More informationWALL ROUGHNESS EFFECTS ON SHOCK BOUNDARY LAYER INTERACTION FLOWS
ISSN (Online) : 2319-8753 ISSN (Print) : 2347-6710 International Journal of Innovative Research in Science, Engineering and Technology An ISO 3297: 2007 Certified Organization, Volume 2, Special Issue
More informationIntroduction and Basic Concepts
Topic 1 Introduction and Basic Concepts 1 Flow Past a Circular Cylinder Re = 10,000 and Mach approximately zero Mach = 0.45 Mach = 0.64 Pictures are from An Album of Fluid Motion by Van Dyke Flow Past
More informationNUMERICAL OPTIMIZATION OF THE SHAPE OF A HOLLOW PROJECTILE
NUMERICAL OPTIMIZATION OF THE SHAPE OF A HOLLOW PROJECTILE Wessam Mahfouz Elnaggar, Zhihua Chen and Hui Zhang Key Laboratory of Transient Physics, Nanjing University of Science and Technology, Nanjing,
More informationDepartment of Mechanical Engineering
Department of Mechanical Engineering AMEE401 / AUTO400 Aerodynamics Instructor: Marios M. Fyrillas Email: eng.fm@fit.ac.cy HOMEWORK ASSIGNMENT #2 QUESTION 1 Clearly there are two mechanisms responsible
More informationHigh Speed Aerodynamics. Copyright 2009 Narayanan Komerath
Welcome to High Speed Aerodynamics 1 Lift, drag and pitching moment? Linearized Potential Flow Transformations Compressible Boundary Layer WHAT IS HIGH SPEED AERODYNAMICS? Airfoil section? Thin airfoil
More informationApplication of the immersed boundary method to simulate flows inside and outside the nozzles
Application of the immersed boundary method to simulate flows inside and outside the nozzles E. Noël, A. Berlemont, J. Cousin 1, T. Ménard UMR 6614 - CORIA, Université et INSA de Rouen, France emeline.noel@coria.fr,
More informationFALLING FILM FLOW ALONG VERTICAL PLATE WITH TEMPERATURE DEPENDENT PROPERTIES
Proceedings of the International Conference on Mechanical Engineering 2 (ICME2) 8-2 December 2, Dhaka, Bangladesh ICME-TH-6 FALLING FILM FLOW ALONG VERTICAL PLATE WITH TEMPERATURE DEPENDENT PROPERTIES
More informationDetailed Numerical Simulation of Liquid Jet in Cross Flow Atomization: Impact of Nozzle Geometry and Boundary Condition
ILASS-Americas 25th Annual Conference on Liquid Atomization and Spray Systems, Pittsburgh, PA, May 23 Detailed Numerical Simulation of Liquid Jet in Cross Flow Atomization: Impact of Nozzle Geometry and
More informationBLUFF-BODY AERODYNAMICS
International Advanced School on WIND-EXCITED AND AEROELASTIC VIBRATIONS OF STRUCTURES Genoa, Italy, June 12-16, 2000 BLUFF-BODY AERODYNAMICS Lecture Notes by Guido Buresti Department of Aerospace Engineering
More informationModeling full-mach-range cavitating flow with sharp interface model
Center for Turbulence Research Proceedings of the Summer Program 28 83 Modeling full-mach-range cavitating flow with sharp interface model By X. Y. Hu, N. A. Adams, E. Johnsen AND G. Iaccarino In this
More informationNumerical Investigation of the Transonic Base Flow of A Generic Rocket Configuration
1 Numerical Investigation of the Transonic Base Flow of A Generic Rocket Configuration A. Henze, C. Glatzer, M. Meinke, W. Schröder Institute of Aerodynamics, RWTH Aachen University, Germany March 21,
More informationAPPLICATION OF HYBRID CFD/CAA TECHNIQUE FOR MODELING PRESSURE FLUCTUATIONS IN TRANSONIC FLOWS
TASK QUARTERLY Vol. 17, Nos 3 4, 2013, pp. 145 154 APPLICATION OF HYBRID CFD/CAA TECHNIQUE FOR MODELING PRESSURE FLUCTUATIONS IN TRANSONIC FLOWS SŁAWOMIR DYKAS, WŁODZIMIERZ WRÓBLEWSKI AND SEBASTIAN RULIK
More informationInterpreting Differential Equations of Transport Phenomena
Interpreting Differential Equations of Transport Phenomena There are a number of techniques generally useful in interpreting and simplifying the mathematical description of physical problems. Here we introduce
More informationNumerical Simulation of Cavitating Flows with Homogeneous Models
Numerical Simulation of Cavitating Flows with Homogeneous Models Eric Goncalvès da Silva, Régiane Fortes Patella To cite this version: Eric Goncalvès da Silva, Régiane Fortes Patella. Numerical Simulation
More informationA NEW MODELING OF SHEET CAVITATION CONSIDERING THE THERMODYNAMIC EFFECTS. Yumiko Sekino Ebara Cooporation, Ltd., Ota-ku, Tokyo JAPAN
Cav3-GS-6-3 Fifth International Symposium on Cavitation (CAV3) Osaka, Japan, November -4, 3 A NEW MODELING OF SHEET CAVITATION CONSIDERING THE THERMODYNAMIC EFFECTS Takashi Tokumasu Institute of Fluid
More informationIntroduction to Heat and Mass Transfer. Week 12
Introduction to Heat and Mass Transfer Week 12 Next Topic Convective Heat Transfer» Heat and Mass Transfer Analogy» Evaporative Cooling» Types of Flows Heat and Mass Transfer Analogy Equations governing
More informationFluid Flow, Heat Transfer and Boiling in Micro-Channels
L.P. Yarin A. Mosyak G. Hetsroni Fluid Flow, Heat Transfer and Boiling in Micro-Channels 4Q Springer 1 Introduction 1 1.1 General Overview 1 1.2 Scope and Contents of Part 1 2 1.3 Scope and Contents of
More informationContents. I Introduction 1. Preface. xiii
Contents Preface xiii I Introduction 1 1 Continuous matter 3 1.1 Molecules................................ 4 1.2 The continuum approximation.................... 6 1.3 Newtonian mechanics.........................
More informationFluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition
Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition C. Pozrikidis m Springer Contents Preface v 1 Introduction to Kinematics 1 1.1 Fluids and solids 1 1.2 Fluid parcels and flow
More informationNonlinear Wave Theory for Transport Phenomena
JOSO 2016 March 9-11 2015 Nonlinear Wave Theory for Transport Phenomena ILYA PESHKOV CHLOE, University of Pau, France EVGENIY ROMENSKI Sobolev Institute of Mathematics, Novosibirsk, Russia MICHAEL DUMBSER
More informationEntropy stable schemes for compressible flows on unstructured meshes
Entropy stable schemes for compressible flows on unstructured meshes Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore deep@math.tifrbng.res.in http://math.tifrbng.res.in/
More informationCFD in COMSOL Multiphysics
CFD in COMSOL Multiphysics Mats Nigam Copyright 2016 COMSOL. Any of the images, text, and equations here may be copied and modified for your own internal use. All trademarks are the property of their respective
More informationNumerical Investigation of Cavitation Bubble Collapsing Behavior
First International Smposium on Marine Propulsors SMP 09, Trondheim, Norwa, June 009 Numerical Investigation of Cavitation Bubble Collapsing Behavior Beong Rog Shin 1 1 Deptartment of Mechanical Engineering,
More informationChapter 1: Basic Concepts
What is a fluid? A fluid is a substance in the gaseous or liquid form Distinction between solid and fluid? Solid: can resist an applied shear by deforming. Stress is proportional to strain Fluid: deforms
More information3. FORMS OF GOVERNING EQUATIONS IN CFD
3. FORMS OF GOVERNING EQUATIONS IN CFD 3.1. Governing and model equations in CFD Fluid flows are governed by the Navier-Stokes equations (N-S), which simpler, inviscid, form is the Euler equations. For
More informationCompressible Flow. Professor Ugur GUVEN Aerospace Engineer Spacecraft Propulsion Specialist
Compressible Flow Professor Ugur GUVEN Aerospace Engineer Spacecraft Propulsion Specialist What is Compressible Flow? Compressible Flow is a type of flow in which the density can not be treated as constant.
More informationFLUID MECHANICS. Chapter 9 Flow over Immersed Bodies
FLUID MECHANICS Chapter 9 Flow over Immersed Bodies CHAP 9. FLOW OVER IMMERSED BODIES CONTENTS 9.1 General External Flow Characteristics 9.3 Drag 9.4 Lift 9.1 General External Flow Characteristics 9.1.1
More information58:160 Intermediate Fluid Mechanics Bluff Body Professor Fred Stern Fall 2014
Professor Fred Stern Fall 04 Chapter 7 Bluff Body Fluid flows are broadly categorized:. Internal flows such as ducts/pipes, turbomachinery, open channel/river, which are bounded by walls or fluid interfaces:
More informationModeling Cavitation over Axisymmetric Bodies:
Modeling Cavitation over Axisymmetric Bodies: VOF Technique versus Boundary Element Method H. Moin, I. Rashidi, Mohammad P. Fard 3, Mahmoud P. Fard 4 and E. Roohi 5,,3,4 Department of Mechanical Engineering,
More informationModelling of single bubble-dynamics and thermal effects
Journal of Physics: Conference Series PAPER OPEN ACCESS Modelling of single bubble-dynamics and thermal effects To cite this article: D Papoulias and M Gavaises 205 J. Phys.: Conf. Ser. 656 02098 View
More informationEigenmode Analysis of Boundary Conditions for the One-dimensional Preconditioned Euler Equations
NASA/CR-1998-208741 ICASE Report No. 98-51 Eigenmode Analysis of Boundary Conditions for the One-dimensional Preconditioned Euler Equations David L. Darmofal Massachusetts Institute of Technology, Cambridge,
More informationRecent advances in numerical methods for compressible two-phase flow with heat & mass transfers
Recent advances in numerical methods for compressible two-phase flow with heat & mass transfers Keh-Ming Shyue Institute of Applied Mathematical Sciences National Taiwan University Joint work with Marica
More informationNumerical Studies of Droplet Deformation and Break-up
ILASS Americas 14th Annual Conference on Liquid Atomization and Spray Systems, Dearborn, MI, May 2001 Numerical Studies of Droplet Deformation and Break-up B. T. Helenbrook Department of Mechanical and
More informationWhere does Bernoulli's Equation come from?
Where does Bernoulli's Equation come from? Introduction By now, you have seen the following equation many times, using it to solve simple fluid problems. P ρ + v + gz = constant (along a streamline) This
More informationNUMERICAL INVESTIGATION OF PERFORMANCE OF KAPLAN TURBINE WITH DRAFT TUBE
NUMERICAL INVESTIGATION OF PERFORMANCE OF KAPLAN TURBINE WITH DRAFT TUBE Mohamed Adel and Nabil H. Mostafa ABSTRACT Mechanical Power Engineering Department, Engineering Faculty Zagazig University, Zagazig,
More informationNumerical Analysis of Unsteady Open Water Characteristics of Surface Piercing Propeller
Third International Symposium on Marine Propulsors smp 13, Launceston, Tasmania, Australia, May 2013 Numerical Analysis of Unsteady Open Water Characteristics of Surface Piercing Propeller Kohei Himei
More information