Introduction to asymptotic techniques for stochastic systems with multiple time-scales
|
|
- Vanessa Conley
- 5 years ago
- Views:
Transcription
1 Introduction to asymptotic techniques for stochastic systems with multiple time-scales Eric Vanden-Eijnden Courant Institute
2 Motivating examples Consider the ODE {Ẋ = Y 3 + sin(πt) + cos( 2πt) X() = x Ẏ = ε 1 (Y X) Y () = y. If ε 1, Y is very fast and it will adjust rapidly to the current value of X, i.e. we will have Y = X + O(ε) at all times. Then the equation for X reduces to (1) Ẋ = X 3 + sin(πt) + cos( 2πt). (2)
3 #)& # ")& "!)&!!!)&!"! " # $ % & ' ( The solution of (1) when ε =.5 and we took X t= = 2, Y t= = 1. X is shown in blue, and Y in green. Also shown in red is the solution of the limiting equation (2).
4 In contrast, consider the SDE {Ẋ = Y 3 + sin(πt) + cos( 2πt), X t= = x Ẏ = ε 1 (Y X) + ε 1/2 Ẇ, Y t= = y. The equation for Y at fixed X = x defines an Ornstein-Uhlenbeck process whose equilibrium probability density is ρ(y x) = e (y x)2 π. The limiting equation is obtained by averaging the right hand-side of the equation for X in (3) with respect to this density: (3) Ẋ = X X + sin(πt) + cos( 2πt), X = x. (4) Note the new term 3 2 α2 X due to the noise in (3).
5 % $ # "!!"!#!$! " # $ % & ' ( The solution of (3) with X t= = 2, Y t= = 1 when ε = 1 3. X is shown in blue, and Y in green. Also shown in red is the solution of the limiting equation (4). Notice how noisy Y is.
6 Singular perturbations techniques for Markov processes Consider the system {Ẋ = f(x, Y ), Xt= = x R n Ẏ = ε 1 g(x, Y ), Y t= = y R m. Assume that the equation for Y at fixed X = x defines an Markov process which is ergodic for every x with respect to the probability distribution dµ x (y) If F (x) = f(x, y)dµ x (y) exists R m then in the limit as ε the evolution for X solution of (5) is governed by (5) Ẋ = F (X) X() = x. (6)
7 Derivation. Let L be the infinitesimal generator of the Markov process generated by (5) 1( lim Ex,y φ(x t, Y t ) φ(x, y) ) = (Lφ)(x, y) (7) t + t Then L = L + ε 1 L 1, and u(x, y, t) = E x,y φ(y, X) satisfies the backward Kolmogorov equation u t = L u + ε 1 L 1 u, u t= = φ (8) Look for a solution in the form of so that lim ε u = u. u = u + εu 1 + O(ε 2 ) (9)
8 Inserting (9) into (8) and equating equal powers in ε leads to the hierarchy of equations L 1 u =, L 1 u 1 = u t L u, L 1 u 2 = (1) The first equation tells that u belong to the null-space of L 1. Assuming that for every x, L 1 is the generator of an ergodic Markov process with equilibrium distribution µ x (y), this null-space is spanned by functions constant in y, i.e. u = u (x, t). Since the null-space of L 1 is non-trivial, the next equations each requires a solvability condition, namely that their right hand-side belongs to the range of L 1.
9 To see what this solvability condition actually is, take the expectation of both sides of the second equation in (1) with respect to dµ x (x). This gives ( u ) = dµ x (y) R t L u (11) m Explicitly, (11) is where u t = F (x) xu (12) F (x) = R m f(x, y)dµ x (y) (13) (12) is the backward Kolmogorov equation of (6).
10 Remark: Computing the expectation with respect to µ x (y) in practice. We have (e L1t φ)(x, y) φ(x, y)dµ x (y) as t R m In other words, if u(x, y, t) satisfies u t = L 1u, u t= = φ so that formally u(x, y, t) = (e L1t φ)(x, y), we have lim u(x, y, t) = φ(x, y)dµ x (y) t R m But since u(x, y, t) = E y φ(y x ) where it follows that Ẏ x = g(x, Y x ) t ) φ(x, y)dµ x (y) = lim E y φ(y x R t m 1 = lim T T T φ(y x t )dt (14)
11 Example: the Lorenz 96 (L96) model L96 consists of K slow variables X k coupled to J K fast variables Y j,k whose evolution is governed by Ẋ k = X k 1 (X k 2 X k+1 ) X k + F x + h J x Y j,k J j=1 (15) Ẏ j,k = 1 ( ) Yj+1,k (Y j+2,k Y j 1,k ) Y j,k + h y X k. ε We will study (15) with F x = 1, h x =.8, h y = 1, K = 9, J = 8, and two values of ε: ε = 1/128 and ε = 1/124.
12 Typical time-series of the slow (black line) and fast (grey line) modes; K = 9, J = 8, ε = 1/128. The subplot displays a typical snapshot of the slow and fast modes at a given time.
13 PDF of the slow variable; K = 9, J = 8, black line: ε = 1/16, grey line: ε = 1/124. The insensitivity in ε of the PDFs indicates that the slow variables have already converged close to their limiting behavior when ε = 1/16.
14 ACFs of the slow (thick line) and fast (thin line) variables; ε = 1/128, grey line: ε = 1/124. The subplot is the zoom-in of the main graph which shows the transient decay of the ACFs of the fast modes becoming faster as ε is decreased: this is the only signature in the ACFs of the fact that the Y j,k s are faster.
15 Typical PDFs of the coupling term (h x /J) J j=1 Z j,k(x) for various values of x. These PDFs are robust against variations in the initial conditions indicating that the dynamics of the fast modes conditional on the slow ones being fixed is ergodic.
16 Comparison between the ACFs and PDFs; black line: limiting dynamics; full grey line: ε = 1/128. Also shown in dashed grey are the corresponding ACF and PDF produced by the truncated dynamics where the coupling of the slow modes X k with the fast ones, Y j,k, is artificially switched off.
17 Black points: scatterplot of the forcing F (x) in the limiting dynamics. Grey points: scatter plot of the bare coupling term (h x /J) J j=1 Y j,k(x) when ε = 1/128.
18 Diffusive time-scales Suppose that F (x) = f(x, y)dµ x (y) =. (16) R m Then the limiting equation on the O(1) time-scale is trivial, Ẋ =, and the interesting dynamics arises on the diffusive time-scale O(ε 1 ). Consider then {Ẋ = ε 1 f(x, Y ), X t= = x R n Ẏ = ε 2 g(x, Y ), Y t= = y R m. (17) Proceeding as before we arrive at the following limiting equation when ε : where and Ẋ = b(x) + σ(x) Ẇ t, (18) ( Lφ)(x) b(x) x φ(x) (σσt )(x) : x x φ(x) ( = dt dµ x (y)f(x, y) x Ey f(x, Yt x ) x φ(x) ) (19) R m Ẏ x = g(x, Y x )
19 Derivation. The backward Kolmogorov equation for u(x, y, t) = E x,y f(x) is now u t = ε 1 L u + ε 2 L 1 u. Inserting the expansion u = u + εu 1 + ε 2 u 2 + O(ε 2 ) (we will have to go one order in ε higher than before) in this equation now gives L 1 u =, L 1 u 1 = L u, L 1 u 2 = u t L u 1, L 1 u 3 = The first equation tells that u (x, y, t) = u (x, t). (2) The solvability condition for the second equation is satisfied by assumption because of (16). Therefore this equation can be formally solved as u 1 = L 1 1 L 2u. Inserting this expression in the third equation in (2) and considering the solvability condition for this equation, we obtain the limiting equation for u : u t = Lu, (21) where L = dµ x (y)l L 1 1 L. R m
20 To see what this equation is explicitly, notice that L 1 1 g(y) is the steady state solution of v t = L 1v + g(y). The solution of this equation with the initial condition v(y, ) = can be represented by Feynman-Kac formula as Therefore L 1 1 g(y) = E y g(yt x )dt, L in the limiting backward Kolmogorov equa- and the operator tion (21) is (19). v(y, t) = E y t g(y x s )ds,
21 Example: the Lorenz 96 (L96) model Consider Ẋ k = ε (X k 1 (X k 2 X k+1 ) + X k ) + h J x ( ) Yj,k+1 Y j,k 1 J j=1 Ẏ j,k = 1 ( ) Yj+1,k (Y j+2,k Y j 1,k ) Y j,k + F y + hy X k, ε (22)
22 ACFs and PDFs (subplot) of the slow variable X k evolving under (22). Grey line: ε = 1/256; full black line: ε = 1/128; dashed black line: ε = 1/128 with time rescaled as t 2t. The near perfect match confirms that evolution of X k converges to some limiting dynamics on the O(1/ε) time-scale.
23 Coupled truncated Burgers-Hopf (TBH). I. Stable periodic orbits Ẋ 1 = 1 ε b 1X 2 Y 1 + ay 1 (R 2 (X1 2 + X2)) 2 bx 2 (α + (X1 2 + X2)), 2 Ẋ 2 = 1 ε b 2X 1 Y 1 + ax 2 (R 2 (X1 2 + X2)) 2 + bx 1 (α + (X1 2 + X2)), 2 Ẏ k = Re ik Up Uq ε b 3δ 1,k X 1 X 2, p+q+k= Ż k = Im ik Up Uq, 2 p+q+k= where U k = Y k + iz k. Truncated system: stable periodic orbit (limit cycle) (X 1 (t), X 2 (t)) = R (cos ωt, sin ωt) with frequency ω = b(α + R 2 ).
24 NB: Truncated Burgers (Majda & Timofeyev, 2) Fourier-Galerkin truncation of the inviscid Burgers-Hopf equation, u t (u2 ) x = : U k = ik 2 U p U q, k < Λ k+p+q= p, q Λ Features common with many complex systems. In particular: Display deterministic chaos. Ergodic on E = k U k 2. Scaling law for the correlation functions with t k O(k 1 ). Here: Used as a model for unresolved modes. Couple truncated Burgers with two resolved variables.
25 Ẋ 1 = 1 ε b 1X 2 Y 1 + ay 1 (R 2 (X1 2 + X2 2 )) bx 2(α + (X1 2 + X2 2 )), Ẋ 2 = 1 ε b 2X 1 Y 1 + ax 2 (R 2 (X1 2 + X2)) 2 + bx 1 (α + (X1 2 + X2)), 2 Ẏ k = Re ik Up Uq ε b 3δ 1,k X 1 X 2, p+q+k= Ż k = Im ik 2 p+q+k= U p U q, Limiting SDEs: Ẋ 1 = b 1 b 2 X 1 + N 1 X 1 X σ 1X 2 Ẇ (t) + ax 1 (R 2 (X X2 2 )) bx 2(α + (X X2 2 )), Ẋ 2 = b 1 b 2 X 2 + N 2 X 2 X σ 2X 1 Ẇ (t) + ax 2 (R 2 (X X2 2 )) + bx 1(α + X X2 2 )),
26 (a) (b) 1 1 x x Contour plots of the joint probability density for the climate variables X 1 and X 2 ; (a) deterministic system with 12 variables; (b) limit SDE. There only remains a ghost of the limit cycle.
27 3 2.5 (a) DNS (b) Reduced model PDF of x PDF of x x x (c) DNS (d) Reduced model PDF of x PDF of x x x 2 Marginal PDFs of X 1 and X 2 for the simulations of the full deterministic system with 12 variables (DNS) and the limit SDE.
28 corr of x 1 (a) t 8 (c) corr of x 2 (b) t 8 (d) cross corr of x 1, x t 8 K 2 (t) t 8 Two-point statistics for X 1 and X 2 ; solid lines - deterministic system with 12 degrees of freedom; dashed lines - limit SDE; (a), (b) correlation functions of X 1 and X 2, respectively; (c) cross-correlation function of X 1 and X 2 ; (d) normalized correlation of energy, K 2 (t) = X 2 2 (t)x2 2 () X X 2 (t)x 2 () 2
29 Coupled truncated Burgers-Hopf (TBH). II. Multiple equilibria Ẋ = 1 ε b 1Y 1 Z 1 + λ(1 αx 2 )x, Ẏ k = Re ik 2 p+q+k= U p U q + 1 ε b 2δ 1,k XZ k, Ż k = Im ik 2 p+q+k= u pu q + 1 ε b 3δ 1,k XY k, Limiting SDE: Ẋ = γx + σẇ (t) + λ(1 αx 2 )X
30 .6 (a) (b) (c).6.6 λ = 1.2 λ =.5 λ=.15 PDF of x x PDF of X for the simulations of the deterministic model (solid lines) and the limit SDE (dashed lines) in three regimes, λ = 1.2,.5,.15.
31 1 (a) 1 (b) 1 (c) corr of x.5 λ =1.2.5 λ =.5.5 λ = t (d) 1.2 (e) 1.2 (f) K(t) λ =1.2 λ = λ = t Two-point statistics of X in three regimes, λ = 1.2,.5,.15 for the deterministic equations (solid lines) and limit SDE (dashed lines) (a), (b), (c) correlation function of x; (d), (e), (f) correlation of energy, K(t), in x.
32 Generalizations Let Z t S be the sample path of a continuous-time Markov process with generator L = L + ε 1 L 1 and assume that L 1 has several ergodic components indexed by x S with equilibrium distribution µ x (z) Then, as ε there exists a limiting on S with generator L = E µx L Similar results available on diffusive time-scales. Can be applied to SDEs, Markov chains (i.e. discrete state-space), deterministic systems (periodic or chaotic), etc.
33 Consider Other situations with limiting dynamics? {Ẋt = f(x t, Y t ), Ẏ t = g(x t, Y t ), (23) and denote by ϕ(x [,t] ) the solution of the equation for Y at time t assuming that X t is known on [, t] Observe that ϕ(x [,t] ) is a functional of {X s, s [, t]}. Then X t satisfies the closed equation Ẋ t = f(x t, ϕ(x [,t] )) (24) In general, X t is not Markov! It becomes Markov when Y t is faster, or...?
34 Other possibility: Weak coupling The system Ẋ t = 1 N Ẏ n t N n=1 = g(x t, Yt n), f(x t, Y n t ), (all the Y n t coupled only via X t ) may have a limit behavior as N which is the same as the limit behavior as Ẋ t = f(x t, Y t ), when ε. Ẏ t = 1 ε g(x t, Y t ),
35 Asymptotic techniques for singularly perturbed Markov processes: R. Z. Khasminsky, On Stochastic Processes Defined by Differential Equations with a Small Parameter, Theory Prob. Applications, 11: , R. Z. Khasminsky, A Limit Theorem for the Solutions of Differential Equations with Random Right-Hand Sides, Theory Prob. Applications, 11:39 46, T. G. Kurtz, Semigroups of conditioned shifts and approximations of Markov processes, Annals of Probability, 3: , G. Papanicolaou, Some probabilistic problems and methods in singular perturbations, Rocky Mountain J. Math, 6: , M. I. Freidlin and A. D. Wentzell, Random perturbations of dynamical systems, 2nd edition, Springer-Verlag, Effective dynamics in L96: E. N. Lorenz, Predictability A problem partly solved, pp in: ECMWF Seminar Proceedings on Predictability, Reading, United Kingdom, ECMWF, C. Rödenbeck, C. Beck, H. Kantz, Dynamical systems with time scale-separation: averaging, stochastic modeling, and central limit theorems, pp in: Stochastic Climate Models, P. Imkeller, J.-S. von Storch eds., Progress in Probability 49, Birkhäuser Verlag, Basel, 21. I. Fatkullin and E. Vanden-Eijnden, A computational strategy for multiscale systems with applications to Lorenz 96 model, J. Comp. Phys. 2:65 638, 24.
36 (Coupled) TBH: A. J. Majda and I. Timofeyev, Remarkable statistical behavior for truncated Burgers- Hopf dynamics, Proc. Nat. Acad. Sci. USA, 97: , 2. A. J. Majda and I. Timofeyev, Statistical mechanics for truncations of the Burgers- Hopf equation: a model for intrinsic stochastic behavior with scaling, Milan Journal of Mathematics, 7(1):39 96, 22. A. J. Majda, I. Timofeyev, and E. Vanden-Eijnden, Systematic strategies for stochastic mode reduction in climate. J. Atmos. Sci., 6(14): , 23. A. J. Majda, I. Timofeyev, and E. Vanden-Eijnden, Stochastic models for selected slow variables in large deterministic systems, Nonlinearity, 19: , 26.
Stochastic Mode Reduction in Large Deterministic Systems
Stochastic Mode Reduction in Large Deterministic Systems I. Timofeyev University of Houston With A. J. Majda and E. Vanden-Eijnden; NYU Related Publications: http://www.math.uh.edu/ ilya Plan Mode-Elimination
More informationLow-dimensional chaotic dynamics versus intrinsic stochastic noise: a paradigm model
Physica D 199 (2004) 339 368 Low-dimensional chaotic dynamics versus intrinsic stochastic noise: a paradigm model A. Majda a,b, I. Timofeyev c, a Courant Institute of Mathematical Sciences, New York University,
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 informationStochastic models for selected slow variables in large deterministic systems
INSTITUTE OF PHYSICS PUBLISHING Nonlinearity 19 (26) 769 794 NONLINEARITY doi:1.188/951-7715/19/4/1 Stochastic models for selected slow variables in large deterministic systems A Majda 1,2, I Timofeyev
More informationA computational strategy for multiscale chaotic systems with applications to Lorenz 96 model
A computational strategy for multiscale chaotic systems with applications to Lorenz 96 model Ibrahim Fatkullin a, Eric Vanden-Eijnden b,1 a California Institute of Technology, Pasadena, CA 91125, USA b
More informationA SIMPLE LINEAR RESPONSE CLOSURE APPROXIMATION FOR SLOW DYNAMICS OF A MULTISCALE SYSTEM WITH LINEAR COUPLING
A SIMPLE LINEAR RESPONSE CLOSURE APPROXIMATION FOR SLOW DYNAMICS OF A MULTISCALE SYSTEM WITH LINEAR COUPLING RAFAIL V. ABRAMOV Abstract. Many applications of contemporary science involve multiscale dynamics,
More informationIntroduction 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 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 informationStochastic Volatility and Correction to the Heat Equation
Stochastic Volatility and Correction to the Heat Equation Jean-Pierre Fouque, George Papanicolaou and Ronnie Sircar Abstract. From a probabilist s point of view the Twentieth Century has been a century
More informationLecture 12: Detailed balance and Eigenfunction methods
Miranda Holmes-Cerfon Applied Stochastic Analysis, Spring 2015 Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility),
More informationFast-slow systems with chaotic noise
Fast-slow systems with chaotic noise David Kelly Ian Melbourne Courant Institute New York University New York NY www.dtbkelly.com May 1, 216 Statistical properties of dynamical systems, ESI Vienna. David
More informationMathematical Theories of Turbulent Dynamical Systems
Mathematical Theories of Turbulent Dynamical Systems Di Qi, and Andrew J. Majda Courant Institute of Mathematical Sciences Fall 2016 Advanced Topics in Applied Math Di Qi, and Andrew J. Majda (CIMS) Mathematical
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 informationLYAPUNOV EXPONENTS AND STABILITY FOR THE STOCHASTIC DUFFING-VAN DER POL OSCILLATOR
LYAPUNOV EXPONENTS AND STABILITY FOR THE STOCHASTIC DUFFING-VAN DER POL OSCILLATOR Peter H. Baxendale Department of Mathematics University of Southern California Los Angeles, CA 90089-3 USA baxendal@math.usc.edu
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationHandbook of Stochastic Methods
Springer Series in Synergetics 13 Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences von Crispin W Gardiner Neuausgabe Handbook of Stochastic Methods Gardiner schnell und portofrei
More informationHandbook of Stochastic Methods
C. W. Gardiner Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences Third Edition With 30 Figures Springer Contents 1. A Historical Introduction 1 1.1 Motivation I 1.2 Some Historical
More informationKolmogorov Equations and Markov Processes
Kolmogorov Equations and Markov Processes May 3, 013 1 Transition measures and functions Consider a stochastic process {X(t)} t 0 whose state space is a product of intervals contained in R n. We define
More informationA priori tests of a stochastic mode reduction strategy
Physica D 170 (2002) 206 252 A priori tests of a stochastic mode reduction strategy A. Majda a,b, I. Timofeyev a,b, E. Vanden-Eijnden a, a Courant Institute of Mathematical Sciences, New Yor University,
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 informationLecture 12: Detailed balance and Eigenfunction methods
Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility), 4.2-4.4 (explicit examples of eigenfunction methods) Gardiner
More informationEdward Lorenz: Predictability
Edward Lorenz: Predictability Master Literature Seminar, speaker: Josef Schröttle Edward Lorenz in 1994, Northern Hemisphere, Lorenz Attractor I) Lorenz, E.N.: Deterministic Nonperiodic Flow (JAS, 1963)
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 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 informationFast-slow systems with chaotic noise
Fast-slow systems with chaotic noise David Kelly Ian Melbourne Courant Institute New York University New York NY www.dtbkelly.com May 12, 215 Averaging and homogenization workshop, Luminy. Fast-slow systems
More informationSPDEs, criticality, and renormalisation
SPDEs, criticality, and renormalisation Hendrik Weber Mathematics Institute University of Warwick Potsdam, 06.11.2013 An interesting model from Physics I Ising model Spin configurations: Energy: Inverse
More informationPhysics Constrained Nonlinear Regression Models for Time Series
Physics Constrained Nonlinear Regression Models for Time Series Andrew J. Majda and John Harlim 2 Department of Mathematics and Center for Atmosphere and Ocean Science, Courant Institute of Mathematical
More informationPath integrals for classical Markov processes
Path integrals for classical Markov processes Hugo Touchette National Institute for Theoretical Physics (NITheP) Stellenbosch, South Africa Chris Engelbrecht Summer School on Non-Linear Phenomena in Field
More informationBIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs
BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits. 2. Equilibria. 3. Periodic orbits and limit cycles. 4. Homoclinic orbits.
More informationA Mathematical Framework for Stochastic Climate Models
A Mathematical Framewor for Stochastic Climate Models ANDREW J. MAJDA ILYA TIMOFEYEV AND ERIC VANDEN EIJNDEN Courant Institute Abstract There has been a recent burst of activity in the atmosphere-ocean
More informationIntroduction to Stochastic SIR Model
Introduction to Stochastic R Model Chiu- Yu Yang (Alex), Yi Yang R model is used to model the infection of diseases. It is short for Susceptible- Infected- Recovered. It is important to address that R
More informationDATA-DRIVEN TECHNIQUES FOR ESTIMATION AND STOCHASTIC REDUCTION OF MULTI-SCALE SYSTEMS
DATA-DRIVEN TECHNIQUES FOR ESTIMATION AND STOCHASTIC REDUCTION OF MULTI-SCALE SYSTEMS A Dissertation Presented to the Faculty of the Department of Mathematics University of Houston In Partial Fulfillment
More informationRandom and Deterministic perturbations of dynamical systems. Leonid Koralov
Random and Deterministic perturbations of dynamical systems Leonid Koralov - M. Freidlin, L. Koralov Metastability for Nonlinear Random Perturbations of Dynamical Systems, Stochastic Processes and Applications
More informationContinuous dependence estimates for the ergodic problem with an application to homogenization
Continuous dependence estimates for the ergodic problem with an application to homogenization Claudio Marchi Bayreuth, September 12 th, 2013 C. Marchi (Università di Padova) Continuous dependence Bayreuth,
More informationStatistical Mechanics for the Truncated Quasi-Geostrophic Equations
Statistical Mechanics for the Truncated Quasi-Geostrophic Equations Di Qi, and Andrew J. Majda Courant Institute of Mathematical Sciences Fall 6 Advanced Topics in Applied Math Di Qi, and Andrew J. Majda
More informationMathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )
Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio (2014-2015) Etienne Tanré - Olivier Faugeras INRIA - Team Tosca November 26th, 2014 E. Tanré (INRIA - Team Tosca) Mathematical
More information1 Elementary probability
1 Elementary probability Problem 1.1 (*) A coin is thrown several times. Find the probability, that at the n-th experiment: (a) Head appears for the first time (b) Head and Tail have appeared equal number
More informationNested Stochastic Simulation Algorithm for KMC with Multiple Time-Scales
Nested Stochastic Simulation Algorithm for KMC with Multiple Time-Scales Eric Vanden-Eijnden Courant Institute Join work with Weinan E and Di Liu Chemical Systems: Evolution of an isothermal, spatially
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 informationUncertainty quantification and systemic risk
Uncertainty quantification and systemic risk Josselin Garnier (Université Paris Diderot) with George Papanicolaou and Tzu-Wei Yang (Stanford University) February 3, 2016 Modeling systemic risk We consider
More informationA MODEL FOR THE LONG-TERM OPTIMAL CAPACITY LEVEL OF AN INVESTMENT PROJECT
A MODEL FOR HE LONG-ERM OPIMAL CAPACIY LEVEL OF AN INVESMEN PROJEC ARNE LØKKA AND MIHAIL ZERVOS Abstract. We consider an investment project that produces a single commodity. he project s operation yields
More informationProf. Krstic Nonlinear Systems MAE281A Homework set 1 Linearization & phase portrait
Prof. Krstic Nonlinear Systems MAE28A Homework set Linearization & phase portrait. For each of the following systems, find all equilibrium points and determine the type of each isolated equilibrium. Use
More informationReconstruction Deconstruction:
Reconstruction Deconstruction: A Brief History of Building Models of Nonlinear Dynamical Systems Jim Crutchfield Center for Computational Science & Engineering Physics Department University of California,
More informationReaction-Diffusion Equations In Narrow Tubes and Wave Front P
Outlines Reaction-Diffusion Equations In Narrow Tubes and Wave Front Propagation University of Maryland, College Park USA Outline of Part I Outlines Real Life Examples Description of the Problem and Main
More information1 Brownian Local Time
1 Brownian Local Time We first begin by defining the space and variables for Brownian local time. Let W t be a standard 1-D Wiener process. We know that for the set, {t : W t = } P (µ{t : W t = } = ) =
More informationAPPPHYS217 Tuesday 25 May 2010
APPPHYS7 Tuesday 5 May Our aim today is to take a brief tour of some topics in nonlinear dynamics. Some good references include: [Perko] Lawrence Perko Differential Equations and Dynamical Systems (Springer-Verlag
More informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
More informationNested Stochastic Simulation Algorithm for Chemical Kinetic Systems with Disparate Rates. Abstract
Nested Stochastic Simulation Algorithm for Chemical Kinetic Systems with Disparate Rates Weinan E Department of Mathematics and PACM, Princeton University, Princeton, NJ 08544, USA Di Liu and Eric Vanden-Eijnden
More informationIntertwinings for Markov processes
Intertwinings for Markov processes Aldéric Joulin - University of Toulouse Joint work with : Michel Bonnefont - Univ. Bordeaux Workshop 2 Piecewise Deterministic Markov Processes ennes - May 15-17, 2013
More informationRandom Dynamical Systems with Non-Gaussian Fluctuations Some Advances & Perspectives
Random Dynamical Systems with Non-Gaussian Fluctuations Some Advances & Perspectives Jinqiao Duan Illinois Institute of Technology Department of Applied Mathematics Chicago, IL 60616, USA www.iit.edu/~duan
More informationPart II. Dynamical Systems. Year
Part II Year 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2017 34 Paper 1, Section II 30A Consider the dynamical system where β > 1 is a constant. ẋ = x + x 3 + βxy 2, ẏ = y + βx 2
More informationStability of Stochastic Differential Equations
Lyapunov stability theory for ODEs s Stability of Stochastic Differential Equations Part 1: Introduction Department of Mathematics and Statistics University of Strathclyde Glasgow, G1 1XH December 2010
More informationElliptic Operators with Unbounded Coefficients
Elliptic Operators with Unbounded Coefficients Federica Gregorio Universitá degli Studi di Salerno 8th June 2018 joint work with S.E. Boutiah, A. Rhandi, C. Tacelli Motivation Consider the Stochastic Differential
More informationOne dimensional Maps
Chapter 4 One dimensional Maps The ordinary differential equation studied in chapters 1-3 provide a close link to actual physical systems it is easy to believe these equations provide at least an approximate
More informationENGI 9420 Lecture Notes 8 - PDEs Page 8.01
ENGI 940 Lecture Notes 8 - PDEs Page 8.01 8. Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives
More information1. (i) Determine how many periodic orbits and equilibria can be born at the bifurcations of the zero equilibrium of the following system:
1. (i) Determine how many periodic orbits and equilibria can be born at the bifurcations of the zero equilibrium of the following system: ẋ = y x 2, ẏ = z + xy, ż = y z + x 2 xy + y 2 + z 2 x 4. (ii) Determine
More informationDynamical systems with Gaussian and Levy noise: analytical and stochastic approaches
Dynamical systems with Gaussian and Levy noise: analytical and stochastic approaches Noise is often considered as some disturbing component of the system. In particular physical situations, noise becomes
More informationData assimilation in high dimensions
Data assimilation in high dimensions David Kelly Courant Institute New York University New York NY www.dtbkelly.com February 12, 2015 Graduate seminar, CIMS David Kelly (CIMS) Data assimilation February
More informationInterest Rate Models:
1/17 Interest Rate Models: from Parametric Statistics to Infinite Dimensional Stochastic Analysis René Carmona Bendheim Center for Finance ORFE & PACM, Princeton University email: rcarmna@princeton.edu
More informationStrong uniqueness for stochastic evolution equations with possibly unbounded measurable drift term
1 Strong uniqueness for stochastic evolution equations with possibly unbounded measurable drift term Enrico Priola Torino (Italy) Joint work with G. Da Prato, F. Flandoli and M. Röckner Stochastic Processes
More informationA Study of the Van der Pol Equation
A Study of the Van der Pol Equation Kai Zhe Tan, s1465711 September 16, 2016 Abstract The Van der Pol equation is famous for modelling biological systems as well as being a good model to study its multiple
More informationAnti-synchronization of a new hyperchaotic system via small-gain theorem
Anti-synchronization of a new hyperchaotic system via small-gain theorem Xiao Jian( ) College of Mathematics and Statistics, Chongqing University, Chongqing 400044, China (Received 8 February 2010; revised
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 informationSuperparameterization and Dynamic Stochastic Superresolution (DSS) for Filtering Sparse Geophysical Flows
Superparameterization and Dynamic Stochastic Superresolution (DSS) for Filtering Sparse Geophysical Flows June 2013 Outline 1 Filtering Filtering: obtaining the best statistical estimation of a nature
More informationMARKOV CHAIN AND TIME-DELAY REDUCED MODELING OF NONLINEAR SYSTEMS
MARKOV CHAIN AND TIME-DELAY REDUCED MODELING OF NONLINEAR SYSTEMS A Dissertation Presented to the Faculty of the Department of Mathematics University of Houston In Partial Fulfillment of the Requirements
More informationENGI 9420 Lecture Notes 8 - PDEs Page 8.01
ENGI 940 ecture Notes 8 - PDEs Page 8.0 8. Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives
More informationData Assimilation in Slow Fast Systems Using Homogenized Climate Models
APRIL 2012 M I T C H E L L A N D G O T T W A L D 1359 Data Assimilation in Slow Fast Systems Using Homogenized Climate Models LEWIS MITCHELL AND GEORG A. GOTTWALD School of Mathematics and Statistics,
More informationStochastic analysis of biochemical reaction networks with absolute concentration robustness
Stochastic analysis of biochemical reaction networks with absolute concentration robustness David F. Anderson anderson@math.wisc.edu Department of Mathematics University of Wisconsin - Madison Boston University
More informationPerformance Evaluation of Generalized Polynomial Chaos
Performance Evaluation of Generalized Polynomial Chaos Dongbin Xiu, Didier Lucor, C.-H. Su, and George Em Karniadakis 1 Division of Applied Mathematics, Brown University, Providence, RI 02912, USA, gk@dam.brown.edu
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 informationMATH 425, FINAL EXAM SOLUTIONS
MATH 425, FINAL EXAM SOLUTIONS Each exercise is worth 50 points. Exercise. a The operator L is defined on smooth functions of (x, y by: Is the operator L linear? Prove your answer. L (u := arctan(xy u
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 information2012 NCTS Workshop on Dynamical Systems
Barbara Gentz gentz@math.uni-bielefeld.de http://www.math.uni-bielefeld.de/ gentz 2012 NCTS Workshop on Dynamical Systems National Center for Theoretical Sciences, National Tsing-Hua University Hsinchu,
More information11 Chaos in Continuous Dynamical Systems.
11 CHAOS IN CONTINUOUS DYNAMICAL SYSTEMS. 47 11 Chaos in Continuous Dynamical Systems. Let s consider a system of differential equations given by where x(t) : R R and f : R R. ẋ = f(x), The linearization
More informationSUB-SAMPLING AND PARAMETRIC ESTIMATION FOR MULTISCALE DYNAMICS
SUB-SAMPLING AND PARAMETRIC ESTIMATION FOR MULTISCALE DYNAMICS R. AZENCOTT, A. BERI, A. JAIN, AND I. TIMOFEYEV Abstract. We study the problem of adequate data sub-sampling for consistent parametric estimation
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationStochastic contraction BACS Workshop Chamonix, January 14, 2008
Stochastic contraction BACS Workshop Chamonix, January 14, 2008 Q.-C. Pham N. Tabareau J.-J. Slotine Q.-C. Pham, N. Tabareau, J.-J. Slotine () Stochastic contraction 1 / 19 Why stochastic contraction?
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 informationDiffusion for a Markov, Divergence-form Generator
Diffusion for a Markov, Divergence-form Generator Clark Department of Mathematics Michigan State University Arizona School of Analysis and Mathematical Physics March 16, 2012 Abstract We consider the long-time
More informationEntrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4.
Entrance Exam, Differential Equations April, 7 (Solve exactly 6 out of the 8 problems). Consider the following initial value problem: { y + y + y cos(x y) =, y() = y. Find all the values y such that the
More informationTHREE DIMENSIONAL SYSTEMS. Lecture 6: The Lorenz Equations
THREE DIMENSIONAL SYSTEMS Lecture 6: The Lorenz Equations 6. The Lorenz (1963) Equations The Lorenz equations were originally derived by Saltzman (1962) as a minimalist model of thermal convection in a
More informationFiltering Sparse Regular Observed Linear and Nonlinear Turbulent System
Filtering Sparse Regular Observed Linear and Nonlinear Turbulent System John Harlim, Andrew J. Majda Department of Mathematics and Center for Atmosphere Ocean Science Courant Institute of Mathematical
More informationIntroduction to Dynamical Systems Basic Concepts of Dynamics
Introduction to Dynamical Systems Basic Concepts of Dynamics A dynamical system: Has a notion of state, which contains all the information upon which the dynamical system acts. A simple set of deterministic
More informationEffective dynamics of many-body quantum systems
Effective dynamics of many-body quantum systems László Erdős University of Munich Grenoble, May 30, 2006 A l occassion de soixantiéme anniversaire de Yves Colin de Verdiére Joint with B. Schlein and H.-T.
More informationLecture 6. Lorenz equations and Malkus' waterwheel Some properties of the Lorenz Eq.'s Lorenz Map Towards definitions of:
Lecture 6 Chaos Lorenz equations and Malkus' waterwheel Some properties of the Lorenz Eq.'s Lorenz Map Towards definitions of: Chaos, Attractors and strange attractors Transient chaos Lorenz Equations
More informationHigh-Gain Observers in Nonlinear Feedback Control. Lecture # 2 Separation Principle
High-Gain Observers in Nonlinear Feedback Control Lecture # 2 Separation Principle High-Gain ObserversinNonlinear Feedback ControlLecture # 2Separation Principle p. 1/4 The Class of Systems ẋ = Ax + Bφ(x,
More informationSTABILITY. Phase portraits and local stability
MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),
More informationKolmogorov equations in Hilbert spaces IV
March 26, 2010 Other types of equations Let us consider the Burgers equation in = L 2 (0, 1) dx(t) = (AX(t) + b(x(t))dt + dw (t) X(0) = x, (19) where A = ξ 2, D(A) = 2 (0, 1) 0 1 (0, 1), b(x) = ξ 2 (x
More informationFinite-time Blowup of Semilinear PDEs via the Feynman-Kac Representation. CENTRO DE INVESTIGACIÓN EN MATEMÁTICAS GUANAJUATO, MEXICO
Finite-time Blowup of Semilinear PDEs via the Feynman-Kac Representation JOSÉ ALFREDO LÓPEZ-MIMBELA CENTRO DE INVESTIGACIÓN EN MATEMÁTICAS GUANAJUATO, MEXICO jalfredo@cimat.mx Introduction and backgrownd
More information1.Introduction: 2. The Model. Key words: Prey, Predator, Seasonality, Stability, Bifurcations, Chaos.
Dynamical behavior of a prey predator model with seasonally varying parameters Sunita Gakkhar, BrhamPal Singh, R K Naji Department of Mathematics I I T Roorkee,47667 INDIA Abstract : A dynamic model based
More informationScaling Limits of Waves in Convex Scalar Conservation Laws under Random Initial Perturbations
Scaling Limits of Waves in Convex Scalar Conservation Laws under Random Initial Perturbations Jan Wehr and Jack Xin Abstract We study waves in convex scalar conservation laws under noisy initial perturbations.
More informationStochastic relations of random variables and processes
Stochastic relations of random variables and processes Lasse Leskelä Helsinki University of Technology 7th World Congress in Probability and Statistics Singapore, 18 July 2008 Fundamental problem of applied
More informationXt i Xs i N(0, σ 2 (t s)) and they are independent. This implies that the density function of X t X s is a product of normal density functions:
174 BROWNIAN MOTION 8.4. Brownian motion in R d and the heat equation. The heat equation is a partial differential equation. We are going to convert it into a probabilistic equation by reversing time.
More informationMetastability for the Ginzburg Landau equation with space time white noise
Barbara Gentz gentz@math.uni-bielefeld.de http://www.math.uni-bielefeld.de/ gentz Metastability for the Ginzburg Landau equation with space time white noise Barbara Gentz University of Bielefeld, Germany
More informationHigh-Gain Observers in Nonlinear Feedback Control. Lecture # 3 Regulation
High-Gain Observers in Nonlinear Feedback Control Lecture # 3 Regulation High-Gain ObserversinNonlinear Feedback ControlLecture # 3Regulation p. 1/5 Internal Model Principle d r Servo- Stabilizing u y
More informationDerivation of Itô SDE and Relationship to ODE and CTMC Models
Derivation of Itô SDE and Relationship to ODE and CTMC Models Biomathematics II April 23, 2015 Linda J. S. Allen Texas Tech University TTU 1 Euler-Maruyama Method for Numerical Solution of an Itô SDE dx(t)
More informationON THE FIRST TIME THAT AN ITO PROCESS HITS A BARRIER
ON THE FIRST TIME THAT AN ITO PROCESS HITS A BARRIER GERARDO HERNANDEZ-DEL-VALLE arxiv:1209.2411v1 [math.pr] 10 Sep 2012 Abstract. This work deals with first hitting time densities of Ito processes whose
More informationContinuous Threshold Policy Harvesting in Predator-Prey Models
Continuous Threshold Policy Harvesting in Predator-Prey Models Jon Bohn and Kaitlin Speer Department of Mathematics, University of Wisconsin - Madison Department of Mathematics, Baylor University July
More informationDETERMINATION OF MODEL VALID PREDICTION PERIOD USING THE BACKWARD FOKKER-PLANCK EQUATION
.4 DETERMINATION OF MODEL VALID PREDICTION PERIOD USING THE BACKWARD FOKKER-PLANCK EQUATION Peter C. Chu, Leonid M. Ivanov, and C.W. Fan Department of Oceanography Naval Postgraduate School Monterey, California.
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 4. Bifurcations Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Local bifurcations for vector fields 1.1 The problem 1.2 The extended centre
More information