등방성 난류

난류 유체역학 전산유체역학 마지막 수정: 2026-08-04 04:52:40

1. 개요[편집]

등방성 난류(homogeneous isotropic turbulence, HIT)는 난류의 통계량이 공간 평행이동에 대해 불변(균질)이고 회전과 반사에 대해서도 불변(등방)인, 현실에는 없지만 이론적으로는 가장 다루기 좋은 이상화된 난류 상태다. 벽도 없고 평균 전단도 없고 대류도 없다. 남는 것은 오직 비선형 상호작용과 점성뿐이라, 난류 연구자들이 난류라는 현상의 뼈대만 뽑아 관찰하는 표준 실험실로 쓴다.

현실의 유동은 죄다 벽이나 전단이나 부력을 달고 있어서 등방성이 아니다. 그럼에도 이 이상화가 반세기 넘게 살아남은 이유는 두 가지다. 첫째, 이 조건에서만 나비에-스토크스로부터 정확한 관계식이 유도된다. 둘째, 콜모고로프의 국소 등방성 가설에 따르면 충분히 작은 스케일은 큰 스케일의 이방성을 잊어버리므로, 소산 영역 근처의 물리는 어떤 유동이든 HIT와 통계적으로 비슷해진다.

이 문서는 이상화의 정의와 그것을 계산기 위에 세우는 방법에 집중한다. 캐스케이드와 $\eta$ 의 정의는 콜모고로프 스케일, 속도 증분의 모멘트는 구조함수, K41 이탈의 상세는 간헐성, DNS의 비용 논리는 직접수치모사에 각각 위임한다.

2. 이상화의 정확한 정의[편집]

균질성은 임의의 $n$ 점 통계가 좌표 원점 선택과 무관하다는 조건이다. 등방성은 그 위에 회전 불변성과 반사 불변성을 더한다. 반사까지 요구하는 것이 중요한데, 반사 대칭을 깨면 평균 헬리시티 $\langle \mathbf{u}\cdot\boldsymbol{\omega}\rangle$ 가 0이 아닌 “비반사 등방 난류”라는 별개의 부류가 된다.

이 세 대칭을 걸면 2점 상관 텐서가 딱 두 개의 스칼라 함수로 압축된다.

R_{ij}(\mathbf{r}) = \langle u_i(\mathbf{x})u_j(\mathbf{x}+\mathbf{r})\rangle = u'^2\left[\frac{f(r)-g(r)}{r^2}r_i r_j + g(r)\,\delta_{ij}\right]

여기서 $f$ 는 종방향, $g$ 는 횡방향 상관함수다. 비압축 조건 $\partial_i R_{ij}=0$ 을 걸면 둘 사이에 관계식 하나가 더 생겨서, 결국 함수 하나로 모든 2점 통계가 결정된다. 텐서 성분 아홉 개짜리 문제가 1변수 함수 하나로 줄어드는 이 압축이 HIT가 이론의 무대가 된 실무적 이유다.

여기서 카르만-하워스 방정식(1938)이 나오고, 고레이놀즈수 극한을 취하면 난류 이론에서 몇 안 되는 엄밀한 결과인 콜모고로프의 4/5 법칙이 떨어진다.

S_3(r) \equiv \langle (\delta u_\parallel)^3\rangle = -\frac{4}{5}\,\varepsilon\, r

계수도 지수도 근사가 아니고, 조정 가능한 상수가 하나도 없다. 난류 이론에서 이런 문장을 쓸 수 있는 결과는 손에 꼽는다.¹

3. 왜 이것이 난류 이론의 실험실인가[편집]

관성 부영역에서 $\varepsilon$ 과 $k$ 만이 통계를 지배한다는 K41의 차원 해석은 그 유명한 $-5/3$ 스펙트럼을 낳는다.

E(k) = C\,\varepsilon^{2/3}k^{-5/3}, \qquad C \approx 1.5\text{–}1.6

이 결과를 오염 없이 검증할 수 있는 유일한 무대가 HIT다. 풍동 격자 난류(콩트-벨로 & 코르신, 1966·1971)나 벽 난류에서는 평균 전단·비등방성·유한 크기 효과가 스펙트럼에 섞여 들어와, 기울기가 $-5/3$ 에서 벗어나도 그게 이론의 문제인지 실험 설정의 문제인지 구분이 안 된다. 주기경계 상자 안의 강제 HIT에는 변명거리가 없다.

같은 이유로 K41의 수정판들도 여기서 검증된다. 소산이 얇은 구조에 몰려 있는 간헐성 때문에 고차 구조함수의 스케일링 지수 $\zeta_p$ 는 $p/3$ 에서 아래로 휘고, 대표적인 두 처방은 다음과 같다.

\text{K62: } \zeta_p = \frac{p}{3} - \frac{\mu}{18}p(p-3), \quad \mu \approx 0.25

\text{She-Lévêque: } \zeta_p = \frac{p}{9} + 2\left[1-\left(\frac{2}{3}\right)^{p/3}\right]

둘 다 $\zeta_3 = 1$ 을 정확히 만족하도록 설계됐다(4/5 법칙이 강제하는 조건이다). $p \le 8$ 쯤까지는 She-Lévêque 쪽이 데이터에 잘 맞고, K62는 큰 $p$ 에서 지수가 다시 내려가는 비물리적 거동을 보인다. 이 구분을 지을 수 있으려면 $p=8$ 이상의 모멘트를 통계적으로 수렴시켜야 하는데, 그게 곧 HIT DNS의 상자 크기와 적분 시간을 결정한다.

4. 표준 DNS 설정[편집]

강제 등방난류 DNS는 이 바닥에서 가장 표준화된 계산이다. 레시피가 거의 고정돼 있다.

도메인: $[0,2\pi]^3$ 삼중 주기 상자(주기경계조건). 균질성을 이산화 수준에서 강제하는 가장 싼 방법이다.
공간 이산화: 의사스펙트럼 방법. 비선형 항만 물리공간에서 계산하고 나머지는 파수공간에서 처리하며, 앨리어싱은 2/3 규칙으로 잘라낸다( $k_{max} = N/3$ ).
시간 적분: 점성 항은 적분인자로 정확히 처리하고 비선형 항만 RK2/RK3로 전진시키는 조합이 표준. 점성 항의 강성을 시간 간격 제약에서 빼내는 것이 목적이다.
압력: 파수공간 사영 연산자 $\left(\delta_{ij} - k_ik_j/k^2\right)$ 한 방으로 비압축성이 기계 정밀도로 만족된다. 포아송 방정식을 반복법으로 풀 필요가 없다는 것이 스펙트럴 HIT의 반칙급 이점이다.

강제(forcing)를 넣는 이유는 단순하다. 아무것도 안 하면 난류는 점성에 먹혀 감쇠하고, 통계를 모을 정상상태가 생기지 않기 때문이다. 두 계열이 쓰인다.

저파수 강제. $k \lesssim 2$ 인 셸의 에너지를 매 스텝 목표값으로 되돌리거나(결정론적), 오른슈타인-울렌벡 과정으로 그 셸에만 랜덤 힘을 넣는 방식(에스와란 & 포프, 1988). 에너지 주입 스케일과 소산 스케일을 확실히 분리해 준다는 장점이 있지만, 강제 방식이 대규모 구조의 통계를 직접 결정하므로 큰 스케일 결과의 해석에 주의가 필요하다.

선형 강제. 물리공간에서 $\mathbf{f} = A\mathbf{u}$ 를 그냥 더한다(룬드그렌 2003, 로살레스 & 메네보 2005). 정상상태에서 주입률이 $\varepsilon = 2A K$ ( $K$ 는 난류 운동에너지)로 딱 떨어지고, 파수공간을 건드리지 않으므로 비주기 격자·유한체적 코드에도 그대로 이식된다. 대신 강제가 모든 스케일에 걸리고 큰 스케일 요동이 커서, 같은 $Re_\lambda$ 를 얻으려면 저파수 강제보다 상자가 더 커야 한다.

5. $Re_\lambda$ 와 해상도 기준[편집]

HIT의 레이놀즈수는 적분 스케일이 아니라 테일러 마이크로스케일 $\lambda$ 로 정의하는 것이 관례다.

\lambda = u'\sqrt{\frac{15\nu}{\varepsilon}}, \qquad Re_\lambda = \frac{u'\lambda}{\nu} = \sqrt{\frac{20}{3}Re_L}

$\lambda$ 자체는 물리적 의미가 뚜렷한 스케일이 아니라 상관함수의 곡률에서 나오는 중간 척도지만, 실험에서 재기 쉽고 코드 사이 비교가 쉬워서 사실상 HIT의 신분증이 되었다. $Re_\lambda \approx 100$ 이면 관성 부영역이 겨우 생기기 시작하고, 200~400이면 스펙트럼에서 $-5/3$ 구간을 눈으로 확인할 수 있으며, 고차 통계로 간헐성을 논하려면 최소 수백 이상이 필요하다.

해상도 판정은 $\eta$ 대비 최대 해상 파수로 한다.

k_{max}\,\eta \gtrsim 1.5

이 값의 근거는 소산 스펙트럼 $2\nu k^2E(k)$ 가 $k\eta \approx 0.2$ 부근에서 최대이고 $k\eta \sim 1$ 이후로 급격히 꺼진다는 사실이다. 즉 $k_{max}\eta = 1.5$ 는 “소산의 대부분을 잡았다”는 뜻이지 “다 잡았다”는 뜻이 아니다. 그래서 에너지 스펙트럼과 $\varepsilon$ 만 볼 거면 $1.5$ 로 충분하지만, 속도 기울기의 고차 모멘트나 엔스트로피 $\langle|\boldsymbol{\omega}|^2\rangle/2$ 의 극단 사건을 볼 거면 $k_{max}\eta \gtrsim 3$ , 심하면 $6$ 까지 요구된다.² 2/3 규칙을 감안하면 $k_{max}\eta = 1.5$ 는 대략 $N\eta \approx 4.5$ , 즉 콜모고로프 스케일당 격자 4~5개에 해당한다.

여기에 반대쪽 제약이 붙는다. 적분 스케일 $L$ 이 상자의 1/5~1/8보다 크면 가장 큰 와가 자기 주기 이미지와 상관되어 통계가 오염된다. 위아래에서 동시에 조이니 $N$ 은 직접수치모사 문서가 말하는 $Re^{9/4}$ 법칙 그대로 폭발한다. 실제로 지구 시뮬레이터의 $4096^3$ ( $Re_\lambda \approx 1200$ , 카네다 외 2003)과 그 후속인 $12288^3$ ( $Re_\lambda \approx 2300$ )이 오랫동안 이 분야의 상한선 역할을 했다.

6. 어디에 쓰는가[편집]

모델 검증대. 대와류 모사의 부격자 모델은 거의 예외 없이 HIT에서 먼저 검증된다. 스마고린스키 상수 $C_s\approx 0.17$ 은 관성 부영역의 에너지 전달률을 맞추도록 HIT에서 유도된 값이고, 동적 모델 역시 여기서 교정된다.
감쇠 난류 벤치마크. 강제를 끄면 $K \sim t^{-n}$ 으로 감쇠하는데, 지수 $n$ 이 초기 스펙트럼의 저파수 거동에 따라 새프먼형 $6/5$ 와 배츨러형 $10/7$ 로 갈린다. 콩트-벨로 & 코르신의 풍동 데이터가 지금도 LES 코드의 통과 의례로 쓰인다.
입자·연소·혼합 문제의 배경장. 관성 입자의 우선 집중(preferential concentration), 스칼라 혼합, 예혼합 화염 주름 연구는 전부 HIT 상자를 배경으로 깐다. 유동 자체를 통제 변수에서 빼기 위해서다.
기계학습 클로저의 학습 데이터. 부격자 응력을 학습하는 신경망 클로저의 절대다수가 HIT DNS 데이터로 훈련된다. 그리고 그 모델이 벽 난류에서 잘 안 된다는 보고 역시 절대다수다.³

이상화가 이만큼 쓰이는 이유는 결국 하나다. 난류에서 유일하게 변수를 통제할 수 있는 설정이기 때문이다. 물론 그 대가로 현실에서 제일 중요한 것들(벽, 전단, 부력)을 전부 빼놓았으니, HIT에서 잘 되는 모델이 실제 유동에서 잘 되리라는 보장은 어디에도 없다.⁴

7. 관련 문서[편집]

4/5 법칙이 특별한 이유는 비선형 항을 근사하지 않고 그대로 다뤘기 때문이다. 대신 조건이 빡세다 — 균질·등방·비압축·정상상태·무한 레이놀즈수. 이 다섯 개를 다 만족하는 실제 유동은 지구상에 없다. 그래도 있는 게 없는 것보단 낫다. ↩
이래서 논문의 $k_{max}\eta$ 값은 그 연구가 무엇을 보려 했는지를 알려주는 지문 같은 숫자다. $1.5$ 면 스펙트럼 논문, $3$ 이상이면 기울기·간헐성 논문일 확률이 높다. ↩
등방성 난류로 학습한 모델을 벽 난류에 던지고 “일반화가 안 되네요”라고 말하는 것은, 사막에서 배운 운전으로 빙판길에 나가는 것과 비슷하다. 문제는 데이터가 사막에만 있다는 것이다. ↩
그럼에도 신규 CFD 코드가 나오면 제일 먼저 돌리는 게 HIT다. 정답 통계가 공개돼 있고, 코드가 틀리면 스펙트럼 꼬리가 즉시 들려 올라가서 숨길 수가 없기 때문이다. 난류 코드계의 리트머스 시험지. ↩