BROWNIAN DYNAMICS SIMULATIONS WITH HYDRODYNAMICS. Abstract
|
|
- Oswin Simmons
- 5 years ago
- Views:
Transcription
1 BROWNIAN DYNAMICS SIMULATIONS WITH HYDRODYNAMICS Juan J. Cerdà 1 1 Institut für Computerphysik, Pfaffenwaldring 27, Universität Stuttgart, Stuttgart, Germany. (Dated: July 21, 2009) Abstract - 1
2 I. MAIN RESULTS WE WILL OBTAIN Under the following assumptions τ < τ a = a2 D o Spherical particles of radius a. On colloidal time scales τ τ B the motion of the solvent can be described by the so-called Creeping flow equations (aka Stokes equations ), i.e. p(r) + η o 2 u(r) = 0 (1) u(r) = 0 (2) where p(r) is the pressure of the fluid at point r, and u(r) is the fluid velocity in such point. This assumption implies that the Reynolds number Re = ρ a v η o (3) is low. Stick boundary conditions between the fluid and the surface of the particles. we will show that: 2
3 It is possible to get the following finite difference algorithm (known as Ermak- McCammon scheme r i (t 0 + τ) = r i (t o ) + vi D (X o )τ + (2τ) 1/2 d(x o ) n + o(τ) (4) where v D i (X o ) = [βd ij (X o ) f j (X o ) + j D ij (X o )] (5) and X o = (r 1 (t o ),...r N (t o )) X = (r 1 (t o + τ),...r N (t o + τ)) β = 1 k B T n is a Gaussian Random vector of independently distributed components of mean zero and variance one. Usually computed using a uniform random number χ, such that n R = 12(χ 0.5) d(x) is square-root matrix of the positive definite matrix D which are 3N 3N. D is the translational hydrodynamic diffusivity tensor (or the product β D is known as the translational hydrodynamic mobility tensor (we will usually refer to them as diffusivity tensors. D(X 0 ) = d(x o ) d(x o ) (6) When there is no hydrodynamics: D ij (X) = D o δ ij 1 When there is hydrodynamics, D can be very complex, and obtaining analytical expressions is very difficult. If we take the approximation of just considering pairwise interactions, i.e. given a pair of spheres interacting i and j the other N 2 we assume do not exert any influence on the interaction between i and j, then we will show that it is possible to write where D ij (X) = D far ij (r ij ) + Dij near (X) (7) 3
4 D far ij (r ij ) is the far-field part, valid when the distances between particles i and j is large. A usual expression for the far-field is the so know as Rotne-Prager (RP) approximation Dij RP (r ij ) given by [ ( Dij RP 3 a (r ij ) = δ ij D o 1+(1 δ ij ) D o 4 we will see later where it comes from. D near ij r ij ) [1 + ˆr ij ˆr ij ] + 1 ( ) 3 a [1 3ˆr ij ˆr ij]] 2 r ij (X) is the near-field part (all terms not accounted in Dij RP (r ij )) plus the lubrication effects. It is important when the distances between particles i and j become smaller. We can neglect it if we have for instance a diluted system of particles that repeal each other. We cannot neglect it for instance when we have attraction between particles, or in the case of hard spheres. Another thing we will learn is that, within the approaches we do, the static properties of a system can be computed without any need of including Hydrodyamic Interactions. So, BD without HI can be still very useful!!! One should notice that in general: HI is not pairwise-additive They are long range, decaying as 1/r. They can even diverge for certain types of motion when particles approach very closely (lubrication effects) Implementation of Ermak scheme is in principle O(N 3 ) bottleneck in BD. Few particles can be simulated, typically 30 < N < 300. Let s proof the different results we have stated... (8) 4
5 II. HOW TO GET THE FAR-FIELD : THE ROTNE-PRAGER (RP) APPROX- IMATION A. Steps we will do (1) Derive expressions for the velocity and pressure in a fluid that obeys the Creeping flow equations due to the action of external forces f ext : the Oseen tensor, and the pressure vector. (2) Suppose that the surface elements of the N colloidal spheres cause the external forces f ext that act on the fluid, and we have stick boundary conditions on the surface of the particles. (3) Simplify the Oseen tensor by assuming that the distance between two particles is much larger than the diameter of the particles 2a: the Oseen approximation for the microscopic diffusion matrices. (4) Assume the two-particle level approach. Get the general form of the diffusion matrices. Show that if we stay cut the series at first order we get the Rotne-Prager (RP) approximation. 5
6 B. Derive expressions for the velocity and pressure in a fluid that obeys the Creeping flow equations due to the action of external forces f ext : the Oseen tensor, and the pressure vector. We assume the fluid obeys the so-called Creeping flow equations (aka Stokes equations ), i.e. p(r) + η o 2 u(r) = 0 (9) u(r) = 0 (10) which are linear dif. equations. If we have an external force acting on point r f ext = f o δ(r r ), where the prefactor f o is the total force acting on the fluid f o = dr f ext (r ). (11) (12) then, due to the linearity of the Creeping flow equations it is possible to show that velocity and pressure in the fluid will be proportional to f o, ie v(r) f o p(r) f o and we can write this as the following equalities v(r) = T (r r )f o p(r) = g(r r )f o where T (r r ) is known as the Oseen Tensor (because it seems that Oseen was the first to derive an expression for this tensor). g(r r ) is known as the pressure vector. 6
7 In short, the Oseen Tensor and the pressure vector are the Green s functions of the creeping flow equations. Let s suppose now that our force is not just acting on r but over the entire fluid, then, once again thanks to the linearity of the Creeping equations the velocity and pressure at one given point r is the sum of the contributions acting on that point that arise from all points of the fluid, i.e., we can write the velocity and pressure at point r as: v(r) = dr T (r r ) f ext (r ) (13) p(r) = dr g(r r ) f ext (r ) (14) By substituting previous expressions of v(r) and p(r) into the creeping flow equations, we get (see B1.25) T (r) = g(r) = 1 4π π η o r [ 1 + rr r 2 ] (15) r r 3 (16) 7
8 C. Suppose that the surface elements of the N colloidal spheres cause the external forces f ext that act on the fluid, and we have stick boundary conditions on the surface of the particles. Let s suppose we have N particles, and V j denotes the surface of the spheres. Therefore the velocity and pressure of the fluid at a point r is given by where v(r) = p(r) = ds T (r r ) f j (r ) (17) V j ds g(r r ) f j (r ) (18) V j f j is the force per unit area that a surface element ds of the sphere j which is located at position r exerts on the fluid. We assume stick boundary conditions, i.e. the velocity of the fluid at the sphere surfaces is the same than the velocity of the surfaces. This mathematically can be expressed as v particle = v fluid (r) when r V i (19) v i + Ω (r r i ) = ds T (r r ) f j (r ) when r V i V j (20) Due to symmetry (if there is no external torques) Ω (r r i ) = 0. And we get v i = V j ds T (r r ) f j (r ) when r V i (21) Now we can integrate the previous expression over the whole surface of particle (r V i ) and we get v i = 1 4πa 2 1 4πa 2 V i ds j i V i ds ds T (r r ) f i (r ) + (22) V i ds T (r r ) f j (r ) (23) V j Since we now T, we know that ds T (r r ) = 1 2a V i 3η 0 for r V i (24) 8
9 and therefore V i ds ds T (r r ) f i (r ) = 1 V i 6πη 0 a F i h (25) Fi h (t) = ds f i (r ) (26) V i i.e., F h i (t) is the total force that the fluid exerts on the particle i. Now still we need to compute, 1 4πa 2 j i V i ds V j ds T (r r ) f j (r ) (27) Let s write it in terms of the position vectors of the centers of the spheres br i and br j and the local vectors R = r r i and R = r r j, (plot B1.27), we can get for each term of the sum 1 ds ds T (R R + r 4πa 2 i r j ) f j (R + r j ) (28) V i V j 9
10 D. Simplify the Oseen tensor by assuming that the distance between two particles is much larger than the diameter of the particles 2a: the Oseen approximation for the microscopic diffusion matrices. If R R < 2a, then T (R R + r i r j ) T (r i r j ) we thus get for each term of the sum that 1 ds ds T (r r ) f 4πa 2 j (r ) T (r i r j ) Fj h (29) V i V j Thus finally we get v i = 1 6πη 0 a F i h j i T (r i r j ) F h j (30) the diffusion matrix elements can be obtained by recalling that D ij must satisfy v i = β D ij (X) F h j (31) which leads to the Oseen approximation for the microscopic diffusion matrices if we just take the leading order expansion respect a/r ij, namely D ii = D o 1 (32) D ij k B T T (r i r j ) = 3 4 D a o [1 + ˆr ij ˆr ij ] r ij when i j (33) where r ij = r i r j, and D o = k B T/(6πη o a) (the Stokes-Einstein diffusion coefficient). Remember this expression is only valid for large values of r ij (i.e. small values of a/r ij ) otherwise we need to add more terms to the expansion. 10
11 E. Assume the two-particle level approach. Get the general form of the diffusion matrices. Show that if we stay cut the series at first order we get the Rotne-Prager (RP) approximation If we assume that in the interaction between particle i and j no other of the N 2 particles is able to modify the interaction (i.e. low densities are assumed), then, in this two-particle approach it is possible to show that the general expressions for the diffusion matrices are N D ii = D o 1 + D o [A s (r ij ) ˆr ij ˆr ij + B s (r ij ) [1 ˆr ij ˆr ij ]] (34), j i D ij = D o [A c (r ij ) ˆr ij ˆr ij + B c (r ij ) [1 ˆr ij ˆr ij ]] (35) where the so called mobility functions (s=self, c=cross) are A s (r ij ) = 15 ( ) 4 a + 11 ( ) ( 6 ( ) ) 8 a a + O 4 r ij 2 r ij r ij A c (r ij ) = + 3 ( ) ( ) 3 a a + 75 ( ) ( 7 ( ) ) 9 a a + O 2 r ij r ij 4 r ij r ij B s (r ij ) = 17 ( ) ( 6 ( ) ) 8 a a + O 16 r ij r ij B c (r ij ) = 3 ( ) a + 1 ( ) ( 3 ( ) ) 9 a a + O 4 2 r ij r ij r ij (36) (37) (38) (39) (40) 11
12 III. HOW TO GET THE ERMAK-MCCAMMON SCHEME A. Steps we will do (1) Derivation of the Smoluchowski equation with HI (2) Recast the Smoluchowski equation in terms of the particle drift velocities vi D (X). (3) Assume that the time step τ τ a and therefore the configuration has changed so little that D(X) D(X o ) and vi D (X) vi D (X o ). Get a solution for the simplified Smoluchowski equation. (4) Derive the finite difference algorithm. 12
13 B. Derivation of the Smoluchowski equation with HI What we want to proof is that given a pdf P (X, t) is governed by the equation P (X, t) dt = Ô (...) Ô P (X, t) (41) [ ri D i,j (X) (β (...) rj Φ(X) + rj (...) )] (42) i, There are different ways of getting the Smoluchowski equation with HI (a) Integration of the Liouville equation P (X, t) dt = ˆL P (X, t) (43) where ˆL (...) [ p ] j m r j (...) + [ rj Φ] pj (...) (44) is the Liouville operator. But integrating the Liouville equation respect the phase space coordinates and the momenta of the colloids is very involved (see JCP 54, 3547, (1971); JCP 57, 2098, (1972).) (b) As we will do... Let be now X = (r 0,..., r N ) (only the positions). Then N(t) = dx P (X, t) (45) W where W is an arbitrary volume of the positional phase space. Therefore, given that only the component of dx/dt can lead to a in-/out- flow of probability across the surface W [ ] dn(t) P (X, t) dx = dx = ds (X, t) (46) dt t dt W W Now using Gauss integral theorem [ [ ]] P (X, t) dx dx + X P (X, t) = 0 (47) t dt W where X is respect the 3N position coordinates. Given W is arbitrary, we get [ ] P (X, t) dx = X P (X, t) t dt 13 (48)
14 which can be rephrased as P (X, t) t = rj [v cg j P (X, t)] (49) where cg means coarse-grained to the diffusive time scale. Ok. Now, what we need is an expression for the v cg j. Let s assume that in the scale we measure things the inertial forces are very small (relaxation has already passed) when compared to Forces arising from the total potential Φ The Brownian Force F brow j forces due to the friction with the solvent F h j induces a fluid flow which affects other particles. which make than when a fluid moves, it Given, we assume Creeping flow equations to be valid, its linearity allows to say that the force acting on a sphere j must be linearly proportional to the velocities of all the other spheres, namely F h j = i=1 B ij v cg j (50) where B is the hydrodynamic friction tensor. In the case there is no hydrodynamics, then we just have B ij = γ 1 δ Kronecker ij (51) i.e, the usual friction law. Therefore, in the diffusive time scale, the sum of all previous forces must be zero rj Φ B ij v cg j + F brown j = 0 (52) which can be rewritten in the 3N momenta space notation X as X Φ B ij v cg + F brown = 0 (53) where v cg (v cg i,..., v cg N ), idem for F brown. If we define the diffusion tensor as D k B T B 1 (54) 14
15 we can get by inverting the previous expression 53, and returning to the notation for each particle, that v cg i = D ij [ ] β rj Φ + βfj brown (55) Thus, we get the following expression for the Smoluchowski equation P (X, t) dt = Ô (...) Ô P (X, t) (56) [ ri D i,j (X) (β (...) rj Φ(X) + rj (...) )] (57) i, aka GSE= Generalized Smoluchowski equation with HI. Keep in mind that the previous equation you can just apply it on the Brownian time and length scales, when both particles and solvent are quasi-inertia-free. 15
16 C. Recast the Smoluchowski equation in terms of the particle drift velocities v D i (X). The Smoluchowski equation P (X, t) dt = Ô (...) Ô P (X, t) (58) [ ri D i,j (X) (β (...) rj Φ(X) + rj (...) )] (59) i, can be recast as where P (X, t) t v D i (X) = i (v i D (X) P (X, t) ) + i=1 i [ j (D i,j P (X, t)] (60) i, [β D ij (X) F j (X) + j D ij (X)] (61) v D i (X) are the drift velocities. A new term arises due to hydrodynamics: j j D ij, which physically acts as a repulsive force (tends to send particles to regions where they have more mobility). Notice that P eq (X) exp( βφ(x)) Z (62) where Z is the partition function, is a stationary solution of the Smoluchowski equation. So what? The last fact implies that Hydrodynamics do no have effect on static equilibrium properties!!! (within the approaches we have taken, keep that also in mind) Your simulations without HI should give the correct equilibrium properties, no need of doing for them the costly HI simulations. 16
17 D. Assume that the time step τ τ a and therefore the configuration has changed so little that D(X) D(X o ) and v D i (X) vd i (X o). Get a solution for the simplified Smoluchowski equation. In that case, the GSE reduces to a diffusion equation with constant coefficients, namely (in the 3N position space coordinates) P (X, τ X o ) t = v D (X o ) P (X, τ X o ) + [ (D(X o ) P (X, τ X o ) (63) is It can be proof that the solution to the previous equation, if P (X, 0 X o ) = δ(x X o ) P sol delta (X, τ X o ) = ( (4πτ) 3N/2 1 det(d(xo )) exp 4τ (64) [X X o τv D (X o )] T D 1 (X 0 ) [X X o τv D (X o )] ) (65) which is a Gaussian. 17
18 E. Derive the finite difference algorithm. We can see that the true solution for our system will be P (X, τ X o ) = P sol delta (X, τ X o ) + o(τ) (66) therefore we get the momentum relations (X X o ) o = vi D (X o )τ + o(τ) (67) (X Xo ) 2 o = 2DD ij (X o )τ + o(τ) (68) which can be proof that what we have then is stochastically equivalent, in the limit τ 0 to the Ermak-McCammon scheme r i (t 0 + τ) = r i (t o ) + v D i (X o )τ + (2τ) 1/2 d(x o ) n + o(τ) (69) where v D i (X o ) = [βd ij (X o ) f j (X o ) + j D ij (X o )] (70) 18
Langevin 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 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 informationCHAPTER V. Brownian motion. V.1 Langevin dynamics
CHAPTER V Brownian motion In this chapter, we study the very general paradigm provided by Brownian motion. Originally, this motion is that a heavy particle, called Brownian particle, immersed in a fluid
More informationRouse vs. Zimm Regime: Hydrodynamic Interactions 1, a)
THEOY OF POLYME DYNAMICS ouse vs. Zimm egime: Hydrodynamic Interactions 1, a) Henri Menke Institute for Computational Physics, Universität Stuttgart, Allmandring 3, D-70569 Stuttgart, Germany (Dated: February
More informationProtein Dynamics, Allostery and Function
Protein Dynamics, Allostery and Function Lecture 3. Protein Dynamics Xiaolin Cheng UT/ORNL Center for Molecular Biophysics SJTU Summer School 2017 1 Obtaining Dynamic Information Experimental Approaches
More informationStatistical Mechanics of Active Matter
Statistical Mechanics of Active Matter Umberto Marini Bettolo Marconi University of Camerino, Italy Naples, 24 May,2017 Umberto Marini Bettolo Marconi (2017) Statistical Mechanics of Active Matter 2017
More informationFokker-Planck Equation with Detailed Balance
Appendix E Fokker-Planck Equation with Detailed Balance A stochastic process is simply a function of two variables, one is the time, the other is a stochastic variable X, defined by specifying: a: the
More informationBrownian Motion: Fokker-Planck Equation
Chapter 7 Brownian Motion: Fokker-Planck Equation The Fokker-Planck equation is the equation governing the time evolution of the probability density of the Brownian particla. It is a second order differential
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 information3.1 Hydrodynamic interactions in a Gaussian chain
Chapter 3 The Zimm model 3. Hydrodynamic interactions in a Gaussian chain In the previous chapter we have focused on the Rouse chain, which gives a good description of the dynamics of unentangled concentrated
More informationCoupling an Incompressible Fluctuating Fluid with Suspended Structures
Coupling an Incompressible Fluctuating Fluid with Suspended Structures Aleksandar Donev Courant Institute, New York University & Rafael Delgado-Buscalioni, UAM Florencio Balboa Usabiaga, UAM Boyce Griffith,
More informationHydrodynamic interactions and wall drag effect in colloidal suspensions and soft matter systems
Hydrodynamic interactions and wall drag effect in colloidal suspensions and soft matter systems Maciej Lisicki" Institute of Theoretical Physics" Faculty of Physics, University of Warsaw" Poland SOMATAI
More informationBrownian motion and the Central Limit Theorem
Brownian motion and the Central Limit Theorem Amir Bar January 4, 3 Based on Shang-Keng Ma, Statistical Mechanics, sections.,.7 and the course s notes section 6. Introduction In this tutorial we shall
More information2.5 Stokes flow past a sphere
Lecture Notes on Fluid Dynamics.63J/.J) by Chiang C. Mei, MIT 007 Spring -5Stokes.tex.5 Stokes flow past a sphere Refs] Lamb: Hydrodynamics Acheson : Elementary Fluid Dynamics, p. 3 ff One of the fundamental
More informationComputer simulation on the collision-sticking dynamics of two colloidal particles in an optical trap
THE JOURNAL OF CHEMICAL PHYSICS 126, 144903 2007 Computer simulation on the collision-sticking dynamics of two colloidal particles in an optical trap Shenghua Xu and Zhiwei Sun a NML, Institute of Mechanics,
More informationPolymer dynamics. Course M6 Lecture 5 26/1/2004 (JAE) 5.1 Introduction. Diffusion of polymers in melts and dilute solution.
Course M6 Lecture 5 6//004 Polymer dynamics Diffusion of polymers in melts and dilute solution Dr James Elliott 5. Introduction So far, we have considered the static configurations and morphologies of
More informationAn electrokinetic LB based model for ion transport and macromolecular electrophoresis
An electrokinetic LB based model for ion transport and macromolecular electrophoresis Raffael Pecoroni Supervisor: Michael Kuron July 8, 2016 1 Introduction & Motivation So far an mesoscopic coarse-grained
More informationModeling and simulation of nanoparticle aggregation in colloidal systems
Retrospective Theses and Dissertations Iowa State University Capstones, Theses and Dissertations 2008 Modeling and simulation of nanoparticle aggregation in colloidal systems Sergiy Markutsya Iowa State
More informationFluid-Particles Interaction Models Asymptotics, Theory and Numerics I
Fluid-Particles Interaction Models Asymptotics, Theory and Numerics I J. A. Carrillo collaborators: T. Goudon (Lille), P. Lafitte (Lille) and F. Vecil (UAB) (CPDE 2005), (JCP, 2008), (JSC, 2008) ICREA
More informationPolymer Solution Thermodynamics:
Polymer Solution Thermodynamics: 3. Dilute Solutions with Volume Interactions Brownian particle Polymer coil Self-Avoiding Walk Models While the Gaussian coil model is useful for describing polymer solutions
More informationFluid Mechanics II Viscosity and shear stresses
Fluid Mechanics II Viscosity and shear stresses Shear stresses in a Newtonian fluid A fluid at rest can not resist shearing forces. Under the action of such forces it deforms continuously, however small
More informationPolymer Dynamics and Rheology
Polymer Dynamics and Rheology 1 Polymer Dynamics and Rheology Brownian motion Harmonic Oscillator Damped harmonic oscillator Elastic dumbbell model Boltzmann superposition principle Rubber elasticity and
More informationAA210A Fundamentals of Compressible Flow. Chapter 5 -The conservation equations
AA210A Fundamentals of Compressible Flow Chapter 5 -The conservation equations 1 5.1 Leibniz rule for differentiation of integrals Differentiation under the integral sign. According to the fundamental
More informationSmoluchowski Diffusion Equation
Chapter 4 Smoluchowski Diffusion Equation Contents 4. Derivation of the Smoluchoswki Diffusion Equation for Potential Fields 64 4.2 One-DimensionalDiffusoninaLinearPotential... 67 4.2. Diffusion in an
More informationComputational Fluid Dynamics Prof. Dr. Suman Chakraborty Department of Mechanical Engineering Indian Institute of Technology, Kharagpur
Computational Fluid Dynamics Prof. Dr. Suman Chakraborty Department of Mechanical Engineering Indian Institute of Technology, Kharagpur Lecture No. # 02 Conservation of Mass and Momentum: Continuity and
More informationActive Matter Lectures for the 2011 ICTP School on Mathematics and Physics of Soft and Biological Matter Lecture 3: Hydrodynamics of SP Hard Rods
Active Matter Lectures for the 2011 ICTP School on Mathematics and Physics of Soft and Biological Matter Lecture 3: of SP Hard Rods M. Cristina Marchetti Syracuse University Baskaran & MCM, PRE 77 (2008);
More informationThe Navier-Stokes equations for an incompressible fluid in indicial notation are. i u i = 0 (9.1)
Chapter 9 Viscous flow The Navier-Stokes equations for an incompressible fluid in indicial notation are i u i = 0 9.1 ρ t + u j j u i = i p + µ 2 j u i 9.2 When these are applied to practical problems,
More informationDiscontinuous shear thickening on dense non-brownian suspensions via lattice Boltzmann method
Discontinuous shear thickening on dense non-brownian suspensions via lattice Boltzmann method Pradipto and Hisao Hayakawa Yukawa Insitute for Theoretical Physics Kyoto University Rheology of disordered
More informationFrom a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction
From a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction Carnegie Mellon University Center for Nonlinear Analysis Working Group, October 2016 Outline 1 Physical Framework 2 3 Free Energy
More information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
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 informationFORCE ENERGY. only if F = F(r). 1 Nano-scale (10 9 m) 2 Nano to Micro-scale (10 6 m) 3 Nano to Meso-scale (10 3 m)
What is the force? CEE 618 Scientific Parallel Computing (Lecture 12) Dissipative Hydrodynamics (DHD) Albert S. Kim Department of Civil and Environmental Engineering University of Hawai i at Manoa 2540
More informationFluid dynamics is the study of flows of fluid. However, before we begin discussing the dynamics of
Hydrodynamics of Newtonian Fluids Cameron Stewart Universität Stuttgart cameron.n.stewart@gmail.com April 5, 2017 Abstract In this review we provide an overview of the hydrodynamics of viscous, Newtonian
More informationEmergence of collective dynamics in active biological systems -- Swimming micro-organisms --
12/08/2015, YITP, Kyoto Emergence of collective dynamics in active biological systems -- Swimming micro-organisms -- Norihiro Oyama John J. Molina Ryoichi Yamamoto* Department of Chemical Engineering,
More informationMesoscale Simulation Methods. Ronojoy Adhikari The Institute of Mathematical Sciences Chennai
Mesoscale Simulation Methods Ronojoy Adhikari The Institute of Mathematical Sciences Chennai Outline What is mesoscale? Mesoscale statics and dynamics through coarse-graining. Coarse-grained equations
More informationPart III. Polymer Dynamics molecular models
Part III. Polymer Dynamics molecular models I. Unentangled polymer dynamics I.1 Diffusion of a small colloidal particle I.2 Diffusion of an unentangled polymer chain II. Entangled polymer dynamics II.1.
More informationFast Brownian Dynamics for Colloidal Suspensions
Fast Brownian Dynamics for Colloidal Suspensions Aleksandar Donev, CIMS and collaborators: Florencio Balboa CIMS) Andrew Fiore and James Swan MIT) Courant Institute, New York University Modeling Complex
More informationON ALGORITHMS FOR BROWNIAN DYNAMICS COMPUTER SIMULATIONS
COMPUTATIONAL METHODS IN SCIENCE AND TECHNOLOGY 4,35-42 (1998) ON ALGORITHMS FOR BROWNIAN DYNAMICS COMPUTER SIMULATIONS ARKADIUSZ C. BRAŃKA Institute of Molecular Physics, Polish Academy of Sciences, Smoluchowskiego
More informationCEE 618 Scientific Parallel Computing (Lecture 12)
1 / 26 CEE 618 Scientific Parallel Computing (Lecture 12) Dissipative Hydrodynamics (DHD) Albert S. Kim Department of Civil and Environmental Engineering University of Hawai i at Manoa 2540 Dole Street,
More informationMicrostructural theory and the rheology of. concentrated colloidal suspensions
Under consideration for publication in J. Fluid Mech. 1 Microstructural theory and the rheology of concentrated colloidal suspensions E H S S A N N A Z O C K D A S T 1 and J E F F R E Y F. M O R R I S
More informationArtificial Intelligence & Neuro Cognitive Systems Fakultät für Informatik. Robot Dynamics. Dr.-Ing. John Nassour J.
Artificial Intelligence & Neuro Cognitive Systems Fakultät für Informatik Robot Dynamics Dr.-Ing. John Nassour 25.1.218 J.Nassour 1 Introduction Dynamics concerns the motion of bodies Includes Kinematics
More informationThe 5 basic equations of semiconductor device physics: We will in general be faced with finding 5 quantities:
6.012 - Electronic Devices and Circuits Solving the 5 basic equations - 2/12/08 Version The 5 basic equations of semiconductor device physics: We will in general be faced with finding 5 quantities: n(x,t),
More informationLecture 19: Biomixing, part 4: Viscous swimmer
Lecture 9: Biomixing, part 4: Viscous swimmer The dumbbell swimmer The Green s function for the Stokes equation p + µ 2 u = fδ(r), u = 0, () is f G(r), where G(r) is the Oseen tensor: G(r) = 8πµ r ( I
More informationBrownian Motion and Langevin Equations
1 Brownian Motion and Langevin Equations 1.1 Langevin Equation and the Fluctuation- Dissipation Theorem The theory of Brownian motion is perhaps the simplest approximate way to treat the dynamics of nonequilibrium
More informationChapter 2: Fluid Dynamics Review
7 Chapter 2: Fluid Dynamics Review This chapter serves as a short review of basic fluid mechanics. We derive the relevant transport equations (or conservation equations), state Newton s viscosity law leading
More informationCircular motion. Aug. 22, 2017
Circular motion Aug. 22, 2017 Until now, we have been observers to Newtonian physics through inertial reference frames. From our discussion of Newton s laws, these are frames which obey Newton s first
More informationTransport (kinetic) phenomena: diffusion, electric conductivity, viscosity, heat conduction...
Transport phenomena 1/16 Transport (kinetic) phenomena: diffusion, electric conductivity, viscosity, heat conduction... Flux of mass, charge, momentum, heat,...... J = amount (of quantity) transported
More information8.6 Drag Forces in Fluids
86 Drag Forces in Fluids When a solid object moves through a fluid it will experience a resistive force, called the drag force, opposing its motion The fluid may be a liquid or a gas This force is a very
More information7 The Navier-Stokes Equations
18.354/12.27 Spring 214 7 The Navier-Stokes Equations In the previous section, we have seen how one can deduce the general structure of hydrodynamic equations from purely macroscopic considerations and
More informationRelativistic Hydrodynamics L3&4/SS14/ZAH
Conservation form: Remember: [ q] 0 conservative div Flux t f non-conservative 1. Euler equations: are the hydrodynamical equations describing the time-evolution of ideal fluids/plasmas, i.e., frictionless
More informationEXTREMELY CHARGED STATIC DUST DISTRIBUTIONS IN GENERAL RELATIVITY
arxiv:gr-qc/9806038v1 8 Jun 1998 EXTREMELY CHARGED STATIC DUST DISTRIBUTIONS IN GENERAL RELATIVITY METÍN GÜRSES Mathematics department, Bilkent University, 06533 Ankara-TURKEY E-mail: gurses@fen.bilkent.edu.tr
More information16. Rotational Dynamics
6. Rotational Dynamics A Overview In this unit we will address examples that combine both translational and rotational motion. We will find that we will need both Newton s second law and the rotational
More information15.3 The Langevin theory of the Brownian motion
5.3 The Langevin theory of the Brownian motion 593 Now, by virtue of the distribution (2), we obtain r(t) =; r 2 (t) = n(r,t)4πr 4 dr = 6Dt t, (22) N in complete agreement with our earlier results, namely
More informationChapter 1. Governing Equations of GFD. 1.1 Mass continuity
Chapter 1 Governing Equations of GFD The fluid dynamical governing equations consist of an equation for mass continuity, one for the momentum budget, and one or more additional equations to account for
More informationSupporting Information
Supporting Information A: Calculation of radial distribution functions To get an effective propagator in one dimension, we first transform 1) into spherical coordinates: x a = ρ sin θ cos φ, y = ρ sin
More informationLecture IX: Field equations, cosmological constant, and tides
Lecture IX: Field equations, cosmological constant, and tides Christopher M. Hirata Caltech M/C 350-17, Pasadena CA 91125, USA (Dated: October 28, 2011) I. OVERVIEW We are now ready to construct Einstein
More informationBiomolecular hydrodynamics
Biomolecular hydrodynamics Chem 341, Fall, 2014 1 Frictional coefficients Consider a particle moving with velocity v under the influence of some external force F (say a graviational or electrostatic external
More informationStatistical Mechanics of Coarse-Graining
Statistical Mechanics of Coarse-Graining Pep Español Departamento de Física Fundamental, UNED, Apartado 6141, 288 Madrid, Spain 1 Introduction Complex fluids are fluids which despite their continuum appearance
More informationChapter 7. Kinematics. 7.1 Tensor fields
Chapter 7 Kinematics 7.1 Tensor fields In fluid mechanics, the fluid flow is described in terms of vector fields or tensor fields such as velocity, stress, pressure, etc. It is important, at the outset,
More informationC C C C 2 C 2 C 2 C + u + v + (w + w P ) = D t x y z X. (1a) y 2 + D Z. z 2
This chapter provides an introduction to the transport of particles that are either more dense (e.g. mineral sediment) or less dense (e.g. bubbles) than the fluid. A method of estimating the settling velocity
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 information18 Green s function for the Poisson equation
8 Green s function for the Poisson equation Now we have some experience working with Green s functions in dimension, therefore, we are ready to see how Green s functions can be obtained in dimensions 2
More informationIn this approach, the ground state of the system is found by modeling a diffusion process.
The Diffusion Monte Carlo (DMC) Method In this approach, the ground state of the system is found by modeling a diffusion process. Diffusion and random walks Consider a random walk on a lattice with spacing
More informationHydrodynamics in Espresso. Kai Grass 1
Hydrodynamics in Espresso Kai Grass 1 Hydrodynamics in Espresso Explicit solvent: Lennard-Jones liquid at constant energy Dissipative Particle Dynamics (DPD): thermostated particle fluid Lattice Boltzmann
More informationCoupling an Incompressible Fluctuating Fluid with Suspended Structures
Coupling an Incompressible Fluctuating Fluid with Suspended Structures Aleksandar Donev Courant Institute, New York University & Rafael Delgado-Buscalioni, UAM Florencio Balboa Usabiaga, UAM Boyce Griffith,
More informationOn the hydrodynamic diffusion of rigid particles
On the hydrodynamic diffusion of rigid particles O. Gonzalez Introduction Basic problem. Characterize how the diffusion and sedimentation properties of particles depend on their shape. Diffusion: Sedimentation:
More informationMathematical Concepts & Notation
Mathematical Concepts & Notation Appendix A: Notation x, δx: a small change in x t : the partial derivative with respect to t holding the other variables fixed d : the time derivative of a quantity that
More informationSalmon: Lectures on partial differential equations
6. The wave equation Of the 3 basic equations derived in the previous section, we have already discussed the heat equation, (1) θ t = κθ xx + Q( x,t). In this section we discuss the wave equation, () θ
More informationDynamics of interacting vortices on trapped Bose-Einstein condensates. Pedro J. Torres University of Granada
Dynamics of interacting vortices on trapped Bose-Einstein condensates Pedro J. Torres University of Granada Joint work with: P.G. Kevrekidis (University of Massachusetts, USA) Ricardo Carretero-González
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF PHYSICS. DOCTORAL GENERAL EXAMINATION PART 1 August 27, 2012 SOLUTIONS
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF PHYSICS EDUCATION OFFICE PHONE (617)-253-4842 ROOM 4-352 FAX (617)-258-8319 DOCTORAL GENERAL EXAMINATION PART 1 August 27, 2012 SOLUTIONS 1 GROUP I Solutions
More informationElectrical Transport. Ref. Ihn Ch. 10 YC, Ch 5; BW, Chs 4 & 8
Electrical Transport Ref. Ihn Ch. 10 YC, Ch 5; BW, Chs 4 & 8 Electrical Transport The study of the transport of electrons & holes (in semiconductors) under various conditions. A broad & somewhat specialized
More informationEffects of Forcing Scheme on the Flow and the Relative Motion of Inertial Particles in DNS of Isotropic Turbulence
Effects of Forcing Scheme on the Flow and the Relative Motion of Inertial Particles in DNS of Isotropic Turbulence Rohit Dhariwal PI: Sarma L. Rani Department of Mechanical and Aerospace Engineering The
More informationLecture 3: From Random Walks to Continuum Diffusion
Lecture 3: From Random Walks to Continuum Diffusion Martin Z. Bazant Department of Mathematics, MIT February 3, 6 Overview In the previous lecture (by Prof. Yip), we discussed how individual particles
More informationOrigins of Mechanical and Rheological Properties of Polymer Nanocomposites. Venkat Ganesan
Department of Chemical Engineering University of Texas@Austin Origins of Mechanical and Rheological Properties of Polymer Nanocomposites Venkat Ganesan $$$: NSF DMR, Welch Foundation Megha Surve, Victor
More informationPhysics 212: Statistical mechanics II Lecture IV
Physics 22: Statistical mechanics II Lecture IV Our program for kinetic theory in the last lecture and this lecture can be expressed in the following series of steps, from most exact and general to most
More informationPHYSICS DEPARTMENT, PRINCETON UNIVERSITY PHYSICS 301 FINAL EXAMINATION. January 13, 2005, 7:30 10:30pm, Jadwin A10 SOLUTIONS
PHYSICS DEPARTMENT, PRINCETON UNIVERSITY PHYSICS 301 FINAL EXAMINATION January 13, 2005, 7:30 10:30pm, Jadwin A10 SOLUTIONS This exam contains five problems. Work any three of the five problems. All problems
More informationarxiv: v1 [cond-mat.soft] 7 Mar 2012
arxiv:1203.1409v1 [cond-mat.soft] 7 Mar 2012 Hydrodynamic effects in driven soft matter Manoel Manghi, 1 Xaver Schlagberger, 2 Yong-Woon Kim, 3 and Roland R. Netz 2 1 Laboratoire de Physique Théorique,
More informationThe Quantum Hydrodynamic Smooth Effective Potential
The Quantum Hydrodynamic Smooth Effective Potential Carl L. Gardner and Christian Ringhofer Department of Mathematics Arizona State University Tempe, AZ 85287-184 carl.gardner@asu.edu Abstract An extension
More information(2.1) Is often expressed using a dimensionless drag coefficient:
1. Introduction Multiphase materials occur in many fields of natural and engineering science, industry, and daily life. Biological materials such as blood or cell suspensions, pharmaceutical or food products,
More informationKerr black hole and rotating wormhole
Kerr Fest (Christchurch, August 26-28, 2004) Kerr black hole and rotating wormhole Sung-Won Kim(Ewha Womans Univ.) August 27, 2004 INTRODUCTION STATIC WORMHOLE ROTATING WORMHOLE KERR METRIC SUMMARY AND
More information4 Evolution of density perturbations
Spring term 2014: Dark Matter lecture 3/9 Torsten Bringmann (torsten.bringmann@fys.uio.no) reading: Weinberg, chapters 5-8 4 Evolution of density perturbations 4.1 Statistical description The cosmological
More informationLecture 4: Hydrodynamics of Bacterial Suspensions. Plan
Lecture 4: Hydrodynamics of Bacterial Suspensions M. Cristina Marchetti Syracuse University Boulder School 2009 Aphrodite Ahmadi (SU SUNY Cortland) Shiladitya Banerjee (SU) Aparna Baskaran (SU Brandeis)
More informationCritical Phenomena under Shear Flow
Critical Phenomena under Shear Flow Pavlik Lettinga, Hao Wang, Jan K.G. Dhont Close to a gas-liquid critical point, effective interactions between particles become very long ranged, and the dynamics of
More informationHow DLS Works: Interference of Light
Static light scattering vs. Dynamic light scattering Static light scattering measures time-average intensities (mean square fluctuations) molecular weight radius of gyration second virial coefficient Dynamic
More informationPHY481 - Outline of solutions to Homework 3
1 PHY481 - Outline of solutions to Homework 3 Problem 3.8: We consider a charge outside a conducting sphere that is neutral. In order that the sphere be neutral, we have to introduce a new image charge
More informationAsk class: what is the Minkowski spacetime in spherical coordinates? ds 2 = dt 2 +dr 2 +r 2 (dθ 2 +sin 2 θdφ 2 ). (1)
1 Tensor manipulations One final thing to learn about tensor manipulation is that the metric tensor is what allows you to raise and lower indices. That is, for example, v α = g αβ v β, where again we use
More informationDresden, September 20-24, 2010 by Hartmut Löwen
Computer simulations of colloidal dispersions Outline 1) Introduction 2) Colloidal sedimentation 3) Lane formation in driven colloids 4) Band formation in oscillatory fields 5) Lane formation in complex
More informationActive microrheology: Brownian motion, strong forces and viscoelastic media
OR1394 DFG Research Unit Nonlinear Response to Probe Vitrification Active microrheology: Brownian motion, strong forces and viscoelastic media Matthias Fuchs Fachbereich Physik, Universität Konstanz New
More informationPhysics 142 Energy in Mechanics Page 1. Energy in Mechanics
Physics 4 Energy in Mechanics Page Energy in Mechanics This set of notes contains a brief review of the laws and theorems of Newtonian mechanics, and a longer section on energy, because of its central
More informationTopics covered so far:
Topics covered so far: Chap 1: The kinetic theory of gases P, T, and the Ideal Gas Law Chap 2: The principles of statistical mechanics 2.1, The Boltzmann law (spatial distribution) 2.2, The distribution
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 informationRigid bodies - general theory
Rigid bodies - general theory Kinetic Energy: based on FW-26 Consider a system on N particles with all their relative separations fixed: it has 3 translational and 3 rotational degrees of freedom. Motion
More informationGeometrized units. Specific energy and specific angular momentum
In this lecture we will continue our discussion of general relativity. We first introduce a convention that allows us to drop the many factors of G and c that appear in formulae, then talk in more detail
More informationEffects of Forcing Scheme on the Flow and the Relative Motion of Inertial Particles in DNS of Isotropic Turbulence
Effects of Forcing Scheme on the Flow and the Relative Motion of Inertial Particles in DNS of Isotropic Turbulence Rohit Dhariwal and Vijaya Rani PI: Sarma L. Rani Department of Mechanical and Aerospace
More informationMicrostructural studies for rheology. Chapter 7: Microstructural studies for rheology. Separation of length scales. Micro & macro views
Chapter 7: Microstructural studies for rheology Microstructural studies for rheology To calculate the flow of complex fluids, need governing equations, in particular, the constitutive equation relating
More informationFinal Review Prof. WAN, Xin
General Physics I Final Review Prof. WAN, Xin xinwan@zju.edu.cn http://zimp.zju.edu.cn/~xinwan/ About the Final Exam Total 6 questions. 40% mechanics, 30% wave and relativity, 30% thermal physics. Pick
More informationLinear Theory of Evolution to an Unstable State
Chapter 2 Linear Theory of Evolution to an Unstable State c 2012 by William Klein, Harvey Gould, and Jan Tobochnik 1 October 2012 2.1 Introduction The simple theory of nucleation that we introduced in
More informationDiscontinuous Shear Thickening
Discontinuous Shear Thickening dynamic jamming transition Ryohei Seto, Romain Mari, Jeffrey F. Morris, Morton M. Denn Levich Institute, City College of New York First experimental data Williamson and Hecker
More informationTemporal Integrators for Langevin Equations with Applications to Fluctuating Hydrodynamics and Brownian Dynamics
Temporal Integrators for Langevin Equations with Applications to Fluctuating Hydrodynamics and Brownian Dynamics by Steven Delong A dissertation submitted in partial fulfillment of the requirements for
More information4. The Green Kubo Relations
4. The Green Kubo Relations 4.1 The Langevin Equation In 1828 the botanist Robert Brown observed the motion of pollen grains suspended in a fluid. Although the system was allowed to come to equilibrium,
More information