Misplaced Pages

von Kármán wind turbulence model

Article snapshot taken from Wikipedia with creative commons attribution-sharealike license. Give it a read and then ask your questions in the chat. We can research this topic together.
(Redirected from Von Kármán Wind Turbulence Model) Mathematical model of continuous gusts

The von Kármán wind turbulence model (also known as von Kármán gusts) is a mathematical model of continuous gusts. It matches observed continuous gusts better than that Dryden Wind Turbulence Model and is the preferred model of the United States Department of Defense in most aircraft design and simulation applications. The von Kármán model treats the linear and angular velocity components of continuous gusts as spatially varying stochastic processes and specifies each component's power spectral density. The von Kármán wind turbulence model is characterized by irrational power spectral densities, so filters can be designed that take white noise inputs and output stochastic processes with the approximated von Kármán gusts' power spectral densities.

History

The von Kármán wind turbulence model first appeared in a 1957 NACA report based on earlier work by Theodore von Kármán.

Further information: Sandia method

Power spectral densities

The von Kármán model is characterized by single-sided power spectral densities for gusts' three linear velocity components (ug, vg, and wg),

Φ u g ( Ω ) = σ u 2 2 L u π 1 ( 1 + ( 1.339 L u Ω ) 2 ) 5 6 Φ v g ( Ω ) = σ v 2 2 L v π 1 + 8 3 ( 2.678 L v Ω ) 2 ( 1 + ( 2.678 L v Ω ) 2 ) 11 6 Φ w g ( Ω ) = σ w 2 2 L w π 1 + 8 3 ( 2.678 L w Ω ) 2 ( 1 + ( 2.678 L w Ω ) 2 ) 11 6 {\displaystyle {\begin{aligned}\Phi _{u_{g}}(\Omega )&=\sigma _{u}^{2}{\frac {2L_{u}}{\pi }}{\frac {1}{\left(1+(1.339L_{u}\Omega )^{2}\right)^{\frac {5}{6}}}}\\\Phi _{v_{g}}(\Omega )&=\sigma _{v}^{2}{\frac {2L_{v}}{\pi }}{\frac {1+{\frac {8}{3}}(2.678L_{v}\Omega )^{2}}{\left(1+(2.678L_{v}\Omega )^{2}\right)^{\frac {11}{6}}}}\\\Phi _{w_{g}}(\Omega )&=\sigma _{w}^{2}{\frac {2L_{w}}{\pi }}{\frac {1+{\frac {8}{3}}(2.678L_{w}\Omega )^{2}}{\left(1+(2.678L_{w}\Omega )^{2}\right)^{\frac {11}{6}}}}\end{aligned}}}

where σi and Li are the turbulence intensity and scale length, respectively, for the ith velocity component, and Ω is a spatial frequency. These power spectral densities give the stochastic process spatial variations, but any temporal variations rely on vehicle motion through the gust velocity field. The speed with which the vehicle is moving through the gust field V allows conversion of these power spectral densities to different types of frequencies,

Ω = ω V Φ i ( Ω ) = V Φ i ( ω ) {\displaystyle {\begin{aligned}\Omega &={\frac {\omega }{V}}\\\Phi _{i}(\Omega )&=V\Phi _{i}\left(\omega \right)\end{aligned}}}

where ω has units of radians per unit time.

The gust angular velocity components (pg, qg, rg) are defined as the variations of the linear velocity components along the different vehicle axes,

p g = ∂ w g ∂ y q g = ∂ w g ∂ x r g = − ∂ v g ∂ x {\displaystyle {\begin{aligned}p_{g}&={\frac {\partial w_{g}}{\partial y}}\\q_{g}&={\frac {\partial w_{g}}{\partial x}}\\r_{g}&=-{\frac {\partial v_{g}}{\partial x}}\end{aligned}}}

though different sign conventions may be used in some sources. The power spectral densities for the angular velocity components are

Φ p g ( ω ) = σ w 2 2 V L w 0.8 ( 2 π L w 4 b ) 1 3 1 + ( 4 b ω π V ) 2 Φ q g ( ω ) = ± ( ω V ) 2 1 + ( 4 b ω π V ) 2 Φ w g ( ω ) Φ r g ( ω ) = ∓ ( ω V ) 2 1 + ( 3 b ω π V ) 2 Φ v g ( ω ) {\displaystyle {\begin{aligned}\Phi _{p_{g}}(\omega )&={\frac {\sigma _{w}^{2}}{2VL_{w}}}{\frac {0.8\left({\frac {2\pi L_{w}}{4b}}\right)^{\frac {1}{3}}}{1+\left({\frac {4b\omega }{\pi V}}\right)^{2}}}\\\Phi _{q_{g}}(\omega )&={\frac {\pm \left({\frac {\omega }{V}}\right)^{2}}{1+\left({\frac {4b\omega }{\pi V}}\right)^{2}}}\Phi _{w_{g}}(\omega )\\\Phi _{r_{g}}(\omega )&={\frac {\mp \left({\frac {\omega }{V}}\right)^{2}}{1+\left({\frac {3b\omega }{\pi V}}\right)^{2}}}\Phi _{v_{g}}(\omega )\end{aligned}}}

The military specifications give criteria based on vehicle stability derivatives to determine whether the gust angular velocity components are significant.

Spectral factorization

The gusts generated by the von Kármán model are not a white noise process and therefore may be referred to as colored noise. Colored noise may, in some circumstances, be generated as the output of a minimum phase linear filter through a process known as spectral factorization. Consider a linear time invariant system with a white noise input that has unit variance, transfer function G(s), and output y(t). The power spectral density of y(t) is

Φ y ( ω ) = | G ( i ω ) | 2 {\displaystyle \Phi _{y}(\omega )=|G(i\omega )|^{2}}

where i = -1. For irrational power spectral densities, such as that of the von Kármán model, a suitable transfer function can be found whose magnitude squared evaluated along the imaginary axis approximates the power spectral density. The MATLAB documentation provides a realization of such a transfer function for von Kármán gusts that is consistent with the military specifications,

G u g ( s ) = σ u 2 L u π V ( 1 + 0.25 L u V s ) 1 + 1.357 L u V s + 0.1987 ( L u V s ) 2 G v g ( s ) = σ v 2 L v π V ( 1 + 2.7478 2 L v V s + 0.3398 ( 2 L v V s ) 2 ) 1 + 2.9958 2 L v V s + 1.9754 ( 2 L v V s ) 2 + 0.1539 ( 2 L v V s ) 3 G w g ( s ) = σ w 2 L w π V ( 1 + 2.7478 2 L w V s + 0.3398 ( 2 L w V s ) 2 ) 1 + 2.9958 2 L w V s + 1.9754 ( 2 L w V s ) 2 + 0.1539 ( 2 L w V s ) 3 G p g ( s ) = σ w 0.8 V ( π 4 b ) 1 6 ( 2 L w ) 1 3 ( 1 + 4 b π V s ) G q g ( s ) = ± s V 1 + 4 b π V s G w g ( s ) G r g ( s ) = ∓ s V 1 + 3 b π V s G v g ( s ) {\displaystyle {\begin{aligned}G_{u_{g}}(s)&={\frac {\sigma _{u}{\sqrt {\frac {2L_{u}}{\pi V}}}\left(1+0.25{\frac {L_{u}}{V}}s\right)}{1+1.357{\frac {L_{u}}{V}}s+0.1987\left({\frac {L_{u}}{V}}s\right)^{2}}}\\G_{v_{g}}(s)&={\frac {\sigma _{v}{\sqrt {\frac {2L_{v}}{\pi V}}}\left(1+2.7478{\frac {2L_{v}}{V}}s+0.3398\left({\frac {2L_{v}}{V}}s\right)^{2}\right)}{1+2.9958{\frac {2L_{v}}{V}}s+1.9754\left({\frac {2L_{v}}{V}}s\right)^{2}+0.1539\left({\frac {2L_{v}}{V}}s\right)^{3}}}\\G_{w_{g}}(s)&={\frac {\sigma _{w}{\sqrt {\frac {2L_{w}}{\pi V}}}\left(1+2.7478{\frac {2L_{w}}{V}}s+0.3398\left({\frac {2L_{w}}{V}}s\right)^{2}\right)}{1+2.9958{\frac {2L_{w}}{V}}s+1.9754\left({\frac {2L_{w}}{V}}s\right)^{2}+0.1539\left({\frac {2L_{w}}{V}}s\right)^{3}}}\\G_{p_{g}}(s)&=\sigma _{w}{\sqrt {\frac {0.8}{V}}}{\frac {\left({\frac {\pi }{4b}}\right)^{\frac {1}{6}}}{(2L_{w})^{\frac {1}{3}}\left(1+{\frac {4b}{\pi V}}s\right)}}\\G_{q_{g}}(s)&={\frac {\pm {\frac {s}{V}}}{1+{\frac {4b}{\pi V}}s}}G_{w_{g}}(s)\\G_{r_{g}}(s)&={\frac {\mp {\frac {s}{V}}}{1+{\frac {3b}{\pi V}}s}}G_{v_{g}}(s)\end{aligned}}}

Driving these filters with independent, unit variance, band-limited white noise yields outputs with power spectral densities that approximate the power spectral densities of the velocity components of the von Kármán model. The outputs can, in turn, be used as wind disturbance inputs for aircraft or other dynamic systems.

Altitude dependence

Further information: Continuous Gusts § Altitude Dependence

The von Kármán model is parameterized by a length scale and turbulence intensity. The combination of these two parameters determine the shape of the power spectral densities and therefore the quality of the model's fit to spectra of observed turbulence. Many combinations of length scale and turbulence intensity give realistic power spectral densities in the desired frequency ranges. The Department of Defense specifications include choices for both parameters, including their dependence on altitude.

See also

Notes

  1. ^ Hoblit 1988, Chap. 4.
  2. ^ MIL-STD-1797A 1990, p. 678.
  3. Diedrich, Franklin W.; Joseph A. Drischler (1957). "Effect of Spanwise Variations in Gust Intensity on the Lift Due to Atmospheric Turbulence": NACA TN 3920. {{cite journal}}: Cite journal requires |journal= (help)
  4. de Kármán, Theodore; Leslie Howarth (1938). "On the Statistical Theory of Isotropic Turbulence". Proceedings of the Royal Society of London. Series A, Mathematical and Physical Sciences. 164 (917): 192–215. Bibcode:1938RSPSA.164..192D. doi:10.1098/rspa.1938.0013.
  5. von Kármán, Theodore (1948). "Progress in the Statistical Theory of Turbulence". Proceedings of the National Academy of Sciences. 34 (11): 530–539. Bibcode:1948PNAS...34..530V. doi:10.1073/pnas.34.11.530. PMC 1079162. PMID 16588830.
  6. von Kármán, T.; Lin, C. C. (1951). "On the Statistical Theory of Isotropic Turbulence". In von Mises, Richard; von Kármán, Theodore (eds.). Advances in Applied Mechanics. Academic Press, Inc. pp. 1–19. ISBN 9780080563800.
  7. Hoblit 1988, p. ***.
  8. ^ "Von Karman Wind Turbulence Model (Continuous)". MATLAB Reference Pages. The MathWorks, Inc. 2010. Retrieved May 24, 2013.
  9. MIL-STD-1797A 1990, p. 680.
  10. Richardson 2013, p. 33.
  11. MIL-STD-1797A 1990, pp. 673, 678–685, 702.

References

Categories: