User talk:Paul Wormer/scratchbook: Difference between revisions
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 | 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 | 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 | 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 Δz and cross section ''O'' filled with air. | {{Image|Pressure on pilbox.png|right|250px|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 air. The pillbox is pushed | 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 | ||
:<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 ρ times the molar mass ''m'' times the volume of the pillbox Δ''z''·''O'' | The total mass in the pillbox is the molar density ρ times the molar mass ''m'' times the volume of the pillbox Δ''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} = - | \frac{1}{p} \frac{\Delta p}{\Delta z} = - \frac{mg}{RT} . | ||
</math> | </math> | ||
Taking the limit of | Taking the limit of infinitesimally small Δ 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} . | \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>) ≡ ''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., | 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 + | T \rightarrow T + L (z-z_0). | ||
</math> | </math> | ||
The modified equation is: | |||
:<math> | :<math> | ||
\int_{p_0}^{p} d\ln\,p' = - \frac{mg}{R} \int_{z_0}^z \frac{1}{T + | \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 + | 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}{ | - \frac{mg}{L R} \int^{T+L(z-z_0)}_T \frac{dy}{y} = | ||
- \frac{mg}{L R} \ln \frac{T+L(z-z_0)}{T} = | |||
- \frac{mg}{ | \ln\;\left[ \frac{T}{T+L(z-z_0)} \right]^{\frac{mg}{L R}}. | ||
\ln\;\left[ \frac{T}{T+ | |||
</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+ | \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
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