User talk:Paul Wormer/scratchbook: Difference between revisions

From Citizendium
Jump to navigation Jump to search
imported>Paul Wormer
No edit summary
imported>Paul Wormer
No edit summary
Line 1: Line 1:
A  '''barometric formula''' gives the pressure in a planetary atmosphere—usually the [[Earth's atmosphere]]—as a function of height. Since the [[temperature]] and [[composition]] of the atmosphere may be a complicated function of height, a  barometric formula may be very intricate.  
A  '''barometric formula''' gives the pressure in a planetary atmosphere—usually the [[Earth's atmosphere|Earth's]]—as a function of height. Since the [[temperature]] and [[composition]] of the atmosphere are  complicated functions of height, and because gravitation is  an inverse function of distance, a completely general barometric formula is very intricate.  


However, two simplified formulas are in common use. The first one assumes constant temperature and composition throughout the atmosphere. It assumes that the atmosphere is an [[ideal gas law|ideal gas]]. The equation is a form of a [[Boltzmann distribution]],
However, two simplified formulas are in common use. The first is based on the assumption that temperature, gravitation,  and composition are constant throughout the atmosphere. Further it is used that the atmosphere is an [[ideal gas law|ideal gas]]. The equation thus obtained is a form of a [[Boltzmann distribution]],
:<math>
:<math>
p(z) = p(z_0) \; \exp[ - \frac{mg}{RT} (z-z_0)]
p(z) = p(z_0) \; \exp\left[ - \frac{mg}{RT} (z-z_0)\right].
</math>
</math>
Here ''p''(''z'')  is the pressure at height ''z''; ''z''<sub>0/sub> is a suitable reference height, often the surface of the planet; ''m'' is the  mass of one mole of the atmospheric gas (note that the equation holds for one substance only, if the atmosphere contains more than one compound a weighted mass is used); ''g'' is the [[gravitational acceleration]], ''R'' is the [[molar gas constant]]; ''T'' is the absolute temperature (assumed independent of ''z'').
Here ''p''(''z'')  is the pressure at height ''z''; ''z''<sub>0</sub> is a suitable reference height, often the surface of the planet and usually referred to as "base"; ''m'' is the  mass of one mole of the atmospheric gas (note that the equation holds for a pure gas, if the atmosphere contains more than one compound a weighted-average mass is used); ''g'' is the constant [[gravitational acceleration]]; ''R'' is the [[molar gas constant]]; ''T'' is the absolute [[temperature]] (assumed to be independent of ''z'').


The second barometric formula is based on the same assumptions.
The second barometric formula is based on the same assumptions as the first, the only generalization being that the temperature is taken to be linear in  height with a slope ''L'', the [[lapse rate]].
==Proof==  
The equation is
:<math>
p(z) = p(z_0) \; \left[ \frac{T}{T + L (z-z_0)}\right]^{mg/(LR)}.
</math>
==Proof==
The two equations follow by integration of a simple differential equation that holds for the pressure of an ideal gas in a constant gravitational field.
===Differential equation===
===Differential equation===
{{Image|Pressure on pilbox.png|right|250px|A pillbox  of height &Delta;z and cross section ''O'' filled with air. Pillbox,  with upward pressure ''p'' at its bottom, is part of a cylinder filled with air.  Total mass ''M'' of air in pillbox is pulled downward by gravitational acceleration ''g'. }}
{{Image|Pressure on pilbox.png|right|250px|A pillbox  of height &Delta;z and cross section ''O'' filled with air. The pillbox,  with upward pressure ''p'' at its bottom, is part of a cylinder filled with air.  Total mass ''M'' of air in pillbox is pulled downward by gravitational acceleration ''g'. }}


In the figure we see a stationary pillbox at arbitrary position ''z'' in a cylinder filled with air. The pillbox is pushed up by a [[pressure]] force ''pO'', and pulled down by gravitation, which gives a force '''G''' equal to the total mass ''M'' of air in the pillbox times the [[gravitational acceleration]] ''g''. Also the air pressure (''p''+ &Delta;''p'')''O'' at the top acts downward. Forces upward are positive, downward negative, and since the total force is zero (no motion of air), we have
In the figure we see a stationary pillbox at arbitrary position ''z'' in a cylinder filled with an atmospheric gas that, for obvious reason,  we  call "air". The pillbox is pushed upward by a [[pressure]] force ''pO'', and pulled down by gravitation that acts with the force '''G''' equal to the total mass ''M'' of air in the pillbox times the [[gravitational acceleration]] ''g''. Also the air pressure (''p''+ &Delta;''p'')''O'' at the top acts downward. Forces upward are positive, downward negative, and since the total force is zero (no motion of air), we have
:<math>
:<math>
pO -Mg - (p+\Delta p)O = 0 \;\Longrightarrow \; \Delta p = - \frac{Mg}{O}
pO -Mg - (p+\Delta p)O = 0 \;\Longrightarrow \; \Delta p = - \frac{Mg}{O}.
</math>
</math>
The total mass in the pillbox is the molar density &rho; times the molar mass ''m'' times the volume of the pillbox &Delta;''z''·''O''
The total mass in the pillbox is the molar density &rho; times the molar mass ''m'' times the volume of the pillbox &Delta;''z''·''O''
:<math>
:<math>
M = \rho m \Delta z O .\;
M = \rho\, m\, \Delta z O .\;
</math>
</math>
Hence
Hence
:<math>
:<math>
\Delta p = -  \rho m g \Delta z \,
\Delta p = -  \rho\, m\, g\, \Delta z .\,
</math>
</math>
Noting that the molar density is the number of moles per volume, we get for an [[ideal gas law|ideal gas]],
Noting that the molar density is the number of moles per volume, we get for an [[ideal gas law|ideal gas]],
:<math>
:<math>
p = \rho R T \;\Longrightarrow\; \rho = \frac{p}{RT}
p = \rho R T \;\Longrightarrow\; \rho = \frac{p}{RT}.
</math>
</math>
Assuming that the gas in the pillbox is ideal and of absolute [[temperature]] ''T'', we get
Assuming that the gas in the pillbox is ideal and of absolute [[temperature]] ''T'', we get
Line 32: Line 37:
:<math>
:<math>
\Delta p = - p \frac{mg}{RT} \Delta z \;\Longrightarrow\;  
\Delta p = - p \frac{mg}{RT} \Delta z \;\Longrightarrow\;  
\frac{1}{p} \frac{\Delta p}{\Delta z} = - p \frac{mg}{RT}
\frac{1}{p} \frac{\Delta p}{\Delta z} = - \frac{mg}{RT} .
</math>  
</math>  
Taking  the limit of infinitesimal small  &Delta;  and
Taking  the limit of infinitesimally small  &Delta;  and
remembering that the derivative of the natural logarithm is given by
remembering that the derivative of the natural logarithm is given by
:<math>
:<math>
Line 41: Line 46:
we arrive at the the differential equation
we arrive at the the differential equation
:<math>
:<math>
  \frac{d\ln\,p}{dz} = - \frac{mg}{RT} . \qquad\qquad\qquad\qquad\qquad\qquad\qquad(1)
  \frac{d\ln\,p(z)}{dz} = - \frac{mg}{RT} . \qquad\qquad\qquad\qquad\qquad\qquad(1)
</math>
</math>


===Constant temperature===
===Constant temperature===
The master equation (1) can be solved by integrating the left hand side from ''p''<sub>0</sub> to  
The master equation (1) can be solved by integrating the left hand side from ''p''(''z''<sub>0</sub>) &equiv; ''p''<sub>0</sub> to  
''p'' and the right hand side from ''z''<sub>0</sub> to ''z''. Hence
''p''(''z'') and the right hand side from ''z''<sub>0</sub> to ''z''. Hence
:<math>
:<math>
\ln\frac{p}{p_0} = - \frac{mg}{RT} (z-z_0) \qquad\qquad\qquad\qquad\qquad(2)
\ln\frac{p}{p_0} = - \frac{mg}{RT} (z-z_0), \qquad\qquad\qquad\qquad\qquad(2)
</math>
</math>
where we recall that  
where we recall that  
:<math>
:<math>
\ln p - \ln p_0 = \ln\frac{p}{p_0}
\ln p - \ln p_0 = \ln\frac{p}{p_0}.
</math>
</math>
Taking the exponent (the inverse function of the natural logarithm) of both sides of equation (2)
Taking the exponent (the inverse function of the natural logarithm) of both sides of equation (2)
Line 59: Line 64:
\frac{p}{p_0} = \frac{\exp\left[-  \frac{mg}{RT} z\right]}{\exp\left[-  \frac{mg}{RT} z_0\right]}
\frac{p}{p_0} = \frac{\exp\left[-  \frac{mg}{RT} z\right]}{\exp\left[-  \frac{mg}{RT} z_0\right]}
\;\Longleftrightarrow\;
\;\Longleftrightarrow\;
p = p_0 \, \exp\left[ - \frac{mg}{RT}(z-z_0) \right]
p = p_0 \, \exp\left[ - \frac{mg}{RT}(z-z_0) \right].
</math>
</math>
This  distribution is known in [[statistical thermodynamics]] as a [[Boltzmann distribution]].
This  distribution is known in [[statistical thermodynamics]] as a [[Boltzmann distribution]].
===Temperature linear in height===
===Temperature linear in height===
We assumed in the previous subsection that ''T'' is constant throughout the cylinder, which makes the equation useful for cases where there is no temperature gradient, i.e., case where the temperature at the top of the cylinder is equal to the temperature at the bottom.  The master equation (1) can be modified for the case of a constant temperature gradient (constant [[lapse rate]]). That is, we  make the substitution
We assumed in the previous subsection that ''T'' is constant throughout the cylinder, which makes the equation useful for cases where there is no temperature gradient, i.e., cases where the temperature at the top of the cylinder is equal to the temperature at the bottom.  The master equation (1) can be modified for the case of a constant temperature gradient (constant [[lapse rate]]). That is, we  make in equation (1) the substitution
:<math>
:<math>
T \rightarrow T + K (z-z_0)
T \rightarrow T + L (z-z_0).
</math>  
</math>  
in equation (1). The modified equation is still easily integrated:
The modified equation is:
:<math>
:<math>
\int_{p_0}^{p} d\ln\,p'  = - \frac{mg}{R} \int_{z_0}^z \frac{1}{T + K (z'-z_0)} dz'   
\int_{p_0}^{p} d\ln\,p'  = - \frac{mg}{R} \int_{z_0}^z \frac{1}{T + L (z'-z_0)} dz'  .
</math>
</math>
Substitute  
Substitute  
:<math>
:<math>
y = T +K(z'-z_0)\; \Longrightarrow \; dz'= \frac{1}{K} dy
y = T +L(z'-z_0)\; \Longrightarrow \; dz'= \frac{1}{L} dy .
</math>
</math>
The integral on the right hand side becomes
The integral on the right hand side becomes
:<math>
:<math>
- \frac{mg}{R} \int_{z_0}^z \frac{1}{T + K (z'-z_0)} dz' =
- \frac{mg}{L R} \int^{T+L(z-z_0)}_T \frac{dy}{y}  =
- \frac{mg}{K R} \int^{T+K(z-z_0)}_T \frac{dy}{y}  =
- \frac{mg}{L R} \ln \frac{T+L(z-z_0)}{T} =  
- \frac{mg}{K R} \ln \frac{T+K(z-z_0)}{T} =  
\ln\;\left[ \frac{T}{T+L(z-z_0)} \right]^{\frac{mg}{L R}}.
\ln\;\left[ \frac{T}{T+K(z-z_0)} \right]^{\frac{mg}{K R}}
</math>
</math>
We already saw that the integral on the left hand side is the logarithm of ''p''/''p''<sub>0</sub>. hence
We already saw that the integral on the left hand side is the logarithm of ''p''/''p''<sub>0</sub>. hence
:<math>
:<math>
\frac{p}{p_0} = \left[ \frac{T}{T+K(z-z_0)} \right]^{\frac{mg}{K R}}.
\frac{p}{p_0} = \left[ \frac{T}{T+L(z-z_0)} \right]^{\frac{mg}{L R}}.
</math>
</math>

Revision as of 03:00, 27 August 2009

A barometric formula gives the pressure in a planetary atmosphere—usually the Earth's—as a function of height. Since the temperature and composition of the atmosphere are complicated functions of height, and because gravitation is an inverse function of distance, a completely general barometric formula is very intricate.

However, two simplified formulas are in common use. The first is based on the assumption that temperature, gravitation, and composition are constant throughout the atmosphere. Further it is used that the atmosphere is an ideal gas. The equation thus obtained is a form of a Boltzmann distribution,

Here p(z) is the pressure at height z; z0 is a suitable reference height, often the surface of the planet and usually referred to as "base"; m is the mass of one mole of the atmospheric gas (note that the equation holds for a pure gas, if the atmosphere contains more than one compound a weighted-average mass is used); g is the constant gravitational acceleration; R is the molar gas constant; T is the absolute temperature (assumed to be independent of z).

The second barometric formula is based on the same assumptions as the first, the only generalization being that the temperature is taken to be linear in height with a slope L, the lapse rate. The equation is

Proof

The two equations follow by integration of a simple differential equation that holds for the pressure of an ideal gas in a constant gravitational field.

Differential equation

PD Image
A pillbox of height Δz and cross section O filled with air. The pillbox, with upward pressure p at its bottom, is part of a cylinder filled with air. Total mass M of air in pillbox is pulled downward by gravitational acceleration g'.

In the figure we see a stationary pillbox at arbitrary position z in a cylinder filled with an atmospheric gas that, for obvious reason, we call "air". The pillbox is pushed upward by a pressure force pO, and pulled down by gravitation that acts with the force G equal to the total mass M of air in the pillbox times the gravitational acceleration g. Also the air pressure (p+ Δp)O at the top acts downward. Forces upward are positive, downward negative, and since the total force is zero (no motion of air), we have

The total mass in the pillbox is the molar density ρ times the molar mass m times the volume of the pillbox Δz·O

Hence

Noting that the molar density is the number of moles per volume, we get for an ideal gas,

Assuming that the gas in the pillbox is ideal and of absolute temperature T, we get the equation

Taking the limit of infinitesimally small Δ and remembering that the derivative of the natural logarithm is given by

we arrive at the the differential equation

Constant temperature

The master equation (1) can be solved by integrating the left hand side from p(z0) ≡ p0 to p(z) and the right hand side from z0 to z. Hence

where we recall that

Taking the exponent (the inverse function of the natural logarithm) of both sides of equation (2) we get

This distribution is known in statistical thermodynamics as a Boltzmann distribution.

Temperature linear in height

We assumed in the previous subsection that T is constant throughout the cylinder, which makes the equation useful for cases where there is no temperature gradient, i.e., cases where the temperature at the top of the cylinder is equal to the temperature at the bottom. The master equation (1) can be modified for the case of a constant temperature gradient (constant lapse rate). That is, we make in equation (1) the substitution

The modified equation is:

Substitute

The integral on the right hand side becomes

We already saw that the integral on the left hand side is the logarithm of p/p0. hence