To send content items to your account,
please confirm that you agree to abide by our usage policies.
If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your account.
Find out more about sending content to .
To send content items to your Kindle, first ensure firstname.lastname@example.org
is added to your Approved Personal Document E-mail List under your Personal Document Settings
on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part
of your Kindle email address below.
Find out more about sending to your Kindle.
Note you can select to send to either the @free.kindle.com or @kindle.com variations.
‘@free.kindle.com’ emails are free but can only be sent to your device when it is connected to wi-fi.
‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.
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
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
is a good example since the standard Birkhoff normal form applies only when
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 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 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 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
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.
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
. For hyperbolic space, the standard sphere, and the two-dimensional torus, we prove that
is enough. For general compact manifolds, due to the lack of a local smoothing effect, it is hard to improve on the bound
from interpolation. We managed to go below 1 for dimension
3. The more interesting thing is that, for a one-dimensional compact manifold,
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
. 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.
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.
A two-layer non-hydrostatic numerical model is proposed to simulate the formation of undular bores by tsunami wave fission. These phenomena could not be produced by a hydrostatic model. Here, we derived a modified Shallow Water Equations with involving hydrodynamic pressure using two layer approach. Staggered finite volume method with predictor corrector step is applied to solve the equation numerically. Numerical dispersion relation is derived from our model to confirm the exact linear dispersion relation for dispersive waves. To illustrate the performance of our non-hydrostatic scheme in case of linear wave dispersion and non-linearity, four test cases of free surface flows are provided. The first test case is standing wave in a closed basin, which test the ability of the numerical scheme in simulating dispersive wave motion with the correct frequency. The second test case is the solitary wave propagation as the examination of owing balance between dispersion and nonlinearity. Regular wave propagation over a submerged bar test by Beji is simulated to show that our non-hydrostatic scheme described well the shoaling process as well as the linear dispersion compared with the experimental data. The last test case is the undular bore propagation.
We describe a numerical model to simulate the non-linear elasto-plastic dynamics of compressible materials. The model is fully Eulerian and it is discretized on a fixed Cartesian mesh. The hyperelastic constitutive law considered is neohookean and the plasticity model is based on a multiplicative decomposition of the inverse deformation tensor. The model is thermodynamically consistent and it is shown to be stable in the sense that the norm of the deviatoric stress tensor beyond yield is non increasing. The multimaterial integration scheme is based on a simple numerical flux function that keeps the interfaces sharp. Numerical illustrations in one to three space dimensions of high-speed multimaterial impacts in air are presented.
For 2D elastic-plastic flows with the hypo-elastic constitutive model and von Mises’ yielding condition, the non-conservative character of the hypo-elastic constitutive model and the von Mises’ yielding condition make the construction of the solution to the Riemann problem a challenging task. In this paper, we first analyze the wave structure of the Riemann problem and develop accordingly a Four-Rarefaction wave approximate Riemann Solver with Elastic waves (FRRSE). In the construction of FRRSE one needs to use an iterative method. A direct iteration procedure for four variables is complex and computationally expensive. In order to simplify the solution procedure we develop an iteration based on two nested iterations upon two variables, and our iteration method is simple in implementation and efficient. Based on FRRSE as a building block, we propose a 2nd-order cell-centered Lagrangian numerical scheme. Numerical results with smooth solutions show that the scheme is of second-order accuracy. Moreover, a number of numerical experiments with shock and rarefaction waves demonstrate the scheme is essentially non-oscillatory and appears to be convergent. For shock waves the present scheme has comparable accuracy to that of the scheme developed by Maire et al., while it is more accurate in resolving rarefaction waves.