2 Dissipative particle dynamics with energy conservation: Heat conduction that is dissipated due to the velocity dependent forces is transformed into
|
|
- Dwayne Hutchinson
- 5 years ago
- Views:
Transcription
1 International Journal of Modern Physics C, fc World Scientic Publishing Company DISSIPATIVE PARTICLE DYNAMICS WITH ENERGY CONSERVATION: HEAT CONDUCTION M. RIPOLL, P. ESPA ~NOL Departamento de Fsica Fundamental, UNED, Apartado 604 Madrid 28080, Spain and M.H. ERNST Institute for Theoretical Physics, Utrecht University, Princetonplein 5, P.O. Box TA Utrecht, The Netherlands Received (received date) Revised (revised date) We study by means of numerical simulations the model of dissipative particle dynamics with energy conservation for the simple case of thermal conduction. It is shown that the model displays correct equilibrium uctuations and reproduces Fourier law. The connection between \mesoscopic coarse-graining" and \resolution" is claried.. Introduction Mesoscopic approaches to the study of complex uids are receiving a great deal of attention in view of their potential to address problems with disparate time-scales. The method of dissipative particle dynamics introduced originally as an hybrid between lattice gas and molecular dynamics simulations, has proven to be a exible and powerful tool in the study of complex uids as colloidal suspensions 2;;4, porous ow, polymer suspensions 5, and multicomponent ows 6. DPD models a uid as a collection of point particles that are interpreted as representing the center of mass of a mesoscopic cluster of the atoms or molecules that constitute the uid. By formulating physically motivated interactions (that is, dissipative forces that conserve momentum ;7 ) between these \lumps" of uid, it is assured that the macroscopic behavior of the system is hydrodynamic 8;9. Most importantly, a random interaction is introduced in order to account for thermal uctuations that occur at the mesoscopic level and are responsible for Brownian eects appearing in complex uids. One of the drawbacks of classic DPD is that the total energy of the system is not conserved because the forces are velocity dependent. This has been recently remedied simultaneously and independently in Refs. 0;, where an internal energy variable and a temperature is introduced for each particle. The mechanical energy pep@sfun.uned.es
2 2 Dissipative particle dynamics with energy conservation: Heat conduction that is dissipated due to the velocity dependent forces is transformed into internal energy of the particles. This viscous heating process is accompanied by a thermal conduction process where exchange of internal energy between neighboring particles occurs when dierences of temperatures exist. The new DPD algorithm with energy conservation opens up the possibility of studying not only rheological aspects of complex uids but also thermal issues. It is therefore important to validate the technique in simple situations for which the dynamical aspects are well-known, either analytically or by means of other numerical simulations techniques. As a rst step towards the understanding of the behavior of DPD with energy conservation we focus in this paper on the simplest problem of conduction in a quiescent uid (or a solid). This is a particular case of the model introduced in Refs. 0;. The particles are at rest and located randomly and represent portions of material that can interchange energy. This simple model can be regarded as a way of solving a uctuating heat conduction equation. In section 2 we review the model that was introduced in Refs. 0;, for the case that only heat conduction takes place. In section we further specify the model functions. In section 4 we present simulation results that show the validity of the simulation method by comparing the numerical results with the theoretical results. Finally, we end up with a section of conclusions. 2. The Model We have introduced in Ref. 0 a mesoscopic model for describing a uctuating viscous and thermally conducting uid. Here we apply this model to a quiescent uid (or solid) within a box. The system is represented by a set of N particles of xed positions r i which are distributed at random within the box. These particles are understood as thermodynamic subsystems of the whole system, that is, small portions of material that have a suciently large number of degrees of freedom. To each particle we associate an internal energy variable i an entropy function s( i ), and a temperature T i? i =@ i. Each particle represents a mesoscopic portion of the material and its internal energy content is subject to thermal uctuations, which in a continuum theory would be described by a random heat ux. 2 Dierences in temperatures between neighboring particles will cause a transfer of internal energy. Therefore, the following stochastic dierential equation for the internal energy of each particle is postulated d i = X j!(r ij ) ij T i? T j dt + X j (2k B!(r ij ) ij ) =2 dw ij () The rst term in the r.h.s. is deterministic and species that a temperature difference causes exchange of energy. The second term is stochastic and takes into account thermally induced uctuations in each particle caused by a random interchange of energy between particles. The notation is as follows. k B is Boltzmann's constant, r ij = jr i? r j j is the distance between particles i; j and the dimensionless weight function!(r) determines the range of inuence between particles. The pa-
3 Dissipative particle dynamics with energy conservation: Heat conduction rameter ij governs the overall amplitude of the deterministic and random parts. It depends in general on the state variables of particles i; j and is symmetric under particle interchange. The random \heat ux" from particle i to j is expressed in terms of the increments of the Wiener process dwij. The increments of the Wiener process are antisymmetric under particle interchange dwij =?dw ji in such a way that the total internal energy of the system governed by Eqn. () is exactly conserved d dt i i = 0. Because ij might depend in general on i ; j we need to P provide a stochastic interpretation of Eqn. (), as the noise is multiplicative. We choose It^o interpretation. It has been shown that, ij ij = 0 y, then Eqn. () is mathematically equivalent to a Fokker-Planck equation which has as equilibrium solution eq () / exp ( k? B X i s( i ) ) P( X k k ) (2) P where P( k k) is a function of the total internal energy which is determined by the initial distribution of total energy. Particular examples are given by the canonical and microcanonical ensembles mic ( ; : : : ; N ) = can ( ; : : : ; N ) = ( (E 0 ; N) exp ( Z(; N) exp k? B k? B X i X i s( i ) ) ( X i s( i )? i ) i? E 0 ) Here, the factors (E 0 ; N) and Z(; N) are obtained by normalization.. Perfect solid and Fourier law In this section we further specify the model by providing the undetermined functions T (), ( + 0 ) and!(r) appearing in Eqn. (). The simplest possible form for the equation of state T () is that of a perfect solid, which is a very good approximation for metals. The perfect solid has the following equation of state T () = () C v (4) where the heat capacity of the mesoscopic particles C v is a constant independent of the energy. It is an extensive property that depends on the \size" of the particles. For mesoscopic particles, the dimensionless heat capacity C v =k B. The entropy s() is obtained by integrating the temperature and, apart from irrelevant constant additive terms, it is given by s() = C v ln. Now, let us assume the following functional form for ij, y This implies that if ij depends only on i ; j then, necessarily, it will be a function of i + j. This can be proved by considering a Taylor expansion of a function of two variables and using the fact that, at each point, both partial derivatives coincide.
4 4 Dissipative particle dynamics with energy conservation: Heat conduction ij = C v~ 2 T 2 i + j = C v~ 2 4 (T 2 i + T j ) 2 (5) where is the average distance between particles (this is = n?= and n is the average density number) and ~ has dimensions of a diusion coecient. In the last equality, Eqn. (4) has been used. Substitution of Eqn. (5) into the dynamical equation for the energy leads to the following equation for the temperature of a particle 2kB ~ dt i = ~ 2 X j!(r ij ) (T i + T j ) 2 4T i T j (T j? T i )dt + C v 2 =2 X j! =2 (r ij ) T i + T j dwij 2 (6) (Ti+Tj )2 The factor 4T it j should be very close to if the temperature dierences between particles within an action sphere are small. We see, then, that for the perfect solid by a proper selection of ij the form of the postulated SDE for the energy leads to a discrete version of the Fourier law of heat conduction. Note that the factor 2 in the rst term of the left hand side of Eqn. (6) can be understood as the \lattice spacing" squared that would arise in a nite dierence discretization of the Laplacian in the heat conduction equation. In a sense, Eqn. (6) can be understood as a nite dierence discretization of the macroscopic heat equation on a random lattice which, in addition, has the correct form for the thermal uctuations. We will assume that the weight function!(r) has the form of the Lucy weight function, which is a smooth bell-shaped function with continuous derivatives. It is given by!(r) = 05 6s + r r c? r r c (7) where r c is range of interaction between particle and s = r c = is the overlapping coecient which gives, essentially, the number of neighbors of a given particle. We have normalized the weight function as Z d d r!(r) = n (8) 4. Simulation results In this section we present some results obtained from numerical simulation of Eqn. (6). The physical system is assumed to be in a three dimensional cubic box of edge length L. Initially the box is seeded with N mesoscopic particles located at random. The density is then n = N=L and the typical interparticle distance is = n?=. Periodic boundary conditions are assumed in the y; z directions and in the x direction we assume either periodic boundary conditions when considering equilibrium situations or we impose the temperature at x =?L=2 and x = L=2. This is achieved by considering two extra layers of particles in each boundary in the
5 Dissipative particle dynamics with energy conservation: Heat conduction 5 x direction lled with particles that are at a constant temperature. These layers act, therefore, as thermal baths. They have a width as large as r c in order that any particle within the system interacts with the same number of particles. The basic model parameters are the following ones: ~, N, L, r c, C v, T 0. Here, T 0 is a reference temperature, that of equilibrium or, in the case of nonequilibrium situations, that of the colder thermal bath. We work in reduced units such that L is the unit of length, L 2 =~ is unit of time, T 0 is the unit of temperature, and k B is taken as the unit of entropy or heat capacity. We introduce the dimensionless heat capacity of mesoscopic particles = C v =k B and the dimensionless heat capacity of the total simulated sample, which is C L = N because heat capacity is an extensive property. On the other hand, C L = L c v =k B, where c v is the heat capacity per unit volume, which is a material constant. For example, c v = : J=Km for copper at room temperature. Therefore, by xing C L we are xing the actual volume of the sample. For copper, a value :60?0 m 0:. In most of the simulations we use C L = 0 8 and this xes the value of depending on the number of particles in the system. Note that by xing C L, large N implies small, in accordance with the idea that the heat capacity of the mesoscopic particles is proportional to its typical size. We note that the overall intensity of the noise in (6) is proportional to (k B =C v ) =2 =?=2 and, therefore, inversely proportional to the square root of the volume of the particles. 4 The higher the number density of mesoscopic particles used to discretize a given sample, the smaller is the size of such particles and, therefore, the larger is the thermal noise. Eqns. (6) are solved with a conventional Euler method. For values of < 0, the stochastic term produces from time to time a temperature which is negative for some particle. This has disastrous consequences because the factor accompanying (T i? T j ) in the deterministic term of (6) is negative and produces a ow of energy from the cold particles to the hot particles. When this happens the system becomes unstable. However, this never occurs for suciently large values of as the ones we typically use and no instabilities are observed. In order to check that the code performs properly, we rst run a simulation with fully periodic boundary conditions. Initially, all the particles are assumed to have the same temperature T =. After some equilibration period the system reaches equilibrium. The energy is exactly conserved (no round-o errors) because we use a Verlet list 5, in such a way that what is gained by a particle is exactly what is lost by the rest of particles. In Fig. we show the probability density distribution of the energy of a particle as obtained from simulation of Eqn. () for two dierent values of = C v =k B, the dimensionless heat capacity of the mesoscopic particles, for a system of N = 000. The larger the larger the average value of the energy and variance. However, the relative uctuations are smaller, proportional to N?=2. Also shown in Fig. are the canonical f can () and microcanonical f mic () theoretical results for the given values of. The analytical expressions can be obtained from () and are given by C L = 0 8 implies a submicron sample size of L = C = L
6 6 Dissipative particle dynamics with energy conservation: Heat conduction f can () = f mic () = M?( + ) () expf?g E 0 E 0? E 0 (N?)(+)? (9) where?(x) is the gamma function. The parameter can be explicitly computed for a perfect solid by requiring that the average energy of a particle is E 0 =N, where E 0 is the total energy of the isolated system. The result is = N( + )=E 0. In the thermodynamic limit the canonical and microcanonical distributions are identical. In practice, they are indistinguishable for values N > 50. The perfect coincidence of the theoretical and simulation results in Fig. gives condence on the implementation of the numerical code used f() Fig.. The probability density distribution f() of the energy of a given particle. Two simulations with = 50 (left) and = 00 (right) are presented. Other parameters are N = 000, s = :5, T eq =. Also shown, although not visible because they are perfectly superimposed, are the theoretical predictions for the microcanonical and canonical distributions for the given values of. Next, we compute the thermal diusivity by a macroscopic measurement. Two heat baths are constructed with temperatures T cold and T hot in such a way that a gradient (T hot? T cold )=L is applied. A steady state is reached for suciently long times. Fig. 2 shows the average of the temperature once the steady state has been reached for dierent number of particles for the case that T cold =, T hot = 2 and the overlapping coecient is s = :5. For N = 00 the spatial noise is considerable. Fluctuations in this graph are not due to statistical noise but to the
7 Dissipative particle dynamics with energy conservation: Heat conduction 7 spatial inhomogeneities that appear for small number of particles and they would smear out by averaging over dierent realizations of the initial spatial conguration of the particles. For larger number of particles a smooth linear prole is achieved, in accordance with our expectation that the discrete model reproduces Fourier law. The temperature eld T (x) is obtained by binning the x axis and averaging the temperatures of the particles within each bin of volume xl T (x) x Fig. 2. Average temperature of the particles in each of the 25 bins in which the x axis of the box has been divided. Solid line is the theoretical Fourier prole. Open circles correspond to N = 00 and solid circles correspond to N = 000 particles within the box. Other parameters are s = :5, C L = 0 8. We have performed a series of simulations with dierent temperature gradients and have computed the macroscopic heat ux in the steady state. Two dierent ways of computing the macroscopic heat ux have been used. In the rst case, the macroscopic heat ux is obtained by computing at each time step the energy gained by the cold bath (which is equal to the energy lost by the hot bath) and dividing it by the time step, this is E=t. The second way follows from the microscopic denition of the heat ux as obtained from projection operator methods or from kinetic theory 6 X Q (0) ij!(r ij ) ij T i? T j rij 2 All transport uxes are caused by collisional transfer. As the particles do no move, there are no kinetic uxes. This second form produces better results because it involves a very large number of pairs of particles and it is therefore adopted here. In Fig. we plot the macroscopic heat ux in terms of the applied temperature gradient. A linear dependence is obtained, which allows to extract the thermal
8 8 Dissipative particle dynamics with energy conservation: Heat conduction diusivity from the slope Q x (rt ) x Fig.. Average x component of the heat ux as a function of the imposed temperature gradient. A linear dependence is obtained with small deviations at large gradients. The slope of the straight line gives the thermal diusivity. When this experiment is performed for dierent number densities and dierent overlapping we can compare with the results from kinetic theory as developed in Ref. 6. Fig. 4 shows the value of the macroscopic thermal diusivity as a function of the overlapping s. Also shown is the kinetic theory prediction for the thermal diusivity 6, reading, D = s2 ~ () 24 We observe that the theoretical prediction improves for large overlapping. The kinetic theory result () has been obtained through the choices Eqns. (4) and (5) and shows no dependence on the density. This is a consequence of introducing the factor?2 in Eqn. (5). We verify in Fig. 5 that simulations performed at dierent number density, with the remaining parameters s; C L held xed, show that the thermal diusivity D is indeed independent of n, as predicted in Eqn. (). 5. Conclusions In this paper we have applied the model introduced in Ref. 0 to the case of heat conduction in a random solid. We have chosen the equation of state s( i ) for a perfect solid model and a particular model for ij, Eqn. (5). It has been shown that the model has the correct equilibrium properties and that it reproduces Fourier's law in non-equilibrium situations. We have also corroborated, by means
9 Dissipative particle dynamics with energy conservation: Heat conduction D Fig. 4. Thermal diusivity coecient for dierent values of the overlapping. Also shown (solid line) is the kinetic theory prediction. s 2 of simulations, the kinetic theory predictions for this model which are presented elsewhere. Kinetic theory predicts a thermal diusivity which depends only on the overlapping and not on the density. By introducing the factor?2 in the denition of ij we have shown that it is possible to interpret n not as the density but as the resolution at which the physical problem is being studied. A similar discussion about resolution was presented for the case of DPD with no energy conservation in Ref. 7. Acknowledgments M.R. acknowledges the support of a F.P.I. grant from Ministerio de Educacion. P.E. has received partial support from DGICYT Project No PB and by E.C. Contract ERB-CHRXCT References. P.J. Hoogerbrugge and J.M.V.A. Koelman, Europhys. Lett. 9, 55 (992). 2. J.M.V.A. Koelman and P.J. Hoogerbrugge, Europhys. Lett. 2, 69 (99).. E.S. Boek, P.V. Coveney, H.N.W. Lekkerkerker, and P. van der Schoot, Phys. Rev. E 55, 24 (997). 4. E.S. Boek, P.V. Coveney, and H.N.W. Lekkerkerker, J. Phys.: Condens. Matter 8, 9509 (997). 5. A.G. Schlijper, P.J. Hoogerbrugge, and C.W. Manke, J. Rheol. 9, 567 (995). 6. P.V. Coveney and K. Novik, Phys. Rev. E 54, 54 (996). 7. P. Espa~nol and P. Warren, Europhys. Lett. 0, 9 (995).
10 0 Dissipative particle dynamics with energy conservation: Heat conduction D n Fig. 5. Thermal diusivity coecient for dierent values of the density. Also shown (solid line) is the kinetic theory prediction for this value of the overlapping s = P. Espa~nol, Phys. Rev. E, 52, 74 (995). 9. C. Marsh, G. Backx, and M.H. Ernst, Europhys. Lett. 8, 4 (997). C. Marsh, G. Backx, and M.H. Ernst, Phys. Rev. E 56, 976 (997). 0. P. Espa~nol, Europhys. Lett. 40, 6 (997).. J. Bonet Avalos and M. Mackie, Europhys. Lett. 40, 4 (997). 2. L.D. Landau and E.M. Lifshitz, Fluid Mechanics (Pergamon Press, 959).. L.B. Lucy, Astron. J. 82, 0 (977). 4. P. Espa~nol, Physica A 248, 77 (997). 5. M.P. Allen and D.J. Tildesley, Computer Simulation of Liquids (Clarendon Press, Oxford, 987). 6. M. Ripoll, M.H. Espa~nol, and M.H. Ernst, preprint. 7. P. Espa~nol, Phys. Rev. E 57, 290 (998).
APMA 2811T. By Zhen Li. Today s topic: Lecture 2: Theoretical foundation and parameterization. Sep. 15, 2016
Today s topic: APMA 2811T Dissipative Particle Dynamics Instructor: Professor George Karniadakis Location: 170 Hope Street, Room 118 Time: Thursday 12:00pm 2:00pm Dissipative Particle Dynamics: Foundation,
More informationMD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
MD Thermodynamics Lecture 1 3/6/18 1 Molecular dynamics The force depends on positions only (not velocities) Total energy is conserved (micro canonical evolution) Newton s equations of motion (second order
More informationApplication of the Dissipative Particle Dynamics Method to Ferromagnetic Colloidal Dispersions
Application of the Dissipative Particle Dynamics Method to Ferromagnetic Colloidal Dispersions AKIRA SATOH Professor, Akita Prefectural University, Japan Roy W. Chantrell Professor, York University, U.
More informationSmoothed Dissipative Particle Dynamics: theory and applications to complex fluids
2015 DPD Workshop September 21-23, 2015, Shanghai University Smoothed Dissipative Particle Dynamics: Dynamics theory and applications to complex fluids Marco Ellero Zienkiewicz Centre for Computational
More informationThe dynamics of small particles whose size is roughly 1 µmt or. smaller, in a fluid at room temperature, is extremely erratic, and is
1 I. BROWNIAN MOTION The dynamics of small particles whose size is roughly 1 µmt or smaller, in a fluid at room temperature, is extremely erratic, and is called Brownian motion. The velocity of such particles
More informationIntroduction to Computer Simulations of Soft Matter Methodologies and Applications Boulder July, 19-20, 2012
Introduction to Computer Simulations of Soft Matter Methodologies and Applications Boulder July, 19-20, 2012 K. Kremer Max Planck Institute for Polymer Research, Mainz Overview Simulations, general considerations
More informationMesoscopic simulation for the structural change of a surfactant solution using dissipative particle dynamics
Korean J. Chem. Eng., 26(6), 1717-1722 (2009) DOI: 10.1007/s11814-009-0235-2 RAPID COMMUNICATION Mesoscopic simulation for the structural change of a surfactant solution using dissipative particle dynamics
More informationUvA-DARE (Digital Academic Repository)
UvA-DARE (Digital Academic Repository) Poiseuille flow to measure the viscosity of particle model fluids. Backer, J.A.; Lowe, C.P.; Hoefsloot, H.C.J.; Iedema, P.D. Published in: Journal of Chemical Physics
More informationwith angular brackets denoting averages primes the corresponding residuals, then eq. (2) can be separated into two coupled equations for the time evol
This paper was published in Europhys. Lett. 27, 353{357, 1994 Current Helicity the Turbulent Electromotive Force N. Seehafer Max-Planck-Gruppe Nichtlineare Dynamik, Universitat Potsdam, PF 601553, D-14415
More informationA comparative study between dissipative particle dynamics and molecular dynamics for simple- and complex-geometry flows
THE JOURNAL OF CHEMICAL PHYSICS 123, 104107 2005 A comparative study between dissipative particle dynamics and molecular dynamics for simple- and complex-geometry flows Eric E. Keaveny, Igor V. Pivkin,
More informationA mathematical introduction to coarse-grained stochastic dynamics
A mathematical introduction to coarse-grained stochastic dynamics Gabriel STOLTZ gabriel.stoltz@enpc.fr (CERMICS, Ecole des Ponts & MATHERIALS team, INRIA Paris) Work also supported by ANR Funding ANR-14-CE23-0012
More informationAbsolute Equilibrium Entropy
To appear in the Journal of Plasma Physics Absolute Equilibrium Entropy By JOHN V. SHEBALIN 1 National Aeronautics and Space Administration, Langley Research Center, Hampton, Virginia 3681, USA The entropy
More informationDynamics that trigger/inhibit cluster formation in a one-dimensional granular gas
Physica A 342 (24) 62 68 www.elsevier.com/locate/physa Dynamics that trigger/inhibit cluster formation in a one-dimensional granular gas Jose Miguel Pasini a;1, Patricio Cordero b; ;2 a Department of Theoretical
More informationCOARSE-GRAINING AND THERMODYNAMICS IN FAR-FROM-EQUILIBRIUM SYSTEMS
Vol. 44 (2013) ACTA PHYSICA POLONICA B No 5 COARSE-GRAINING AND THERMODYNAMICS IN FAR-FROM-EQUILIBRIUM SYSTEMS J. Miguel Rubí, A. Pérez-Madrid Departament de Física Fonamental, Facultat de Física, Universitat
More informationAn alternative approach to dissipative particle dynamics
EUROPHYSICS LETTERS 15 July 1999 Europhys. Lett., 47 (2), pp. 145-151 (1999) An alternative approach to dissipative particle dynamics C. P. Lowe Computational Physics, Delft University of Technology Lorentzweg
More informationOrdering periodic spatial structures by non-equilibrium uctuations
Physica A 277 (2000) 327 334 www.elsevier.com/locate/physa Ordering periodic spatial structures by non-equilibrium uctuations J.M.G. Vilar a;, J.M. Rub b a Departament de F sica Fonamental, Facultat de
More informationChE 503 A. Z. Panagiotopoulos 1
ChE 503 A. Z. Panagiotopoulos 1 STATISTICAL MECHANICAL ENSEMLES 1 MICROSCOPIC AND MACROSCOPIC ARIALES The central question in Statistical Mechanics can be phrased as follows: If particles (atoms, molecules,
More informationAsymptotic solution of the Boltzmann equation for the shear ow of smooth inelastic disks
Physica A 275 (2000) 483 504 www.elsevier.com/locate/physa Asymptotic solution of the Boltzmann equation for the shear ow of smooth inelastic disks V. Kumaran Department of Chemical Engineering, Indian
More informationCoarse-graining limits in open and wall-bounded dissipative particle dynamics systems
THE JOURNAL OF CHEMICAL PHYSICS 124, 184101 2006 Coarse-graining limits in open and wall-bounded dissipative particle dynamics systems Igor V. Pivkin and George E. Karniadakis a Division of Applied Mathematics,
More informationContents. 1 Introduction and guide for this text 1. 2 Equilibrium and entropy 6. 3 Energy and how the microscopic world works 21
Preface Reference tables Table A Counting and combinatorics formulae Table B Useful integrals, expansions, and approximations Table C Extensive thermodynamic potentials Table D Intensive per-particle thermodynamic
More informationMesoscale fluid simulation of colloidal systems
Mesoscale fluid simulation of colloidal systems Mingcheng Yang Institute of Physics, CAS Outline (I) Background (II) Simulation method (III) Applications and examples (IV) Summary Background Soft matter
More informationFluctuations of intensive variables and non-equivalence of thermodynamic ensembles
Fluctuations of intensive variables and non-equivalence of thermodynamic ensembles A. Ya. Shul'man V.A. Kotel'nikov Institute of Radio Engineering and Electronics of the RAS, Moscow, Russia 7th International
More informationUniversity of Groningen. Mathematical modeling of complex systems Maurits, N.M.
University of Groningen Mathematical modeling of complex systems Maurits, N.M. MPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if you wish to cite from it. Please check
More informationThe goal of equilibrium statistical mechanics is to calculate the diagonal elements of ˆρ eq so we can evaluate average observables < A >= Tr{Â ˆρ eq
Chapter. The microcanonical ensemble The goal of equilibrium statistical mechanics is to calculate the diagonal elements of ˆρ eq so we can evaluate average observables < A >= Tr{Â ˆρ eq } = A that give
More informationDissipative nuclear dynamics
Dissipative nuclear dynamics Curso de Reacciones Nucleares Programa Inter universitario de Fisica Nuclear Universidad de Santiago de Compostela March 2009 Karl Heinz Schmidt Collective dynamical properties
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Statistical Physics I Spring Term 2013 Notes on the Microcanonical Ensemble
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department 8.044 Statistical Physics I Spring Term 2013 Notes on the Microcanonical Ensemble The object of this endeavor is to impose a simple probability
More informationMonte Carlo Simulation of the 2D Ising model
Monte Carlo Simulation of the 2D Ising model Emanuel Schmidt, F44 April 6, 2 Introduction Monte Carlo methods are a powerful tool to solve problems numerically which are dicult to be handled analytically.
More informationDiffusive Transport Enhanced by Thermal Velocity Fluctuations
Diffusive Transport Enhanced by Thermal Velocity Fluctuations Aleksandar Donev 1 Courant Institute, New York University & Alejandro L. Garcia, San Jose State University John B. Bell, Lawrence Berkeley
More informationChapter 4: Going from microcanonical to canonical ensemble, from energy to temperature.
Chapter 4: Going from microcanonical to canonical ensemble, from energy to temperature. All calculations in statistical mechanics can be done in the microcanonical ensemble, where all copies of the system
More informationarxiv:comp-gas/ v1 28 Apr 1993
Lattice Boltzmann Thermohydrodynamics arxiv:comp-gas/9304006v1 28 Apr 1993 F. J. Alexander, S. Chen and J. D. Sterling Center for Nonlinear Studies and Theoretical Division Los Alamos National Laboratory
More informationAPMA 2811T. By Zhen Li. Today s topic: Lecture 3: New Methods beyond traditional DPD. Sep. 22, Division of Applied Mathematics, Brown University
Today s topic: APMA 2811T Dissipative Particle Dynamics Instructor: Professor George Karniadakis Location: 170 Hope Street, Room 118 Time: Thursday 12:00pm 2:00pm Dissipative Particle Dynamics: Foundation,
More informationOn angular momentum balance for particle systems with periodic boundary conditions
On angular momentum balance for particle systems with periodic boundary conditions Vitaly A. Kuzkin 1 Institute for Problems in Mechanical Engineering RAS Saint Petersburg State Polytechnical University
More informationRemoving the mystery of entropy and thermodynamics. Part 3
Removing the mystery of entropy and thermodynamics. Part 3 arvey S. Leff a,b Physics Department Reed College, Portland, Oregon USA August 3, 20 Introduction In Part 3 of this five-part article, [, 2] simple
More informationParticle-Simulation Methods for Fluid Dynamics
Particle-Simulation Methods for Fluid Dynamics X. Y. Hu and Marco Ellero E-mail: Xiangyu.Hu and Marco.Ellero at mw.tum.de, WS 2012/2013: Lectures for Mechanical Engineering Institute of Aerodynamics Technical
More informationDissipative particle dynamics: Bridging the gap between atomistic and mesoscopic simulation
Dissipative particle dynamics: Bridging the gap between atomistic and mesoscopic simulation Robert D. Groot and Patrick B. Warren Unilever Research Port Sunlight Laboratory, Quarry Road East, Bebington,
More informationin order to insure that the Liouville equation for f(?; t) is still valid. These equations of motion will give rise to a distribution function f(?; t)
G25.2651: Statistical Mechanics Notes for Lecture 21 Consider Hamilton's equations in the form I. CLASSICAL LINEAR RESPONSE THEORY _q i = @H @p i _p i =? @H @q i We noted early in the course that an ensemble
More informationMolecular Dynamics Simulations
Molecular Dynamics Simulations Dr. Kasra Momeni www.knanosys.com Outline Long-range Interactions Ewald Sum Fast Multipole Method Spherically Truncated Coulombic Potential Speeding up Calculations SPaSM
More informationStochastic Particle Methods for Rarefied Gases
CCES Seminar WS 2/3 Stochastic Particle Methods for Rarefied Gases Julian Köllermeier RWTH Aachen University Supervisor: Prof. Dr. Manuel Torrilhon Center for Computational Engineering Science Mathematics
More informationModeling of Micro-Fluidics by a Dissipative Particle Dynamics Method. Justyna Czerwinska
Modeling of Micro-Fluidics by a Dissipative Particle Dynamics Method Justyna Czerwinska Scales and Physical Models years Time hours Engineering Design Limit Process Design minutes Continious Mechanics
More informationEntropy production fluctuation theorem and the nonequilibrium work relation for free energy differences
PHYSICAL REVIEW E VOLUME 60, NUMBER 3 SEPTEMBER 1999 Entropy production fluctuation theorem and the nonequilibrium work relation for free energy differences Gavin E. Crooks* Department of Chemistry, University
More informationMapping Closure Approximation to Conditional Dissipation Rate for Turbulent Scalar Mixing
NASA/CR--1631 ICASE Report No. -48 Mapping Closure Approximation to Conditional Dissipation Rate for Turbulent Scalar Mixing Guowei He and R. Rubinstein ICASE, Hampton, Virginia ICASE NASA Langley Research
More informationFinite-temperature magnetism of ultrathin lms and nanoclusters PhD Thesis Booklet. Levente Rózsa Supervisor: László Udvardi
Finite-temperature magnetism of ultrathin lms and nanoclusters PhD Thesis Booklet Levente Rózsa Supervisor: László Udvardi BME 2016 Background of the research Magnetic materials continue to play an ever
More information2 Introduction Standard conservative schemes will admit non-physical oscillations near material interfaces, which can be avoided by the application of
A Conservative Approach to the Multiphase Euler Equations without Spurious Pressure Oscillations Xu-Dong Liu Ronald P. Fedkiw Stanley Osher June 8, 998 Abstract In [4], the pressure evolution equation
More informationPhase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany
Phase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals Phase Equilibria Phase diagrams and classical thermodynamics
More informationThermodynamics of Three-phase Equilibrium in Lennard Jones System with a Simplified Equation of State
23 Bulletin of Research Center for Computing and Multimedia Studies, Hosei University, 28 (2014) Thermodynamics of Three-phase Equilibrium in Lennard Jones System with a Simplified Equation of State Yosuke
More informationMACROSCOPIC VARIABLES, THERMAL EQUILIBRIUM. Contents AND BOLTZMANN ENTROPY. 1 Macroscopic Variables 3. 2 Local quantities and Hydrodynamics fields 4
MACROSCOPIC VARIABLES, THERMAL EQUILIBRIUM AND BOLTZMANN ENTROPY Contents 1 Macroscopic Variables 3 2 Local quantities and Hydrodynamics fields 4 3 Coarse-graining 6 4 Thermal equilibrium 9 5 Two systems
More informationNon-equilibrium Stefan Boltzmann law
J. Non-Equilib. Thermodyn. 35 (2010), 279 288 DOI 10.1515/JNETDY.2010.17 de Gruyter 2010 Non-equilibrium Stefan Boltzmann law Agustin Pérez-Madrid, J. Miguel Rubí and Luciano C. Lapas Abstract We study
More informationEnsemble equivalence for non-extensive thermostatistics
Physica A 305 (2002) 52 57 www.elsevier.com/locate/physa Ensemble equivalence for non-extensive thermostatistics Raul Toral a;, Rafael Salazar b a Instituto Mediterraneo de Estudios Avanzados (IMEDEA),
More informationSmoothed Dissipative Particle Dynamics Model for Predicting Self-Assembled Nano-Cellulose Fibre Structures
Smoothed Dissipative Particle Dynamics Model for Predicting Self-Assembled Nano-Cellulose Fibre Structures David Vidal and Tetsu Uesaka FPInnovations, Pointe-Claire, Québec, CANADA Nano-cellulose fibres
More informationWhat is Classical Molecular Dynamics?
What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential functions Newton s equations of motion are integrated
More informationThermodynamics. Basic concepts. Thermal equilibrium and temperature
hermodynamics Basic concepts hermodynamics is a phenomenological description of macroscopic systems with many degrees of freedom. Its microscopic justication is provided by statistical mechanics. Equilibrium
More informationEfficient numerical integrators for stochastic models
Efficient numerical integrators for stochastic models G. De Fabritiis a, M. Serrano b P. Español b P.V.Coveney a a Centre for Computational Science, Department of Chemistry, University College London,
More informationSliding Friction in the Frenkel-Kontorova Model
arxiv:cond-mat/9510058v1 12 Oct 1995 to appear in "The Physics of Sliding Friction", Ed. B.N.J. Persson (Kluwer Academic Publishers), 1995 Sliding Friction in the Frenkel-Kontorova Model E. Granato I.N.P.E.,
More informationAn Extended van der Waals Equation of State Based on Molecular Dynamics Simulation
J. Comput. Chem. Jpn., Vol. 8, o. 3, pp. 97 14 (9) c 9 Society of Computer Chemistry, Japan An Extended van der Waals Equation of State Based on Molecular Dynamics Simulation Yosuke KATAOKA* and Yuri YAMADA
More informationThus, the volume element remains the same as required. With this transformation, the amiltonian becomes = p i m i + U(r 1 ; :::; r N ) = and the canon
G5.651: Statistical Mechanics Notes for Lecture 5 From the classical virial theorem I. TEMPERATURE AND PRESSURE ESTIMATORS hx i x j i = kt ij we arrived at the equipartition theorem: * + p i = m i NkT
More informationExplicit approximations to estimate the perturbative diusivity in the presence of convectivity and damping (Part 1): Semi-innite slab approximations
Explicit approximations to estimate the perturbative diusivity in the presence of convectivity and damping (Part 1): Semi-innite slab approximations M. van Berkel 1,2,3,4, H.J. Zwart 5,6, G.M.D. Hogeweij
More informationNPTEL
NPTEL Syllabus Nonequilibrium Statistical Mechanics - Video course COURSE OUTLINE Thermal fluctuations, Langevin dynamics, Brownian motion and diffusion, Fokker-Planck equations, linear response theory,
More informationThermodynamics of nuclei in thermal contact
Thermodynamics of nuclei in thermal contact Karl-Heinz Schmidt, Beatriz Jurado CENBG, CNRS/IN2P3, Chemin du Solarium B.P. 120, 33175 Gradignan, France Abstract: The behaviour of a di-nuclear system in
More informationBrownian Motion and The Atomic Theory
Brownian Motion and The Atomic Theory Albert Einstein Annus Mirabilis Centenary Lecture Simeon Hellerman Institute for Advanced Study, 5/20/2005 Founders Day 1 1. What phenomenon did Einstein explain?
More informationTemperature and Pressure Controls
Ensembles Temperature and Pressure Controls 1. (E, V, N) microcanonical (constant energy) 2. (T, V, N) canonical, constant volume 3. (T, P N) constant pressure 4. (T, V, µ) grand canonical #2, 3 or 4 are
More informationPrancing Through Quantum Fields
November 23, 2009 1 Introduction Disclaimer Review of Quantum Mechanics 2 Quantum Theory Of... Fields? Basic Philosophy 3 Field Quantization Classical Fields Field Quantization 4 Intuitive Field Theory
More informationStochastic equations for thermodynamics
J. Chem. Soc., Faraday Trans. 93 (1997) 1751-1753 [arxiv 1503.09171] Stochastic equations for thermodynamics Roumen Tsekov Department of Physical Chemistry, University of Sofia, 1164 Sofia, ulgaria The
More informationStatistical mechanics of classical systems
Statistical mechanics of classical systems States and ensembles A microstate of a statistical system is specied by the complete information about the states of all microscopic degrees of freedom of the
More informationDepartment of Mechanical Engineering ME 96. Free and Forced Convection Experiment. Revised: 25 April Introduction
CALIFORNIA INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering ME 96 Free and Forced Convection Experiment Revised: 25 April 1994 1. Introduction The term forced convection refers to heat transport
More informationCHAPTER 4. Basics of Fluid Dynamics
CHAPTER 4 Basics of Fluid Dynamics What is a fluid? A fluid is a substance that can flow, has no fixed shape, and offers little resistance to an external stress In a fluid the constituent particles (atoms,
More informationIntroduction Statistical Thermodynamics. Monday, January 6, 14
Introduction Statistical Thermodynamics 1 Molecular Simulations Molecular dynamics: solve equations of motion Monte Carlo: importance sampling r 1 r 2 r n MD MC r 1 r 2 2 r n 2 3 3 4 4 Questions How can
More informationIntroduction to Simulation - Lectures 17, 18. Molecular Dynamics. Nicolas Hadjiconstantinou
Introduction to Simulation - Lectures 17, 18 Molecular Dynamics Nicolas Hadjiconstantinou Molecular Dynamics Molecular dynamics is a technique for computing the equilibrium and non-equilibrium properties
More informationComputer simulations as concrete models for student reasoning
Computer simulations as concrete models for student reasoning Jan Tobochnik Department of Physics Kalamazoo College Kalamazoo MI 49006 In many thermal physics courses, students become preoccupied with
More informationU U B U P x
Smooth Quantum Hydrodynamic Model Simulation of the Resonant Tunneling Diode Carl L. Gardner and Christian Ringhofer y Department of Mathematics Arizona State University Tempe, AZ 8587-184 Abstract Smooth
More informationMonte Carlo methods for phase equilibria of uids 2 1. Introduction The phase behavior of uids and their mixtures is of central importance to many tech
TOPICAL REVIEW Monte Carlo methods for phase equilibria of uids Athanassios Z. Panagiotopoulos z Institute for Physical Science and Technology and Department of Chemical Engineering, University of Maryland,
More informationChaotic dynamics is one of the most extensively developing subjects in physics. Although
CHAOS vs THERMALIZATION IN THE NUCLEAR SHELL MODEL Mihai Horoi 1;3, Vladimir Zelevinsky 1;2;4 and B. Alex Brown 1;2 1 National Superconducting Cyclotron Laboratory, Michigan State University, East Lansing,
More informationNon equilibrium thermodynamics: foundations, scope, and extension to the meso scale. Miguel Rubi
Non equilibrium thermodynamics: foundations, scope, and extension to the meso scale Miguel Rubi References S.R. de Groot and P. Mazur, Non equilibrium Thermodynamics, Dover, New York, 1984 J.M. Vilar and
More informationLangevin Methods. Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 10 D Mainz Germany
Langevin Methods Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 1 D 55128 Mainz Germany Motivation Original idea: Fast and slow degrees of freedom Example: Brownian motion Replace
More informationto satisfy the large number approximations, W W sys can be small.
Chapter 12. The canonical ensemble To discuss systems at constant T, we need to embed them with a diathermal wall in a heat bath. Note that only the system and bath need to be large for W tot and W bath
More informationA Study of the Thermal Properties of a One. Dimensional Lennard-Jones System
A Study of the Thermal Properties of a One Dimensional Lennard-Jones System Abstract In this study, the behavior of a one dimensional (1D) Lennard-Jones (LJ) system is simulated. As part of this research,
More informationMigration Behavior of Bead-spring Dumbbell Models under Microchannel Flow from Dissipative Particle Dynamics Simulations
2426 Bull. Korean Chem. Soc. 2007, Vol. 28, No. 12 Kwang Jin Oh Migration Behavior of Bead-spring Dumbbell Models under Microchannel Flow from Dissipative Particle Dynamics Simulations Kwang Jin Oh Supercomputing
More information(a). t = 0 (b). t = 0.276
10 (a). t = 0 (b). t = 0.276 Figure 8: Snapshots of the sedimentation of 5040 disks of diameter 10/192 cm in 2D (W =8 cm). The initial lattice is disordered. from t = 0:026 to t = 0:476. In Figure 6, a
More informationII Relationship of Classical Theory to Quantum Theory A Quantum mean occupation number
Appendix B Some Unifying Concepts Version 04.AppB.11.1K [including mostly Chapters 1 through 11] by Kip [This appendix is in the very early stages of development] I Physics as Geometry A Newtonian Physics
More informationNon-equilibrium phenomena and fluctuation relations
Non-equilibrium phenomena and fluctuation relations Lamberto Rondoni Politecnico di Torino Beijing 16 March 2012 http://www.rarenoise.lnl.infn.it/ Outline 1 Background: Local Thermodyamic Equilibrium 2
More informationHydrodynamics. Stefan Flörchinger (Heidelberg) Heidelberg, 3 May 2010
Hydrodynamics Stefan Flörchinger (Heidelberg) Heidelberg, 3 May 2010 What is Hydrodynamics? Describes the evolution of physical systems (classical or quantum particles, fluids or fields) close to thermal
More informationChE 210B: Advanced Topics in Equilibrium Statistical Mechanics
ChE 210B: Advanced Topics in Equilibrium Statistical Mechanics Glenn Fredrickson Lecture 1 Reading: 3.1-3.5 Chandler, Chapters 1 and 2 McQuarrie This course builds on the elementary concepts of statistical
More informationDissipative Particle Dynamics
CHAPTER 2 Dissipative Particle Dynamics Igor V. Pivkin, a,b Bruce Caswell, a and George Em Karniadakis a a Division of Applied Mathematics, Brown University, Providence, Rhode Island b Department of Materials
More informationEntropy, extropy and information potential in stochastic systems far from equilibrium
Physica A 305 (2002) 445 466 www.elsevier.com/locate/physa Entropy, extropy and information potential in stochastic systems far from equilibrium B. Gaveau a;, K. Martinas b, M. Moreau c,j.toth d;1 a UFR
More informationn ~ k d ~ k =2 ~ k 2 d ~ k 3 (2) 2 ( (2) ~ k + ~ k 6 ~ k2 ~ k3 )2(k 2 + k 2 k 2 2 k 2 ) 3 n ~k n ~k n ~k2 n ~k3 n ~k + n ~k n ~k2 n ~k3! : (2
Kolmogorov and Nonstationary Spectra of Optical Turbulence Vladimir M. Malkin Budker Institute of Nuclear Physics Novosibirsk, 630090, Russia () Abstract Optical turbulence, which develops in the transverse
More informationTemperature Fluctuations and Entropy Formulas
Temperature Fluctuations and Entropy Formulas Does it or does not? T.S. Biró G.G. Barnaföldi P. Ván MTA Heavy Ion Research Group Research Centre for Physics, Budapest February 1, 2014 J.Uffink, J.van Lith:
More informationExamples of the Dirac approach to dynamics of systems with constraints
Physica A 290 (2001) 431 444 www.elsevier.com/locate/physa Examples of the Dirac approach to dynamics of systems with constraints Sonnet Nguyen, Lukasz A. Turski Center for Theoretical Physics, College
More information7 To solve numerically the equation of motion, we use the velocity Verlet or leap frog algorithm. _ V i n = F i n m i (F.5) For time step, we approxim
69 Appendix F Molecular Dynamics F. Introduction In this chapter, we deal with the theories and techniques used in molecular dynamics simulation. The fundamental dynamics equations of any system is the
More information18.13 Review & Summary
5/2/10 10:04 PM Print this page 18.13 Review & Summary Temperature; Thermometers Temperature is an SI base quantity related to our sense of hot and cold. It is measured with a thermometer, which contains
More informationStatistical Thermodynamics and Monte-Carlo Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany
Statistical Thermodynamics and Monte-Carlo Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals From Micro to Macro Statistical Mechanics (Statistical
More informationUniqueness of thermodynamic projector and kinetic basis of molecular individualism
Available online at www.sciencedirect.com Physica A 336 (2004) 391 432 www.elsevier.com/locate/physa Uniqueness of thermodynamic projector and kinetic basis of molecular individualism Alexander N. Gorban
More informationQuantum measurement theory and micro-macro consistency in nonequilibrium statistical mechanics
Nagoya Winter Workshop on Quantum Information, Measurement, and Quantum Foundations (Nagoya, February 18-23, 2010) Quantum measurement theory and micro-macro consistency in nonequilibrium statistical mechanics
More informationTurbulent penetrative convection with an internal heat source
Fluid Dynamics Research 8 (001) 397 41 Turbulent penetrative convection with an internal heat source J.R. Chasnov, K.L. Tse 1 Department of Mathematics, The Hong Kong University of Science and Technology,
More informationMultistage pulse tubes
Cryogenics 40 (2000) 459±464 www.elsevier.com/locate/cryogenics Multistage pulse tubes A.T.A.M. de Waele *, I.A. Tanaeva, Y.L. Ju Department of Physics, Eindhoven University of Technology, P.O. Box 513,
More informationKinetic Theory. Motivation - Relaxation Processes Violent Relaxation Thermodynamics of self-gravitating system
Kinetic Theory Motivation - Relaxation Processes Violent Relaxation Thermodynamics of self-gravitating system negative heat capacity the gravothermal catastrophe The Fokker-Planck approximation Master
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 11 Dec 2002
arxiv:cond-mat/0212257v1 [cond-mat.stat-mech] 11 Dec 2002 International Journal of Modern Physics B c World Scientific Publishing Company VISCOUS FINGERING IN MISCIBLE, IMMISCIBLE AND REACTIVE FLUIDS PATRICK
More informationDerivation of Van der Waal s equation of state in microcanonical ensemble formulation
arxiv:180.01963v1 [physics.gen-ph] 9 Nov 017 Derivation of an der Waal s equation of state in microcanonical ensemble formulation Aravind P. Babu, Kiran S. Kumar and M. Ponmurugan* Department of Physics,
More informationAdaptive resolution molecular-dynamics simulation: Changing the degrees of freedom on the fly
THE JOURNAL OF CHEMICAL PHYSICS 123, 224106 2005 Adaptive resolution molecular-dynamics simulation: Changing the degrees of freedom on the fly Matej Praprotnik, a Luigi Delle Site, and Kurt Kremer b Max-Planck-Institut
More informationHeat and Mass Transfer over Cooled Horizontal Tubes 333 x-component of the velocity: y2 u = g sin x y : (4) r 2 The y-component of the velocity eld is
Scientia Iranica, Vol., No. 4, pp 332{338 c Sharif University of Technology, October 2004 Vapor Absorption into Liquid Films Flowing over a Column of Cooled Horizontal Tubes B. Farhanieh and F. Babadi
More informationJavier Junquera. Statistical mechanics
Javier Junquera Statistical mechanics From the microscopic to the macroscopic level: the realm of statistical mechanics Computer simulations Thermodynamic state Generates information at the microscopic
More information