27.3 난류 에너지 스펙트럼과 에너지 캐스케이드

1. 난류 에너지 스펙트럼의 개념

난류 에너지 스펙트럼(turbulence energy spectrum) $E(k)$ 는 난류 운동 에너지를 공간 파수 $k$ 의 함수로 분포시킨 함수이다. 파수 $k$ 는 공간 주파수이며, 단위는 rad/m이다. 에너지 스펙트럼은 다음의 적분 관계를 만족한다.

$\int_0^{\infty} E(k) \, dk = \frac{1}{2} \overline{u_i' u_i'} = k_{TKE}$

여기서 $k_{TKE}$ 는 난류 운동 에너지이다. 에너지 스펙트럼은 난류의 스케일별 에너지 분포를 기술하는 핵심 도구이며, 난류 이론과 실용 응용의 중심 개념이다.

에너지 캐스케이드의 개념

에너지 캐스케이드(energy cascade)는 난류에서 에너지가 큰 스케일 와류에서 점진적으로 작은 스케일 와류로 전달되는 과정을 지칭한다. 이 개념은 Richardson의 1922년 시 “Big whirls have little whirls that feed on their velocity, and little whirls have lesser whirls and so on to viscosity“로 유명하며, Richardson의 “Weather Prediction by Numerical Process”(Cambridge University Press, 1922)에서 제시되었다. 에너지는 가장 큰 스케일에서 외부 입력(바람 전단, 부력, 대규모 기상 등)에 의해 공급되며, 와류 불안정성을 통해 점진적으로 작은 와류로 분할된다. 최종적으로는 Kolmogorov 소산 스케일에서 점성에 의해 열로 변환된다.

에너지 스펙트럼의 세 영역

난류 에너지 스펙트럼은 다음 세 주요 영역으로 구분된다. 첫째, 에너지 생성 영역(energy-containing range)은 큰 스케일 와류가 평균 유동 또는 외부에서 에너지를 공급받는 영역이다. 이 영역에서 스펙트럼은 적분 길이 스케일 $L$ 에 의해 특성화된다. 둘째, 관성 부영역(inertial subrange)은 에너지 생성이나 소산이 지배적이지 않고 에너지가 스케일 간 전달되는 영역이다. 이 영역에서 Kolmogorov의 $-5/3$ 법칙이 성립한다. 셋째, 소산 영역(dissipation range)은 Kolmogorov 스케일 근방에서 점성 소산이 지배적인 영역이다. 이 영역에서 스펙트럼은 지수 함수적으로 감소한다.

Kolmogorov의 5/3 법칙

관성 부영역에서의 에너지 스펙트럼은 다음의 보편적 형태를 가진다.

$E(k) = C_K \varepsilon^{2/3} k^{-5/3}$

여기서 $C_K \approx 1.5$ 는 Kolmogorov 상수, $\varepsilon$ 은 난류 에너지 소산율이다. 이 관계는 Kolmogorov의 1941년 논문 “The Local Structure of Turbulence in Incompressible Viscous Fluid for Very Large Reynolds Numbers”(Doklady Akademii Nauk SSSR, 1941)에서 도출되었다. 5/3 법칙은 대기 난류의 실험적 측정에서 광범위하게 확인되었으며, 난류 이론의 핵심 결과로 평가된다.

2. 차원 해석에 의한 유도

Kolmogorov의 5/3 법칙은 차원 해석으로 유도될 수 있다. 관성 부영역에서 에너지 스펙트럼 $E(k)$ 는 파수 $k$ 와 에너지 소산율 $\varepsilon$ 에만 의존한다고 가정하면(국지적 등방성과 관성 부영역의 가정), 차원 고려에 의해

$[E(k)] = \frac{\text{m}^3}{\text{s}^2}, \quad [k] = \frac{1}{\text{m}}, \quad [\varepsilon] = \frac{\text{m}^2}{\text{s}^3}$

이므로 $E(k) \propto \varepsilon^{2/3} k^{-5/3}$ 이 유일한 차원 조합이다. 이러한 단순한 차원 해석이 실제 난류의 정량적 특성을 예측할 수 있다는 점은 Kolmogorov 이론의 깊은 통찰을 보여준다.

적분 스케일과 에너지 공급

에너지 생성 영역은 적분 길이 스케일 $L$ 에 의해 특성화된다. 이 스케일에서 에너지는 평균 유동의 운동 에너지에서 난류 운동 에너지로 전달된다. 스펙트럼의 피크 파수는 약 $k_{peak} \sim 1/L$ 이며, 이 영역에서의 에너지 밀도는 난류 변동 분산의 대부분을 차지한다. 대기 경계층의 경우 적분 길이 스케일은 수십 미터에서 수백 미터 범위이며, 이에 해당하는 피크 주파수는 드론 운용 속도에서 수 헤르츠에서 수십 헤르츠 범위이다.

Kolmogorov 소산 스케일

에너지가 점성에 의해 열로 변환되는 최소 스케일은 Kolmogorov 스케일이라 불린다. 이 스케일의 특성량은 다음과 같이 정의된다.

$\eta = \left(\frac{\nu^3}{\varepsilon}\right)^{1/4}, \quad \tau_\eta = \left(\frac{\nu}{\varepsilon}\right)^{1/2}, \quad u_\eta = (\nu \varepsilon)^{1/4}$

여기서 $\nu$ 는 운동 점성 계수, $\varepsilon$ 은 에너지 소산율이다. 대기 조건에서 $\eta$ 는 약 1 mm 수준이며, $\tau_\eta$ 는 약 0.1초 수준이다. 이 스케일에서 Reynolds 수가 1에 가까워지며 점성이 지배적이 된다.

3. 스펙트럼의 실험적 검증

Kolmogorov 스펙트럼은 다양한 실험 환경에서 검증되었다. Grant, Stewart, and Moilliet의 “Turbulence Spectra from a Tidal Channel”(Journal of Fluid Mechanics, 1962)는 해양 유동에서 5/3 법칙을 검증한 초기 연구이다. 대기 경계층에서도 수많은 측정이 수행되어 관성 부영역의 존재와 Kolmogorov 법칙의 성립이 확인되었다. 다만 대기의 경우 안정도에 따라 소규모 편차가 관찰될 수 있다.

4. 공간 파수에서 주파수로의 변환

실제 측정에서는 고정 위치에서 시간에 따라 변동하는 속도가 측정되므로 주파수 $\omega$ 에서의 스펙트럼이 얻어진다. Taylor의 고립 가설(frozen turbulence hypothesis)을 적용하면 평균 속도 $U$ 로 이류되는 난류의 시간 스펙트럼 $\Phi_{uu}(\omega)$ 와 공간 스펙트럼 $E(k)$ 가 다음의 관계를 가진다.

$\Phi_{uu}(\omega) = \frac{E(k)}{U}, \quad \text{with} \quad k = \omega/U$

이 관계는 대기 난류의 실용적 분석에서 필수적이다.

3차원 스펙트럼 대 1차원 스펙트럼

에너지 스펙트럼은 3차원 등방 스펙트럼 $E(k)$ 과 1차원 스펙트럼 $E_{ij}(k_1)$ 으로 구분된다. 등방 난류에서는 이들 간에 다음과 같은 관계가 성립한다.

$E_{11}(k_1) = \int_{k_1}^{\infty} \frac{E(k)}{k^3}(k^2 - k_1^2) dk$

여기서 $E_{11}$ 은 종방향 속도 성분 $u$ 의 1차원 스펙트럼이다. 실제 측정에서는 1차원 스펙트럼이 직접 얻어지며, 3차원 스펙트럼은 등방성 가정하에 역산된다.

5. 이방성 난류의 스펙트럼

대기 경계층의 난류는 엄밀한 의미에서 등방성이 아니며, 종방향, 횡방향, 수직 방향 스펙트럼이 서로 다르다. 일반적으로 종방향 변동이 가장 크고, 수직 변동이 가장 작으며, 횡방향은 중간 수준이다. 이러한 이방성은 적분 길이 스케일의 방향성에도 반영된다. Dryden 모델 및 von Kármán 모델에서는 이러한 방향별 특성이 명시적으로 반영되어 있다.

6. 에너지 전달률

관성 부영역에서의 에너지 전달률(energy transfer rate)은 파수 공간에서 일정하며, 이는 에너지 소산율 $\varepsilon$ 과 같다. 수학적으로는 각 파수 $k$ 에서의 에너지 전달은 3차 속도 상관(triadic interaction)을 통해 이루어진다. 이러한 에너지 전달의 상세 분석은 Batchelor의 “The Theory of Homogeneous Turbulence”(Cambridge University Press, 1953)와 같은 교과서에서 체계적으로 다루어진다.

7. 드론 응용에서의 의의

드론의 크기는 일반적으로 수십 센티미터에서 수 미터 범위이며, 이는 대기 난류의 관성 부영역 또는 소산 영역 시작 부분에 해당한다. 따라서 드론에 작용하는 공력 외란은 주로 관성 부영역 이하의 스케일 난류에 의해 결정되며, Kolmogorov 법칙에 기반한 스펙트럼 모델이 효과적으로 적용될 수 있다. 드론 크기와 난류 길이 스케일의 비율은 공력 응답 해석의 중요한 비차원 수이다.

8. 출처

Richardson, L. F., “Weather Prediction by Numerical Process,” Cambridge University Press, 1922.
Kolmogorov, A. N., “The Local Structure of Turbulence in Incompressible Viscous Fluid for Very Large Reynolds Numbers,” Doklady Akademii Nauk SSSR, Vol. 30, 1941.
Grant, H. L., Stewart, R. W., and Moilliet, A., “Turbulence Spectra from a Tidal Channel,” Journal of Fluid Mechanics, Vol. 12, 1962.
Batchelor, G. K., “The Theory of Homogeneous Turbulence,” Cambridge University Press, 1953.
Pope, S. B., “Turbulent Flows,” Cambridge University Press, 2000.
Frisch, U., “Turbulence: The Legacy of A. N. Kolmogorov,” Cambridge University Press, 1995.

9. 버전

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