Virginie Ehrlacher

Amina Benaceur, Alexandre Ern, Virginie Ehrlacher et al.

We investigate new developments of the combined Reduced-Basis and Empirical Interpolation Methods (RB-EIM) for parametrized nonlinear parabolic problems. In many situations, the cost of the EIM in the offline stage turns out to be prohibitive since a significant number of nonlinear time-dependent problems need to be solved using the high-fidelity (or full-order) model. In the present work, we develop a new methodology, the Progressive RB-EIM (PREIM) method for nonlinear parabolic problems.The purpose is to reduce the offline cost while maintaining the accuracy of the RB approximation in the online stage. The key idea is a progressive enrichment of both the EIM approximation and the RB space, in contrast to the standard approach where the EIM approximation and the RB space are built separately. PREIM uses high-fidelity computations whenever available and RB computationsotherwise. Another key feature of each PREIM iteration is to select twice the parameter in a greedy fashion, the second selection being made after computing the high-fidelity solution for the firstly selected value of the parameter. Numerical examples are presented on nonlinear heat transfer problems.

Eric Cances, Virginie Ehrlacher, Tony Lelievre

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.

Eric Cancès, Virginie Ehrlacher, Frederic Legoll et al.

This contribution is the numerically oriented companion article of the work [E. Cancès, V. Ehrlacher, F. Legoll, B. Stamm and S. Xiang, arxiv preprint 1807.05131]. We focus here on the numerical resolution of the embedded corrector problem introduced in [E. Cancès, V. Ehrlacher, F. Legoll and B. Stamm, CRAS 2015; E. Cancès, V. Ehrlacher, F. Legoll, B. Stamm and S. Xiang, arxiv preprint 1807.05131] in the context of homogenization of diffusion equations. Our approach consists in considering a corrector-type problem, posed on the whole space, but with a diffusion matrix which is constant outside some bounded domain. In [E. Cancès, V. Ehrlacher, F. Legoll, B. Stamm and S. Xiang, arxiv preprint 1807.05131], we have shown how to define three approximate homogenized diffusion coefficients on the basis of the embedded corrector problems. We have also proved that these approximations all converge to the exact homogenized coefficients when the size of the bounded domain increases. We show here that, under the assumption that the diffusion matrix is piecewise constant, the corrector problem to solve can be recast as an integral equation. In case of spherical inclusions with isotropic materials, we explain how to efficiently discretize this integral equation using spherical harmonics, and how to use the fast multipole method (FMM) to compute the resulting matrix-vector products at a cost which scales only linearly with respect to the number of inclusions. Numerical tests illustrate the performance of our approach in various settings.

Eric Cancès, Virginie Ehrlacher, Frederic Legoll et al.

This article is the first part of a two-fold study, the objective of which is the theoretical analysis and numerical investigation of new approximate corrector problems in the context of stochastic homogenization. We present here three new alternatives for the approximation of the homogenized matrix for diffusion problems with highly-oscillatory coefficients. These different approximations all rely on the use of an embedded corrector problem (that we previously introduced in [Cancès, Ehrlacher, Legoll and Stamm, C. R. Acad. Sci. Paris, 2015]), where a finite-size domain made of the highly oscillatory material is embedded in a homogeneous infinite medium whose diffusion coefficients have to be appropriately determined. The motivation for considering such embedded corrector problems is made clear in the companion article [Cancès, Ehrlacher, Legoll, Stamm and Xiang, in preparation], where a very efficient algorithm is presented for the resolution of such problems for particular heterogeneous materials. In the present article, we prove that the three different approximations we introduce converge to the homogenized matrix of the medium when the size of the embedded domain goes to infinity.

Virginie Ehrlacher, Damiano Lombardi

A numerical method is proposed to solve the full-Eulerian time-dependent Vlasov-Poisson system in high dimension. The algorithm relies on the construction of a tensor decomposition of the solution whose rank is adapted at each time step. This decomposition is obtained through the use of an efficient modified Progressive Generalized Decomposition (PGD) method, whose convergence is proved. We suggest in addition a symplectic time-discretization splitting scheme that preserves the Hamiltonian properties of the system. This scheme is naturally obtained by considering the tensor structure of the approximation. The efficiency of our approach is illustrated through time-dependent 2D-2D numerical examples.

Mi-Song Dupuy, Virginie Ehrlacher, Clément Guillot

We prove in this paper that the solution of the time-dependent Schr{ö}dinger equation can be expressed as the solution of a global space-time quadratic minimization problem that is amenable to Galerkin time-space discretization schemes, using an appropriate least-square formulation. The present analysis can be applied to the electronic many-body time-dependent Schr{ö}dinger equation with an arbitrary number of electrons and interaction potentials with Coulomb singularities. We motivate the interest of the present approach with two goals: first, the design of Galerkin space-time discretization methods; second, the definition of dynamical low-rank approximations following a variational principle different from the classical Dirac-Frenkel principle, and for which it is possible to prove the global-in-time existence of solutions.

Giulia Sambataro, Virginie Ehrlacher

This work investigates model order reduction for time-dependent parametrized variational inequalities, with a focus on discrete contact problems. As a prototypical example, we consider an agent-based crowd model [Maury et al., 2011] in which agent velocities are obtained at each time step from a constrained least-squares problem. Geometric parameter variations induce significant variability in both agent positions and contact forces, leading to a slowly decaying Kolmogorov $n$-width of the solution manifold. We propose a nonlinear approach that combines a linear reduced-order model with a deep-learning-based correction. The method utilizes a greedy index selection (gIS) algorithm for compressing Lagrange multipliers and Proper Orthogonal Decomposition (POD) applied to velocity snapshots. Additionally, we explore hyper-reduction techniques, comparing the Empirical Interpolation Method (EIM) and the Empirical Quadrature (EQ) procedure from both computational complexity and accuracy perspectives. Finally, we demonstrate the applicability of the methodology in a complex scenario involving many agents in a highly congested geometric configuration. This work represents the first attempt to apply model order reduction to a discrete contact problem of the type introduced in [Maury et al., 2011] and paves the way for future advancements in nonlinear MOR specifically for this class of problems.

Elise Bonnet-Weill, Virginie Ehrlacher, Luca Nenna

In this work, we aim at efficiently solving a parametrized family of optimal transport problems by using model order reduction methods. We propose a reduced-order model by adding to the primal (respectively dual) version of the high-fidelity model the additional constraint to live in a non negative sub cone (resp. in subspaces) of small dimension. The reduced-order model then reads as a linear program with a small number of degrees of freedom and constraints. We identify explicit conditions under which this reduced-order model has at least one solution. We propose two a posteriori error estimations that bounds the error between the optimal values of the high-fidelity problem and the reduced-order model. As one of these estimations requires the computation of non linear terms (with respect to the reduction of dimension), we use an Empirical Interpolation Method (EIM) (see e.g. \cite{maday2007general} or \cite{barrault2004empirical}) to numerically efficiently compute this estimation. We apply the whole methodology on a simple 1D example and on a problem of color transfer between images, and compare its performances to Sinkhorn algorithm.

Aurélien Alfonsi, Rafaël Coyaud, Virginie Ehrlacher et al.

Optimal Transport (OT) problems arise in a wide range of applications, from physics to economics. Getting numerical approximate solution of these problems is a challenging issue of practical importance. In this work, we investigate the relaxation of the OT problem when the marginal constraints are replaced by some moment constraints. Using Tchakaloff's theorem, we show that the Moment Constrained Optimal Transport problem (MCOT) is achieved by a finite discrete measure. Interestingly, for multimarginal OT problems, the number of points weighted by this measure scales linearly with the number of marginal laws, which is encouraging to bypass the curse of dimension. This approximation method is also relevant for Martingale OT problems. We show the convergence of the MCOT problem toward the corresponding OT problem. In some fundamental cases, we obtain rates of convergence in $O(1/n)$ or $O(1/n^2)$ where $n$ is the number of moments, which illustrates the role of the moment functions. Last, we present algorithms exploiting the fact that the MCOT is reached by a finite discrete measure and provide numerical examples of approximations.

Thomas Boiveau, Virginie Ehrlacher, Alexandre Ern et al.

We devise a space-time tensor method for the low-rank approximation of linear parabolic evolution equations. The proposed method is a stable Galerkin method, uniformly in the discretization parameters, based on a Minimal Residual formulation of the evolution problem in Hilbert--Bochner spaces. The discrete solution is sought in a trial space composed of tensors of discrete functions in space and in time and is characterized as the unique minimizer of a discrete functional where the dual norm of the residual is evaluated in a space semi-discrete test space. The resulting global space-time linear system is solved iteratively by a greedy algorithm. Numerical results are presented to illustrate the performances of the proposed method on test cases including non-selfadjoint and time-dependent differential operators in space. The results are also compared to those obtained using a fully discrete Petrov--Galerkin setting to evaluate the dual residual norm.

Eric Cances, Virginie Ehrlacher, Frederic Legoll et al.

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.