View source for Brunt–Väisälä frequency

{{Short description|Measure of fluid stability against vertical displacement}}
{{Other uses|Vaisala (disambiguation)}}
In [[atmospheric dynamics]], [[oceanography]], [[asteroseismology]] and [[geophysics]], the '''Brunt–Väisälä frequency''', or '''[[buoyancy]] frequency''', is a measure of the stability of a fluid to vertical displacements such as those caused by [[convection]]. More precisely it is the frequency at which a vertically displaced parcel will oscillate within a statically stable environment. It is named after [[David Brunt]] and [[Vilho Väisälä]]. It can be used as a measure of atmospheric stratification.

==Derivation for a general fluid==
Consider a parcel of water or gas that has density <math>\rho_0</math>. This parcel is in an environment of other water or gas particles where the density of the environment is a function of height: <math>\rho = \rho (z)</math>. If the parcel is displaced by a small vertical increment <math>z'</math>, ''and it maintains its original density, so that its volume does not change,'' it will be subject to an extra gravitational force against its surroundings of:

:<math>\rho_0 \frac{\partial^2 z'}{\partial t^2} = - g \left[\rho (z)-\rho (z+z')\right]</math>

where <math>g</math> is the gravitational acceleration, and is defined to be positive. We make a [[linear approximation]] to <math>\rho (z+z') - \rho (z) = \frac{\partial \rho (z)}{\partial z} z'</math>, and move <math>\rho_0</math> to the RHS:

:<math>\frac{\partial^2 z'}{\partial t^2} = \frac{g}{\rho_0} \frac{\partial \rho (z)}{\partial z} z' </math>

The above second-order [[differential equation]] has the following solution:

:<math>z' = z'_0 e^{i N t}\!</math>

where the Brunt–Väisälä frequency <math>N</math> is:<ref>{{Cite book|title=Atmospheric and oceanic fluid dynamics: fundamentals and large-scale circulation|last=Vallis|first=Geoffrey K.|isbn=9781107588417|edition=2nd|location=Cambridge|publisher=[[Cambridge University Press]]|date=2017|oclc=990033511|doi=10.1017/9781107588417|bibcode=2017aofd.book.....V |s2cid=115699889 }}</ref>

:<math>N = \sqrt{- \frac{g}{\rho_0} \frac{\partial \rho (z)}{\partial z}}</math>

For negative <math>\frac{\partial \rho (z)}{\partial z}</math>, the displacement <math>z'</math> has oscillating solutions (and N gives our angular frequency). If it is positive, then there is run away growth – i.e. the fluid is statically unstable.

==In meteorology and astrophysics==
For a gas parcel, the density will only remain fixed as assumed in the previous derivation if the pressure, <math>P</math>, is constant with height, which is not true in an atmosphere confined by gravity. Instead, the parcel will expand adiabatically as the pressure declines. Therefore a more general formulation used in meteorology is:

: <math>N \equiv \sqrt{\frac{g}{\theta}\frac{d\theta}{dz}}</math>, where <math>\theta</math> is [[potential temperature]], <math>g</math> is the local acceleration of [[gravity]], and <math>z</math> is [[geometric height]].<ref>{{Cite book|title=Atmospheric Convection|last=Emmanuel|first=K.A.|publisher=Oxford University Press|year=1994|isbn=0195066308|doi=10.1002/joc.3370150709}}</ref>

Since <math>\theta = T (P_0/P)^{R/c_P}</math>, where <math>P_0</math> is a constant reference pressure, for a perfect gas this expression is equivalent to:
: <math>N^2 \equiv g\left\{\frac{1}{T}\frac{dT}{dz}-\frac{R}{c_P}\frac{1}{P}\frac{dP}{dz}\right\}=g\left\{\frac{1}{T}\frac{dT}{dz}-\frac{\gamma - 1}{\gamma}\frac{1}{P}\frac{dP}{dz}\right\}</math>,
where in the last form  <math>\gamma = c_P/c_V</math>, the [[adiabatic index]]. Using the [[ideal gas law]], we can eliminate the temperature to express <math>N^2</math> in terms of pressure and density: 
: <math>N^2 \equiv g\left\{\frac{1}{\gamma}\frac{1}{P}\frac{dP}{dz}-\frac{1}{\rho}\frac{d\rho}{dz}\right\}=g\left\{\frac{1}{\gamma}\frac{d\ln P}{dz}-\frac{d\ln\rho}{dz}\right\}</math>.

This version is in fact more general than the first, as it applies when the chemical composition of the gas varies with height, and also for imperfect gases with variable adiabatic index, in which case <math>\gamma \equiv \gamma_{01}= (\partial \ln P / \partial \ln \rho)_{S}</math>, i.e. the derivative
is taken at constant [[entropy]], <math>S</math>.<ref>{{Citation|first=Jørgen|last=Christensen-Dalsgaard|year=2014|title=Lecture Notes on Stellar Oscillations|edition=5th|url=https://users-phys.au.dk/~jcd/oscilnotes/Lecture_Notes_on_Stellar_Oscillations.pdf}}</ref>

If a gas parcel is pushed up and <math>N^2>0</math>, the air parcel will move up and down around the height where the density of the parcel matches the density of the surrounding air. If the air parcel is pushed up and <math>N^2=0</math>, the air parcel will not move any further. If the air parcel is pushed up and <math>N^2<0</math>, (i.e. the Brunt–Väisälä frequency is imaginary), then the air parcel will rise and rise unless <math>N^2</math> becomes positive or zero again further up in the atmosphere. In practice this leads to convection, and hence the [[Schwarzschild criterion]]  for stability against convection (or the [[Paul Ledoux#Ledoux criterion|Ledoux criterion]] if there is compositional stratification) is equivalent to the statement that <math>N^2</math> should be positive.

The Brunt–Väisälä frequency commonly appears in the thermodynamic equations for the atmosphere and in the structure of stars.

[[File:Brunt-Vaisala N-1.png|thumb|Trajectory of a parcel of fluid displaced by 1m in an unstably stratified fluid of Brunt-Väisälä frequency N² = -1 /s²]]
[[File:Brunt-Vaisala N1.png|thumb|Oscillations of a parcel of fluid initially displaced by 1m in a stably stratified fluid with Brunt-Väisälä frequency N=1/s.]]



==In oceanography==
In the '''ocean''' where [[salinity]] is important, or in fresh water lakes near freezing, where density is not a linear function of temperature:<math display="block">N\equiv \sqrt{-{g\over{\rho}} {d\rho\over{dz}}}</math>where <math>\rho</math>, the [[potential density]], depends on both temperature and salinity. An example of Brunt–Väisälä oscillation in a density stratified liquid can be observed in the 'Magic Cork' movie [https://www.youtube.com/watch?v=_DqOyvWwqWI here] .

==Context==
The concept derives from Newton's Second Law when applied to a fluid parcel in the presence of a background stratification (in which the density changes in the vertical - i.e. the density can be said to have multiple vertical layers). The parcel, perturbed vertically from its starting position, experiences a vertical acceleration. If the acceleration is back towards the initial position, the stratification is said to be stable and the parcel oscillates vertically. In this case, {{math|<var>N</var><sup>2</sup> > 0}} and the [[angular frequency]] of oscillation is given {{math|<var>N</var>}}. If the acceleration is away from the initial position ({{math|<var>N</var><sup>2</sup> < 0}}), the stratification is unstable. In this case, overturning or convection generally ensues.

The Brunt–Väisälä frequency relates to [[internal gravity waves]]: it is the frequency when the waves propagate horizontally; and it provides a useful description of atmospheric and oceanic stability.

==See also==
*[[Buoyancy]]
*[[Bénard cell]]

== References ==
{{DEFAULTSORT:Brunt-Vaisala Frequency}}
[[Category:Atmospheric thermodynamics]]
[[Category:Atmospheric dynamics]]
[[Category:Fluid dynamics]]
[[Category:Oceanography]]
[[Category:Buoyancy]]
<references />