[논문 리뷰] Velocity‑vorticity‑pressure formulation for the Oseen problem with variable viscosity

Verónica Anaya, Rubén Caraballo외 3명|arXiv (Cornell University)|2022. 04. 20.

Advanced Numerical Methods in Computational Mathematics참고 문헌 57인용 수 14

한 줄 요약

이 논문은 비정상 점성과 비선형 대류항을 고려한 Oseen 문제에 대해 속도-渇도-압력 혼합 유한요소 형식을 제안하며, 최소제곱 안정화 기법을 사용하여 비정상 점성과 대류항을 처리한다. 이 방법은 어떤 스토크스 호환 속도-압력 쌍과 임의의 순서의 난류 공간에 대해 inf-sup 안정성을 보장하며, 적응 메쉬 구분에서 최적 수렴률과 강인한 성능을 달성한다. 복잡한 기하구조와 변동 점성계수 분포를 포함한 수치 실험을 통해 검증되었다.

ABSTRACT

We propose and analyse an augmented mixed finite element method for the Oseen equations written in terms of velocity, vorticity, and pressure with non-constant viscosity and homogeneous Dirichlet boundary condition for the velocity. The weak formulation includes least-squares terms arising from the constitutive equation and from the incompressibility condition, and we show that it satisfies the hypotheses of the Babuška-Brezzi theory. Repeating the arguments of the continuous analysis, the stability and solvability of the discrete problem are established. The method is suited for any Stokes inf-sup stable finite element pair for velocity and pressure, while for vorticity any generic discrete space (of arbitrary order) can be used. A priori and a posteriori error estimates are derived using two specific families of discrete subspaces. Finally, we provide a set of numerical tests illustrating the behaviour of the scheme, verifying the theoretical convergence rates, and showing the performance of the adaptive algorithm guided by residual a posteriori error estimation

연구 동기 및 목표

• 비정상 점성으로 인해 기존의 난류 형식이 무너지는 상황에서, 비정상 점성계수를 가진 Oseen 문제에 대해 안정적이고 수렴하는 혼합 유한요소 방법을 개발한다.
• 이전에 Brinkman 및 스토크스 방정식에 사용된 증강 혼합 형식을, 대류항과 비정상 점성계수를 포함한 Oseen 문제로 확장한다.
• 점성계수가 공간적으로 변할 경우에도 Babuška-Brezzi 이론에 기반해 이산 문제의 안정성과 해의 존재성을 보장한다.
• 두 종류의 유한요소 하위공간에 대해 사전 및 사후 오차 추정식을 유도하여, 적응 메쉬 구분 전략을 가능하게 한다.

제안 방법

• 탄성 법칙과 비압축성 조건에서 유도된 최소제곱 항을 포함하여, 속도, 난류, 압력을 변수로 하는 Oseen 문제의 변분 형식을 정의한다.
• 강한 형태에서 유도된 안정화 항을 포함한 증강 혼합 유한요소 방법을 사용하여, 혼합 방법론의 Babuška-Brezzi 이론과의 호환성을 확보한다.
• 속도장의 컬을 통해 유도된 난류 방정식을 약한 형식으로 변환하며, 변형률 텐서 ε(u) = 1/2(∇u + ∇u^T)를 통해 비정상 점성계수에 적응한다.
• 가우스-레이스칼 최소제곱 접근법을 적용하여, 비정상 점성계수와 대류항 (β·∇)u로 인해 발생하는 도전 과제를 안정화한다.
• 잔차 기반 사후 오차 추정식을 사용하여 메쉬 구분을 유도하며, 재입구 모서리나 경계층과 같은 고오차 영역에 집중한다.
• 오차 지표에 따라 유도되는 블루-그린 메쉬 구분 전략을 적용하여, 특이점이나 급격한 기울기를 포함한 상황에서도 최적 수렴을 달성한다.

실험 결과

연구 질문

• RQ1 비정상 점성계수를 가진 Oseen 문제에 대해 속도, 난류, 압력을 주요 변수로 사용하는 안정적이고 수렴하는 혼합 유한요소 방법을 구성할 수 있는가?
• RQ2 비정상 점성계수에서는 항등식 curl(curl v) = −Δv + ∇(div v)가 직접 적용되지 않기 때문에, 기존의 난류 기반 형식을 어떻게 비정상 점성계수에 적응시킬 수 있는가?
• RQ3 점성계수가 공간적으로 변할 때와 대류항이 존재할 경우, inf-sup 안정성과 최적 수렴률을 보장하는 안정화 기법은 무엇인가?
• RQ4 사후 오차 추정이 특이 도메인(예: L자형 도메인)에서 최적 수렴률을 회복하기 위해 효과적인 적응 메쉬 구분을 이끌 수 있는가?
• RQ5 환자 맞춤형 abdominal arch 모델과 같은 실제 복잡한 기하구조에서, 비정상 점성계수 조건에서 이 방법은 어떻게 성능을 발휘하는가?

주요 결과

• 제안된 증강 혼합 유한요소 방법은 어떤 스토크스 안정 속도-압력 쌍과 임의의 순서의 난류 근사에도 대해 inf-sup 조건을 만족하며, 이산 문제의 유일한 해 존재성을 보장한다.
• 사전 오차 추정식은 모든 변수에서 최적 수렴률을 확인한다: 속도는 H1 노름, 난류는 L2 노름, 압력은 L2 노름 기준으로, 정렬 및 비정렬 메쉬에서 모두 검증되었다.
• 사후 오차 추정식의 효과성 지수는 약 1.12~1.16로 근접하여 신뢰할 수 있는 오차 추정과 효과적인 적응 메쉬 구분을 나타낸다.
• 재입구 모서리가 있는 L자형 도메인에서 적응 메쉬 구분이 최적 수렴률을 회복했으며, 균일 메쉬 구분은 특이성으로 인해 비최적 수렴을 보였다.
• 3D 대동맥 아치 시뮬레이션에서는 경계층과 난류 구조와 같은 복잡한 유동 특징을 잘 포착했으며, 적응 메쉬 구분은 내림내림 대동맥과 벽면 근처에 요소를 집중시키며 실제 혈역학 응용에서의 강인함을 확인했다.