[논문 리뷰] On Euler's equation and `EPDiff'

David Mumford, Peter W. Michor|arXiv (Cornell University)|2012. 09. 28.

Nonlinear Waves and Solitons참고 문헌 6인용 수 3

한 줄 요약

이 논문은 비압축성 비점성 유체 흐름에 대한 오일러 방정식의 해가 매개변수 $ \varepsilon, \eta \geq 0 $를 가진 정규화 연산자 $ L_{\varepsilon,\eta} $를 사용하여 미분형식군 위의 지오데식 방정식 가족인 EPDiff(에일러-포이아르 방정식)의 해의 극한임을 규명한다. 핵심 결과는 EPDiff 해와 오일러 방정식 사이의 국소 시간 내 수렴 바ounds이며, 유동에 따른 운동량 보존이 안정성과 컴팩트 지지의 전파를 보장한다.

ABSTRACT

We study a family of approximations to Euler's equation depending on two parameters $\varepsilon,η\ge 0$. When $\varepsilon=η=0$ we have Euler's equation and when both are positive we have instances of the class of integro-differential equations called EPDiff in imaging science. These are all geodesic equations on either the full diffeomorphism group $\operatorname{Diff}_{H^\infty}(\mathbb R^n)$ or, if $\varepsilon = 0$, its volume preserving subgroup. They are defined by the right invariant metric induced by the norm on vector fields given by $$ \|v\|_{\varepsilon,η} = \int_{\mathbb R^n} dx $$ where $L_{\varepsilon,η} = (I- frac{η^2}{p} riangle)^p \circ (I- frac1{\varepsilon^2} abla \circ ÷)$. All geodesic equations are locally well-posed, and the $L_{\varepsilon,η}$-equation admits solutions for all time if $η>0$ and $p\ge (n+3)/2$. We tie together solutions of all these equations by estimates which, however, are only local in time. This approach leads to a new notion of momentum which is transported by the flow and serves as a generalization of vorticity. We also discuss how delta distribution momenta lead to "vortex-solitons", also called "landmarks" in imaging science, and to new numeric approximations to fluids.

연구 동기 및 목표

• 오일러 방정식과 미분형식군 위의 지오데식 방정식 가족 간의 엄밀한 연결 고리를 확립하기 위해.
• 유체 흐름 방정식에 대한 연산자 $ L_{\varepsilon,\eta} = (I - \frac{\eta^2}{p}\Delta)^p \circ (I - \frac{1}{\varepsilon^2}\nabla\circ\text{div}) $ 의 정규화 효과를 분석하기 위해.
• EPDiff 해가 $ \eta > 0 $ 이고 $ p \geq (n+3)/2 $ 라면 전역적으로 잘 정의되어 있으며, 오일러 방정식의 근사화를 가능하게 함을 보여주기 위해.
• 유동에 의해 운반되는 일반화된 운동량 개념을 도입하여 비선형성의 개념을 확장하기 위해.
• 델타-분포 운동량에서 유래된 비틀림-솔리톤(랜드마크)의 출현과 수치적 유체 근사에서의 역할을 탐색하기 위해.

제안 방법

• 무한 차수의 미분형식군 $ \text{Diff}_{H^\infty}(\mathbb{R}^n) $ 에서 오른쪽 불변 계량을 정의하는 데 사용된 노름 $ \|v\|_{\varepsilon,\eta}^2 = \int \langle L_{\varepsilon,\eta}v, v \rangle \, dx $ 를 사용하며, 여기서 $ L_{\varepsilon,\eta} $ 는 스무딩과 발산 보정을 통합한다.
• 지오데식 방정식으로서의 EPDiff 방정식 $ \partial_t m = -(v\cdot\nabla)m - \text{div}(v)m - m\cdot(Dv)^T $ 의 유도를 통해, 여기서 $ m = L_{\varepsilon,\eta}v $ 이다.
• 유동 $ \varphi $ 에 대해 운동량 1-form $ \widetilde{m} $ 이 불변임을 증명한다. 즉, $ \widetilde{m}(\cdot,t) = \varphi(\cdot,t)_* \widetilde{m}(\cdot,0) $ 이며, 이는 지지의 보존과 감쇠 성질의 유지를 의미한다.
• EPDiff 해가 $ \varepsilon, \eta > 0 $ 인 경우와 오일러 방정식 해 사이의 국소 시간 내 추정을 수립함으로써, $ \varepsilon, \eta \to 0 $ 일 때의 수렴성을 확립한다.
• 유동을 초기 운동량으로 표현하기 위해 라그랑주 형식 $ \partial_t \varphi = K^{\varphi} \ast (\varphi_* \widetilde{m}(\cdot,0)) $ 를 사용하며, 여기서 $ K $ 는 $ L_{\varepsilon,\eta} $ 의 그린 함수이다.
• $ \delta $-함수 운동량을 통한 랜드마크 지오데식의 수치적 탐색을 통해, 각운동량과 에너지에 따라 상호작용이 반발력 또는 흡인력으로 나타남을 보였다.

실험 결과

연구 질문

• RQ1 정규화 매개변수 $ \varepsilon, \eta > 0 $ 를 가진 EPDiff 방정식의 해가 $ \varepsilon, \eta \to 0 $ 일 때 오일러 방정식의 해로 수렴하는 방식은 어떠한가?
• RQ2 일반화된 운동량 $ m = L_{\varepsilon,\eta}v $ 는 유동 하에서 어떤 구조를 유지하는 데 기여하며, 이는 비선형성의 개념을 어떻게 일반화하는가?
• RQ3 EPDiff 방정식이 전역 해를 갖는 조건은 무엇이며, 이는 오일러 방정식의 잘 정의됨과 어떻게 관련되는가?
• RQ4 델타-분포 운동량은 어떻게 비틀림-솔리톤(랜드마크)을 유도하며, 어떤 역학적 행동을 나타내는가?
• RQ5 EPDiff 프레임워크는 오일러 방정식의 안정적인 수치 근사로 기능할 수 있는가, 특히 $ \eta \to 0 $ 일 때?

주요 결과

• 오일러 방정식의 해는 $ L_{\varepsilon,\eta} $ 를 사용하여 $ \varepsilon, \eta \to 0 $ 일 때 EPDiff 방정식의 해의 극한으로 나타나며, 국소 시간 내 수렴 바ounds 가 확립되었다.
• $ \eta > 0 $ 이고 $ p \geq (n+3)/2 $ 인 경우, EPDiff 방정식은 전역적으로 잘 정의되어 있으며, 모든 시간 동안 유일한 해를 가지며 정규화된 근사화를 보장한다.
• 운동량 1-form $ \widetilde{m} $ 는 유동 $ \varphi $ 에 대해 불변이므로 컴팩트 지지와 급격한 감쇠 성질이 유지되며, 이는 구조적 안정성을 의미한다.
• 델타-분포 운동량에서 유래된 비틀림-솔리톤(랜드마크)는 에너지와 각운동량에 따라 산산이 흩어지거나 포착되는 역학적 행동을 보이며, 임계 에너지 한계는 $ E = \frac{8}{5}|\omega|^2 $ 에서 발생한다.
• 핵심 함수 $ K_{0,1} $ 의 도함수 근처에서 유도된 벡터장은 $ (x_1,x_2) $-평면에서 나선형 유선을 보이며 $ x_3 $-축을 따라 외부 유동을 나타내어 국소적인 비틀림-링 모양의 행동을 보인다.
• 랜드마크 지오데식의 수치적 통합 결과, 배경 운동량 $ \overline{m} $ 의 변화에 따라 반발력과 흡인력 상호작용이 모두 나타나며, 이는 유체 불안정성의 모델링에 유용함을 시사한다.