136 research outputs found
Quelques modèles effectifs pour la modélisation de matériaux à l'échelle atomistique à température finie
National audienceSee http://hal.archives-ouvertes.fr/docs/00/59/28/49/ANNEX/r_X1U8C757.pd
Contributions à l'étude mathématique et numérique de quelques modèles en simulation multi-échelle des matériaux
The first part of the manuscript is concerned with various works in molecular simulation. We consider systems composed of point particles (that typically represent atoms in a molecular system), that interact through a potential energy. The degrees of freedom of the system are the position and the momentum of each particle. The very large number of degrees of freedom makes the problem challenging: many systems considered by the applied community are composed of more than 10^5 atoms, and the largest systems that can today be simulated contain a few tens of billions of atoms. A first question is to sample the Boltzmann-Gibbs measure associated to the potential energy of the system. Several methods have been proposed in the literature, based on Markov chains, Markov processes (typically solutions of stochastic differential equations), and ordinary differential equations. We briefly review these methods, and present our results on the non-ergodicity of some deterministic methods (namely, the Nosé-Hoover method). One has hence to be cautious when using this method to sample the Boltzmann-Gibbs measure. We next consider questions related to the construction of effective dynamics. Typically, the potential energy has many local minima, separated by high barriers. However, it is not always needed to know all the degrees of freedom of the system to know in which conformation the system is. The knowledge of a function of the degrees of freedom is often sufficient enough. Starting from the microscopic description of the system (where the state is described by a vector X, the dimension of which is large) we then introduce a macroscopic description, based on a function of the microscopic state, that we denote here xi(X). This function (that we assume to be scalar-valued) is our quantity of interest. We propose a strategy to design a dynamics that approximates the dynamics xi(X_t) of the quantity of interest, when the system, at the microscopic scale, evolves according to the overdamped Langevin equation. Under a time scale separation assumption, we show that the dynamics we propose is an accurate approximation of xi(X_t). The works described above have been performed in a finite temperature setting. We have also considered isolated systems, the evolution of which is given by Newton laws, that we write as a Hamiltonian dynamical system. Simulating this dynamics is challenging, due to the gap between the fast time scales present in the system (vibrations of some chemical bonds occur with a typical period of the order of the femtosecond), and the slow time scales, which are of the order of the microsecond or the millisecond. Such simulation times are essentially out of reach if one uses a standard integrator, that would use time steps of the order of the femtosecond. We have addressed this problem in a two-fold manner. First, we have proposed numerical integrators for highly oscillatory Hamiltonian dynamics (for which the time step is not limited by the fastest characteristic time scales), following a homogenization (in time) approach. We have also proposed a variant of the parareal scheme (that uses parallel computations to simulate a dynamics) that is better adapted to the Hamiltonian framework than the original parareal algorithm. In the second part of the manuscript, considering solids described at the atomistic scale, we derive coarse-grained models (written at the continuum scale) and couple these two models. The fine scale model is similar to models considered in the first part of the manuscript, whereas the coarse-grained model is written at a much more macroscopic scale. During the PhD, we have addressed this question in a variational setting (corresponding to a vanishing temperature). The configuration we are after is a (sometimes unique) minimizer of some energy. We have then studied how to couple the atomistic model with the corresponding continuum model. Recently, we have turned to a finite temperature setting, in which macroscopic quantities of interest are averages of functions depending of the microscopic state of the system (the position of all the atoms) with respect to the Boltzmann measure. In this setting, we have derived coarse-grained models, where the temperature is a parameter, following thermodynamic limit approaches. In the third part of the manuscript, we turn to materials modelled at the continuum scale by (linear) elliptic partial differential equations. We consider random materials, the properties of which oscillate at a fast scale. Random homogenization is a well developed theory to handle such problems. However, even of such simple cases, the associated numerical methods today available (such as methods to compute the homogenized matrix) generally lead to very expensive computations. We address this problem from a numerical viewpoint: we consider an equation that is simple from the theoretical viewpoint, and wish to design more efficient numerical strategies. We have followed two directions. The first one is concerned with variance reduction. In the setting we consider, even though the exact homogenized matrix is deterministic, it turns out that, due to the numerical procedure that is standardly used, one can only have access to a random approximation of the homogenized matrix. We have proposed variance reduction techniques, to typically obtain an approximation with a smaller confidence interval, in the framework of Monte Carlo methods. A second direction consists in considering weakly stochastic materials. Indeed, real materials are rarely periodic, but they are not always strongly random. The case when the randomness comes as a small perturbation of a periodic model is thus relevant, in practice, to model a large class of materials. We then show that, in this case, we can compute, at the first orders, the homogenized matrix with a workload comparable to that encountered in periodic homogenization, and hence much smaller than the one required by generic stochastic homogenization.La première partie du mémoire résume des travaux en simulation moléculaire. On s'intéresse à des systèmes de particules ponctuelles (représentant typiquement les noyaux des atomes d'un système moléculaire), qui interagissent via une énergie potentielle. Les degrés de liberté du système sont la position et l'impulsion de chaque particule. La complexité du problème vient du nombre de degrés de liberté en jeu, qui peut atteindre (et dépasser!) plusieurs centaines de milliers d'atomes pour les systèmes d'intérêt pratique. Les questions étudiées portent sur l'échantillonnage de la mesure de Boltzmann-Gibbs (avec des résultats concernant la non-ergodicité de certains systèmes dynamiques proposés dans la littérature), et sur la construction de dynamiques effectives: supposant que le système suit une dynamique X_t régie par l'équation de Langevin amortie, et se donnant une variable scalaire macroscopique xi(X), lente en un certain sens, nous proposons une dynamique mono-dimensionnelle fermée qui approche xi(X_t), et dont la précision est estimée à l'aide de méthodes d'entropie relative. Une autre partie du travail consiste à développer de nouveaux schémas numériques pour des problèmes Hamiltoniens hautement oscillants (souvent rencontrés en simulation moléculaire), en suivant une démarche d'homogénéisation en temps. Nous avons aussi proposé une adaptation au contexte Hamiltonien de l'algorithme pararéel, permettant d'obtenir la solution d'un problème d'évolution par des méthodes de calcul parallèle. La seconde partie du mémoire présente des travaux sur la dérivation de modèles à l'échelle du continuum à partir de modèles discrets (à l'échelle atomistique), pour les solides, et sur le couplage de ces deux modèles, discret et continu. Une première approche consiste à poser le problème sous forme variationnelle (modélisation à température nulle). Nous nous sommes aussi intéressés au cas de systèmes à température finie, modélisés dans le cadre de la mécanique statistique. Dans certains cas, nous avons obtenu des modèles réduits, macroscopiques, où la température est un paramètre, en suivant des approches de type limite thermodynamique. La troisième partie du mémoire s'intéresse à des questions d'homogénéisation stochastique, pour des équations aux dérivées partielles elliptiques linéaires. Les matériaux sont donc modélisés à l'échelle du continuum. Le constat qui motive notre travail est le fait que, même dans les cas les plus simples sur le plan théorique, les méthodes numériques à ce jour disponibles en homogénéisation stochastique conduisent à des calculs très lourds. Nous avons travaillé dans deux directions. La première consiste à réduire la variance des quantités aléatoires effectivement calculées, seules accessibles en pratique pour approcher la matrice homogénéisée. La seconde est d'étudier le cas de problèmes faiblement stochastiques, en partant du constat que les matériaux hétérogènes, rarement périodiques, ne sont pas pour autant systématiquement fortement aléatoires. Le cas d'un matériau aléatoire pour lequel cet aléa n'est qu'une petite perturbation autour d'un modèle périodique est donc intéressant, et peut se traiter avec un coût calcul beaucoup plus abordable
Derivation of Langevin Dynamics in a Nonzero Background Flow Field
We propose a derivation of a nonequilibrium Langevin dynamics for a large
particle immersed in a background flow field. A single large particle is placed
in an ideal gas heat bath composed of point particles that are distributed
consistently with the background flow field and that interact with the large
particle through elastic collisions. In the limit of small bath atom mass, the
large particle dynamics converges in law to a stochastic dynamics. This
derivation follows the ideas of [D. D\"urr, S. Goldstein, and J. L. Lebowitz,
1981 and 1983; P. Calderoni, D. D\"urr, and S. Kusuoka, 1989] and provides
extensions to handle the nonzero background flow. The derived nonequilibrium
Langevin dynamics is similar to the dynamics in [M. McPhie, et al., 2001]. Some
numerical experiments illustrate the use of the obtained dynamic to simulate
homogeneous liquid materials under flow.Comment: Minor revisions, refined discussion of the laminar bath approach and
non-Hamiltonian dynamics approach in Section 2. 41 pages, 8 figure
An accurate scheme to solve cluster dynamics equations using a Fokker-Planck approach
We present a numerical method to accurately simulate particle size
distributions within the formalism of rate equation cluster dynamics. This
method is based on a discretization of the associated Fokker-Planck equation.
We show that particular care has to be taken to discretize the advection part
of the Fokker-Planck equation, in order to avoid distortions of the
distribution due to numerical diffusion. For this purpose we use the
Kurganov-Noelle-Petrova scheme coupled with the monotonicity-preserving
reconstruction MP5, which leads to very accurate results. The interest of the
method is highlighted on the case of loop coarsening in aluminum. We show that
the choice of the models to describe the energetics of loops does not
significantly change the normalized loop distribution, while the choice of the
models for the absorption coefficients seems to have a significant impact on
it
Integrators for highly oscillatory Hamiltonian systems: an homogenization approach
International audienceWe introduce a systematic way to construct symplectic schemes for the numerical integration of a large class of highly oscillatory Hamiltonian systems. The bottom line of our construction is to consider the Hamilton-Jacobi form of the Newton equations of motion, and to perform a two-scale expansion of the solution, for small times and high frequencies. The approximation obtained for the solution is then used as a generating function, from which the numerical scheme is derived. Several options for the derivation are presented, some of them also giving rise to non symplectic variants. The various integrators obtained are tested and compared to several existing algorithms. The numerical results demonstrate their efficiency
Asymptotic Behavior of Green Functions of Divergence form Operators with Periodic Coefficients
International audienceThis article is concerned with the asymptotic behavior, at infinity and at the origin, of Green functions of operators of the form Lu= −div(A∇u), where A is a periodic, coercive, and bounded matrix
Réduction de variance en homogénéisation stochastique : l'exemple des variables antithétiques
We show that we can reduce the variance in a simple problem of stochastic homogenization using the classical technique of antithetic variables. The setting, and the presentation, are deliberately kept elementary. We point out the main issues, show some illustrative results, and demonstrate, both theoretically and numerically, the efficiency of the approach on simple cases
Stochastic homogenization for a scalar viscoelastic model exhibiting stress–strain hysteresis
Motivated by rate-independent stress-strain hysteresis observed in filled rubber, this article considers a scalar viscoelastic model in which the constitutive law is random and varies on a lengthscale which is small relative to the overall size of the solid. Using stochastic two-scale convergence as introduced by Bourgeat, Mikelic and Wright, we obtain the homogenized limit of the evolution, and demonstrate that under certain hypotheses, the homogenized model exhibits hysteretic behaviour which persists under asymptotically slow loading. These results are illustrated by means of numerical simulations in a particular one-dimensional instance of the model
Finite-temperature coarse-graining of one-dimensional models: mathematical analysis and computational approaches
Published in "Journal of Nonlinear Science", 2010, vol. 20, no 2, p. 241-275 http://hal-enpc.archives-ouvertes.fr/hal-00625882/fr/International audienceWe present a possible approach for the computation of free energies and ensemble averages of one-dimensional coarse-grained models in materials science. The approach is based upon a thermodynamic limit process, and makes use of ergodic theorems and large deviations theory. In addition to providing a possible efficient computational strategy for ensemble averages, the approach allows for assessing the accuracy of approximations commonly used in practice
- …