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

Frederic Legoll

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

ResearcherUnclaimed source record
math.NAmath.PRcond-mat.stat-mechcond-mat.mtrl-scicond-mat.othermath.APmath.DS

Catalog footprint

What is connected

22works
7topics
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

22 published item(s)

preprint2016arXiv

Examples of computational approaches to accommodate randomness in elliptic PDEs

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Pathwise estimates for an effective dynamics

Starting from the overdamped Langevin dynamics in $\mathbb{R}^n$, $$ dX_t = -\nabla V(X_t) dt + \sqrt{2 β^{-1}} dW_t, $$ we consider a scalar Markov process $ξ_t$ which approximates the dynamics of the first component $X^1_t$. In the previous work [F. Legoll, T. Lelievre, Nonlinearity 2010], the fact that $(ξ_t)_{t \ge 0}$ is a good approximation of $(X^1_t)_{t \ge 0}$ is proven in terms of time marginals, under assumptions quantifying the timescale separation between the first component and the other components of $X_t$. Here, we prove an upper bound on the trajectorial error $\mathbb{E} \left( \sup_{0 \leq t \leq T} \left| X^1_t - ξ_t \right| \right)$, for any $T > 0$, under a similar set of assumptions. We also show that the technique of proof can be used to obtain quantitative averaging results.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A numerical comparison of some Multiscale Finite Element approaches for convection-dominated problems in heterogeneous media

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Coupling a reactive potential with a harmonic approximation for atomistic simulations of material failure

Molecular dynamics (MD) simulations involving reactive potentials can be used to model material failure. The empirical potentials which are used in such simulations are able to adapt to the atomic environment, at the expense of a significantly higher computational cost than non-reactive potentials. However, during a simulation of failure, the reactive ability is needed only in some limited parts of the system, where bonds break or form and the atomic environment changes. Therefore, simpler non-reactive potentials can be used in the remainder of the system, provided that such potentials reproduce correctly the behavior of the reactive potentials in this region, and that seamless coupling is ensured at the interface between the reactive and non-reactive regions. In this article, we propose a methodology to combine a reactive potential with a non-reactive approximation thereof, made of a set of harmonic pair and angle interactions and whose parameters are adjusted to predict the same energy, geometry and Hessian in the ground state of the potential. We present a methodology to construct the non-reactive approximation of the reactive potential, and a way to couple these two potentials. We also propose a criterion for on-the-fly substitution of the reactive potential by its non-reactive approximation during a simulation. We illustrate the correctness of this hybrid technique for the case of MD simulation of failure in two-dimensional graphene originally modeled with REBO potential.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Some variance reduction methods for numerical stochastic homogenization

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Special Quasirandom Structures: a selection approach for stochastic homogenization

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A control variate approach based on a defect-type theory for variance reduction in stochastic homogenization

We consider a variance reduction approach for the stochastic homogenization of divergence form linear elliptic problems. Although the exact homogenized coefficients are deterministic, their practical approximations are random. We introduce a control variate technique to reduce the variance of the computed approximations of the homogenized coefficients. Our approach is based on a surrogate model inspired by a defect-type theory, where a perfect periodic material is perturbed by rare defects. This model has been introduced in [A. Anantharaman and C. Le Bris, CRAS 2010] in the context of weakly random models. In this work, we address the fully random case, and show that the perturbative approaches proposed in [A. Anantharaman and C. Le Bris, CRAS 2010, MMS 2011] can be turned into an efficient control variable. We theoretically demonstrate the efficiency of our approach in simple cases. We next provide illustrating numerical results and compare our approach with other variance reduction strategies. We also show how to use the Reduced Basis approach proposed in [C. Le Bris and F. Thomines, Chinese Ann. Math. 2012] so that the cost of building the surrogate model remains limited.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Periodic homogenization using the Lippmann--Schwinger formalism

When homogenizing elliptic partial differential equations, the so-called corrector problem is pivotal to compute the macroscale effective coefficients from the microscale information. To solve this corrector problem in the periodic setting, Moulinec and Suquet introduced in the mid-nineties a numerical strategy based on the reformulation of that problem as an integral equation (known as the Lippmann--Schwinger equation), which is then suitably discretized. This results in an iterative, matrix-free method, which is of particular interest for complex microstructures. Since the seminal work of Moulinec and Suquet, several variants of their scheme have been proposed. The aim of this contribution is twofold. First, we provide an overview of these methods, recast in the language of the applied mathematics community. These methods are presented as asymptotically consistent Galerkin discretizations of the Lippmann--Schwinger equation. The bilinear form arising in the weak form of this integral equation is indeed the sum of a local and a non-local term. We show that most of the variants proposed in the literature correspond to alternative approximations of this non-local term. Second, we propose a mathematical analysis of the discretized problem. In particular, we prove under mild hypotheses the convergence of these numerical schemes with respect to the grid-size. We also provide a priori error estimates on the solution. The article closes on a three-dimensional numerical application within the framework of linear elasticity.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A micro-macro parareal algorithm: application to singularly perturbed ordinary differential equations

We introduce a micro-macro parareal algorithm for the time-parallel integration of multiscale-in-time systems. The algorithm first computes a cheap, but inaccurate, solution using a coarse propagator (simulating an approximate slow macroscopic model), which is iteratively corrected using a fine-scale propagator (accurately simulating the full microscopic dynamics). This correction is done in parallel over many subintervals, thereby reducing the wall-clock time needed to obtain the solution, compared to the integration of the full microscopic model. We provide a numerical analysis of the algorithm for a prototypical example of a micro-macro model, namely singularly perturbed ordinary differential equations. We show that the computed solution converges to the full microscopic solution (when the parareal iterations proceed) only if special care is taken during the coupling of the microscopic and macroscopic levels of description. The convergence rate depends on the modeling error of the approximate macroscopic model. We illustrate these results with numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

An MsFEM type approach for perforated domains

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Effective dynamics for a kinetic Monte-Carlo model with slow and fast time scales

We consider several multiscale-in-time kinetic Monte Carlo models, in which some variables evolve on a fast time scale, while the others evolve on a slow time scale. In the first two models we consider, a particle evolves in a one-dimensional potential energy landscape which has some small and some large barriers, the latter dividing the state space into metastable regions. In the limit of infinitely large barriers, we identify the effective dynamics between these macro-states, and prove the convergence of the process towards a kinetic Monte Carlo model. We next consider a third model, which consists of a system of two particles. The state of each particle evolves on a fast time-scale while conserving their respective energy. In addition, the particles can exchange energy on a slow time scale. Considering the energy of the first particle, we identify its effective dynamics in the limit of asymptotically small ratio between the characteristic times of the fast and the slow dynamics. For all models, our results are illustrated by representative numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Multi-Level Monte Carlo approaches for numerical homogenization

In this article, we study the application of Multi-Level Monte Carlo (MLMC) approaches to numerical random homogenization. Our objective is to compute the expectation of some functionals of the homogenized coefficients, or of the homogenized solutions. This is accomplished within MLMC by considering different levels of representative volumes (RVE), and, when it comes to homogenized solutions, different levels of coarse-grid meshes. Many inexpensive computations with the smallest RVE size and the largest coarse mesh are combined with fewer expensive computations performed on larger RVEs and smaller coarse meshes. We show that, by carefully selecting the number of realizations at each level, we can achieve a speed-up in the computations in comparison to a standard Monte Carlo method. Numerical results are presented both for one-dimensional and two-dimensional test-cases.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Variance reduction using antithetic variables for a nonlinear convex stochastic homogenization problem

We consider a nonlinear convex stochastic homogenization problem, in a stationary setting. In practice, the deterministic homogenized energy density can only be approximated by a random apparent energy density, obtained by solving the corrector problem on a truncated domain. We show that the technique of antithetic variables can be used to reduce the variance of the computed quantities, and thereby decrease the computational cost at equal accuracy. This leads to an efficient approach for approximating expectations of the apparent homogenized energy density and of related quantities. The efficiency of the approach is numerically illustrated on several test cases. Some elements of analysis are also provided.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A numerical strategy for coarse-graining two-dimensional atomistic models at finite temperature: the membrane case

We present a numerical strategy to compute ensemble averages of coarse-grained two-dimensional membrane-like models. The approach consists in generalizing to these two-dimensional models a one-dimensional strategy exposed in [Blanc, Le Bris, Legoll, Patz, JNLS 2010], which is based on applying the ergodic theorem to Markov chains. This may be considered as a first step towards computing the constitutive law associated to such models, in the thermodynamic limit.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Negative thermal conductivity of chains of rotors with mechanical forcing

We consider chains of rotors subjected to both thermal and mechanical forcings,in a nonequilibrium steady-state. Unusual nonlinear profiles of temperature and velocities are observed in the system. In particular, the temperature is maximal in the center, which is an indication of the nonlocal behavior of the system. In spite of that, local equilibrium holds for long enough chains. Our numerical results also show that, when the mechanical forcing is strong enough, the energy current can be increased by an inverse temperature gradient. This counterintuitive result again reveals the complexity of nonequilibrium states.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Free energy calculations: An efficient adaptive biasing potential method

We develop an efficient sampling and free energy calculation technique within the adaptive biasing potential (ABP) framework. By mollifying the density of states we obtain an approximate free energy and an adaptive bias potential that is computed directly from the population along the coordinates of the free energy. Because of the mollifier, the bias potential is "nonlocal" and its gradient admits a simple analytic expression. A single observation of the reaction coordinate can thus be used to update the approximate free energy at every point within a neighborhood of the observation. This greatly reduces the equilibration time of the adaptive bias potential. This approximation introduces two parameters: strength of mollification and the zero of energy of the bias potential. While we observe that the approximate free energy is a very good estimate of the actual free energy for a large range of mollification strength, we demonstrate that the errors associated with the mollification may be removed via deconvolution. The zero of energy of the bias potential, which is easy to choose, influences the speed of convergence but not the limiting accuracy. This method is simple to apply to free energy or mean force computation in multiple dimensions and does not involve second derivatives of the reaction coordinates, matrix manipulations nor on-the-fly adaptation of parameters. For the alanine dipeptide test case, the new method is found to gain as much as a factor of ten in efficiency as compared to two common adaptive biasing force formulations and it is shown to be as efficient as well-tempered metadynamics with the post-process deconvolution giving a clear advantage to the mollified density of states method.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Some remarks on free energy and coarse-graining

We present recent results on coarse-graining techniques for thermodynamic quantities (canonical averages) and dynamical quantities (averages of path functionals over solutions of overdamped Langevin equations). The question is how to obtain reduced models to compute such quantities, in the specific case when the functional to be averaged only depends on a few degrees of freedom. We mainly review, numerically illustrate and extend results from [3,18], concerning the computation of the stress-strain relation for one-dimensional chains of atoms, and the construction of an effective dynamics for a scalar coarse-grained variable when the complete system evolves according to the overdamped Langevin equation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Symplectic schemes for highly oscillatory Hamiltonian systems: the homogenization approach beyond the constant frequency case

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Effective dynamics using conditional expectations

The question of coarse-graining is ubiquitous in molecular dynamics. In this article, we are interested in deriving effective properties for the dynamics of a coarse-grained variable $ξ(x)$, where $x$ describes the configuration of the system in a high-dimensional space $\R^n$, and $ξ$ is a smooth function with value in $\R$ (typically a reaction coordinate). It is well known that, given a Boltzmann-Gibbs distribution on $x \in \R^n$, the equilibrium properties on $ξ(x)$ are completely determined by the free energy. On the other hand, the question of the effective dynamics on $ξ(x)$ is much more difficult to address. Starting from an overdamped Langevin equation on $x \in \R^n$, we propose an effective dynamics for $ξ(x) \in \R$ using conditional expectations. Using entropy methods, we give sufficient conditions for the time marginals of the effective dynamics to be close to the original ones. We check numerically on some toy examples that these sufficient conditions yield an effective dynamics which accurately reproduces the residence times in the potential energy wells. We also discuss the accuracy of the effective dynamics in a pathwise sense, and the relevance of the free energy to build a coarse-grained dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Thermal conductivity of the Toda lattice with conservative noise

We study the thermal conductivity of the one dimensional Toda lattice perturbed by a stochastic dynamics preserving energy and momentum. The strength of the stochastic noise is controlled by a parameter $γ$. We show that heat transport is anomalous, and that the thermal conductivity diverges with the length $n$ of the chain according to $κ(n) \sim n^α$, with $0 < α\leq 1/2$. In particular, the ballistic heat conduction of the unperturbed Toda chain is destroyed. Besides, the exponent $α$ of the divergence depends on $γ$.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Non-ergodicity of Nose-Hoover dynamics

The numerical integration of the Nose-Hoover dynamics gives a deterministic method that is used to sample the canonical Gibbs measure. The Nose-Hoover dynamics extends the physical Hamiltonian dynamics by the addition of a "thermostat" variable, that is coupled nonlinearly with the physical variables. The accuracy of the method depends on the dynamics being ergodic. Numerical experiments have been published earlier that are consistent with non-ergodicity of the dynamics for some model problems. The authors recently proved the non-ergodicity of the Nose-Hoover dynamics for the one-dimensional harmonic oscillator. In this paper, this result is extended to non-harmonic one-dimensional systems. It is also shown for some multidimensional systems that the averaged dynamics for the limit of infinite thermostat "mass" have many invariants, thus giving theoretical support for either non-ergodicity or slow ergodization. Numerical experiments for a two-dimensional central force problem and the one-dimensional pendulum problem give evidence for non-ergodicity.

0 citations0 reviewsNo code linkedOpen access