[논문 리뷰] Finite element methods respecting the discrete maximum principle for convection-diffusion equations

Gabriel R. Barrenechea, Volker John외 1명|arXiv (Cornell University)|2022. 04. 15.

Advanced Numerical Methods in Computational Mathematics인용 수 6

한 줄 요약

이 논문은 비선형성 우세 영역에서의 대류-확산-반응 방정식에 대해 이산 최대원칙(DMP)을 만족하는 유한요소법을 조사한다. 여기서는 대칭화된 안정화 방법과 FEM-FCT, AFC와 같은 비선형 방법이 현재까지 유일하게 DMP를 유지하면서도 정확한 해를 제공하는 방법임을 밝히며, 단조성과 고차원 이산화의 최근 발전을 강조한다.

ABSTRACT

Convection-diffusion-reaction equations model the conservation of scalar quantities. From the analytic point of view, solution of these equations satisfy under certain conditions maximum principles, which represent physical bounds of the solution. That the same bounds are respected by numerical approximations of the solution is often of utmost importance in practice. The mathematical formulation of this property, which contributes to the physical consistency of a method, is called Discrete Maximum Principle (DMP). In many applications, convection dominates diffusion by several orders of magnitude. It is well known that standard discretizations typically do not satisfy the DMP in this convection-dominated regime. In fact, in this case, it turns out to be a challenging problem to construct discretizations that, on the one hand, respect the DMP and, on the other hand, compute accurate solutions. This paper presents a survey on finite element methods, with a main focus on the convection-dominated regime, that satisfy a local or a global DMP. The concepts of the underlying numerical analysis are discussed. The survey reveals that for the steady-state problem there are only a few discretizations, all of them nonlinear, that at the same time satisfy the DMP and compute reasonably accurate solutions, e.g., algebraically stabilized schemes. Moreover, most of these discretizations have been developed in recent years, showing the enormous progress that has been achieved lately. Methods based on algebraic stabilization, nonlinear and linear ones, are currently as well the only finite element methods that combine the satisfaction of the global DMP and accurate numerical results for the evolutionary equations in the convection-dominated situation.

연구 동기 및 목표

• 대류 우세 대류-확산 문제에서 이산 최대원칙(DMP)을 엄격히 만족하는 유한요소법을 특정하고 분석하는 것.
• 안정화된 유한요소법에서 DMP 만족도와 수치 정확도 사이의 상충 관계를 평가하는 것.
• 선형과 비선형 안정화 기법이 해의 물리적 범위를 유지하는 데 어떻게 기여하는지 평가하는 것.
• 정상 상태 및 시간에 의존하는 문제에 대해 대칭화된 안정화 및 유량 보정 유한요소법의 최근 진전을 요약하는 것.
• 대류 우세 흐름에서 정규성 보존과 비물리적 진동을 방지하는 데 기여하는 방법들을 종합적으로 개괄하는 것.

제안 방법

• 대칭화된 안정화를 통해 DMP를 강제하는 유한요소법에 중점을 두며, 인공 확산 및 국소 투영 기법을 포함한다.
• 정상 상태 문제에 대해 Mizukami–Hughes 방법, Burman–Ern 방법, 그리고 대칭화된 유량 보정(AFC)과 같은 비선형 방법을 분석한다.
• 시간에 의존하는 문제에 대해 FEM-FCT와 단조성 보장 국소 투영 안정화(LPS) 방법을 검토하며, DMP 유지 능력을 강조한다.
• 상향 스무딩 및 가장자리 평균 유한요소와 같은 선형 방법을 검토하며, 대류 우세 조건에서 DMP를 만족하는 데 한계가 있음을 지적한다.
• 한계기법과 기울기 근사치가 단조성과 선형성 보존을 유지하면서 정확도를 향상시키는 데 수행하는 역할을 평가한다.
• 다양한 유한요소 유형 간의 방법 비교를 수행하며, $\mathbb{Q}_1$, 고차원 $H^1$-일치, Crouzeix–Raviart, 그리고 비연속 갈레르킨 요소를 포함한다.

실험 결과

연구 질문

• RQ1 대류 우세 영역에서 대류-확산 방정식에 대해 이산 최대원칙(DMP)을 엄격히 만족하는 유한요소법은 무엇인가?
• RQ2 대칭화된 안정화 방법은 전통적인 안정화 방법에 비해 DMP 만족도와 해 정확도 측면에서 어떻게 비교되는가?
• RQ3 유량 보정 및 국소 투영과 같은 비선형 안정화 메커니즘이 단조성 유지와 진동 방지를 위해 수행하는 역할은 무엇인가?
• RQ4 선형 유한요소 방법은 고차원 정확도를 유지하면서도 전역 DMP를 만족하도록 수정될 수 있는가?
• RQ5 3차원에서 DMP를 유지하는 방법을 구현할 경우 실용적 상충 관계는 무엇이며, 특히 계산 비용과 비선형 솔버 효율성에 미치는 영향은 어떠한가?

주요 결과

• 정상 상태 방법 중에서 유일하게 비선형 방법—특히 대칭화된 안정화 및 유량 보정 방법—이 DMP를 만족하면서도 정확한 해를 제공한다.
• Mizukami–Hughes 방법은 유일하게 DMP를 만족함을 증명한 전통적 안정화 방법이지만, 대부분의 다른 방법들은 대류 우세 조건에서 실패한다.
• 최근 개발된 대칭화된 안정화 방법, 예를 들어 수정된 LPS와 기울기 기반 한계기법은 더 높은 정확도를 확보하면서도 과도한 퍼짐 없이 테라스 현상을 감소시킨다.
• 시간에 의존하는 문제에 대해 현재까지 DMP를 전역적으로 유지하면서도 양호한 정확도를 유지하는 오직 하나의 방법은 대칭화된 안정화 FEM-FCT 및 AFC 방법뿐이다.
• 대칭화된 안정화의 선형 변형은 정확도와 계산 효율성 사이의 유리한 균형을 제공하므로 실용적 응용에 적합하다.
• 반연속적 스킴에 단조성 보장 국소 투영 연산자를 추가하면 테라스 현상이 크게 감소하고 해의 품질이 향상되며, DMP 유지에 영향을 주지 않는다.