skip to main content


Title: Orbital stability of smooth solitary waves for the Degasperis-Procesi equation
The Degasperis-Procesi (DP) equation is an integrable Camassa-Holm-type model which is an asymptotic approximation for the unidirectional propagation of shallow water waves. This work establishes the orbital stability of localized smooth solitary waves to the DP equation on the real line, extending our previous work on their spectral stability [J. Math. Pures Appl. (9) 142 (2020), pp. 298–314]. The main difficulty stems from the fact that the natural energy space is a subspace of L 3 L^3 , but the translation symmetry for the DP equation gives rise to a conserved quantity equivalent to the L 2 L^2 -norm, resulting in L 3 L^3 higher-order nonlinear terms in the augmented Hamiltonian. But the usual coercivity estimate is in terms of L 2 L^2 norm for DP equation, which cannot be used to control the L 3 L^3 higher order term directly. The remedy is to observe that, given a sufficiently smooth initial condition satisfying some mild constraint, the L ∞ L^\infty orbital norm of the perturbation is bounded above by a function of its L 2 L^2 orbital norm, yielding the higher order control and the orbital stability in the L 2 ∩ L ∞ L^2\cap L^\infty space.  more » « less
Award ID(s):
1815079
NSF-PAR ID:
10382651
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Proceedings of the American Mathematical Society
Volume:
151
Issue:
763
ISSN:
0002-9939
Page Range / eLocation ID:
151 to 160
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We consider the linear third order (in time) PDE known as the SMGTJ-equation, defined on a bounded domain, under the action of either Dirichlet or Neumann boundary control \begin{document}$ g $\end{document}. Optimal interior and boundary regularity results were given in [1], after [41], when \begin{document}$ g \in L^2(0, T;L^2(\Gamma)) \equiv L^2(\Sigma) $\end{document}, which, moreover, in the canonical case \begin{document}$ \gamma = 0 $\end{document}, were expressed by the well-known explicit representation formulae of the wave equation in terms of cosine/sine operators [19], [17], [24,Vol Ⅱ]. The interior or boundary regularity theory is however the same, whether \begin{document}$ \gamma = 0 $\end{document} or \begin{document}$ 0 \neq \gamma \in L^{\infty}(\Omega) $\end{document}, since \begin{document}$ \gamma \neq 0 $\end{document} is responsible only for lower order terms. Here we exploit such cosine operator based-explicit representation formulae to provide optimal interior and boundary regularity results with \begin{document}$ g $\end{document} "smoother" than \begin{document}$ L^2(\Sigma) $\end{document}, qualitatively by one unit, two units, etc. in the Dirichlet boundary case. To this end, we invoke the corresponding results for wave equations, as in [17]. Similarly for the Neumann boundary case, by invoking the corresponding results for the wave equation as in [22], [23], [37] for control smoother than \begin{document}$ L^2(0, T;L^2(\Gamma)) $\end{document}, and [44] for control less regular in space than \begin{document}$ L^2(\Gamma) $\end{document}. In addition, we provide optimal interior and boundary regularity results when the SMGTJ equation is subject to interior point control, by invoking the corresponding wave equations results [42], [24,Section 9.8.2].

     
    more » « less
  2. Abstract

    This work analyzes the forward and inverse scattering series for scalar waves based on the Helmholtz equation and the diffuse waves from the time-independent diffusion equation, which are important partial differential equations (PDEs) in various applications. Different from previous works, which study the radius of convergence for the forward and inverse scattering series, the stability, and the approximation error of the series under theLpnorms, we study these quantities under the SobolevHsnorm, which associates with a general class ofL2-based function spaces. TheHsnorm has a natural spectral bias based on its definition in the Fourier domain: the cases < 0 biases towards the lower frequencies, while the cases > 0 biases towards the higher frequencies. We compare the stability estimates using differentHsnorms for both the parameter and data domains and provide a theoretical justification for the frequency weighting techniques in practical inversion procedures. We also provide numerical inversion examples to demonstrate the differences in the inverse scattering radius of convergence under different metric spaces.

     
    more » « less
  3. The existence and stability of the Landau equation (1936) in a general bounded domain with a physical boundary condition is a long-outstanding open problem. This work proves the global stability of the Landau equation with the Coulombic potential in a general smooth bounded domain with the specular reflection boundary condition for initial perturbations of the Maxwellian equilibrium states. The highlight of this work also comes from the low-regularity assumptions made for the initial distribution. This work generalizes the recent global stability result for the Landau equation in a periodic box (Kim et al. in Peking Math J, 2020). Our methods consist of the generalization of the wellposedness theory for the Fokker–Planck equation (Hwang et al. SIAM J Math Anal 50(2):2194–2232, 2018; Hwang et al. Arch Ration Mech Anal 214(1):183–233, 2014) and the extension of the boundary value problem to a whole space problem, as well as the use of a recent extension of De Giorgi–Nash–Moser theory for the kinetic Fokker–Planck equations (Golse et al. Ann Sc Norm Super Pisa Cl Sci 19(1):253–295, 2019) and the Morrey estimates (Bramanti et al. J Math Anal Appl 200(2):332–354, 1996) to further control the velocity derivatives, which ensures the uniqueness. Our methods provide a new understanding of the grazing collisions in the Landau theory for an initial-boundary value problem. 
    more » « less
  4. Abstract The distribution of natural frequencies of the Euler–Bernoulli beam subject to fully non-dissipative boundary conditions is investigated. The beam is clamped at the left end and equipped with a 4-parameter ($\alpha ,\beta ,k_1,k_2$) linear boundary feedback law at the right end. The $2 \times 2$ boundary feedback matrix relates the control input (a vector of velocity and its spatial derivative at the right end), to the output (a vector of shear and moment at the right end). The initial boundary value problem describing the dynamics of the beam has been reduced to the first order in time evolution equation in the state Hilbert space equipped with the energy norm. The dynamics generator has a purely discrete spectrum (the vibrational modes) denoted by $\{\nu _n\}_{n\in \mathbb {Z}^{\prime}}$. The role of the control parameters is examined and the following results have been proven: (i) when $\beta \neq 0$, the set of vibrational modes is asymptotically close to the vertical line on the complex $\nu$-plane given by the equation $\Re \nu = \alpha + (1-k_1k_2)/\beta$; (ii) when $\beta = 0$ and the parameter $K = (1-k_1 k_2)/(k_1+k_2)$ is such that $\left |K\right |\neq 1$ then the following relations are valid: $\Re (\nu _n/n) = O\left (1\right )$ and $\Im (\nu _n/n^2) = O\left (1\right )$ as $\left |n\right |\to \infty$; (iii) when $\beta =0$, $|K| = 1$, and $\alpha = 0$, then the following relations are valid: $\Re (\nu _n/n^2) = O\left (1\right )$ and $\Im (\nu _n/n) = O\left (1\right )$ as $\left |n\right |\to \infty$; (iv) when $\beta =0$, $|K| = 1$, and $\alpha>0$, then the following relations are valid: $\Re (\nu _n/\ln \left |n\right |) = O\left (1\right )$ and $\Im (\nu _n/n^2) = O\left (1\right )$ as $\left |n\right |\to \infty$. 
    more » « less
  5. null (Ed.)
    In this paper we propose and analyze a finite difference numerical scheme for the Poisson-Nernst-Planck equation (PNP) system. To understand the energy structure of the PNP model, we make use of the Energetic Variational Approach (EnVarA), so that the PNP system could be reformulated as a non-constant mobility H − 1 H^{-1} gradient flow, with singular logarithmic energy potentials involved. To ensure the unique solvability and energy stability, the mobility function is explicitly treated, while both the logarithmic and the electric potential diffusion terms are treated implicitly, due to the convex nature of these two energy functional parts. The positivity-preserving property for both concentrations, n n and p p , is established at a theoretical level. This is based on the subtle fact that the singular nature of the logarithmic term around the value of 0 0 prevents the numerical solution reaching the singular value, so that the numerical scheme is always well-defined. In addition, an optimal rate convergence analysis is provided in this work, in which many highly non-standard estimates have to be involved, due to the nonlinear parabolic coefficients. The higher order asymptotic expansion (up to third order temporal accuracy and fourth order spatial accuracy), the rough error estimate (to establish the ℓ ∞ \ell ^\infty bound for n n and p p ), and the refined error estimate have to be carried out to accomplish such a convergence result. In our knowledge, this work will be the first to combine the following three theoretical properties for a numerical scheme for the PNP system: (i) unique solvability and positivity, (ii) energy stability, and (iii) optimal rate convergence. A few numerical results are also presented in this article, which demonstrates the robustness of the proposed numerical scheme. 
    more » « less