Source author record

Virginie Ehrlacher

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

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

16 published item(s)

preprint2023arXiv

A Wasserstein-type metric for generic mixture models, including location-scatter and group invariant measures

In this article, we study Wasserstein-type metrics and corresponding barycenters for mixtures of a chosen subset of probability measures called atoms hereafter. In particular, this works extends what was proposed by Delon and Desolneux [A Wasserstein-Type Distance in the Space of Gaussian Mixture Models. SIAM J. Imaging Sci. 13, 936-970 (2020)] for mixtures of gaussian measures to other mixtures. We first prove in a general setting that for a set of atoms equipped with a metric that defines a geodesic space, the set of mixtures based on this set of atoms is also geodesic space for the defined modified Wasserstein metric. We then focus on two particular cases of sets of atoms: (i) the set of location-scatter atoms and (ii) the set of measures that are invariant with respect to some symmetry group. Both cases are particularly relevant for various applications among which electronic structure calculations. Along the way, we also prove some sparsity and symmetry properties of optimal transport plans between measures that are invariant under some well-chosen symmetries.

preprint2022arXiv

Waserstein model reduction approach for parametrized flow problems in porous media

The aim of this work is to build a reduced-order model for parametrized porous media equations. The main challenge of this type of problems is that the Kolmogorov width of the solution manifold typically decays quite slowly and thus makes usual linear model-order reduction methods inappropriate. In this work, we investigate an adaptation of the methodology proposed in a previous work, based on the use of Wasserstein barycenters, to the case of non-conservative problems. Numerical examples in one-dimensional test cases illustrate the advantages and limitations of this approach and suggest further research directions that we intend to explore in the future.

preprint2021arXiv

Constrained overdamped Langevin dynamics for symmetric multimarginal optimal transportation

The Strictly Correlated Electrons (SCE) limit of the Levy-Lieb functional in Density Functional Theory (DFT) gives rise to a symmetric multi-marginal optimal transport problem with Coulomb cost, where the number of marginal laws is equal to the number of electrons in the system, which can be very large in relevant applications. In this work, we design a numerical method, built upon constrained overdamped Langevin processes to solve Moment Constrained Optimal Transport (MCOT) relaxations (introduced in A. Alfonsi, R. Coyaud, V. Ehrlacher and D. Lombardi, Math. Comp. 90, 2021, 689--737) of symmetric multi-marginal optimal transport problems with Coulomb cost. Some minimizers of such relaxations can be written as discrete measures charging a low number of points belonging to a space whose dimension, in the symmetrical case, scales linearly with the number of marginal laws. We leverage the sparsity of those minimizers in the design of the numerical method and prove that any local minimizer to the resulting problem is actually a \emph{global} one. We illustrate the performance of the proposed method by numerical examples which solves MCOT relaxations of 3D systems with up to 100 electrons.

preprint2021arXiv

Multi-center decomposition of molecular densities: a mathematical perspective

The aim of this paper is to analyze from a mathematical perspective some existing schemes to partition a molecular density into several atomic contributions, with a specific focus on Iterative Stockholder Atom (ISA) methods. We provide a unified mathematical framework to describe the latter family of methods and propose a new scheme, named L-ISA (for linear approximation of ISA). We prove several important mathematical properties of the ISA and L-ISA minimization problems and show that the so-called ISA algorithms can be viewed as alternating minimization schemes, which in turn enables us to obtain new convergence results for these numerical methods. Specific mathematical properties of the ISA decomposition for diatomic systems are also presented. We also review the basis-space oriented Distributed Multipole Analysis method, the mathematical formulation of which is also clarified. Different schemes are numerically compared on different molecules and we discuss the advantages and drawbacks of each approach.

preprint2020arXiv

Adaptive force biasing algorithms: new convergence results and tensor approximations of the bias

A modification of the Adaptive Biasing Force method is introduced, in which the free energy is approximated by a sum of tensor products of one-dimensional functions. This enables to handle a larger number of reaction coordinates than the classical algorithm. We prove the algorithm is well-defined and prove the long-time convergence toward a regularized version of the free energy for an idealized version of the algorithm. Numerical experiments demonstrate that the method is able to capture correlations between reaction coordinates.

preprint2020arXiv

Existence of weak solutions to a cross-diffusion Cahn-Hilliard type system

The aim of this article is to study a Cahn-Hilliard model for a multicomponent mixture with cross-diffusion effects, degenerate mobility and where only one of the species does separate from the others. We define a notion of weak solution adapted to possible degeneracies and our main result is (global in time) existence. In order to overcome the lack of a-priori estimates, our proof uses the formal gradient flow structure of the system and an extension of the boundedness by entropy method which involves a careful analysis of an auxiliary variational problem. This allows to obtain solutions to an approximate, time-discrete system. Letting the time step size go to zero, we recover the desired weak solution where, due to their low regularity, the Cahn-Hilliard terms require a special treatment.

preprint2020arXiv

Finite volumes for the Stefan-Maxwell cross-diffusion system

The aim of this work is to propose a provably convergent finite volume scheme for the so-called Stefan-Maxwell model, which describes the evolution of the composition of a multi-component mixture and reads as a cross-diffusion system. The scheme proposed here relies on a two-point flux approximation, and preserves at the discrete level some fundamental theoretical properties of the continuous models, namely the non-negativity of the solutions, the conservation of mass and the preservation of the volume-filling constraints. In addition, the scheme satisfies a discrete entropy-entropy dissipation relation, very close to the relation which holds at the continuous level. In this article, we present this scheme together with its numerical analysis, and finally illustrate its behaviour with some numerical results.

preprint2015arXiv

Global existence of bounded weak solutions to degenerate cross-diffusion equations in moving domain

The aim of this note is to present preliminary existence results for a system of cross-diffusion equations defined on a domain with moving boundaries, which model the evolution of the concentrations of different chemical species in a solid during a Chemical Vapor Deposition process. The system of equations, when the domain remains fixed over time, can be seen formally as a gradient flow system which can be analyzed using the boundedness by entropy method introduced by Burger and Jüngel. Preliminary existence results are presented in the case of a one-dimensional moving boundary domain.

preprint2015arXiv

Statistical methods for critical scenarios in aeronautics

We present numerical results obtained on the CEMRACS project Predictive SMS proposed by Safety Line. The goal of this work was to elaborate a purely statistical method in order to reconstruct the deceleration profile of a plane during landing under normal operating conditions, from a database containing around $1500$ recordings. The aim of Safety Line is to use this model to detect malfunctions of the braking system of the plane from deviations of the measured deceleration profile of the plane to the one predicted by the model. This yields to a multivariate nonparametric regression problem, which we chose to tackle using a Bayesian approach based on the use of gaussian processes. We also compare this approach with other statistical methods.

preprint2014arXiv

An embedded corrector problem to approximate the homogenized coefficients of an elliptic equation

We consider a diffusion equation with highly oscillatory coefficients that admits a homogenized limit. As an alternative to standard corrector problems, we introduce here an embedded corrector problem, written as a diffusion equation in the whole space in which the diffusion matrix is uniform outside some ball of radius $R$. Using that problem, we next introduce three approximations of the homogenized coefficients. These approximations, which are variants of the standard approximations obtained using truncated (supercell) corrector problems, are shown to converge when $R \to \infty$. We also discuss efficient numerical methods to solve the embedded corrector problem.

preprint2013arXiv

Greedy algorithms for high-dimensional eigenvalue problems

In this article, we present two new greedy algorithms for the computation of the lowest eigenvalue (and an associated eigenvector) of a high-dimensional eigenvalue problem, and prove some convergence results for these algorithms and their orthogonalized versions. The performance of our algorithms is illustrated on numerical test cases (including the computation of the buckling modes of a microstructured plate), and compared with that of another greedy algorithm for eigenvalue problems introduced by Ammar and Chinesta.

preprint2012arXiv

Greedy algorithms for high-dimensional non-symmetric linear problems

In this article, we present a family of numerical approaches to solve high-dimensional linear non-symmetric problems. The principle of these methods is to approximate a function which depends on a large number of variates by a sum of tensor product functions, each term of which is iteratively computed via a greedy algorithm. There exists a good theoretical framework for these methods in the case of (linear and nonlinear) symmetric elliptic problems. However, the convergence results are not valid any more as soon as the problems considered are not symmetric. We present here a review of the main algorithms proposed in the literature to circumvent this difficulty, together with some new approaches. The theoretical convergence results and the practical implementation of these algorithms are discussed. Their behaviors are illustrated through some numerical examples.

preprint2012arXiv

Non-consistent approximations of self-adjoint eigenproblems: Application to the supercell method

In this article, we introduce a general theoretical framework to analyze non-consistent approximations of the discrete eigenmodes of a self-adjoint operator. We focus in particular on the discrete eigenvalues laying in spectral gaps. We first provide a priori error estimates on the eigenvalues and eigenvectors in the absence of spectral pollution. We then show that the supercell method for perturbed periodic Schrödinger operators falls into the scope of our study. We prove that this method is spectral pollution free, and we derive optimal convergence rates for the planewave discretization method, taking numerical integration errors into account. Some numerical illustrations are provided.

preprint2011arXiv

Convergence of a greedy algorithm for high-dimensional convex nonlinear problems

In this article, we present a greedy algorithm based on a tensor product decomposition, whose aim is to compute the global minimum of a strongly convex energy functional. We prove the convergence of our method provided that the gradient of the energy is Lipschitz on bounded sets. The main interest of this method is that it can be used for high-dimensional nonlinear convex problems. We illustrate this method on a prototypical example for uncertainty propagation on the obstacle problem.

preprint2011arXiv

Periodic Schrödinger operators with local defects and spectral pollution

This article deals with the numerical calculation of eigenvalues of perturbed periodic Schrödinger operators located in spectral gaps. Such operators are encountered in the modeling of the electronic structure of crystals with local defects, and of photonic crystals. The usual finite element Galerkin approximation is known to give rise to spectral pollution. In this article, we give a precise description of the corresponding spurious states. We then prove that the supercell model does not produce spectral pollution. Lastly, we extend results by Lewin and Séré on some no-pollution criteria. In particular, we prove that using approximate spectral projectors enables one to eliminate spectral pollution in a given spectral gap of the reference periodic Schödinger operator.

preprint2010arXiv

Local defects are always neutral in the Thomas-Fermi-von Weiszäcker theory of crystals

The aim of this article is to propose a mathematical model describing the electronic structure of crystals with local defects in the framework of the Thomas-Fermi-von Weizsäcker (TFW) theory. The approach follows the same lines as that used in {\it E. Cancès, A. Deleurence and M. Lewin, Commun. Math. Phys., 281 (2008), pp. 129--177} for the reduced Hartree-Fock model, and is based on thermodynamic limit arguments. We prove in particular that it is not possible to model charged defects within the TFW theory of crystals. We finally derive some additional properties of the TFW ground state electronic density of a crystal with a local defect, in the special case when the host crystal is modelled by a homogeneous medium.