[논문 리뷰] Metastability and rapid convergence to quasi-stationary bar states for the 2D Navier-Stokes Equations

Margaret Beck, C. Eugene Wayne|arXiv (Cornell University)|2011. 08. 17.

Fluid Dynamics and Turbulent Flows참고 문헌 17인용 수 4

한 줄 요약

이 논문은 토러스 위 2차원 나비에-스토크스 방정식에서 바 상태(bar states)의 준안정성에 대해 동역학계 이론을 제시한다. 바 상태 주변의 선형화된 연산자를 근사하고, 속도가 $\mathcal{O}(e^{-\sqrt{\nu}t})$ 속도로 급격히 감쇠됨을 보여줌으로써, 느린 점성 감쇠($\mathcal{O}(e^{-\nu t})$)에도 불구하고 이러한 상태가 오랫동안 지속되는 이유를 설명한다. 수치적 증거는 전체 선형 연산자에 대해 이 감쇠 속도가 최적임을 확인한다.

ABSTRACT

Quasi-stationary, or metastable, states play an important role in two-dimensional turbulent fluid flows where they often emerge on time-scales much shorter than the viscous time scale, and then dominate the dynamics for very long time intervals. In this paper we propose a dynamical systems explanation of the metastability of an explicit family of solutions, referred to as bar states, of the two-dimensional incompressible Navier-Stokes equation on the torus. These states are physically relevant because they are associated with certain maximum entropy solutions of the Euler equations, and they have been observed as one type of metastable state in numerical studies of two-dimensional turbulence. For small viscosity (high Reynolds number), these states are quasi-stationary in the sense that they decay on the slow, viscous timescale. Linearization about these states leads to a time-dependent operator. We show that if we approximate this operator by dropping a higher-order, non-local term, it produces a decay rate much faster than the viscous decay rate. We also provide numerical evidence that the same result holds for the full linear operator, and that our theoretical results give the optimal decay rate in this setting.

연구 동기 및 목표

• 고레이놀즈 수에서 2차원 비압축성 나비에-스토크스 유동에서 바 상태의 장수 및 준정적 행동을 설명하기 위해.
• 바 상태 주변의 선형화된 동역학을 분석하고, 이러한 상태로의 빠른 수렴 메커니즘을 규명하기 위해.
• 선형화된 연산자의 근사가 점성 시간 스케일보다 훨씬 빠른 감쇠 속도를 제공함을 보여주기 위해.
• 이론적 감쇠 속도가 전체 선형화된 연산자에 대해 최적임을 수치적으로 증명하기 위해.
• 비선형 나비에-스토크스 방정식으로의 선형 분석 확장을 위한 기초를 마련하기 위해.

제안 방법

• 고차항이면서 국소적이지 않은 항을 생략함으로써 바 상태 주변의 시간에 의존하는 선형화된 연산자를 근사하여, 분석 가능한 스펙트럼 성질을 갖는 단순화된 연산자를 도출한다.
• 스펙트럼 분석을 통해 단순화된 연산자가 $\mathcal{O}(e^{-\sqrt{\nu}t})$ 속도로 감쇠됨을 보이며, 이는 점성 감쇠 속도 $\mathcal{O}(e^{-\nu t})$보다 훨씬 빠르다는 것을 입증한다.
• 변수변환 $w = \sqrt{1 + \Delta^{-1}} u$를 적용하여 선형 연산자를 비대칭 성분을 갖는 형태로 변환함으로써 스펙트럼 분석을 용이하게 한다.
• 전체 선형화된 연산자의 고유값을 수치적으로 계산하여 근사로부터 도출된 이론적 감쇠 속도와 비교한다.
• 푸리에 분석과 비에트-사바르 법칙을 사용하여 코시누스 모드를 통해 회전율과 속도를 연결하고, 바 상태 해를 명시적으로 정의한다.
• 바 상태가 올리어 방정식의 최대 엔트로피 해임을 이용하여, 2차원 난류에서 물리적으로 관련된 장수 상태와 연결한다.

실험 결과

연구 질문

• RQ1 왜 바 상태는 점성 시간 스케일에서는 불안정함에도 불구하고 2차원 나비에-스토크스 유동에서 오랫동안 지속되는가?
• RQ2 동역학계 접근법이 수치 시뮬레이션에서 관측된 바 상태로의 빠른 수렴을 설명할 수 있는가?
• RQ3 선형화된 근사에서의 감쇠 속도 $\mathcal{O}(e^{-\sqrt{\nu}t})$는 전체 선형화된 연산자에 대해 최적인가?
• RQ4 비자기적 성질(non-self-adjointness)은 선형화된 동역학에서 빠른 감쇠 모드를 가능하게 하는 역할을 하는가?
• RQ5 바 상태에 대한 분석 프레임워크를 비선형 나비에-스토크스 방정식으로 확장할 수 있는가?

주요 결과

• 바 상태 주변의 단순화된 선형화된 연산자는 $\mathcal{O}(e^{-\sqrt{\nu}t})$ 속도로 감쇠됨을 보이며, 이는 점성 감쇠 속도 $\mathcal{O}(e^{-\nu t})$보다 훨씬 빠르다.
• 전체 선형화된 연산자의 고유값을 수치적으로 계산한 결과, 감쇠 속도 $\mathcal{O}(e^{-\sqrt{\nu}t})$가 근사에 국한되지 않고 실제 시스템에 대해서도 최적임을 확인하였다.
• 근사된 선형 연산자에 대한 불변 부분공간은 해가 바 상태로 빠르게 수렴하도록 보장하며, 준안정성의 메커니즘을 제공한다.
• 이 빠른 감쇠의 메커니즘은 버거스 방정식의 경우와 다름을 보이며, 이는 큰 계수의 고유함수 전개에서 비롯되지 않고 $\sqrt{\nu}$ 비례하는 작은 고유값에서 기인한다.
• 선형화된 연산자는 비자기적임을 보이며, 이는 빠른 감쇠 모드 존재와 시간 스케일 분리의 근본 원인이다.
• 결과는 비선형 나비에-스토크스 방정식으로의 분석 확장을 위한 길을 제시하지만, 선형화된 연산자가 시간에 따라 변하기 때문에 표준 스펙트럼 투영 기법을 적용하는 것이 복잡하다.