27.12 폰 카르만 모델의 파워 스펙트럼 밀도 함수

1. 폰 카르만 PSD의 역할

폰 카르만 모델의 파워 스펙트럼 밀도(Power Spectral Density, PSD) 함수는 난류 속도 변동의 주파수별 에너지 분포를 기술하는 핵심 함수이다. 이 PSD는 비행체의 난류 응답 해석, 제어 시스템 설계, 비행 시뮬레이션에서 핵심적 역할을 한다. 본 절에서는 폰 카르만 PSD의 수학적 형태, 각 속도 성분별 특성, 그리고 실용적 구현을 상세히 서술한다.

2. 종방향 속도 PSD

비행체 속도 $U$ 기준의 종방향 속도 변동 $u_g$ 의 PSD는 다음과 같이 표현된다.

$\Phi_u(\Omega) = \frac{2 \sigma_u^2 L_u}{\pi} \cdot \frac{1}{[1 + (1.339 L_u \Omega)^2]^{5/6}}$

여기서 $\Omega = \omega / U$ 는 공간 파수이고, $\omega$ 는 시간 각주파수, $L_u$ 는 종방향 적분 길이 스케일, $\sigma_u$ 는 종방향 난류 분산의 제곱근이다. 상수 $1.339$ 는 폰 카르만 모델의 특징적 수치이며, 적분 길이 스케일의 정의와 관련된다.

주파수 기반 표현

실제 응용에서는 주파수 $f = \omega/(2\pi)$ 또는 각주파수 $\omega$ 기반의 PSD가 더 자주 사용된다. 주파수 기반 폰 카르만 PSD는 다음과 같다.

$\Phi_u(\omega) = \sigma_u^2 \frac{2 L_u / U}{\pi} \cdot \frac{1}{[1 + (1.339 L_u \omega / U)^2]^{5/6}}$

이 식은 비행체 고정 좌표계에서 측정되는 난류 속도 변동의 주파수 스펙트럼을 기술한다. 총 분산은 $\int_0^\infty \Phi_u(\omega) d\omega = \sigma_u^2$ 의 관계를 만족한다.

3. 횡방향 속도 PSD

횡방향 속도 변동 $v_g$ 의 PSD는 종방향과 다른 형태를 가진다.

$\Phi_v(\omega) = \sigma_v^2 \frac{L_v / U}{\pi} \cdot \frac{1 + (8/3)(1.339 L_v \omega / U)^2}{[1 + (1.339 L_v \omega / U)^2]^{11/6}}$

분자의 추가 항은 횡방향 난류의 공간 상관 구조에서 비롯된다. 등방성 난류에서 $L_v = L_u / 2$ 의 관계가 성립하므로, 실제 적용에서는 $L_u$ 만 지정하고 $L_v$ 를 계산하여 사용한다.

수직 방향 속도 PSD

수직 방향 속도 변동 $w_g$ 의 PSD는 횡방향과 동일한 함수 형태를 가진다.

$\Phi_w(\omega) = \sigma_w^2 \frac{L_w / U}{\pi} \cdot \frac{1 + (8/3)(1.339 L_w \omega / U)^2}{[1 + (1.339 L_w \omega / U)^2]^{11/6}}$

등방성 난류에서 $L_w = L_u / 2$ 이고 $\sigma_w = \sigma_u$ 이다. 실제 대기에서는 이방성으로 인해 이 관계가 정확히 성립하지 않으며, 표준 규정에서 고도별 값이 제공된다.

4. PSD의 특성 주파수

각 PSD는 특성 주파수 $\omega_c = U / L$ 로 정의되는 전이점을 가진다. $\omega < \omega_c$ 에서 PSD는 거의 일정한 값(저주파 플랫 영역)이며, $\omega > \omega_c$ 에서 고주파 감쇠가 시작된다. 종방향 PSD의 경우 고주파 영역에서 $\Phi_u \propto \omega^{-5/3}$ 로 감쇠하며, 이는 Kolmogorov의 관성 부영역 법칙과 일치한다.

5. 파라미터의 고도 의존성

폰 카르만 모델의 적분 길이 스케일과 난류 강도는 고도에 따라 변화한다. MIL-HDBK-1797은 다음과 같은 규정을 제공한다. 저고도(10 ft 초과 1000 ft 이하)에서

$L_w = h, \quad L_u = 2 L_v = \frac{h}{(0.177 + 0.000823 h)^{1.2}}$

고고도(2000 ft 이상)에서

$L_u = 2 L_v = 2 L_w = 1750 \, \text{ft}$

이러한 고도 의존성은 대기 경계층 내에서 와류 크기가 지표로부터의 거리에 제약된다는 물리적 관찰을 반영한다.

6. 관성 부영역의 검증

폰 카르만 PSD가 관성 부영역에서 Kolmogorov 법칙을 만족하는지 확인하기 위해 고주파 극한을 살펴본다. $L_u \omega / U \gg 1$ 에서 종방향 PSD는

$\Phi_u(\omega) \approx \sigma_u^2 \frac{2 L_u / U}{\pi} \cdot (1.339 L_u \omega / U)^{-5/3}$

$= \frac{2 \sigma_u^2}{\pi (1.339)^{5/3}} \cdot \frac{L_u^{-2/3}}{U^{-2/3}} \cdot \omega^{-5/3}$

Kolmogorov의 관성 부영역 PSD $\Phi_u \propto \varepsilon^{2/3} (\omega/U)^{-5/3}$ 와 비교하면, 폰 카르만 모델이 $-5/3$ 의존성을 정확히 재현함을 확인할 수 있다. 난류 에너지 소산율 $\varepsilon$ 과 모델 파라미터( $\sigma_u$ , $L_u$ ) 사이의 관계도 이로부터 유도된다.

7. 수치 계산의 실무적 고려

폰 카르만 PSD를 수치적으로 계산할 때는 다음 사항을 고려한다. 첫째, 주파수 범위는 저주파 플랫 영역과 고주파 감쇠 영역을 모두 포함해야 한다. 일반적으로 0.01 Hz부터 100 Hz까지의 범위가 사용된다. 둘째, 로그 주파수 간격이 자연스럽다. 셋째, 수치 적분 시 고주파 꼬리에서의 적분 오차를 제어하기 위해 적절한 상한 주파수 설정이 필요하다. 넷째, 분산 정규화를 확인하여 모델의 일관성을 검증한다.

8. 유리 함수 근사

폰 카르만 PSD는 $[...]^{5/6}$ 또는 $[...]^{11/6}$ 지수를 포함하여 유리 함수로 정확히 표현되지 않는다. 시간 영역 시뮬레이션을 위해서는 유리 함수 근사가 필요하다. 대표적 근사는 Reeves의 제안이며, 낮은 차수의 유리 함수로 폰 카르만 PSD를 근사한다. 예시로 종방향 PSD의 2차 근사는

$H_u(s) \approx \sigma_u \sqrt{\frac{2 L_u}{\pi U}} \cdot \frac{1 + 0.25 (L_u/U) s}{1 + 1.357 (L_u/U) s + 0.1987 ((L_u/U) s)^2}$

이러한 근사는 관심 주파수 범위에서 폰 카르만 PSD를 잘 재현한다.

PSD의 실험적 검증

폰 카르만 PSD는 수많은 실측으로 검증되었다. 초음파 풍속계로 측정한 지표층 풍속 변동의 스펙트럼은 관성 부영역에서 Kolmogorov 법칙을 정확히 보이며, 에너지 생성 영역에서도 폰 카르만 모델의 예측과 잘 일치한다. 대류권 상층과 성층권에서도 유사한 결과가 얻어진다. 이러한 광범위한 실증적 지지는 폰 카르만 모델의 학술적 권위를 뒷받침한다.

공간 PSD와 시간 PSD의 변환

폰 카르만 모델은 공간 PSD 형태로 유도되며, 비행체의 이동에 의해 시간 PSD로 변환된다. Taylor의 고립 가설에 의해 비행체 속도 $U$ 로 이동할 때 공간 파수 $k$ 와 시간 각주파수 $\omega$ 는 $\omega = U k$ 의 관계를 가진다. PSD의 변환은

$\Phi(\omega) = \frac{E(k)}{U}\bigg|_{k = \omega/U}$

이러한 변환은 정지 비행체에서의 측정 스펙트럼과 이동 비행체에서의 스펙트럼을 연관시킨다.

9. 드론의 응용에서의 스펙트럼 사용

드론의 난류 응답 해석에서 폰 카르만 PSD는 다음과 같이 사용된다. 첫째, 난류 속도 변동이 드론의 자세 및 위치에 미치는 주파수별 영향을 정량화한다. 둘째, 제어 시스템의 주파수 응답 특성과 결합하여 최종 출력의 주파수 스펙트럼을 예측한다. 셋째, 제어기 설계에서 난류 외란 감쇠 성능을 평가한다. 넷째, 센서 노이즈와 난류 외란의 분리 설계에 활용된다. 이러한 응용에서 정확한 PSD의 사용은 드론 성능 평가와 개선에 필수적이다.

10. 출처

von Kármán, T., “Progress in the Statistical Theory of Turbulence,” Proceedings of the National Academy of Sciences, Vol. 34, No. 11, 1948.
U.S. Department of Defense, “Flying Qualities of Piloted Aircraft,” MIL-HDBK-1797, 1997.
Hoblit, F. M., “Gust Loads on Aircraft: Concepts and Applications,” AIAA Education Series, 1988.
Etkin, B., “Dynamics of Atmospheric Flight,” Dover Publications, 2005.
Reeves, P. M., Joshi, V. M., and Campbell, G. S., “Development of a Non-Gaussian Atmospheric Turbulence Model for Use in Flight Simulators,” NASA CR-2237, 1974.

11. 버전

문서 버전: v1.0
작성 기준일: 2026-04-17