Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA linearly implicit structure-preserving scheme for the Camassa-Holm equation based on multiple scalar auxiliary variables approach
62
0
0.0
(
0
)
Added by
Chaolong Jiang
Publication date
2019
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
In this paper, we present a linearly implicit energy-preserving scheme for the Camassa-Holm equation by using the multiple scalar auxiliary variables approach, which is first developed to construct efficient and robust energy stable schemes for gradient systems. The Camassa-Holm equation is first reformulated into an equivalent system by utilizing the multiple scalar auxiliary variables approach, which inherits a modified energy. Then, the system is discretized in space aided by the standard Fourier pseudo-spectral method and a semi-discrete system is obtained, which is proven to preserve a semi-discrete modified energy. Subsequently, the linearized Crank-Nicolson method is applied for the resulting semi-discrete system to arrive at a fully discrete scheme. The main feature of the new scheme is to form a linear system with a constant coefficient matrix at each time step and produce numerical solutions along which the modified energy is precisely conserved, as is the case with the analytical solution. Several numerical results are addressed to confirm accuracy and efficiency of the proposed scheme.
rate research
Read More
A novel class of high-order linearly implicit energy-preserving exponential integrators are proposed for the nonlinear Schrodinger equation. We firstly done that the original equation is reformulated into a new form with a modified quadratic energy by the scalar auxiliary variable approach. The spatial derivatives of the system are then approximated with the standard Fourier pseudo-spectral method. Subsequently, we apply the extrapolation technique to the nonlinear term of the semi-discretized system and a linearized system is obtained. Based on the Lawson transformation, the linearized system is rewritten as an equivalent one and we further apply the symplectic Runge-Kutta method to the resulting system to gain a fully discrete scheme. We show that the proposed scheme can produce numerical solutions along which the modified energy is precisely conserved, as is the case with the analytical solution and is extremely efficient in the sense that only linear equations with constant coefficients need to be solved at every time step. Numerical results are addressed to demonstrate the remarkable superiority of the proposed schemes in comparison with other high-order structure-preserving method.
We present two semidiscretizations of the Camassa-Holm equation in periodic domains based on variational formulations and energy conservation. The first is a periodic version of an existing conservative multipeakon method on the real line, for which we propose efficient computation algorithms inspired by works of Camassa and collaborators. The second method, and of primary interest, is the periodic counterpart of a novel discretization of a two-component Camassa-Holm system based on variational principles in Lagrangian variables. Applying explicit ODE solvers to integrate in time, we compare the variational discretizations to existing methods over several numerical examples.
In this paper, we develop a new class of high-order energy-preserving schemes for the Korteweg-de Vries equation based on the quadratic auxiliary variable technique, which can conserve the original energy of the system. By introducing a quadratic auxiliary variable, the original system is reformulated into an equivalent form with a modified quadratic energy, where the way of the introduced variable naturally produces a quadratic invariant of the new system. A class of Runge-Kutta methods satisfying the symplectic condition is applied to discretize the reformulated model in time, arriving at arbitrarily high-order schemes, which naturally conserve the modified quadratic energy and the produced quadratic invariant. Under the consistent initial condition, the proposed methods are rigorously proved to maintain the original energy conservation law of the Korteweg-de Vries equation. In order to match the high order precision of time, the Fourier pseudo-spectral method is employed for spatial discretization, resulting in fully discrete energy-preserving schemes. To solve the proposed nonlinear schemes effectively, we present a very efficient practically-structure-preserving iterative technique, which not only greatly saves the calculation cost, but also achieves the purpose of practically preserving structure. Ample numerical results are addressed to confirm the expected order of accuracy, conservative property and efficiency of the proposed schemes. This new class of numerical strategies is rather general so that they can be readily generalized for any conservative systems with a polynomial energy.
We put forward and analyze an explicit finite difference scheme for the Camassa-Holm shallow water equation that can handle general $H^1$ initial data and thus peakon-antipeakon interactions. Assuming a specified condition restricting the time step in terms of the spatial discretization parameter, we prove that the difference scheme converges strongly in $H^1$ towards a dissipative weak solution of Camassa-Holm equation.
Systems driven by multiple physical processes are central to many areas of science and engineering. Time discretization of multiphysics systems is challenging, since different processes have different levels of stiffness and characteristic time scales. The multimethod approach discretizes each physical process with an appropriate numerical method; the methods are coupled appropriately such that the overall solution has the desired accuracy and stability properties. The authors developed the general-structure additive Runge-Kutta (GARK) framework, which constructs multimethods based on Runge-Kutta schemes. This paper constructs the new GARK-ROS/GARK-ROW families of multimethods based on linearly implicit Rosenbrock/Rosenbrock-W schemes. For ordinary differential equation models, we develop a general order condition theory for linearly implicit methods with any number of partitions, using exact or approximate Jacobians. We generalize the order condition theory to two-way partitioned index-1 differential-algebraic equations. Applications of the framework include decoupled linearly implicit, linearly implicit/explicit, and linearly implicit/implicit methods. Practical GARK-ROS and GARK-ROW schemes of order up to four are constructed.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
