Long time dynamics of the smoothed step initial value problem or dispersive Riemann problem for the Benjamin‐Bona‐Mahony (BBM) equation
We consider a higher‐dimensional version of the Benjamin‐Ono (HBO) equation in the 2D setting:
- Award ID(s):
- 1927258
- NSF-PAR ID:
- 10362243
- Publisher / Repository:
- Wiley-Blackwell
- Date Published:
- Journal Name:
- Studies in Applied Mathematics
- Volume:
- 148
- Issue:
- 2
- ISSN:
- 0022-2526
- Page Range / eLocation ID:
- p. 498-542
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract are studied using asymptotic methods and numerical simulations. The catalog of solutions of the dispersive Riemann problem for the BBM equation is much richer than for the related, integrable, Korteweg‐de Vries equation . The transition width of the initial smoothed step is found to significantly impact the dynamics. Narrow width gives rise to rarefaction and dispersive shock wave (DSW) solutions that are accompanied by the generation of two‐phase linear wavetrains, solitary wave shedding, and expansion shocks. Both narrow and broad initial widths give rise to two‐phase nonlinear wavetrains or DSW implosion and a new kind of dispersive Lax shock for symmetric data. The dispersive Lax shock is described by an approximate self‐similar solution of the BBM equation whose limit as is a stationary, discontinuous weak solution. By introducing a slight asymmetry in the data for the dispersive Lax shock, the generation of an incoherent solitary wavetrain is observed. Further asymmetry leads to the DSW implosion regime that is effectively described by a pair of coupled nonlinear Schrödinger equations. The complex interplay between nonlocality, nonlinearity, and dispersion in the BBM equation underlies the rich variety of nonclassical dispersive hydrodynamic solutions to the dispersive Riemann problem. -
more » « less
Abstract We consider the mapping properties of the integral operator arising in nonlocal slender body theory (SBT) for the model geometry of a straight, periodic filament. It is well known that the classical singular SBT integral operator suffers from high wavenumber instabilities, making it unsuitable for approximating the
slender body inverse problem , where the fiber velocity is prescribed and the integral operator must be inverted to find the force density along the fiber. Regularizations of the integral operator must therefore be used instead. Here, we consider two regularization methods: spectral truncation and the‐regularization of Tornberg and Shelley (2004). We compare the mapping properties of these approximations to the underlying partial differential equation (PDE) solution, which for the inverse problem is simply the Stokes Dirichlet problem with data constrained to be constant on cross sections. For the straight‐but‐periodic fiber with constant radius , we explicitly calculate the spectrum of the operator mapping fiber velocity to force for both the PDE and the approximations. We prove that the spectrum of the original SBT operator agrees closely with the PDE operator at low wavenumbers but differs at high frequencies, allowing us to define a truncated approximation with a wavenumber cutoff . For both the truncated and ‐regularized approximations, we obtain rigorous ‐based convergence to the PDE solution as : A fiber velocity with regularity gives convergence, while a fiber velocity with at least regularity yields convergence. Moreover, we determine the dependence of the ‐regularized error estimate on the regularization parameter . -
more » « less
Abstract The Whitham equation was proposed as a model for surface water waves that combines the quadratic flux nonlinearity
of the Korteweg–de Vries equation and the full linear dispersion relation of unidirectional gravity water waves in suitably scaled variables. This paper proposes and analyzes a generalization of Whitham's model to unidirectional nonlinear wave equations consisting of a general nonlinear flux function and a general linear dispersion relation . Assuming the existence of periodic traveling wave solutions to this generalized Whitham equation, their slow modulations are studied in the context of Whitham modulation theory. A multiple scales calculation yields the modulation equations, a system of three conservation laws that describe the slow evolution of the periodic traveling wave's wavenumber, amplitude, and mean. In the weakly nonlinear limit, explicit, simple criteria in terms of general and establishing the strict hyperbolicity and genuine nonlinearity of the modulation equations are determined. This result is interpreted as a generalized Lighthill–Whitham criterion for modulational instability. -
more » « less
Abstract In this work, we investigate the existence of nonmonotone traveling wave solutions to a reaction‐diffusion system modeling social outbursts, such as rioting activity, originally proposed in Berestycki et al (
Netw Heterog Media . 2015;10(3):443–475). The model consists of two scalar values, the level of unrestand a tension field . A key component of the model is a bandwagon effect in the unrest, provided the tension is sufficiently high. We focus on the so‐called tension‐inhibitive regime, characterized by the fact that the level of unrest has a negative feedback on the tension. This regime has been shown to be physically relevant for the spatiotemporal spread of the 2005 French riots. We use Geometric Singular Perturbation Theory to study the existence of such solutions in two situations. The first is when both and diffuse at a very small rate. Here, the time scale over which the bandwagon effect is observed plays a key role. The second case we consider is when the tension diffuses at a much slower rate than the level of unrest. In this case, we are able to deduce that the driving dynamics are modeled by the well‐known Fisher–Kolmogorov‐Petrovsky‐Piskunov (KPP) equation. -
more » « less
Abstract Weakly nonlinear, bi‐periodic patterns of waves that propagate in the
‐direction with amplitude variation in the ‐direction are generated in a laboratory. The amplitude variation in the ‐direction is studied within the framework of the vector (vNLSE) and scalar (sNLSE) nonlinear Schrödinger equations using the uniform‐amplitude, Stokes‐like solution of the vNLSE and the Jacobi elliptic sine function solution of the sNLSE. The wavetrains are generated using the Stokes‐like solution of vNLSE; however, a comparison of both predictions shows that while they both do a reasonably good job of predicting the observed amplitude variation in , the comparison with the elliptic function solution of the sNLSE has significantly less error when the ratio of ‐wavenumber to the two‐dimensional wavenumber is less than about 0.25. For ratios between about 0.25 and 0.30 (the limit of the experiments), the two models have comparable errors. When the ratio is less than about 0.17, agreement with the vNLSE solution requires a third‐harmonic term in the ‐direction, obtained from a Stokes‐type expansion of interacting, symmetric wavetrains. There is no evidence of instability growth in the ‐direction, consistent with the work of Segur and colleagues, who showed that dissipation stabilizes the modulational instability. Finally, there is some extra amplitude variation in , which is examined via a qualitative stability calculation that allows symmetry breaking in that direction.
