A Diagnostic Model for Initial Winds in Primitive

Total Page:16

File Type:pdf, Size:1020Kb

A Diagnostic Model for Initial Winds in Primitive 1 1 A DIAGNOSTIC MODEL FOR INITIAL WINDS IN PRIMITIVE EQUATIONS FORECASTS by Richard Asselin Department of Meteorology, Ph. D. ABSTRACT The set of primitive meteorologicai equations reduces to a single second order partial differential cquation when the accel- eration is replaced by the derivative of the geostrophic wind in the equation.of motion. From this equation a diagnostic three-dimensional wind field can be computed when the pressure distribution is known. Calculations for a barotropic atmosphere are performed fir st; the rotational part of the wind is found to be similar to that from the classical balance equation and the divergent part agrees well with synoptic experience. A forecast thus initialized contaills Lese; gravi- tational activity than one initialized with the non-divergent wind of the balance equation. With a five-level baroclinic atmosphere. the effects of mountains, surface drag, and release of latent heat are studi~d separately; the results arc consistent wlth classical concepts and other interesting features are revealed. T!;e diagnostic surfa:::e pres­ sure tendency is fo~nd to be in very good agreement with the observed pressure changes. A DIAGNOSTIC MODEL FOR INITIAL WINDS IN PRIMITIVE EQUATiONS FORECASTS by RI CHARD ASSE LIN . A the sis submitted to the Faculty of Graduate Studies and Re search in partial fulfillrnent of the requirernents for the degree of Doctor of Philosophy. Departrnent of Meteorology McGill Univer sity Montreal, Canada July. 1970 1971 A DL~GNOSTIC MODEL FOR INITIAL WINDS IN PRIMITIVE EQUATIONS FORECASTS by RI CHARD ASSE LIN . A thesis submitted ta the Faculty of Graduate Studies and Research in partial fulfillment of the requirements. for the degree of Doctor of Philosophy. Department of Meteorology McGill Univer sity Montreal, Canada July, 1970 @ l-J.cha l'Cl L.s :::;e1in 1971 ACKNOWLEDGEMENTS The author is indebted to his research supervisor, Professor A. J. Robert, for many useful discussions and for his encouragement throughout this work, as weU as to Professor B. W. Boville, who provided helpful advice and editorial assistance. l am grateful to Dr. M. Kwizak for generating my initial incentive in these studies and for his continuing suppurt, particularly in arranging for me to use the computer and other facili­ ties of the Central Analysis Office (CAO). The programs of the barotropic forecast model as well as a diagram were borrowed from his work. Miss S. Boville is responsible for the programming of the experiment with idealized data. Mr. W. S. Creswick did the graphical computations and helped to evaluate the re~u1ts. Mr. P. Sarrazin drafted the figures and Misses T. Savoie and J. Sutherland typed the manuscript. Finally, l wish to acknowledge the financial support of the Canadian Meteorological Servi.ce and of the National Research Council of Canada. ii TABLE OF CONTENTS Acknowledgements ii List of Plates; Tables and Figures vi List of Symbols ix Résumé xi Abstract xii 1. Introduction 1 2. Review of Diagnostic Methods 6 2. 1. Balance -W- Equations 6 2. 2. Geostrophic Diagnostic (G ) 8 2 2. 3. Miyakoda-Moyer Diagnostic 9 2.4. Nitta-Hovermale Method 9 3. Theory of Initialization 11 3. 1. The Simple st Model 11 3. 2. Numerical Compatibility 15 4. Historical R~.view of the G Geos.trophic Method 16 2 Part A: BAROTROPIC EXPERIMENTS· 19 5. Basic Equations and Theory 19 6. Numerical Solution of the Diagnostic Model 25 6. 1. Relaxation 25 6. 2. Ellipticity 26 7. Hig~er Order G Approximation 30 2 iii 8. Numerical Procedures for Initialization of Forecasts 31 8. 1. Compatibility Procedure 31 8. 2. Imposing New Boundary Conditions 34 9. Discussion of Results 37 9. 1. Divergence and Tendency 37 9. 2. Vorticity 43 9.3. Height Changes 46 9.4. Forecasts 50 Part B: BAROCLINIC EXPERIMENTS 54 10. Basic Primitive Equations 54 10.1. Precipitation Model 56 Il. Diagnostic Equations 59 12. Solution of Baroclinic System 64 12. 1. Discretization 64 12.2. Boundary Condition~, Vertical Structure. 65 12.3. Relaxation 67 13. Data Preparation 69 13.1. Interpolation to (j-SurfaceEV 69 ~----- .13. 2. Smoothing 70 14. Ellipticity 72 15. _ Discussion of Results 75 15.1. Effects of Orography and Drag 75 15.2. Effect of Pressure Systems, Ellipticity 85 iv 15.3. Effect of Latent Heat 96 15.4. Combined Effect from All Sources 102 16. Conclusions 109 References III Appendix A: Finitc Differences 115 Appendix B: Special Operations for Equatorial Latitudes 119 Appendix C: 123 v LIST OF PLATES. TABLES. AND FIGURES pagè Plate 1. The baroclinic diagnostic equations in full 60 Table 1. Input data for each baroclinic experiment 76 Table 2. Statistics from each baroclinic experiment 95 Figure 1. The 51. x 55 point grid 5. 2. Two graphs of the tangential wind speed in a station­ ary circular system 23 3. 500 mb analysis for OOz 21 Feb 69 and diagnosed tendency 38 4. Summaryof 25 day integration by Kwizak (1970, Energy. vorticity, divergence • 39 5. Time variation of total available potential energy in barotropic forecasts with bala~ce and G initia1i­ 2 zations 41 6. Time variation of squared divergence in barotropic forecasts with balance and G initializations 42 2 7. Absolute vorticity from balance equation 44 8. Absolute vorticity from G geostrophic initialization 45 2 9. Total height changes ma.de by bal1:\.nce ilütialization 47 10. Total height changes ~ade by G initialization 48 D 2 vi Figure 11. 72-hour forecast from balance initialization 51 12. 72-hour forecast from G initialization 52 2 13. 72-hour verifying analysis, OOZ 24 Feb 1969 53 14. The discrete structure of the baroclinic mode1 66 15. a) response of short wave filter b) response of 10w latitude smoother as functions of wave1ength 71 16. Idealized data used for experiments 7 to 12 78 17. Orientation error pattern in baroclinic model 79 18. Basic results of experiments with idealized mountain 81 19. Resu1ts of Cressman' s (1960) experiments with orography 84 20. Diagnosed surface tendency produced by rea1 oro­ graphy a10ne 86 21. Diagnosed surface tendency produced by orography and 'drag combined 87 22. 300 mb ana1ysis for 122 10 Feb 1970 88 23. 500 mb ana1ysis for 12Z 10 Feb 1970 89 24. 1000 mb analysis for 122 10 Feb. 1970 90 25. 300 mb height change~ resulting from preliminary ellipticization 92 26. 500 mb height changes resulting from preliminary ellipticizati?n 93 vii Figure 27. Results from basic baroc1inic experiment: vertical m.otion, vorticity and tendency for a small volume of space 94 28. Surface chart for 062 10 Feb 1970 97 29. Surface chart for 122 10 Feb 1970 98 30. Surface chart for 182 10 Feb 1970 99 31. Diagnostic vertical motion accompanying observed pressure distribution assuming saturated troposphere 100 32. Diagnostic vertical motion accompanying observed pressure distribution alone 101 33. Diagnostic vertical inotion accompanying observed pressure and moisture distributions with mountain and friction effects also considered 103 34. Vertical motion from the w-equaÜon inc1uding mountain and friction effects . 104 35. Surface pressure tendency from most complete diagnostic model 105 36. Observed 12-hour surface pressure change 107 37. Exact and modified functions of the Coriolis 122 parameter used in baroclinic experiments viii LIST OF SYMBOLS a element of matrix A A' 2 x 2 or 3 x 3 matr ix b element of matrix A B 3 x 3 matrix c element of matrix A C condensation rate - cp specifie heat at const. pre 55. CD surfa.ce drag coefficient d element of matrix A, D horizontal wind divergence distance between grid points F unspecHiedfunction e water vapour pressure f friction f orc e f = 2.Qsin (l, Coriolis parame G time averaged geopotential, g accelei-ation of gravit y, Geostrophic, Gravit y subscript: geostrophic H mean depth of nuid, h height of free surface heating rate unit imaginary numbcr, K kinetic energy, Kelvin index or unit vector for x L latent heat of vapourization j index or unit vector for y pep* ) pressure (at the grcu!'\d), k index or unit vertical vcctor polynomial m map factor Q = f + f absolute vorticity p element of matrix A R radhls of curvature, residual, q mixing ratio Rossby, Gas constant r clement of matrix. A 5 water vapvl:.r SOUl-ce or sink, t time subscript: stationary u velocity in x direction T temperature v velocity in y direction, u velocity in x directio!l subscript: virtual v velo city in y direction w weight W' vector wind Ct x, y grid cartesian coord. W vertical velocity z vertical direction,- height ·of constant pressure surface ix Q wave number, relaxation factor, arbitrary factor f3 arbitrary factor y arbitrary factor ô arbitrary factor, Kronecker delta A difference e: small number used as limit, arbitrary factor ~ relative vorticity about vertical axis 6 latitude k = R/cp 1\ eigen value V iteration inde.x, smoothing coefficient 1T product sign cr _ p Ip*, vertical coordinate é1: E dtr/dt, vertical velocity ~ summation sign 1> = gh or gz, geopotential ~ velocity potential \fi tendency, stream function (,J :.: dP Idt, vertical velocity Il.. earth' s angular velocity DIFFERENTIAL OPERA TORS A. total differential a partial differ ential V = i dlèJ x + j dIa y, gradient 2 ",2= dl Iè. xl + a Id yl, Laplacian J(Fl,Fl) = èiF1/dx· dFl/O y - dFl/ë)y· àFZldX, Jacobian x RESUME ~ , L'ensemble des lois gouvernant la meteorologie dyna- mique se condense en une seule équation differentielle du second ordre . " , .. " si on substitue la derivee du vent geostrophique a l'acceleration dans la loi du mouvement. On peut ensuite se servir de cette équation pour d~terminer le champ des trois composantes du vent à partir de la distribution de la.
Recommended publications
  • Lecture 5: Saturn, Neptune, … A. P. Ingersoll [email protected] 1
    Lecture 5: Saturn, Neptune, … A. P. Ingersoll [email protected] 1. A lot a new material from the Cassini spacecraft, which has been in orbit around Saturn for four years. This talk will have a lot of pictures, not so much models. 2. Uranus on the left, Neptune on the right. Uranus spins on its side. The obliquity is 98 degrees, which means the poles receive more sunlight than the equator. It also means the sun is almost overhead at the pole during summer solstice. Despite these extreme changes in the distribution of sunlight, Uranus is a banded planet. The rotation dominates the winds and cloud structure. 3. Saturn is covered by a layer of clouds and haze, so it is hard to see the features. Storms are less frequent on Saturn than on Jupiter; it is not just that clouds and haze makes the storms less visible. Saturn is a less active planet. Nevertheless the winds are stronger than the winds of Jupiter. 4. The thickness of the atmosphere is inversely proportional to gravity. The zero of altitude is where the pressure is 100 mbar. Uranus and Neptune are so cold that CH forms clouds. The clouds of NH3 , NH4SH, and H2O form at much deeper levels and have not been detected. 5. Surprising fact: The winds increase as you move outward in the solar system. Why? My theory is that power/area is smaller, small‐scale turbulence is weaker, dissipation is less, and the winds are stronger. Jupiter looks more turbulent, and it has the smallest winds.
    [Show full text]
  • ATMOSPHERIC and OCEANIC FLUID DYNAMICS Fundamentals and Large-Scale Circulation
    ATMOSPHERIC AND OCEANIC FLUID DYNAMICS Fundamentals and Large-Scale Circulation Geoffrey K. Vallis Contents Preface xi Part I FUNDAMENTALS OF GEOPHYSICAL FLUID DYNAMICS 1 1 Equations of Motion 3 1.1 Time Derivatives for Fluids 3 1.2 The Mass Continuity Equation 7 1.3 The Momentum Equation 11 1.4 The Equation of State 14 1.5 The Thermodynamic Equation 16 1.6 Sound Waves 29 1.7 Compressible and Incompressible Flow 31 1.8 * More Thermodynamics of Liquids 33 1.9 The Energy Budget 39 1.10 An Introduction to Non-Dimensionalization and Scaling 43 2 Effects of Rotation and Stratification 51 2.1 Equations in a Rotating Frame 51 2.2 Equations of Motion in Spherical Coordinates 55 2.3 Cartesian Approximations: The Tangent Plane 66 2.4 The Boussinesq Approximation 68 2.5 The Anelastic Approximation 74 2.6 Changing Vertical Coordinate 78 2.7 Hydrostatic Balance 80 2.8 Geostrophic and Thermal Wind Balance 85 2.9 Static Instability and the Parcel Method 92 2.10 Gravity Waves 98 v vi Contents 2.11 * Acoustic-Gravity Waves in an Ideal Gas 100 2.12 The Ekman Layer 104 3 Shallow Water Systems and Isentropic Coordinates 123 3.1 Dynamics of a Single, Shallow Layer 123 3.2 Reduced Gravity Equations 129 3.3 Multi-Layer Shallow Water Equations 131 3.4 Geostrophic Balance and Thermal wind 134 3.5 Form Drag 135 3.6 Conservation Properties of Shallow Water Systems 136 3.7 Shallow Water Waves 140 3.8 Geostrophic Adjustment 144 3.9 Isentropic Coordinates 152 3.10 Available Potential Energy 155 4 Vorticity and Potential Vorticity 165 4.1 Vorticity and Circulation 165
    [Show full text]
  • Chapter 7 Balanced Flow
    Chapter 7 Balanced flow In Chapter 6 we derived the equations that govern the evolution of the at- mosphere and ocean, setting our discussion on a sound theoretical footing. However, these equations describe myriad phenomena, many of which are not central to our discussion of the large-scale circulation of the atmosphere and ocean. In this chapter, therefore, we focus on a subset of possible motions known as ‘balanced flows’ which are relevant to the general circulation. We have already seen that large-scale flow in the atmosphere and ocean is hydrostatically balanced in the vertical in the sense that gravitational and pressure gradient forces balance one another, rather than inducing accelera- tions. It turns out that the atmosphere and ocean are also close to balance in the horizontal, in the sense that Coriolis forces are balanced by horizon- tal pressure gradients in what is known as ‘geostrophic motion’ – from the Greek: ‘geo’ for ‘earth’, ‘strophe’ for ‘turning’. In this Chapter we describe how the rather peculiar and counter-intuitive properties of the geostrophic motion of a homogeneous fluid are encapsulated in the ‘Taylor-Proudman theorem’ which expresses in mathematical form the ‘stiffness’ imparted to a fluid by rotation. This stiffness property will be repeatedly applied in later chapters to come to some understanding of the large-scale circulation of the atmosphere and ocean. We go on to discuss how the Taylor-Proudman theo- rem is modified in a fluid in which the density is not homogeneous but varies from place to place, deriving the ‘thermal wind equation’. Finally we dis- cuss so-called ‘ageostrophic flow’ motion, which is not in geostrophic balance but is modified by friction in regions where the atmosphere and ocean rubs against solid boundaries or at the atmosphere-ocean interface.
    [Show full text]
  • Dynamics of Jupiter's Atmosphere
    6 Dynamics of Jupiter's Atmosphere Andrew P. Ingersoll California Institute of Technology Timothy E. Dowling University of Louisville P eter J. Gierasch Cornell University GlennS. Orton J et P ropulsion Laboratory, California Institute of Technology Peter L. Read Oxford. University Agustin Sanchez-Lavega Universidad del P ais Vasco, Spain Adam P. Showman University of A rizona Amy A. Simon-Miller NASA Goddard. Space Flight Center Ashwin R . V asavada J et Propulsion Laboratory, California Institute of Technology 6.1 INTRODUCTION occurs on both planets. On Earth, electrical charge separa­ tion is associated with falling ice and rain. On Jupiter, t he Giant planet atmospheres provided many of the surprises separation mechanism is still to be determined. and remarkable discoveries of planetary exploration during The winds of Jupiter are only 1/ 3 as strong as t hose t he past few decades. Studying Jupiter's atmosphere and of aturn and Neptune, and yet the other giant planets comparing it with Earth's gives us critical insight and a have less sunlight and less internal heat than Jupiter. Earth broad understanding of how atmospheres work that could probably has the weakest winds of any planet, although its not be obtained by studying Earth alone. absorbed solar power per unit area is largest. All the gi­ ant planets are banded. Even Uranus, whose rotation axis Jupiter has half a dozen eastward jet streams in each is tipped 98° relative to its orbit axis. exhibits banded hemisphere. On average, Earth has only one in each hemi­ cloud patterns and east- west (zonal) jets.
    [Show full text]
  • ATM 500: Atmospheric Dynamics End-Of-Semester Review, December 6 2017 Here, in Bullet Points, Are the Key Concepts from the Second Half of the Course
    ATM 500: Atmospheric Dynamics End-of-semester review, December 6 2017 Here, in bullet points, are the key concepts from the second half of the course. They are not meant to be self-contained notes! Refer back to your course notes for notation, derivations, and deeper discussion. Static stability The vertical variation of density in the environment determines whether a parcel, disturbed upward from an initial hydrostatic rest state, will accelerate away from its initial position, or oscillate about that position. For an incompressible or Boussinesq fluid, the buoyancy frequency is g dρ~ N 2 = − ρ~dz For a compressible ideal gas, the buoyancy frequency is g dθ~ N 2 = θ~dz (the expressions are different because density is not conserved in a small upward displacement of a compressible fluid, but potential temperature is conserved) Fluid is statically stable if N 2 > 0; otherwise it is statically unstable. Typical value for troposphere: N = 0:01 s−1, with corresponding oscillation period of 10 minutes. Dry adiabatic lapse rate g −1 Γd = = 9:8 K km cp The rate at which temperature decreases during adiabatic ascent (for which Dθ Dt = 0) \Dry" here means that there is no latent heating from condensation of water vapor. Static stability criteria for a dry atmosphere Environment is stable to dry con- vection if dθ~ dT~ > 0 or − < Γ dz dz d and otherwise unstable. 1 Buoyancy equation for Boussinesq fluid with background stratification Db0 + N 2w = 0 Dt where we have separated the full buoyancy into a mean vertical part ^b(z) and 0 2 d^b 2 everything else (b ), and N = dz .
    [Show full text]
  • Ch4: Basics of Stratified Fluid ∂
    AOS611Ch.4, Z. Liu, 02/16/2016 1 Ch4: Basics of Stratified Fluid Sec. 4.1: Basic Equations Stratification will introduce new physics. We first derive the equations. Denoting rotation vector as Ω , gravity potential , u u,v, w, i x j y k z . The momentum equations are du 2Ω u p F . (4.1.1) dt The mass equation is d u 0 or u 0 (4.1.2) dt t The equation of state in general is p, T, S... (4.1.3) In the ocean, neglecting salinity, the equation of state is m 1T T0 (4.1.4) 1 where is the coefficient of thermal expansion and is a constant m T P reference density, which can be chosen as the average density. In the atmosphere, using the perfect gas, the equation of state is P (4.1.5) RT where R is the gas constant. The thermodynamic equation describes the internal energy change. For the ocean, which is incompressible, the thermodynamic equation is dT ~ C Q k 2T (4.1.6) p dt ~ where k is the thermal conductivity, Q is the heating rate per unit mass. This can be rewritten as dT J k 2T (4.1.6a) dt ~ Q k where J and k . Cp Cp Copyright 2014, Zhengyu Liu AOS611Ch.4, Z. Liu, 02/16/2016 2 The atmosphere is compressible, so the thermodynamic equation is dT 1 dp J k 2T (4.1.7) dt C p dt p 0 R , Define the potential temperature as T , where C we have p p p0 ln lnT ln lnT ln p const p d dT dp T p T p0 T d dT dp dT dp T p p p d p0 2 J k T (4.1.8) dt p T RT 1 where we have used (4.1.7) and the ideal gas law (4.1.5), such that .
    [Show full text]
  • Understanding Basic Concepts of Numerical Weather Prediction Using Single Level Atmospheric Model
    American Journal of Engineering Research (AJER) 2017 American Journal of Engineering Research (AJER) e-ISSN: 2320-0847 p-ISSN : 2320-0936 Volume-6, Issue-3, pp-57-62 www.ajer.org Research Paper Open Access Understanding Basic Concepts of Numerical Weather Prediction Using Single Level Atmospheric Model Ahmed Z. Abdulmajed, Kais J. Al-Jumaily Department of Atmospheric Sciences, College of Science, Al-Mustansiriyah University, Baghdad, Iraq ABSTRACT: Weather impacts all daily lives from a very hot summer day to a cold winter storm. Predicting these events and other weather phenomena has lead scientist to develop the most accurate way to forecast weather, numerical weather prediction. This system developed over the last 100 years has created more accurate and longer forecast periods, as well as, created a highly complex system of calculations designed to simulate atmospheric motions. This simulation requires the use of large computing systems for the computation but through the use of a simplified computer program one can gain a full understanding of the process of numerical weather prediction and the complexity of its parts. The aim of this research is to use a single level atmospheric model to understand some basic concepts of numerical weather prediction. Reference geopotential height of 5500 meters which represent the height of 500 mb pressure level was used as a single level. The model results showed that initial geopotential has two distinct centres of low and high geopotentials and the high centre is situated close to reference point. The results also indicated that the change in the time step interval of the integration does not affect the solution of shallow water equations and increasing the forecast length would result in the deformation in the geopotential patterns.
    [Show full text]
  • The Problem of the Professional Training of Meteorological Personnel of All Grades in the Less-Developed Countries
    WORLD METEOROLOGICAL ORGANIZATION TECHNICAL NOTE No. 50 THE PROBLEM OF THE PROFESSIONAL TRAINING OF METEOROLOGICAL PERSONNEL OF ALL GRADES IN THE LESS-DEVELOPED COUNTRIES PROFESSOR J. VAN MIEGHEM Consultant on Meteorological Training to the Secretary-General of WMO Reprinted in 1967 PRICE: Sw. fro 6.- I WMO· No. 132.TP. 59 I Secretariat of the World Meteorological Organization - Geneva - Switzerland 1963 NOTE The designations employed and the presentation of the material in this publication do not imply the expression of any opinion whatsoever on the part of the Secretariat of the World Meteorological Organization concerning the legal status of any country or territory or of its authorities, or concerning the delimitation of its frontiers. TABLE OF CONTENTS Foreword • .. .. .. .. .. .. .. .. .. .. .. .. • .. .. .. Summaries (English, French, Russian, Spanish). VII 1. Introduction .. .. .. 1 2. Collaboration with UNESCO and ICAO 1 3· Qualifications and recruitment of meteorological personnel 2 4. Equivalence of curricula and diplomas, standardization of Syllabi. 4 The diversity of personnel in a meteorological service and the need for uniform grading .. .. .. 5 6. Short-course training 7 7· Grading of meteorological personnel into four classes I duties devolving on each class and standards of knowledge required. 7 7.1 Class I meteorological personnel • 8 7.2 Class II meteorological personnel 8 7.3 Class III meteorological personnel 9 7.4 Class IV meteorological personnel 9 7.5 Syllabi of the courses • • • • • • . 9 7.6 Duration of meteorological courses 10 8. Plan for developing professional meteorological training 11 8.1 Fundamentals .. •• " •• • • •• • .. • • •••••• .. • 11 8.2 Instruction and tTaining of Class IV meteorological personnel 12 8.3 Instruction and training of Class III and Class II meteorological personnel.
    [Show full text]
  • Circulation and Vorticity
    Lecture 4: Circulation and Vorticity • Circulation •Bjerknes Circulation Theorem • Vorticity • Potential Vorticity • Conservation of Potential Vorticity ESS227 Prof. Jin-Yi Yu Measurement of Rotation • Circulation and vorticity are the two primary measures of rotation in a fluid. • Circul ati on, whi ch i s a scal ar i nt egral quantit y, i s a macroscopic measure of rotation for a finite area of the fluid. • Vorticity, however, is a vector field that gives a microscopic measure of the rotation at any point in the fluid. ESS227 Prof. Jin-Yi Yu Circulation • The circulation, C, about a closed contour in a fluid is defined as thliihe line integral eval uated dl along th e contour of fh the component of the velocity vector that is locally tangent to the contour. C > 0 Î Counterclockwise C < 0 Î Clockwise ESS227 Prof. Jin-Yi Yu Example • That circulation is a measure of rotation is demonstrated readily by considering a circular ring of fluid of radius R in solid-body rotation at angular velocity Ω about the z axis . • In this case, U = Ω × R, where R is the distance from the axis of rotation to the ring of fluid. Thus the circulation about the ring is given by: • In this case the circulation is just 2π times the angular momentum of the fluid ring about the axis of rotation. Alternatively, note that C/(πR2) = 2Ω so that the circulation divided by the area enclosed by the loop is just twice thlhe angular spee dfifhid of rotation of the ring. • Unlike angular momentum or angular velocity, circulation can be computed without reference to an axis of rotation; it can thus be used to characterize fluid rotation in situations where “angular velocity” is not defined easily.
    [Show full text]
  • Moist Versus Dry Barotropic Instability in a Shallow-Water Model of the Atmosphere with Moist Convection
    1234 JOURNAL OF THE ATMOSPHERIC SCIENCES VOLUME 68 Moist versus Dry Barotropic Instability in a Shallow-Water Model of the Atmosphere with Moist Convection JULIEN LAMBAERTS,GUILLAUME LAPEYRE, AND VLADIMIR ZEITLIN Laboratoire de Me´te´orologie Dynamique/IPSL, Ecole Normale Supe´rieure/CNRS/UPMC, Paris, France (Manuscript received 28 April 2010, in final form 6 January 2011) ABSTRACT Dynamical influence of moist convection upon development of the barotropic instability is studied in the rotating shallow-water model. First, an exhaustive linear ‘‘dry’’ stability analysis of the Bickley jet is per- formed, and the most unstable mode identified in this way is used to initialize simulations to compare the development and the saturation of the instability in dry and moist configurations. High-resolution numerical simulations with a well-balanced finite-volume scheme reveal substantial qualitative and quantitative dif- ferences in the evolution of dry and moist-convective instabilities. The moist effects affect both balanced and unbalanced components of the flow. The most important differences between dry and moist evolution are 1) the enhanced efficiency of the moist-convective instability, which manifests itself by the increase of the growth rate at the onset of precipitation, and by a stronger deviation of the end state from the initial one, measured with a number of different norms; 2) a pronounced cyclone–anticyclone asymmetry during the nonlinear evolution of the moist-convective instability, which leads to an additional, with respect to the dry case, geostrophic adjustment, and the modification of the end state; and 3) an enhanced ageostrophic activity in the precipitation zones but also in the nonprecipitating areas because of the secondary geostrophic adjustment.
    [Show full text]
  • An Introduction to Dynamic Meteorology.Pdf
    January 24, 2004 12:0 Elsevier/AID aid An Introduction to Dynamic Meteorology FOURTH EDITION January 24, 2004 12:0 Elsevier/AID aid This is Volume 88 in the INTERNATIONAL GEOPHYSICS SERIES A series of monographs and textbooks Edited by RENATADMOWSKA, JAMES R. HOLTON and H.THOMAS ROSSBY A complete list of books in this series appears at the end of this volume. January 24, 2004 12:0 Elsevier/AID aid AN INTRODUCTION TO DYNAMIC METEOROLOGY Fourth Edition JAMES R. HOLTON Department of Atmospheric Sciences University of Washington Seattle,Washington Amsterdam Boston Heidelberg London New York Oxford Paris San Diego San Francisco Singapore Sydney Tokyo January 24, 2004 12:0 Elsevier/AID aid Senior Editor, Earth Sciences Frank Cynar Editorial Coordinator Jennifer Helé Senior Marketing Manager Linda Beattie Cover Design Eric DeCicco Composition Integra Software Services Printer/Binder The Maple-Vail Manufacturing Group Cover printer Phoenix Color Corp Elsevier Academic Press 200 Wheeler Road, Burlington, MA 01803, USA 525 B Street, Suite 1900, San Diego, California 92101-4495, USA 84 Theobald’s Road, London WC1X 8RR, UK This book is printed on acid-free paper. Copyright c 2004, Elsevier Inc. All rights reserved. No part of this publication may be reproduced or transmitted in any form or by any means, electronic or mechanical, including photocopy, recording, or any information storage and retrieval system, without permission in writing from the publisher. Permissions may be sought directly from Elsevier’s Science & Technology Rights Department in
    [Show full text]
  • Local and Nonlocal Advected Invariants and Helicities in MHD and Gas Dynamics I: Lie Dragging Approach, J
    Local and Nonlocal Advected Invariants and Helicities in Magnetohydrodynamics and Gas Dynamics II: Noether’s Theorems and Casimirs G.M. Webb1, B. Dasgupta1, J. F. McKenzie1,3, Q. Hu1,2 and G. P. Zank1,2 1 CSPAR, The University of Alabama Huntsville in Huntsville AL 35805, USA 2 Dept. of Physics, The University of Alabama in Huntsville, Huntsville AL 35899, USA 3Department of Mathematics and Statistics, Durban University of Technology, Steve Biko Campus, Durban South Africa and School of Mathematical Sciences, University of Kwa-Zulu, Natal, Durban, South Africa E-mail: [email protected] Abstract. Conservation laws in ideal gas dynamics and magnetohydrodynamics (MHD) associated with fluid relabelling symmetries are derived using Noether’s first and second theorems. Lie dragged invariants are discussed in terms of the MHD Casimirs. A nonlocal conservation law for fluid helicity applicable for a non-barotropic fluid involving Clebsch variables is derived using Noether’s theorem, in conjunction with a fluid relabelling symmetry and a gauge transformation. A nonlocal cross helicity conservation law involving Clebsch potentials, and the MHD energy conservation law are derived by the same method. An Euler Poincar´evariational approach is also used to derive conservation laws associated with fluid relabelling symmetries using Noether’s second theorem. arXiv:1307.1038v2 [math-ph] 8 Jan 2014 PACS numbers: 95.30.Qd,47.35.Tv,52.30.Cv,45.20Jj,96.60.j,96.60.Vg Submitted to: J. Phys. A., Math. and Theor.,28 June, 2013, Revised 4 June 2018 Keywords: magnetohydrodynamics, advection, invariants, Hamiltonian Advected Invariants in MHD and Fluids: Noether’s Theorems and Casimirs 2 1.
    [Show full text]