[논문 리뷰] Enhancing stability of correction procedure via reconstruction using summation-by-parts operators II: Modal filtering

Jan Glaubitz, Hendrik Ranocha외 2명|arXiv (Cornell University)|2016. 06. 03.

Medical Imaging Techniques and Applications참고 문헌 12인용 수 5

한 줄 요약

이 논문은 초순수 보존법에 대한 수정 절차를 통한 복구(CPR) 방법의 안정성을 향상하기 위해 합성부분합(SBP) 프레임워크 내에서 모드 필터링을 도입함으로써, 시간에 대한 일계 정확도까지 모드 필터링과 스펙트럼 점성 간의 등가성을 확립하고, 적응형 필터 강도 전략을 제안하며, 명시적 오일러 시간 이산화를 사용하는 완전 이산화된 방법에 대해 증명 가능한 안정성을 입증한다. 이는 이전의 인위적 소산 방법에서 관찰된 시간 단계 제약을 극복한다.

ABSTRACT

A recently introduced framework of semidiscretisations for hyperbolic conservation laws known as correction procedure via reconstruction (CPR, also known as flux reconstruction) is considered in the extended setting of summation-by-parts (SBP) operators using simultaneous approximation terms (SATs). This reformulation can yield stable semidiscretisations for linear advection and Burgers' equation as model problems. In order to enhance these properties, modal filters are introduced to this framework. As a second part of a series, the results of Ranocha, Glaubitz, Öffner, and Sonar ("Enhancing stability of correction procedure via reconstruction using summation-by-parts operators I: Artificial dissipation", 2016) concerning artificial dissipation / spectral viscosity are extended, yielding fully discrete stable schemes. Additionally, a new adaptive strategy to compute the filter strength is introduced and different possible applications of modal filters are compared both theoretically and numerically.

연구 동기 및 목표

• SBP-SAT 프레임워크 내에서 모드 필터링을 적용하여 초순수 보존법에 대한 수정 절차를 통한 복구(CPR) 방법의 안정성을 향상시키는 것.
• 제1부에서 다룬 인위적 소산의 안정성 결과를 모드 필터링으로 확장하여 완전 이산화된 안정성을 보장하는 것.
• 사용자 입력을 줄이고 강건성을 향상시키기 위해 자동으로 필터 강도를 선택하는 새로운 적응형 전략을 개발하는 것.
• 다양한 모드 필터링 적용 전략(해의 필터링, 시간 도함수의 필터링, 연산자 분할 접근법)을 비교하고 분석하는 것.
• SBP-CPR 스킴의 맥락에서 모드 필터링과 스펙트럼 점성 간의 이론적 및 수치적 등가성을 확립하는 것.

제안 방법

• 선형 이동 및 버거스 방정식에 대해 보존적이고 안정적인 부분 이산화를 보장하기 위해, 동시 approximation term(SATs)을 갖는 합성부분합(SBP) 연산자를 사용하여 CPR 방법을 재구성한다.
• 해의 또는 시간 도함수에 대해 지수 필터를 통해 모드 필터링을 적용하며, 필터 행렬은 $\underline{\underline{F}}$로 표기된다.
• 연산자 분할 접근법을 사용하여 모드 필터링과 스펙트럼 점성 간의 등가성을 분석하며, 시간 도함수를 필터링하는 것이 $\underline{\underline{M}}\underline{\underline{F}}^{-1}$에 의해 유도되는 노름에서 안정적인 부분 이산화에 해당함을 보여준다.
• $(\Delta t)^2$ 오차 항을 보상하기 위해 새로운 적응형 필터 강도 전략을 도입하여, 명시적 오일러 시간 이산화를 사용할 때도 증명 가능한 안정성을 확보한다.
• 시간 이산화에는 명시적 오일러 스킴을 사용하며, 필터링된 시간 도함수는 $\partial_t \underline{u} \mapsto \partial_t \underline{\underline{F}} \underline{u}$로 적용되어 $\underline{u}_{+} = \underline{u} + \Delta t \, \partial_t \underline{\underline{F}} \underline{u}$로 이어진다.
• 수치 실험을 통해 해의 필터링, 시간 도함수의 필터링, 연산자 분할 접근법을 비교하며, 에너지 노름을 사용하여 안정성을 평가한다.

실험 결과

연구 질문

• RQ1 모드 필터링은 SBP-CPR 스킴에 효과적으로 적용되어 증명 가능한 안정성과 완전 이산화된 방법을 달성할 수 있는가?
• RQ2 안정성과 정확성 측면에서 모드 필터링은 인위적 소산 및 스펙트럼 점성과 어떻게 비교될 수 있는가?
• RQ3 해의 필터링과 시간 도함수의 필터링이 수치적 안정성과 진동 제어에 미치는 영향은 어떠한가?
• RQ4 시간 단계와 해의 행동에 따라 자동으로 필터 강도를 조정할 수 있는 적응형 필터 강도 전략을 설계할 수 있는가?
• RQ5 제안된 필터링 전략은 이전 인위적 소산 방법에서 관찰된 시간 단계 제약을 제거하는가?

주요 결과

• 모드 필터링은 시간에 대해 일계 정확도까지 스펙트럼 점성과 이론적으로 등가적이며, 인위적 소산의 안정적인 대체 방법을 제공한다.
• 적응형 필터 강도 전략은 $(\Delta t)^2$ 오차 항을 성공적으로 보상하여, 명시적 오일러 시간 이산화를 사용할 때도 증명 가능한 안정성과 완전 이산화된 스킴을 가능하게 한다.
• 시간 도함수를 필터링하는 것은 $\underline{\underline{M}}\underline{\underline{F}}^{-1}$에 의해 유도되는 노름에서 안정적인 부분 이산화를 이끌지만, 전체 시간 단계에서 역 필터 적용 시 이 안정성이 파괴되어, 해를 필터링하는 것보다 열등한 방법이 된다.
• 수치 실험 결과, 해를 필터링하는 것이 진동을 감소시키고 에너지 행동을 개선하는 반면, 시간 도함수를 필터링하는 것은 진동하는 에너지 노름을 유도함으로써 이론적 예측을 확인한다.
• 불연속적인 초기 자료를 가진 선형 이동 문제에서는 스킴이 안정적이고 정확하게 유지되지만, 충격 형성이 일어나는 비선형 케이스에서는 진동이 지속된다.
• 제안된 프레임워크는 SBP-CPR 스킴의 적용 가능성을 비선형 문제로 확장하며, 특히 적응형 필터링 메커니즘 덕분에 강건성이 향상된다.