27.19 공간적 상관을 고려한 난류 모델링

1. 공간적 상관의 개념과 중요성

난류는 본질적으로 3차원 공간 현상이며, 공간적으로 인접한 두 지점에서의 속도 변동은 서로 상관된다. 이러한 공간적 상관(spatial correlation)은 난류 와류의 특성 크기에 의해 결정되며, 두 지점 간 거리가 적분 길이 스케일보다 훨씬 작으면 상관이 높고, 훨씬 크면 상관이 거의 없다. 공간적 상관을 고려한 난류 모델링은 대형 비행체의 공력 해석, 풍력 터빈 설계, 그리고 다중 드론 군집 비행의 정확한 시뮬레이션에 필수적이다.

2. 공간 상관 함수의 정의

등방성 균질 난류에서 두 지점 $\mathbf{x}$ 와 $\mathbf{x} + \mathbf{r}$ 의 속도 변동 간 공간 상관은 다음과 같이 정의된다.

$R_{ij}(\mathbf{r}) = \overline{u_i'(\mathbf{x}) u_j'(\mathbf{x} + \mathbf{r})}$

여기서 $i, j$ 는 속도 성분 지수이다. 등방성 난류에서 이 상관 텐서는 두 개의 스칼라 함수로 표현된다. 종방향 상관 함수 $f(r) = R_{LL}(r)/\sigma^2$ 와 횡방향 상관 함수 $g(r) = R_{NN}(r)/\sigma^2$ 이다. 여기서 $L$ 은 $\mathbf{r}$ 방향의 속도 성분, $N$ 은 $\mathbf{r}$ 에 수직인 성분이다.

연속 방정식으로부터의 제약

비압축성 연속 방정식 $\nabla \cdot \mathbf{u} = 0$ 으로부터 종방향과 횡방향 상관 함수 간의 관계가 유도된다.

$g(r) = f(r) + \frac{r}{2} \frac{df(r)}{dr}$

이 관계는 등방성 난류 상관 구조의 근본적 제약이다. 이로부터 종방향 상관이 알려지면 횡방향 상관이 자동으로 결정된다. 이러한 관계는 비압축성 가정 하에서 모든 등방성 난류 모델(드라이든, 폰 카르만 등)이 만족해야 하는 기본 조건이다.

3. 적분 길이 스케일의 정의

적분 길이 스케일은 공간 상관 함수의 적분으로 정의된다.

$L_{ii} = \int_0^{\infty} \frac{R_{ii}(r)}{\sigma_i^2} dr$

이 길이 스케일은 속도 변동이 유의하게 상관되는 최대 공간 거리를 나타내며, 난류의 특성 크기이다. 대기 경계층에서 $L_u$ 는 수십 미터에서 수백 미터 수준이다. 적분 길이 스케일은 난류 와류의 전형적 크기로 해석될 수 있다.

비행체 규모와 공간 상관

비행체가 경험하는 난류의 공간 상관 구조는 비행체 규모와 적분 길이 스케일의 비율에 의존한다. 비행체 규모가 $L$ 보다 훨씬 작으면 비행체 전체에 걸쳐 거의 균일한 난류가 작용하므로 공간 상관을 무시할 수 있다. 비행체 규모가 $L$ 에 비교할 만하거나 크면 비행체의 다른 부분에서 상이한 난류가 작용하며, 공간 상관의 반영이 필수적이다.

소형 드론에서의 공간 상관

소형 드론(크기 수십 cm에서 수 m)은 일반적으로 적분 길이 스케일보다 훨씬 작으므로 점 근사(point approximation)가 타당하다. 드론 전체에 동일한 난류 속도가 작용한다고 가정할 수 있다. 그러나 드론 프로펠러 직경 스케일에서는 작은 스케일 난류의 영향이 있을 수 있으며, 프로펠러 간의 상관도 고려될 수 있다. 이러한 소규모 공간 상관은 대부분 무시되나, 특수 해석에서는 반영될 수 있다.

대형 항공기와 UAM의 공간 상관

대형 항공기의 날개 폭이 적분 길이 스케일에 비교할 만할 때 날개 양단에서 경험하는 난류가 부분적으로 상관된다. 이 상관은 비행체에 롤링 모멘트를 유발한다. 전체 날개에 걸쳐 균일한 상승 돌풍은 수직 가속도만을 유발하지만, 공간적으로 변동하는 돌풍은 비대칭 하중을 생성한다. UAM 기체의 다중 로터 구성에서도 로터 간 상관이 기체 자세 응답에 영향을 미친다.

다변량 모델링

공간적 상관을 반영한 난류 모델링은 다변량 확률 과정으로 접근한다. $M$ 개의 공간 지점에서의 속도 변동 $u_1(t), u_2(t), \ldots, u_M(t)$ 는 $M$ -차원 다변량 확률 과정이다. 이 과정의 특성은 파워 스펙트럼 밀도 행렬 $\Phi_{ij}(\omega)$ 로 기술되며, 대각선 요소는 각 지점의 PSD, 비대각선 요소는 지점 간 크로스-PSD이다. 다변량 모델링은 이 행렬을 바탕으로 상관된 시계열을 생성한다.

Cholesky 분해 기반 합성

다변량 상관 시계열 생성의 대표적 방법은 Cholesky 분해 기반 합성이다. 각 주파수 $\omega_n$ 에서 PSD 행렬 $\mathbf{\Phi}(\omega_n)$ 을 Cholesky 분해한다.

$\mathbf{\Phi}(\omega_n) = \mathbf{L}(\omega_n) \mathbf{L}^H(\omega_n)$

여기서 $\mathbf{L}$ 은 하삼각 행렬, $\mathbf{L}^H$ 는 에르미트 공액 전치이다. 독립 가우시안 랜덤 벡터 $\mathbf{w}$ 에 $\mathbf{L}$ 을 곱하면 PSD 행렬에 해당하는 다변량 스펙트럼이 얻어진다. 이를 역 푸리에 변환하여 시간 시계열을 생성한다.

4. Mann 모델

Mann의 “Wind Field Simulation”(Probabilistic Engineering Mechanics, 1998)에서 제시된 Mann 모델은 3차원 공간의 상관 난류 필드를 생성하는 대표적 방법이다. 이 모델은 Rapid Distortion Theory에 기반하여 대기 경계층의 이방성 난류를 재현한다. 풍력 터빈 하중 계산 및 대기 시뮬레이션에서 표준으로 사용된다. 국제 전기 기술 위원회(IEC) 표준 IEC 61400-1에서도 이 모델이 권장된다.

5. Kaimal 모델

Kaimal 모델은 대기 경계층 난류의 스펙트럼과 공간 상관을 경험적으로 기술하는 모델이다. Kaimal, Wyngaard, Izumi, and Coté의 “Spectral Characteristics of Surface-Layer Turbulence”(Quarterly Journal of the Royal Meteorological Society, 1972)에서 제시되었다. 이 모델은 단순성과 실측 데이터와의 좋은 일치로 풍력 공학에서 광범위하게 사용된다. IEC 61400-1의 대안 모델로도 채택되어 있다.

6. 공간 주파수 표현

3차원 공간에서 정의된 난류의 주파수 기술은 공간 파수 벡터 $\mathbf{k} = (k_1, k_2, k_3)$ 를 사용한다. 에너지 스펙트럼 텐서 $\Phi_{ij}(\mathbf{k})$ 는 3차원 공간 스펙트럼을 기술한다. 이로부터 1차원 스펙트럼 및 특정 방향의 공간 상관이 유도된다. 3차원 FFT를 사용하여 공간 난류 필드가 효율적으로 합성된다.

7. 비등방성 공간 모델

실제 대기 경계층 난류는 완전한 등방성을 가지지 않는다. 종방향, 횡방향, 수직 방향 난류가 서로 다른 강도와 길이 스케일을 가진다. 비등방성 모델은 이러한 방향별 차이를 명시적으로 반영한다. Mann 모델은 Reynolds 응력 모델과 결합하여 비등방성을 재현한다. 이러한 모델은 정확한 공력 해석에 필수적이다.

8. 도심 환경의 공간 상관

도심 환경에서 난류의 공간 상관은 일반적인 균질 난류와 크게 다르다. 건물 후류, 빌딩 캐니언 와류, 국지적 경계층 변형 등이 공간 상관을 비균질적으로 만든다. 이러한 환경에서 정확한 난류 모델링은 전산유체역학 또는 현장 측정에 기반한 경험 모델이 필요하다. UAM 기체의 도심 운용 시뮬레이션에서 이러한 복잡한 공간 구조가 반영되어야 한다.

9. 수치 구현과 계산 비용

공간적 상관을 고려한 난류 생성은 균질 단일 지점 모델보다 계산 비용이 크다. 3차원 FFT 기반 생성은 격자 크기의 세제곱에 비례하는 계산 시간을 요구한다. 메모리 사용량도 상당하다. 따라서 응용의 요구 정확도와 계산 자원에 따라 적절한 모델과 해상도가 선택된다. 실시간 시뮬레이션에서는 단순화된 모델이 선호되며, 오프라인 해석에서는 고충실도 모델이 사용된다.

10. 검증과 활용

공간 상관 난류 모델의 검증은 실측 또는 고충실도 CFD 시뮬레이션과의 비교로 수행된다. 검증 항목은 PSD, 공간 상관, 코히어런스 함수 등이다. 검증된 모델은 항공기 구조 설계, 풍력 터빈 피로 해석, 대형 건물 풍하중 평가 등에서 활용된다. 드론 및 UAM 분야에서도 공간 상관 난류 모델의 활용이 확대되고 있다.

11. 출처

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.
Batchelor, G. K., “The Theory of Homogeneous Turbulence,” Cambridge University Press, 1953.
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.
Shinozuka, M., and Deodatis, G., “Simulation of Stochastic Processes by Spectral Representation,” Applied Mechanics Reviews, Vol. 44, No. 4, 1991.

12. 버전

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