[논문 리뷰] An unfitted discontinuous Galerkin scheme for conservation laws on evolving surfaces

Christian Engwer, Thomas Ranner외 1명|arXiv (Cornell University)|2016. 02. 02.

Advanced Mathematical Modeling in Engineering참고 문헌 6인용 수 3

한 줄 요약

이 논문은 표면에 맞는 메시가 필요 없이 변화하는 표면에서 보존법칙을 해결하기 위한 비일치 불연속 갈레르킨(DG) 방법을 제시한다. 표면은 레벨셋 표현을 사용하고, 고정된 배경 메시에서 표준 이동 문제를 해결함으로써, 복잡한 기하구조에서도 질량 보존과 최적 수렴률(1차 수렴률: L¹ 및 L²)을 달성한다. 수치 결과는 수축하는 원형 표면에서의 안정성과 정확성을 확인한다.

ABSTRACT

Motivated by considering partial differential equations arising from conservation laws posed on evolving surfaces, a new numerical method for an advection problem is developed and simple numerical tests are performed. The method is based on an unfitted discontinuous Galerkin approach where the surface is not explicitly tracked by the mesh which means the method is extremely flexible with respect to geometry. Furthermore, the discontinuous Galerkin approach is well-suited to capture the advection driven by the evolution of the surface without the need for a space-time formulation, back-tracking trajectories or streamline diffusion. The method is illustrated by a one-dimensional example and numerical results are presented that show good convergence properties for a simple test problem.

연구 동기 및 목표

• 표면에 맞는 메시가 필요 없이 시간에 따라 변화하는 표면에서 광범위하게 이동하는 보존법칙을 해결하기 위한 강력한 수치적 방법을 개발하기 위해.
• 공간-시간 수식, 역경로 추적, 스트림라인 확산을 피하기 위해 문제를 비일치 영역에서의 표준 이동 방정식으로 재구성함으로써, 복잡한 수식을 피하기 위해.
• 수치적 스킴에서 기계 정밀도 수준까지 이산 질량 보존을 보장하기 위해.
• 특히 표면 운동에 의해 이동하는 문제에 대해, 변화하는 표면을 포함한 기준 문제에서 수렴성과 정확성을 입증하기 위해.
• 미래의 확장에 기초를 마련하기 위해, 시간에 따라 변화하는 문제와 비영인 플럭스 항(q ≠ 0)을 포함한 문제에 적용 가능하도록 하기 위해.

제안 방법

• 표면는 레벨셋 함수 Φ(x,t)를 통해 암묵적으로 표현되며, 변화하는 표면에 맞는 메시를 맞출 필요 없이, 비일치 불연속 갈레르킨 수식을 적용한다.
• 표면의 진화는 미분형 매개변수화 G(·,t)의 시간 도함수를 통해 정의된 속도장 w를 통해 기술되며, 표면 운동을 기술할 수 있다.
• 기본 방정식은 물질 도함수 문제로 재구성된다: ∂•u + u∇Γ·w + ∇Γ·q = 0이며, 제시된 경우 q = 0이므로 표면 운동에 의한 순수 이동 문제로 축소된다.
• 고정된 배경 메시에서 불연속 갈레르킨 변분 수식을 적용하며, 요소 간의 불연속성을 다루고 안정성을 확보하기 위해 요소 경계에서 수치적 플럭스를 계산한다.
• 레벨셋 함수를 통한 표면 측도의 일관성 있는 통합을 사용하여, 재메시 없이 잘린 요소에서 정확한 적분을 가능하게 한다.
• 스킴은 다항식 기반 형상 함수(k=0이 보고된 테스트에서)를 사용하며, 요소 경계를 넘는 해의 연속성은 약한 방식으로 강제하여 질량 보존을 확보한다.

실험 결과

연구 질문

• RQ1 비일치 불연속 갈레르킨 방법이 표면에 맞는 메시가 필요 없이 변화하는 표면에서 안정적이고 정확한 해를 달성할 수 있는가?
• RQ2 표면이 배경 메시를 가로질러 변화하더라도, 이산 설정에서 기계 정밀도 수준까지 질량을 보존하는가?
• RQ3 수축하는 원형 표면에서 모델 이동 문제에 대해 L¹, L², L∞ 노름에서 달성할 수 있는 수렴률은 무엇인가?
• RQ4 비연속 초기 자료에 대해, 특히 수치적 확산과 불연속성의 날카움 측면에서 방법의 성능은 어떠한가?
• RQ5 미래 작업에서 시간에 따라 변화하는 문제와 비영인 플럭스 항(q ≠ 0)을 포함한 방정식으로 확장할 수 있는가?

주요 결과

• 메시 세분화(h-refinement) 조건에서 L¹ 및 L² 노름에서 약 1의 최적 수렴 순서를 달성하며, 충분히 미세한 메시에서 수렴률이 1에 가까워진다.
• L∞ 노름에서는 굵은 해상도에서 수렴률이 일관되지 않지만, h < 0.01일 경우 순서 1에 수렴하는 것으로 나타나, 메시 세분화에 따라 정확도가 향상됨을 시사한다.
• 이산 질량은 기계 정밀도 수준에서 보존된다: 모든 세분화에서 총 질량은 6.283185 ± 1e-6 범위에 있으며, 해석적 값 2π와 매우 가까이 일치한다.
• 수축하는 원형 표면에서 일정한 초기 농도(u₀ ≡ 1)가 최종 농도(u ≡ 2)로 이동하는 이동 현상을 효과적으로 포착하며, 수치적 확산이 최소화된다.
• 급격한 점프가 있는 비정상적인 초기 자료(이진 분포)에 대해서도 불연속성의 날카움을 상당히 잘 유지하지만, 일부 수치적 확산이 관찰되며, 이는 고차원 요소와 플럭스 리미터를 사용하여 완화될 수 있다.
• 표면이 배경 메시를 가로질러도 방법은 안정적이고 정확하며, 재메시 없이도 복잡한 변화하는 기하구조를 효과적으로 처리함을 보여주며, 강건성을 입증한다.