엔스트로피

편집 역사 토론
난류 유체역학 마지막 수정: 2026-08-05 04:55:09

1. 개요[편집]

엔스트로피
Enstrophy
정의$\mathcal{E} = \tfrac{1}{2}\int |\boldsymbol{\omega}|^2\,dV$
재료와도 $\boldsymbol{\omega} = \nabla\times\mathbf{u}$
3D에서와도 신장으로 생성됨, 비보존
2D에서비점성 극한에서 보존 → 이중 캐스케이드
소산과의 관계$\varepsilon = 2\nu\langle\mathcal{E}\rangle$ (균질 난류)
수치적 의의엔스트로피 보존 이산화(아라카와)가 장기 적분을 살림

에너지는 “얼마나 세게 흐르는가”, 엔스트로피는 “얼마나 심하게 비틀리는가”.

엔스트로피(enstrophy)는 유동장의 와도 크기 제곱을 전체 부피에 대해 적분한 양으로, 보통 다음과 같이 정의한다.

E=12∫V∣ω∣2 dV,ω=∇×u\mathcal{E} = \frac{1}{2}\int_V |\boldsymbol{\omega}|^2\,dV, \qquad \boldsymbol{\omega} = \nabla\times\mathbf{u}

운동에너지가 속도의 L2L^2 노름이라면 엔스트로피는 와도의 L2L^2 노름, 즉 속도 기울기의 크기를 재는 양이다. 계수 12\frac{1}{2}를 붙이지 않는 관례도 흔하니 논문을 읽을 때는 정의부터 확인하는 게 안전하다.1 이름은 그리스어 어근에서 왔고, 열역학의 엔트로피와는 아무 관계도 없다.

에너지가 유동의 “총량”을 말한다면, 엔스트로피는 그 에너지가 어느 스케일에 앉아 있는지를 말해준다. 스펙트럼으로 쓰면 에너지가 ∫E(k)dk\int E(k)dk일 때 엔스트로피는 ∫k2E(k) dk\int k^2 E(k)\,dk라, 작은 스케일(큰 kk)에 훨씬 무겁게 가중된다. 이 k2k^2 하나가 이 문서의 거의 모든 내용을 설명한다.

2. 소산과의 관계[편집]

균질 비압축성 난류에서 단위 부피당 운동에너지 소산율은

ε=ν⟨∣ω∣2⟩=2ν⟨E⟩\varepsilon = \nu\langle|\boldsymbol{\omega}|^2\rangle = 2\nu\langle\mathcal{E}\rangle

로 엔스트로피와 직결된다. 즉 점성이 에너지를 갉아먹는 속도는 그 순간의 엔스트로피가 결정한다. 콜모고로프 스케일에서 ε\varepsilon이 관성영역을 통과한 에너지 플럭스와 같아진다는 사실과 합치면, “작은 소용돌이가 얼마나 많이 만들어졌는가”가 곧 “얼마나 빨리 열로 바뀌는가”다.

여기서 유명한 반전이 있다. 점성 ν\nu를 줄여도 ε\varepsilon은 0으로 가지 않는다 — 대신 엔스트로피가 1/ν1/\nu 꼴로 커져서 곱이 유한하게 남는다. 이것이 소산 이상(dissipative anomaly)이며, 난류가 점성의 크기와 무관하게 에너지를 소모하는 이유다.

3. 3차원 — 와도 신장[편집]

와도 수송 방정식은 (비압축·바로트로픽 가정에서)

DωDt=(ω⋅∇)u+ν∇2ω\frac{D\boldsymbol{\omega}}{Dt} = (\boldsymbol{\omega}\cdot\nabla)\mathbf{u} + \nu\nabla^2\boldsymbol{\omega}

이고, 엔스트로피 수지는

dEdt=∫ω⋅S⋅ω dV−ν∫∣∇ω∣2 dV\frac{d\mathcal{E}}{dt} = \int \boldsymbol{\omega}\cdot\mathbf{S}\cdot\boldsymbol{\omega}\,dV - \nu\int|\nabla\boldsymbol{\omega}|^2\,dV

가 된다(S\mathbf{S}는 변형률 텐서). 첫 항이 와도 신장(vortex stretching)으로, 와관이 잡아 늘여지면 각운동량 보존 때문에 단면이 줄고 회전이 빨라져 엔스트로피가 생성된다. 피겨 스케이팅 선수가 팔을 모으는 그 장면이 유체 버전으로 온 공간에서 동시에 벌어지는 셈.

그래서 3D 난류에서 엔스트로피는 보존량이 아니라 성장하는 양이고, 이것이 에너지가 큰 스케일에서 작은 스케일로 내려가는 순방향 캐스케이드의 동력이다. 덧붙여 3차원 나비에-스토크스 방정식의 미해결 정칙성 문제도 정확히 이 양의 문제다. 엔스트로피가 유한하게 유지되면 해는 매끄럽게 남고(레이-라디젠스카야 계열 결과), 유한 시간에 발산한다면 특이점이 생긴다. 비일-카토-마이다 판정 기준이 와도 최댓값의 시간 적분으로 붕괴 여부를 규정하는 것도 같은 맥락이다.2

4. 2차원 — 이중 캐스케이드[편집]

2차원 유동에서는 와도가 평면에 수직한 스칼라라 (ω⋅∇)u=0(\boldsymbol{\omega}\cdot\nabla)\mathbf{u} = 0이다. 와도 신장이 원리적으로 존재하지 않는다. 비점성 극한에서 와도는 유체 입자를 따라 그냥 실려 가고, 따라서 에너지와 엔스트로피가 둘 다 보존된다.

보존량이 두 개라는 사실만으로 캐스케이드 방향이 강제된다. 피외르토프트(Fjørtoft) 논증이 그것인데, 에너지 ∫E(k)dk\int E(k)dk와 엔스트로피 ∫k2E(k)dk\int k^2E(k)dk를 동시에 지키려면 중간 파수에 주입된 에너지의 대부분은 작은 kk(큰 스케일)로 가고 엔스트로피는 큰 kk(작은 스케일)로 갈 수밖에 없다. 크라이크넌(1967), 리스(1968), 배철러(1969)가 정리한 결론이 이중 캐스케이드다.

영역플럭스스펙트럼
주입 파수보다 큰 스케일에너지 역캐스케이드E(k)∼ε2/3k−5/3E(k)\sim \varepsilon^{2/3}k^{-5/3}
주입 파수보다 작은 스케일엔스트로피 순캐스케이드E(k)∼η2/3k−3E(k)\sim \eta^{2/3}k^{-3}

여기서 η\eta는 엔스트로피 소산율이다. k−3k^{-3} 영역에는 크라이크넌이 지적한 로그 보정이 붙는다. 역캐스케이드는 작은 소용돌이가 합쳐져 더 큰 소용돌이가 되는 과정이라, 2D 난류에서는 큰 와류가 저절로 자라 오래 살아남는다. 목성의 줄무늬와 대적점, 지구 대기의 저기압성 소용돌이가 오래 유지되는 이유가 여기 있다.

5. 준2차원 지구유체[편집]

지구 대기·해양은 정확히 2차원은 아니지만, 수평 규모가 연직 규모보다 압도적으로 크고 회전(코리올리)과 성층이 연직 운동을 억제해 준2차원으로 행동한다. 원시방정식 기반 모델과 준지균 근사에서는 와도 대신 퍼텐셜 와도(PV)가 보존되는 양으로 등장하고, 대응하는 이차량이 퍼텐셜 엔스트로피다.

관측된 대기 스펙트럼(나스트롬-게이지)은 종관 규모에서 k−3k^{-3}, 중규모(수백 km 이하)에서 k−5/3k^{-5/3}로 꺾이는데, 이 꺾임의 해석 자체가 지금도 논쟁거리다. 어느 쪽이든 어느 스케일에서 엔스트로피가 어디로 흐르는가가 논쟁의 축이다.

6. 수치해석 — 엔스트로피를 지키는 이산화[편집]

여기가 시뮬레이션 하는 사람에게 실질적으로 중요한 부분이다. 연속 방정식이 두 개의 이차 보존량을 갖는다고 해서 이산화된 방정식이 그걸 물려받지는 않는다.

  • 중심차분 이류는 에너지를 그럭저럭 지켜도 엔스트로피를 지키지 못해, 작은 스케일에 에너지가 쌓이는 가짜 캐스케이드를 만든다. 격자 스케일에서 에너지가 축적되면 결국 발산한다.
  • 상류 차분으로 안정화하면 발산은 막히지만, 이번엔 수치 소산과 분산이 엔스트로피를 물리보다 훨씬 빠르게 갉아먹는다. 소용돌이가 며칠 만에 뭉개지는 대기 모델이 이렇게 태어난다.

이 딜레마의 고전적 해법이 아라카와 야코비안(Arakawa, 1966)이다. 2D 와도-유선함수 형식의 이류항 야코비안 J(ψ,ζ)J(\psi,\zeta)를 세 가지 방식으로 이산화하면 각각 다른 것을 보존하는데, 이들을 J=(J1+J2+J3)/3J = (J_1+J_2+J_3)/3로 평균하면 에너지와 엔스트로피를 동시에 이산 수준에서 보존한다. 그 결과 격자 스케일로의 가짜 에너지 유입이 구조적으로 차단되어, 인공 점성 없이도 장기 적분이 안정해진다. 오늘날 대기 모델의 역학코어가 “보존 성질”을 그렇게까지 따지는 관행의 출발점이 이 논문이다.3

관련 실무 지식 몇 가지.

  • 대와류 모사: 3D LES의 스마고린스키 모델은 부격자 항이 엔스트로피 소산을 대신한다. 그런데 이 모델을 그대로 2D에 가져다 쓰면 과도하게 소산적이라 역캐스케이드를 죽여버린다. 2D·준2D에서는 와도 기울기 기반의 리스(Leith) 점성 νe∝∣∇ζ∣\nu_e \propto |\nabla\zeta| 계열을 쓴다.
  • 초점성(hyperviscosity): ν∇2\nu\nabla^2 대신 (−1)p+1νp∇2p(-1)^{p+1}\nu_p\nabla^{2p}를 써서 큰 스케일은 건드리지 않고 격자 근처의 엔스트로피만 걷어내는 처방. 스펙트럴 코드의 국룰이다.
  • 진단량으로서의 엔스트로피: 계산이 이상해질 때 총 에너지·엔스트로피 시계열을 찍어 보면 원인이 대개 보인다. 에너지는 멀쩡한데 엔스트로피가 폭증하고 있다면 격자 스케일에서 뭔가 쌓이고 있다는 뜻이고, 격자를 조밀하게 하기 전에 도식을 먼저 의심해야 한다.

7. 관련 문서[편집]

8. Footnotes[편집]

  1. 12\frac{1}{2}가 붙은 정의는 에너지 12∫∣u∣2dV\frac{1}{2}\int|\mathbf{u}|^2dV와 짝을 맞추려는 것이고, 안 붙인 정의는 ε=ν⟨∣ω∣2⟩\varepsilon = \nu\langle|\boldsymbol{\omega}|^2\rangle 공식을 깔끔하게 쓰려는 것이다. 논문 사이에서 인수 2가 안 맞아 반나절을 태웠다면, 축하한다. 통과의례다. ↩

  2. 밀레니엄 문제 100만 달러가 걸린 3D 나비에-스토크스 정칙성 문제를 한 줄로 줄이면 “엔스트로피가 유한 시간에 터질 수 있는가”다. 수치적으로 후보 초기 조건을 잡아 폭발 직전까지 밀어붙이는 연구가 계속 나오지만, 아직 아무도 못 터뜨렸다(혹은 격자가 먼저 터졌다). ↩

  3. 1966년 논문 하나가 60년째 현역인 이유는, 그게 “정확도를 한 차수 올린 도식”이 아니라 “장기 적분이 애초에 가능하게 만든 도식”이었기 때문이다. 정확도는 격자로 살 수 있지만 보존성은 못 산다. ↩