The study of spherically symmetric motion is important for the theory of explosion waves. In this paper, we consider a ‘spherical piston’ problem for the relativistic Euler equations, which describes the wave motion produced by a sphere expanding into an infinite surrounding medium. We use the reflected characteristics method to construct a global piecewise smooth solution with a single shock of this spherical piston problem, provided that the speed of the sphere is a small perturbation of a constant speed.

In this note, we study the hyperbolic stochastic damped sine-Gordon equation (SdSG), with a parameter β2 > 0, and its associated Gibbs dynamics on the two-dimensional torus. After introducing a suitable renormalization, we first construct the Gibbs measure in the range 0 < β2 < 4π via the variational approach due to Barashkov-Gubinelli (2018). We then prove almost sure global well-posedness and invariance of the Gibbs measure under the hyperbolic SdSG dynamics in the range 0 < β2 < 2π. Our construction of the Gibbs measure also yields almost sure global well-posedness and invariance of the Gibbs measure for the parabolic sine-Gordon model in the range 0 < β2 < 4π.

We prove that any simple planar travelling wave solution to the membrane equation in spatial dimension $d\geqslant 3$ with bounded spatial extent is globally nonlinearly stable under sufficiently small compactly supported perturbations, where the smallness depends on the size of the support of the perturbation as well as on the initial travelling wave profile. The main novelty of the argument is the lack of higher order peeling in our vector-field-based method. In particular, the higher order energies (in fact, all energies at order $2$ or higher) are allowed to grow polynomially (but in a controlled way) in time. This is in contrast with classical global stability arguments, where only the ‘top’ order energies used in the bootstrap argument exhibit growth, and reflects the fact that the background travelling wave solution has ‘infinite energy’ and the coefficients of the perturbation equation are not asymptotically Lorentz invariant. Nonetheless, we can prove that the perturbation converges to zero in $C^{2}$ by carefully analysing the nonlinear interactions and exposing a certain ‘vestigial’ null structure in the equations.

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.

The aim of the paper is to introduce and investigate a dynamical system which consists of a variational–hemivariational inequality of hyperbolic type combined with a non-linear evolution equation. Such a dynamical system arises in studies of complicated contact problems in mechanics. Existence, uniqueness and regularity of a global solution to the system are established. The approach is based on a new semi-discrete approximation with an application of a surjectivity result for a pseudomonotone perturbation of a maximal monotone operator. A new dynamic viscoelastic frictional contact model with adhesion is studied as an application, in which the contact boundary condition is described by a generalised normal damped response condition with unilateral constraint and a multivalued frictional contact law.

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.

We present a model for a class of non-local conservation laws arising in traffic flow modelling at road junctions. Instead of a single velocity function for the whole road, we consider two different road segments, which may differ for their speed law and number of lanes (hence their maximal vehicle density). We use an upwind type numerical scheme to construct a sequence of approximate solutions, and we provide uniform L∞ and total variation estimates. In particular, the solutions of the proposed model stay positive and below the maximum density of each road segment. Using a Lax–Wendroff type argument and the doubling of variables technique, we prove the well-posedness of the proposed model. Finally, some numerical simulations are provided and compared with the corresponding (discontinuous) local model.

The numerical entropy production (NEP) for shallow water equations (SWE) is discussed and implemented as a smoothness indicator. We consider SWE in three different dimensions, namely, one-dimensional, one-and-a-half-dimensional, and two-dimensional SWE. An existing numerical entropy scheme is reviewed and an alternative scheme is provided. We prove the properties of these two numerical entropy schemes relating to the entropy steady state and consistency with the entropy equality on smooth regions. Simulation results show that both schemes produce NEP with the same behaviour for detecting discontinuities of solutions and perform similarly as smoothness indicators. An implementation of the NEP for an adaptive numerical method is also demonstrated.

The study of radially symmetric motion is important for the theory of explosion waves. We construct rigorously self-similar entropy solutions to Riemann initial-boundary value problems for the radially symmetric relativistic Euler equations. We use the assumption of self-similarity to reduce the relativistic Euler equations to a system of nonlinear ordinary differential equations, from which we obtain detailed structures of solutions besides their existence. For the ultra-relativistic Euler equations, we also obtain the uniqueness of the self-similar entropy solution to the Riemann initial-boundary value problems.

The present paper concerns the system ut + [ϕ(u)]x = 0, vt + [ψ(u)v]x = 0 having distributions as initial conditions. Under certain conditions, and supposing ϕ, ψ: ℝ → ℝ functions, we explicitly solve this Cauchy problem within a convenient space of distributions u,v. For this purpose, a consistent extension of the classical solution concept defined in the setting of a distributional product (not constructed by approximation processes) is used. Shock waves, δ-shock waves, and also waves defined by distributions that are not measures are presented explicitly as examples. This study is carried out without assuming classical results about conservation laws. For reader's convenience, a brief survey of the distributional product is also included.

This paper is concerned with the periodic (in time) solutions to an one-dimensional semilinear wave equation with x-dependent coefficients. Such a model arises from the forced vibrations of a nonhomogeneous string and propagation of seismic waves in nonisotropic media. By combining variational methods with saddle point reduction technique, we obtain the existence of at least three periodic solutions whenever the period is a rational multiple of the length of the spatial interval. Our method is based on a delicate analysis for the asymptotic character of the spectrum of the wave operator with x-dependent coefficients, and the spectral properties play an essential role in the proof.

We investigate the Cauchy problem of the viscous liquid-gas two-phase flow model in ℝ3. Under the assumption that the initial data is close to the constant equilibrium state in the framework of Sobolev space H2(ℝ3), the Cauchy problem is shown to be globally well-posed by an energy method. If additionally, for 1 ⩽ p < 6/5, Lp-norm of the initial perturbation is bounded, the optimal convergence rates of the solutions in Lq-norm with 2 ⩽ q ⩽ 6 and optimal convergence rates of their spatial derivatives in L2-norm are also obtained by combining spectral analysis with energy methods.

We demonstrate the global existence of weak solutions to a class of semilinear strongly damped wave equations possessing nonlinear hyperbolic dynamic boundary conditions. The associated linear operator is $(-\unicode[STIX]{x1D6E5}_{W})^{\unicode[STIX]{x1D703}}\unicode[STIX]{x2202}_{t}u$, where $\unicode[STIX]{x1D703}\in [\frac{1}{2},1)$ and $\unicode[STIX]{x1D6E5}_{W}$ is the Wentzell–Laplacian. A balance condition is assumed to hold between the nonlinearity defined on the interior of the domain and the nonlinearity on the boundary. This allows for arbitrary (supercritical) polynomial growth of each potential, as well as mixed dissipative/antidissipative behaviour.

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 prove the stability with respect to the flux of solutions to initial – boundary value problems for scalar non autonomous conservation laws in one space dimension. Key estimates are obtained through a careful construction of the solutions.

In this article we are interested in the rigorous construction of WKB expansions for hyperbolic boundary value problems in the strip $\mathbb{R}^{d-1}\times [0,1]$. In this geometry, a new inversibility condition has to be imposed to construct the WKB expansion. This new condition is due to selfinteraction phenomenon which naturally appear when several boundary conditions are imposed. More precisely, by selfinteraction we mean that some rays can regenerated themselves after some rebounds against the sides of the strip. This phenomenon is not new and has already been studied in Benoit (Geometric optics expansions for hyperbolic corner problems, I: self-interaction phenomenon, Anal. PDE9(6) (2016), 1359–1418), Sarason and Smoller (Geometrical optics and the corner problem, Arch. Rat. Mech. Anal.56 (1974/75), 34–69) in the corner geometry. In this framework the existence of such selfinteracting rays is linked to specific geometries of the characteristic variety and may seem to be somewhat anecdotal. However for the strip geometry such rays become generic. The new inversibility condition, used to construct the WKB expansion, is a microlocalized version of the one characterizing the uniform in time strong well-posedness (Benoit, Lower exponential strong well-posedness of hyperbolic boundary value problems in a strip (preprint)). It is interesting to point here that such a situation already occurs in the half space geometry (Kreiss, Initial boundary value problems for hyperbolic systems, Comm. Pure Appl. Math.23 (1970), 277–298).

In this paper, we study the initial boundary value problem for a class of fourth order damped wave equations with arbitrary positive initial energy. In the framework of the energy method, we further exploit the properties of the Nehari functional. Finally, the global existence and finite time blow-up of solutions are obtained.

For bounded domains Ω, we prove that the Lp-norm of a regular function with compact support is controlled by weighted Lp-norms of its gradient, where the weight belongs to a class of symmetric non-negative definite matrix-valued functions. The class of weights is defined by regularity assumptions and structural conditions on the degeneracy set, where the determinant vanishes. In particular, the weight A is assumed to have rank at least 1 when restricted to the normal bundle of the degeneracy set S. This generalization of the classical Poincaré inequality is then applied to develop a robust theory of first-order Lp-based Sobolev spaces with matrix-valued weight A. The Poincaré inequality and these Sobolev spaces are then applied to produce various results on existence, uniqueness and qualitative properties of weak solutions to boundary-value problems for degenerate elliptic, degenerate parabolic and degenerate hyperbolic partial differential equations (PDEs) of second order written in divergence form, where A is calibrated to the matrix of coefficients of the second-order spatial derivatives. The notion of weak solution is variational: the spatial states belong to the matrix-weighted Sobolev spaces with p = 2. For the degenerate elliptic PDEs, the Dirichlet problem is treated by the use of the Poincaré inequality and Lax–Milgram theorem, while the treatment of the Cauchy–Dirichlet problem for the degenerate evolution equations relies only on the Poincaré inequality and the parabolic and hyperbolic counterparts of the Lax–Milgram theorem.

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.

In this study, we propose a high order well-balanced weighted compact nonlinear (WCN) scheme for the gas dynamic equations under gravitational fields. The proposed scheme is an extension of the high order WCN schemes developed in (S. Zhang, S. Jiang, C.-W Shu, J. Comput. Phys. 227 (2008) 7294-7321). For the purpose of maintaining the exact steady state solution, the well-balanced technique in (Y. Xing, C.-W Shu, J. Sci. Comput. 54 (2013) 645-662) is employed to split the source term into two terms. The proposed scheme can maintain the isothermal equilibrium solution exactly, genuine high order accuracy and resolve small perturbations of the hydrostatic balance state on the coarse meshes. Furthermore, in order to capture the strong discontinuities and large gradients, the fifth-order upwind weighted nonlinear interpolations together with the fourth/sixth order cell-centered compact schemes with local characteristic projections are used to construct different WCN schemes. Several representative one- and two-dimensional examples are simulated to demonstrate the good performance of the proposed schemes.