Related papers: Efficient boundary corrected Strang splitting
Exponential Runge-Kutta methods are a well-established tool for the numerical integration of parabolic evolution equations. However, these schemes are typically developed under the assumption of homogeneous boundary conditions. In this…
In this paper we consider splitting methods for nonlinear ordinary differential equations in which one of the (partial) flows that results from the splitting procedure can not be computed exactly. Instead, we insert a well-chosen state…
The implementation of multi-stage splitting integrators is essentially the same as the implementation of the familiar Strang/Verlet method. Therefore multi-stage formulas may be easily incorporated into software that now uses the…
We extend to the N-level Bloch model the splitting scheme which use exact numerical solutions of sub-equations. These exact solutions involve matrix exponentials which we want to avoid to calculate at each time step. We use Newton…
Fractional-order elliptic problems are investigated in case of inhomogeneous Dirichlet boundary data. The boundary integral form is proposed as a suitable mathematical model. The corresponding theory is completed by sharpening the mapping…
We propose and analyze a second-order Strang splitting method for a class of stiff matrix differential equations with Sylvester-type structure. The method splits the dynamics into a stiff linear part, treated exactly via matrix…
We provide general product formulas for the solutions of non-autonomous abstract Cauchy problems. The main technical tool is the application of evolution semigroup methods, allowing the direct application of existing results on autonomous…
A rigorous convergence analysis of the Strang splitting algorithm for Vlasov-type equations in the setting of abstract evolution equations is provided. It is shown that under suitable assumptions the convergence is of second order in the…
In this paper, we present a class of high-order and efficient compact difference schemes for nonlinear convection diffusion equations, which can preserve both bounds and mass. For the one-dimensional problem, we first introduce a high-order…
We present a parametric family of semi-implicit second order accurate numerical methods for non-conservative and conservative advection equation for which the numerical solutions can be obtained in a fixed number of forward and backward…
A practical and efficient scheme for the higher order integration of the Landau-Lifschitz-Gilbert (LLG) equation is presented. The method is based on extrapolation of the two-step explicit midpoint rule and incorporates adaptive time step…
We derive a numerical method, based on operator splitting, to abstract parabolic semilinear boundary coupled systems. The method decouples the linear components which describe the coupling and the dynamics in the bulk and on the surface,…
We propose a boundary value correction approach for cases when curved boundaries are approximated by straight lines (planes) and Lagrange multipliers are used to enforce Dirichlet boundary conditions. The approach allows for optimal order…
The error behavior of exponential operator splitting methods for nonlinear Schr{\"o}dinger equations in the semiclassical regime is studied. For the Lie and Strang splitting methods, the exact form of the local error is determined and the…
For systems of the form $\dot q = M^{-1} p$, $\dot p = -Aq+f(q)$, common in many applications, we analyze splitting integrators based on the (linear/nonlinear) split systems $\dot q = M^{-1} p$, $\dot p = -Aq$ and $\dot q = 0$, $\dot p =…
An equation containing a fractional power of an elliptic operator of second order is studied for Dirichlet boundary conditions. Finite difference approximations in space are employed. The proposed numerical algorithm is based on solving an…
We study the Strang splitting scheme for quasilinear Schr\"odinger equations. We establish the convergence of the scheme for solutions with small initial data. We analyze the linear instability of the numerical scheme, which explains the…
Classical and new numerical schemes are generated using evolutionary computing. Differential Evolution is used to find the coefficients of finite difference approximations of function derivatives, and of single and multi-step integration…
We present a number of new contributions to the topic of constructing efficient higher-order splitting methods for the numerical integration of evolution equations. Particular schemes are constructed via setup and solution of polynomial…
Splitting methods constitute a widely used class of numerical integrators for ordinary and partial differential equations, particularly well suited to problems that can be decomposed into simpler subproblems. High-order splitting schemes…