BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Philippe G. LeFloch

Philippe G. LeFloch appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.APgr-qcmath.NAmath.DGNumerical Analysismath.CA

Catalog footprint

What is connected

58works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

58 published item(s)

preprint2023arXiv

Optimal localization for the Einstein constraints

We consider asymptotically Euclidean, initial data sets for Einstein's field equations and solve the localization problem at infinity, also called gluing problem. We achieve optimal gluing and optimal decay, in the sense that we encompass solutions with possibly arbitrarily low decay at infinity and establish (super-)harmonic estimates within possibly arbitrarily narrow conical domains. In the localized seed-to-solution method (as we call it), we define a variational projection operator which associates the solution to the Einstein constraints that is closest to any given localized seed data set (as we call it). Our main contribution concerns the derivation of harmonic estimates for the linearized Einstein operator and its formal adjoint which, in particular, includes new analysis on the linearized scalar curvature operator. The statement of harmonic estimates requires the notion of energy-momentum modulators (as we call them), which arise as correctors to the localized seed data sets. For the Hamiltonian and momentum operators, we introduce a notion of harmonic-spherical decomposition and we uncover stability conditions on the localization function, which are localized Poincare and Hardy-type inequalities and, for instance, hold for arbitrarily narrow gluing domains. Our localized seed-to-solution method builds upon the gluing techniques pioneered by Carlotto, Chrusciel, Corvino, Delay, Isenberg, Maxwell, and Schoen, while providing a proof of a conjecture by Carlotto and Schoen on the localization problem and generalize P. LeFloch and Nguyen's theorem on the asymptotic localization problem.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Singularity scattering laws for bouncing cosmologies: a brief overview

For contracting/expanding bouncing cosmologies, the formulation of junction conditions at a bouncing was recently revisited by the author in collaboration with B. Le Floch and G. Veneziano. The regime of interest here is the so-called quiescent regime, in which a monotone behavior of the metric is observed and asymptotic expansions can be derived. Here, we overview our new methodology based on the notion of {\sl singularity scattering maps} and {\sl cyclic spacetimes}, and we present our main conclusions. In particular, we provide a classification of all allowed bouncing junction conditions, including three universal laws.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Formulation and convergence of the finite volume method for conservation laws on spacetimes with boundary

We study nonlinear hyperbolic conservation laws posed on a differential (n+1)-manifold with boundary referred to as a spacetime, and defined from a prescribed flux field of n-forms depending on a parameter (the unknown variable), a class of equations proposed by LeFloch and Okutmustur in 2008. Our main result is a proof of the convergence of the finite volume method for weak solutions satisfying suitable entropy inequalities. A main difference with previous work is that we allow for slices with a boundary and, in addition, introduce a new formulation of the finite volume method involving the notion of total flux functions. Under a natural global hyperbolicity condition on the flux field and the spacetime and by assuming that the spacetime admits a foliation by compact slices with boundary, we establish an existence and uniqueness theory for the initial and boundary value problem, and we prove a contraction property in a geometrically natural L1-type distance.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Compensated compactness and corrector stress tensor for the Einstein equations in T2 symmetry

We consider the Einstein equations in T2 symmetry, either for vacuum spacetimes or coupled to the Euler equations for a compressible fluid, and we introduce the notion of T2 areal flows on T3 with finite total energy. By uncovering a hidden structure of the Einstein equations, we establish a compensated compactness framework and solve the global evolution problem for vacuum spacetimes as well as for self-gravitating compressible fluids. We study the stability and instability of such flows and prove that, when the initial data are well-prepared, any family of T2 areal flows is sequentially compact in a natural topology. In order to handle general initial data we propose a relaxed notion of T2 areal flows endowed with a corrector stress tensor (as we call it) which is a bounded measure generated by geometric oscillations and concentrations propagating at the speed of light. This generalizes a result for vacuum spacetimes in: Le Floch B. and LeFloch P.G., Arch. Rational Mech. Anal. 233 (2019), 45-86. In addition, we determine the global geometry of the corresponding future Cauchy developments and we prove that the area of the T2 orbits generically approaches infinity in the future-expanding regime. In the future-contracting regime, the volume of the T3 spacelike slices approaches zero and, for generic initial data, the area of the orbits of symmetry approaches zero in Gowdy symmetric matter spacetimes and in T2 vacuum spacetimes.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Global evolution of the U(1) Higgs Boson: nonlinear stability and uniform energy bounds

Relying on the hyperboloidal foliation method, we establish the nonlinear stability of the ground state of the $U(1)$ standard model of electroweak interactions. This amounts to establishing a global-in-time theory for the initial value problem for a nonlinear wave-Klein-Gordon system that couples (Dirac, scalar, gauge) massive equations together. In particular, we investigate here the Dirac equation and consider a new energy functional for this field defined with respect to the hyperboloidal foliation of Minkowski spacetime. We provide a novel decay result for the Dirac equation which is uniform in the mass coefficient, and thus allows for the Dirac mass coefficient to be arbitrarily small. Furthermore we obtain energy bounds for the Higgs fields and gauge bosons that are uniform with respect to the hyperboloidal time variable.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Mesh-free error integration in arbitrary dimensions: a numerical study of discrepancy functions

We are interested in mesh-free formulas based on the Monte-Carlo methodology for the approximation of multi-dimensional integrals, and we investigate their accuracy when the functions belong to a reproducing-kernel space. A kernel typically captures regularity and qualitative properties of functions "beyond" the standard Sobolev regularity class. We are interested in the issue whether quantitative error bounds can be a priori guaranteed in applications (e.g. mathematical finance but also scientific computing and machine learning). Our main contribution is a numerical study of the error discrepancy function based on a comparison between several numerical strategies, when one varies the choice of the kernel, the number of approximation points, and the dimension of the problem. We consider two strategies in order to localize to a bounded set the standard kernels defined in the whole Euclidian space (exponential, multiquadric, Gaussian, truncated), namely, on one hand the class of periodic kernels defined via a discrete Fourier transform on a lattice and, on the other hand, a class of transport-based kernels. First of all, relying on a Poisson formula on a lattice, together with heuristic arguments, we discuss the derivation of theoretical bounds for the discrepancy function of periodic kernels. Second, for each kernel of interest, we perform the numerical experiments that are required in order to generate the optimal distributions of points and the discrepancy error functions. Our numerical results allow us to validate our theoretical observations and provide us with quantitative estimates for the error made with a kernel-based strategy as opposed to a purely random strategy.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

An algorithm (CoDeFi) for overcoming the curse of dimensionality in mathematical finance

We present an algorithm (CoDeFi) which overcomes the curse of dimensionality (CoD) in scientific computations and, especially, in mathematical finance (Fi). Our method applies a broad class of partial differential equations such as Kolmogorov-type equations and, for instance, the Black and Scholes equation. As a main feature, our algorithm allows one to solve partial differential equations in large dimensions and provides a general framework for stochastic problems. In insurance or finance applications, the number of dimensions corresponds to the number of risk sources and it is crucial to have a numerical method that remains robust and reliable in large dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Boundedness of the total energy of relativistic membranes evolving in a curved spacetime

We establish a global existence theory for the equation governing the evolution of a relativistic membrane in a (possibly curved) Lorentzian manifold, when the spacetime metric is a perturbation of the Minkowski metric. Relying on the Hyperboloidal Foliation Method introduced by LeFloch and Ma in 2014, we revisit a theorem established earlier by Lindblad (who treated membranes in the flat Minkowski spacetime) and we provide a simpler proof of existence, which is also valid in a curved spacetime and, most importantly, leads to the important property that the total energy of the membrane is globally bounded in time.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Future asymptotics and geodesic completeness of polarized T2-symmetric spacetimes

We investigate the late-time asymptotics of future expanding, polarized vacuum Einstein spacetimes with T2-symmetry on T3, which, by definition, admit two spacelike Killing fields. Our main result is the existence of a stable asymptotic regime within this class, that is, we provide here a full description of the late-time asymptotics of the solutions to the Einstein equations when the initial data set is close to the asymptotic regime. Our proof is based on several energy functionals with lower order corrections (as is standard for such problems) and the derivation of a simplified model which we exhibit here. Roughly speaking, the Einstein equations in the symmetry class under consideration consists of a system of wave equations coupled to constraint equations plus a system of ordinary differential equations. The unknowns involved in the system of ordinary equations are blowing up in the future timelike directions. One of our main contributions is the derivation of novel effective equations for suitably renormalized unknowns. Interestingly, this renormalization is not performed with respect to a fixed background, but does involve the energy of the coupled system of wave equations. In addition, we construct an open set of initial data which are arbitrarily close to the expected asymptotic behavior. We emphasize that, in comparison, the class of Gowdy spacetimes exhibits a very different dynamical behavior to the one we uncover in the present work for general polarized T2-symmetric spacetimes. Furthermore, all the conclusions of this paper are valid within the framework of weakly T2-symmetric spacetimes previously introduced by the authors.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Nonlocal conservation laws. I. A new class of monotonicity-preserving models

We introduce a new class of nonlocal nonlinear conservation laws in one space dimension that allow for nonlocal interactions over a finite horizon. The proposed model, which we refer to as the nonlocal pair interaction model, inherits at the continuum level the unwinding feature of finite difference schemes for local hyperbolic conservation laws, so that the maximum principle and certain monotonicity properties hold and, consequently, the entropy inequalities are naturally satisfied. We establish a global-in-time well-posedness theory for these models which covers a broad class of initial data. Moreover, in the limit when the horizon parameter approaches zero, we are able to prove that our nonlocal model reduces to the conventional class of local hyperbolic conservation laws. Furthermore, we propose a numerical discretization method adapted to our nonlocal model, which relies on a monotone numerical flux and a uniform mesh, and we establish that these numerical solutions converge to a solution, providing as by-product both the existence theory for the nonlocal model and the convergence property relating the nonlocal regime and the asymptotic local regime.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The global nonlinear stability of Minkowski space for the Einstein equations in presence of a massive field

We provide a significant extension of the Hyperboloidal Foliation Method introduced by the authors in 2014 in order to establish global existence results for systems of quasilinear wave equations posed on a curved space, when wave equations and Klein-Gordon equations are coupled. This method is based on a (3+1) foliation (of the interior of a future light cone in Minkowski spacetime) by spacelike and asymptotically hyperboloidal hypersurfaces. In the new formulation of the method, we succeed to cover wave-Klein-Gordon systems containing "strong interaction" terms at the level of the metric, and then generalize our method in order to establish a new existence theory for the Einstein equations of general relativity. Following pioneering work by Lindblad and Rodnianski on the Einstein equations in wave coordinates, we establish the nonlinear stability of Minkowski spacetime for self-gravitating massive scalar fields.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The global nonlinear stability of Minkowski spacetime for self-gravitating massive fields. A brief overview

This is a short review of a series of papers which, in collaboration with Yue Ma, establish several novel existence results for systems of coupled wave-Klein-Gordon equation. Our method, the Hyperbolic Hyperboloidal Method, has allowed us to address the global evolution problem for the Einstein equations of general relativity and investigate the global geometry of matter spacetimes that are initially close to Minkowski spacetime.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Weakly regular fluid flows with bounded variation on the domain of outer communication of a Schwarzschild black hole spacetime. A numerical study

We study the dynamical behavior of compressible fluids evolving on the outer domain of communication of a Schwarzschild background. To this end, we design several numerical methods which take the Schwarzschild geometry into account and we treat, both, the relativistic Burgers equation and the relativistic Euler system under the assumption that the flow is spherically symmetric. All the schemes we construct are proven to be well-balanced and therefore to preserve the family of steady state solutions for both models. They enable us to study the nonlinear stability of fluid equilibria, and in particular to investigate the behavior of the fluid near the blackhole horizon. We state and numerically demonstrate several conjectures about the late-time behavior of perturbations of steady solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A central-upwind geometry-preserving method for hyperbolic conservation laws on the sphere

We introduce a second-order, central-upwind finite volume method for the discretization of nonlinear hyperbolic conservation laws posed on the two-dimensional sphere. The semi-discrete version of the proposed method is based on a technique of local propagation speeds and it is free of any Riemann solver. The main advantages of our scheme are the high resolution of discontinuous solutions, its low numerical dissipation, and its simplicity for the implementation. The proposed scheme does not use any splitting approach, which is applied in some cases to upwind schemes in order to simplify the resolution of Riemann problems. The semi-discrete form of the scheme is strongly linked to the analytical properties of the nonlinear conservation law and to the geometry of the sphere. The curved geometry is treated here in an analytical way so that the semi-discrete form of the proposed scheme is consistent with a geometric compatibility property. Furthermore, the time evolution is carried out by using a total-variation-diminishing Runge-Kutta method. A rich family of (discontinuous) stationary solutions is available for the problem under consideration when the flux is nonlinear and foliated (as identified by the author in an earlier work). We present here a series of numerical examples, obtained by considering non-trivial steady state solutions and this leads us to a good validation of the accuracy and efficiency of the proposed central-upwind finite volume method. Our numerical tests confirm the stability of the proposed scheme and clearly show its ability to capture accurately discontinuous steady state solutions to nonlinear hyperbolic conservation laws posed on the sphere.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Lagrangian formulation and a priori estimates for relativistic fluid flows with vacuum

We study the evolution of a compressible fluid surrounded by vacuum and introduce a new symmetrization in Lagrangian coordinates that allows us to encompass both relativistic and non-relativistic fluid flows. The problem under consideration is a free boundary problem of central interest in compressible fluid dynamics and, from the mathematical standpoint, the main challenge to be overcome lies in the loss of regularity in the fluid variables near the free boundary. Based on our Lagrangian formulation, we establish the necessary a priori estimates in weighted Sobolev spaces which are adapted to this loss of regularity.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Structure-preserving shock-capturing methods: late-time asymptotics, curved geometry, small-scale dissipation, and nonconservative products

We present some recent developments on shock capturing methods for nonlinear hyperbolic systems of balance laws, whose prototype is the Euler system of compressible fluid flows, and especially discuss {structure-preserving} techniques. The problems under consideration arise with complex fluids in realistic applications when friction terms, geometrical terms, viscosity and capillarity effects, etc., need to be taken into account in order to achieve a proper description of the physical phenomena. For these problems, it is necessary to design numerical methods that are not only consistent with the given partial differential equations, but remain accurate and robust in certain {asymptotic regimes} of physical interest. That is, certain structural properties of these hyperbolic problems (conservation or balance law, equilibrium state, monotonicity properties, etc.) are essential in many applications, and one seeks that the numerical solutions preserve these properties, which is often a very challenging task.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The global nonlinear stability of Minkowski space for self-gravitating massive fields. The wave-Klein-Gordon model

The Hyperboloidal Foliation Method (introduced by the authors in 2014) is extended here and applied to the Einstein equations of general relativity. Specifically, we establish the nonlinear stability of Minkowski spacetime for self-gravitating massive scalar fields, while existing methods only apply to massless scalar fields. First of all, by analyzing the structure of the Einstein equations in wave coordinates, we exhibit a nonlinear wave-Klein-Gordon model defined on a curved background, which is the focus of the present paper. For this model, we prove here the existence of global-in-time solutions to the Cauchy problem, when the initial data have sufficiently small Sobolev norms. A major difficulty comes from the fact that the class of conformal Killing fields of Minkowski space is significantly reduced in presence of a massive scalar field, since the scaling vector field is not conformal Killing for the Klein-Gordon operator. Our method relies on the foliation (of the interior of the light cone) of Minkowski spacetime by hyperboloidal hypersurfaces and uses Lorentz-invariant energy norms. We introduce a frame of vector fields adapted to the hyperboloidal foliation and we establish several key properties: Sobolev and Hardy-type inequalities on hyperboloids, as well as sup-norm estimates which correspond to the sharp time decay for the wave and the Klein-Gordon equations. These estimates allow us to control interaction terms associated with the curved geometry and the massive field, by distinguishing between two levels of regularity and energy growth and by a successive use of our key estimates in order to close a bootstrap argument.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The global nonlinear stability of self-gravitating irrotational Chaplygin fluids in a FRW geometry

We analyze the global nonlinear stability of FRW (Friedmann-Robertson-Walker) spacetimes in presence of an irrotational perfect fluid. We assume that the fluid is governed by the so-called (generalized) Chaplygin equation of state relating the pressure to the mass-energy density. We express the Einstein equations in wave gauge as a systems of coupled nonlinear wave equations and by performing a suitable conformal transformation, we are able to analyze the global behavior of solutions in future timelike directions. We establish that the (3+1)-spacetime metric and the mass density and velocity vector describing the evolution of the fluid remain globally close to a reference FRW solution, under small initial data perturbations. Our analysis provides also the precise asymptotic behavior of the perturbed solutions in the future directions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The relativistic Burgers equation on a FLRW background and its finite volume approximation

A relativistic generalization of the inviscid Burgers equation was proposed by LeFloch, Makhlof, and Okutmustur and then investigated on a Schwarzschild background. Here, we extend their analysis to a Friedmann-Lemaitre-Robertson-Walker (FLRW) background. This problem is more challenging due to the existence of non-trivial spatially homogeneous solutions. First, we derive the relativistic Burgers model of interest and determine its spatially homogeneous solutions. Second, we design a numerical scheme based on the finite volume methodology, which is well-balanced in the sense that spatially homogeneous solutions are preserved at the discrete level of approximation. Numerical experiments demonstrate the efficiency of the proposed method for weak solutions containing shock waves.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Weakly regular fluid flows with bounded variation on a Schwarzschild background

We study the global dynamics of isothermal fluids evolving in the domain of outer communication of a Schwarzschild black hole. We first formulate the initial value problem within a class of weak solutions with bounded variation (BV), possibly containing shock waves. We then introduce a version of the random choice method and establish a global-in-time existence theory for the initial value problem within the proposed class of weakly regular fluid flows. The initial data may have arbitrary large bounded variation and can possibly blow up near the horizon of the black hole. Furthermore, we study the class of possibly discontinuous, equilibrium solutions and design a version of the random choice method in which these fluid equilibria are exactly preserved. This leads us to a nonlinear stability property for fluid equilibria under small perturbations with bounded variation. Furthermore, we can also encompass several limiting regimes (stiff matter, non-relativistic flows, extremal black hole) by letting the physical parameters (mass of the black hole, light speed, sound speed) reach extremal values.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Mathematical Validity of the f(R) Theory of Modified Gravity

We establish a well-posedness theory for the f(R) theory of modified gravity, which is a generalization of Einstein's theory of gravitation. The scalar curvature R of the spacetime, which arises in the integrand of the Einstein-Hilbert functional, is replaced by a function f=f(R). The field equations involve up to fourth-order derivatives of the spacetime metric, and the challenge is to understand the structure of these high-order terms. We propose a formulation of the initial value problem in modified gravity when the initial data are prescribed on a hypersurface. In addition to the induced metric and second fundamental form of the initial slice and the initial matter content, an initial data set must also provide the spacetime scalar curvature and its time-derivative. We introduce an augmented conformal formulation, as we call it, in which the spacetime curvature is an independent variable. In the so-called wave gauge, the field equations of modified gravity reduce a coupled system of nonlinear wave-Klein-Gordon equations with defocusing potential, whose unknowns are the conformally-transformed metric and the scalar curvature, as well as the matter fields. We are able to establish the existence of maximal globally hyperbolic developments when the matter is represented by a scalar field. We analyze the so-called Jordan coupling and work with the Einstein metric, which is conformally equivalent to the physical metric --the conformal factor depending upon the scalar curvature. Our analysis leads us to a rigorous validation of the theory of modified gravity. We derive quantitative estimates, which are uniform in terms of the nonlinearity f(R), and we prove that asymptotically flat spacetimes of modified gravity are `close' to Einstein spacetimes, when the function f(R) in the action functional of modified gravity is `close' to the Einstein-Hilbert integrand R.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the area of the symmetry orbits in weakly regular Einstein-Euler spacetimes with Gowdy symmetry

This paper establishes novel bounds for Gowdy-symmetric Einstein-Euler spacetimes and completes the analysis, initiated by LeFloch and Rendall, of the global areal foliation for these spacetimes. We thus consider the initial value problem for the Einstein-Euler equations under the assumption of Gowdy symmetry. We establish that, for the maximal Cauchy development of future contracting initial data, the area of the group orbits approaches zero toward the future. This property holds as one approaches the future boundary of the spacetime, provided a geometry invariant associated with the Gowdy symmetry property is initially non-vanishing. Our condition is sharp within the class of spatially homogeneous spacetimes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Revisiting the method of characteristics via a convex hull algorithm

We revisit the method of characteristics for shock wave solutions to nonlinear hyperbolic problems and we describe a novel numerical algorithm - the convex hull algorithm (CHA) - in order to compute, both, entropy dissipative solutions (satisfying all relevant entropy inequalities) and entropy conservative (or multivalued) solutions to nonlinear hyperbolic conservation laws. Our method also applies to Hamilton-Jacobi equations and other problems endowed with a method of characteristics. From the multivalued solutions determined by the method of characteristic, our algorithm "extracts" the entropy dissipative solutions, even after the formation of shocks. It applies to, both, convex or non-convex flux/Hamiltonians. We demonstrate the relevance of the proposed approach with a variety of numerical tests including a problem from fluid dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The equations of elastostatics in a Riemannian manifold

To begin with, we identify the equations of elastostatics in a Riemannian manifold, which generalize those of classical elasticity in the three-dimensional Euclidean space. Our approach relies on the principle of least energy, which asserts that the deformation of the elastic body arising in response to given loads minimizes over a specific set of admissible deformations the total energy of the elastic body, defined as the difference between the strain energy and the potential of the loads. Assuming that the strain energy is a function of the metric tensor field induced by the deformation, we first derive the principle of virtual work and the associated nonlinear boundary value problem of nonlinear elasticity from the expression of the total energy of the elastic body. We then show that this boundary value problem possesses a solution if the loads are sufficiently small (in a sense we specify).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The hyperboloidal foliation method

The Hyperboloidal Foliation Method presented in this monograph is based on a (3+1)-foliation of Minkowski spacetime by hyperboloidal hypersurfaces. It allows us to establish global-in-time existence results for systems of nonlinear wave equations posed on a curved spacetime and to derive uniform energy bounds and optimal rates of decay in time. We are also able to encompass the wave equation and the Klein-Gordon equation in a unified framework and to establish a well-posedness theory for nonlinear wave-Klein-Gordon systems and a large class of nonlinear interactions. The hyperboloidal foliation of Minkowski spacetime we rely upon in this book has the advantage of being geometric in nature and, especially, invariant under Lorentz transformations. As stated, our theory applies to many systems arising in mathematical physics and involving a massive scalar field, such as the Dirac-Klein-Gordon system. As it provides uniform energy bounds and optimal rates of decay in time, our method appears to be very robust and should extend to even more general systems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The positive mass theorem for manifolds with distributional curvature

We formulate and prove a positive mass theorem for n-dimensional spin manifolds whose metrics have only the Sobolev regularity $C^0 \cap W^{1,n}$. At this level of regularity, the curvature of the metric is defined in the distributional sense only, and we propose here a (generalized) notion of ADM mass for such a metric. Our main theorem establishes that if the manifold is asymptotically flat and has non-negative scalar curvature distribution, then its (generalized) ADM mass is well-defined and non-negative, and vanishes only if the manifold is isometric to Euclidian space. Prior applications of Witten's spinor method by Lee and Parker and by Bartnik required the much stronger regularity $W^{2,2}$. Our proof is a generalization of Witten's arguments, in which we must treat the Dirac operator and its associated Lichnerowicz-Weitzenbock identity in the distributional sense and cope with certain averages of first-order derivatives of the metric over annuli that approach infinity. Finally, we observe that our arguments are not specific to scalar curvature and also allow us to establish a universal positive mass theorem.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Weakly regular T2 symmetric spacetimes. The future causal geometry of Gowdy spaces

We investigate the future asymptotic behavior of Gowdy spacetimes on T3, when the metric satisfies weak regularity conditions, so that the metric coefficients (in suitable coordinates) are only in the Sobolev space H1 or have even weaker regularity. The authors recently introduced this class of spacetimes in the broader context of T2 symmetric spacetimes and established the existence of a global foliation by spacelike hypersurfaces when the time function is chosen to be the area of the surfaces of symmetry. In the present paper, we identify the global causal geometry of these spacetimes and, in particular, establish that weakly regular Gowdy spacetimes are future causally geodesically complete. This result extends a theorem by Ringström for metrics with sufficiently high regularity. We emphasize that our proof of the energy decay is based on an energy functional inspired by the Gowdy-to-Ernst transformation. In order to establish the geodesic completeness property, we prove a higher regularity property concerning the metric coefficients along timelike curves and we provide a novel analysis of the geodesic equation for Gowdy spacetimes, which does not require high-order regularity estimates. Even when sufficient regularity is assumed, our proof provides an alternative and shorter proof of the energy decay and of the geodesic completeness property for Gowdy spacetimes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Numerical methods with controlled dissipation for small-scale dependent shocks

We provide a `user guide' to the literature of the past twenty years concerning the modeling and approximation of discontinuous solutions to nonlinear hyperbolic systems that admit small-scale dependent shock waves. We cover several classes of problems and solutions: nonclassical undercompressive shocks, hyperbolic systems in nonconservative form, boundary layer problems. We review the relevant models arising in continuum physics and describe the numerical methods that have been proposed to capture small-scale dependent solutions. In agreement with the general well-posedness theory, small-scale dependent solutions are characterized by a kinetic relation, a family of paths, or an admissible boundary set. We provide a review of numerical methods (front tracking schemes, finite difference schemes, finite volume schemes), which, at the discrete level, reproduce the effect of the physically-meaningful dissipation mechanisms of interest in the applications. An essential role is played by the equivalent equation associated with discrete schemes, which is found to be relevant even for solutions containing shock waves.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

High-order asymptotic-preserving methods for fully nonlinear relaxation problems

We study solutions to nonlinear hyperbolic systems with fully nonlinear relaxation terms in the limit of, both, infinitely stiff relaxation and arbitrary late time. In this limit, the dynamics is governed by effective systems of parabolic type, with possibly degenerate and/or fully nonlinear diffusion terms. For this class of problems, we develop here an implicit-explicit method based on Runge-Kutta discretization in time, and we use this method in order to investigate several examples of interest in compressible fluid dynamics. Importantly, we impose here a realistic stability condition on the time-step and we demonstrate that solutions in the hyperbolic-to-parabolic regime can be computed numerically with high robustness and accuracy, even in presence of fully nonlinear relaxation terms.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Plane-symmetric spacetimes with positive cosmological constant. The case of stiff fluids

We consider plane-symmetric spacetimes satisfying Einstein's field equations with positive cosmological constant, when the matter is a fluid whose pressure is equal to its mass-energy density (i.e. a so-called stiff fluid). We study the initial-value problem for the associated Einstein equations and establish a global existence result. The late-time asymptotics of solutions is also rigorously derived, and we conclude that the spacetime approaches the de Sitter spacetime while the matter disperses asymptotically. A technical difficulty dealt with here lies in the fact that solutions may contain vacuum states as well as velocities approaching the speed of light, both possibilities leading to singular behavior in the evolution equations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Relativistic Burgers equations on curved spacetimes. Derivation and finite volume approximation

Within the class of nonlinear hyperbolic balance laws posed on a curved spacetime (endowed with a volume form), we identify a hyperbolic balance law that enjoys the same Lorentz invariance property as the one satisfied by the Euler equations of relativistic compressible fluids. This model is unique up to normalization and converges to the standard inviscid Burgers equation in the limit of infinite light speed. Furthermore, from the Euler system of relativistic compressible flows on a curved background, we derive, both, the standard inviscid Burgers equation and our relativistic generalizations. The proposed models are referred to as relativistic Burgers equations on curved spacetimes and provide us with simple models on which numerical methods can be developed and analyzed. Next, we introduce a finite volume scheme for the approximation of discontinuous solutions to these relativistic Burgers equations. Our scheme is formulated geometrically and is consistent with the natural divergence form of the balance laws under consideration. It applies to weak solutions containing shock waves and, most importantly, is well-balanced in the sense that it preserves steady solutions. Numerical experiments are presented which demonstrate the convergence of the proposed finite volume scheme and its relevance for computing entropy solutions on a curved background.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The finite energy method for compressible fluids. The Navier-Stokes-Korteweg model

This is the first of a series of papers devoted to the initial value problem for the Euler system of compressible fluids and augmented versions containing higher-order terms. We encompass solutions that have finite total energy and enjoy a certain symmetry (for instance, plane symmetry); these solutions may have unbounded amplitude and contain cavitation regions in which the mass density vanishes. In the present paper, we are interested in dispersive shock waves and analyze the zero viscosity-capillarity limit associated with the Navier-Stokes-Korteweg system. Specifically, we establish the existence of finite energy solutions as well as their convergence toward entropy solutions to the Euler system. We encompass a broad class of nonlinear Navier-Stokes-Korteweg constitutive laws, which is determined by two main conditions relating the viscosity and capillarity coefficients, that is, on one hand the strong coercivity condition (as we call it) which provides a favorable sign for the integrated dissipation associated with an effective energy, and on the other hand the tame capillarity condition (as we call it), which restricts pointwise the strength of the capillarity relatively to the viscosity. Rather mild conditions on the growth of the constitutive functions are aso imposed, which are required in order to define finite energy weak solutions to the Navier-Stokes-Korteweg system, even in the presence of cavitation. Our method of proof relies on fine algebraic properties of the Euler system and combines together energy and effective energy estimates, dissipation and effective dissipation estimates, a nonlinear Sobolev inequality, high-integrability properties for the mass density and for the velocity, and compactness properties based on entropies.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The geometric finite volume method for compressible fluid flows on Schwarzschild spacetime

We consider the relativistic Euler equations governing spherically symmetric, perfect fluid flows on the outer domain of communication of Schwarzschild spacetime, and we introduce a version of the finite volume method which is formulated geometrically (without choosing coordinates a priori) and is well--balanced, in the sense that it preserves steady solutions to the Euler equations on the curved geometry under consideration. In order to formulate our method, we first derive a closed formula describing all steady and spherically symmetric solutions to the Euler equations posed on Schwarzschild spacetime. Second, we describe a finite volume method which is formulated geometrically from the family of steady solutions to the Euler system. Our scheme is second--order accurate and, as required, preserves the family of steady solutions at the discrete level. Numerical experiments are presented which demonstrate the efficiency and robustness of the proposed method even for solutions containing shock waves and nonlinear interacting wave patterns. As an application, we investigate the late--time asymptotics of perturbed steady solutions and demonstrate its convergence for late time toward another steady solution, taking the overall effect of the perturbation into account.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The linear stability of shock waves for the nonlinear Schrödinger-Inviscid Burgers system

We investigate the coupling between the nonlinear Schrödinger equation and the inviscid Burgers equation, a system which models interactions between short and long waves, for instance in fluids. Well-posedness for the associated Cauchy problem remains a difficult open problem, and we tackle it here via a linearization technique. Namely, we establish a linearized stability theorem for the Schrödinger--Burgers system, when the reference solution is an entropy--satisfying shock wave to Burgers equation. Our proof is based on suitable energy estimates and on properties of hyperbolic equations with discontinuous coefficients. Numerical experiments support and expand our theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A framework for late-time/stiff relaxation asymptotics

We consider solutions to nonlinear hyperbolic systems of balance laws with stiff relaxation and formally derive a parabolic-type effective system describing the late-time asymptotics of these solutions. We show that many examples from continuous physics fall into our framework, including the Euler equations with (possibly nonlinear) friction. We then turn our attention to the discretization of these stiff problems and introduce a new finite volume scheme which preserves the late-time asymptotic regime. Importantly, our scheme requires only the classical CFL (Courant, Friedrichs, Lewy) condition associated with the hyperbolic system under consideration, rather than the more restrictive, parabolic-type stability condition.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A geometric approach to error estimates for conservation laws posed on a spacetime

We consider a hyperbolic conservation law posed on an (N+1)-dimensional spacetime, whose flux is a field of differential forms of degree N. Generalizing the classical Kuznetsov's method, we derive an L1 error estimate which applies to a large class of approximate solutions. In particular, we apply our main theorem and deal with two entropy solutions associated with distinct flux fields, as well as with an entropy solution and an approximate solution. Our framework encompasses, for instance, equations posed on a globally hyperbolic Lorentzian manifold.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A global foliation of Einstein-Euler spacetimes with Gowdy-symmetry on T3

We investigate the initial value problem for the Einstein-Euler equations of general relativity under the assumption of Gowdy symmetry on T3, and we construct matter spacetimes with low regularity. These spacetimes admit, both, impulsive gravitational waves in the metric (for instance, Dirac mass curvature singularities propagating at light speed) and shock waves in the fluid (i.e., discontinuities propagating at about the sound speed). Given an initial data set, we establish the existence of a future development and we provide a global foliation in terms of a globally and geometrically defined time-function, closely related to the area of the orbits of the symmetry group. The main difficulty lies in the low regularity assumed on the initial data set which requires a distributional formulation of the Einstein-Euler equations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Godunov-type method for the shallow water equations with discontinuous topography in the resonant regime

We investigate the Riemann problem for the shallow water equations with variable and (possibly) discontinuous topography and provide a complete description of the properties of its solutions: existence; uniqueness in the non-resonant regime; multiple solutions in the resonant regime. This analysis leads us to a numerical algorithm that provides one with a Riemann solver. Next, we introduce a Godunov-type scheme based on this Riemann solver, which is well-balanced and of quasi-conservative form. Finally, we present numerical experiments which demonstrate the convergence of the proposed scheme even in the resonance regime, except in the limiting situation when Riemann data precisely belong to the resonance hypersurface.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Late-time/stiff relaxation asymptotic-preserving approximations of hyperbolic equations

We investigate the late-time asymptotic behavior of solutions to nonlinear hyperbolic systems of conservation laws containing stiff relaxation terms. First, we introduce a Chapman-Enskog-type asymptotic expansion and derive an effective system of equations describing the late-time/stiff relaxation singular limit. The structure of this new system is discussed and the role of a mathematical entropy is emphasized. Second, we propose a new finite volume discretization which, in late-time asymptotics, allows us to recover a discrete version of the same effective asymptotic system. This is achieved provided we suitably discretize the relaxation term in a way that depends on a matrix-valued free-parameter, chosen so that the desired asymptotic behavior is obtained. Our results are illustrated with several models of interest in continuum physics, and numerical experiments demonstrate the relevance of the proposed theory and numerical strategy.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Null injectivity estimate under an upper bound on the curvature

We establish a uniform estimate for the injectivity radius of the past null cone of a point in a general Lorentzian manifold foliated by spacelike hypersurfaces and satisfying an upper curvature bound. Precisely, our main assumptions are, on one hand, upper bounds on the null curvature of the spacetime and the lapse function of the foliation, and sup-norm bounds on the deformation tensors of the foliation. Our proof is inspired by techniques from Riemannian geometry, and it should be noted that we impose no restriction on the size of the curvature or deformation tensors, and allow for metrics that are "far" from the Minkowski one. The relevance of our estimate is illustrated with a class of plane-symmetric spacetimes which satisfy our assumptions but admit no uniform lower bound on the curvature not even in the L2 norm. The conditions we put forward, therefore, lead to a uniform control of the spacetime geometry and should be useful in the context of general relativity.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Second-order hyperbolic Fuchsian systems and applications

We introduce a new class of singular partial differential equations, referred to as the second-order hyperbolic Fuchsian systems, and we investigate the associated initial value problem when data are imposed on the singularity. First, we establish a general existence theory of solutions with asymptotic behavior prescribed on the singularity, which relies on a new approximation scheme, suitable also for numerical purposes. Second, this theory is applied to the (vacuum) Einstein equations for Gowdy spacetimes, and allows us to recover, by more direct arguments, well-posedness results established earlier by Rendall and collaborators. Another main contribution in this paper is the proposed approximation scheme, which we refer to as the Fuchsian numerical algorithm and is shown to provide highly accurate numerical approximations to the singular initial value problem. For the class of Gowdy spacetimes, the numerical experiments presented here show the interest and efficiency of the proposed method and demonstrate the existence of a class of Gowdy spacetimes containing a smooth, incomplete, and non-compact Cauchy horizon.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Second-order hyperbolic Fuchsian systems. Asymptotic behavior of geodesics in Gowdy spacetimes

Recent work by the authors led to the development of a mathematical theory dealing with `second--order hyperbolic Fuchsian systems', as we call them. In the present paper, we adopt a physical standpoint and discuss the implications of this theory which provides one with a new tool to tackle the Einstein equations of general relativity (under certain symmetry assumptions). Specifically, we formulate the `Fuchsian singular initial value problem' and apply our general analysis to the broad class of vacuum Gowdy spacetimes with spatial toroidal topology. Our main focus is on providing a detailed description of the asymptotic geometry near the initial singularity of these inhomogeneous cosmological spacetimes and, especially, analyzing the asymptotic behavior of timelike geodesics ---which represent the trajectories of freely falling observers --- and null geodesics. In particular, we numerically construct Gowdy spacetimes which contain a black hole--like region together with a flat Minkowski--like region. By using the Fuchsian technique, we investigate the effect of the gravitational interaction between these two regions and we study the unexpected behavior of geodesic trajectories within the intermediate part of the spacetime limited by these two regions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Second-order hyperbolic Fuchsian systems. General theory

We introduce a class of singular partial differential equations, the second-order hyperbolic Fuchsian systems, and we investigate the associated initial value problem when data are imposed on the singularity. First of all, we analyze a class of equations in which hyperbolicity is not assumed and we construct asymptotic solutions of arbitrary order. Second, for the proposed class of second-order hyperbolic Fuchsian systems, we establish the existence of solutions with prescribed asymptotic behavior on the singularity. Our proof is based on a new scheme which is also suitable to design numerical approximations. Furthermore, as shown in a follow-up paper, the second-order Fuchsian framework is appropriate to handle Einstein's field equations for Gowdy symmetric spacetimes and allows us to recover (and slightly generalize) earlier results by Rendall and collaborators, while providing a direct approach leading to accurate numerical solutions. The proposed framework is also robust enough to encompass matter models arising in general relativity.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Second-order hyperbolic Fuchsian systems. Gowdy spacetimes and the Fuchsian numerical algorithm

This is the second part of a series devoted to the singular initial value problem for second-order hyperbolic Fuchsian systems. In the first part, we defined and investigated this general class of systems, and we established a well-posedness theory in weighted Sobolev spaces. This theory is applied here to the vacuum Einstein equations for Gowdy spacetimes admitting, by definition, two Killing fields satisfying certain geometric conditions. We recover, by more direct and simpler arguments, the well-posedness results established earlier by Rendall and collaborators. In addition, in this paper we introduce a natural approximation scheme, which we refer to as the Fuchsian numerical algorithm and is directly motivated by our general theory. This algorithm provides highly accurate, numerical approximations of the solution to the singular initial value problem. In particular, for the class of Gowdy spacetimes under consideration, various numerical experiments are presented which show the interest and efficiency of the proposed method. Finally, as an application, we numerically construct Gowdy spacetimes containing a smooth, incomplete, non-compact Cauchy horizon.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The characteristic initial value problem for plane symmetric spacetimes with weak regularity

We investigate the existence and the global causal structure of plane symmetric spacetimes with weak regularity when the matter consists of an irrotational perfect fluid with pressure equal to its mass-energy density. Our theory encompasses the class of weakly regular spacetimes whose metric coefficients have square-integrable first-order derivatives and whose curvature must be understood in the sense of distributions. We formulate the characteristic initial value problem with data posed on two null hypersurfaces intersecting along a two-plane. Relying on Newman-Penrose's formalism and expressing our weak regularity conditions in terms of the Newman-Penrose scalars, we arrive at a fully geometrical formulation in which, along each initial hypersurface, two scalar fields describing the incoming radiation must be prescribed. To analyze the future boundary of such a spacetime and identify its global causal structure, we introduce a gauge that reduces the Einstein equations to a coupled system of wave equations and ordinary differential equations for well-chosen unknowns. We prove that, within the weak regularity class under consideration and for generic initial data, a true spacetime singularity forms in finite proper time. Our formulation is robust enough so that propagating discontinuities in the curvature or in the matter variables do not prevent us from constructing a spacetime whose curvature generically blows-up on the future boundary. Earlier work on the problem studied here was restricted to sufficiently regular and vacuum spacetimes.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Weakly regular T2 symmetric spacetimes. The global geometry of future developments

Under weak regularity assumptions, only, we develop a fully geometric theory of vacuum Einstein spacetimes with T2 symmetry, establish the global well-posedness of the initial value problem for Einstein's field equations, and investigate the global causal structure of the constructed spacetimes. Our weak regularity assumptions are the minimal ones allowing to give a meaning to the Einstein equations under the assumed symmetry and to solve the initial value problem. First of all, we introduce a frame adapted to the symmetry in which each Christoffel symbol can be checked to belong to some Lp space. We identify certain cancellation properties taking place in the expression of the Riemann and Ricci curvatures, and this leads us to a reformulation of the initial value problem for the Einstein field equations when the initial data set has weak regularity. Second, we investigate the future development of a weakly regular initial data set. We check that the area R of the orbits of symmetry must grow to infinity in the future timelike directions, and we establish the existence of a global foliation by the level sets of R. Our weak regularity assumptions only require that R is Lipschitz continuous while the metric coefficients describing the initial geometry of the orbits of symmetry are in the Sobolev space H1 and the remaining coefficients have even weaker regularity. We develop here the compactness arguments required to cover the natural level of regularity associated with the energy of the system of partial differential equations determined from Einstein's field equations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Why many theories of shock waves are necessary. Kinetic relations for nonconservative systems

For a class of nonconservative hyperbolic systems of partial differential equations endowed with a strictly convex mathematical entropy, we formulate the initial value problem by supplementing the equations with a kinetic relation prescribing the rate of entropy dissipation across shock waves. Our condition can be regarded as a generalization to nonconservative systems of a similar concept introduced by Abeyaratne, Knowles, and Truskinovsky for subsonic phase transitions and by LeFloch for nonclassical undercompressive shocks to nonlinear hyperbolic systems. The proposed kinetic relation for nonconservative systems turns out to be equivalent, for the class of systems under consideration at least, to Dal Maso, LeFloch, and Murat's definition based on a prescribed family of Lipschitz continuous paths. In agreement with previous theories, the kinetic relation should be derived from a phase plane analysis of traveling solutions associated with an augmented version of the nonconservative system. We illustrate with several examples that nonconservative systems arising in the applications fit in our framework, and for a typical model of turbulent fluid dynamics, we provide a detailed analysis of the existence and properties of traveling waves which yields the corresponding kinetic function.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Existence of naked singularities in Brans-Dicke theory of gravitation. An analytical and numerical study

Within the framework of the scalar-tensor models of gravitation and by relying on analytical and numerical techniques, we establish the existence of a class of spherically symmetric spacetimes containing a naked singularity. Our result relies on and extends a work by Christodoulou on the existence of naked singularities for the Einstein-scalar field equations. We establish that a key parameter in Christodoulou's construction couples to the Brans-Dicke field and becomes a dynamical variable, which enlarges and modifies the phase space of solutions. We recover analytically many properties first identified by Christodoulou, in particular the loss of regularity (especially at the center), and then investigate numerically the properties of these spacetimes.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Global geometry of T2 symmetric spacetimes with weak regularity

We define the class of weakly regular spacetimes with T2 symmetry, and investigate their global geometry structure. We formulate the initial value problem for the Einstein vacuum equations with weak regularity, and establish the existence of a global foliation by the level sets of the area R of the orbits of symmetry, so that each leaf can be regarded as an initial hypersurface. Except for the flat Kasner spacetimes which are known explicitly, R takes all positive values. Our weak regularity assumptions only require that the gradient of R is continuous while the metric coefficients belong to the Sobolev space H1 (or have even less regularity).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Hyperbolic conservation laws on spacetimes

We present a generalization of Kruzkov's theory to manifolds. Nonlinear hyperbolic conservation laws are posed on a differential (n+1)-manifold, called a spacetime, and the flux field is defined as a field of n-forms depending on a parameter. The entropy inequalities take a particularly simple form as the exterior derivative of a family of n-form fields. Under a global hyperbolicity condition on the spacetime, which allows arbitrary topology for the spacelike hypersurfaces of the foliation, we establish the existence and uniqueness of an entropy solution to the initial value problem, and we derive a geometric version of the standard L1 semi-group property. We also discuss an alternative framework in which the flux field consists of a parametrized family of vector fields.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Kinetic relations for undercompressive shock waves. Physical, mathematical, and numerical issues

Kinetic relations are required in order to characterize nonclassical undercompressive shock waves and formulate a well-posed initial value problem for nonlinear hyperbolic systems of conservation laws. Such nonclassical waves arise in weak solutions of a large variety of physical models: phase transitions, thin liquid films, magnetohydrodynamics, Camassa-Holm model, martensite-austenite materials, semi-conductors, combustion theory, etc. This review presents the research done in the last fifteen years which led the development of the theory of kinetic relations for undercompressive shocks and has now covered many physical, mathematical, and numerical issues. The main difficulty overcome here in our analysis of nonclassical entropy solutions comes from their lack of monotonicity with respect to initial data. Undercompressive shocks of hyperbolic conservation laws turn out to exhibit features that are very similar to shocks of nonconservative hyperbolic systems, who were investigated earlier by the author.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Uniqueness of Kottler spacetime and Besse conjecture

We establish a black hole uniqueness theorem for Schwarzschild-de Sitter spacetime, also called Kottler spacetime, which satisfies Einstein's field equations of general relativity with positive cosmological constant. Our result concerns the class of static vacuum spacetimes with compact spacelike slices and regular maximal level set of the lapse function. We provide a characterization of the interior domain of communication of the Kottler spacetime, which surrounds an inner horizon and is surrounded by a cosmological horizon. The proposed proof combines arguments from the theory of partial differential equations and differential geometry, and is centered on a detailed study of a possibly singular foliation. We also apply our technique in the Riemannian setting, and establish the validity of the so-called Besse conjecture.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Hyperbolic conservation laws on the sphere. A geometry-compatible finite volume scheme

We consider entropy solutions to the initial value problem associated with scalar nonlinear hyperbolic conservation laws posed on the two-dimensional sphere. We propose a finite volume scheme which relies on a web-like mesh made of segments of longitude and latitude lines. The structure of the mesh allows for a discrete version of a natural geometric compatibility condition, which arose earlier in the well-posedness theory established by Ben-Artzi and LeFloch. We study here several classes of flux vectors which define the conservation law under consideration. They are based on prescribing a suitable vector field in the Euclidean three-dimensional space and then suitably projecting it on the sphere's tangent plane; even when the flux vector in the ambient space is constant, the corresponding flux vector is a non-trivial vector field on the sphere. In particular, we construct here "equatorial periodic solutions", analogous to one-dimensional periodic solutions to one-dimensional conservation laws, as well as a wide variety of stationary (steady state) solutions. We also construct "confined solutions", which are time-dependent solutions supported in an arbitrarily specified subdomain of the sphere. Finally, representative numerical examples and test-cases are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Local foliations and optimal regularity of Einstein spacetimes

We investigate the local regularity of pointed spacetimes, that is, time-oriented Lorentzian manifolds in which a point and a future-oriented, unit timelike vector (an observer) are selected. Our main result covers the class of Einstein vacuum spacetimes. Under curvature and injectivity bounds only, we establish the existence of a local coordinate chart defined in a ball with definite size in which the metric coefficients have optimal regularity. The proof is based on quantitative estimates, on one hand, for a constant mean curvature (CMC) foliation by spacelike hypersurfaces defined locally near the observer and, on the other hand, for the metric in local coordinates that are spatially harmonic in each CMC slice. The results and techniques in this paper should be useful in the context of general relativity for investigating the long-time behavior of solutions to the Einstein equations.

0 citations0 reviewsNo code linkedOpen access