21.30 돌풍과 난류 모델링 (Gust and Turbulence Modeling)

1. 대기 난류의 본질과 공기역학적 의미

대기 난류(atmospheric turbulence)는 대기 중에서 발생하는 불규칙하고 비정상적인 유동 변동으로서, 비행체에 작용하는 공기역학적 힘과 모멘트를 시간적으로 변동시키는 주요 외란(disturbance)이다. 대기 난류는 태양 복사에 의한 열적 대류(thermal convection), 지형에 의한 기계적 난류(mechanical turbulence), 풍속 전단(wind shear), 대류 활동(convective activity) 등 다양한 기상학적 원인에 의해 발생한다.

비행 로봇의 관점에서 대기 난류는 받음각과 옆미끄럼각의 순간적 변동, 비행 속도의 변화, 비대칭 공력 하중의 발생 등을 초래하며, 이는 비행 경로의 이탈, 구조적 피로 하중, 탑재 센서의 교란, 최악의 경우 비행체의 제어 상실로 이어질 수 있다. 따라서 대기 난류의 정확한 모델링은 비행체의 구조 설계, 비행 제어 시스템 설계, 비행 시뮬레이션에서 필수적이다.

2. 난류의 통계적 기술

대기 난류는 확률적(stochastic) 현상으로서, 결정론적(deterministic) 예측이 불가능하며 통계적 방법으로 기술된다. 난류의 통계적 특성을 기술하는 주요 매개변수는 다음과 같다:

난류 강도(turbulence intensity, \sigma_u, \sigma_v, \sigma_w): 각 방향 속도 변동의 표준 편차로서, 난류의 세기를 나타낸다. 수직 방향 난류 강도 \sigma_w가 비행체의 받음각 변동에 가장 직접적으로 영향을 미친다.

자기 상관 함수(autocorrelation function): 난류 속도 변동의 시간적 상관 구조를 나타낸다. 자기 상관 함수의 적분 스케일(integral scale)은 난류의 특성 시간 스케일에 해당한다.

파워 스펙트럼 밀도(Power Spectral Density, PSD): 난류 에너지의 주파수 분포를 나타내며, 난류 모델의 핵심 특성이다. 콜모고로프의 관성 부영역(inertial subrange) 이론에 의하면, 등방성 난류의 PSD는 고주파 영역에서 -5/3 지수 법칙을 따른다.

난류 길이 스케일(turbulence length scale, L): 난류 와류의 특성 공간 스케일로서, 난류의 공간적 상관 구조를 특성화한다.

3. 드라이든 난류 모델

드라이든 난류 모델(Dryden turbulence model)은 대기 난류를 유리(rational) 전달 함수 형태의 파워 스펙트럼으로 기술하는 모델이다. 드라이든 모델의 종방향(u) 파워 스펙트럼 밀도는 다음과 같다:

\Phi_u(\Omega) = \sigma_u^2 \frac{2L_u}{\pi} \frac{1}{1 + (L_u \Omega)^2}

여기서 \Omega = \omega/V는 공간 주파수(spatial frequency), \omega는 각주파수, V는 비행 속도, L_u는 종방향 난류 길이 스케일이다.

수직 방향(w)의 드라이든 스펙트럼은:

\Phi_w(\Omega) = \sigma_w^2 \frac{L_w}{\pi} \frac{1 + 3(L_w \Omega)^2}{[1 + (L_w \Omega)^2]^2}

드라이든 모델의 장점은 유리 전달 함수로 표현되므로, 백색 잡음(white noise)을 입력으로 하는 성형 필터(shaping filter)로 시간 영역에서 효율적으로 구현할 수 있다는 것이다. 이는 비행 시뮬레이션과 비행 제어 시스템의 설계에 직접 활용된다. 미국 군사 규격 MIL-HDBK-1797에 드라이든 모델의 매개변수가 고도별로 규정되어 있다 (MIL-HDBK-1797, 1997).

4. 폰 카르만 난류 모델

폰 카르만 난류 모델(von Kármán turbulence model)은 드라이든 모델보다 고주파 영역에서의 난류 스펙트럼을 더 정확하게 기술하는 모델이다. 종방향 파워 스펙트럼 밀도는:

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

수직 방향의 폰 카르만 스펙트럼은:

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

폰 카르만 모델은 고주파 영역에서 -5/3 지수 법칙을 더 정확하게 재현하며, 등방성 난류의 이론적 스펙트럼과 더 잘 일치한다. 그러나 비유리(irrational) 전달 함수 형태이므로, 시간 영역 시뮬레이션에서의 구현에는 유리 함수 근사가 필요하다. 국제 규격 ISO 2533과 ESDU에서 폰 카르만 모델의 사용을 권고하고 있다.

5. 이산 돌풍 모델

이산 돌풍(discrete gust)은 단일의 결정론적 속도 변화로 모델링되며, 비행체의 과도 응답(transient response) 해석에 사용된다. 대표적인 이산 돌풍 모델은 다음과 같다:

1-코사인 돌풍(1-minus-cosine gust): 가장 널리 사용되는 이산 돌풍 모델로서, 돌풍 속도가 시간(또는 거리)에 따라 다음과 같이 변화한다:

w_g(x) = \frac{w_0}{2}\left[1 - \cos\left(\frac{2\pi x}{2H}\right)\right], \quad 0 \leq x \leq 2H

여기서 w_0는 돌풍의 최대 속도, H는 돌풍의 구배 거리(gradient distance)이다. 이 모델은 항공기 인증 규격(FAR/CS 25.341)에서 구조적 돌풍 하중의 산정에 사용된다.

날카로운 돌풍(sharp-edged gust): 돌풍 속도가 순간적으로 최대값에 도달하는 이상화된 모델로서, 퀴슈너 함수(Küssner function)에 의한 비정상 공력 응답의 해석에 사용된다.

21.30.6 돌풍에 대한 비정상 공력 응답

비행체가 돌풍을 조우하면 받음각이 순간적으로 변화하지만, 양력의 발생에는 시간 지연이 존재한다. 이러한 비정상 공력 응답은 퀴슈너 함수(Küssner function, \psi)에 의해 기술된다. 퀴슈너 함수는 날카로운 돌풍에 대한 양력의 시간적 성장을 나타내는 무차원 함수이다.

무차원 시간(반시위 이동 수, s = 2Vt/c)에 대한 퀴슈너 함수의 근사식은:

\psi(s) \approx 1 - 0.5e^{-0.13s} - 0.5e^{-s}

이 함수는 돌풍 진입 직후 양력이 정상 값의 약 50%에서 시작하여, 시간이 경과함에 따라 점근적으로 정상 값에 도달함을 나타낸다. 양력의 시간 지연은 후류 와류의 발달과 유동의 비정상적 조정(unsteady adjustment)에 기인한다.

임의의 돌풍 프로파일 w_g(t)에 대한 비정상 양력은 뒤하멜 적분(Duhamel integral)으로 계산된다:

\Delta C_L(t) = \frac{2\pi}{V}\int_0^t \frac{dw_g(\tau)}{d\tau}\psi(t - \tau) \, d\tau

21.30.7 대기 경계층의 난류 구조

지표면 부근의 대기 경계층(Atmospheric Boundary Layer, ABL)은 비행 로봇의 주요 운용 고도에 해당하며, 이 영역의 난류 구조는 비행 성능과 안전성에 직접적 영향을 미친다.

대기 경계층의 평균 풍속 분포는 고도에 따라 증가하며, 대수 법칙(logarithmic law)에 의해 기술된다:

\bar{u}(z) = \frac{u_*}{\kappa}\ln\left(\frac{z}{z_0}\right)

여기서 u_*는 마찰 속도(friction velocity), \kappa \approx 0.41은 카르만 상수, z_0는 표면 거칠기 길이(surface roughness length)이다. 난류 강도는 일반적으로 지표면에 가까울수록 높으며, 고도의 증가에 따라 감소한다.

소형 비행 로봇이 운용되는 저고도(지상 10~150 m)에서의 난류 특성은 고고도 항공기의 운용 환경과 현저히 다르며, 지형과 구조물에 의한 기계적 난류, 열적 대류, 풍속 전단 등이 복합적으로 작용한다 (Watkins et al., 2006).

6. 난류 시뮬레이션 기법

비행 시뮬레이션에서 난류 시간 이력(time history)을 생성하기 위한 주요 기법은 다음과 같다:

성형 필터법(shaping filter method): 백색 잡음을 드라이든 또는 폰 카르만 모델의 전달 함수를 통해 통과시켜 난류 시간 이력을 생성한다. 실시간 시뮬레이션에 적합하며, 드라이든 모델은 직접 구현이 가능하고, 폰 카르만 모델은 유리 함수 근사 후 구현한다.

공간 난류장 생성(spatial turbulence field generation): 다수의 비행체가 동일 난류장을 경험하는 편대 비행 시뮬레이션이나, 비행체의 스팬 방향으로의 난류 공간 상관을 고려하는 경우에는 공간적으로 상관된 난류장의 생성이 필요하다. 스펙트럼 방법(spectral method)이나 회전 합성법(vortex particle method)이 사용된다.

대와류 시뮬레이션(LES): 대기 경계층의 난류를 직접 수치적으로 해석하여 난류 유동장을 생성한다. 가장 물리적으로 정확하지만 계산 비용이 높다.

7. 소형 비행 로봇에 대한 난류의 상대적 영향

소형 비행 로봇은 크기와 질량이 작으므로, 동일한 대기 난류에 대한 응답이 대형 항공기보다 현저히 크다. 이는 다음의 물리적 요인에 기인한다:

크기 대 난류 스케일의 비: 소형 비행 로봇의 크기(스팬 0.5~3 m)가 대기 난류의 에너지 함유 와류(energy-containing eddy)의 스케일(수 미터~수백 미터)보다 작으므로, 비행체는 난류 와류에 의해 강체적으로 운반(bodily transport)된다. 대형 항공기에서는 비행체 크기가 난류 스케일과 유사하거나 더 커서 공간적 평균화(spatial averaging)에 의해 난류 응답이 완화된다.

낮은 관성: 소형 비행 로봇의 낮은 질량과 관성 모멘트는 외란에 대한 각속도 응답을 크게 만든다.

낮은 날개 하중: 낮은 날개 하중(W/S)은 돌풍에 대한 하중 배수 변동을 증가시킨다.

이러한 이유로 소형 비행 로봇의 설계에서 대기 난류에 대한 강건성(robustness)은 핵심적 설계 요구 사항이다 (Watkins et al., 2006).

8. 로봇 공학에서의 돌풍 경감 기술

비행 로봇의 난류 환경에서의 안정적 운용을 위한 돌풍 경감(gust alleviation) 기술은 활발히 연구되고 있다:

수동적 돌풍 경감: 날개의 탄성 변형(aeroelastic deformation)에 의한 하중 완화, 날개 비틀림에 의한 유효 받음각 변동의 감소 등 구조적·공기역학적 수동 메커니즘이다.

능동적 돌풍 경감: 센서(받음각 센서, 가속도계, 피토관 등)에 의해 돌풍을 감지하고, 제어면(에일러론, 엘리베이터)이나 로터 회전 속도를 조정하여 돌풍에 의한 공력 변동을 보상하는 피드백 제어 시스템이다.

전방 감지(look-ahead sensing): 라이다(LiDAR)나 초음파 풍속계를 이용하여 비행체 전방의 돌풍을 사전에 감지하고, 피드포워드(feedforward) 제어를 통해 돌풍 도달 전에 선제적으로 대응하는 기술이다. 이 접근법은 피드백 제어의 시간 지연 문제를 근본적으로 극복할 수 있는 가능성을 제공한다 (Mohamed et al., 2014).

참고 문헌

• Etkin, B. (1972). Dynamics of Atmospheric Flight. John Wiley & Sons.
• MIL-HDBK-1797 (1997). Flying Qualities of Piloted Aircraft. U.S. Department of Defense.
• Mohamed, A., Massey, K., Watkins, S., & Clothier, R. (2014). The attitude control of fixed-wing MAVs in turbulent environments. Progress in Aerospace Sciences, 66, 37–48.
• Watkins, S., Milbank, J., Loxton, B. J., & Melbourne, W. H. (2006). Atmospheric winds and their implications for microair vehicles. AIAA Journal, 44(11), 2591–2600.

버전: 2026-04-16-v1.0