Computation of some dispersive equations through their iterated linearisation

Guannan Chen,Arieh Iserles,Karolina Kropielnicka,Pranav Singh
2024-05-08
Abstract:It is often the case that, while the numerical solution of the non-linear dispersive equation $\mathrm{i}\partial_t u(t)=\mathcal{H}(u(t),t)u(t)$ represents a formidable challenge, it is fairly easy and cheap to solve closely related linear equations of the form $\mathrm{i}\partial_t u(t)=\mathcal{H}_1(t)u(t)+\widetilde{\mathcal H}_2(t)u(t)$, where $\mathcal{H}_1(t)+\mathcal{H}_2(v,t)=\mathcal{H}(v,t)$. In that case we advocate an iterative linearisation procedure that involves fixed-point iteration of the latter equation to solve the former. A typical case is when the original problem is a nonlinear Schrödinger or Gross--Pitaevskii equation, while the `easy' equation is linear Schrödinger with time-dependent potential.
Numerical Analysis,Computational Physics
What problem does this paper attempt to address?
This paper mainly discusses a new method for solving nonlinear dispersion equations - by iteratively linearizing them. Nonlinear dispersion equations are usually difficult to solve, but the associated linear equations are relatively easy and cheap. The authors propose an iterative scheme to approximate the solution of the nonlinear equation by continuously linearizing the original equation. Specifically, they write the nonlinear equation in the form of a linear part and a nonlinear part, and then solve it step by step using the fixed point iteration method. In each iteration, they solve a linear dispersion equation with a time-variable potential to approximate the solution of the nonlinear equation. The paper analyzes the convergence, speed, and preservation of some geometric properties of the original equation, such as Hamiltonian property or L2 norm conservation, of the iterative scheme. They effectively approximate the linearized but non-autonomous equation by using a combination of Magnus expansion and Hermite integration. The paper also demonstrates the application of this method using the cubic nonlinear Schrödinger equation and the Gross-Pitaevskii equation, and compares it with the splitting method and the Krylov subspace method through numerical experiments, proving the effectiveness of the proposed method. The main contribution of the paper is the introduction of an iterative linearization strategy, which approximates the solution of the nonlinear dispersion equation by repeatedly solving linear equations, while analyzing the convergence and structural properties of this method. This method can be used in combination with other methods for approximating matrix exponentials, providing a new approach for solving nonlinear dispersion equations.