ترغب بنشر مسار تعليمي؟ اضغط هنا

A new class of high-order energy-preserving schemes for the Korteweg-de Vries equation based on the quadratic auxiliary variable (QAV) approach

135   0   0.0 ( 0 )
 نشر من قبل Yuezheng Gong
 تاريخ النشر 2021
  مجال البحث الهندسة المعلوماتية
والبحث باللغة English




اسأل ChatGPT حول البحث

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.


قيم البحث

اقرأ أيضاً

88 - Kai Yang 2021
This paper proposes a new class of arbitarily high-order conservative numerical schemes for the generalized Korteweg-de Vries (KdV) equation. This approach is based on the scalar auxiliary variable (SAV) method. The equation is reformulated into an e quivalent system by introducing a scalar auxiliary variable, and the energy is reformulated into a sum of two quadratic terms. Therefore, the quadratic preserving Runge-Kutta method will preserve both the mass and the reformulated energy in the discrete time flow. With the Fourier pseudo-spectral spatial discretization, the scheme conserves the first and third invariant quantities (momentum and energy) exactly in the fully discrete sense. The discrete mass possesses the precision of the spectral accuracy.
In this paper, we design a novel linearized and momentum-preserving Fourier pseudo-spectral scheme to solve the Rosenau-Korteweg de Vries equation. With the aid of a new semi-norm equivalence between the Fourier pseudo-spectral method and the finite difference method, a prior bound of the numerical solution in discrete $L^{infty}$-norm is obtained from the discrete momentum conservation law. Subsequently, based on the energy method and the bound of the numerical solution, we show that, without any restriction on the mesh ratio, the scheme is convergent with order $O(N^{-s}+tau^2)$ in discrete $L^infty$-norm, where $N$ is the number of collocation points used in the spectral method and $tau$ is the time step. Numerical results are addressed to confirm our theoretical analysis.
In this paper, we propose a novel family of high-order numerical schemes for the gradient flow models based on the scalar auxiliary variable (SAV) approach, which is named the high-order scalar auxiliary variable (HSAV) method. The newly proposed sch emes could be shown to reach arbitrarily high order in time while preserving the energy dissipation law without any restriction on the time step size (i.e., unconditionally energy stable). The HSAV strategy is rather general that it does not depend on the specific expression of the effective free energy, such that it applies to a class of thermodynamically consistent gradient flow models arriving at semi-discrete high-order energy-stable schemes. We then employ the Fourier pseudospectral method for spatial discretization. The fully discrete schemes are also shown to be unconditionally energy stable. Furthermore, we present several numerical experiments on several widely-used gradient flow models, to demonstrate the accuracy, efficiency and unconditionally energy stability of the HSAV schemes. The numerical results verify that the HSAV schemes can reach the expected order of accuracy, and it allows a much larger time step size to reach the same accuracy than the standard SAV schemes.
In this paper, a class of arbitrarily high-order linear momentum-preserving and energy-preserving schemes are proposed, respectively, for solving the regularized long-wave equation. For the momentum-preserving scheme, our key ideas mainly follow the extrapolation/prediction-correction technique and symplectic Runge-Kutta (RK) methods in time combined with the standard Fourier pseudo-spectral method in space. We show that it is uniquely solvable, unconditionally stable and can exactly preserve the momentum of the system. Subsequently, based on the energy quadratization approach and the analogous linearized idea used in the construction of the linear momentum-preserving scheme, the energy-preserving scheme is presented and it is proven to preserve both the discrete mass and quadratic energy. Numerical results are addressed to demonstrate the accuracy and efficiency of the schemes.
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 gradi ent 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.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا