Scale analysis (mathematics)

From Wikipedia, the free encyclopedia
  (Redirected from Order of magnitude analysis)
Jump to: navigation, search

Scale analysis is a powerful tool used in the mathematical sciences for the simplification of equations with many terms. First the approximate magnitude of individual terms in the equations is determined. Then some negligibly small terms may be ignored.

Example: vertical momentum in synoptic-scale meteorology[edit]

Consider for example the momentum equation of the Navier–Stokes equations in the vertical coordinate direction of the atmosphere

{{\partial w }\over{\partial t }} + u {\frac{\partial w}{\partial x}} + v {\frac{\partial w}{\partial y}} + w {\frac{\partial w}{\partial z}} - {\frac{u^2 + v^2}{R}}= - { { \frac{1}{\varrho}}{\frac{\partial p}{\partial z}}} - g +2{\Omega u \cos \varphi} + \nu \left({\frac{\partial^2 w}{\partial x^2}}+{\frac{\partial^2 w}{\partial y^2}}+{\frac{\partial^2 w}{\partial z^2}}\right),\qquad(1)

where R is Earth radius, Ω is frequency of rotation of the Earth, g is gravitational acceleration, φ is latitude ρ is density of air and ν is kinematic viscosity of air (we can neglect turbulence in free atmosphere).

In synoptic scale we can expect horizontal velocities about U = 101 m.s−1 and vertical about W = 10−2 m.s−1. Horizontal scale is L = 106 m and vertical scale is H = 104 m. Typical time scale is T = L/U = 105 s. Pressure differences in troposphere are ΔP = 104 Pa and density of air ρ = 100 kg·m−3. Other physical properties are approximately:

R = 6.378 × 106 m;
Ω = 7.292 × 10-5 rad·s−1;
ν = 1.46 × 10−5 m2·s−1;
g = 9.81 m·s−2.

Estimates of the different terms in equation (1) can be made using their scales:


\begin{align}
{{\partial w }\over{\partial t }} &\sim \frac{W}{T} \\[1.2ex]
u {\frac{\partial w}{\partial x}} &\sim U\frac{W}{L} &\qquad
v {\frac{\partial w}{\partial y}} &\sim U\frac{W}{L} &\qquad
w {\frac{\partial w}{\partial z}} &\sim W\frac{W}{H} \\[1.2ex]
{\frac{u^2}{R}} &\sim \frac{U^2}{R} &\qquad
{\frac{v^2}{R}} &\sim \frac{U^2}{R} \\[1.2ex]
\frac{1}{\varrho}\frac{\partial p}{\partial z} &\sim \frac{1}{\varrho}\frac{\Delta P}{H} &\qquad
\Omega u \cos \varphi &\sim \Omega U \\[1.2ex]
\nu \frac{\partial^2 w}{\partial x^2} &\sim \nu \frac{W}{L^2} &\qquad
\nu \frac{\partial^2 w}{\partial y^2} &\sim \nu \frac{W}{L^2} &\qquad
\nu \frac{\partial^2 w}{\partial z^2} &\sim \nu \frac{W}{H^2}
\end{align}

Now we can introduce these scales and their values into equation (1):


{\frac{10^{-2}}{10^5}}+10{\frac{10^{-2}}{10^6}}
+10{\frac{10^{-2}}{10^6}}
+10^{-2}{\frac{10^{-2}}{10^4}}
-{\frac{10^2+10^2}{10^6}}



= - {{\frac{1}{1}} {\frac{10^4}{10^4}} } - 10 + 2 \times 10^{-4} \times 10 + 10^{-5} \left({\frac{10^{-2}}{10^{12}}} + {\frac{10^{-2}}{10^{12}}} + {\frac{10^{-2}}{10^{8}}}  \right).
\qquad (2)

We can see that all terms — except the first and second on the right-hand side — are negligibly small. Thus we can simplify the vertical momentum equation to the hydrostatic equilibrium equation:

{ { \frac{1}{\varrho}}{\frac{\partial p}{\partial z}}} = - g. \qquad (3)

See also[edit]

References[edit]

External links[edit]