Primitive equations
From Wikipedia, the free encyclopedia
The primitive equations are a version of the Navier-Stokes equations which describe hydrodynamical flow on the sphere under the assumptions that vertical motion is much smaller than horizontal motion (hydrostasis) and that the fluid layer depth is small compared to the radius of the sphere. Thus, they are a good approximation of global atmospheric flow and are used in most atmospheric models. The primitive equations may be linearized to yield Laplace's tidal equations, an eigenvalue problem from which the analytical solution to the latitudinal structure of the flow may be determined.
In general, nearly all forms of the primitive equations relate the five variables (u,v,ω,T,φ), and their evolution over space and time.
Contents |
[edit] Definitions
- u is the zonal velocity (velocity in the east/west direction tangent to the sphere).
- v is the meridional velocity (velocity in the north/south direction tangent to the sphere).
- ω is the vertical velocity
- T is the temperature
- φ is the geopotential
- f is the term corresponding to the Coriolis force, and is equal to 2Ωsin(φ), where Ω is the angular rotation rate of the Earth (2π / 24 radians/hour), and φ is the latitude.
- R is the gas constant
- p is the pressure
- cp is the specific heat
- J is the heat flow per unit time per unit mass
- π is the exner function
- θ is the potential temperature
[edit] Forms of the primitive equations
The precise form of the primitive equations depends on the vertical coordinate system chosen, such as pressure coordinates, log pressure coordinates, or sigma coordinates. Furthermore, the velocity, temperature, and geopotential variables may be decomposed into mean and perturbation components using Reynolds decomposition.
[edit] Vertical pressure, cartesian tangential plane
In this form pressure is selected as the vertical coordinate and the horizontal coordinates are written for the cartesian tangential plane (i.e. a plane tangent to some point on the surface of the Earth). This form does not take the curvature of the Earth into account, but is useful for visualizing some of the physical processes involved in formulating the equations due to its relative simplicity.
Note that the capital derivatives are the material derivatives.
- the geostrophic momentum equations
- the hydrostatic equation, a special case of the vertical momentum equation in which there is no background vertical acceleration.
- the continuity equation, connecting horizontal divergence/convergence to vertical motion
- and the Thermodynamic Energy equation, a consequence of the first law of thermodynamics
When a statement of the conservation of water vapor substance is included, these six equations form the basis for any numerical weather prediction scheme.
[edit] Primitive equations using sigma coordinate system, polar stereographic projection
- According to the National Weather Service Handbook No. 1 - Facsimile Products, the primitive equations can be simplified into the following equations:
- Temperature: ∂T/∂t = u (∂Tx/∂X) + v (∂Ty/∂Y) + w (∂Tz/∂Z)
- Wind in E-W direction: ∂u/∂t = ηv - ∂Φ/∂x – Cp θ (∂π/∂x) – z (∂u/∂σ) – [∂(u2 + y) / 2] / ∂x
- Wind in N-S direction: ∂v/∂t = -η(u/v) - ∂Φ/∂y – Cp θ (∂π/∂y) – z (∂v/∂σ) – [∂(u2 + y) / 2] / ∂y
- Precipitable water: ∂W/∂t = u (∂Wx/∂X) + v (∂Wy/∂Y) + z (∂Wz/∂Z)
- Pressure Thickness: ∂(∂p/∂σ)/∂t = u [(∂p/∂σ)x /∂X] + v [(∂p/∂σ)y /∂Y] + z [(∂p/∂σ)z /∂Z]
- These simplifications make it much easier to understand what is happening in the model. Things like the temperature (potential temperature), precipitable water, and to an extent the pressure thickness simply move from one spot on the grid to another with the wind. The wind is forecasted slightly differently. It uses geopotential, specific heat, the exner function π, and change in sigma coordinate.
[edit] Solution to the primitive equations
The analytic solution to the primitive equations involves a sinusoidal oscillation in time and longitude, modulated by coefficients related to height and latitude.
s and σ are the zonal wavenumber and angular frequency, respectively. The solution represents atmospheric waves and tides.
When the coefficients are separated into their height and latitude components, the height dependence takes the form of propagating or evanescent waves (depending on conditions), while the latitude dependence is given by the Hough functions.
This analytic solution is only possible when the primitive equations are linearized and simplified. Unfortunately many of these simplifications (i.e. no dissipation, isothermal atmosphere) do not correspond to conditions in the actual atmosphere. As a result, a numerical solution which takes these factors into account is often calculated using general circulation models and climate models.
[edit] How the wind vectors, u and v, are calculated
- Normally, wind consists of two components: speed and direction
- In the model, wind still consists the speed in N-S direction and speed in E-W direction
- Speed in E-W direction is the variable u.
- Speed in N-S direction is the variable v.
- Wind the E-W direction is:
- u = W (cos d)
- Wind in the N-S direction is:
- v = W (sin d)
- Where W is the wind speed and d is the direction in degrees
- It should be noted that although the wind in meteorology is measured in degrees from the north (0 is north, 90 is east, 180 is south, 270 is west), in trigonometry degrees are measured starting at east at 0, north is 90, west is 180, south is 270. Therefore the wind direction must be converted to this standard to be calculated on a computer or calulator.
- This can be corrected by using the equation -θ + 90. If the resulting number is below 0, add 360. If the resulting number is above 360, subtract 360. You should get a number between 0 and 360.
[edit] Forecasting temperature
In this form of computer model, the change in temperature at a given location is the easiest thing to forecast. The temperature, given in potential temperature, is simply moved from one spot to another. How?
- If temperature changes by 2 degrees every 1 mile north and by 1 degree every 1 mile east:
- If you multiply the rate of change in temperature per mile north by the wind speed in the north-south direction in mph, you get the change in temperature per hour contributed from this direction.
- If you multiply the rate of change in temperature per mile east by the wind speed in the east-west direction in mph, you get the change in temperature per hour contributed from that direction.
- Adding the change given from both directions gives us the total change in temperature per hour. This total change in temperature added to the original temperature will give us our forecast temperature.
[edit] Forecasting precipitable water and pressure thickness
- The format explained above for calculating the change in temperature is exactly the way the change in precipitable water and pressure thickness can be forecasted. The most difficult calculations involve forecasting the change in wind between run intervals.
[edit] Useful references
- Beniston, Martin. From Turbulence to Climate: Numerical Investigations of the Atmosphere with a Hierarchy of Models. Berlin: Springer, 1998.
- Firth, Robert. Mesoscale and Microscale Meteorological Model Grid Construction and Accuracy. LSMSA, 2006.
- Thompson, Philip. Numerical Weather Analysis and Prediction. New York: The Macmillan Company, 1961.
- Pielke, Roger A. Mesoscale Meteorological Modeling. Orlando: Academic Press, Inc., 1984.
- U.S. Department of Commerce, National Oceanic and Atmospheric Administration, National Weather Service. National Weather Service Handbook No. 1 - Facsimile Products. Washington, DC: Department of Commerce, 1979.
