Introduction to multiscale modeling and simulation. Almost every system is multiscale.
|
|
- Kimberly Matthews
- 6 years ago
- Views:
Transcription
1 Introduction to multiscale modeling and simulation Giovanni Samaey, Scientific Computing Dept. of Computer Science, K.U.Leuven Lecture 1: Course introduction Almost every system is multiscale. We are interested in behavior on large space and long time scales. The phenomenon is influenced by effects taking place on smaller spatial scales and shorter time scales.
2 We don t always notice this multiscale behavior. Microscopically, molecules are subject to Newton s equations of motion. Macroscopically, partial differential equations often describe evolution of density, momentum, energy to sufficient accurately. The microscopic description is then unnecessary! These phenomenological models often predate a microscopic understanding, and were only later derived from first principles. Multiscale problems cannot always be reduced. Bacterial chemotaxis : biased motion towards chemoattractant ( food ) Motion often steered by internal chemical networks source: Wikipedia taken from S. Lubkin microscopic level modeling level detailed individual-based models macroscopic level observation/simulation simplified advection-diffusion PDE
3 Why are multiscale problems difficult to simulate? Effects of different time scales - Micro: small time-steps necessary - Macro: large time intervals Effects of different spatial scales - Micro: fine grid necessary - Macro: large spatial domains y(t) t Effects of different modeling scales - Micro: many degrees of freedom - Macro: averaged effect Often all these effects simultaneously! What will we do in this course? Multiscale model problems (3 lectures) - Relation between microscopic and macroscopic levels of description - Analytic reduction to macroscopic model and reconstruction of microscopic state Multiscale numerical techniques (6 lectures) - Exploit scale separation to reduce computational cost when analytical reduction is not possible - Accurately transfer information from micro to macro and back again Evaluation : read, summarize and evaluate a research paper in the field
4 Disclaimer The literature on analytical and numerical multiscale techniques is - enormous - spread out over many disconnected disciplines - in continuous evolution In this course, you will learn - a particular way of thinking about multiscale problems - that you can use as a starting point to dive into the literature In particular, this course is very much biased by my own research expertise! What will we do in this lecture? Present a number of applications in which multiple scales arise Isolate a single multiscale feature that is relevant in that application The main goals of this lecture are to : - illustrate that multiscale problems are ubiquitous in science and engineering - identify prior knowledge (or gaps therein) in this diverse audience (to adjust the exposition in the remaining lectures) - give an overview of the type of problems that will be tackled (first analytically, then numerically) throughout the remainder of the course
5 Time-scale separation Application 1 : Chemical kinetics Application 2 : A mechanical system (Dahlquist s clock) Application 1 : Chemical kinetics Chemical reactions and mass balance laws Kinetics of enzyme-substrate reaction (Michaelis-Menten reaction) k1 k 3 E+P E + S ES k2 Mass action laws describe evolution of concentration via set of ordinary differential equations (ODEs) d[e]/dt = d[s]/dt = d[es]/dt = d[p ]/dt = k1 [E][S] + k2 [ES] + k3 [ES] k1 [E][S] + k2 [ES] k1 [E][S] k2 [ES] k3[ ES] k3 [ES]
6 Reduction and nondimensionalization Due to mass conservation, we have So that the system can be reduced to a system of 2 ODEs Now we nondimensionalize... t = t/(k 1 E tot ) x =[S]/S tot y =[ES]/E tot... and introduce the nondimensional parameters to obtain [E]+[ES] E tot [S]+[ES]+[P ] S tot d[s]/dt = k 1 E tot [S]+(k 1 [S]+k 2 )[ES] d[es]/dt = k 1 E tot [S] (k 1 [S]+k 2 + k 3 )[ES] E + S k 1 ES k 3 E + P k2 = E tot /S tot λ = k 3 /(k 1 S tot ) κ(k 2 + k 3 )/(k 1 S tot ) x = x +(x + κ λ)y y = x (x + κ)y Time evolution with separated time scales E + S k 1 ES k 3 E + P k2 1 x = x +(x + κ λ)y y = x (x + κ)y 1 := E tot S tot x(t) 0.6 y(t) t Equation for y evolves on a much faster time scale than equation for x Small time-steps are needed to simulate the y-dynamics After a fast initial transient, the system has forgotten its initial condition in y (all behavior in y is very close to each other) t
7 Phase portrait and slaving An informal derivation of a macroscopic model y x = x +(x + κ λ)y y = x (x + κ)y κ =1 λ =1/ x We see that y is uniquely determined by x after a transient. y(x) x x +1 This phenomenon is called slaving, and implies that an approximate macroscopic model exists X X = 2(X + 1) Because the small parameter disappeared, larger time-steps can be taken. Application 2 : Dahlquist s clock A very informal averaging argument Very fast vibrations Slow movement (towards the edge of the table...) Equations of motion x = v v = g(x, v) Instead of approaching a fixed point for each x, v oscillates quickly. The averaged effect of this velocity is felt on long time scales. We would again like to obtain a macroscopic model of the form X = V (X)
8 Time-scale separation : two types of problems Model problem 1 : ODE with fast slaving/invariant manifold structure - Fast variables slaved to (uniquely determined by) slow variables - Slaving relation is used to obtain a model for the slow variables alone Model problem 2 : averaging/homogenization for ODEs - Fast variables behave oscillatory for every value of the slow variables - Appropriately averaged fast behavior drives slow evolution Later on in the course : - derivation of approximate macroscopic models in idealized situations (lecture 2) - numerical methods that exploit scale separation when deriving the macroscopic model analytically is not possible (second part of the course) Space-scale separation Application 1 : Flow through porous media (groundwater, oil reservoir, nuclear waste, CO2 storage,...) Application 2 : Deformation of solid materials with local defects (cracks) The examples are equilibrium problems (no time dependence). If time dependent, the problems of time-scale separation are added.
9 Application 1 : Flow through porous media Darcy s law and groundwater flow equation Darcy s law : flow velocity is proportional to the pressure gradient q(x) = κ p(x) The proportionality constant is called the permeability... κ κ(x; )... and it depends on space and a small-scale parameter. We also have conservation of mass (fluid is incompressible) div q(x) =f(x) Combined, this gives the groundwater flow equation div κ(x; ) p(x) =f(x) source term Heterogeneity of the permeability div κ(x; ) p(x) =f(x) The heterogeneous permeability cannot be measured everywhere. It is therefore modeled as a random field : denotes an average over all possible realizations - Decompose into (space-independent) average + fluctuations κ(x; ) =κ + δκ(x; ) - The fluctuations are correlated with some fine-scale correlation length δκ(x; ) δκ(x ; ) = w(x x ; ) - Generalization to anisotropic fields is also possible One needs to simulate on a very fine grid to sample the heterogeneity.
10 div κ(x; ) p(x) =f(x) Homogenized equation with effective permeability To reduce computational cost, we want to solve on a coarser mesh. First idea : Underresolve (evaluate the heterogeneous permeability on the coarse mesh). This clearly does not converge. Second idea : Homogenize the heterogeneous coefficient. Two options : div κ p(x) =f(x) Replace the permeability coefficient by its mean. This does not converge either (but it is -at this point- not so clear why). Perform a mathematical homogenization to obtain an effective permeability. Application 2 : Solid materials Continuum modeling of deformation Deformed position is original position plus a displacement ϕ(x) =x + u(x) x ϕ(x) Equilibrium of stresses (ignoring boundary conditions) div σ(x) =f(x) System is closed by a constitutive law σ(x) σ(x, u(x),...) For elastic materials, stress only depends on deformation gradient F (x) = ϕ(x) σ σ(f ) If we have hyperelasticity, the stress is derived from an energy σ(f )= W(F )/ F The problem then becomes an energy minimization problem inf W ( ϕ)+...
11 Localized breakdown of continuum description Intendation experiment (large forces locally) Fracture propagation Some problems require an incredibly fine mesh locally. The mesh becomes commensurate with the size of individual atoms. In that region, the continuum description breaks down (cannot be used). Locally, an atomistic description is necessary. Atomistic modeling of deformation Microscopically, a solid material consists of atoms, aligned on a grid. Each atom is displaced Y i = ϕ(x i ), i I There is an interaction potential between every two atoms V (r ij )=V ( Y i Y j ) The energy of the sample is then 1 V (r ij ) 2 i,j, i=j
12 Space-scale separation : two types of problems Model problem 1 : homogenization of PDEs - Partial differential equation with highly heterogeneous coefficients - Goal : obtain/simulate an effective equation with homogenized coefficient Model problem 2 : localized breakdown of continuum descriptions - Macroscopic description: partial differential equation - Microscopic description: particle (atomistic) model - Goal : relate the two level of description/simulate using domain decomposition Later on in the course : - derivation of approximate macroscopic models in idealized situations (lecture 3) - numerical methods that exploit scale separation to simulate when deriving the macroscopic model analytically is not possible (second part of the course) Model-scale separation Application 1 : Complex fluids (blood, liquid crystals, polymers,...) Application 2 : Kinetic models (plasma physics, fluid mechanics, chemotaxis)
13 Before proceeding to application 1... A brief introduction to stochastic differential equations Consider the movement of a heavy particle in a solvent The heavy particle is bombarded with solvent molecules. y The effect is an excitation into random fluctuations + dissipation due to friction. The resulting motion is very irregular (nowhere differentiable). We need a mathematical construction to model such motion. x(t) x t Definition of Brownian motion Aim: model random movement of a particle (drunk man s walk). Consider a process W(t) with the following properties : - W(t) is continuous as a function of time - Increments on a time interval [s,t] are Gaussian variables with variance t-s W (t) W (s) t sn(0, 1) - Non-overlapping increments are independent cov (W (t) W (s),w(u) W (v)) = 0 0 t<s u<v T The resulting Brownian motion is nowhere differentiable, because W (t + t) W (t) t is a random variable with mean 0 and variance 1/ t
14 A different viewpoint: probability distributions Consider starting a Brownian motion at 0 : W (0) = 0 Then, at each point in time, the process evaluates to a Gaussian random variable with variance t and mean 0. W (t) = tn(0, 1) The probability of finding the particle at position x at time t is then 4 P(x W (t) <x+ dx) =ρ(x, t) = exp( x2 /2t) 2πt 2 x(t) t Expectation and Monte Carlo simulation Given a random variable X with a density ρ(x, t) We are interested in the expected value of a function of X (the average value over all possible Brownian paths). f(t) =f(x(t)) Numerically, this value can be approximated by taking many samples f(t) 1 M M f(x m (t)) m=1 Analytically, this corresponds to taking f(x(t)) = f(x)ρ(x, t)dx
15 Definition of stochastic differential equations (SDEs) A particle subject to pure Brownian motion, satisfies formally an SDE... dx(t) =dw (t), X(0) = X 0... which corresponds to saying that t X(t) =X 0 + dw (s) =X(0) + W (t) 0 We can add a drift term, which will make the expectation of X nonzero dx(t) =a(x) dt + bdw(t) This is a lot more to say about this type of equations, especially when the coefficients a and b depend on X. We will say some of these things in lecture 2... The Fokker-Planck equation A relation with a partial differential equation (PDE) Consider a stochastic differential equation dx(t) =a(x) dt + bdw(t) The probability distribution satisfies the Fokker-Planck equation ρ t + (a(x) ρ) x = b2 2 2 ρ x 2 This equation clearly shows that - the drift term in the SDE shows up as an advection term for the density - the Brownian motion in the SDE shows up as a diffusion term for the density
16 Example SDE and corresponding Fokker-Planck Heavy particle in a fluid The SDE for the heavy particle is (Langevin equation) dv = av dt + dw - Friction causes particle to slow down - Collisions with solvent molecules modeled as a Brownian motion of velocity The corresponding Fokker-Planck equation is ρ(v, t) t [avρ(v, t)] v = ρ(v, t) v 2 - Density is Gaussian at all times if initial condition is Gaussian v(t) ρ(v, t) t v Application 1 : Complex fluids Micro-macro simulation of dilute polymer solutions Macroscopic part : Navier-Stokes equation for the solvent ρ u t + u u = η s u p +div(τ p ) div(u) =0 Coupling : non-newtonian stress (Kramers expression) τ p = nx F (X) nk B T Id Microscopic part : SDE for dissolved polymers dx + u Xdt = ux 2ζ F (X) dt + σ 2 ζ dw t,
17 Concrete model problem: FENE dumbbells One-dimensional case / Known velocity / Lagrangian framework / FENE force dx = ux 2ζ F (X) dt + σ 2 X dw t, F(X) = ζ 1 X 2 /b Consider a time-stepper for an ensemble of such polymers X k+1 = s X (X k, u(t),δt), X R N The ensemble represents a distribution of polymer configurations δ = 0.50 δ = 1.00 Example: u(t) =2 ϕ( X ) X realizations - Polymer configuration distributions at two different moments in time Closure approximations Aim : Avoid computationally expensive micro/macro simulation Means : Derive a closed equation for a few macroscopic state variables Questions : U t U i = f i (X) + u U = H(U, u) τ p = T (U) - How to choose the macroscopic state variables? - How to determine the closure relations? Many closures have been proposed: FENE-L, FENE-S, FENE-P, FENE- LS, FENE-DS, { FENE-D }, quasi-equilibrium,...
18 Application 2 : Kinetic models Velocity-jump process for bacterial chemotaxis Each bacteria moves with a constant speed in a given direction. dx t dt = V t Direction changes at jump times, generated via a Poisson process. Tn+1 λ(t)dt = log(u n ), U n = U(0, 1) T n The modeling effort goes into the definition of the intensity of the Poisson process (turning rate). Turning rate for velocity jump processes Constant turning rate (leads to advection-diffusion) λ(t) λ 0 Turning rate biased by chemoattractant concentration λ(t) λ 0 g ( x S(X(t))) T V (t) - Direct sensing of chemoattractant gradient - Reduced turning rate if moving in favorable direction Turning rate biased by internal state (small bacteria) λ(t) =λ(y(t)), dy t dt = F (Y t,s(x t )) - Intensity of Poisson process depends on internal state - Evolution of internal state depends on the environment (memory! )
19 Example: cartoon internal dynamics Excitation/adaptation mechanism y Alignment with desired direction t λ(y) :=λ 0 by 2 λ(y) > λ 0 for y 2 < 0 λ Bacteria has a higher jump rate if it moves in an unfavorable direction. y 2 The nondimensionalized model Non-dimensionalize: t x : time scale associated with environmental changes t λ : typical time between two jumps t y : time scale associated with evolution of internal state := t λ /t x τ := t y /t λ dx t = V t, dt V t =1 V t = V n t [T n,t n+1 ) Tn+1 T n λ(y t )dt = θ n dy t dt = F (Y t,s(x t )) τ 1 <K
20 Kinetic description Equivalent continuum description of the same model in terms of the distribution p(x, v, y, t) Several ways to simulate this kinetic equation deterministically - direct discretization - computationally expensive due to high number of dimensions (x,v,y) stochastically t p + v x p + y (F (y, S(x))p) = Q(p) advection internal dynamics collisions - via an ensemble of individual velocity-jump processes - introduction of numerical noise (variance) due to finite number of particles reduction to an approximate macroscopic model Advection-diffusion limit for internal state model Define the density n(x, t) = Y In the limit 0, and scaling to long time scales t t/ 2, we get 1 t n = x 2 xn T (x)n The advection term depends on the details of the internal dynamics T (x) =b T V p(x, v, y, t)dvdy τ λ 0 τ + I xs(x) For nonzero, the difference of the density with respect to the limiting density is of order. Convergence is pathwise: trajectories converge to the corresponding SDE.
21 Model-scale separation : two types of problems Model problem 1 : approximation of Monte Carlo SDE simulation by moments - Microscopic description : Stochastic differential equation with fast evolution - Macroscopic level : Slowly evolving function of interest (expectation) Model problem 2 : approximation of kinetic equations on long time scales - Microscopic description: kinetic (phase-space) model - Macroscopic description: partial differential equation in space-time domain Later on in the course : - derivation of approximate macroscopic models in idealized situations (lecture 4) - numerical methods that exploit scale separation to simulate when deriving the macroscopic model analytically is not possible (second part of the course) Summary / Take-home message for today Problems with scale separation are everywhere. The multiscale nature can arise by having - multiple time scales - multiple spatial scales - multiple levels of description We have seen a number of applications and model problems for each. In the remainder of the course, we will look into - analytical techniques to obtain macroscopic models - conditions under which these can be derived - numerical techniques for when these conditions are not met
Introduction to multiscale modeling and simulation. Explicit methods for ODEs : forward Euler. y n+1 = y n + tf(y n ) dy dt = f(y), y(0) = y 0
Introduction to multiscale modeling and simulation Lecture 5 Numerical methods for ODEs, SDEs and PDEs The need for multiscale methods Two generic frameworks for multiscale computation Explicit methods
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 information2. Molecules in Motion
2. Molecules in Motion Kinetic Theory of Gases (microscopic viewpoint) assumptions (1) particles of mass m and diameter d; ceaseless random motion (2) dilute gas: d λ, λ = mean free path = average distance
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 informationDirect Modeling for Computational Fluid Dynamics
Direct Modeling for Computational Fluid Dynamics Kun Xu February 20, 2013 Computational fluid dynamics (CFD) is new emerging scientific discipline, and targets to simulate fluid motion in different scales.
More informationThe Polymers Tug Back
Tugging at Polymers in Turbulent Flow The Polymers Tug Back Jean-Luc Thiffeault http://plasma.ap.columbia.edu/ jeanluc Department of Applied Physics and Applied Mathematics Columbia University Tugging
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 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 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 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 1 Direct Modeling for Computational Fluid Dynamics
Chapter 1 Direct Modeling for Computational Fluid Dynamics Computational fluid dynamics (CFD) is a scientific discipline, which aims to capture fluid motion in a discretized space. The description of the
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
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 informationChapter 2. Brownian dynamics. 2.1 Introduction
Chapter 2 Brownian dynamics 2.1 Introduction The equations of physics are deterministic rather than stochastic in nature. Stochastic differential equations are used to approximate reality. They are introduced
More informationHomogenization with stochastic differential equations
Homogenization with stochastic differential equations Scott Hottovy shottovy@math.arizona.edu University of Arizona Program in Applied Mathematics October 12, 2011 Modeling with SDE Use SDE to model system
More informationTime-Dependent Statistical Mechanics 1. Introduction
Time-Dependent Statistical Mechanics 1. Introduction c Hans C. Andersen Announcements September 24, 2009 Lecture 1 9/22/09 1 Topics of concern in the course We shall be concerned with the time dependent
More informationTwo recent works on molecular systems out of equilibrium
Two recent works on molecular systems out of equilibrium Frédéric Legoll ENPC and INRIA joint work with M. Dobson, T. Lelièvre, G. Stoltz (ENPC and INRIA), A. Iacobucci and S. Olla (Dauphine). CECAM workshop:
More informationSystematic Closure Approximations for Multiscale Simulations
Systematic Closure Approximations for Multiscale Simulations Qiang Du Department of Mathematics/Materials Sciences Penn State University http://www.math.psu.edu/qdu Joint work with C. Liu, Y. Hyon and
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 informationNumerical methods in molecular dynamics and multiscale problems
Numerical methods in molecular dynamics and multiscale problems Two examples T. Lelièvre CERMICS - Ecole des Ponts ParisTech & MicMac project-team - INRIA Horizon Maths December 2012 Introduction The aim
More informationChapter 6 - Random Processes
EE385 Class Notes //04 John Stensby Chapter 6 - Random Processes Recall that a random variable X is a mapping between the sample space S and the extended real line R +. That is, X : S R +. A random process
More informationIntroduction to Aspects of Multiscale Modeling as Applied to Porous Media
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media Part III Todd Arbogast Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and
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 informationStochastic differential equation models in biology Susanne Ditlevsen
Stochastic differential equation models in biology Susanne Ditlevsen Introduction This chapter is concerned with continuous time processes, which are often modeled as a system of ordinary differential
More information16. Working with the Langevin and Fokker-Planck equations
16. Working with the Langevin and Fokker-Planck equations In the preceding Lecture, we have shown that given a Langevin equation (LE), it is possible to write down an equivalent Fokker-Planck equation
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 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 informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationExact Simulation of Diffusions and Jump Diffusions
Exact Simulation of Diffusions and Jump Diffusions A work by: Prof. Gareth O. Roberts Dr. Alexandros Beskos Dr. Omiros Papaspiliopoulos Dr. Bruno Casella 28 th May, 2008 Content 1 Exact Algorithm Construction
More informationNumerical methods for kinetic equations
Numerical methods for kinetic equations Lecture 6: fluid-kinetic coupling and hybrid methods Lorenzo Pareschi Department of Mathematics and Computer Science University of Ferrara, Italy http://www.lorenzopareschi.com
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 informationDiffusion of a density in a static fluid
Diffusion of a density in a static fluid u(x, y, z, t), density (M/L 3 ) of a substance (dye). Diffusion: motion of particles from places where the density is higher to places where it is lower, due to
More informationAnomalous Lévy diffusion: From the flight of an albatross to optical lattices. Eric Lutz Abteilung für Quantenphysik, Universität Ulm
Anomalous Lévy diffusion: From the flight of an albatross to optical lattices Eric Lutz Abteilung für Quantenphysik, Universität Ulm Outline 1 Lévy distributions Broad distributions Central limit theorem
More informationMonte Carlo methods for kinetic equations
Monte Carlo methods for kinetic equations Lecture 4: Hybrid methods and variance reduction Lorenzo Pareschi Department of Mathematics & CMCS University of Ferrara Italy http://utenti.unife.it/lorenzo.pareschi/
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 informationLecture 14. Turbulent Combustion. We know what a turbulent flow is, when we see it! it is characterized by disorder, vorticity and mixing.
Lecture 14 Turbulent Combustion 1 We know what a turbulent flow is, when we see it! it is characterized by disorder, vorticity and mixing. In a fluid flow, turbulence is characterized by fluctuations of
More informationAnomalous Transport and Fluctuation Relations: From Theory to Biology
Anomalous Transport and Fluctuation Relations: From Theory to Biology Aleksei V. Chechkin 1, Peter Dieterich 2, Rainer Klages 3 1 Institute for Theoretical Physics, Kharkov, Ukraine 2 Institute for Physiology,
More informationWind and turbulence experience strong gradients in vegetation. How do we deal with this? We have to predict wind and turbulence profiles through the
1 2 Wind and turbulence experience strong gradients in vegetation. How do we deal with this? We have to predict wind and turbulence profiles through the canopy. 3 Next we discuss turbulence in the canopy.
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 informationLecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Approximations
Lecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Approximations Simo Särkkä Aalto University, Finland November 18, 2014 Simo Särkkä (Aalto) Lecture 4: Numerical Solution of SDEs November
More informationBridging the Gap between Center and Tail for Multiscale Processes
Bridging the Gap between Center and Tail for Multiscale Processes Matthew R. Morse Department of Mathematics and Statistics Boston University BU-Keio 2016, August 16 Matthew R. Morse (BU) Moderate Deviations
More informationGillespie s Algorithm and its Approximations. Des Higham Department of Mathematics and Statistics University of Strathclyde
Gillespie s Algorithm and its Approximations Des Higham Department of Mathematics and Statistics University of Strathclyde djh@maths.strath.ac.uk The Three Lectures 1 Gillespie s algorithm and its relation
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 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 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 informationEffective dynamics for the (overdamped) Langevin equation
Effective dynamics for the (overdamped) Langevin equation Frédéric Legoll ENPC and INRIA joint work with T. Lelièvre (ENPC and INRIA) Enumath conference, MS Numerical methods for molecular dynamics EnuMath
More informationInverse Langevin approach to time-series data analysis
Inverse Langevin approach to time-series data analysis Fábio Macêdo Mendes Anníbal Dias Figueiredo Neto Universidade de Brasília Saratoga Springs, MaxEnt 2007 Outline 1 Motivation 2 Outline 1 Motivation
More information1. Introduction Systems evolving on widely separated time-scales represent a challenge for numerical simulations. A generic example is
COMM. MATH. SCI. Vol. 1, No. 2, pp. 385 391 c 23 International Press FAST COMMUNICATIONS NUMERICAL TECHNIQUES FOR MULTI-SCALE DYNAMICAL SYSTEMS WITH STOCHASTIC EFFECTS ERIC VANDEN-EIJNDEN Abstract. Numerical
More informationAn Overview of Fluid Animation. Christopher Batty March 11, 2014
An Overview of Fluid Animation Christopher Batty March 11, 2014 What distinguishes fluids? What distinguishes fluids? No preferred shape. Always flows when force is applied. Deforms to fit its container.
More informationLecture 4: Introduction to stochastic processes and stochastic calculus
Lecture 4: Introduction to stochastic processes and stochastic calculus Cédric Archambeau Centre for Computational Statistics and Machine Learning Department of Computer Science University College London
More informationAnisotropic fluid dynamics. Thomas Schaefer, North Carolina State University
Anisotropic fluid dynamics Thomas Schaefer, North Carolina State University Outline We wish to extract the properties of nearly perfect (low viscosity) fluids from experiments with trapped gases, colliding
More informationBefore we consider two canonical turbulent flows we need a general description of turbulence.
Chapter 2 Canonical Turbulent Flows Before we consider two canonical turbulent flows we need a general description of turbulence. 2.1 A Brief Introduction to Turbulence One way of looking at turbulent
More informationHomogenization Theory
Homogenization Theory Sabine Attinger Lecture: Homogenization Tuesday Wednesday Thursday August 15 August 16 August 17 Lecture Block 1 Motivation Basic Ideas Elliptic Equations Calculation of Effective
More informationIntroduction. Statistical physics: microscopic foundation of thermodynamics degrees of freedom 2 3 state variables!
Introduction Thermodynamics: phenomenological description of equilibrium bulk properties of matter in terms of only a few state variables and thermodynamical laws. Statistical physics: microscopic foundation
More informationStochastic differential equations in neuroscience
Stochastic differential equations in neuroscience Nils Berglund MAPMO, Orléans (CNRS, UMR 6628) http://www.univ-orleans.fr/mapmo/membres/berglund/ Barbara Gentz, Universität Bielefeld Damien Landon, MAPMO-Orléans
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 informationIntroduction to Marine Hydrodynamics
1896 1920 1987 2006 Introduction to Marine Hydrodynamics (NA235) Department of Naval Architecture and Ocean Engineering School of Naval Architecture, Ocean & Civil Engineering First Assignment The first
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationStochastic continuity equation and related processes
Stochastic continuity equation and related processes Gabriele Bassi c Armando Bazzani a Helmut Mais b Giorgio Turchetti a a Dept. of Physics Univ. of Bologna, INFN Sezione di Bologna, ITALY b DESY, Hamburg,
More informationDifferential Equations and Modeling
Differential Equations and Modeling Preliminary Lecture Notes Adolfo J. Rumbos c Draft date: March 22, 2018 March 22, 2018 2 Contents 1 Preface 5 2 Introduction to Modeling 7 2.1 Constructing Models.........................
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 informationTurbulence Modeling I!
Outline! Turbulence Modeling I! Grétar Tryggvason! Spring 2010! Why turbulence modeling! Reynolds Averaged Numerical Simulations! Zero and One equation models! Two equations models! Model predictions!
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 informationStochastic Modelling in Climate Science
Stochastic Modelling in Climate Science David Kelly Mathematics Department UNC Chapel Hill dtbkelly@gmail.com November 16, 2013 David Kelly (UNC) Stochastic Climate November 16, 2013 1 / 36 Why use stochastic
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 informationLANGEVIN EQUATION AND THERMODYNAMICS
LANGEVIN EQUATION AND THERMODYNAMICS RELATING STOCHASTIC DYNAMICS WITH THERMODYNAMIC LAWS November 10, 2017 1 / 20 MOTIVATION There are at least three levels of description of classical dynamics: thermodynamic,
More informationSDE Coefficients. March 4, 2008
SDE Coefficients March 4, 2008 The following is a summary of GARD sections 3.3 and 6., mainly as an overview of the two main approaches to creating a SDE model. Stochastic Differential Equations (SDE)
More informationBrownian Motion. 1 Definition Brownian Motion Wiener measure... 3
Brownian Motion Contents 1 Definition 2 1.1 Brownian Motion................................. 2 1.2 Wiener measure.................................. 3 2 Construction 4 2.1 Gaussian process.................................
More informationModeling of 1D Anomalous Diffusion In Fractured Nanoporous Media
LowPerm2015 Colorado School of Mines Low Permeability Media and Nanoporous Materials from Characterisation to Modelling: Can We Do It Better? IFPEN / Rueil-Malmaison - 9-11 June 2015 CSM Modeling of 1D
More informationLocal vs. Nonlocal Diffusions A Tale of Two Laplacians
Local vs. Nonlocal Diffusions A Tale of Two Laplacians Jinqiao Duan Dept of Applied Mathematics Illinois Institute of Technology Chicago duan@iit.edu Outline 1 Einstein & Wiener: The Local diffusion 2
More informationActive and Driven Soft Matter Lecture 3: Self-Propelled Hard Rods
A Tutorial: From Langevin equation to Active and Driven Soft Matter Lecture 3: Self-Propelled Hard Rods M. Cristina Marchetti Syracuse University Boulder School 2009 A Tutorial: From Langevin equation
More informationExplaining and modelling the rheology of polymeric fluids with the kinetic theory
Explaining and modelling the rheology of polymeric fluids with the kinetic theory Dmitry Shogin University of Stavanger The National IOR Centre of Norway IOR Norway 2016 Workshop April 25, 2016 Overview
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 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 informationDiffusion - The Heat Equation
Chapter 6 Diffusion - The Heat Equation 6.1 Goal Understand how to model a simple diffusion process and apply it to derive the heat equation in one dimension. We begin with the fundamental conservation
More informationFluid Animation. Christopher Batty November 17, 2011
Fluid Animation Christopher Batty November 17, 2011 What distinguishes fluids? What distinguishes fluids? No preferred shape Always flows when force is applied Deforms to fit its container Internal forces
More information2 Equations of Motion
2 Equations of Motion system. In this section, we will derive the six full equations of motion in a non-rotating, Cartesian coordinate 2.1 Six equations of motion (non-rotating, Cartesian coordinates)
More informationAn equation-free computational approach for extracting population-level behavior from individual-based models of biological dispersal
An equation-free computational approach for extracting population-level behavior from individual-based models of biological dispersal Radek Erban Ioannis G. Kevrekidis Hans G. Othmer November 6, 4 Abstract:
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 informationLecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Process Approximations
Lecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Process Approximations Simo Särkkä Aalto University Tampere University of Technology Lappeenranta University of Technology Finland November
More informationThe concentration of a drug in blood. Exponential decay. Different realizations. Exponential decay with noise. dc(t) dt.
The concentration of a drug in blood Exponential decay C12 concentration 2 4 6 8 1 C12 concentration 2 4 6 8 1 dc(t) dt = µc(t) C(t) = C()e µt 2 4 6 8 1 12 time in minutes 2 4 6 8 1 12 time in minutes
More informationModelling across different time and length scales
18/1/007 Course MP5 Lecture 1 18/1/007 Modelling across different time and length scales An introduction to multiscale modelling and mesoscale methods Dr James Elliott 0.1 Introduction 6 lectures by JAE,
More informationMultiscale, multiphysics modeling of turbulent transport and heating in collisionless, magnetized plasmas
Multiscale, multiphysics modeling of turbulent transport and heating in collisionless, magnetized plasmas Michael Barnes Plasma Science & Fusion Center Massachusetts Institute of Technology Collaborators:
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 informationIntroduction to nonequilibrium physics
Introduction to nonequilibrium physics Jae Dong Noh December 18, 2016 Preface This is a note for the lecture given in the 2016 KIAS-SNU Physics Winter Camp which is held at KIAS in December 17 23, 2016.
More informationMathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine
Lecture 2 The wave equation Mathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine V1.0 28/09/2018 1 Learning objectives of this lecture Understand the fundamental properties of the wave equation
More informationViscoelasticity. Basic Notions & Examples. Formalism for Linear Viscoelasticity. Simple Models & Mechanical Analogies. Non-linear behavior
Viscoelasticity Basic Notions & Examples Formalism for Linear Viscoelasticity Simple Models & Mechanical Analogies Non-linear behavior Viscoelastic Behavior Generic Viscoelasticity: exhibition of both
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationEuler equation and Navier-Stokes equation
Euler equation and Navier-Stokes equation WeiHan Hsiao a a Department of Physics, The University of Chicago E-mail: weihanhsiao@uchicago.edu ABSTRACT: This is the note prepared for the Kadanoff center
More informationMalliavin Calculus in Finance
Malliavin Calculus in Finance Peter K. Friz 1 Greeks and the logarithmic derivative trick Model an underlying assent by a Markov process with values in R m with dynamics described by the SDE dx t = b(x
More informationWeek 9 Generators, duality, change of measure
Week 9 Generators, duality, change of measure Jonathan Goodman November 18, 013 1 Generators This section describes a common abstract way to describe many of the differential equations related to Markov
More informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More information(Super) Fluid Dynamics. Thomas Schaefer, North Carolina State University
(Super) Fluid Dynamics Thomas Schaefer, North Carolina State University Hydrodynamics Hydrodynamics (undergraduate version): Newton s law for continuous, deformable media. Fluids: Gases, liquids, plasmas,...
More informationChapter 1 Introduction
Chapter 1 Introduction This thesis is concerned with the behaviour of polymers in flow. Both polymers in solutions and polymer melts will be discussed. The field of research that studies the flow behaviour
More informationNumerical treatment of stochastic delay differential equations
Numerical treatment of stochastic delay differential equations Evelyn Buckwar Department of Mathematics Heriot-Watt University, Edinburgh Evelyn Buckwar Roma, 16.10.2007 Overview An example: machine chatter
More informationSeveral forms of the equations of motion
Chapter 6 Several forms of the equations of motion 6.1 The Navier-Stokes equations Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationt = no of steps of length s
s t = no of steps of length s Figure : Schematic of the path of a diffusing molecule, for example, one in a gas or a liquid. The particle is moving in steps of length s. For a molecule in a liquid the
More informationWhy Complexity is Different
Why Complexity is Different Yaneer Bar-Yam (Dated: March 21, 2017) One of the hardest things to explain is why complex systems are actually different from simple systems. The problem is rooted in a set
More informationNONLOCALITY AND STOCHASTICITY TWO EMERGENT DIRECTIONS FOR APPLIED MATHEMATICS. Max Gunzburger
NONLOCALITY AND STOCHASTICITY TWO EMERGENT DIRECTIONS FOR APPLIED MATHEMATICS Max Gunzburger Department of Scientific Computing Florida State University North Carolina State University, March 10, 2011
More information5 Applying the Fokker-Planck equation
5 Applying the Fokker-Planck equation We begin with one-dimensional examples, keeping g = constant. Recall: the FPE for the Langevin equation with η(t 1 )η(t ) = κδ(t 1 t ) is = f(x) + g(x)η(t) t = x [f(x)p
More information