Emergent Mind

High-Order Schemes of Exponential Time Differencing for Stiff Systems with Nondiagonal Linear Part

(2208.14292)
Published Aug 30, 2022 in math.NA , cs.NA , and nlin.PS

Abstract

Exponential time differencing methods is a power toll for high-performance numerical simulation of computationally challenging problems in condensed matter and chemical physics, where mathematical models often possess fast oscillating or decaying modes -- in other words, are stiff systems. Practical implementation of these methods for the systems with nondiagonal linear part of equations is exacerbated by infeasibility of an analytical calculation of the exponential of a nondiagonal linear operator; in this case, the coefficients of the exponential time differencing scheme cannot be calculated analytically. We suggest an approach, where these coefficients are numerically calculated with auxiliary problems. We rewrite the high-order Runge--Kutta type schemes in terms of the solutions to these auxiliary problems and practically examine the accuracy and computational performance of these methods for a heterogeneous Cahn--Hilliard equation and a sixth-order spatial derivative equation governing pattern formation in the presence of an additional conservation law.

We're not able to analyze this paper right now due to high demand.

Please check back later (sorry!).

Generate a summary of this paper on our Pro plan:

We ran into a problem analyzing this paper.

Newsletter

Get summaries of trending comp sci papers delivered straight to your inbox:

Unsubscribe anytime.