skip to main content


Title: Well-posedness and dynamics of solutions to the generalized KdV with low power nonlinearity
Abstract We consider two types of the generalized Korteweg–de Vries equation, where the nonlinearity is given with or without absolute values, and, in particular, including the low powers of nonlinearity, an example of which is the Schamel equation. We first prove the local well-posedness of both equations in a weighted subspace of H 1 that includes functions with polynomial decay, extending the result of Linares et al (2019 Commun. Contemp. Math. 21 1850056) to fractional weights. We then investigate solutions numerically, confirming the well-posedness and extending it to a wider class of functions that includes exponential decay. We include a comparison of solutions to both types of equations, in particular, we investigate soliton resolution for the positive and negative data with different decay rates. Finally, we study the interaction of various solitary waves in both models, showing the formation of solitons, dispersive radiation and even breathers, all of which are easier to track in nonlinearities with lower power.  more » « less
Award ID(s):
2050971
NSF-PAR ID:
10401344
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Nonlinearity
Volume:
36
Issue:
1
ISSN:
0951-7715
Page Range / eLocation ID:
584 to 635
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract This paper is concerned with two classes of cubic quasilinear equations, which can be derived as asymptotic models from shallow-water approximation to the 2D incompressible Euler equations. One class of the models has homogeneous cubic nonlinearity and includes the integrable modified Camassa–Holm (mCH) equation and Novikov equation, and the other class encompasses both quadratic and cubic nonlinearities. It is demonstrated here that both these models possess localized peaked solutions. By constructing a Lyapunov function, these peaked waves are shown to be dynamically stable under small perturbations in the natural energy space $H^1$, without restriction on the sign of the momentum density. In particular, for the homogeneous cubic nonlinear model, we are able to further incorporate a higher-order conservation law to conclude orbital stability in $H^1\cap W^{1,4}$. Our analysis is based on a strong use of the conservation laws, the introduction of certain auxiliary functions, and a refined continuity argument. 
    more » « less
  2. null (Ed.)
    Recent equations of motion for the large deflections of a cantilevered elastic beam are analyzed. In the traditional theory of beam (and plate) large deflections, nonlinear restoring forces are due to the effect of stretching on bending; for an inextensible cantilever, the enforcement of arc-length preservation leads to quasilinear stiffness effects and inertial effects that are both nonlinear and nonlocal. For this model, smooth solutions are constructed via a spectral Galerkin approach. Additional compactness is needed to pass to the limit, and this is obtained through a complex procession of higher energy estimates. Uniqueness is obtained through a non-trivial decomposition of the nonlinearity. The confounding effects of nonlinear inertia are overcome via the addition of structural (Kelvin–Voigt) damping to the equations of motion. Local well-posedness of smooth solutions is shown first in the absence of nonlinear inertial effects, and then shown with these inertial effects present, taking into account structural damping. With damping in force, global-in-time, strong well-posedness result is obtained by achieving exponential decay for small data. 
    more » « less
  3. We study the focusing NLS equation in $R\mathbb{R}^N$ in the mass-supercritical and energy-subcritical (or intercritical ) regime, with $H^1$ data at the mass-energy threshold $\mathcal{ME}(u_0)=\mathcal{ME}(Q)$, where Q is the ground state. Previously, Duyckaerts–Merle studied the behavior of threshold solutions in the $H^1$-critical case, in dimensions $N = 3, 4, 5$, later generalized by Li–Zhang for higher dimensions. In the intercritical case, Duyckaerts–Roudenko studied the threshold problem for the 3d cubic NLS equation. In this paper, we generalize the results of Duyckaerts–Roudenko for any dimension and any power of the nonlinearity for the entire intercritical range. We show the existence of special solutions, $Q^\pm$, besides the standing wave $e^{it}Q$, which exponentially approach the standing wave in the positive time direction, but differ in its behavior for negative time. We classify solutions at the threshold level, showing either blow-up occurs in finite (positive and negative) time, or scattering in both time directions, or the solution is equal to one of the three special solutions above, up to symmetries. Our proof extends to the $H^1$-critical case, thus, giving an alternative proof of the Li–Zhang result and unifying the critical and intercritical cases. These results are obtained by studying the linearized equation around the standing wave and some tailored approximate solutions to the NLS equation. We establish important decay properties of functions associated to the spectrum of the linearized Schrödinger operator, which, in combination with modulational stability and coercivity for the linearized operator on special subspaces, allows us to use a fixed-point argument to show the existence of special solutions. Finally, we prove the uniqueness by studying exponentially decaying solutions to a sequence of linearized equations. 
    more » « less
  4. The thermal radiative transfer (TRT) equations form an integro-differential system that describes the propagation and collisional interactions of photons. Computing accurate and efficient numerical solutions TRT are challenging for several reasons, the first of which is that TRT is defined on a high-dimensional phase space that includes the independent variables of time, space, and velocity. In order to reduce the dimensionality of the phase space, classical approaches such as the P$_N$ (spherical harmonics) or the S$_N$ (discrete ordinates) ansatz are often used in the literature. In this work, we introduce a novel approach: the hybrid discrete (H$^T_N$) approximation to the radiative thermal transfer equations. This approach acquires desirable properties of both P$_N$ and S$_N$, and indeed reduces to each of these approximations in various limits: H$^1_N$ $\equiv$ P$_N$ and H$^T_0$ $\equiv$ S$_T$. We prove that H$^T_N$ results in a system of hyperbolic partial differential equations for all $T\ge 1$ and $N\ge 0$. Another challenge in solving the TRT system is the inherent stiffness due to the large timescale separation between propagation and collisions, especially in the diffusive (i.e., highly collisional) regime. This stiffness challenge can be partially overcome via implicit time integration, although fully implicit methods may become computationally expensive due to the strong nonlinearity and system size. On the other hand, explicit time-stepping schemes that are not also asymptotic-preserving in the highly collisional limit require resolving the mean-free path between collisions, making such schemes prohibitively expensive. In this work we develop a numerical method that is based on a nodal discontinuous Galerkin discretization in space, coupled with a semi-implicit discretization in time. In particular, we make use of a second order explicit Runge-Kutta scheme for the streaming term and an implicit Euler scheme for the material coupling term. Furthermore, in order to solve the material energy equation implicitly after each predictor and corrector step, we linearize the temperature term using a Taylor expansion; this avoids the need for an iterative procedure, and therefore improves efficiency. In order to reduce unphysical oscillation, we apply a slope limiter after each time step. Finally, we conduct several numerical experiments to verify the accuracy, efficiency, and robustness of the H$^T_N$ ansatz and the numerical discretizations. 
    more » « less
  5. In this note we establish hypocoercivity and exponential relaxation to the Maxwellian for a class of kinetic Fokker-Planck-Alignment equations arising in the studies of collective behavior. Unlike previously known results in this direction that focus on convergence near Maxwellian, our result is global for hydrodynamically dense flocks, which has several consequences. In particular, if communication is long-range, the convergence is unconditional. If communication is local then all nearly aligned flocks quantified by smallness of the Fisher information relax to the Maxwellian. In the latter case the class of initial data is stable under the vanishing noise limit, i.e. it reduces to a non-trivial and natural class of traveling wave solutions to the noiseless Vlasov-Alignment equation.

    The main novelty in our approach is the adaptation of a mollified Favre filtration of the macroscopic momentum into the communication protocol. Such filtration has been used previously in large eddy simulations of compressible turbulence and its new variant appeared in the proof of the Onsager conjecture for inhomogeneous Navier-Stokes system. A rigorous treatment of well-posedness for smooth solutions is provided. Lastly, we prove that in the limit of strong noise and local alignment solutions to the Fokker-Planck-Alignment equation Maxwellialize to solutions of the macroscopic hydrodynamic system with the isothermal pressure.

     
    more » « less