27.9 드라이든(Dryden) 난류 모델의 이론적 기초

1. 드라이든 모델의 개요와 역사적 배경

드라이든 난류 모델(Dryden turbulence model)은 비행체의 난류 응답 해석을 위한 고전적 통계 모델이다. 이 모델은 Hugh L. Dryden이 20세기 초 NACA에서 수행한 난류 연구를 바탕으로 체계화되었으며, Dryden과 Kuethe의 “The Measurement of Fluctuations of Air Speed by the Hot-Wire Anemometer”(NACA Report No. 320, 1929) 이후의 연구에서 기본 원리가 확립되었다. 이 모델은 비행체 인증과 조종 특성 평가에 표준으로 사용되며, 군용 표준 MIL-F-8785C와 MIL-HDBK-1797에 공식 정의되어 있다.

2. 모델의 기본 가정

드라이든 모델은 다음 기본 가정에 기반한다. 첫째, 난류가 정상 확률 과정(stationary stochastic process)이다. 통계적 특성이 시간에 따라 변화하지 않는다. 둘째, 난류가 가우시안(Gaussian) 분포를 따른다. 이는 중심 극한 정리에 따른 근사이다. 셋째, 난류가 등방성 및 균질성(isotropy and homogeneity)을 가진다. 방향성 및 공간적 차이가 없다고 가정된다. 넷째, 비행체가 이동 중인 경우 Taylor의 고립 가설(frozen turbulence hypothesis)이 성립한다. 이 가정들은 실제 대기 난류의 근사이나, 공학적 해석에 충분한 정확도를 제공한다.

3. 드라이든 자기 상관 함수

드라이든 모델에서 난류 변동의 자기 상관 함수(autocorrelation function)는 다음과 같이 표현된다. 종방향 속도 변동 $u'$ 의 자기 상관은

$R_{uu}(\xi) = \sigma_u^2 \exp\left(-\frac{|\xi|}{L_u}\right)$

여기서 $\sigma_u$ 는 종방향 변동의 표준편차, $L_u$ 는 종방향 적분 길이 스케일, $\xi$ 는 공간 이동 거리이다. 횡방향 및 수직 방향 속도 변동은 약간 다른 형태의 자기 상관을 가진다.

$R_{vv}(\xi) = \sigma_v^2 \left(1 - \frac{|\xi|}{2L_v}\right)\exp\left(-\frac{|\xi|}{L_v}\right)$

$R_{ww}(\xi) = \sigma_w^2 \left(1 - \frac{|\xi|}{2L_w}\right)\exp\left(-\frac{|\xi|}{L_w}\right)$

이러한 함수 형태는 난류 공간 상관을 지수 감쇠 형태로 근사하는 단순하면서도 유용한 수학적 구조이다.

드라이든 파워 스펙트럼 밀도

자기 상관 함수의 푸리에 변환으로부터 드라이든 파워 스펙트럼 밀도(Power Spectral Density, PSD)가 유도된다. 종방향 속도 성분의 PSD는 다음과 같다.

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

여기서 $\omega$ 는 각주파수, $U$ 는 비행체 속도이다. 횡방향 성분은

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

수직 방향 성분 $\Phi_w(\omega)$ 도 유사한 형태를 가진다. 이러한 PSD는 비행체의 공력 응답 해석과 제어 시스템 설계에 직접 사용된다.

드라이든 형성 필터

시간 영역 시뮬레이션에서는 드라이든 PSD를 재현하는 형성 필터(forming filter)가 사용된다. 백색 잡음(white noise)을 입력으로 하여 드라이든 PSD를 갖는 출력을 생성하는 전달 함수는 유리 함수(rational function) 형태로 표현된다. 종방향 성분의 경우

$H_u(s) = \sigma_u \sqrt{\frac{2 L_u}{\pi U}} \cdot \frac{1}{1 + (L_u/U) s}$

여기서 $s$ 는 라플라스 변수이다. 이 필터는 1차 저역 통과 필터 형태로서 구현이 단순하다. 횡방향 및 수직 성분은 2차 필터로 표현된다.

4. 난류 강도 파라미터

드라이든 모델의 난류 강도 파라미터 $\sigma_u$ , $\sigma_v$ , $\sigma_w$ 는 대기 조건에 따라 달라진다. 등방성 가정 하에서는 세 성분이 동일하지만, 실제 대기에서는 이방성을 가진다. 군용 표준 MIL-HDBK-1797은 고도에 따른 난류 강도를 제시한다. 저고도(750 ft 이하)에서

$\sigma_w = 0.1 W_{20}$

$\frac{\sigma_u}{\sigma_w} = \frac{\sigma_v}{\sigma_w} = \frac{1}{(0.177 + 0.000823 h)^{0.4}}$

여기서 $W_{20}$ 은 지상 20 ft 높이에서의 풍속, $h$ 는 고도(ft)이다. 이러한 관계는 실험적 측정을 반영한다.

5. 길이 스케일 파라미터

난류 적분 길이 스케일은 난류 구조의 특성 크기를 나타낸다. 드라이든 모델의 길이 스케일은 고도와 대기 조건에 따라 지정된다. 저고도에서

$L_w = h$

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

이 관계는 경계층 내부에서 와류 크기가 고도에 제약된다는 물리적 통찰을 반영한다. 고고도에서는 $L_u = L_v = L_w = 1750 \, \text{ft}$ 로 일정하게 설정된다.

6. 등방성 가정의 한계

드라이든 모델은 등방성 난류를 가정하나 실제 대기 경계층 난류는 이방성을 가진다. 일반적으로 종방향 변동이 크고, 수직 방향 변동이 작다. 모델에서는 이러한 이방성을 파라미터 조정을 통해 부분적으로 반영한다. 또한 대기 안정도에 따른 이방성 변화도 표준 표(lookup table)를 통해 반영된다. 그럼에도 불구하고 모델은 완전한 이방성 재현에 한계가 있다.

7. 관성 부영역에서의 편차

드라이든 모델의 PSD는 $\omega \to \infty$ 에서 $\omega^{-2}$ 로 감쇠하나, Kolmogorov 이론에 따른 실제 관성 부영역의 스펙트럼은 $\omega^{-5/3}$ 로 감쇠한다. 따라서 드라이든 모델은 고주파 영역에서 실제보다 에너지를 과소 예측한다. 그러나 이러한 편차는 공학적 응용에서는 대체로 무시할 수 있는 수준이다. 더 정확한 스펙트럼 형태가 필요한 경우 폰 카르만 모델이 사용된다.

8. 비행체의 공력 응답

드라이든 모델의 PSD를 입력으로 하여 비행체의 공력 응답을 해석할 수 있다. 비행체의 전달 함수 $G(s)$ 가 알려진 경우 출력의 PSD는 다음과 같다.

$\Phi_{out}(\omega) = |G(i\omega)|^2 \Phi_{in}(\omega)$

이 관계는 난류가 비행체 반응(자세, 가속도, 양력 등)에 미치는 영향을 주파수 영역에서 해석하는 기본 도구이다. 이를 통해 난류에 대한 비행체의 감도, 응답 스펙트럼, 그리고 과도 하중 등을 평가할 수 있다.

드라이든 모델의 활용

드라이든 모델은 항공 산업 및 관련 분야에서 다음과 같이 활용된다. 첫째, 새로운 항공기의 난류 응답 예측 및 설계 평가에 사용된다. 둘째, 조종 특성 평가 및 인증 시험의 기준으로 사용된다. 셋째, 비행 시뮬레이터에서 난류 효과를 재현하는 데 사용된다. 넷째, 제어 시스템의 외란 제거 성능 검증에 사용된다. 드론 및 UAM 기체 설계에서도 이 모델이 필수적으로 활용된다.

시뮬레이션 구현

드라이든 모델의 시뮬레이션 구현은 다음 단계로 수행된다. 첫째, 현재 고도와 대기 조건에 맞는 난류 강도와 길이 스케일을 결정한다. 둘째, 형성 필터를 백색 잡음으로 여기시켜 난류 속도 시계열을 생성한다. 셋째, 생성된 시계열을 비행체 공력 모델에 외란으로 입력한다. 넷째, 비행체의 응답을 시간 영역에서 적분한다. MATLAB/Simulink의 Aerospace Toolbox, X-Plane, FlightGear 등의 상용 및 오픈소스 시뮬레이터는 이러한 구현을 제공한다.

모델의 제한 사항과 대안

드라이든 모델은 공학적 응용에 유용하나 다음 제한이 있다. 첫째, 고주파 영역에서의 스펙트럼 형태 오차. 둘째, 이방성의 제한적 반영. 셋째, 비가우시안 극단 현상의 미반영. 넷째, 공간적 상관 구조의 단순화. 이러한 제한을 극복하기 위해 폰 카르만 모델, 수정 드라이든 모델, 그리고 물리 기반 고해상도 모델 등이 대안으로 사용된다. 그럼에도 불구하고 드라이든 모델은 단순성과 표준화로 인해 여전히 널리 사용되고 있다.

출처

Dryden, H. L., and Kuethe, A. M., “The Measurement of Fluctuations of Air Speed by the Hot-Wire Anemometer,” NACA Report No. 320, 1929.
U.S. Department of Defense, “Flying Qualities of Piloted Airplanes,” MIL-F-8785C, 1980.
U.S. Department of Defense, “Flying Qualities of Piloted Aircraft,” MIL-HDBK-1797, 1997.
Etkin, B., “Dynamics of Atmospheric Flight,” Dover Publications, 2005.
Hoblit, F. M., “Gust Loads on Aircraft: Concepts and Applications,” AIAA Education Series, 1988.

버전

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