DOCSLIB.ORG

Course Outline

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Course Outline These lecture notes are constantly being updated, extended and improved This is Version F2021:v0.1 Last Modified: Aug 14, 2021 2 Course Outline Astronomy is the study of everything in our Universe outside of the Earth’s atmo- sphere, including the Universe as a whole (‘cosmology’). As such, it involves all of physics: quantum mechanics, nuclear and particle physics, classical mechanics, special and general relativity, electromagnetism, statistical physics, hydrodynamics, plasma physics, and even solid state physics. As we will see, though, with the excep- tion of rocky planets and asteroids, all objects in the Universe can be characterized as some kind of fluid. In addition, almost all information we receive from these ob- jects reaches us in the form of radiation (photons)1. Consequently, in this course on physical processes in astronomy we will focus almost exclusively on fluid dynam- ics, including neutral fluids, charged fluids (aka plasmas), and collisionless fluids, and radiative processes. We start in Part I with standard hydrodynamics, which applies mainly to neutral, collisional fluids. We derive the fluid equations from kinetic theory: starting with the Liouville theorem we derive the Boltzmann equation, from which in turn we derive the continuity, momentum and energy equations. Next we discuss a variety of different flows; vorticity, incompressible barotropic flow, viscous flow, accretion flow, and turbulent flow, before addressing fluid instabilities and shocks. Next, in Part II, we briefly focus on collisionless dynamics, highlighting the subtle differences between the Jeans equations and the Navier-Stokes equations. We also give a brief account of orbit theory, and discuss the Virial theorem. In Part III we turn our attention to plasmas. We discuss how the relevant equations (Vlasov and Lenard-Balescu) derive from kinetic theory, and then discuss plasma orbit theory, magnetohydrodynamics (MHD) and several applications. Finally, in Part IV we discuss radiative processes and the interaction of light with matter (scattering & absorption). We end with a brief discussion of the cooling function, which describes how baryonic matter can dissipate its internal energy, leading to the formation of galaxies, stars and planets. It is assumed that the student is familiar with vector calculus, with curvi-linear coordinate systems, and with some radiative essential. A brief overview of these topics, highlighting the most important essentials, is provided in Appendices A-E and K. The other appendices present detailed background information that is provided for the interested student, but which is not considered part of the course material. 1Other astrophysical messengers include neutrinos, cosmic rays and gravitational waves 3 CONTENTS Part I: Fluid Dynamics 1: Introduction to Fluids and Plasmas ......................... ............9 2: Dynamical Treatments of Fluids ........................... ............16 3: Hydrodynamic Equations for Ideal Fluids .................... ..........25 4: Viscosity, Conductivity & The Stress Tensor . ..........28 5: Hydrodynamic Equations for Non-Ideal Fluids . .........34 6: Equations of State .................................... .................38 7: Kinetic Theory: from Liouville to Boltzmann . ......45 8: Kinetic Theory: from Boltzmann to Navier-Stokes . ........60 9: Vorticity & Circulation . .............71 10: Hydrostatics and Steady Flows ........................... .............79 11: Viscous Flow and Accretion Flow ............................ .........91 12: Turbulence .......................................... ................102 13: Sound Waves ......................................... ...............110 14: Shocks ............................................. ..................118 15: Fluid Instabilities . .............127 Part II: Collisionless Dynamics 16: Jeans equations & Dynamical Modeling ....................... .......138 17: Orbit Theory ......................................... ...............147 18: Virial Theorem & Gravothermal Catastrophe ................ ........151 19: Collisions & Encounters of Collisionless Systems . .... 156 Part III: Plasma Physics 20: Characteristic time & length scales . ..........167 21: Plasma Orbit Theory .................................... ............178 22: Plasma Kinetic Theory ................................... ...........186 23: Vlasov Equation & Two-Fluid Model .......................... ...... 192 24: Magnetohydrodynamics ................................ ..............199 Part IV: Radiative Processes 25: The Interaction of Light with Matter. I - Scattering ........... .......212 26: The Interaction of Light with Matter. II- Absorption .......... .......222 27: The Interaction of Light with Matter. III - Extinction . ......225 28: Radiative Transfer .................................... ...............229 29: Continuum Emission Processes ............................ ...........237 30: The Cooling Function.................................... ............249 4 APPENDICES Appendix A: Vector Calculus....................................... .........251 Appendix B: Conservative Vector Fields .............................. .......256 Appendix C: Integral Theorems................................... ...........257 Appendix D: Curvi-Linear Coordinate Systems ......................... .....258 Appendix E: The Levi-Civita Symbol .................................. .....268 Appendix F: The Viscous Stress Tensor .............................. .......269 Appendix G: Poisson Brackets ..................................... .........272 Appendix H: The BBGKY Hierarchy ................................... .... 274 Appendix I: Derivation of the Energy equation ......................... .....279 Appendix J: The Chemical Potential ................................. .......283 Appendix K: Radiation Essentials .................................... .......286 Appendix L: Velocity Moments of the Collision Integral . ..292 5 LITERATURE The material covered and presented in these lecture notes has relied heavily on a number of excellent textbooks listed below. The Physics of Fluids and Plasmas • by A. Choudhuri (ISBN-0-521-55543) The Physics of Astrophysics–I. Radiation • by F. Shu (ISBN-0-935702-64-4) The Physics of Astrophysics–II. Gas Dynamics • by F. Shu (ISBN-0-935702-65-2) Modern Fluid Dynamics for Physics and Astrophysics • by O. Regev, O. Umurhan & P. Yecko (ISBN-978-1-4939-3163-7) Principles of Astrophysical Fluid Dynamics • by C. Clarke & B.Carswell (ISBN-978-0-470-01306-9) Introduction to Plasma Theory • by D. Nicholson (ISBN-978-0-894-6467705) Introduction to Modern Magnetohydrodynamics • by S. Galtier (ISBN-978-1-316-69247-9) Astrophysics: Decoding the Cosmos • by J. Irwin (ISBN-978-0-470-01306-9) Theoretical Astrophysics • by M. Bartelmann (ISBN-978-3-527-41004-0) Radiative Processes in Astrophysics • by G Rybicki & A Lightman (ISBN-978-0-471-82759-7) Galactic Dynamics • by J. Binney & S. Tremaine (ISBN-978-0-691-13027-9) Modern Classical Physics • by K. Thorne & R. Blandford (ISBN-978-0-691-15902-7) 6 7 Part I: Fluid Dynamics Almost everything we encounter in the Universe, from gas planets to stars, and from the interstellar medium to galaxies, can be categorized as some kind of fluid. Hence, understanding astrophysical processes requires a solid understanding of fluid dynamics. The following chapters present fairly detailed description of the dynamics of fluids with an application to astrophysics. Fluid dynamics is a rich topic, and one could easily devote an entire course to it. The following chapters therefore only scratch the surface of this rich topic. Readers who want to get more indepth information are referred to the following excellent textbooks - The Physics of Fluids and Plasmas by A. Choudhuri - Modern Fluid Dynamics for Physics and Astrophysics by O.Regev et al. - The Physics of Astrophysics II. Gas Dynamics by F. Shu - Principles of Astrophysical Fluid Dynamics by C. Clarke & B. Carswell - Modern Classical Physics by K.Thorne & R. Blandford 8 CHAPTER 1 Introduction to Fluids & Plasmas What is a fluid? A fluid is a substance that can flow, has no fixed shape, and offers little resistance to an external stress In a fluid the constituent particles (atoms, ions, molecules, stars) can ‘freely’ • move past one another. Fluids take on the shape of their container (if not self-gravitating). • A fluid changes its shape at a steady rate when acted upon by a stress force. • What is a plasma? A plasma is a fluid in which (some of) the consistituent particles are electrically charged, such that the interparticle force (Coulomb force) is long-range in nature. Fluid Demographics: All fluids are made up of large numbers of constituent particles, which can be molecules, atoms, ions, dark matter particles or even stars. Different types of fluids mainly differ in the nature of their interparticle forces. Examples of inter-particle forces are the Coulomb force (among charged particles in a plasma), vanderWaals forces (among molecules in a neutral fluid) and gravity (among the stars in a galaxy). Fluids can be both collisional or collisionless, where we define a collision as an interaction between constituent particles that causes the trajectory of at least one of these particles to be deflected ‘noticeably’. Collisions among particles drive the system towards thermodynamic equilibrium (at least locally) and the velocity distribution towards a Maxwell-Boltzmann distribution. In neutral fluids the particles only interact with each other on very small scales. Typically the inter-particle
Load more

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]
  • A Brief Tour of Vector Calculus
    A BRIEF TOUR OF VECTOR CALCULUS A. HAVENS Contents 0 Prelude ii 1 Directional Derivatives, the Gradient and the Del Operator 1 1.1 Conceptual Review: Directional Derivatives and the Gradient........... 1 1.2 The Gradient as a Vector Field............................ 5 1.3 The Gradient Flow and Critical Points ....................... 10 1.4 The Del Operator and the Gradient in Other Coordinates*............ 17 1.5 Problems........................................ 21 2 Vector Fields in Low Dimensions 26 2 3 2.1 General Vector Fields in Domains of R and R . 26 2.2 Flows and Integral Curves .............................. 31 2.3 Conservative Vector Fields and Potentials...................... 32 2.4 Vector Fields from Frames*.............................. 37 2.5 Divergence, Curl, Jacobians, and the Laplacian................... 41 2.6 Parametrized Surfaces and Coordinate Vector Fields*............... 48 2.7 Tangent Vectors, Normal Vectors, and Orientations*................ 52 2.8 Problems........................................ 58 3 Line Integrals 66 3.1 Defining Scalar Line Integrals............................. 66 3.2 Line Integrals in Vector Fields ............................ 75 3.3 Work in a Force Field................................. 78 3.4 The Fundamental Theorem of Line Integrals .................... 79 3.5 Motion in Conservative Force Fields Conserves Energy .............. 81 3.6 Path Independence and Corollaries of the Fundamental Theorem......... 82 3.7 Green's Theorem.................................... 84 3.8 Problems........................................ 89 4 Surface Integrals, Flux, and Fundamental Theorems 93 4.1 Surface Integrals of Scalar Fields........................... 93 4.2 Flux........................................... 96 4.3 The Gradient, Divergence, and Curl Operators Via Limits* . 103 4.4 The Stokes-Kelvin Theorem..............................108 4.5 The Divergence Theorem ...............................112 4.6 Problems........................................114 List of Figures 117 i 11/14/19 Multivariate Calculus: Vector Calculus Havens 0.
    [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 1 VECTOR ANALYSIS ̂ ̂ Distribution Rule Vector Components
    8/29/2016 CHAPTER 1 1. VECTOR ALGEBRA VECTOR ANALYSIS Addition of two vectors 8/24/16 Communicative Chap. 1 Vector Analysis 1 Vector Chap. Associative Definition Multiplication by a scalar Distribution Dot product of two vectors Lee Chow · , & Department of Physics · ·· University of Central Florida ·· 2 Orlando, FL 32816 • Cross-product of two vectors Cross-product follows distributive rule but not the commutative rule. 8/24/16 8/24/16 where , Chap. 1 Vector Analysis 1 Vector Chap. Analysis 1 Vector Chap. In a strict sense, if and are vectors, the cross product of two vectors is a pseudo-vector. Distribution rule A vector is defined as a mathematical quantity But which transform like a position vector: so ̂ ̂ 3 4 Vector components Even though vector operations is independent of the choice of coordinate system, it is often easier 8/24/16 8/24/16 · to set up Cartesian coordinates and work with Analysis 1 Vector Chap. Analysis 1 Vector Chap. the components of a vector. where , , and are unit vectors which are perpendicular to each other. , , and are components of in the x-, y- and z- direction. 5 6 1 8/29/2016 Vector triple products · · · · · 8/24/16 8/24/16 · ( · Chap. 1 Vector Analysis 1 Vector Chap. Analysis 1 Vector Chap. · · · All these vector products can be verified using the vector component method. It is usually a tedious process, but not a difficult process. = - 7 8 Position, displacement, and separation vectors Infinitesimal displacement vector is given by The location of a point (x, y, z) in Cartesian coordinate ℓ as shown below can be defined as a vector from (0,0,0) 8/24/16 In general, when a charge is not at the origin, say at 8/24/16 to (x, y, z) is given by (x’, y’, z’), to find the field produced by this Chap.
    [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]
  • Vector Calculus Operators
    Vector Calculus Operators July 6, 2015 1 Div, Grad and Curl The divergence of a vector quantity in Cartesian coordinates is defined by: −! @ @ @ r · ≡ x + y + z (1) @x @y @z −! where = ( x; y; z). The divergence of a vector quantity is a scalar quantity. If the divergence is greater than zero, it implies a net flux out of a volume element. If the divergence is less than zero, it implies a net flux into a volume element. If the divergence is exactly equal to zero, then there is no net flux into or out of the volume element, i.e., any field lines that enter the volume element also exit the volume element. The divergence of the electric field and the divergence of the magnetic field are two of Maxwell's equations. The gradient of a scalar quantity in Cartesian coordinates is defined by: @ @ @ r (x; y; z) ≡ ^{ + |^ + k^ (2) @x @y @z where ^{ is the unit vector in the x-direction,| ^ is the unit vector in the y-direction, and k^ is the unit vector in the z-direction. Evaluating the gradient indicates the direction in which the function is changing the most rapidly (i.e., slope). If the gradient of a function is zero, the function is not changing. The curl of a vector function in Cartesian coordinates is given by: 1 −! @ @ @ @ @ @ r × ≡ ( z − y )^{ + ( x − z )^| + ( y − x )k^ (3) @y @z @z @x @x @y where the resulting function is another vector function that indicates the circulation of the original function, with unit vectors pointing orthogonal to the plane of the circulation components.
    [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]
  • Divergence, Gradient and Curl Based on Lecture Notes by James
    Divergence, gradient and curl Based on lecture notes by James McKernan One can formally define the gradient of a function 3 rf : R −! R; by the formal rule @f @f @f grad f = rf = ^{ +| ^ + k^ @x @y @z d Just like dx is an operator that can be applied to a function, the del operator is a vector operator given by @ @ @ @ @ @ r = ^{ +| ^ + k^ = ; ; @x @y @z @x @y @z Using the operator del we can define two other operations, this time on vector fields: Blackboard 1. Let A ⊂ R3 be an open subset and let F~ : A −! R3 be a vector field. The divergence of F~ is the scalar function, div F~ : A −! R; which is defined by the rule div F~ (x; y; z) = r · F~ (x; y; z) @f @f @f = ^{ +| ^ + k^ · (F (x; y; z);F (x; y; z);F (x; y; z)) @x @y @z 1 2 3 @F @F @F = 1 + 2 + 3 : @x @y @z The curl of F~ is the vector field 3 curl F~ : A −! R ; which is defined by the rule curl F~ (x; x; z) = r × F~ (x; y; z) ^{ |^ k^ = @ @ @ @x @y @z F1 F2 F3 @F @F @F @F @F @F = 3 − 2 ^{ − 3 − 1 |^+ 2 − 1 k:^ @y @z @x @z @x @y Note that the del operator makes sense for any n, not just n = 3. So we can define the gradient and the divergence in all dimensions. However curl only makes sense when n = 3. Blackboard 2. The vector field F~ : A −! R3 is called rotation free if the curl is zero, curl F~ = ~0, and it is called incompressible if the divergence is zero, div F~ = 0.
    [Show full text]
  • 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.
    [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]