BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Jean-François Coulombel appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We prove a sharp uniform generalized Gaussian bound for the Green's function of the Lax-Wendroff and Beam-Warming schemes. Our bound highlights the spatial region that leads to the well-known (rather weak) instability of these schemes in the maximum norm. We also recover uniform bounds in the maximum norm when these schemes are applied to initial data of bounded variation.
In this article, we consider a class of finite rank perturbations of Toeplitz operators that have simple eigenvalues on the unit circle. Under a suitable assumption on the behavior of the essential spectrum, we show that such operators are power bounded. The problem originates in the approximation of hyperbolic partial differential equations with boundary conditions by means of finite difference schemes. Our result gives a positive answer to a conjecture by Trefethen, Kreiss and Wu that only a weak form of the so-called Uniform Kreiss-Lopatinskii Condition is sufficient to imply power boundedness.
We explore in this article the possibilities and limitations of the so-called energy method for analyzing the stability of finite difference approximations to the transport equation with extrapolation numerical boundary conditions at the outflow boundary. We first show that for the most simple schemes, namely the explicit schemes with a three point stencil, the energy method can be applied for proving stability estimates when the scheme is implemented with either the first or second order extrapolation boundary condition. We then examine the case of five point stencils and give several examples of schemes and second order extrapolation numerical boundary conditions for which the energy method produces stability estimates. However, we also show that for the standard first or second order translatory extrapolation boundary conditions, the energy method cannot be applied for proving stability of the classical fourth order scheme originally proposed by Strang. This gives a clear limitation of the energy method with respect to the more general approach based on the normal mode decomposition.
The aim of this article is to propose a systematic study of transparent boundary conditions for finite difference approximations of evolution equations. We try to keep the discussion at the highest level of generality in order to apply the theory to the broadest class of problems. We deal with two main issues. We first derive transparent numerical boundary conditions, that is, we exhibit the relations satisfied by the solution to the pure Cauchy problem when the initial condition vanishes outside of some domain. Our derivation encompasses discretized transport, diffusion and dispersive equations with arbitrarily wide stencils. The second issue is to prove sharp stability estimates for the initial boundary value problem obtained by enforcing the boundary conditions derived in the first step. We focus here on discretized transport equations. Under the assumption that the numerical boundary is non-characteristic, our main result characterizes the class of numerical schemes for which the corresponding transparent boundary conditions satisfy the so-called Uniform Kreiss-Lopatinskii Condition introduced in [GKS72]. Adapting some previous works to the non-local boundary conditions considered here, our analysis culminates in the derivation of trace and semigroup estimates for such transparent numerical boundary conditions. Several examples and possible extensions are given.
Among hyperbolic Initial Boundary Value Problems (IBVP), those coming from a variational principle 'generically' admit linear surface waves, as was shown by Serre [J. Funct. Anal. 2006]. At the weakly nonlinear level, the behavior of surface waves is expected to be governed by an amplitude equation that can be derived by means of a formal asymptotic expansion. Amplitude equations for weakly nonlinear surface waves were introduced by Lardner [Int. J. Engng Sci. 1983], Parker and co-workers [J. Elasticity 1985] in the framework of elasticity, and by Hunter [Contemp. Math. 1989] for abstract hyperbolic problems. They consist of nonlocal evolution equations involving a complicated, bilinear Fourier multiplier in the direction of propagation along the boundary. It was shown by the authors in an earlier work [Arch. Ration. Mech. Anal. 2012] that this multiplier, or kernel, inherits some algebraic properties from the original IBVP. These properties are crucial for the (local) well-posedness of the amplitude equation, as shown together with Tzvetkov [Adv. Math., 2011]. Properties of amplitude equations are revisited here in a somehow simpler way, for surface waves in a variational setting. Applications include various physical models, from elasticity of course to the director-field system for liquid crystals introduced by Saxton [Contemp. Math. 1989] and studied by Austria and Hunter [Commun. Inf. Syst. 2013]. Similar properties are eventually shown for the amplitude equation associated with surface waves at reversible phase boundaries in compressible fluids, thus completing a work initiated by Benzoni-Gavage and Rosini [Comput. Math. Appl. 2009].
This technical note is a complement to an earlier paper [Benzoni-Gavage \& Rosini, Comput. Math. Appl. 2009], which aims at a deeper understanding of a basic model for propagating phase boundaries that was proved to admit surface waves [Benzoni-Gavage, Nonlinear Anal. 1998]. The amplitude equation governing the evolution of weakly nonlinear surface waves for that model is computed explicitly, and is eventually found to have enough symmetry properties for the associated Cauchy problem to be locally well-posed.
Leray and Gårding have developed a multiplier technique for deriving a priori estimates for solutions to scalar hyperbolic equations in either the whole space or the torus. In particular, the arguments in Leray and Gårding's work provide with at least one local multiplier and one local energy functional that is controlled along the evolution. The existence of such a local multiplier is the starting point of the argument by Rauch for the derivation of semigroup estimates for hyperbolic initial boundary value problems. In this article, we explain how this multiplier technique can be adapted to the framework of finite difference approximations of transport equations. The technique applies to numerical schemes with arbitrarily many time levels, and encompasses a somehow magical trick that has been known for a long time for the leapfrog scheme. More importantly, the existence and properties of the local multiplier enable us to derive optimal semigroup estimates for fully discrete hyperbolic initial boundary value problems, which answers a problem raised by Trefethen, Kreiss and Wu.
The stability theory for hyperbolic initial boundary value problems relies most of the time on the Laplace transform with respect to the time variable. For technical reasons, this usually restricts the validity of stability estimates to the case of zero initial data. In this article, we consider the class of non-glancing finite difference approximations to the hyperbolic operator. We show that the maximal stability estimates that are known for zero initial data and nonzero boundary source term extend to the case of nonzero initial data in â 2 . The main novelty of our approach is to cover finite difference schemes with an arbitrary number of time levels. As an easy corollary of our main trace estimate, we recover former stability results in the semigroup sense by Kreiss [Kre68] and Osher [Osh69b].