[논문 리뷰] Residual equilibrium schemes for time dependent partial differential equations

Lorenzo Pareschi, Thomas Rey|arXiv (Cornell University)|2016. 02. 08.

Computational Fluid Dynamics and AerodynamicsEngineering참고 문헌 38인용 수 19

한 줄 요약

이 논문은 시간에 의존하는 PDE에 대해 평형 상태 해를 유지하면서 고차 정확도를 달성하는 잔여 평형 계획을 소개한다. 평형 기반 보정과 유량 제한 기법을 사용해 수치적 유량을 수정함으로써, 이는 대류-확산, 쌍곡 평형 법칙, 그리고 운동 방정식에 대해 잘 균형 잡힌 행동을 보장하며, 정확도를 희생시키지 않은 채 안정성과 단조성을 유지한다.

ABSTRACT

Many applications involve partial differential equations which admits nontrivial steady state solutions. The design of schemes which are able to describe correctly these equilibrium states may be challenging for numerical methods, in particular for high order ones. In this paper, inspired by micro-macro decomposition methods for kinetic equations, we present a class of schemes which are capable to preserve the steady state solution and achieve high order accuracy for a class of time dependent partial differential equations including nonlinear diffusion equations and kinetic equations. Extension to systems of conservation laws with source terms are also discussed.

연구 동기 및 목표

• 시간에 의존하는 PDE에서 표준 방법이 소스항 균형으로 인해 실패할 때 비트레이스러운 평형 상태 해를 유지하는 수치적 계획을 개발하는 것.
• 운동 방정식에서의 미세-거시 분해 기법을 비선형 확산과 쌍곡 평형 법칙을 포함한 더 넓은 PDE 클래스로 확장하는 것.
• 소스항이 존재하는 상황에서 고차 정확도를 유지하면서 잘 균형 잡힌 및 TVD(Total Variation Diminishing) 성질을 확보하는 것.
• 특정 문제에 맞는 유량 설계 없이도 평형 상태 유지가 가능한 구조적으로 단순하고 광범위하게 적용 가능한 프레임워크를 제공하는 것.

제안 방법

• 수치적 유량이 평형 상태 해로부터 유도된 평형 유량을 빼는 잔여 분포 접근법을 도입한다.
• 균형 지표 r = ΔUⁿᵢ / ΔUᵉᑫᵢ에 따라 변하는 유량 제한 함수 φ(r)를 도입하여 보정을 제어하고 단조성을 유지한다.
• 스키마는 ∂uᵢ/∂t = −(𝒰ᵢ₊₁/₂ − 𝒰ᵢ₋₁/₂)/Δx + (uᵢ − uᵢᵉᑫ)로 표현되며, 𝒰ᵢ±₁/₂ = Uᵘᵐᵢ±₁/₂ − φᵢUᵉᑫᵢ±₁/₂로 정의되어 평형 상태 유지 보장을 받는다.
• 선형 문제의 경우, uᵢ → uᵢ − uᵢᵉᑫ의 변환은 잔여 스킴을 기본 방법과 동일하게 만들며, 그 성질을 유지한다.
• 비선형 스킴의 경우, 0 ≤ φ(r)/r ≤ 1이 되도록 유량 제한 함수 φ(r)를 설계함으로써 적절한 시간 간격 조건 하에서 TVD 안정성을 확보한다.
• 대류-확산 방정식, 소스항이 있는 쌍곡 시스템, 운동 방정식에 대한 검증을 통해 수치 결과가 정확도와 잘 균형 잡힌 행동을 확인한다.

실험 결과

연구 질문

• RQ1 특정 문제에 맞는 유량 설계 없이도 시간에 의존하는 PDE에서 평형 상태 해를 유지할 수 있는 일반적 프레임워크를 개발할 수 있는가?
• RQ2 소스항이 존재하는 비선형 PDE에 대해 고차 정확도를 유지하면서 잘 균형 잡힌 스킴을 어떻게 확보할 수 있는가?
• RQ3 비선형 스킴에서 평형 상태 유지와 TVD(Total Variation Diminishing) 안정성을 동시에 확보할 수 있는 유량 제한 전략은 무엇인가?
• RQ4 미세-거시 분해 기법은 운동 방정식을 초월하여 다른 PDE 클래스로 얼마나 일반화될 수 있는가?

주요 결과

• 잔여 평형 스킴은 대류-확산 방정식에서 정확히 평형 상태 해 uᵉᑫ(𝑥)를 유지한다. 이를 반경분리 형태 ∂uᵢ/∂t = −(𝒰ᵢ₊₁/₂ − 𝒰ᵢ₋₁/₂)/Δx + (uᵢ − uᵢᵉᑫ)로 확인할 수 있다.
• uᵉᑫ(𝑥) = u_B exp(−𝑥/a)인 선형 대류-확산 문제에서, 메쉬 해상도나 시간 간격에 관계없이 유량 보정이 적용된 한 평형 상태가 유지된다.
• r ≤ 1일 땐 r^α, r > 1일 땐 r^−α로 정의된 유량 제한 함수 φ(r)는 0 ≤ φ(r)/r ≤ 1 조건 하에서 잘 균형 잡힌 행동과 TVD 안정성을 모두 보장한다 (α > 1).
• 수치 결과는 강한 소스항과 비연속적 초깃값이 존재하는 상황에서도 스킴이 고차 정확도를 유지하고 임의의 진동을 피하는 것으로 확인된다.
• TVD 성질은 불등식 TV(uⁿ⁺¹) ≤ (1 − Δt)TV(uⁿ) + ΔtTV(uᵉᑫ)를 통해 보장되며, 이는 TV(uᵉᑫ) ≤ TV(uⁿ)일 경우에 성립한다.
• 이 방법은 소스항이 있는 보존 법칙 시스템과 운동 방정식으로도 확장 가능하며, 고려된 원형 문제를 초월해 광범위한 적용 가능성을 보여준다.