회전·성층성 보수스키 흐름에서 에너지와 퍼텐셜 엔스트로피의 동시 다운스케일 전이

초록

본 연구는 고해상도 DNS 데이터를 이용해 회전·성층성 보수스키 방정식에 대한 코스그레이닝 기법을 적용, 에너지와 퍼텐셜 엔스트로피(잠재 엔스트로피)의 공간·스케일 동시 예산을 도출한다. 강한 회전·중간 성층, 중간 회전·강한 성층, 회전·성층이 동등하게 강한 세 경우에서 모두 큰‑스케일에서 작은‑스케일로의 일정한 전이(플럭스)가 관측되었으며, 이는 퍼텐셜 엔스트로피의 관성 구간 존재를 최초로 실증한다. 결과는 2차원 난류와는 다른 3차원 전이 메커니즘을 보여준다.

상세 분석

이 논문은 회전·성층성 보수스키 흐름에서 비선형 스케일 상호작용을 정량화하기 위해 ‘코스그레이닝(저역통과 필터)’ 접근법을 채택하였다. 필터링 연산 (\bar{a}\ell(\mathbf{x})=\int G\ell(\mathbf{r})a(\mathbf{x+r})d\mathbf{r}) 를 통해 물리량을 길이 스케일 (\ell) 보다 큰 부분만 남기고, 필터링되지 않은 작은 스케일의 효과는 ‘서브그리드 응력’ (\tau_\ell(f,g)=\overline{fg}-\bar f_\ell\bar g_\ell) 으로 나타낸다. 이 절차를 원래 보수스키 방정식(운동량 (1)·스칼라 (2))에 적용하면, 대규모 에너지와 퍼텐셜 엔스트로피의 보존 법칙에 대한 정확한 예산식이 얻어진다. 특히 에너지 예산식(7)에서 정의된 SGS 플럭스 (\Pi_\ell(\mathbf{x})=-\partial_j\bar u_i,\tau_\ell(u_i,u_j)-\partial_j\bar\theta,\tau_\ell(\theta,u_j)) 는 갈릴레오 불변성을 만족하고, (\ell) 가 최대 파수에 도달하면 자동으로 0이 되는 물리적 요구조건을 충족한다. 퍼텐셜 엔스트로피 예산식(11) 역시 동일한 구조를 가지며, 플럭스 (\Pi^Q_\ell(\mathbf{x})) 가 대규모에서 소규모로의 전이를 담당한다.

시뮬레이션은 640³ 격자, 8차 초점성 점성, 대규모 강제(k_f≈4) 조건에서 수행되었다. 세 가지 파라미터 조합(Rs, rS, RS)은 각각 f/N≫1, f/N≪1, f/N≈1을 의미하며, 모두 ‘선형 PV’ (q=f\partial_z\theta-N\omega_z) 가 전체 퍼텐셜 엔스트로피 (Q=\frac12 q^2) 의 99 % 이상을 차지하도록 설계되었다. 이는 KSW06이 제시한 ‘잠재 엔스트로피가 2차 형태로 근사되는’ 조건과 일치한다.

결과적으로, 모든 경우에서 (\langle\Pi_\ell\rangle)와 (\langle\Pi^Q_\ell\rangle)가 (\ell) 범위(≈ k = 6 ~ 30)에서 거의 일정하고 양의 값을 보였다. 이는 에너지와 퍼텐셜 엔스트로피가 동시에 다운스케일 전이를 겪으며, 각각의 관성 구간이 존재함을 의미한다. 특히 퍼텐셜 엔스트로피 플럭스가 일정하다는 사실은 KSW06이 예측한 ‘2/3 법칙’의 실증이며, 기존 연구(예: Herring et al.)가 제시한 ‘퍼텐셜 엔스트로피는 전 스케일에 걸쳐 점성에 의해 소멸한다’는 주장과는 대조적이다. 또한, 에너지 플럭스가 일정함은 2차원 난류에서와 달리 3차원 회전·성층성 흐름이 양방향 전이가 아닌 순수한 전방 전이를 보인다는 중요한 물리적 의미를 갖는다.

이러한 발견은 (1) 회전·성층성 대기·해양 흐름의 스케일 분광 모델링에 새로운 제약조건을 제공하고, (2) LES와 같은 대규모 시뮬레이션에서 서브그리드 모델링을 설계할 때 에너지와 퍼텐셜 엔스트로피의 동시 전이를 고려해야 함을 시사한다. 또한, 고차 점성에 의한 인위적 열화가 관성 구간에 미치는 영향을 최소화했음에도 불구하고 플럭스가 견고하게 유지된 점은 결과의 신뢰성을 높인다.