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

High order conservative schemes for the generalized Benjamin-Ono equation in the unbounded domain

111   0   0.0 ( 0 )
 نشر من قبل Kai Yang
 تاريخ النشر 2021
  مجال البحث الهندسة المعلوماتية
والبحث باللغة English
 تأليف Kai Yang




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

This paper proposes a new class of mass or energy conservative numerical schemes for the generalized Benjamin-Ono (BO) equation on the whole real line with arbitrarily high-order accuracy in time. The spatial discretization is achieved by the pseudo-spectral method with the rational basis functions, which can be implemented by the Fast Fourier transform (FFT) with the computational cost $mathcal{O}( Nlog(N))$. By reformulating the spatial discretized system into the different equivalent forms, either the spatial semi-discretized mass or energy can be preserved exactly under the continuous time flow. Combined with the symplectic Runge-Kutta, with or without the scalar auxiliary variable reformulation, the fully discrete energy or mass conservative scheme can be constructed with arbitrarily high-order temporal accuracy, respectively. Our numerical results show the conservation of the proposed schemes, and also the superior accuracy and stability to the non-conservative (Leap-frog) scheme.

قيم البحث

اقرأ أيضاً

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.
We consider the generalized Benjamin-Ono (gBO) equation on the real line, $ u_t + partial_x (-mathcal H u_{x} + tfrac1{m} u^m) = 0, x in mathbb R, m = 2,3,4,5$, and perform numerical study of its solutions. We first compute the ground state solution to $-Q -mathcal H Q^prime +frac1{m} Q^m = 0$ via Petviashvilis iteration method. We then investigate the behavior of solutions in the Benjamin-Ono ($m=2$) equation for initial data with different decay rates and show decoupling of the solution into a soliton and radiation, thus, providing confirmation to the soliton resolution conjecture in that equation. In the mBO equation ($m=3$), which is $L^2$-critical, we investigate solutions close to the ground state mass, and, in particular, we observe the formation of stable blow-up above it. Finally, we focus on the $L^2$-supercritical gBO equation with $m=4,5$. In that case we investigate the global vs finite time existence of solutions, and give numerical confirmation for the dichotomy conjecture, in particular, exhibiting blow-up phenomena in the supercritical setting.
We consider a higher-dimensional version of the Benjamin-Ono (HBO) equation in the 2D setting: $u_t- mathcal{R}_1 Delta u + frac{1}{2}(u^2)_x=0, (x,y) in mathbb{R}^2$, which is $L^2$-critical, and investigate properties of solutions both analytically and numerically. For a generalized equation (fractional 2D gKdV) after deriving the Pohozaev identities, we obtain non-existence conditions for solitary wave solutions, then prove uniform bounds in the energy space or conditional global existence, and investigate the radiation region, a specific wedge in the negative $x$-direction. We then introduce our numerical approach in a general context, and apply it to obtain the ground state solution in the 2D critical HBO equation, then show that its mass is a threshold for global vs. finite time existing solutions, which is typical in the focusing (mass-)critical dispersive equations. We also observe that globally existing solutions tend to disperse completely into the radiation in this nonlocal equation. The blow-up solutions travel in the positive $x$-direction with the rescaled ground state profile while also radiating dispersive oscillations into the radiative wedge. We conclude with examples of different interactions of two solitary wave solutions, including weak and strong interactions.
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 design a novel class of arbitrarily high-order structure-preserving numerical schemes for the time-dependent Gross-Pitaevskii equation with angular momentum rotation in three dimensions. Based on the idea of the scalar auxiliary var iable approach which is proposed in the recent papers [J. Comput. Phys., 416 (2018) 353-407 and SIAM Rev., 61(2019) 474-506] for developing energy stable schemes for gradient flow systems, we firstly reformulate the Gross-Pitaevskii equation into an equivalent system with a modified energy conservation law. The reformulated system is then discretized by the Gauss collocation method in time and the standard Fourier pseudo-spectral method in space, respectively. We show that the proposed schemes can preserve the discrete mass and modified energy exactly. Numerical results are addressed to verify the efficiency and high-order accuracy of the proposed schemes.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

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