[논문 리뷰] Kernel Based High Order "Explicit" Unconditionally-Stable Scheme for Nonlinear Degenerate Advection-Diffusion Equations

Andrew Christlieb, Wei Guo외 1명|arXiv (Cornell University)|2017. 07. 28.

Differential Equations and Numerical Methods참고 문헌 25인용 수 8

한 줄 요약

이 논문은 선형이 아닌 탈형상(advection-diffusion) 방정식을 위한 고차수, 조건부 안정성 없는 수치적 방법을 제안한다. 이 방법은 선형 변환 기반의 방법론(transpose)에 기반한 커널 기반 수식을 사용하며, 연속적인 콘볼루션, 고차수 WENO 재구성, 그리고 명시적 강안정성 보존(Strong-Stability-Preserving) 룬게-쿠타 시간 적분을 조합함으로써 최대 제3차 정확도를 달성하고, 표준 명시적 방법보다 훨씬 큰 시간 간격을 허용한다. 이는 안정성과 불필요한 진동을 피하면서도 계산 비용을 크게 절감할 수 있도록 한다.

ABSTRACT

In this paper, we present a novel numerical scheme for solving a class of nonlinear degenerate parabolic equations with non-smooth solutions. The proposed method relies on a special kernel based formulation of the solutions found in our early work on the method of lines transpose and successive convolution. In such a framework, a high order weighted essentially non-oscillatory (WENO) methodology and a nonlinear filter are further employed to avoid spurious oscillations. High order accuracy in time is realized by using the high order explicit strong-stability-preserving (SSP) Runge-Kutta method. Moreover, theoretical investigations of the kernel based formulation combined with an explicit SSP method indicates that the combined scheme is unconditionally stable and up to third order accuracy. Evaluation of the kernel based approach is done with a fast $\mathcal{O}(N)$ summation algorithm. The new method allows for much larger time step evolution compared with other explicit schemes with the same order accuracy, leading to remarkable computational savings.

연구 동기 및 목표

• 비연속성이나 급격한 전면을 가진 해가 존재하는 비선형 탈형상 포물형 방정식을 위한 고차수 수치적 방법을 개발한다.
• 특히 확산이 지배적인 문제에서 Δt ∝ Δx²와 같이 엄격한 시간 간격 제약을 초래하는 표준 명시적 방법의 문제를 해결한다.
• 공간 및 시간에서 고차수 정확도를 유지하면서도 조건부 안정성 없는 시간 적분을 달성한다.
• 고차수 WENO 재구성과 비선형 필터를 사용하여 비연속성 근처의 진동을 억제한다.
• 커널 기반 수식의 빠른 O(N) 합산 알고리즘을 통해 효율적인 계산을 가능하게 한다.

제안 방법

• 방법론은 선형 변환 기반의 방법론(transpose)에서 유도된 커널 기반 적분 수식을 사용하여, 각 시간 단계에서 편미분방정식(PDE)을 경계값 문제로 변환한다.
• 연속적인 콘볼루션(해결책 전개)을 사용하여 해를 고차수 공간 도함수를 포함하는 급수 형태로 표현함으로써 고차수 정확도를 달성한다.
• 비연속성 근처에서 뚜렷한 해상도를 유지하기 위해 고차수 가중 평균 비진동(WENO) 재구성을 공간 도함수에 적용한다.
• 급격한 기울기나 비연속성이 존재하는 영역에서 불필요한 진동을 억제하기 위해 비선형 필터를 적용한다.
• 안정성과 정확도를 유지하기 위해 고차수 강안정성 보존(SSP) 룬게-쿠타 방법을 사용하여 시간 적분을 수행한다.
• 커널 기반 적분을 평가하기 위해 빠른 O(N) 합산 알고리즘을 사용하여 계산 비용을 O(N²)에서 O(N)으로 감소시킨다.

실험 결과

연구 질문

• RQ1 MOL T 프레임워크 내에서 커널 기반 수식이 비선형 탈형상 대류-확산 방정식에 대해 조건부 안정성 없는 시간 적분을 달성할 수 있는가?
• RQ2 명시적 시간 스텝 방법을 사용함에도 불구하고 공간 및 시간에서 최대 제3차 정확도를 달성할 수 있는가, 안정성 제약 없이?
• RQ3 WENO 재구성과 비선형 필터가 비연속성 근처의 진동을 효과적으로 억제하면서도 고차수 정확도를 유지할 수 있는가?
• RQ4 제안된 방법이 표준 명시적 방법보다 훨씬 큰 시간 간격을 허용할 수 있는가, 유사한 정확도 단계를 유지하면서?
• RQ5 빠른 O(N) 합산 알고리즘이 커널 기반 수식의 정확도와 안정성을 유지하면서도 높은 효율성을 유지할 수 있는가?

주요 결과

• 제안된 방법은 이론적으로 조건부 안정성 없이 고차수 정확도(공간 및 시간에서 최대 제3차)를 보장한다.
• 표준 명시적 방법보다 훨씬 큰 시간 간격을 허용하여 계산 비용 절감이 가능하다.
• 수치 실험 결과는 설계된 고차수 정확도와 비연속성 및 급격한 전면 근처에서 비진동 행동을 확인한다.
• 커널 기반 수식과 함께 빠른 O(N) 합산 알고리즘을 사용함으로써 계산 효율성이 향상되어 복잡도를 O(N²)에서 O(N)으로 감소시킨다.
• WENO와 비선형 필터의 사용은 급격한 기울기 또는 비연속성이 존재하는 해에서 불필요한 진동을 성공적으로 억제한다.
• 이 방법은 다공성 매체 유동, 플라즈마 이동, 복사 전달을 모델링하는 방정식 등 다양한 방정식에 적용 가능하며, 그 안정성과 강건성 덕분이다.