[논문 리뷰] Convergence and error estimates of a viscosity-splitting finite-element scheme for the Primitive Equations

Francisco Guillén‐González, M. V. Redondo-Neble|arXiv (Cornell University)|2014. 11. 19.

Advanced Numerical Methods in Computational Mathematics참고 문헌 8인용 수 4

한 줄 요약

이 논문은 해양의 3차원 기본방정식에 대한 점성분리 유한요소 스킴을 제안하며, 약한 해로의 무조건적 안정성과 수렴성을 증명한다. $k \leq h^2$ 조건 하에 속도 및 압력에 대해 $O(k + h^l)$ 순서의 최적 오차 추정치를 확립한다. 여기서 $l=1$ 또는 $2$는 유한요소 근사 순서에 따라 달라지며, 시간분할 기법과 인프-서프 안정 유한요소를 조합하고, 코리olis 힘은 안정화 보정 전략을 통해 처리한다.

ABSTRACT

The purpose of this paper is the numerical analysis of a first order fractional-step time-scheme, using decomposition of theviscosity, and "inf-sup" stable finite element space-approximations for the Primitive Equations of the Ocean. The aim of the paper is twofold. Firstly, we prove that the scheme is unconditionally stable and convergent towards weak solutions of the Primitive Equations. Secondly, optimal error estimates for velocity and pressure are provided of order $O(k+h^l)$ for $l=1$ or $l=2$ when either first or second order finite-element approximations are considered ($k$ and $h$ being the time step and the mesh size, respectively). In both cases, these error estimates are obtained under the same constraint $k\le h^2$.

연구 동기 및 목표

• 해양 동역학의 3차원 기본방정식에 대해 안정적이고 수렴하는 수치적 스킴을 개발한다.
• 유한요소 근사에서 속도 및 압력에 대한 엄밀한 오차 추정치를 확립한다.
• 시간단계 제약 조건 $k \leq h^2$ 하에서 안정성과 수렴성을 확보한다.
• 수치적 구현 및 정확도 측면에서 비교 분석을 위해 문제의 두 가지 재구성 — 적분미분형 (Q) 및 완전 미분형 (R) — 을 고려한다.
• 지수적 시간 성장 방지를 위해 코리olis 항을 보정 전략을 통해 구현하고 안정성 추정치를 유지한다.

제안 방법

• 일阶 분할단계 시간 스킴을 사용하여 점성을 분리함으로써 운동량 방정식과 압력 방정식을 분리한다.
• 인프-서프 안정 유한요소 공간을 사용하여 공간 이산화를 수행하며, $h$-재분할을 통해 속도 및 압력 근사에 활용한다.
• 문제를 두 형태로 재구성한다: (Q) 적분미분계로, (R) 완전 미분계로, 둘 다 수직 메쉬 구조를 사용한다.
• 코리olis 항 $\mathbf{b}(\mathbf{u}) = f\mathbf{u}^\perp$ 은 중간 속도 단계에서 명시적으로 처리하고, 두 번째 부분단계에서 보정하여 안정성을 유지한다.
• 이sovotropic 추정 및 이산 그론발의 부등식을 사용하여 오차 전파를 통제하고 안정성 한계에서 지수적 성장 방지를 한다.
• 부분단계 2에서 수직 스토크스 단계를 해결하며, $A_h^{-1}$ 은 이산 역연산자를 나타내며, 보간 및 역연산 부등식을 통해 $W^{1,6}$-노름 추정치를 유도한다.

실험 결과

연구 질문

• RQ1 점성분리 유한요소 스킴이 3차원 기본방정식에 대해 무조건적 안정성과 수렴성을 달성할 수 있는가?
• RQ2 시간단계 $k$ 와 메쉬 크기 $h$ 에 대해 스킴의 최적 수렴 속도는 무엇인가?
• RQ3 적분미분형 (Q) 과 완전 미분형 (R) 재구성 간의 선택이 정확도 및 구현 복잡도에 어떤 영향을 미치는가?
• RQ4 코리olis 항이 안정성 추정치에서 지수적 시간 성장을 방지할 수 있는 방식으로 처리될 수 있는가?
• RQ5 $k \leq h^2$ 제약 조건이 스킴의 전체 정확도 및 실용성에 어떤 영향을 미치는가?

주요 결과

• 스킴은 약한 해로의 무조건적 안정성과 수렴성을 확보한다.
• 속도 및 압력에 대해 $O(k + h^l)$ 순서의 최적 오차 추정치를 달성하며, $l=1$은 선형 요소, $l=2$는 2차 요소에 해당한다.
• $k \leq h^2$ 조건은 두 재구성형 (Q) 및 (R) 모두에서 이러한 최적 오차 한계를 확보하기 위해 필수적이고 충분하다.
• (R)-스킴은 분석적으로 다소 열악하지만, 더 단순한 구현을 제공하며, $l=2$일 경우 속도 및 압력에 대해 $H^1 \times L^2$-노름에서 $O(k + h^2)$ 정확도를 달성한다.
• 코리olis 항은 지수적 시간 성장 방지를 위해 보정 전략을 통해 처리되어 최적 수렴성을 유지한다.
• 이산 역연산자 $A_h^{-1}$ 에 대한 이sovotropic $W^{1,6}$-노름 추정치가 유도되었으며, 보간 및 역연산 부등식을 통해 오차 분석을 지원한다.