[논문 리뷰] A viscosity-independent error estimate of a pressure-stabilized Lagrange-Galerkin scheme for the Oseen problem

Shinya Uchiumi|arXiv (Cornell University)|2017. 12. 12.

Advanced Numerical Methods in Computational Mathematics참고 문헌 42인용 수 10

한 줄 요약

이 논문은 동일한 순서의 $P_k/P_k$ 유한요소를 사용하는 일시적 Oseen 문제에 대해 압력 안정화된 라그랑주-갈레르킨 유한요소 스킴을 제안한다. 이 방법은 점성에 독립적인 오차 추정을 달성하며, 속도의 $L^2$-노름과 $ u^{1/2}$로 스케일링된 기울기의 $L^2$-노름에서, 점성 $ u$에 관계없이 $ u$-독립적인 수렴 속도 $ abla t + h^2 + h^k$를 확보한다. 수치 결과는 작은 $ u$와 큰 압력 조건에서 높은 정확도를 보이며, 타일러후드 요소를 능가한다.

ABSTRACT

We consider a pressure-stabilized Lagrange-Galerkin scheme for the transient Oseen problem with small viscosity. In the scheme we use the equal-order approximation of order $k$ for both the velocity and pressure, and add a symmetric pressure stabilization term. We show an error estimate for the velocity with a constant independent of the viscosity if the exact solution is sufficiently smooth. Numerical examples show high accuracy of the scheme for problems with small viscosity.

연구 동기 및 목표

• 작은 점성에서 안정적이고 고차수의 유한요소 스킴을 개발하기 위해.
• 대류 지배성과 압력 근사 오차로 인한 동일 순서 유한요소 방법의 불안정성 문제를 해결하기 위해.
• 속도와 그 기울기에 대해 점성 $ abla$에 독립적인 오차 추정을 확보하기 위해.
• 대칭 안정화를 사용하는 압력 안정화 라그랑주-갈레르킨 방법을 고차수 $P_k/P_k$ 요소로 확장하기 위해.
• 작은 $ abla$와 큰 압력 기울기 조건에서 스킴의 안정성과 정확도를 수치적으로 검증하기 위해.

제안 방법

• 속도와 압력에 대해 동일한 순서의 $P_k/P_k$ 유한요소($k \geq 1$)를 사용하여 압력 안정화 라그랑주-갈레르킨 스킴을 수립한다.
• Burman의 방법에 기반한 대칭 압력 안정화를 적용하여 동일 순서 근사의 안정성을 확보한다.
• 물질 도함수와 대류 지배 흐름을 다루기 위해 특성 기반 라그랑주-갈레르킨 방법과 시간 이산화된 변분 공식을 융합한다.
• 핵심 기법은 점성 $ abla$에 독립적인 오차 한계를 갖는 정확한 속도 해의 투영을 포함한다.
• 정밀한 오차 분해에 기반하여 기울기-안정화를 요구하지 않고도 분석을 수행한다.
• 시간 이산화에는 후진 오일러 스킴을 사용하며, 속도 장의 특성 기반 업데이트를 수행한다.

실험 결과

연구 질문

• RQ1 압력 안정화 라그랑주-갈레르킨 스킴이 $P_k/P_k$ 요소를 사용할 때 Oseen 문제에 대해 점성에 독립적인 오차 추정을 달성할 수 있는가?
• RQ2 대칭 압력 안정화 항은 인프-서프 안정성 조건을 요구하지 않고도 고차수 요소에 대해 안정성과 최적 수렴을 보장하는가?
• RQ3 작은 점성과 큰 압력 기울기 조건에서 이 스킴은 타일러후드 요소와 비교해 어떻게 성능을 발휘하는가?
• RQ4 속도와 그 기울기에 대한 오차 한계가 $L^2$-노름에서 점성 $ abla$에 독립적인가?
• RQ5 이 스킴은 비선형 문제에 대해 비슷한 안정성과 수렴 성질을 갖는 나비에-스토크스 문제로 확장 가능한가?

주요 결과

• 스킴은 속도의 $L^2$-노름과 속도 기울기의 $ u^{1/2}$-가중 $L^2$-노름에서 점성에 독립적인 오차 추정을 달성한다.
• 속도와 $ u^{1/2}$-가중 기울기의 수렴 속도는 $ abla t + h^2 + h^k$이며, 이는 점성 $ abla$에 관계없이 일정하다.
• 수치 실험 결과, $ abla = 10^{-4}$ 조건에서도 스킴이 높은 정확도를 유지하며, 큰 압력 기울기 조건에서 타일러후드 요소 스킴을 능가한다.
• 특히 작은 $ abla$ 조건에서 압력 안정화 스킴은 비안정화 스킴에 비해 속도 장의 진동을 크게 감소시킨다.
• 스킴은 본질적으로 조건부 안정이며, 시간 간격 $ abla t$와 메쉬 크기 $h$를 독립적으로 선택할 수 있다.
• 나비에-스토크스 문제에 대한 수치 결과는 스킴의 강건성과 비선형 문제에 대한 잠재적 적용 가능성을 보여준다.