구조함수

난류 유체역학 통계 마지막 수정: 2026-07-28 05:06:33

1. 개요[편집]

구조함수(structure function)는 난류 속도장에서 거리 $r$ 만큼 떨어진 두 점의 속도 차 — 속도 증분 — 의 $p$ 차 모멘트다. 종방향 증분을

\delta u_\parallel(r) = \left[\mathbf{u}(\mathbf{x}+\mathbf{r}) - \mathbf{u}(\mathbf{x})\right]\cdot \hat{\mathbf{r}}

로 정의하면 $p$ 차 종방향 구조함수는

S_p(r) = \left\langle \left(\delta u_\parallel(r)\right)^p \right\rangle

이다. 별것 아닌 정의처럼 보이지만, 이 한 줄이 난류 통계 이론의 절반을 떠받친다. 왜냐하면 구조함수는 평균 유동을 자동으로 지워버리기 때문이다. 두 점의 속도를 빼면 균일한 이송 성분이 상쇄되고, 남는 것은 크기 $r$ 정도의 소용돌이가 만드는 요동뿐이다. 그래서 $S_p(r)$ 을 $r$ 의 함수로 보면 스케일별 통계가 곧바로 읽힌다.

콜모고로프 스케일이 “어느 크기까지 소용돌이가 존재하는가”를 다루고 간헐성이 “그 통계가 어떻게 K41에서 벗어나는가”를 다룬다면, 이 문서는 그 논의의 측정 도구 자체 — 정의, 정확한 결과, 그리고 실제로 재려 할 때 무엇이 사람을 괴롭히는가 — 를 다룬다.

2. 종방향과 횡방향[편집]

증분을 분리 벡터 방향으로 투영하면 종방향, 수직 방향으로 투영하면 횡방향 $\delta u_\perp$ 다. 둘은 독립이 아니다. 등방·비압축 난류에서 연속방정식이 두 2차 구조함수를 묶는다.

S_2^\perp(r) = S_2^\parallel(r) + \frac{r}{2}\frac{dS_2^\parallel}{dr}

관성 부영역에서 $S_2^\parallel \propto r^{2/3}$ 을 넣으면 $S_2^\perp / S_2^\parallel = 4/3$ 이라는 깔끔한 값이 나온다. 실험에서 이 비가 $4/3$ 에 얼마나 가까운지가 등방성이 얼마나 잘 성립하는지를 재는 1차 체크리스트다. 벽 근처나 전단이 강한 유동에서는 이 값이 한참 벗어나고, 그러면 K41류 해석 자체가 위태로워진다.¹

3. 세 영역과 K41[편집]

구조함수는 $r$ 에 따라 세 구간을 보인다.

구간	범위	$S_2$ 거동	물리
소산영역	$r \lesssim \eta$	$\propto r^2$	속도장이 매끄러워 $\delta u \approx r\,\partial u/\partial x$
관성영역	$\eta \ll r \ll L$	$\propto r^{2/3}$	점성도 강제도 안 보이는 캐스케이드
적분 스케일	$r \gtrsim L$	포화 $\to 2\langle u'^2\rangle$	두 점의 속도가 무상관

관성영역에서 K41은 차원 해석만으로 $S_2(r) = C_2 (\varepsilon r)^{2/3}$ 을 준다( $C_2 \approx 2.0$ ). 그런데 K41에는 차원 해석이 아니라 나비에-스토크스에서 직접 유도되는 결과가 하나 있다. 그 유명한 4/5 법칙이다.

S_3(r) = -\frac{4}{5}\,\varepsilon r

균질·등방·정상 상태이고 레이놀즈수가 무한대인 극한에서, 폰 카르만-하워스 방정식으로부터 어떤 모형 가정도 없이 나온다. 난류 이론에서 이 급의 정확한 결과는 손에 꼽는다. 부호가 음수라는 것도 중요한데, 이는 에너지가 큰 스케일에서 작은 스케일로 흐른다는 캐스케이드의 방향을 통계량 하나로 못 박은 것이다.² 더 일반적으로는 $\langle \delta u_\parallel |\delta \mathbf{u}|^2 \rangle = -\tfrac{4}{3}\varepsilon r$ 형태의 4/3 법칙이 있고, 스칼라 수송에도 대응물이 존재한다.

4. 에너지 스펙트럼과의 대응[편집]

$S_2(r)$ 과 에너지 스펙트럼 $E(k)$ 는 상관함수를 매개로 푸리에 쌍으로 엮여 있다. $E(k) \propto k^{-5/3}$ 이 $S_2 \propto r^{2/3}$ 에 대응하는 것이 그 예다. 다만 대응이 항상 되는 건 아니다. 스펙트럼 지수가 $-3$ 보다 가파르면 $S_2$ 는 $r^2$ 에서 포화해 버려 더 이상 지수를 구분하지 못한다. 매끄러운 장에서 $S_2 \propto r^2$ 이 나오는 것과 같은 이유다. 그래서 스펙트럼이 급한 계(예: 2D 난류의 엔스트로피 영역)에서는 구조함수로 지수를 읽으려 하면 안 된다.

반대 방향의 장점도 분명하다. 구조함수는 주기 경계도, 균일 격자도, 창 함수도 필요 없다. 불규칙하게 샘플링된 실험 데이터나 비주기 도메인에서도 두 점만 있으면 계산되기 때문에, 대기 경계층·해양·태양풍 관측처럼 FFT를 쓸 수 없는 데이터에서는 사실상 유일한 스케일 분석 도구다.

5. 스케일링 지수와 확장 자기유사성[편집]

관성영역에서 $S_p(r) \propto r^{\zeta_p}$ 로 놓으면, K41은 $\zeta_p = p/3$ 을 예측한다. 실제 측정치는 다르다.

$p$	K41 $p/3$	측정 $\zeta_p$	She-Lévêque
2	0.667	0.70	0.696
3	1.000	1.00	1.000
4	1.333	1.28	1.280
6	2.000	1.78	1.778
8	2.667	2.17	2.211

고차로 갈수록 $\zeta_p$ 가 직선 아래로 휘어지는데, 이 휘어짐이 곧 간헐성의 정량적 정의다. 1994년 She와 Lévêque가 제안한

\zeta_p = \frac{p}{9} + 2\left[1 - \left(\frac{2}{3}\right)^{p/3}\right]

는 조절 파라미터 없이 측정치를 놀랄 만큼 잘 맞힌다(정의상 $\zeta_3 = 1$ 을 정확히 만족한다). 소산 구조가 1차원 실 형태(로그-푸아송 계층)라는 가정에서 나온 결과다. 다중프랙탈 형식론은 같은 곡선을 특이도 스펙트럼의 르장드르 변환으로 서술한다.

문제는 실험실 레이놀즈수에서 관성영역 자체가 반 자릿수도 안 나온다는 것. 로그-로그 그래프에 직선 구간이 없는데 기울기를 읽으라니 곤란하다. 여기서 나온 트릭이 확장 자기유사성(ESS, Benzi 등 1993)이다. $S_p$ 를 $r$ 이 아니라 $|S_3|$ 에 대해 그리면 훨씬 넓은 범위에서 깨끗한 멱법칙이 나타나고, 그 기울기가 $\zeta_p/\zeta_3 = \zeta_p$ 를 준다. 저Re DNS와 실험에서 지수를 뽑는 표준 절차가 됐다. 다만 ESS는 “왜 되는지”에 대한 완전한 이론적 정당화가 없는 경험적 처방에 가깝고, 유한Re 보정과 비등방성이 남아 있으면 여전히 편향된다.

6. 실제로 재는 법과 수렴 요건[편집]

실험에서는 보통 핫와이어 한 개의 시간 신호를 쓴다. 테일러 동결 가설( $r = U\tau$ )로 시간 지연을 공간 거리로 바꾸는데, 이 근사는 난류 강도 $u'/U$ 가 충분히 작아야 성립한다. 벽 근처나 후류처럼 $u'/U$ 가 20 %를 넘으면 동결 가설 자체가 오차원이 된다. **직접수치모사**에서는 3차원 장이 통째로 있으니 여러 방향·여러 시점에 대해 직접 평균할 수 있어 이 문제가 없다. 대신 도메인 크기가 적분 스케일의 몇 배밖에 안 되면 큰 $r$ 쪽 통계가 상자에 갇힌다.

가장 자주 간과되는 게 고차 모멘트의 표본 요구량이다. $S_p$ 는 PDF의 꼬리에 지배되고, 간헐성 때문에 그 꼬리는 가우스보다 훨씬 두껍다. 실무 규칙은 이렇다.

적분항을 직접 그려본다. $\delta u^p\, P(\delta u)$ 곡선이 양 끝에서 확실히 0으로 떨어지지 않으면 그 $p$ 차 모멘트는 수렴하지 않은 것이다. 지수만 뽑아 놓고 이 그림을 안 그리면 리뷰어에게 잡힌다.
유한 표본은 꼬리를 잘라먹으므로 $\zeta_p$ 를 체계적으로 과대 추정하는 방향(K41 쪽으로)으로 편향된다. 즉 표본이 부족하면 간헐성이 실제보다 약해 보인다.
홀수 차수는 양·음 기여가 상쇄돼 신호 대 잡음이 더 나쁘다. 그래서 $|\delta u|^p$ 의 절댓값 구조함수를 따로 쓰는 관행이 있다.
아무리 좋은 데이터라도 신뢰 한계는 대체로 $p \lesssim 8$ 정도. 그 위 값을 보고하는 논문은 표본 수를 반드시 명시해야 한다.³

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

모델 검증: 대와류 모사 하위격자 모델이나 난류 모델링 결과가 관성영역 스케일링을 재현하는지 보는 표준 지표. 평균장은 맞는데 $S_p$ 기울기가 틀린 모델이 흔하다.
수치 소산 진단: 격자 스케일 근처에서 $S_2$ 가 이론보다 빨리 꺾이면 그건 물리적 소산이 아니라 수치 소산과 분산이다.
관측 데이터 분석: 대기 경계층, 해양 혼합층, 태양풍 자기유체 난류. 압축성 난류에서는 밀도 가중 증분 $\rho^{1/3}\mathbf{u}$ 를 써야 $-4/5$ 류 관계가 회복된다.
실험 설계: $\zeta_p$ 를 신뢰도 있게 뽑으려면 몇 개의 표본이 필요한지 역산하는 데 위 수렴 규칙이 그대로 쓰인다.

8. 관련 문서[편집]

실험 논문에서 “등방성이 잘 만족되었다”는 문장 뒤에 이 $4/3$ 검증 그래프가 없으면, 그 문장은 대체로 희망 사항이다. ↩
4/5 법칙은 난류 이론에서 몇 안 되는 “정확한 결과”라 종종 성지 취급을 받는다. 다만 조건이 무한 레이놀즈수라, 실제 데이터에서 $-S_3/(\varepsilon r)$ 를 그려보면 $4/5$ 에 도달하지 못하고 $0.6$ 언저리에서 멈추는 경우가 흔하다. 유한Re 보정을 넣으면 설명이 되지만, 처음 그려본 대학원생은 대개 코드부터 의심하며 밤을 새운다. ↩
“12차 구조함수를 측정했다”는 주장은 곧 “꼬리에서 열 몇 개의 표본이 결과를 다 결정했다”는 뜻일 수 있다. 이 바닥에서 고차 모멘트는 통계량이 아니라 신앙 고백에 가까워지는 지점이 있다. ↩