BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Claude Le Bris appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We consider an homogenization problem for the second order elliptic equation $- Δu^{\varepsilon} + \dfrac{1}{\varepsilon} V(./\varepsilon) u^{\varepsilon} + νu^{\varepsilon} =f$ when the highly oscillatory potential $V$ belongs to a particular class of non-periodic potentials. We show the existence of an adapted corrector and prove the convergence of $u^{\varepsilon}$ to its homogenized limit.
We numerically investigate, and improve upon, a computational approach originally introduced in [Cottereau, IJNME 2013] which aims at evaluating the effective coefficient of a medium modelled by a highly oscillatory coefficient. This computational approach is based on a Arlequin type coupling. It combines the original fine-scale description of the medium (modelled by an oscillatory coefficient) with an effective description (modelled by a constant coefficient) and optimizes upon the coefficient of the effective medium to best fit the response of the actual heterogeneous medium using a purely homogeneous medium. We present here a mathematical formalization of the approach along with various improvements of the algorithms, in order to obtain a procedure as efficient as possible. Representative numerical results demonstrate the added value of our approach in comparison to the original approach.
We overview a series of recent works addressing numerical simulations of partial differential equations in the presence of some elements of randomness. The specific equations manipulated are linear elliptic, and arise in the context of multiscale problems, but the purpose is more general. On a set of prototypical situations, we investigate two critical issues present in many settings: variance reduction techniques to obtain sufficiently accurate results at a limited computational cost when solving PDEs with random coefficients, and finite element techniques that are sufficiently flexible to carry over to geometries with random fluctuations. Some elements of theoretical analysis and numerical analysis are briefly mentioned. Numerical experiments, although simple, provide convincing evidence of the efficiency of the approaches.
The purpose of this work is to investigate the behavior of Multiscale Finite Element type methods for advection-diffusion problems in the advection-dominated regime. We present, study and compare various options to address the issue of the simultaneous presence of both heterogeneity of scales and strong advection. Classical MsFEM methods are compared with adjusted MsFEM methods, stabilized versions of the methods, and a splitting method that treats the multiscale diffusion and the strong advection separately.
We present a twofold contribution to the numerical simulation of Lindblad equations. First, an adaptive numerical approach to approximate Lindblad equations using low-rank dynamics is described: a deterministic low-rank approximation of the density operator is computed, and its rank is adjusted dynamically, using an on-the-fly estimator of the error committed when reducing the dimension. On the other hand, when the intrinsic dimension of the Lindblad equation is too high to allow for such a deterministic approximation, we combine classical ensemble averages of quantum Monte Carlo trajectories and a denoising technique. Specifically, a variance reduction method based upon the consideration of a low-rank dynamics as a control variate is developed. Numerical tests for quantum collapse and revivals show the efficiency of each approach, along with the complementarity of the two approaches.
We overview a series of recent works devoted to variance reduction techniques for numerical stochastic homogenization. Numerical homogenization requires solving a set of problems at the micro scale, the so-called corrector problems. In a random environment, these problems are stochastic and therefore need to be repeatedly solved, for several configurations of the medium considered. An empirical average over all configurations is then performed using the Monte-Carlo approach, so as to approximate the effective coefficients necessary to determine the macroscopic behavior. Variance severely affects the accuracy and the cost of such computations. Variance reduction approaches, borrowed from other contexts of the engineering sciences, can be useful. Some of these variance reduction techniques are presented, studied and tested here.
We adapt and study a variance reduction approach for the homogenization of elliptic equations in divergence form. The approach, borrowed from atomistic simulations and solid-state science [von Pezold et al, Physical Review B 2010; Wei et al, Physical Review B 1990; Zunger et al, Physical Review Letters 1990], consists in selecting random realizations that best satisfy some statistical properties (such as the volume fraction of each phase in a composite material) usually only obtained asymptotically. We study the approach theoretically in some simplified settings (one-dimensional setting, perturbative setting in higher dimensions), and numerically demonstrate its efficiency in more general cases.
We study a one-dimensional equation arising in the multiscale modeling of some non-Newtonian fluids. At a given shear rate, the equation provides the instantaneous mesoscopic response of the fluid, allowing to compute the corresponding stress. In a simple setting, we study the well-posedness of the equation and next the long-time behavior of its solution. In the limit of a response of the fluid much faster than the time variations of the ambient shear rate, we derive some equivalent macroscopic differential equations that relate the shear rate and the stress. Our analytical conclusions are confronted to some numerical experiments. The latter quantitatively confirm our derivations.
We provide in this article a new proof of the uniqueness of the flow solution to ordinary differential equations with $BV$ vector-fields that have divergence in $L^\infty$ (or in $L^1$) and that are nearly incompressible (see the text for the definition of this term). The novelty of the proof lies in the fact it does not use the associated transport equation.
We follow up on our previous work [C. Le Bris, F. Legoll and A. Lozinski, Chinese Annals of Mathematics 2013] where we have studied a multiscale finite element (MsFEM) type method in the vein of the classical Crouzeix-Raviart finite element method that is specifically adapted for highly oscillatory elliptic problems. We adapt the approach to address here a multiscale problem on a perforated domain. An additional ingredient of our approach is the enrichment of the multiscale finite element space using bubble functions. We first establish a theoretical error estimate. We next show that, on the problem we consider, the approach we propose outperforms all dedicated existing variants of MsFEM we are aware of.
A systematic numerical approach to approximate high dimensional Lindblad equations is described. It is based on a deterministic rank m approximation of the density operator, the rank m being the only parameter to adjust. From a known initial value, this rank m approximation gives at each time-step an estimate of the largest m eigen-values with their eigen-vectors of the density operator. A numerical scheme is proposed. Its numerical efficiency in the case of a rank 12 approximation is demonstrated for oscillation revivals of 50 atoms interacting resonantly with a slightly damped coherent quantized field of 200 photons. The approach may be employed for other similar equations. We in particularly show how to adapt such low-rank approximation for Riccati differential equations appearing in Kalman filtering.
We study mathematically a system of partial differential equations arising in the modelling of an aging fluid, a particular class of non Newtonian fluids. We prove well-posedness of the equations in appropriate functional spaces and investigate the longtime behaviour of the solutions.
We study the long-time behaviour of a nonlinear Fokker-Planck equation, which models the evolution of rigid polymers in a given flow, after a closure approximation. The aim of this work is twofold: first, we propose a microscopic derivation of the classical Doi closure, at the level of the kinetic equation ; second, we prove the convergence of the solution to the Fokker-Planck equation to periodic solutions in the long-time limit.
We present in this paper an approach for computing the homogenized behavior of a medium that is a small random perturbation of a periodic reference material. The random perturbation we consider is, in a sense made precise in our work, a rare event at the microscopic level. It however affects the macroscopic properties of the material, and we indeed provide a method to compute the first and second-order corrections. To this end, we formally establish an asymptotic expansion of the macroscopic properties. Our perturbative approach shares common features with a defect-type theory of solid state physics. The computational efficiency of the approach is demonstrated.
This work is a follow-up to our previous work "A numerical approach related to defect-type theories for some weakly random problems in homogenization" (preprint available on this archive). It extends and complements, both theoretically and experimentally, the results presented there. Under consideration is the homogenization of a model of a weakly random heterogeneous material. The material consists of a reference periodic material randomly perturbed by another periodic material, so that its homogenized behavior is close to that of the reference material. We consider laws for the random perturbations more general than in our previous work cited above. We prove the validity of an asymptotic expansion in a certain class of settings. We also extend the formal approach introduced in our former work. Our perturbative approach shares common features with a defect-type theory of solid state physics. The computational efficiency of the approach is demonstrated.
We report here on the recent application of a now classical general reduction technique, the Reduced-Basis approach initiated in [C. Prud'homme, D. Rovas, K. Veroy, Y. Maday, A. T. Patera, and G. Turinici. Reliable real-time solution of parametrized partial differential equations: Reduced-basis output bounds methods. Journal of Fluids Engineering, 124(1):7080, 2002.], to the specific context of differential equations with random coefficients. After an elementary presentation of the approach, we review two contributions of the authors: [S. Boyaval, C. Le Bris, Y. Maday, N.C. Nguyen, and A.T. Patera. A reduced basis approach for variational problems with stochastic parameters: Application to heat conduction with variable Robin co-efficient. Computer Methods in Applied Mechanics and Engineering, 198(4144):3187-3206, 2009.], which presents the application of the RB approach for the discretization of a simple second order elliptic equation supplied with a random boundary condition, and [S. Boyaval and T. Lelièvre, A variance reduction method for parametrized stochastic differential equations using the reduced basis paradigm with T. Lelièvre, Commun. Math. Sci. 8, special Issue "Mathematical Issue on Complex Fluids" P. Zhang ed., to appear, 2010, ARXIV preprint arXiv:0906.3600], which uses a RB type approach to reduce the variance in the Monte-Carlo simulation of a stochastic differential equation. We conclude the review with some general comments and also discuss possible tracks for further research in the direction.
We follow up on our previous works which presented a possible approach for deriving symplectic schemes for a certain class of highly oscillatory Hamiltonian systems. The approach considers the Hamilton-Jacobi form of the equations of motion, formally homogenizes it and infers an appropriate symplectic integrator for the original system. In our previous work, the case of a system exhibiting a single constant fast frequency was considered. The present work successfully extends the approach to systems that have either one varying fast frequency or several constant frequencies. Some related issues are also examined.