Choosing ${\kappa }$ (horizontal ordinate of the saddle point associated to the homoclinic orbit) as bifurcation parameter, bifurcations of the travelling wave solutions is studied in a perturbed $(1 + 1)$-dimensional dispersive long wave equation. The solitary wave solution exists at a suitable wave speed $c$ for the bifurcation parameter ${\kappa }\in \left (0,1-\frac {\sqrt 3}{3}\right )\cup \left (1+\frac {\sqrt 3}{3},2\right )$, while the kink and anti-kink wave solutions exist at a unique wave speed $c^*=\sqrt {15}/3$ for $\kappa =0$ or $\kappa =2$. The methods are based on the geometric singular perturbation (GSP, for short) approach, Melnikov method and invariant manifolds theory. Interestingly, not only the explicit analytical expression of the complicated homoclinic Melnikov integral is directly obtained for the perturbed long wave equation, but also the explicit analytical expression of the limit wave speed is directly given. Numerical simulations are utilized to verify our mathematical results.

Dispersive and Strichartz estimates are obtained for solutions to the wave equation with a Laguerre potential in spatial dimension three. To obtain the desired dispersive estimate, based on the spectral properties of the Schrödinger operator involved, we subsequently prove the dispersive estimate for the corresponding Schrödinger semigroup, obtain a Gaussian-type upper bound, establish Bernstein-type inequalities, and finally pass to the Müller–Seeger’s subordination formula. The desired Strichartz estimates follow by the established dispersive estimate and the standard argument of Keel–Tao.

Let $G$ be a compact Lie group. In this article, we investigate the Cauchy problem for a nonlinear wave equation with the viscoelastic damping on $G$. More precisely, we investigate some $L^2$-estimates for the solution to the homogeneous nonlinear viscoelastic damped wave equation on $G$ utilizing the group Fourier transform on $G$. We also prove that there is no improvement of any decay rate for the norm $\|u(t,\,\cdot )\|_{L^2(G)}$ by further assuming the $L^1(G)$-regularity of initial data. Finally, using the noncommutative Fourier analysis on compact Lie groups, we prove a local in time existence result in the energy space $\mathcal {C}^1([0,\,T],\,H^1_{\mathcal {L}}(G)).$

We establish local-in-time Strichartz estimates for solutions of the model case Dirichlet wave equation inside cylindrical convex domains $\Omega \subset \mathbb {R}^ 3$ with smooth boundary $\partial \Omega \neq \emptyset $ . The key ingredients to prove Strichartz estimates are dispersive estimates, energy estimates, interpolation and $TT^*$ arguments. Strichartz estimates for waves inside an arbitrary domain $\Omega $ have been proved by Blair, Smith and Sogge [‘Strichartz estimates for the wave equation on manifolds with boundary’, Ann. Inst. H. Poincaré Anal. Non Linéaire 26 (2009), 1817–1829]. We provide a detailed proof of the usual Strichartz estimates from dispersive estimates inside cylindrical convex domains for a certain range of the wave admissibility.

We study in $\mathbb {R}^{3+1}$ a system of nonlinearly coupled Klein-Gordon equations under the null condition, with (possibly vanishing) mass varying in the interval $[0, 1]$. Our goal is three-fold, which extends the results in the earlier work of [5, 3]: 1) we want to establish the global well-posedness result to the system that is uniform in terms of the mass parameter (i.e., the smallness of the initial data is independent of the mass parameter); 2) we want to obtain a unified pointwise decay result for the solution to the system, in the sense that the solution decays more like a wave component (independent of the mass parameter) in a certain range of time, while the solution decays as a Klein-Gordon component with a factor depending on the mass parameter in the other part of the time range; 3) the solution to the Klein-Gordon system converges to the solution to the corresponding wave system in a certain sense when the mass parameter goes to 0. In order to achieve these goals, we will rely on both the flat and hyperboloidal foliation of the spacetime and prove a mass-independent $L^2$–type energy estimate for the Klein-Gordon equations with possibly vanishing mass. In addition, the case of the Klein-Gordon equations with certain restricted large data is discussed.

The article studies inverse problems of determining unknown coefficients in various semi-linear and quasi-linear wave equations given the knowledge of an associated source-to-solution map. We introduce a method to solve inverse problems for nonlinear equations using interaction of three waves that makes it possible to study the inverse problem in all globally hyperbolic spacetimes of the dimension $n+1\geqslant 3$ and with partial data. We consider the case when the set $\Omega _{\mathrm{in}}$ , where the sources are supported, and the set $\Omega _{\mathrm{out}}$ , where the observations are made, are separated. As model problems we study both a quasi-linear equation and a semi-linear wave equation and show in each case that it is possible to uniquely recover the background metric up to the natural obstructions for uniqueness that is governed by finite speed of propagation for the wave equation and a gauge corresponding to change of coordinates. The proof consists of two independent components. In the geometric part of the article we introduce a novel geometrical object, the three-to-one scattering relation. We show that this relation determines uniquely the topological, differential and conformal structures of the Lorentzian manifold in a causal diamond set that is the intersection of the future of the point $p_{in}\in \Omega _{\mathrm{in}}$ and the past of the point $p_{out}\in \Omega _{\mathrm{out}}$ . In the analytic part of the article we study multiple-fold linearisation of the nonlinear wave equation using Gaussian beams. We show that the source-to-solution map, corresponding to sources in $\Omega _{\mathrm{in}}$ and observations in $\Omega _{\mathrm{out}}$ , determines the three-to-one scattering relation. The methods developed in the article do not require any assumptions on the conjugate or cut points.

We prove that solutions to the quintic semilinear wave equation with variable coefficients in ${{\mathbb {R}}}^{1+3}$ scatter to a solution to the corresponding linear wave equation. The coefficients are small and decay as $|x|\to \infty$, but are allowed to be time dependent. The proof uses local energy decay estimates to establish the decay of the $L^{6}$ norm of the solution as $t\to \infty$.

Von Neumann’s original proof of the ergodic theorem is revisited. A uniform convergence rate is established under the assumption that one can control the density of the spectrum of the underlying self-adjoint operator when restricted to suitable subspaces. Explicit rates are obtained when the bound is polynomial, with applications to the linear Schrödinger and wave equations. In particular, decay estimates for time averages of solutions are shown.

We consider the nonlinear wave equation (NLW) on the $d$-dimensional torus $\mathbb{T}^{d}$ with a smooth nonlinearity of order at least 2 at the origin. We prove that, for almost any mass, small and smooth solutions of high Sobolev indices are stable up to arbitrary long times with respect to the size of the initial data. To prove this result, we use a normal form transformation decomposing the dynamics into low and high frequencies with weak interactions. While the low part of the dynamics can be put under classical Birkhoff normal form, the high modes evolve according to a time-dependent linear Hamiltonian system. We then control the global dynamics by using polynomial growth estimates for high modes and the preservation of Sobolev norms for the low modes. Our general strategy applies to any semilinear Hamiltonian Partial Differential Equations (PDEs) whose linear frequencies satisfy a very general nonresonance condition. The (NLW) equation on $\mathbb{T}^{d}$ is a good example since the standard Birkhoff normal form applies only when $d=1$ while our strategy applies in any dimension.

Let $(M^{n},g)$ be a Riemannian manifold without boundary. We study the amount of initial regularity required so that the solution to a free Schrödinger equation converges pointwise to its initial data. Assume the initial data is in $H^{\unicode[STIX]{x1D6FC}}(M)$. For hyperbolic space, the standard sphere, and the two-dimensional torus, we prove that $\unicode[STIX]{x1D6FC}>\frac{1}{2}$ is enough. For general compact manifolds, due to the lack of a local smoothing effect, it is hard to improve on the bound $\unicode[STIX]{x1D6FC}>1$ from interpolation. We managed to go below 1 for dimension ${\leqslant}$ 3. The more interesting thing is that, for a one-dimensional compact manifold, $\unicode[STIX]{x1D6FC}>\frac{1}{3}$ is sufficient.

We consider the finite-time blow-up of solutions for the following two kinds of nonlinear wave equation in de Sitter spacetime:

This proof is based on a new blow-up criterion, which generalizes that by Sideris. Furthermore, we give the lifespan estimate of solutions for the problems.

We improve a previous result about the local energy decay for the damped wave equation on $\mathbb{R}^{d}$. The problem is governed by a Laplacian associated with a long-range perturbation of the flat metric and a short-range absorption index. Our purpose is to recover the decay ${\mathcal{O}}(t^{-d+\unicode[STIX]{x1D700}})$ in the weighted energy spaces. The proof is based on uniform resolvent estimates, given by an improved version of the dissipative Mourre theory. In particular, we have to prove the limiting absorption principle for the powers of the resolvent with inserted weights.

In this paper a three-step scheme is applied to solve the Camassa-Holm (CH) shallow water equation. The differential order of the CH equation has been reduced in order to facilitate development of numerical scheme in a comparatively smaller grid stencil. Here a three-point seventh-order spatially accurate upwinding combined compact difference (CCD) scheme is proposed to approximate the firstorder derivative term. We conduct modified equation analysis on the CCD scheme and eliminate the leading discretization error terms for accurately predicting unidirectional wave propagation. The Fourier analysis is carried out as well on the proposed numerical scheme to minimize the dispersive error. For preserving Hamiltonians in Camassa- Holm equation, a symplecticity conserving time integrator has been employed. The other main emphasis of the present study is the use of u–P–α formulation to get nondissipative CH solution for peakon-antipeakon and soliton-anticuspon head-on wave collision problems.

We prove a family of sharp bilinear space–time estimates for the half-wave propagator $\text{e}^{\text{i}t\sqrt{-\unicode[STIX]{x1D6E5}}}$ . As a consequence, for radially symmetric initial data, we establish sharp estimates of this kind for a range of exponents beyond the classical range.

We present a new method for investigating the Lp-type pullback attractors (2 ≤ p < ∞) of a semilinear heat equation on a time-varying domain under quite general assumptions on the nonlinear and forcing terms. The existing approach does not appear applicable here as it is impossible to show the existence of a pullback absorbing set in Lp space when p is large. A new asymptotic decomposition scheme for a non-autonomous pullback attractor has been introduced. The abstract results and preliminary lemmas are also of independent interest and applicable to other systems.

An algebraic multilevel method is proposed for efficiently simulating linear wave propagation using higher-order numerical schemes. This method is used in conjunction with the Finite Volume Time Domain (FVTD) technique for the numerical solution of the time-domain Maxwell’s equations in electromagnetic scattering problems. In the multilevel method the solution is cycled through spatial operators of varying orders of accuracy, while maintaining highest-order accuracy at coarser approximation levels through the use of the relative truncation error as a forcing function. Higher-order spatial accuracy can be enforced using the multilevel method at a fraction of the computational cost incurred in a conventional higher-order implementation. The multilevel method is targeted towards electromagnetic scattering problems at large electrical sizes which usually require long simulation times due to the use of very fine meshes dictated by point-per-wavelength requirements to accurately model wave propagation over long distances.

We couple a node-centered finite volume method to a high order finite difference method to simulate dynamic earthquake ruptures along nonplanar faults in two dimensions. The finite volume method is implemented on an unstructured mesh, providing the ability to handle complex geometries. The geometric complexities are limited to a small portion of the overall domain and elsewhere the high order finite difference method is used, enhancing efficiency. Both the finite volume and finite difference methods are in summation-by-parts form. Interface conditions coupling the numerical solution across physical interfaces like faults, and computational ones between structured and unstructured meshes, are enforced weakly using the simultaneous-approximation-term technique. The fault interface condition, or friction law, provides a nonlinear relation between fields on the two sides of the fault, and allows for the particle velocity field to be discontinuous across it. Stability is proved by deriving energy estimates; stability, accuracy, and efficiency of the hybrid method are confirmed with several computational experiments. The capabilities of the method are demonstrated by simulating an earthquake rupture propagating along the margins of a volcanic plug.

A new weak boundary procedure for hyperbolic problems is presented. We consider high order finite difference operators of summation-by-parts form with weak boundary conditions and generalize that technique. The new boundary procedure is applied near boundaries in an extended domain where data is known. We show how to raise the order of accuracy of the scheme, how to modify the spectrum of the resulting operator and how to construct non-reflecting properties at the boundaries. The new boundary procedure is cheap, easy to implement and suitable for all numerical methods, not only finite difference methods, that employ weak boundary conditions. Numerical results that corroborate the analysis are presented.

In this paper, we investigate the existence of global weak solutions to an integrable two-component Camassa–Holm shallow-water system, provided the initial data u0(x) and ρ0(x) have end states u± and ρ±, respectively. By perturbing the Cauchy problem of the system around rarefaction waves of the well-known Burgers equation, we obtain a global weak solution for the system under the assumptions u− ≤ u+ and ρ− ≤ ρ+.

This paper is concerned with the initial boundary value problem of a class of nonlinear wave equations and reaction–diffusion equations with several nonlinear source terms of different signs. For the initial boundary value problem of the nonlinear wave equations, we derive a blow up result for certain initial data with arbitrary positive initial energy. For the initial boundary value problem of the nonlinear reaction–diffusion equations, we discuss some probabilities of the existence and nonexistence of global solutions and give some sufficient conditions for the global and nonglobal existence of solutions at high initial energy level by employing the comparison principle and variational methods.