An official website of the United States government
Abstract The inverse scattering transform allows explicit construction of solutions to many physically significant nonlinear wave equations. Notably, this method can be extended to fractional nonlinear evolution equations characterized by anomalous dispersion using completeness of suitable eigenfunctions of the associated linear scattering problem. In anomalous diffusion, the mean squared displacement is proportional to t α , α > 0, while in anomalous dispersion, the speed of localized waves is proportional to A α , where A is the amplitude of the wave. Fractional extensions of the modified Korteweg–deVries (mKdV), sine-Gordon (sineG) and sinh-Gordon (sinhG) and associated hierarchies are obtained. Using symmetries present in the linear scattering problem, these equations can be connected with a scalar family of nonlinear evolution equations of which fractional mKdV (fmKdV), fractional sineG (fsineG), and fractional sinhG (fsinhG) are special cases. Completeness of solutions to the scalar problem is obtained and, from this, the nonlinear evolution equation is characterized in terms of a spectral expansion. In particular, fmKdV, fsineG, and fsinhG are explicitly written. One-soliton solutions are derived for fmKdV and fsineG using the inverse scattering transform and these solitons are shown to exhibit anomalous dispersion.
more »
« less
Christov expansion method for nonlocal nonlinear evolution equations
Christov functions are a complete orthonormal set of functions on L^2(-∞,∞) that allow us to expand derivatives, nonlinear products, and nonlocal (integro-differential) terms back into the same basis. These properties are beneficial when solving nonlinear evolution equations using Galerkin spectral methods. In this work, we demonstrate such a “Christov expansion method” for the Benjamin–Ono (BO) equation. In the BO equation, the dispersion term is nonlocal, given by the Hilbert transform of the second spatial derivative of the unknown function. The Hilbert transform of the Christov functions can be computed using complex integration and Cauchy’s residue theorem to obtain simple relations. Then, a Galerkin spectral expansion can be used to the solve the BO equation. Time integration is performed using a Crank–Nicolson-type scheme. Importantly, the Christov expansion method yields a banded matrix for the spatial discretization, even though the spatial terms are nonlocal. To demonstrate the approach and its implementation, we perform numerical experiments showing the steady propagation of single and the overtaking interaction of multiple BO solitary waves.
more »
« less
- Award ID(s):
- 2029540
- PAR ID:
- 10481631
- Editor(s):
- Todorov, M D
- Publisher / Repository:
- IOP Publishing
- Date Published:
- Journal Name:
- Journal of Physics: Conference Series
- Volume:
- 2675
- ISSN:
- 1742-6588
- Page Range / eLocation ID:
- 012022
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
In this work, we investigate the two-component modified Korteweg-de Vries (mKdV) equation, which is a complete integrable system, and accepts a generalization of 4 × 4 matrix Ablowitz–Kaup–Newell-Segur (AKNS)-type Lax pair. By using of the unified transform approach, the initial-boundary value (IBV) problem of the two-component mKdV equation associated with a 4 × 4 matrix Lax pair on the half-line will be analyzed. Supposing that the solution {u1(x, t), u2(x, t)} of the two-component mKdV equation exists, we will show that it can be expressed in terms of the unique solution of a 4 × 4 matrix Riemann–Hilbert problem formulated in the complex λ-plane. Moreover, we will prove that some spectral functions s(λ) and S(λ) are not independent of each other but meet the global relationship.more » « less
-
In this paper, we propose an efficient exponential integrator finite element method for solving a class of semilinear parabolic equations in rectangular domains. The proposed method first performs the spatial discretization of the model equation using the nite element approximation with continuous multilinear rectangular basis functions, and then takes the explicit exponential Runge-Kutta approach for time integration of the resulting semi-discrete system to produce fully-discrete numerical solution. Under certain regularity assumptions, error estimates measured in H1-norm are successfully derived for the proposed schemes with one and two RK stages. More remarkably, the mass and coefficient matrices of the proposed method can be simultaneously diagonalized with an orthogonal matrix, which provides a fast solution process based on tensor product spectral decomposition and fast Fourier transform. Various numerical experiments in two and three dimensions are also carried out to validate the theoretical results and demonstrate the excellent performance of the proposed method.more » « less
-
Abstract Integrable standard and nonlocal derivative nonlinear Schrödinger equations are investigated. The direct and inverse scattering are constructed for these equations; included are both the Riemann–Hilbert and Gel’fand–Levitan–Marchenko approaches and soliton solutions. As a typical application, it is shown how these derivative NLS equations can be obtained as asymptotic limits from a nonlinear Klein–Gordon equation.more » « less
-
Summary This paper presents an approach for efficient uncertainty analysis (UA) using an intrusive generalized polynomial chaos (gPC) expansion. The key step of the gPC‐based uncertainty quantification(UQ) is the stochastic Galerkin (SG) projection, which can convert a stochastic model into a set of coupled deterministic models. The SG projection generally yields a high‐dimensional integration problem with respect to the number of random variables used to describe the parametric uncertainties in a model. However, when the number of uncertainties is large and when the governing equation of the system is highly nonlinear, the SG approach‐based gPC can be challenging to derive explicit expressions for the gPC coefficients because of the low convergence in the SG projection. To tackle this challenge, we propose to use a bivariate dimension reduction method (BiDRM) in this work to approximate a high‐dimensional integral in SG projection with a few one‐ and two‐dimensional integrations. The efficiency of the proposed method is demonstrated with three different examples, including chemical reactions and cell signaling. As compared to other UA methods, such as the Monte Carlo simulations and nonintrusive stochastic collocation (SC), the proposed method shows its superior performance in terms of computational efficiency and UA accuracy.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1088/1742-6596/2675/1/012022
