[논문 리뷰] Discontinuous Galerkin methods for hyperbolic and advection-dominated problems on surfaces

Andreas Dedner, Pravin Madhavan|arXiv (Cornell University)|2015. 05. 25.

Advanced Numerical Methods in Computational Mathematics참고 문헌 19인용 수 5

한 줄 요약

이 논문은 음의 표면에서 이질성과 이류 지배 문제에 대한 불연속 갈레르킨(DG) 방법을 확장하며, 이산 표면에서 점프 안정화된 업윈드 플럭스를 사용하여 증명된 최적의 사전 오차 추정치를 제시한다. 수치적으로 Raviart-Thomas 유형의 속도 재구성으로 검증되었으며, 비정상적 진동이 없고 안정성을 보여준다.

ABSTRACT

We extend the discontinuous Galerkin (DG) framework to the analysis of first-order hyperbolic and advection-dominated problems posed on implicitly defined surfaces. The focus will be on the hyperbolic part, which is discretised using a "discrete surface" generalisation of the jump-stabilised upwind flux. A key issue arising in the analysis (which does not appear in the planar setting) is the treatment of the discrete velocity field, choices of which play an important role in the stability of the scheme. We then prove optimal error estimates in an appropriate norm given a number of assumptions on the discrete velocity field, which are then investigated and discussed in more detail. The theoretical results are verified numerically for a number of test problems exhibiting advection-dominated behaviour.

연구 동기 및 목표

• 암묵적으로 정의된 표면에서 1차 히퍼볼릭 및 이류 지배 문제에 대한 안정적인 불연속 갈레르킨 방법을 개발하기 위해.
• 표면이 곡선일 경우 표준 유한요소 방법이 비정상적 진동을 유도할 수 있는 이류 지배 흐름을 안정화하는 데 도전하는 문제를 다루기 위해.
• 평면에서의 점프 안정화된 업윈드 플럭스를 표면 기반 DG 설정으로 확장하기 위해.
• 이산 속도장의 근사가 방법의 안정성과 수렴성에 미치는 영향을 분석하기 위해.
• 예측 가능한 날카러운 층과 곡면 기하학을 포함한 벤치마크 문제에서 수치적으로 방법을 검증하기 위해.

제안 방법

• 배경 카르테시안 메쉬에서의 레벨셋 투영을 통해 연속 다각형 Γ의 이산 표면 근사를 구성한다.
• 표면 DG 설정에 점프 안정화된 업윈드 플럭스를 적용하여 히퍼볼릭 항을 다루며, Brezzi 등(2004)의 접근을 표면으로 일반화한다.
• 접선성과 발산 오차 감소를 위해 표면 메쉬에 Raviart-Thomas 유형의 보간을 사용해 이산 속도장을 도입한다.
• 특히 정규 성분과 발산에 관해 가정된 이산 속도장 조건 하에 L² 및 DG 노름에서 최적의 사전 오차 추정치를 유도한다.
• 라플라스-벨트라미 연산자를 위해 표면 유한요소 방법을 사용하고, 이와 함께 이류 항에 대해 DG를 결합한다.
• 이산 표면에서의 수치 적분을 수행하고, 요소 경계를 넘는 해와 플럭스의 점프를 고려한 일致한 이차형식을 구성한다.

실험 결과

연구 질문

• RQ1 어떻게 불연속 갈레르킨 방법을 암묵적으로 정의된 표면에서 이질성 및 이류 지배 문제에 확장할 수 있는가?
• RQ2 평면 영역과 비교해 곡면에서 이산 속도장이 초래하는 주요 안정성 도전 과제는 무엇인가?
• RQ3 이산 속도장에 대한 가정이 표면 DG 방법의 수렴성과 안정성에 어떤 영향을 미치는가?
• RQ4 속도장에 대해 Raviart-Thomas 유형의 보간이 비제로 발산 및 정규 성분 점프 문제를 완화시킬 수 있는가?
• RQ5 표면에서 이류 지배 문제에 대해 실질적으로 얻을 수 있는 수렴 속도와 정확도는 무엇인가?

주요 결과

• 점프 안정화된 업윈드 플럭스를 사용한 표면 DG 방법은 이류 지배 문제에서 L² 노름에서 최적의 O(h²) 수렴과 DG 노름에서 O(h) 수렴을 달성한다.
• 수치 실험 결과, 안정화되지 않은 표면 FEM은 전역적인 비정상적 진동을 보였지만, 표면 IP/UP 방법은 안정적이고 진동이 없는 것으로 나타났다.
• 이산 속도장에 Raviart-Thomas 유형의 보간을 사용하면 조각별 상수 투영 대비 오차를 약 0.5배로 줄일 수 있다.
• ε가 매우 작을 경우(예: ε = 10⁻⁶)에도 최적의 수렴 속도를 유지하여, 이류 지배 영역에서의 강인함을 확인했다.
• 이론적 오차 추정치는 수치적으로 검증되었으며, 적절한 속도장 가정 하에 수렴 차수는 예측된 점근적 비율에 가까워졌다.
• 이 방법은 단위 구와 같이 복잡한 기하학에서도 안정성과 정확도를 유지하며 효과적으로 작동한다.