1 / 22

Richard B. Rood (Room 2525, SRB) rbrood@umich 734-647-3530 Cell: 301-526-8572

AOSS 401 Geophysical Fluid Dynamics: Atmospheric Dynamics Prepared: 20131010 Equations of Motion in Pressure Coordinates / Thermal Wind. Richard B. Rood (Room 2525, SRB) rbrood@umich.edu 734-647-3530 Cell: 301-526-8572. Class News. Ctools site ( AOSS 401 001 F13 )

Télécharger la présentation

Richard B. Rood (Room 2525, SRB) rbrood@umich 734-647-3530 Cell: 301-526-8572

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. AOSS 401Geophysical Fluid Dynamics:Atmospheric DynamicsPrepared: 20131010Equations of Motion in Pressure Coordinates / Thermal Wind Richard B. Rood (Room 2525, SRB) rbrood@umich.edu 734-647-3530 Cell: 301-526-8572

  2. Class News • Ctools site (AOSS 401 001 F13) • Second Examination on December 10, 2013 • Homework • Watch for Ctools / Email

  3. Weather • National Weather Service • Model forecasts: • Weather Underground • Model forecasts: • NCAR Research Applications Program

  4. Outline • Equations of Motion in Pressure Coordinates

  5. Return to Pressure Coordinates • In the coming weeks, we will work mostly with the equations of motion in pressure coordinates • We have already derived and used extensively the horizontal momentum equations in p-coordinates • Remaining equations: • Thermodynamic equation • Continuity equation • Vertical momentum equation (hydrostatic balance) • Plus, what to do with the material derivative?

  6. Our approximated horizontal momentum equations (in p coordinates) No viscosity, no metric terms, no cos-Coriolis terms Subscript h: horizontal Subscript p: constant p surfaces!Sometimes subscript is omitted,  tells you that this is on p surfaces

  7. Geostrophic Wind(in p coordinates) Component form: Vector form: Here we implicitly assume that the partial derivatives in x and y direction are computed on constant p surfaces. Advantage of the p coordinates: Simplicity of the equation!Density no longer appears explicitly.

  8. What do we do with the material derivative? By definition: Example: Total derivative of the temperature on constant z surfaces (subscript omitted)

  9. What do we do with the material derivative when using p in the vertical? By definition: Total derivative DT/Dt on constant pressure surfaces: (subscript omitted)

  10. Think about vertical motion in pressure coordinates • For upward motion, what is the sign of ‘w’ (vertical velocity in height coordinates)? • What is the sign of ω for upward vertical motion? Remember: • And pressure decreases with height… • What is the sign of Dp/Dt (change in pressure following the motion) for upward motion? For downward motion?

  11. Think about vertical motion in pressure coordinates • Something to keep in the back of your mind… • ω is negative for upward motion • ω is positive for downward motion

  12. Thermodynamic equation(in p coordinates) From Holton, we can Derive two forms of the thermodynamic equation Expand the material derivative and use

  13. Thermodynamic equation(in p coordinates) Collect terms in ω and divide through by cp Use the equation of state (ideal gas law)

  14. Thermodynamic equation(in p coordinates) Define the static stability parameterSp What is static stability?

  15. Adiabatic rising/sinking Diabatic heating (radiation, condensation) Thermodynamic equation(in p coordinates) If there is no horizontal advection, then the time rate of change of temperature is due to…? http://www.wunderground.com/modelmaps/maps.asp?model=NAM&domain=US

  16. Continuity equation in z coordinates: in p coordinates: We could try to derive the pressure-coordinate version from the height-coordinate equation (as we did with the pressure gradient force in the horizontal momentum equations), but…

  17. z y x Re-deriving the continuity equationin p coordinates It turns out to be easier to re-derive it from mass conservation Start with an air parcel withvolume:V= x y z Apply the hydrostatic equation p= -gz to express the volume element as V= - x y p/(g) The mass of this fluid element is:M = V= -  x y p/(g) = - x y p/g

  18. z y x Re-deriving the continuity equation in p coordinates (M = - x y p/g) Use the product rule Differential calculus Take the limit Continuity equation inp coordinates

  19. Continuity equation(in p coordinates) This form of the continuity equation contains no reference to the density field and does not involve time derivatives. The simplicity of this equation is one of the chief advantages of the isobaric system. Ease of computing vertical motion from convergence/divergence is another…

  20. Hydrostatic equation (in p coordinates) Hydrostatic equation in z-coords Rearrange (and assume g constant) Apply the equation of state Replaces the vertical momentum equation

  21. Approximated equations of motion in pressure coordinates

  22. In the derivation of the equations in pressure coordinates: • Have used conservation principles. • Have relied heavily on the hydrostatic assumption. • Required that conservation principles hold in all coordinate systems. • Plus we did some implicit scaling (metric terms, cosine coriolis terms dropped out).

More Related