Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA Fourier Pseudospectral Method for the Good Boussinesq Equation with Second Order Temporal Accuracy
102
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
In this paper, we discuss the nonlinear stability and convergence of a fully discrete Fourier pseudospectral method coupled with a specially designed second order time-stepping for the numerical solution of the good Boussinesq equation. Our analysis improves the existing results presented in earlier literature in two ways. First, an $l_infty(0, T^*; H2)$ convergence for the solution and $l_infty(0, T^*; l_2)$ convergence for the time-derivative of the solution are obtained in this paper, instead of the $l_infty(0, T^*; l_2)$ convergence for the solution and the $l_infty(0, T^*; H^{-2})$ convergence for the time-derivative, given in [17]. In addition, the stability and convergence of this method is shown to be unconditional for the time step in terms of the spatial grid size, compared with a severe restriction time step restriction $Delta t leq Ch^2$ reported in [17].
rate research
Read More
Recent theoretical and experimental advances show that the inertia of magnetization emerges at sub-picoseconds and contributes to the ultrafast magnetization dynamics which cannot be captured intrinsically by the LLG equation. Therefore, as a generalization, the inertial Landau-Lifshitz-Gilbert (iLLG) equation is proposed to model the ultrafast magnetization dynamics. Mathematically, the LLG equation is a nonlinear system of parabolic type with (possible) degeneracy. However, the iLLG equation is a nonlinear system of mixed hyperbolic-parabolic type with degeneracy, and exhibits more complicated structures. It behaves like a hyperbolic system at the sub-picosecond scale while behaves like a parabolic system at larger timescales. Such hybrid behaviors impose additional difficulties on designing numerical methods for the iLLG equation. In this work, we propose a second-order semi-implicit scheme to solve the iLLG equation. The second temporal derivative of magnetization is approximated by the standard centered difference scheme and the first derivative is approximated by the midpoint scheme involving three time steps. The nonlinear terms are treated semi-implicitly using one-sided interpolation with the second-order accuracy. At each step, the unconditionally unique solvability of the unsymmetric linear system of equations in the proposed method is proved with a detailed discussion on the condition number. Numerically, the second-order accuracy in both time and space is verified. Using the proposed method, the inertial effect of ferromagnetics is observed in micromagnetics simulations at small timescales, in consistency with the hyperbolic property of the model at sub-picoseconds. For long time simulations, the results of the iLLG model are in nice agreements with those of the LLG model, in consistency with the parabolic feature of the iLLG model at larger timescales.
We develop an inverse scattering transform formalism for the good Boussinesq equation on the line. Assuming that the solution exists, we show that it can be expressed in terms of the solution of a $3 times 3$ matrix Riemann-Hilbert problem. The Riemann-Hilbert problem is formulated in terms of two reflection coefficients whose definitions involve only the initial data, and it has a form which makes it suitable for the evaluation of long-time asymptotics via Deift-Zhou steepest descent arguments.
In this paper, we study a novel second-order energy stable Backward Differentiation Formula (BDF) finite difference scheme for the epitaxial thin film equation with slope selection (SS). One major challenge for the higher oder in time temporal discretization is how to ensure an unconditional energy stability and an efficient numerical implementation. We propose a general framework for designing the higher order in time numerical scheme with unconditional energy stability by using the BDF method with constant coefficient stabilized terms. Based on the unconditional energy stability property, we derive an $L^infty_h (0,T; H_{h}^2)$ stability for the numerical solution and provide an optimal the convergence analysis. To deal with the 4-Laplacian solver in an $L^{2}$ gradient flow at each time step, we apply an efficient preconditioned steepest descent algorithm and preconditioned nonlinear conjugate gradient algorithm to solve the corresponding nonlinear system. Various numerical simulations are present to demonstrate the stability and efficiency of the proposed schemes and slovers.
We consider the initial-value problem for the ``good Boussinesq equation on the line. Using inverse scattering techniques, the solution can be expressed in terms of the solution of a $3 times 3$-matrix Riemann-Hilbert problem. We establish formulas for the long-time asymptotics of the solution by performing a Deift-Zhou steepest descent analysis of a regularized version of this Riemann-Hilbert problem.
In this paper, we analyse a new exponential-type integrator for the nonlinear cubic Schrodinger equation on the $d$ dimensional torus $mathbb T^d$. The scheme has recently also been derived in a wider context of decorated trees in [Y. Bruned and K. Schratz, arXiv:2005.01649]. It is explicit and efficient to implement. Here, we present an alternative derivation, and we give a rigorous error analysis. In particular, we prove second-order convergence in $H^gamma(mathbb T^d)$ for initial data in $H^{gamma+2}(mathbb T^d)$ for any $gamma > d/2$. This improves the previous work in [Knoller, A. Ostermann, and K. Schratz, SIAM J. Numer. Anal. 57 (2019), 1967-1986]. The design of the scheme is based on a new method to approximate the nonlinear frequency interaction. This allows us to deal with the complex resonance structure in arbitrary dimensions. Numerical experiments that are in line with the theoretical result complement this work.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
