NADec 28, 2015
Variance Reduction using Nonreversible Langevin SamplersA. B. Duncan, T. Lelievre, G. A. Pavliotis
A standard approach to computing expectations with respect to a given target measure is to introduce an overdamped Langevin equation which is reversible with respect to the target distribution, and to approximate the expectation by a time-averaging estimator. As has been noted in recent papers, introducing an appropriately chosen nonreversible component to the dynamics is beneficial, both in terms of reducing the asymptotic variance and of speeding up convergence to the target distribution. In this paper we present a detailed study of the dependence of the asymptotic variance on the deviation from reversibility. Our theoretical findings are supported by numerical simulations.
NANov 4, 2008
Results and questions on a nonlinear approximation approach for solving high-dimensional partial differential equationsC. Le Bris, T. Lelievre, Y. Maday
We investigate mathematically a nonlinear approximation type approach recently introduced in [A. Ammar et al., J. Non-Newtonian Fluid Mech., 2006] to solve high dimensional partial differential equations. We show the link between the approach and the greedy algorithms of approximation theory studied e.g. in [R.A. DeVore and V.N. Temlyakov, Adv. Comput. Math., 1996]. On the prototypical case of the Poisson equation, we show that a variational version of the approach, based on minimization of energies, converges. On the other hand, we show various theoretical and numerical difficulties arising with the non variational version of the approach, consisting of simply solving the first order optimality equations of the problem. Several unsolved issues are indicated in order to motivate further research.
NAApr 9, 2008
Generalized Navier Boundary Condition and Geometric Conservation Law for surface tensionJ. -F. Gerbeau, T. Lelievre
We consider two-fluid flow problems in an Arbitrary Lagrangian Eulerian (ALE) framework. The purpose of this work is twofold. First, we address the problem of the moving contact line, namely the line common to the two fluids and the wall. Second, we perform a stability analysis in the energy norm for various numerical schemes, taking into account the gravity and surface tension effects. The problem of the moving contact line is treated with the so-called Generalized Navier Boundary Conditions. Owing to these boundary conditions, it is possible to circumvent the incompatibility between the classical no-slip boundary condition and the fact that the contact line of the interface on the wall is actually moving. The energy stability analysis is based in particular on an extension of the Geometry Conservation Law (GCL) concept to the case of moving surfaces. This extension is useful to study the contribution of the surface tension. The theoretical and computational results presented in this paper allow us to propose a strategy which offers a good compromise between efficiency, stability and artificial diffusion.