Stochastic Particle Dynamics for Unresolved Degrees of Freedom. Sebastian Reich in collaborative work with Colin Cotter (Imperial College)
|
|
- Vivien Austin
- 6 years ago
- Views:
Transcription
1 Stochastic Particle Dynamics for Unresolved Degrees of Freedom Sebastian Reich in collaborative work with Colin Cotter (Imperial College)
2 1. The Motivation Classical molecular dynamics (MD) as well as particle methods for numerical weather prediction (NWP) lead to large systems of conservative Newtonian equations of motion.
3 1. The Motivation Classical molecular dynamics (MD) as well as particle methods for numerical weather prediction (NWP) lead to large systems of conservative Newtonian equations of motion. Underresolved simulations may occur for two reasons:
4 1. The Motivation Classical molecular dynamics (MD) as well as particle methods for numerical weather prediction (NWP) lead to large systems of conservative Newtonian equations of motion. Underresolved simulations may occur for two reasons: certain degrees of freedom are ignored and absent in the model,
5 1. The Motivation Classical molecular dynamics (MD) as well as particle methods for numerical weather prediction (NWP) lead to large systems of conservative Newtonian equations of motion. Underresolved simulations may occur for two reasons: certain degrees of freedom are ignored and absent in the model, the equations of motion are integrated with time steps large compared to the fastest time-scales present in the system.
6 1. The Motivation Classical molecular dynamics (MD) as well as particle methods for numerical weather prediction (NWP) lead to large systems of conservative Newtonian equations of motion. Underresolved simulations may occur for two reasons: certain degrees of freedom are ignored and absent in the model, the equations of motion are integrated with time steps large compared to the fastest time-scales present in the system. One might want to model the effect of unresolved motion using stochastic models.
7 2. The Plan We start with the stiff spring pendulum and compare
8 2. The Plan We start with the stiff spring pendulum and compare the conservative model,
9 2. The Plan We start with the stiff spring pendulum and compare the conservative model, the conservative model plus Langevin dynamics,
10 2. The Plan We start with the stiff spring pendulum and compare the conservative model, the conservative model plus Langevin dynamics, the dissipative particle dynamics (DPD) model,
11 2. The Plan We start with the stiff spring pendulum and compare the conservative model, the conservative model plus Langevin dynamics, the dissipative particle dynamics (DPD) model, the extended DPD model.
12 2. The Plan We start with the stiff spring pendulum and compare the conservative model, the conservative model plus Langevin dynamics, the dissipative particle dynamics (DPD) model, the extended DPD model. We then discuss these algorithms in the context of large time step methods for MD, unresolved wave motion in NWP.
13 3. The Stiff Spring Pendulum We use the stiff spring pendulum equations ṙ = p, ṗ = g K r ( r L) r as our model problem. Here r, p R 3, g = (0, 0, 1) T, K = 1000 is the spring constant, and L = 1 is the length of the pendulum.
14 3. The Stiff Spring Pendulum We use the stiff spring pendulum equations ṙ = p, ṗ = g K r ( r L) r as our model problem. Here r, p R 3, g = (0, 0, 1) T, K = 1000 is the spring constant, and L = 1 is the length of the pendulum. The stiff pendulum equations are Hamiltonian with conserved energy where we introduced r = r. E = 1 2 pt p + K 2 (r L)2 r T g,
15 Furthermore, because of the relatively large spring constant K = 1000, the (oscillatory) energy in the spring is an adiabatic invariant, i.e., J = 1 2 p2 r + K 2 (r L)2 with p r = r T p/r, is approximately conserved over long intervals of time.
16 Furthermore, because of the relatively large spring constant K = 1000, the (oscillatory) energy in the spring is an adiabatic invariant, i.e., J = 1 2 p2 r + K 2 (r L)2 with p r = r T p/r, is approximately conserved over long intervals of time. In the following figure we plot E(t), J(t), and the slow rotational energy E s (t) = E(t) J(t) from a simulation using the symplectic Störmer-Verlet method.
17
18 4. The Langevin Model for the Stiff Spring Pendulum We next replace the Hamiltonian equations of motion by the Langevin model ṙ = p, ṗ = g K r r (r L) γp + γẇ(t), with γ = 0.1, σ = 0.2, k B T = 1, and W(t) a vector of independent Wiener processes.
19 4. The Langevin Model for the Stiff Spring Pendulum We next replace the Hamiltonian equations of motion by the Langevin model ṙ = p, ṗ = g K r (r L) γp + γẇ(t), r with γ = 0.1, σ = 0.2, k B T = 1, and W(t) a vector of independent Wiener processes. It is apparent from the following figure that none of the energy contributions is conserved (even approximately) any longer. Note that, by definition, J 0, E r 1, and E s r 1.
20
21 5. Dissipative Particle Dynamics (DPD) The DPD model (Hoogerbrugge & Koelman, Espanol & Warren) leads to a stochastic coupling only in the direction of the stiff spring and the equations of motion become ṙ = p, ṗ = g r [ K(r L)+γpr + γẇ (t) ], r with the parameters values as before and W (t) a scalar Wiener process.
22 5. Dissipative Particle Dynamics (DPD) The DPD model (Hoogerbrugge & Koelman, Espanol & Warren) leads to a stochastic coupling only in the direction of the stiff spring and the equations of motion become ṙ = p, ṗ = g r [ K(r L)+γpr + γẇ (t) ], r with the parameters values as before and W (t) a scalar Wiener process. Indeed, the simulation results, as displayed in the following figure, reveal that the slow energy is now approximately conserved while both the total and the oscillatory energy drift.
23 5. Dissipative Particle Dynamics (DPD) The DPD model (Hoogerbrugge & Koelman, Espanol & Warren) leads to a stochastic coupling only in the direction of the stiff spring and the equations of motion become ṙ = p, ṗ = g r [ K(r L)+γpr + γẇ (t) ], r with the parameters values as before and W (t) a scalar Wiener process. Indeed, the simulation results, as displayed in the following figure, reveal that the slow energy is now approximately conserved while both the total and the oscillatory energy drift. This demonstrates the capability of the DPD model to thermostate degrees of freedom in a rather targeted manner.
24
25 6. The Extended DPD Model (Cotter & Reich, 2003) Finally, we replace the DPD model by the extended version ṙ = p, ṗ = g r [K(r L) + s], r εṡ = s + γp r + σẇ (t) for various values of the correlation parameter ε. The parameters γ and σ are left as before and W (t) is again a Wiener process.
26 6. The Extended DPD Model (Cotter & Reich, 2003) Finally, we replace the DPD model by the extended version ṙ = p, ṗ = g r [K(r L) + s], r εṡ = s + γp r + σẇ (t) for various values of the correlation parameter ε. The parameters γ and σ are left as before and W (t) is again a Wiener process. The fluctuation-dissipation variable s satisfies an Ornstein- Uhlenbeck process, which has exponential decay of correlation K(τ) ε 1 e τ /ε.
27 6. The Extended DPD Model (Cotter & Reich, 2003) Finally, we replace the DPD model by the extended version ṙ = p, ṗ = g r [K(r L) + s], r εṡ = s + γp r + σẇ (t) for various values of the correlation parameter ε. The parameters γ and σ are left as before and W (t) is again a Wiener process. The fluctuation-dissipation variable s satisfies an Ornstein- Uhlenbeck process, which has exponential decay of correlation K(τ) ε 1 e τ /ε. The extended model reduces to the standard DPD model for ε 0.
28 The rapid change in energy behavior shown in the following figure for parameter values from an interval 1 ε 0.01 is quite striking and indicates that the model s response to added fluctuationdissipation can depend crucially on the decay rate ε 1 of the associated auto-covariance matrix.
29
30 7. Conservative Large Time-Step Methods Consider a spiced-up stiff spring pendulum with Hamiltonian H = 1 2 pt p + K 2 (r L)2 + V short (r) + V long (r) where V short is some short-ranged potential energy function while V long is long-ranged and expensive to compute.
31 7. Conservative Large Time-Step Methods Consider a spiced-up stiff spring pendulum with Hamiltonian H = 1 2 pt p + K 2 (r L)2 + V short (r) + V long (r) where V short is some short-ranged potential energy function while V long is long-ranged and expensive to compute. Then symplectic multiple-time-stepping is a good idea but prone to numerical resonance instabilities. It was suggested by Garcia- Archilla, Sanz-Serna & Skeel to mollify the long range potential energy V long to avoid these instabilities.
32 A simple mollifier is Equilibrium where H = 1 2 pt p + K 2 (r L)2 + V short (r) + V long ( r) and r(r) is defined by the equations r = (1 + λ)r, r = L.
33 A simple mollifier is Equilibrium where H = 1 2 pt p + K 2 (r L)2 + V short (r) + V long ( r) and r(r) is defined by the equations r = (1 + λ)r, r = L. It has been shown by Izaguirre, Reich & Skeel that Equilibrium allows step-size of up to 7 fs for flexible water molecules. The energy drift observed for step-sizes between 7 fs and 15 fs is rather mild. This suggests to...
34 8. Stochastic Large Time-Step Methods... add a Langevin term to the equations of motion (see Izaguirre, Catarello, Wozniak & Skeel). But that approach couples noise to all degrees of freedom and in all directions.
35 8. Stochastic Large Time-Step Methods... add a Langevin term to the equations of motion (see Izaguirre, Catarello, Wozniak & Skeel). But that approach couples noise to all degrees of freedom and in all directions. More recently Izaguirre & Ma suggested Targeted MOLLY (TM), which combines mollification of the long range potential with a DPD-type fluctuation-dissipation on the highly oscillatory degrees of freedom.
36 8. Stochastic Large Time-Step Methods... add a Langevin term to the equations of motion (see Izaguirre, Catarello, Wozniak & Skeel). But that approach couples noise to all degrees of freedom and in all directions. More recently Izaguirre & Ma suggested Targeted MOLLY (TM), which combines mollification of the long range potential with a DPD-type fluctuation-dissipation on the highly oscillatory degrees of freedom. For flexible water this approach allow for step-sizes up to fs without affecting diffusion constants etc.
37 9. Particle Methods for NWP The rotating two-dimensional shallow-water equations (SWEs) describe a shallow layer of fluid subject to gravity and rotation.
38 9. Particle Methods for NWP The rotating two-dimensional shallow-water equations (SWEs) describe a shallow layer of fluid subject to gravity and rotation. The Lagrangian formulation of the SWEs is ẍ = fẋ c 2 0 xh(x),
39 9. Particle Methods for NWP The rotating two-dimensional shallow-water equations (SWEs) describe a shallow layer of fluid subject to gravity and rotation. The Lagrangian formulation of the SWEs is ẍ = fẋ c 2 0 xh(x), where x = (x, y) T are the particle positions (a continuous function of both space and time),
40 9. Particle Methods for NWP The rotating two-dimensional shallow-water equations (SWEs) describe a shallow layer of fluid subject to gravity and rotation. The Lagrangian formulation of the SWEs is ẍ = fẋ c 2 0 xh(x), where x = (x, y) T are the particle positions (a continuous function of both space and time),c 0 = gh, g is the gravitational constant, H is the mean layer-depth,
41 9. Particle Methods for NWP The rotating two-dimensional shallow-water equations (SWEs) describe a shallow layer of fluid subject to gravity and rotation. The Lagrangian formulation of the SWEs is ẍ = fẋ c 2 0 xh(x), where x = (x, y) T are the particle positions (a continuous function of both space and time),c 0 = gh, g is the gravitational constant, H is the mean layer-depth,f is twice the (constant) angular velocity of the fluid, and ẋ = ( ẏ, ẋ) T.
42 Frank, Gottwald & Reich suggested a Hamiltonian particle-mesh (HPM) method for the solution of the two-dimensional shallowwater equations.
43 Frank, Gottwald & Reich suggested a Hamiltonian particle-mesh (HPM) method for the solution of the two-dimensional shallowwater equations. The HPM method may be viewed as an accurate numerical discretisation of the regularised fluid equations: ẍ = fẋ c 2 0 x [A h(x)] where A is a smoothing operator with some smoothing length α.
44 Frank, Gottwald & Reich suggested a Hamiltonian particle-mesh (HPM) method for the solution of the two-dimensional shallowwater equations. The HPM method may be viewed as an accurate numerical discretisation of the regularised fluid equations: ẍ = fẋ c 2 0 x [A h(x)] where A is a smoothing operator with some smoothing length α. We denote the numerically unresolved part of the layer-depth by η = h A h.
45 The HPM method uses a regular grid x kl = (k x, l y) T, moving particles x i = (x i, y i ) T, and grid-centred basis functions ψ kl (X).
46 The HPM method uses a regular grid x kl = (k x, l y) T, moving particles x i = (x i, y i ) T, and grid-centred basis functions ψ kl (X). The finite-dimensional Hamiltonian equations of motion are given by ẍ i = fẋ i c 2 0 k,l xi ψ kl (x i )ĥ kl.
47 The HPM method uses a regular grid x kl = (k x, l y) T, moving particles x i = (x i, y i ) T, and grid-centred basis functions ψ kl (X). The finite-dimensional Hamiltonian equations of motion are given by ẍ i = fẋ i c 2 0 k,l xi ψ kl (x i )ĥ kl. The numerically unresolved gravity waves η can be modelled by a generalized Langevin process.
48 The HPM method uses a regular grid x kl = (k x, l y) T, moving particles x i = (x i, y i ) T, and grid-centred basis functions ψ kl (X). The finite-dimensional Hamiltonian equations of motion are given by ẍ i = fẋ i c 2 0 k,l xi ψ kl (x i )ĥ kl. The numerically unresolved gravity waves η can be modelled by a generalized Langevin process. This idea can be mathematically motivated by representing η as the solution of a linear wave equation coupled to the particle system and subsequent reduction following the Kac-Zwanzig approach.
49 This leads to the extended DPD equations given by ẍ i = fẋ i c 2 ] 0 [ĥkl + η kl Xi ψ kl (X i ), k,l ε η kl = η kl + γḣ kl + σẇ kl.
50 This leads to the extended DPD equations given by ẍ i = fẋ i c 2 ] 0 [ĥkl + η kl Xi ψ kl (X i ), k,l ε η kl = η kl + γḣ kl + σẇ kl. It should be noted that h kl = ĥ kl + η kl is ment to be a better approximation to the true layer-depth h at the grid point x kl than ĥ kl alone.
51 This leads to the extended DPD equations given by ẍ i = fẋ i c 2 ] 0 [ĥkl + η kl Xi ψ kl (X i ), k,l ε η kl = η kl + γḣ kl + σẇ kl. It should be noted that h kl = ĥ kl + η kl is ment to be a better approximation to the true layer-depth h at the grid point x kl than ĥ kl alone. The unresolved waves η kl are modeled by Ornstein-Uhlenbeck processes. Whether this is reasonable in terms of real fluid motion remains an open question. But it can be tested numerically (time consuming...).
52 10. The HPM Method on the Sphere Lately we had great fun implementating the particle method for spherical geometry. A common test case is provided by a five day forecast for 21st December 1978 (500 mb geopotential height field).
53
54
55
56 11. Conclusion
57 11. Conclusion We started by looking at a simple stiff spring pendulum and the effect of various stochastic modifications on the dynamics.
58 11. Conclusion We started by looking at a simple stiff spring pendulum and the effect of various stochastic modifications on the dynamics. DPD and its extended version appear as more sensitive than straightforward Langevin dynamics which is still widely used in MD.
59 11. Conclusion We started by looking at a simple stiff spring pendulum and the effect of various stochastic modifications on the dynamics. DPD and its extended version appear as more sensitive than straightforward Langevin dynamics which is still widely used in MD. Two particular examples for which this is certainly true are (i) stablized long time step methods and (ii) particle methods for fluid dynamics.
60 11. Conclusion We started by looking at a simple stiff spring pendulum and the effect of various stochastic modifications on the dynamics. DPD and its extended version appear as more sensitive than straightforward Langevin dynamics which is still widely used in MD. Two particular examples for which this is certainly true are (i) stablized long time step methods and (ii) particle methods for fluid dynamics. The material presented in this talk is part of a survey article written by us for a handbook on nanoscale computation techniques (American Scientific Publisher).
61 11. Conclusion We started by looking at a simple stiff spring pendulum and the effect of various stochastic modifications on the dynamics. DPD and its extended version appear as more sensitive than straightforward Langevin dynamics which is still widely used in MD. Two particular examples for which this is certainly true are (i) stablized long time step methods and (ii) particle methods for fluid dynamics. The material presented in this talk is part of a survey article written by us for a handbook on nanoscale computation techniques (American Scientific Publisher). THE END
A Hamiltonian Numerical Scheme for Large Scale Geophysical Fluid Systems
A Hamiltonian Numerical Scheme for Large Scale Geophysical Fluid Systems Bob Peeters Joint work with Onno Bokhove & Jason Frank TW, University of Twente, Enschede CWI, Amsterdam PhD-TW colloquium, 9th
More informationLarge-scale atmospheric circulation, semi-geostrophic motion and Lagrangian particle methods
Large-scale atmospheric circulation, semi-geostrophic motion and Lagrangian particle methods Colin Cotter (Imperial College London) & Sebastian Reich (Universität Potsdam) Outline 1. Hydrostatic and semi-geostrophic
More informationThe Hamiltonian Particle-Mesh Method for the Spherical Shallow Water Equations
The Hamiltonian Particle-Mesh Method for the Spherical Shallow Water Equations Jason Fran CWI P.O. Box 94079, 1090 GB Amsterdam, The Netherlands Sebastian Reich Department of Mathematics, Imperial College
More informationHamiltonian particle-mesh simulations for a non-hydrostatic vertical slice model
Hamiltonian particle-mesh simulations for a non-hydrostatic vertical slice model Seoleun Shin Sebastian Reich May 6, 29 Abstract A Lagrangian particle method is developed for the simulation of atmospheric
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 informationM2A2 Problem Sheet 3 - Hamiltonian Mechanics
MA Problem Sheet 3 - Hamiltonian Mechanics. The particle in a cone. A particle slides under gravity, inside a smooth circular cone with a vertical axis, z = k x + y. Write down its Lagrangian in a) Cartesian,
More informationMolecular Dynamics Lecture 3
Molecular Dynamics Lecture 3 Ben Leimkuhler the problem of the timestep in MD constraints - SHAKE/RATTLE multiple timestepping stochastic methods for holonomic constraints stochastic multiple timestepping
More informationThe Two Ball Newton s Cradle. Glendinning, Paul. MIMS EPrint: Manchester Institute for Mathematical Sciences School of Mathematics
The Two Ball Newton s Cradle Glendinning, Paul 11 MIMS EPrint: 11.65 Manchester Institute for Mathematical Sciences School of Mathematics The University of Manchester Reports available from: And by contacting:
More informationLANGEVIN THEORY OF BROWNIAN MOTION. Contents. 1 Langevin theory. 1 Langevin theory 1. 2 The Ornstein-Uhlenbeck process 8
Contents LANGEVIN THEORY OF BROWNIAN MOTION 1 Langevin theory 1 2 The Ornstein-Uhlenbeck process 8 1 Langevin theory Einstein (as well as Smoluchowski) was well aware that the theory of Brownian motion
More informationAPMA 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 informationA mollified ensemble Kalman filter
A mollified ensemble Kalman filter Kay Bergemann Sebastian Reich November 2, 2018 arxiv:1002.3091v2 [math.na] 22 Jun 2010 Abstract It is well recognized that discontinuous analysis increments of sequential
More informationThe Truth about diffusion (in liquids)
The Truth about diffusion (in liquids) Aleksandar Donev Courant Institute, New York University & Eric Vanden-Eijnden, Courant In honor of Berni Julian Alder LLNL, August 20th 2015 A. Donev (CIMS) Diffusion
More informationPhysical Dynamics (SPA5304) Lecture Plan 2018
Physical Dynamics (SPA5304) Lecture Plan 2018 The numbers on the left margin are approximate lecture numbers. Items in gray are not covered this year 1 Advanced Review of Newtonian Mechanics 1.1 One Particle
More informationA mollified ensemble Kalman filter
Quarterly Journal of the Royal Meteorological Society Q. J. R. Meteorol. Soc. 136: 1636 1643, July 1 Part B A mollified ensemble Kalman filter Kay Bergemann and Sebastian Reich* Universität Potsdam, Potsdam,
More informationHypoelliptic multiscale Langevin diffusions and Slow fast stochastic reaction diffusion equations.
Hypoelliptic multiscale Langevin diffusions and Slow fast stochastic reaction diffusion equations. Wenqing Hu. 1 (Joint works with Michael Salins 2 and Konstantinos Spiliopoulos 3.) 1. Department of Mathematics
More informationPhysical Dynamics (PHY-304)
Physical Dynamics (PHY-304) Gabriele Travaglini March 31, 2012 1 Review of Newtonian Mechanics 1.1 One particle Lectures 1-2. Frame, velocity, acceleration, number of degrees of freedom, generalised coordinates.
More informationHamiltonian. March 30, 2013
Hamiltonian March 3, 213 Contents 1 Variational problem as a constrained problem 1 1.1 Differential constaint......................... 1 1.2 Canonic form............................. 2 1.3 Hamiltonian..............................
More informationData assimilation as an optimal control problem and applications to UQ
Data assimilation as an optimal control problem and applications to UQ Walter Acevedo, Angwenyi David, Jana de Wiljes & Sebastian Reich Universität Potsdam/ University of Reading IPAM, November 13th 2017
More informationGeometric Numerical Integration
Geometric Numerical Integration (Ernst Hairer, TU München, winter 2009/10) Development of numerical ordinary differential equations Nonstiff differential equations (since about 1850), see [4, 2, 1] Adams
More informationLecture 1: Pragmatic Introduction to Stochastic Differential Equations
Lecture 1: Pragmatic Introduction to Stochastic Differential Equations Simo Särkkä Aalto University, Finland (visiting at Oxford University, UK) November 13, 2013 Simo Särkkä (Aalto) Lecture 1: Pragmatic
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 informationEstimating Accuracy in Classical Molecular Simulation
Estimating Accuracy in Classical Molecular Simulation University of Illinois Urbana-Champaign Department of Computer Science Institute for Mathematics and its Applications July 2007 Acknowledgements Ben
More informationFinite element exterior calculus framework for geophysical fluid dynamics
Finite element exterior calculus framework for geophysical fluid dynamics Colin Cotter Department of Aeronautics Imperial College London Part of ongoing work on UK Gung-Ho Dynamical Core Project funded
More informationCP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS. Prof. N. Harnew University of Oxford TT 2017
CP1 REVISION LECTURE 3 INTRODUCTION TO CLASSICAL MECHANICS Prof. N. Harnew University of Oxford TT 2017 1 OUTLINE : CP1 REVISION LECTURE 3 : INTRODUCTION TO CLASSICAL MECHANICS 1. Angular velocity and
More informationOnsager theory: overview
Onsager theory: overview Pearu Peterson December 18, 2006 1 Introduction Our aim is to study matter that consists of large number of molecules. A complete mechanical description of such a system is practically
More informationModelling Interactions Between Weather and Climate
Modelling Interactions Between Weather and Climate p. 1/33 Modelling Interactions Between Weather and Climate Adam Monahan & Joel Culina monahana@uvic.ca & culinaj@uvic.ca School of Earth and Ocean Sciences,
More informationNumerical Methods for ODEs. Lectures for PSU Summer Programs Xiantao Li
Numerical Methods for ODEs Lectures for PSU Summer Programs Xiantao Li Outline Introduction Some Challenges Numerical methods for ODEs Stiff ODEs Accuracy Constrained dynamics Stability Coarse-graining
More informationJ10M.1 - Rod on a Rail (M93M.2)
Part I - Mechanics J10M.1 - Rod on a Rail (M93M.2) J10M.1 - Rod on a Rail (M93M.2) s α l θ g z x A uniform rod of length l and mass m moves in the x-z plane. One end of the rod is suspended from a straight
More informationM3-4-5 A16 Notes for Geometric Mechanics: Oct Nov 2011
M3-4-5 A16 Notes for Geometric Mechanics: Oct Nov 2011 Text for the course: Professor Darryl D Holm 25 October 2011 Imperial College London d.holm@ic.ac.uk http://www.ma.ic.ac.uk/~dholm/ Geometric Mechanics
More informationthe EL equation for the x coordinate is easily seen to be (exercise)
Physics 6010, Fall 2016 Relevant Sections in Text: 1.3 1.6 Examples After all this formalism it is a good idea to spend some time developing a number of illustrative examples. These examples represent
More informationClassical Equations of Motion
3 Classical Equations of Motion Several formulations are in use Newtonian Lagrangian Hamiltonian Advantages of non-newtonian formulations more general, no need for fictitious forces better suited for multiparticle
More informationG : Statistical Mechanics
G25.2651: Statistical Mechanics Notes for Lecture 15 Consider Hamilton s equations in the form I. CLASSICAL LINEAR RESPONSE THEORY q i = H p i ṗ i = H q i We noted early in the course that an ensemble
More informationMolecular Dynamics Lecture 2
Molecular Dynamics Lecture 2 Ben Leimkuhler the problem of the timestep in md splitting methods in Langevin dynamics constraints multiple timestepping Spring School, WIAS, Berlin 2017 fruit fly alanine
More informationLegendre Transforms, Calculus of Varations, and Mechanics Principles
page 437 Appendix C Legendre Transforms, Calculus of Varations, and Mechanics Principles C.1 Legendre Transforms Legendre transforms map functions in a vector space to functions in the dual space. From
More informationGeometric Mechanics and Global Nonlinear Control for Multi-Body Dynamics
Geometric Mechanics and Global Nonlinear Control for Multi-Body Dynamics Harris McClamroch Aerospace Engineering, University of Michigan Joint work with Taeyoung Lee (George Washington University) Melvin
More informationarxiv: v1 [cond-mat.stat-mech] 8 Jan 2019
Investigation of dissipative particle dynamics with colored noise arxiv:1901.02451v1 [cond-mat.stat-mech] 8 Jan 2019 Morgane Borreguero (a), Marco Ellero (b), and N.A. Adams (a) (a) Chair of Aerodynamics
More informationEvaluation of three spatial discretization schemes with the Galewsky et al. test
Evaluation of three spatial discretization schemes with the Galewsky et al. test Seoleun Shin Matthias Sommer Sebastian Reich Peter Névir February 22, 2 Abstract We evaluate the Hamiltonian Particle Methods
More informationFrom Molecular Dynamics to hydrodynamics a novel Galilean invariant thermostat
To appear in: Journal of Chemical Physics, 08 APR 005 From Molecular Dynamics to hydrodynamics a novel Galilean invariant thermostat Simeon D. Stoyanov and Robert D. Groot Unilever Research Vlaardingen,
More informationStochastic Stokes drift of a flexible dumbbell
Stochastic Stokes drift of a flexible dumbbell Kalvis M. Jansons arxiv:math/0607797v4 [math.pr] 22 Mar 2007 Department of Mathematics, University College London, Gower Street, London WC1E 6BT, UK January
More informationExponential integrators for oscillatory second-order differential equations
Exponential integrators for oscillatory second-order differential equations Marlis Hochbruck and Volker Grimm Mathematisches Institut Heinrich Heine Universität Düsseldorf Cambridge, March 2007 Outline
More informationVIII.B Equilibrium Dynamics of a Field
VIII.B Equilibrium Dynamics of a Field The next step is to generalize the Langevin formalism to a collection of degrees of freedom, most conveniently described by a continuous field. Let us consider the
More informationExamination paper for TMA4195 Mathematical Modeling
Department of Mathematical Sciences Examination paper for TMA4195 Mathematical Modeling Academic contact during examination: Elena Celledoni Phone: 48238584, 73593541 Examination date: 11th of December
More informationSemi-implicit methods, nonlinear balance, and regularized equations
ATMOSPHERIC SCIENCE LETTERS Atmos. Sci. Let. 8: 1 6 (7 Published online 9 January 7 in Wiley InterScience (www.interscience.wiley.com.1 Semi-implicit methods, nonlinear balance, and regularized equations
More informationPhysics 6010, Fall 2016 Constraints and Lagrange Multipliers. Relevant Sections in Text:
Physics 6010, Fall 2016 Constraints and Lagrange Multipliers. Relevant Sections in Text: 1.3 1.6 Constraints Often times we consider dynamical systems which are defined using some kind of restrictions
More informationF r (t) = 0, (4.2) F (t) = r= r a (t)
Chapter 4 Stochastic Equations 4.1 Langevin equation To explain the idea of introduction of the Langevin equation, let us consider the concrete example taken from surface physics: motion of an atom (adatom)
More informationAdaptive timestepping for SDEs with non-globally Lipschitz drift
Adaptive timestepping for SDEs with non-globally Lipschitz drift Mike Giles Wei Fang Mathematical Institute, University of Oxford Workshop on Numerical Schemes for SDEs and SPDEs Université Lille June
More informationDifferent approaches to model wind speed based on stochastic differential equations
1 Universidad de Castilla La Mancha Different approaches to model wind speed based on stochastic differential equations Rafael Zárate-Miñano Escuela de Ingeniería Minera e Industrial de Almadén Universidad
More informationTime Stepping Algorithms for Classical Molecular Dynamics
Time Stepping Algorithms for Classical Molecular Dynamics Colin John Cotter an Sebastian Reich Department of Mathematics, Imperial College 180 Queen s Gate, Lonon, SW7 2AZ (Date: March 27, 2004) Electronic
More information4.1 Important Notes on Notation
Chapter 4. Lagrangian Dynamics (Most of the material presented in this chapter is taken from Thornton and Marion, Chap. 7) 4.1 Important Notes on Notation In this chapter, unless otherwise stated, the
More informationGaussian processes for inference in stochastic differential equations
Gaussian processes for inference in stochastic differential equations Manfred Opper, AI group, TU Berlin November 6, 2017 Manfred Opper, AI group, TU Berlin (TU Berlin) inference in SDE November 6, 2017
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 informationAlgorithms for Ensemble Control. B Leimkuhler University of Edinburgh
Algorithms for Ensemble Control B Leimkuhler University of Edinburgh Model Reduction Across Disciplines Leicester 2014 (Stochastic) Ensemble Control extended ODE/SDE system design target ensemble SDE invariant
More informationThe Shallow Water Equations
If you have not already done so, you are strongly encouraged to read the companion file on the non-divergent barotropic vorticity equation, before proceeding to this shallow water case. We do not repeat
More informationBACKGROUND IN SYMPLECTIC GEOMETRY
BACKGROUND IN SYMPLECTIC GEOMETRY NILAY KUMAR Today I want to introduce some of the symplectic structure underlying classical mechanics. The key idea is actually quite old and in its various formulations
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 informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationLecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3.12 in Boas)
Lecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3 in Boas) As suggested in Lecture 8 the formalism of eigenvalues/eigenvectors has many applications in physics, especially in
More informationLecture 4: High oscillations
Geometric Numerical Integration TU München Ernst Hairer January February Lecture 4: High oscillations Table of contents A Fermi Pasta Ulam type problem Application of classical integrators 4 3 Exponential
More informationME8230 Nonlinear Dynamics
ME8230 Nonlinear Dynamics Lecture 1, part 1 Introduction, some basic math background, and some random examples Prof. Manoj Srinivasan Mechanical and Aerospace Engineering srinivasan.88@osu.edu Spring mass
More informationCOMPLETE ALL ROUGH WORKINGS IN THE ANSWER BOOK AND CROSS THROUGH ANY WORK WHICH IS NOT TO BE ASSESSED.
BSc/MSci EXAMINATION PHY-304 Time Allowed: Physical Dynamics 2 hours 30 minutes Date: 28 th May 2009 Time: 10:00 Instructions: Answer ALL questions in section A. Answer ONLY TWO questions from section
More informationChapter 14: Periodic motion
Chapter 14: Periodic motion Describing oscillations Simple harmonic motion Energy of simple harmonic motion Applications of simple harmonic motion Simple pendulum & physical pendulum Damped oscillations
More informationMassachusetts Institute of Technology Department of Physics. Final Examination December 17, 2004
Massachusetts Institute of Technology Department of Physics Course: 8.09 Classical Mechanics Term: Fall 004 Final Examination December 17, 004 Instructions Do not start until you are told to do so. Solve
More informationSignatures of Trans-Planckian Dissipation in Inflationary Spectra
Signatures of Trans-Planckian Dissipation in Inflationary Spectra 3. Kosmologietag Bielefeld Julian Adamek ITPA University Würzburg 8. May 2008 Julian Adamek 1 / 18 Trans-Planckian Dissipation in Inflationary
More informationSolar Sailing near a collinear point
Solar Sailing near a collinear point 6th AIMS International Conference on Dynamical Systems and Differential Equations 26 Ariadna Farrès & Àngel Jorba Departament de Matemàtica Aplicada i Anàlisi Universitat
More informationA path integral approach to the Langevin equation
A path integral approach to the Langevin equation - Ashok Das Reference: A path integral approach to the Langevin equation, A. Das, S. Panda and J. R. L. Santos, arxiv:1411.0256 (to be published in Int.
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 informationI. Dissipative Dynamics
I. Dissipative Dynamics I.A Brownian Motion of a Particle Observations under a microscope indicate that a dust particle in a liquid drop undergoes a random jittery motion. This is because of the random
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 information[#1] R 3 bracket for the spherical pendulum
.. Holm Tuesday 11 January 2011 Solutions to MSc Enhanced Coursework for MA16 1 M3/4A16 MSc Enhanced Coursework arryl Holm Solutions Tuesday 11 January 2011 [#1] R 3 bracket for the spherical pendulum
More informationAB-267 DYNAMICS & CONTROL OF FLEXIBLE AIRCRAFT
FLÁIO SILESTRE DYNAMICS & CONTROL OF FLEXIBLE AIRCRAFT LECTURE NOTES LAGRANGIAN MECHANICS APPLIED TO RIGID-BODY DYNAMICS IMAGE CREDITS: BOEING FLÁIO SILESTRE Introduction Lagrangian Mechanics shall be
More informationGSHMC: An efficient Markov chain Monte Carlo sampling method. Sebastian Reich in collaboration with Elena Akhmatskaya (Fujitsu Laboratories Europe)
GSHMC: An efficient Markov chain Monte Carlo sampling method Sebastian Reich in collaboration with Elena Akhmatskaya (Fujitsu Laboratories Europe) 1. Motivation In the first lecture, we started from a
More informationIntegrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows
Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows Alexander Chesnokov Lavrentyev Institute of Hydrodynamics Novosibirsk, Russia chesnokov@hydro.nsc.ru July 14,
More informationTreecodes for Cosmology Thomas Quinn University of Washington N-Body Shop
Treecodes for Cosmology Thomas Quinn University of Washington N-Body Shop Outline Motivation Multipole Expansions Tree Algorithms Periodic Boundaries Time integration Gravitational Softening SPH Parallel
More informationPhysics 5153 Classical Mechanics. Canonical Transformations-1
1 Introduction Physics 5153 Classical Mechanics Canonical Transformations The choice of generalized coordinates used to describe a physical system is completely arbitrary, but the Lagrangian is invariant
More informationGraceful exit from inflation for minimally coupled Bianchi A scalar field models
Graceful exit from inflation for minimally coupled Bianchi A scalar field models Florian Beyer Reference: F.B. and Leon Escobar (2013), CQG, 30(19), p.195020. University of Otago, Dunedin, New Zealand
More informationP321(b), Assignement 1
P31(b), Assignement 1 1 Exercise 3.1 (Fetter and Walecka) a) The problem is that of a point mass rotating along a circle of radius a, rotating with a constant angular velocity Ω. Generally, 3 coordinates
More informationPendulum with a vibrating base
Pendulum with a vibrating base Chennai, March 14, 2006 1 1 Fast perturbations- rapidly oscillating perturbations We consider perturbations of a Hamiltonian perturbed by rapid oscillations. Later we apply
More informationBoundary conditions for molecular dynamics
Xiantao Li Department of Mathematics Penn State University Collaborators Weinan E and Jerry Z. Yang Princeton University AtC Workshop: Sandia National Lab, Mar 21, 26 Motivations Boundary Reflection Fracture
More informationJ07M.1 - Ball on a Turntable
Part I - Mechanics J07M.1 - Ball on a Turntable J07M.1 - Ball on a Turntable ẑ Ω A spherically symmetric ball of mass m, moment of inertia I about any axis through its center, and radius a, rolls without
More informationLangevin Dynamics of a Single Particle
Overview Langevin Dynamics of a Single Particle We consider a spherical particle of radius r immersed in a viscous fluid and suppose that the dynamics of the particle depend on two forces: A drag force
More informationSketchy Notes on Lagrangian and Hamiltonian Mechanics
Sketchy Notes on Lagrangian and Hamiltonian Mechanics Robert Jones Generalized Coordinates Suppose we have some physical system, like a free particle, a pendulum suspended from another pendulum, or a field
More informationHETEROGENEOUS MULTISCALE METHODS FOR MECHANICAL SYSTEMS WITH VIBRATIONS
HETEROGENEOUS MULTISCALE METHODS FOR MECHANICAL SYSTEMS WITH VIBRATIONS M. P. CALVO AND J. M. SANZ-SERNA Abstract. We are concerned with heterogeneous multiscale methods for the numerical integration of
More information06. Lagrangian Mechanics II
University of Rhode Island DigitalCommons@URI Classical Dynamics Physics Course Materials 2015 06. Lagrangian Mechanics II Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative Commons License
More informationControlled Diffusions and Hamilton-Jacobi Bellman Equations
Controlled Diffusions and Hamilton-Jacobi Bellman Equations Emo Todorov Applied Mathematics and Computer Science & Engineering University of Washington Winter 2014 Emo Todorov (UW) AMATH/CSE 579, Winter
More informationStructural Dynamics Lecture 4. Outline of Lecture 4. Multi-Degree-of-Freedom Systems. Formulation of Equations of Motions. Undamped Eigenvibrations.
Outline of Multi-Degree-of-Freedom Systems Formulation of Equations of Motions. Newton s 2 nd Law Applied to Free Masses. D Alembert s Principle. Basic Equations of Motion for Forced Vibrations of Linear
More information= 0. = q i., q i = E
Summary of the Above Newton s second law: d 2 r dt 2 = Φ( r) Complicated vector arithmetic & coordinate system dependence Lagrangian Formalism: L q i d dt ( L q i ) = 0 n second-order differential equations
More informationThis is a Gaussian probability centered around m = 0 (the most probable and mean position is the origin) and the mean square displacement m 2 = n,or
Physics 7b: Statistical Mechanics Brownian Motion Brownian motion is the motion of a particle due to the buffeting by the molecules in a gas or liquid. The particle must be small enough that the effects
More informationAPPENDIX B. The primitive equations
APPENDIX B The primitive equations The physical and mathematical basis of all methods of dynamical atmospheric prediction lies in the principles of conservation of momentum, mass, and energy. Applied to
More informationSummary of the new Modelling Vocabulary
Summary of the new Modelling Vocabulary These two pages attempts to summarise in a concise manner the Modelling Vocabulary. What are Models? What are Simulations? Materials Models consist of Physics or
More informationLatest thoughts on stochastic kinetic energy backscatter - good and bad
Latest thoughts on stochastic kinetic energy backscatter - good and bad by Glenn Shutts DARC Reading University May 15 2013 Acknowledgments ECMWF for supporting this work Martin Leutbecher Martin Steinheimer
More informationEQUATION LANGEVIN. Physics, Chemistry and Electrical Engineering. World Scientific. With Applications to Stochastic Problems in. William T.
SHANGHAI HONG WorlrfScientific Series krtonttimfjorary Chemical Physics-Vol. 27 THE LANGEVIN EQUATION With Applications to Stochastic Problems in Physics, Chemistry and Electrical Engineering Third Edition
More informationFysikens matematiska metoder week: Semester: Spring 2015 (1FA121) Lokal: Bergsbrunnagatan 15, Sal 2 Tid: 08:00-13:00 General remark: Material al
Fysikens matematiska metoder week: 13-23 Semester: Spring 2015 (1FA121) Lokal: Bergsbrunnagatan 15, Sal 2 Tid: 08:00-13:00 General remark: Material allowed to be taken with you to the exam: Physics Handbook,
More informationState Space Representation of Gaussian Processes
State Space Representation of Gaussian Processes Simo Särkkä Department of Biomedical Engineering and Computational Science (BECS) Aalto University, Espoo, Finland June 12th, 2013 Simo Särkkä (Aalto University)
More information6.12 MECHANICS 4, M4 (4764) A2
6.1 MECHANICS 4, M4 (4764) A Objectives To prepare students for more advanced courses at university by extending the use of calculus in mechanics. Students will be expected to be technically competent
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 informationEnergy and Forces in DFT
Energy and Forces in DFT Total Energy as a function of nuclear positions {R} E tot ({R}) = E DF T ({R}) + E II ({R}) (1) where E DF T ({R}) = DFT energy calculated for the ground-state density charge-density
More informationGaussian Process Approximations of Stochastic Differential Equations
Gaussian Process Approximations of Stochastic Differential Equations Cédric Archambeau Centre for Computational Statistics and Machine Learning University College London c.archambeau@cs.ucl.ac.uk CSML
More information11/17/10. Chapter 14. Oscillations. Chapter 14. Oscillations Topics: Simple Harmonic Motion. Simple Harmonic Motion
11/17/10 Chapter 14. Oscillations This striking computergenerated image demonstrates an important type of motion: oscillatory motion. Examples of oscillatory motion include a car bouncing up and down,
More informationParallel-in-time integrators for Hamiltonian systems
Parallel-in-time integrators for Hamiltonian systems Claude Le Bris ENPC and INRIA Visiting Professor, The University of Chicago joint work with X. Dai (Paris 6 and Chinese Academy of Sciences), F. Legoll
More informationStochastic-Dynamical Thermostats for Constraints and Stiff Restraints
EPJ manuscript No. (will be inserted by the editor) Stochastic-Dynamical Thermostats for Constraints and Stiff Restraints Janis Bajars 1,a, Jason Frank 1,b, and Ben Leimkuhler 1,c 1 Centrum Wiskunde &
More information