27.18 주파수 영역 기반 난류 생성 기법

1. 주파수 영역 접근법의 개요

주파수 영역 기반 난류 생성 기법은 난류의 파워 스펙트럼 밀도(PSD)로부터 직접 시간 시계열을 합성하는 접근법이다. 이 방법은 각 주파수 성분의 진폭을 PSD로부터 결정하고 무작위 위상을 부여하여 합성한다. 시간 영역 형성 필터 접근과 비교하여 스펙트럼 형태를 정확히 재현할 수 있다는 장점이 있다. 본 절에서는 주파수 영역 생성 기법의 이론적 기초, 구현 방법, 그리고 실용적 고려사항을 서술한다.

2. 푸리에 급수 표현

난류 속도 시계열은 푸리에 급수로 표현될 수 있다.

$u(t) = \sum_{n=1}^{N} A_n \cos(\omega_n t + \phi_n)$

여기서 $A_n$ 은 $n$ 번째 주파수 성분의 진폭, $\omega_n$ 은 각주파수, $\phi_n$ 은 위상이다. 무작위 위상 가정하에서 이 급수의 PSD는 각 성분의 $A_n^2 / 2$ 의 합으로 표현된다. 주어진 PSD $\Phi(\omega)$ 를 재현하기 위해서는

$A_n = \sqrt{2 \Phi(\omega_n) \Delta \omega}$

로 진폭을 설정한다. $\Delta \omega$ 는 주파수 간격이다.

3. 무작위 위상 가정

주파수 영역 생성의 핵심은 위상 $\phi_n$ 을 $[0, 2\pi]$ 범위에서 균일 분포로 무작위 생성하는 것이다. 이 무작위 위상으로 합성된 시계열은 정상 확률 과정(stationary stochastic process)의 특성을 가진다. 중심 극한 정리에 따라 충분히 많은 $N$ 에 대해 합성 시계열은 가우시안 분포에 수렴한다. 이는 대기 난류의 가우시안 분포 가정을 자연스럽게 만족시킨다.

4. 이산 푸리에 변환 기반 구현

빠른 푸리에 변환(Fast Fourier Transform, FFT)을 이용하면 주파수 영역 생성이 계산 효율적으로 구현된다. 절차는 다음과 같다. 첫째, $N$ 개의 주파수 $\omega_n = n \Delta \omega$ 를 정의한다. 둘째, 각 주파수에 대해 PSD 값을 계산한다. 셋째, 무작위 위상을 생성하고 복소 주파수 스펙트럼 $X(\omega_n) = \sqrt{N \Phi(\omega_n) \Delta \omega} e^{i \phi_n}$ 을 구성한다. 넷째, 역 FFT를 적용하여 시간 시계열을 얻는다. 이 방법은 대규모 시계열 생성에 적합하다.

5. 주파수 범위와 해상도

주파수 영역 생성에서 주파수 범위 $[0, \omega_{max}]$ 와 해상도 $\Delta \omega$ 의 선택이 중요하다. 최대 주파수 $\omega_{max}$ 는 나이퀴스트 주파수( $\omega_{max} = \pi / T_s$ , $T_s$ 는 샘플링 주기)이며, 관심 주파수 대역을 모두 포함해야 한다. 최소 주파수는 시뮬레이션 총 시간의 역수이며, 긴 시계열을 생성할수록 저주파 해상도가 향상된다. 실용적으로는 최소 주파수가 $1/T_{sim}$ 이 되며, 여기서 $T_{sim}$ 은 시뮬레이션 기간이다.

6. 장점과 단점

주파수 영역 기반 생성의 장점은 다음과 같다. 첫째, PSD 형태를 정확히 재현할 수 있다. 둘째, 임의 PSD에 대해 적용 가능하여 드라이든, 폰 카르만, 또는 사용자 정의 스펙트럼을 모두 처리한다. 셋째, FFT 기반 구현은 $O(N \log N)$ 의 계산 복잡도로 효율적이다. 넷째, 통계적 특성이 명확하게 제어된다.

단점도 존재한다. 첫째, 전체 시계열을 한 번에 생성해야 하므로 실시간 시뮬레이션에 부적합하다. 둘째, 메모리 요구량이 크다. 셋째, 주파수 해상도가 시뮬레이션 총 시간에 의해 제약된다. 넷째, 시계열의 주기성(순환성)이 수치적 아티팩트를 유발할 수 있다.

7. 시간 영역 형성 필터 방법과의 비교

시간 영역 형성 필터 방법은 백색 잡음에 필터를 적용하여 난류 시계열을 점진적으로 생성한다. 이 방법은 실시간 시뮬레이션에 적합하며 메모리 사용이 적다. 그러나 PSD 형태를 유리 함수로 근사해야 하므로 폰 카르만 모델과 같은 비유리 PSD는 근사 오차가 발생한다.

주파수 영역 방법과 시간 영역 방법의 선택은 응용에 따라 다르다. 오프라인 해석에서는 주파수 영역 방법의 정확성이 선호되며, 실시간 시뮬레이션에서는 시간 영역 방법의 효율성이 선호된다.

8. 다변량 난류 생성

공간적으로 상관된 다변량 난류(예: 여러 지점의 난류)를 생성하려면 공간 상관 구조를 포함하는 다변량 PSD 행렬이 필요하다. 크로스-스펙트럼 밀도(Cross-Spectral Density, CSD) $\Phi_{ij}(\omega)$ 가 각 쌍의 지점 간 상관을 기술한다. 다변량 주파수 영역 생성은 각 주파수에서 다변량 가우시안 랜덤 벡터를 생성하여 구성한다. 이러한 방법은 Mann의 “Wind Field Simulation”(Probabilistic Engineering Mechanics, 1998) 등에서 체계적으로 제시되었다.

9. 차원 공간 난류 생성

풍력 터빈이나 대형 기체의 해석에서는 3차원 공간 난류 필드가 필요하다. Mann 모델 또는 Kaimal 모델 등이 표준으로 사용된다. 3차원 공간에서 정의된 PSD 텐서로부터 3차원 FFT를 사용하여 속도 필드를 합성한다. 이러한 필드는 비행체 또는 구조물의 각 지점에서 시간에 따른 난류 속도를 제공한다. 풍력 공학에서는 IEC 표준에 준거한 3차원 난류 필드 생성이 광범위하게 사용된다.

10. 스펙트럼 가우시안 합성 방법

비가우시안 특성을 가진 난류를 생성하기 위해 수정된 주파수 영역 방법이 제안되었다. Reeves, Joshi, and Campbell의 “Development of a Non-Gaussian Atmospheric Turbulence Model for Use in Flight Simulators”(NASA CR-2237, 1974)는 비가우시안 대기 난류 모델을 제시한다. 이 모델은 극단적 돌풍의 발생 확률을 가우시안 모델보다 정확하게 재현하여 비행 안전성 평가에 유용하다.

11. 중첩 시뮬레이션

실제 대기 조건에서는 연속 난류와 이산 돌풍이 공존한다. 주파수 영역 방법으로 생성된 연속 난류에 1-코사인 돌풍과 같은 이산 돌풍을 중첩하여 복합 환경을 재현한다. 이러한 중첩 시뮬레이션은 비행체가 다양한 대기 조건에서 어떻게 반응하는지 종합적으로 평가할 수 있게 한다.

12. 인증 시뮬레이션에서의 활용

항공기 인증을 위한 연속 난류 시뮬레이션에서 주파수 영역 방법은 표준 도구 중 하나이다. CS-25와 FAR Part 25의 연속 난류 해석은 폰 카르만 또는 드라이든 PSD를 입력으로 하는 확률적 응답 해석을 요구한다. 이러한 해석에서 주파수 영역 생성은 정확한 스펙트럼 재현을 보장하여 인증 요건을 만족시킨다.

13. 드론과 UAM에의 적용

드론과 UAM 기체의 고충실도 시뮬레이션에서 주파수 영역 기반 난류 생성이 활용된다. 특히 도심 환경의 복잡한 난류 조건을 정확히 모사하기 위해 사용자 정의 PSD를 사용하는 경우가 많다. 전산유체역학 또는 풍동 시험에서 얻은 스펙트럼을 기반으로 현실적 난류 환경을 재현할 수 있다.

14. 구현 도구

주파수 영역 난류 생성을 위한 도구는 여러 환경에서 제공된다. MATLAB의 Aerospace Toolbox는 드라이든 및 폰 카르만 모델을 지원하며, 주파수 영역 방법도 선택 가능하다. TurbSim 소프트웨어는 풍력 공학에서 3차원 공간 난류 생성의 표준 도구이다. 오픈소스 라이브러리로는 NumPy와 SciPy 기반의 구현이 널리 사용된다.

15. 검증과 품질 보증

생성된 난류의 통계적 특성은 모델 규정과 비교 검증되어야 한다. 주요 검증 항목은 평균, 표준편차, PSD, 그리고 자기 상관 함수이다. 충분히 긴 시계열에서 통계적 수렴이 확인되며, 앙상블 평균을 통해 통계적 신뢰성이 확보된다. 이러한 품질 보증 절차는 시뮬레이션 결과의 신뢰성에 필수적이다.

16. 출처

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.
Mann, J., “Wind Field Simulation,” Probabilistic Engineering Mechanics, Vol. 13, No. 4, 1998.
Kaimal, J. C., Wyngaard, J. C., Izumi, Y., and Coté, O. R., “Spectral Characteristics of Surface-Layer Turbulence,” Quarterly Journal of the Royal Meteorological Society, Vol. 98, 1972.
Shinozuka, M., and Jan, C.-M., “Digital Simulation of Random Processes and Its Applications,” Journal of Sound and Vibration, Vol. 25, No. 1, 1972.
Veers, P. S., “Three-Dimensional Wind Simulation,” Sandia National Laboratories Report SAND88-0152, 1988.
International Electrotechnical Commission, “Wind Turbines—Part 1: Design Requirements,” IEC 61400-1, 2019.

17. 버전

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