DOCSLIB.ORG

CHAPTER 2 Viscous Compressible Fluid Dynamics

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
CHAPTER 2 Viscous Compressible Fluid Dynamics CHAPTER 2 Viscous compressible fluid dynamics L. Škerget & M. Hriberšek Faculty of Mechanical Engineering, University of Maribor, Smetanova, Maribor, Slovenia. Abstract The development of boundary element methods for computation of fluid flow is predominantly focused on the flow of incompressible fluids. This assumption pro- vides a good approximation for a wide variety of engineering applications, but to achieve a better representation of physical phenomena, a compressible fluid flow must be taken into account. The goal of the chapter is, first, to give an in depth expla- nation for the transformation of Navier–Stokes equations for compressible fluids into its velocity–vorticity formulation equivalent. In addition, the pressure equation for the velocity–vorticity formulation is derived. This is followed by integral trans- formations of the governing equations based on the use of a parabolic diffusion fundamental solution. Section 8 explains discrete models, including the subdomain type discretisation. As a test example, natural convection in a differentially heated tall enclosure is presented. 1 Introduction Most of the studies dealing with transport phenomena are based on presuming that the fluid is incompressible and viscous, where the mass density is a constant quantity, and the velocity does not depend on the mass density. Pressure in the incompressible fluid flow model is not a thermodynamic state variable, but simply a force in the linear momentum balance equation. Such an easy rheological model for the fluid is suitable for modelling of slow flows, or flows with small pressure and temperature gradients or no chemical reaction and where, therefore, the mass density differences may be neglected. WIT Transactions on State of the Art in Science and Engineering, Vol 14, © 2007 WIT Press www.witpress.com, ISSN 1755-8336 (on-line) doi:10.2495/978-1-84564-100-9/02 34 Domain Decomposition Techniques for Boundary Elements In this chapter we will deal with the real compressible viscous fluid with restric- tion to the subsonic flows where the difference in mass density significantly influ- ences the vector velocity field. However, there are no shock waves and no sudden sharp changes in the values of the field functions. In the model of compressible fluid the pressure is a thermodynamic quantity p = p(ρ , T). 2 Conservation equations The analytical description of the motion of a continuous viscous compressible media is based on the conservation of mass, momentum, and energy with asso- ciated rheological models and equations of state [1]. The present development is focused on the laminar flow of compressible isotropic fluid in solution region bounded by boundary . The field functions of interest are the velocity vector field vi(rj, t), the scalar pressure field p(rj, t), the temperature field T(rj, t) and the field of mass density ρ(rj, t). The mass, momentum, and energy conservation equations are given by ∂ρ ∂ρvj + = 0, (1) ∂t ∂xj ∂ρvi ∂ρvivj ∂σij + = + ρgi,(2) ∂t ∂xj ∂xj ∂ρT ∂ρTvj ∂qj cp + =− + S + , (3) ∂t ∂xj ∂xj in the Cartesian frame xi, where ρ and cp denote variable fluid mass density and isobaric specific heat capacity per unit mass, t is time, gi is gravitational acceleration vector, σij represents the components of the total stress tensor, qi is specific heat diffusion flux, while S stands for the heat source term and is a Rayleigh viscous dissipation function which stands for the conversion of mechanical energy to heat and acts as an additional heat source. With the definition of Stokes material derivative of the variable (·), as given by D (·) /Dt = ∂ (·) /∂t + vk∂ (·) /∂xk, the Navier–Stokes equations are given by ∂vj 1 Dρ + = 0, (4) ∂xj ρ Dt Dvi ∂σij ρ = + ρgi, (5) Dt ∂xj DT ∂qj c =− + S + , (6) Dt ∂xj where c denotes specific heat per unit volume, c = cp ρ. If the incompressible fluid model is applied, than the total rate of mass density variation is identically zero Dρ =0, (7) Dt WIT Transactions on State of the Art in Science and Engineering, Vol 14, © 2007 WIT Press www.witpress.com, ISSN 1755-8336 (on-line) Viscous Compressible Fluid Dynamics 35 and the mass conservation is satisfied with a simple restriction condition for the solenoidal velocity vector ∂vj =0. (8) ∂xj The set of eqns (4)–(6) represents an unclosed system of partial differential equations that has to be closed and solved in conjuction with appropriate rheological equations, equations of state and boundary, as well as initial conditions of the problem. The Cauchy total stress tensor σij can be decomposed into a pressure contribution p plus an extra deviatoric stress tensor field function τij σij =−pδij + τij, (9) where δij is the Kronecker delta. The Rayleigh viscous dissipation term is given by ∂vi = τij . (10) ∂xj 3 Linear gradient type of constitutive models In the general Reiner–Rivlins rheological model of viscous shear fluid the deviatoric stress tensor reads as τij = αδij + βε˙ij + γ ε˙ikε˙kj. (11) By considering the equalities α =−2ηε˙ii/3 and β = 2η, the following form of the constitutive model for compressible viscous shear fluid may be written 2 τij = 2ηε˙ij − η Dδij, (12) 3 where D = div v =˙εii represents the divergence of the velocity field. Considering the eqn (10) we may write ∂v ∂v ∂vj ∂v 2 = η i i + i − ηD2. (13) ∂xj ∂xj ∂xi ∂xj 3 In the case of intensive unsteady heat transfer, it is important to take into account a terminal velocity of a moving temperature front, namely ∂T ∂qi qi =−k − λ , (14) ∂xi ∂t where material constants k and λ are the heat conductivity and the heat relaxation time. For most heat transfer problems of practical importance, the simplification known as the Fourier law of heat diffusion is accurate enough, namely ∂T qi =−k . (15) ∂xi WIT Transactions on State of the Art in Science and Engineering, Vol 14, © 2007 WIT Press www.witpress.com, ISSN 1755-8336 (on-line) 36 Domain Decomposition Techniques for Boundary Elements 4 Primitive variables formulation Combining constitutive models for stress tensor or momentum flux tensor and heat diffusion flux, eqns (12) and (15) in conservation eqns (5) and (6), the following system of nonlinear equations is developed ∂vj 1 Dρ = D =− , (16) ∂xj ρ Dt Dvi ∂ 2 ∂p ρ = 2η ε˙ij − η Dδij − + ρgi, (17) Dt ∂xj 3 ∂xi DT ∂ ∂T c = k + S + . (18) Dt ∂xj ∂xj Because of analytical reasons in developing the velocity–vorticity formulation of governing equations, the operator div τ is worth writing in an extended form ∂ 2 div τ = η Lij + Lji − Dδij ∂xj 3 ∂η ∂Lij ∂Lji 2 ∂ηD = 2 ε˙ij + η + η − , (19) ∂xj ∂xj ∂xj 3 ∂xi where the term ∂Lji/∂xj is equal to the grad D, due to the continuity equation, div v = D. Thus the following relation is valid ∂Lji ∂ ∂vj ∂D = = . (20) ∂xj ∂xi ∂xj ∂xi ˙ By considering the equality between tensor field functions ε˙ij = Lij − ij, one can derive an expression ∂η ∂η ˙ ∂Lij ∂D 2 ∂ηD div τ = 2 Lij − 2 ij + η + η − . (21) ∂xj ∂xj ∂xj ∂xi 3 ∂xi Finally, substituting the equalities 2 ∂Lij ∂ vi ∂η ˙ ∂η = and 2 ij =−eijk ωk, (22) ∂xj ∂xj∂xj ∂xj ∂xj the first extended form of div τ may be formulated as follows 2 div τ = η v + grad η × ω + 2 grad v · grad η + η grad D − grad (ηD). (23) 3 By the derivation of the second extended form we also take into account the equalities div (grad v ) = grad D − rot ω = v, (24) rot (ηω ) = η rot ω − ω × grad η, (25) WIT Transactions on State of the Art in Science and Engineering, Vol 14, © 2007 WIT Press www.witpress.com, ISSN 1755-8336 (on-line) Viscous Compressible Fluid Dynamics 37 obtaining the second extended form 4 div τ =−rot(ηω ) + 2grad η × ω + 2grad v · grad η + grad (ηD) 3 − 2D grad η. (26) Considering the second extended form of the term div τ , eqn (26), the momentum equation (17) may be written in a form suitable to derive the velocity–vorticity formulation, e.g. in a vector form Dv 4 ρ =−rot(ηω ) + 2grad η × ω + 2grad v · grad η + grad (ηD) − 2D grad η Dt 3 − grad p + ρg , (27) or in Cartesian tensor formulation Dvi ∂ηωk ∂η ∂η ∂vi 4 ∂ηD ∂η ρ =−eijk + 2eijk ωk + 2 + − 2D Dt ∂xj ∂xj ∂xj ∂xj 3 ∂xi ∂xi ∂p − + ρgi. (28) ∂xi Representing the material properties of the fluid, e.g. the dynamic viscosity η, heat conductivity k, the specific heat per unit volume c, and the mass density ρ,as a sum of a constant and a variable part ! η = ηo + !η, k = ko + k, c = co +!c, ρ = ρo + ρ!, (29) then the momentum and energy equations (27) and (18) can be written in analogy to the basic conservation equations formulated for the constant material properties Dv 1 ρ 1 m =−νo∇× ω − ∇ p + g + f , (30) Dt ρo ρo ρo DT S Sm = ao T + + + , (31) Dt co co co with the pseudo-body force term f m and the pseudo-heat source term Sm introduced into the momentum equation (30) and into the energy equation (31), respectively, capturing the variable material property effects, and given by the expressions 4 f m =−∇× (!ηω ) + 2∇ η × ω + 2 ∇ v · ∇ η + ∇ (ηD) − 2D∇ η − ρ!a , (32) 3 or in tensor notation form as ! ! m =− ∂ηωk + ∂η + ∂η ∂vi + 4 ∂ηD − ∂η − ! fi eijk 2eijk ωk 2 2D ρai, (33) ∂xj ∂xj ∂xj ∂xj 3 ∂xi ∂xi WIT Transactions on State of the Art in Science and Engineering, Vol 14, © 2007 WIT Press www.witpress.com, ISSN 1755-8336 (on-line) 38 Domain Decomposition Techniques for Boundary Elements which for plane flow problems reduces to the relation ! ! m =− ∂ηω + ∂η + ∂η ∂vi + 4 ∂ηD − ∂η − ! fi eij 2eij ω 2 2D ρai, (34) ∂xj ∂xj ∂xj ∂xj 3 ∂xi ∂xi while the pseudo-heat source term is given by an expression DT Sm = ∇ (!k ∇ T) −!c , (35) Dt in which the kinematic viscosity is νo = ηo/ρo, the heat diffusivity ao = ko/co, and the inertia acceleration vector is a = Dv/ Dt.
Load more

Recommended publications
  • Che 253M Experiment No. 2 COMPRESSIBLE GAS FLOW
    Rev. 8/15 AD/GW ChE 253M Experiment No. 2 COMPRESSIBLE GAS FLOW The objective of this experiment is to familiarize the student with methods for measurement of compressible gas flow and to study gas flow under subsonic and supersonic flow conditions. The experiment is divided into three distinct parts: (1) Calibration and determination of the critical pressure ratio for a critical flow nozzle under supersonic flow conditions (2) Calculation of the discharge coefficient and Reynolds number for an orifice under subsonic (non- choked) flow conditions and (3) Determination of the orifice constants and mass discharge from a pressurized tank in a dynamic bleed down experiment under (choked) flow conditions. The experimental set up consists of a 100 psig air source branched into two manifolds: the first used for parts (1) and (2) and the second for part (3). The first manifold contains a critical flow nozzle, a NIST-calibrated in-line digital mass flow meter, and an orifice meter, all connected in series with copper piping. The second manifold contains a strain-gauge pressure transducer and a stainless steel tank, which can be pressurized and subsequently bled via a number of attached orifices. A number of NIST-calibrated digital hand held manometers are also used for measuring pressure in all 3 parts of this experiment. Assorted pressure regulators, manual valves, and pressure gauges are present on both manifolds and you are expected to familiarize yourself with the process flow, and know how to operate them to carry out the experiment. A process flow diagram plus handouts outlining the theory of operation of these devices are attached.
    [Show full text]
  • Laws of Similarity in Fluid Mechanics 21
    Laws of similarity in fluid mechanics B. Weigand1 & V. Simon2 1Institut für Thermodynamik der Luft- und Raumfahrt (ITLR), Universität Stuttgart, Germany. 2Isringhausen GmbH & Co KG, Lemgo, Germany. Abstract All processes, in nature as well as in technical systems, can be described by fundamental equations—the conservation equations. These equations can be derived using conservation princi- ples and have to be solved for the situation under consideration. This can be done without explicitly investigating the dimensions of the quantities involved. However, an important consideration in all equations used in fluid mechanics and thermodynamics is dimensional homogeneity. One can use the idea of dimensional consistency in order to group variables together into dimensionless parameters which are less numerous than the original variables. This method is known as dimen- sional analysis. This paper starts with a discussion on dimensions and about the pi theorem of Buckingham. This theorem relates the number of quantities with dimensions to the number of dimensionless groups needed to describe a situation. After establishing this basic relationship between quantities with dimensions and dimensionless groups, the conservation equations for processes in fluid mechanics (Cauchy and Navier–Stokes equations, continuity equation, energy equation) are explained. By non-dimensionalizing these equations, certain dimensionless groups appear (e.g. Reynolds number, Froude number, Grashof number, Weber number, Prandtl number). The physical significance and importance of these groups are explained and the simplifications of the underlying equations for large or small dimensionless parameters are described. Finally, some examples for selected processes in nature and engineering are given to illustrate the method. 1 Introduction If we compare a small leaf with a large one, or a child with its parents, we have the feeling that a ‘similarity’ of some sort exists.
    [Show full text]
  • Aerodynamics Material - Taylor & Francis
    CopyrightAerodynamics material - Taylor & Francis ______________________________________________________________________ 257 Aerodynamics Symbol List Symbol Definition Units a speed of sound ⁄ a speed of sound at sea level ⁄ A area aspect ratio ‐‐‐‐‐‐‐‐ b wing span c chord length c Copyrightmean aerodynamic material chord- Taylor & Francis specific heat at constant pressure of air · root chord tip chord specific heat at constant volume of air · / quarter chord total drag coefficient ‐‐‐‐‐‐‐‐ , induced drag coefficient ‐‐‐‐‐‐‐‐ , parasite drag coefficient ‐‐‐‐‐‐‐‐ , wave drag coefficient ‐‐‐‐‐‐‐‐ local skin friction coefficient ‐‐‐‐‐‐‐‐ lift coefficient ‐‐‐‐‐‐‐‐ , compressible lift coefficient ‐‐‐‐‐‐‐‐ compressible moment ‐‐‐‐‐‐‐‐ , coefficient , pitching moment coefficient ‐‐‐‐‐‐‐‐ , rolling moment coefficient ‐‐‐‐‐‐‐‐ , yawing moment coefficient ‐‐‐‐‐‐‐‐ ______________________________________________________________________ 258 Aerodynamics Aerodynamics Symbol List (cont.) Symbol Definition Units pressure coefficient ‐‐‐‐‐‐‐‐ compressible pressure ‐‐‐‐‐‐‐‐ , coefficient , critical pressure coefficient ‐‐‐‐‐‐‐‐ , supersonic pressure coefficient ‐‐‐‐‐‐‐‐ D total drag induced drag Copyright material - Taylor & Francis parasite drag e span efficiency factor ‐‐‐‐‐‐‐‐ L lift pitching moment · rolling moment · yawing moment · M mach number ‐‐‐‐‐‐‐‐ critical mach number ‐‐‐‐‐‐‐‐ free stream mach number ‐‐‐‐‐‐‐‐ P static pressure ⁄ total pressure ⁄ free stream pressure ⁄ q dynamic pressure ⁄ R
    [Show full text]
  • A Semi-Hydrostatic Theory of Gravity-Dominated Compressible Flow
    Generated using version 3.2 of the official AMS LATEX template 1 A semi-hydrostatic theory of gravity-dominated compressible flow 2 Thomas Dubos ∗ IPSL-Laboratoire de Météorologie Dynamique, Ecole Polytechnique, Palaiseau, France Fabrice voitus CNRM-Groupe d’étude de l’Atmosphere Metéorologique, Météo-France, Toulouse, France ∗Corresponding author address: Thomas Dubos, LMD, École Polytechnique, 91128 Palaiseau, France. E-mail: [email protected] 1 3 Abstract 4 From Hamilton’s least action principle, compressible equations of motion with density diag- 5 nosed from potential temperature through hydrostatic balance are derived. Slaving density 6 to potential temperature suppresses the degrees of freedom supporting the propagation of 7 acoustic waves and results in a sound-proof system. The linear normal modes and dispersion 8 relationship for an isothermal state of rest on f- and β- planes are accurate from hydrostatic 9 to non-hydrostatic scales, except for deep internal gravity waves. Especially the Lamb wave 10 and long Rossby waves are not distorted, unlike with anelastic or pseudo-incompressible 11 systems. 12 Compared to similar equations derived by Arakawa and Konor (2009), the semi-hydrostatic 13 system derived here possesses an additional term in the horizontal momentum budget. This 14 term is an apparent force resulting from the vertical coordinate not being the actual height 15 of an air parcel, but its hydrostatic height, i.e. the hypothetical height it would have after 16 the atmospheric column it belongs to has reached hydrostatic balance through adiabatic 17 vertical displacements of air parcels. The Lagrange multiplier λ introduced in Hamilton’s 18 principle to slave density to potential temperature is identified as the non-hydrostatic ver- 19 tical displacement, i.e.
    [Show full text]
  • Modeling Conservation of Mass
    Modeling Conservation of Mass How is mass conserved (protected from loss)? Imagine an evening campfire. As the wood burns, you notice that the logs have become a small pile of ashes. What happened? Was the wood destroyed by the fire? A scientific principle called the law of conservation of mass states that matter is neither created nor destroyed. So, what happened to the wood? Think back. Did you observe smoke rising from the fire? When wood burns, atoms in the wood combine with oxygen atoms in the air in a chemical reaction called combustion. The products of this burning reaction are ashes as well as the carbon dioxide and water vapor in smoke. The gases escape into the air. We also know from the law of conservation of mass that the mass of the reactants must equal the mass of all the products. How does that work with the campfire? mass – a measure of how much matter is present in a substance law of conservation of mass – states that the mass of all reactants must equal the mass of all products and that matter is neither created nor destroyed If you could measure the mass of the wood and oxygen before you started the fire and then measure the mass of the smoke and ashes after it burned, what would you find? The total mass of matter after the fire would be the same as the total mass of matter before the fire. Therefore, matter was neither created nor destroyed in the campfire; it just changed form. The same atoms that made up the materials before the reaction were simply rearranged to form the materials left after the reaction.
    [Show full text]
  • Derivation of Fluid Flow Equations
    TPG4150 Reservoir Recovery Techniques 2017 1 Fluid Flow Equations DERIVATION OF FLUID FLOW EQUATIONS Review of basic steps Generally speaking, flow equations for flow in porous materials are based on a set of mass, momentum and energy conservation equations, and constitutive equations for the fluids and the porous material involved. For simplicity, we will in the following assume isothermal conditions, so that we not have to involve an energy conservation equation. However, in cases of changing reservoir temperature, such as in the case of cold water injection into a warmer reservoir, this may be of importance. Below, equations are initially described for single phase flow in linear, one- dimensional, horizontal systems, but are later on extended to multi-phase flow in two and three dimensions, and to other coordinate systems. Conservation of mass Consider the following one dimensional rod of porous material: Mass conservation may be formulated across a control element of the slab, with one fluid of density ρ is flowing through it at a velocity u: u ρ Δx The mass balance for the control element is then written as: ⎧Mass into the⎫ ⎧Mass out of the ⎫ ⎧ Rate of change of mass⎫ ⎨ ⎬ − ⎨ ⎬ = ⎨ ⎬ , ⎩element at x ⎭ ⎩element at x + Δx⎭ ⎩ inside the element ⎭ or ∂ {uρA} − {uρA} = {φAΔxρ}. x x+ Δx ∂t Dividing by Δx, and taking the limit as Δx approaches zero, we get the conservation of mass, or continuity equation: ∂ ∂ − (Aρu) = (Aφρ). ∂x ∂t For constant cross sectional area, the continuity equation simplifies to: ∂ ∂ − (ρu) = (φρ) . ∂x ∂t Next, we need to replace the velocity term by an equation relating it to pressure gradient and fluid and rock properties, and the density and porosity terms by appropriate pressure dependent functions.
    [Show full text]
  • Law of Conversation of Energy
    Law of Conservation of Mass: "In any kind of physical or chemical process, mass is neither created nor destroyed - the mass before the process equals the mass after the process." - the total mass of the system does not change, the total mass of the products of a chemical reaction is always the same as the total mass of the original materials. "Physics for scientists and engineers," 4th edition, Vol.1, Raymond A. Serway, Saunders College Publishing, 1996. Ex. 1) When wood burns, mass seems to disappear because some of the products of reaction are gases; if the mass of the original wood is added to the mass of the oxygen that combined with it and if the mass of the resulting ash is added to the mass o the gaseous products, the two sums will turn out exactly equal. 2) Iron increases in weight on rusting because it combines with gases from the air, and the increase in weight is exactly equal to the weight of gas consumed. Out of thousands of reactions that have been tested with accurate chemical balances, no deviation from the law has ever been found. Law of Conversation of Energy: The total energy of a closed system is constant. Matter is neither created nor destroyed – total mass of reactants equals total mass of products You can calculate the change of temp by simply understanding that energy and the mass is conserved - it means that we added the two heat quantities together we can calculate the change of temperature by using the law or measure change of temp and show the conservation of energy E1 + E2 = E3 -> E(universe) = E(System) + E(Surroundings) M1 + M2 = M3 Is T1 + T2 = unknown (No, no law of conservation of temperature, so we have to use the concept of conservation of energy) Total amount of thermal energy in beaker of water in absolute terms as opposed to differential terms (reference point is 0 degrees Kelvin) Knowns: M1, M2, T1, T2 (Kelvin) When add the two together, want to know what T3 and M3 are going to be.
    [Show full text]
  • Chapter 5 Dimensional Analysis and Similarity
    Chapter 5 Dimensional Analysis and Similarity Motivation. In this chapter we discuss the planning, presentation, and interpretation of experimental data. We shall try to convince you that such data are best presented in dimensionless form. Experiments which might result in tables of output, or even mul- tiple volumes of tables, might be reduced to a single set of curves—or even a single curve—when suitably nondimensionalized. The technique for doing this is dimensional analysis. Chapter 3 presented gross control-volume balances of mass, momentum, and en- ergy which led to estimates of global parameters: mass flow, force, torque, total heat transfer. Chapter 4 presented infinitesimal balances which led to the basic partial dif- ferential equations of fluid flow and some particular solutions. These two chapters cov- ered analytical techniques, which are limited to fairly simple geometries and well- defined boundary conditions. Probably one-third of fluid-flow problems can be attacked in this analytical or theoretical manner. The other two-thirds of all fluid problems are too complex, both geometrically and physically, to be solved analytically. They must be tested by experiment. Their behav- ior is reported as experimental data. Such data are much more useful if they are ex- pressed in compact, economic form. Graphs are especially useful, since tabulated data cannot be absorbed, nor can the trends and rates of change be observed, by most en- gineering eyes. These are the motivations for dimensional analysis. The technique is traditional in fluid mechanics and is useful in all engineering and physical sciences, with notable uses also seen in the biological and social sciences.
    [Show full text]
  • Sensitivity Analysis of Non-Linear Steep Waves Using VOF Method
    Tenth International Conference on ICCFD10-268 Computational Fluid Dynamics (ICCFD10), Barcelona, Spain, July 9-13, 2018 Sensitivity Analysis of Non-linear Steep Waves using VOF Method A. Khaware*, V. Gupta*, K. Srikanth *, and P. Sharkey ** Corresponding author: [email protected] * ANSYS Software Pvt Ltd, Pune, India. ** ANSYS UK Ltd, Milton Park, UK Abstract: The analysis and prediction of non-linear waves is a crucial part of ocean hydrodynamics. Sea waves are typically non-linear in nature, and whilst models exist to predict their behavior, limits exist in their applicability. In practice, as the waves become increasingly steeper, they approach a point beyond which the wave integrity cannot be maintained, and they 'break'. Understanding the limits of available models as waves approach these break conditions can significantly help to improve the accuracy of their potential impact in the field. Moreover, inaccurate modeling of wave kinematics can result in erroneous hydrodynamic forces being predicted. This paper investigates the sensitivity of non-linear wave modeling from both an analytical and a numerical perspective. Using a Volume of Fluid (VOF) method, coupled with the Open Channel Flow module in ANSYS Fluent, sensitivity studies are performed for a variety of non-linear wave scenarios with high steepness and high relative height. These scenarios are intended to mimic the near-break conditions of the wave. 5th order solitary wave models are applied to shallow wave scenarios with high relative heights, and 5th order Stokes wave models are applied to short gravity waves with high wave steepness. Stokes waves are further applied in the shallow regime at high wave steepness to examine the wave sensitivity under extreme conditions.
    [Show full text]
  • Key Concepts: Conservation of Mass, Momentum, Energy Fluid: a Material
    Key concepts: Conservation of mass, momentum, energy Fluid: a material that deforms continuously and permanently under the application of a shearing stress, no matter how small. Fluids are either gases or liquids. (Under very specialized conditions, a phase of intermediate properties can be stable, but we won’t consider that possibility.) In liquids, the molecules are relatively closely spaced, allowing the magnitude of their (attractive, electrically-based) interaction energy to be of the same magnitude as their kinetic energy. As a result, they exist as a loose collection of clusters. In gases, the molecules are much more widely separated, so the kinetic energy (at a given temperature, identical to that in the liquid) is far greater than the interaction energy (much less than in the liquid), and molecule do not form clusters. In a liquid, the molecules themselves typically occupy a few percent of the total space available; in a gas, they occupy a few thousandths of a percent. Nevertheless, for our purposes, all fluids are considered to be continua (no voids or holes).The absence of significant intermolecular attraction allows gases to fill whatever volume is available to them, whereas the presence of such attraction in liquids prevents them from doing so. The attractive forces in liquid water are unusually strong, compared to other liquids. Properties of Fluids: m Density is mass/volume: ρ = . The density of liquid water is V 3 o −3 3 ~1.0 kg/m ; that of air at 20 C is ~1.2x10 kg/m . mg W Specific weight is weight/volume: γ = = = ρg V V C:\Adata\CLASNOTE\342\Class Notes\Key concepts_Topic 1.doc 1 γ Specific gravity is density normalized to the density of water: sg..= i γ w V 1 Specific volume is volume/mass: V = = m ρ Bulk modulus or modulus of elasticity is the pressure change per dp dp fractional change in volume or density: E =− = .
    [Show full text]
  • Waves and Structures
    WAVES AND STRUCTURES By Dr M C Deo Professor of Civil Engineering Indian Institute of Technology Bombay Powai, Mumbai 400 076 Contact: [email protected]; (+91) 22 2572 2377 (Please refer as follows, if you use any part of this book: Deo M C (2013): Waves and Structures, http://www.civil.iitb.ac.in/~mcdeo/waves.html) (Suggestions to improve/modify contents are welcome) 1 Content Chapter 1: Introduction 4 Chapter 2: Wave Theories 18 Chapter 3: Random Waves 47 Chapter 4: Wave Propagation 80 Chapter 5: Numerical Modeling of Waves 110 Chapter 6: Design Water Depth 115 Chapter 7: Wave Forces on Shore-Based Structures 132 Chapter 8: Wave Force On Small Diameter Members 150 Chapter 9: Maximum Wave Force on the Entire Structure 173 Chapter 10: Wave Forces on Large Diameter Members 187 Chapter 11: Spectral and Statistical Analysis of Wave Forces 209 Chapter 12: Wave Run Up 221 Chapter 13: Pipeline Hydrodynamics 234 Chapter 14: Statics of Floating Bodies 241 Chapter 15: Vibrations 268 Chapter 16: Motions of Freely Floating Bodies 283 Chapter 17: Motion Response of Compliant Structures 315 2 Notations 338 References 342 3 CHAPTER 1 INTRODUCTION 1.1 Introduction The knowledge of magnitude and behavior of ocean waves at site is an essential prerequisite for almost all activities in the ocean including planning, design, construction and operation related to harbor, coastal and structures. The waves of major concern to a harbor engineer are generated by the action of wind. The wind creates a disturbance in the sea which is restored to its calm equilibrium position by the action of gravity and hence resulting waves are called wind generated gravity waves.
    [Show full text]
  • Continuity Equation in Pressure Coordinates
    Continuity Equation in Pressure Coordinates Here we will derive the continuity equation from the principle that mass is conserved for a parcel following the fluid motion (i.e., there is no flow across the boundaries of the parcel). This implies that δxδyδp δM = ρ δV = ρ δxδyδz = − g is conserved following the fluid motion: 1 d(δM ) = 0 δM dt 1 d()δM = 0 δM dt g d ⎛ δxδyδp ⎞ ⎜ ⎟ = 0 δxδyδp dt ⎝ g ⎠ 1 ⎛ d(δp) d(δy) d(δx)⎞ ⎜δxδy +δxδp +δyδp ⎟ = 0 δxδyδp ⎝ dt dt dt ⎠ 1 ⎛ dp ⎞ 1 ⎛ dy ⎞ 1 ⎛ dx ⎞ δ ⎜ ⎟ + δ ⎜ ⎟ + δ ⎜ ⎟ = 0 δp ⎝ dt ⎠ δy ⎝ dt ⎠ δx ⎝ dt ⎠ Taking the limit as δx, δy, δp → 0, ∂u ∂v ∂ω Continuity equation + + = 0 in pressure ∂x ∂y ∂p coordinates 1 Determining Vertical Velocities • Typical large-scale vertical motions in the atmosphere are of the order of 0. 01-01m/s0.1 m/s. • Such motions are very difficult, if not impossible, to measure directly. Typical observational errors for wind measurements are ~1 m/s. • Quantitative estimates of vertical velocity must be inferred from quantities that can be directly measured with sufficient accuracy. Vertical Velocity in P-Coordinates The equivalent of the vertical velocity in p-coordinates is: dp ∂p r ∂p ω = = +V ⋅∇p + w dt ∂t ∂z Based on a scaling of the three terms on the r.h.s., the last term is at least an order of magnitude larger than the other two. Making the hydrostatic approximation yields ∂p ω ≈ w = −ρgw ∂z Typical large-scale values: for w, 0.01 m/s = 1 cm/s for ω, 0.1 Pa/s = 1 μbar/s 2 The Kinematic Method By integrating the continuity equation in (x,y,p) coordinates, ω can be obtained from the mean divergence in a layer: ⎛ ∂u ∂v ⎞ ∂ω ⎜ + ⎟ + = 0 continuity equation in (x,y,p) coordinates ⎝ ∂x ∂y ⎠ p ∂p p2 p2 ⎛ ∂u ∂v ⎞ ∂ω = − ⎜ + ⎟ ∂p rearrange and integrate over the layer ∫p ∫ ⎜ ⎟ 1 ∂x ∂y p1⎝ ⎠ p ⎛ ∂u ∂v ⎞ ω(p )−ω(p ) = (p − p )⎜ + ⎟ overbar denotes pressure- 2 1 1 2 ⎜ ⎟ weighted vertical average ⎝ ∂x ∂y ⎠ p To determine vertical motion at a pressure level p2, assume that p1 = surface pressure and there is no vertical motion at the surface.
    [Show full text]