Related papers: Higher order operator splitting Fourier spectral m…
An efficient geometric integrator is proposed for solving the perturbed Kepler motion. This method is stable and accurate over long integration time, which makes it appropriate for treating problems in astrophysics, like solar system…
In this paper, we propose and analyze a time-stepping method for the time fractional Allen-Cahn equation. The key property of the proposed method is its unconditional stability for general meshes, including the graded mesh commonly used for…
Operator splitting schemes have been successfully used in computational sciences to reduce complex problems into a series of simpler subproblems. Since 1950s, these schemes have been widely used to solve problems in PDE and control.…
The Cahn-Hilliard-Navier-Stokes (CHNS) equations represent the fundamental building blocks of hydrodynamic phase-field models for multiphase fluid flow dynamics. Due to the coupling between the Navier-Stokes equation and the Cahn-Hilliard…
We consider a class of second-order Strang splitting methods for Allen-Cahn equations with polynomial or logarithmic nonlinearities. For the polynomial case both the linear and the nonlinear propagators are computed explicitly. We show that…
A set of semi-analytical techniques based on Fourier analysis is used to solve wave scattering problems in variously shaped waveguides with varying normal admittance boundary conditions. Key components are newly developed conformal mapping…
The precise description of quantum nuclear fluctuations in atomistic modelling is possible by employing path integral techniques, which involve a considerable computational overhead due to the need of simulating multiple replicas of the…
A practical and simple stable method for calculating Fourier integrals is proposed, effective both at low and at high frequencies. An approach based on the fruitful idea of Levin, to use of the collocation method to approximate the slowly…
Computing solutions to partial differential equations using the fast Fourier transform can lead to unwanted oscillatory behavior. Due to the periodic nature of the discrete Fourier transform, waves that leave the computational domain on one…
Prior to the recent development of symplectic integrators, the time-stepping operator $\e^{h(A+B)}$ was routinely decomposed into a sum of products of $\e^{h A}$ and $\e^{hB}$ in the study of hyperbolic partial differential equations. In…
Given an unconditionally stable algorithm for solving the Cahn-Hilliard equation, we present a general calculation for an analytic time step $\d \tau$ in terms of an algorithmic time step $\dt$. By studying the accumulative multi-step error…
The explicit split-operator algorithm is often used for solving the linear and nonlinear time-dependent Schr\"{o}dinger equations. However, when applied to certain nonlinear time-dependent Schr\"{o}dinger equations, this algorithm loses…
We consider the Cauchy problem for a second-order evolution equation, in which the problem operator is the sum of two self-adjoint operators. The main feature of the problem is that one of the operators is represented in the form of the…
We consider Hamiltonian PDEs that can be split into a linear unbounded operator and a regular non linear part. We consider abstract splitting methods associated with this decomposition where no discretization in space is made. We prove a…
We study operator-splitting schemes for approximating Koopman generators of linear semigroups induced by nonlinear flows, a framework originating with Dorroh and Neuberger. Building on ideas of Lie, Kowalewski, and Gr\"{o}bner, we analyze…
We propose a new fourth-order compact time-splitting ($S_\text{4c}$) Fourier pseudospectral method for the Dirac equation by splitting the Dirac equation into two parts together with using the double commutator between them to integrate the…
The object of the present work is to present the new classes of third-order and fourth-order iterative methods for solving nonlinear equations. Our third-order method includes methods of Weerakoon \cite{Weerakoon}, Homeier \cite{Homeier2},…
In this article we provide a splitting method for solving monotone inclusions in a real Hilbert space involving four operators: a maximally monotone, a monotone-Lipschitzian, a cocoercive, and a monotone-continuous operator. The proposed…
We develop a fourth order simulation algorithm for solving the stochastic Langevin equation. The method consists of identifying solvable operators in the Fokker-Planck equation, factorizing the evolution operator for small time steps to…
In this work, we present a second-order nonuniform time-stepping scheme for the time-fractional Allen-Cahn equation. We show that the proposed scheme preserves the discrete maximum principle, and by using the convolution structure of…