A pressure-robust HHO method for the solution of the incompressible Navier-Stokes equations on general meshes

기술 요약: 일반 메쉬 상의 비압축성 나비에-스토크스 방정식을 위한 압력 강건한 HHO 방법

문제 정의
본 논문은 비압축성 나비에-스토크스 방정식의 수치 근사를 다루며, 특히 **압력 강건성(pressure robustness)**을 달리는 데 초점을 맞춥니다. 표준 이산화 방식에서는 속도장의 오차 추정치가 압력에 의존하게 되어, 체력(body force)에 큰 회전 성분 없는 성분(즉, 스칼라 포텐셜의 기울기)이 포함될 경우 성능이 저하되는 문제가 발생합니다. 표준 유한 요소 및 일치하는 심플리셜 메쉬(simplicial meshes)를 위한 압력 강건한 방법들은 개발되어 왔으나, 이를 일반 다면체 메쉬(general polytopal meshes)(임의의 다면체 요소 및 비일치 인터페이스를 지원하는 메쉬)로 확장하는 것은 여전히 중요한 과제로 남아 있습니다. 저자들은 이전 연구 [11]에서 제기된 열린 질문, 즉 일반 다면체 메쉬에 대해 압력 강건한 하이브리드 고차(H/H/O) 방법을 구축할 수 있는지에 대한 해답을 제시하고자 합니다.

방법론
제안된 방법은 2차원 및 3차원에서 유효한 완전한 압력 강건 HHO 스킴입니다. 핵심 혁신은 각 메쉬 요소 $T$의 서브트라이앵귤레이션(subtriangulation) 상에 정의된 국소 발산 보존 속도 재구성(local divergence-preserving velocity reconstruction) ($R^k_T$)의 구축에 있습니다.

1. 메쉬 및 서브메쉬: 본 방법은 다면체 메쉬 $\mathcal{M}_h$ 상에서 작동합니다. 결정적으로, 각 요소 $T$ 내에 일치하는 심플리셜 서브메쉬 $\mathfrak{T}_T$가 존재한다고 가정합니다. 모든 심플렉스가 공통 정점($x_T$)을 공유한다는 특정 기하학적 가정이 적용되며, 이를 통해 피라미드형 또는 비피라미드형 서브메쉬를 구성할 수 있습니다.
2. 발산 보존 재구성: 주어진 이산 속도 미지수 $\mathbf{v}_T$에 대하여, 재구성 $R^k_T \mathbf{v}_T$는 서브메쉬 $\mathfrak{T}_T$ 상의 국소 혼합 문제(local mixed problem)를 풀어 얻습니다. 이 문제는 $RT_k(\mathfrak{T}_T)$ 공간 내의 함수를 찾는 과정으로, 다음을 만족해야 합니다:
   • 요소 경계에서의 이산 속도의 법선 성분과 일치함.
   • 이산 발산을 보존함 (이산 발산 연산자 $D^k_T$와 교환 가능).
   • 코슐 공간(Koszul space) $G_{c, k-1}(T)$ 및 요소 내부와의 직교 조건을 만족함.
     이러한 구축 방식은 제한적인 면(face) 기하 구조를 요구하거나 기울기 유계(gradient bounds)에 관한 이론적 문제를 야기했던 기존 접근 방식(예: Wachspress 좌표계 의존성)의 한계를 극복합니다.
3. 이산 정식화:
   • 점성 항(Viscous Term): 서브메쉬 상의 표준 HHO 기울기 재구성을 사용합니다.
   • 체력(Body Force): 체력 $\mathbf{f}$는 가공되지 않은 이산 미지수가 아닌, 재구성된 속도 $R^k_h \mathbf{v}_h$를 테스트 함수로 사용하여 이산화됩니다. 이는 발산 보존 특성을 활용하여 회전 성분이 없는 힘(기울기)이 속도 해에 영향을 미치지 않도록 보장합니다.
   • 대류 항(Convective Term): 재구성된 속도 $R^k_h$와 확장된 기울기 재구성 $G^{2(k+1)}_{\mathfrak{T}_T}$를 사용하는 새로운 삼선형 형식을 도입합니다. 이 정식화는 비소산성(non-dissipativity)과 일관성을 보장하면서도 압력 강건성을 유지합니다.

주요 기여

• 일반 메쉬 지원: 본 논문은 특정 심플리셜 또는 특정 다면체 기하 구조에 국한되었던 이전의 제한 사항을 극복하고, 일반 다면체 요소로 구성된 메쉬에 대한 압력 강건한 HHO 방법의 확장에 대해 긍정적인 답변을 제공합니다.
• 새로운 재구성: 서브트라이앵귤레이션 상의 국소 혼합 문제에 기반한 발산 보존 속도 재구성을 도입했습니다. 이 접근 방식은 서브메쉬가 공통 정점 가정을 만족하는 경우, 임의의 다항식 차수 $k \ge 0$ 및 일반 다면체에 대해 적절히 정의됨을 입증하였습니다.
• 이론적 분석: 저자들은 엄격한 오차 추정치를 확립했습니다. 속도 오차가 압력 및 체력의 회전 성분 없는 부분(Hodge 분해의 파라미터 $\lambda$)으로부터 독립적임을 증명했습니다. 표준 정칙성 가정 하에, 에너지 노름에서의 수렴 속도는 $O(h^{k+1})$이며, 속도의 $L^2$ 노름에서는 $O(h^{k+2})$입니다.
• 강건성 증명: 이론적 결과는 제안된 방법이 거대한 회전 성분 없는 체력에 대해서도 강건함을 보여주며, 이는 표준 HHO 방식이 갖지 못한 특성입니다.

수치 결과
이론적 결과는 2차원에서의 종합적인 수치 테스트를 통해 검증되었습니다:

1. Kovasznay Flow: Cartesian, Hexagonal, Kershaw 메쉬에서 수렴 차수를 검증했습니다. 결과는 제안된 방법과 표준 HHO 모두에 대해 최적의 수렴 차수($h^{k+1}$ (에너지 노름), $h^{k+2}$ ($L^2$ 속도 오차))를 확인시켜 줍니다.
2. 강건성 테스트: 거대한 회전 성분 없는 체력($\lambda = 10^6$)을 가진 문제를 해결했습니다. 제안된 방법은 (k=1인 경우) 정확한 해를 재현하거나 $\lambda$의 크기에 관계없이 빠르게 수렴합니다. 반면, [8]의 표준 HHO 방법은 메쉬 세밀화에도 불구하고 줄어들지 않는 큰 오차를 보이며 압력 강건성이 결여되었음을 확인시켜 주었습니다.
3. Lid-Driven Cavity: 레이놀즈 수 $Re=1000$인 고전적인 lid-driven cavity 문제에 본 방법을 적용했습니다. 결과는 문헌의 참조 해와 일치합니다. 또한, 회전 성분 없는 체력이 추가되었을 때 제안된 방법의 속도 프로파일은 변하지 않지만, 표준 방식은 잘못된 다른 해로 수렴하는 것을 확인했습니다.

의의 및 주장
본 논문은 다면체 요소로 구성된 일반 메쉬 상에서 작동하는 완전한 압력 강건 HHO 방법을 제공한다고 주장합니다. 그 의의는 다면체 방법의 유연성(복잡한 기하 구조 및 비일치 인터페이스 처리 능력)과 압력 강건한 스킴의 물리적 충실도(회전 성분 없는 힘을 정확하게 처리하는 능력) 사이의 간극을 메우는 데 있습니다. 저자들은 서브트라이앵귤레이션 상의 국소 혼합 문제에 기반한 구축 방식이 기하학적 제한(예: 삼각형 면만 허용)이나 이전의 다면체 접근 방식에서 문제가 되었던 이론적 이슈(예: Wachspress 좌표계의 유계 문제) 없이 필요한 연속성 및 근사 유계(approximation bounds)를 도출할 수 있음을 강조합니다. 이 연구는 임의 차수의 정확도를 유지하면서 압력 강건한 이산화의 적용 범위를 더 넓은 클래스의 메쉬로 확장합니다.