A Semi-Lagrangian discretization of non linear fokker Planck equations
|
|
- Rodger Nicholson
- 5 years ago
- Views:
Transcription
1 A Semi-Lagrangian discretization of non linear fokker Planck equations E. Carlini Università Sapienza di Roma joint works with F.J. Silva RICAM, Linz November 2-25, 206 / 3
2 Outline A Semi-Lagrangian scheme for a nonlinear Fokker-Planck equation 2 Convergence Analysis 3 Numerical simulation 2 / 3
3 A nonlinear Fokker-Planck equation The nonlinear FP equation { t m 2 i,j x i x j (a ij [m](x, t)m) + div(b[m](x, t) m) = 0 in R d R + m(, 0) = m 0 ( ) in R d, () describes the evolution of the law of the diffusion process X(t) where dx(s ) = b[m](x(s ), s ))ds + σ[m](x(s ), s )dw (s ), X(s) = x, b[m] : R d R + R R d, given a vector field, depending non locally on m; σ[m] : R d R + R R d r, (A[m](x, t)) i,j = a ij [m](x, t) := (σ[m](x, t)σ[m](x, t) ) ij, is the diffusion matrix (possible degenerate),depending non locally on m; the density of the initial law is given by m 0 ; m 0 0 and R d m 0 (x)dx =. (2) 3 / 3
4 Some applications Non local interactions due to collective phenomena (biophysics, social behavior) Mean Filed Games: b[m](x, t) = DH(Dv[m](x, t)) where v[m] is the solution of a backward HJB { t v σ2 2 v + H(Dv) = f(x, m(t)) v(x, T ) = g(x, m(t ) Since b[m](x, t) depends on m(s) in all time 0 s T we gets an Implicit scheme: Fixed Point Iterations needed to solve it Hughes model b[m](x, t) = f 2 (m(x, t))dv[m](x, t) where v[m] is the solution of a stationary HJB Dv = f(m(x, t)) Since b[m](x, t) depends on m(s) only in time past 0 s t we gets an Explicit scheme 4 / 3
5 Weak solution For many problems of interest, () has only a formal meaning, henceforth we mean m is a weak solution of () if for all φ Cc (R d ) we have that d φ(x)dm(t)(x) dt R d = R d 2 a ij [m](x, t) xi x j φ(x) + b i [m](x, t)φ xi (x) dm(t)(x) i,j i [ = 2 tr(a[m](x, t)d2 φ(x)) + b[m](x, t) Dφ(x)] dm(t)(x) R d (3) 5 / 3
6 Numerical approximation of FP Some references Linear Chang and Cooper (970): finite difference scheme s.t. preserves positivity, equilibrium states and mass of the distribution (explicit version is stable under a parabolic type CFL condition) Kushner (976): finite difference via probabilistic method Naess and Johnsen (993) : Path Integration Method (time integration of the evolution probability density function, where the transition probability density is approximated by a Gaussian) Jordan, Kinderlehrer,Otto (998): variational scheme for FokkerPlanck equations for which the drift term is given by the gradient of a potential (JKO) (It preserves positivity, and mass of the distribution, costly in high dimension). Achdou,Camilli,Capuzzo Dolcetta (202): implicit finite difference scheme s.t. preserves positivity, and mass of the distribution. Non Linear Drozdov, Morillo (995): finite difference scheme s.t. preserves equilibrium states and mass of the distribution, (high order). Benamou, Carlier, Laborde (205): semi implicit variant of JKO 6 / 3
7 Numerical approximation based on SL for FP We propose a Semi-Lagrangian (SL) scheme for non linear FP equation s.t. it is first order accurate (numerically) it is explicit and it allows for large time steps it preserves the positivity of the density and conserves its integral equals to Ref. Silva, Camilli (203), Carlini and Silva (204,205) The scheme has been applied to numerically compute the solution of Mean Field Games model (with F. Silva) a regularized Hughes model for pedestrian flow (with A.Festa, F.Silva, M.T. Wolfram) Hughes model for pedestrian flow with different congestion penalty function (with A.Festa, F.Silva) 7 / 3
8 A Semi-discrete in time SL for a nonlinear d Fokker-Planck equation Given t, we define t k = k t, k = 0,..., N. We multiply the first equation in () by a smooth test function φ with compact support and integrate by parts to get: φ(x)dm(t k+ ) = φ(x)dm(t k ) (4) R R tk+ + [b[m](x, t)dφ(x) + 2 σ2 [m](x, t)d 2 φ(x)]dm(t)dt. t k R 8 / 3
9 Semi-discrete in time We first approximate (4) as φ(x)dm(t k+ ) = R [φ(x) + tb[m k ](x, t k )Dφ(x) + t 2 σ2 [m k ](x, t k )D 2 φ(x)]dm(t k+ ). R Note that the right hand side corresponds to a Taylor expansion. Hence we write φ(x)dm(t k+ ) = R [φ(x + tb[m k ](x, t k ) + tσ 2 [m k ](x, t k )dm(t k )+ 2 R [φ(x + tb[m k ](x, t k ) tσ 2 [m k ](x, t k )]dm(t k ). 2 R 9 / 3
10 Fully-discrete in space Given x > 0 consider a space grid, for i =,..., M E i = [x i 2 x, x i + 2 x], m i,k := E i dm(t k ). By the standard rectangular quadrature formula, we get φ(x j )m j,k+ = j Z 2 φ(φ + j [m k)])m j,k + φ(φ j 2 [m k)])m j,k, j Z j Z (5) (6) where, for µ R M, j Z, k = 0,..., N, we have defined Φ ± j [µ] := x j + t b[µ](x j, t k ) ± tσ[µ](x j, t k ). (7) 0 / 3
11 Interpretation of the scheme by means of characteristics Note that Φ defined in (7) (with µ j = m(x j, t k )) can be interpreted as a single Euler step approximation of dx(s) = b[m] (X(s), s)) ds + σ[m] (X(s), s)) dw (s), s (t k, t k+ ), X(t k ) = x j, (8) with a random walk discretization of the Brownian motion W ( ). Indeed, considering a random value Z in R such that f P(Z = ) = P(Z = ) = 2 (9) the function Φ ± j [m k] corresponds to one realization of x j + tb[m k ](x j, t k ) ± 2 tσ[m k ](x j, t k )Z. / 3
12 2 -,5 - -0,5 0 0,5,5 2 Fully-discrete in space Figure: P-basis function, β i Given i Z and setting in (6) φ = β i we have for all i Z m i,k+ = G(m k, i, k) = ( ) 2 j Z β i (Φ + j [m k]) + β i (Φ j [m k]) m i,0 = E i m 0 (x)dx m j,k (0) Non-negative : m i,k 0 for k = 0,..., N, i Z Mass conservative : x i m i,k = for k = 0,..., N Generalizable to any dimension Generalizable to handle Dirichlet and Neumann Boundary conditions Generalizable to handle degeneracy of the diffusion matrix 2 / 3
13 Boundary Condition Let us consider () in Ω R 2, we denote by ˆn(x) the outward normal to Ω at x Ω. Homogeneous Neumann condition on Γ N Ω J[m](x, t) ˆn = 0 on Γ N, () where J i [m](x, t) = b i [m](x, t)m k x k (a i,k [m](x, t)m) Dirichlet boundary condition on Γ D := Ω \ Γ N m = 0 on Γ D. (2) The idea is to reflect the discrete characteristics when they intersect Γ N and to cut them at Γ D. 3 / 3
14 Boundary Condition: Dirichlet, d = 2 Let us define, for l =, 2 ĥ l i,± = inf{γ > 0 ; x + γb[m k ](x i, t k ) ± 2dγσ l [m k ](x i, t k ) Γ D } h, and redefine Φ l i,± as Φ l i,± := x i + ĥl,± i,k b[m k](x i, t k ) + Ref. Falcone, Ferretti(203) 2dĥl i,± σ l[m k ](x i, t k ). 4 / 3
15 Boundary Condition: Neumann, d = 2 In order to consider the reflection when exiting at Γ N, we use a symmetrized Euler scheme: we project the discrete characteristics on Ω using the operator P : R d Ω, defined as z, if z Ω, P (z) := 2w z, if z / Ω, (3) where w := argmin z w. w Ω Ref. M. Bossy and E. Gobet and D.Talay, (2004) 5 / 3
16 Unstructured grids Let T be a triangulation of Ω, with x the maximum diameter of the triangles. We call {β i ; i =,..., M} the set of affine shape functions related to the triangular mesh, such that β i (x j ) = δ i,j and i β i(x) = for each x Ω. Median dual control volume E i := E i,t, T T :x i T E i,t := {x T : β j (x) β i (x) j i}. We approximate the solution m by a sequence {m k } k, such that m i,k dm(t k ) E i 6 / 3
17 Consistency Let us denote m x, t a proper extension in R d [0, T ] of the discrete approximation m i,k and P a propability space, metrized by the Kantorowich-Rubinstein distance: { } d (µ, µ 2 ) := sup f(x)d(µ (x) µ 2 (x)) ; Df. R d Proposition (Weak consistency) Assuming that b and σ are continuous and Lipschitz w.r. to x, for every φ C0 (R) if ( x n, t n ) is s.t. as n ( x n, t n ) 0 and xn t n 0, and m xn, t n m in C(0, T ; P (R d )), then for k n such that t k n t lim φ(x) [dm xn, t n n (t kn +) dm xn, t n (0)] R d t [ = 2 tr(a[m](x, s)d2 φ) + b[m](x, s) Dφ ] dm(s). R 0 7 / 3
18 Stability Proposition Suppose that ( x) 2 = O( t), b and σ continuous functions s.t. there exists C > 0 b(x, t, µ) C x σ(x, t, µ) C x, for any x, t, µ. Then, there exists a constant C > 0 (independent of ( x, t)) such that: Proposition R d x 2 dm x, t (t) C t [0, T ]. (4) Suppose that ( x) 2 = O( t). Then, there exists a constant C > 0 (independent of ( x, t)) such that: d (m x, t (t ), m x, t (t 2 )) C t t 2 2 t, t 2 [0, T ]. (5) 8 / 3
19 Convergence Theorem Suppose that ( x) 2 = O( t), b and σcontinuous functions Lipschtiz w.r. to x and there exists C > 0 such that for any x, t, µ b(x, t, µ) C x σ(x, t, µ) C x, then m x, t m in C([0, T ], P ), where m is solution of () and m x, t is a sequence of solution of (0) (extended to a general dimension d). 9 / 3
20 Interacting species Let m and m 2 be the densities of two interacting species t m = div(m ( E (m ) + U (m, m 2 ))) t m 2 = div(m 2 ( E 2 (m 2 ) + U 2 (m, m 2 ))) b.c. internal energy E (m) = E 2 (m) = 3 2 m2 K ε (repulsive self-interactions) cross interaction potential U (m, m 2 ) = W m + W 2 m 2 U (m, m 2 ) = W 2 m + W 22 m 2 Ref. J.D.Benamou, G. Carlier, M.Laborde 20 / 3
21 Interacting species Let us set W = 3 x 2 2,W 2 = 3 x 2 2 W 2 = x 2 2,W 22 = x 2 2 we solve the system on Ω = [.5,.5] [.5,.5] and [0, T ] = [0, 3], we set h = 0.03 and t = h/2 Homogenous Neumann Condition 2 / 3
22 Interacting species t = 0 t = 0.6 t = t =.8 t = 2.4 t = 3 22 / 3
23 Interacting species without self-interactions (E = E 2 = 0 Figure: t = 3 23 / 3
24 Interacting species without selfinteractions (E = E 2 = 0) t = 0 t = 0.6 t = t =.8 t = 2.4 t = 3 24 / 3
25 Interacting species without selfinteractions (E = E 2 = 0 Figure: t = 3 25 / 3
26 The Hughes model In 2002 R. Hughes proposed a macroscopic model for pedestrian dynamics given by t m(x, t) div(m(x, t) f 2 (m(x, t))du(x, t)) = 0, Du(x, t) = f(m(x, t)), (6) where x Ω, t (0, T ], T R +, m corresponds to the density of a crowd trying to exit a domain as fast as possible u the weighted shortest distance to a target, for example an exit. Typical choice for f: f(m) = m where corresponds to the maximum scaled pedestrian density. Few analytic results are available, all of them restricted to spatial dimension one. In this case the dependence of the drift on m is local!, and the convergence analysis does not apply. 26 / 3
27 A modified version of the Hughes model (Joint work with A.Festa,F.J Silva, M.T. Wolframm) A modified version corresponds to t m(x, t) ε m(x, t) div(m(x, t) f 2 (m(x, t))du(x, t)) = 0, ε u(x, t) + 2 Du(x, t) 2 = 2f 2 (m(x, t)) + δ, (7) in Ω (0, T ). δ > 0 prevents the blow-up of the cost when approaching the maximum density one. ε > 0 represents a small diffusion. Ref: Di Francesco M, Markowich P A, Pietschmann F P, Wolfram M T (20) 27 / 3
28 Exit scenario with two doors Figure: Density contour lines at time t = 0.3, ε = The exits are marked with the letter E. The position of the two exits induces an asymmetric splitting of the crowd 28 / 3
29 Exit scenario with two doors Figure: Density contour lines at time t = 0.3 (left) and t =.2 (right). M 0 = 0.7 and ε = Initially, a large part of the crowd chooses to move to the right exit. After a while this exit gets congested, inducing a part of the population to change objective. 29 / 3
30 References Conclusions: we have presented a new SL scheme for a non linear FP main advantage: it is explicit, it allows large time steps, it preserves the positivity and it conserves the mass On going work convergence analysis SL scheme for non linear FP convergence analysis SL scheme for HJB on bounded domain with Neumann end Dirichlet condition Open questions high order SL scheme for FP 30 / 3
31 References E.Carlini, F.J. Silva, A Semi-Lagrangian scheme for a degenerate second order Mean Field Game system, DCDS-A, ( 205) E.Carlini, F.J. Silva, A Semi-Lagrangian scheme for the Fokker-Planck equation, 2nd IFAC Workshop on Control of Systems Governed E.Carlini, A. Festa, F.J. Silva, M.T. Wolfram, A Semi-Lagrangian scheme for a modified version of the Hughes model for pedestrian flow, Dynamic Games and Applications, Springer, (206) 3 / 3
32 References E.Carlini, F.J. Silva, A Semi-Lagrangian scheme for a degenerate second order Mean Field Game system, DCDS-A, ( 205) E.Carlini, F.J. Silva, A Semi-Lagrangian scheme for the Fokker-Planck equation, 2nd IFAC Workshop on Control of Systems Governed E.Carlini, A. Festa, F.J. Silva, M.T. Wolfram, A Semi-Lagrangian scheme for a modified version of the Hughes model for pedestrian flow, Dynamic Games and Applications, Springer, (206) THANK YOU 3 / 3
A Semi-Lagrangian scheme for a regularized version of the Hughes model for pedestrian flow
A Semi-Lagrangian scheme for a regularized version of the Hughes model for pedestrian flow Adriano FESTA Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy of Sciences
More informationMean field games and related models
Mean field games and related models Fabio Camilli SBAI-Dipartimento di Scienze di Base e Applicate per l Ingegneria Facoltà di Ingegneria Civile ed Industriale Email: Camilli@dmmm.uniroma1.it Web page:
More informationarxiv: v2 [math.na] 2 Jan 2018
ON THE DISCRETIZATION OF SOME NONLINEAR FOKKER-PLANCK-KOLMOGOROV EQUATIONS AND APPLICATIONS ELISABETTA CARLINI AND FRANCISCO J. SILVA arxiv:1708.004v [math.na] Jan 018 Abstract. In this work, we consider
More informationWeak solutions to Fokker-Planck equations and Mean Field Games
Weak solutions to Fokker-Planck equations and Mean Field Games Alessio Porretta Universita di Roma Tor Vergata LJLL, Paris VI March 6, 2015 Outlines of the talk Brief description of the Mean Field Games
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
More informationGradient Flows: Qualitative Properties & Numerical Schemes
Gradient Flows: Qualitative Properties & Numerical Schemes J. A. Carrillo Imperial College London RICAM, December 2014 Outline 1 Gradient Flows Models Gradient flows Evolving diffeomorphisms 2 Numerical
More informationNumerical methods for Mean Field Games: additional material
Numerical methods for Mean Field Games: additional material Y. Achdou (LJLL, Université Paris-Diderot) July, 2017 Luminy Convergence results Outline 1 Convergence results 2 Variational MFGs 3 A numerical
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 informationFiltered scheme and error estimate for first order Hamilton-Jacobi equations
and error estimate for first order Hamilton-Jacobi equations Olivier Bokanowski 1 Maurizio Falcone 2 2 1 Laboratoire Jacques-Louis Lions, Université Paris-Diderot (Paris 7) 2 SAPIENZA - Università di Roma
More informationMean Field Games on networks
Mean Field Games on networks Claudio Marchi Università di Padova joint works with: S. Cacace (Rome) and F. Camilli (Rome) C. Marchi (Univ. of Padova) Mean Field Games on networks Roma, June 14 th, 2017
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationTopology-preserving diffusion equations for divergence-free vector fields
Topology-preserving diffusion equations for divergence-free vector fields Yann BRENIER CNRS, CMLS-Ecole Polytechnique, Palaiseau, France Variational models and methods for evolution, Levico Terme 2012
More informationAn adaptive RBF-based Semi-Lagrangian scheme for HJB equations
An adaptive RBF-based Semi-Lagrangian scheme for HJB equations Roberto Ferretti Università degli Studi Roma Tre Dipartimento di Matematica e Fisica Numerical methods for Hamilton Jacobi equations in optimal
More informationarxiv: v2 [math.na] 2 May 2018
A FULLY-DISCRETE SCHEME FOR SYSTEMS OF NONLINEAR FOKKER-PLANCK-KOLMOGOROV EQUATIONS ELISABETTA CARLINI AND FRANCISCO J. SILVA arxiv:7.6v [math.na May 8 Abstract. We consider a system of Fokker-Planck-Kolmogorov
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 informationStationary mean-field games Diogo A. Gomes
Stationary mean-field games Diogo A. Gomes We consider is the periodic stationary MFG, { ɛ u + Du 2 2 + V (x) = g(m) + H ɛ m div(mdu) = 0, (1) where the unknowns are u : T d R, m : T d R, with m 0 and
More informationVariational approach to mean field games with density constraints
1 / 18 Variational approach to mean field games with density constraints Alpár Richárd Mészáros LMO, Université Paris-Sud (based on ongoing joint works with F. Santambrogio, P. Cardaliaguet and F. J. Silva)
More informationON MEAN FIELD GAMES. Pierre-Louis LIONS. Collège de France, Paris. (joint project with Jean-Michel LASRY)
Collège de France, Paris (joint project with Jean-Michel LASRY) I INTRODUCTION II A REALLY SIMPLE EXAMPLE III GENERAL STRUCTURE IV TWO PARTICULAR CASES V OVERVIEW AND PERSPECTIVES I. INTRODUCTION New class
More informationControlled Diffusions and Hamilton-Jacobi Bellman Equations
Controlled Diffusions and Hamilton-Jacobi Bellman Equations Emo Todorov Applied Mathematics and Computer Science & Engineering University of Washington Winter 2014 Emo Todorov (UW) AMATH/CSE 579, Winter
More information2.20 Fall 2018 Math Review
2.20 Fall 2018 Math Review September 10, 2018 These notes are to help you through the math used in this class. This is just a refresher, so if you never learned one of these topics you should look more
More informationMean-Field Games with non-convex Hamiltonian
Mean-Field Games with non-convex Hamiltonian Martino Bardi Dipartimento di Matematica "Tullio Levi-Civita" Università di Padova Workshop "Optimal Control and MFGs" Pavia, September 19 21, 2018 Martino
More informationStochastic Differential Equations.
Chapter 3 Stochastic Differential Equations. 3.1 Existence and Uniqueness. One of the ways of constructing a Diffusion process is to solve the stochastic differential equation dx(t) = σ(t, x(t)) dβ(t)
More informationThe HJB-POD approach for infinite dimensional control problems
The HJB-POD approach for infinite dimensional control problems M. Falcone works in collaboration with A. Alla, D. Kalise and S. Volkwein Università di Roma La Sapienza OCERTO Workshop Cortona, June 22,
More informationIntroduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes Heat equation The simple parabolic PDE with the initial values u t = K 2 u 2 x u(0, x) = u 0 (x) and some boundary conditions
More informationp 1 ( Y p dp) 1/p ( X p dp) 1 1 p
Doob s inequality Let X(t) be a right continuous submartingale with respect to F(t), t 1 P(sup s t X(s) λ) 1 λ {sup s t X(s) λ} X + (t)dp 2 For 1 < p
More informationFast Marching Methods for Front Propagation Lecture 1/3
Outline for Front Propagation Lecture 1/3 M. Falcone Dipartimento di Matematica SAPIENZA, Università di Roma Autumn School Introduction to Numerical Methods for Moving Boundaries ENSTA, Paris, November,
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
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 informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More informationPH.D. PRELIMINARY EXAMINATION MATHEMATICS
UNIVERSITY OF CALIFORNIA, BERKELEY Dept. of Civil and Environmental Engineering FALL SEMESTER 2014 Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem
More informationSébastien Chaumont a a Institut Élie Cartan, Université Henri Poincaré Nancy I, B. P. 239, Vandoeuvre-lès-Nancy Cedex, France. 1.
A strong comparison result for viscosity solutions to Hamilton-Jacobi-Bellman equations with Dirichlet condition on a non-smooth boundary and application to parabolic problems Sébastien Chaumont a a Institut
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationGaussian processes for inference in stochastic differential equations
Gaussian processes for inference in stochastic differential equations Manfred Opper, AI group, TU Berlin November 6, 2017 Manfred Opper, AI group, TU Berlin (TU Berlin) inference in SDE November 6, 2017
More informationOn nonlocal Hamilton-Jacobi equations related to jump processes, some recent results
On nonlocal Hamilton-Jacobi equations related to jump processes, some recent results Daria Ghilli Institute of mathematics and scientic computing, Graz, Austria 25/11/2016 Workshop "Numerical methods for
More informationFinite Difference and Finite Element Methods
Finite Difference and Finite Element Methods Georgy Gimel farb COMPSCI 369 Computational Science 1 / 39 1 Finite Differences Difference Equations 3 Finite Difference Methods: Euler FDMs 4 Finite Element
More informationDiscretization of SDEs: Euler Methods and Beyond
Discretization of SDEs: Euler Methods and Beyond 09-26-2006 / PRisMa 2006 Workshop Outline Introduction 1 Introduction Motivation Stochastic Differential Equations 2 The Time Discretization of SDEs Monte-Carlo
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationApproximation of BSDEs using least-squares regression and Malliavin weights
Approximation of BSDEs using least-squares regression and Malliavin weights Plamen Turkedjiev (turkedji@math.hu-berlin.de) 3rd July, 2012 Joint work with Prof. Emmanuel Gobet (E cole Polytechnique) Plamen
More informationA Semi-Lagrangian Scheme for the Open Table Problem in Granular Matter Theory
A Semi-Lagrangian Scheme for the Open Table Problem in Granular Matter Theory M. Falcone and S. Finzi Vita Abstract We introduce and analyze a new scheme for the approximation of the two-layer model proposed
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationHomogeniza*ons in Perforated Domain. Ki Ahm Lee Seoul Na*onal University
Homogeniza*ons in Perforated Domain Ki Ahm Lee Seoul Na*onal University Outline 1. Perforated Domain 2. Neumann Problems (joint work with Minha Yoo; interes*ng discussion with Li Ming Yeh) 3. Dirichlet
More informationApproximating diffusions by piecewise constant parameters
Approximating diffusions by piecewise constant parameters Lothar Breuer Institute of Mathematics Statistics, University of Kent, Canterbury CT2 7NF, UK Abstract We approximate the resolvent of a one-dimensional
More informationCOURSE Iterative methods for solving linear systems
COURSE 0 4.3. Iterative methods for solving linear systems Because of round-off errors, direct methods become less efficient than iterative methods for large systems (>00 000 variables). An iterative scheme
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 informationDensities for the Navier Stokes equations with noise
Densities for the Navier Stokes equations with noise Marco Romito Università di Pisa Universitat de Barcelona March 25, 2015 Summary 1 Introduction & motivations 2 Malliavin calculus 3 Besov bounds 4 Other
More informationSome Aspects of Universal Portfolio
1 Some Aspects of Universal Portfolio Tomoyuki Ichiba (UC Santa Barbara) joint work with Marcel Brod (ETH Zurich) Conference on Stochastic Asymptotics & Applications Sixth Western Conference on Mathematical
More informationMATH 819 FALL We considered solutions of this equation on the domain Ū, where
MATH 89 FALL. The D linear wave equation weak solutions We have considered the initial value problem for the wave equation in one space dimension: (a) (b) (c) u tt u xx = f(x, t) u(x, ) = g(x), u t (x,
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 informationNumerical schemes of resolution of stochastic optimal control HJB equation
Numerical schemes of resolution of stochastic optimal control HJB equation Elisabeth Ottenwaelter Journée des doctorants 7 mars 2007 Équipe COMMANDS Elisabeth Ottenwaelter (INRIA et CMAP) Numerical schemes
More informationIntroduction to Random Diffusions
Introduction to Random Diffusions The main reason to study random diffusions is that this class of processes combines two key features of modern probability theory. On the one hand they are semi-martingales
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
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 informationA Narrow-Stencil Finite Difference Method for Hamilton-Jacobi-Bellman Equations
A Narrow-Stencil Finite Difference Method for Hamilton-Jacobi-Bellman Equations Xiaobing Feng Department of Mathematics The University of Tennessee, Knoxville, U.S.A. Linz, November 23, 2016 Collaborators
More informationPart I: Mean field game models in pedestrian dynamics
Part I: Mean field game models in pedestrian dynamics M.T. Wolfram University of Warwick, TU Munich and RICAM Graduate Summer School on Mean field games and applications Mean field games and applications
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 informationWeek 6 Notes, Math 865, Tanveer
Week 6 Notes, Math 865, Tanveer. Energy Methods for Euler and Navier-Stokes Equation We will consider this week basic energy estimates. These are estimates on the L 2 spatial norms of the solution u(x,
More informationThe method of lines (MOL) for the diffusion equation
Chapter 1 The method of lines (MOL) for the diffusion equation The method of lines refers to an approximation of one or more partial differential equations with ordinary differential equations in just
More informationStochastic optimal control with rough paths
Stochastic optimal control with rough paths Paul Gassiat TU Berlin Stochastic processes and their statistics in Finance, Okinawa, October 28, 2013 Joint work with Joscha Diehl and Peter Friz Introduction
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 informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationLECTURE 5: THE METHOD OF STATIONARY PHASE
LECTURE 5: THE METHOD OF STATIONARY PHASE Some notions.. A crash course on Fourier transform For j =,, n, j = x j. D j = i j. For any multi-index α = (α,, α n ) N n. α = α + + α n. α! = α! α n!. x α =
More informationSuboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids
Suboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids Jochen Garcke joint work with Axel Kröner, INRIA Saclay and CMAP, Ecole Polytechnique Ilja Kalmykov, Universität
More informationA Hybrid Method for the Wave Equation. beilina
A Hybrid Method for the Wave Equation http://www.math.unibas.ch/ beilina 1 The mathematical model The model problem is the wave equation 2 u t 2 = (a 2 u) + f, x Ω R 3, t > 0, (1) u(x, 0) = 0, x Ω, (2)
More informationTD 1: Hilbert Spaces and Applications
Université Paris-Dauphine Functional Analysis and PDEs Master MMD-MA 2017/2018 Generalities TD 1: Hilbert Spaces and Applications Exercise 1 (Generalized Parallelogram law). Let (H,, ) be a Hilbert space.
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 informationThe first order quasi-linear PDEs
Chapter 2 The first order quasi-linear PDEs The first order quasi-linear PDEs have the following general form: F (x, u, Du) = 0, (2.1) where x = (x 1, x 2,, x 3 ) R n, u = u(x), Du is the gradient of u.
More informationDiscretization of Stochastic Differential Systems With Singular Coefficients Part II
Discretization of Stochastic Differential Systems With Singular Coefficients Part II Denis Talay, INRIA Sophia Antipolis joint works with Mireille Bossy, Nicolas Champagnat, Sylvain Maire, Miguel Martinez,
More informationConservation Laws & Applications
Rocky Mountain Mathematics Consortium Summer School Conservation Laws & Applications Lecture V: Discontinuous Galerkin Methods James A. Rossmanith Department of Mathematics University of Wisconsin Madison
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationA positivity-preserving high order discontinuous Galerkin scheme for convection-diffusion equations
A positivity-preserving high order discontinuous Galerkin scheme for convection-diffusion equations Sashank Srinivasan a, Jonathan Poggie a, Xiangxiong Zhang b, a School of Aeronautics and Astronautics,
More information(Received 10 December 2011, accepted 15 February 2012) x x 2 B(x) u, (1.2) A(x) +
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol13(212) No4,pp387-395 Numerical Solution of Fokker-Planck Equation Using the Flatlet Oblique Multiwavelets Mir Vahid
More informationGradient Schemes for an Obstacle Problem
Gradient Schemes for an Obstacle Problem Y. Alnashri and J. Droniou Abstract The aim of this work is to adapt the gradient schemes, discretisations of weak variational formulations using independent approximations
More informationA path integral approach to the Langevin equation
A path integral approach to the Langevin equation - Ashok Das Reference: A path integral approach to the Langevin equation, A. Das, S. Panda and J. R. L. Santos, arxiv:1411.0256 (to be published in Int.
More 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 informationFIXED POINT ITERATIONS
FIXED POINT ITERATIONS MARKUS GRASMAIR 1. Fixed Point Iteration for Non-linear Equations Our goal is the solution of an equation (1) F (x) = 0, where F : R n R n is a continuous vector valued mapping in
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationPDEs in Image Processing, Tutorials
PDEs in Image Processing, Tutorials Markus Grasmair Vienna, Winter Term 2010 2011 Direct Methods Let X be a topological space and R: X R {+ } some functional. following definitions: The mapping R is lower
More informationCourse Summary Math 211
Course Summary Math 211 table of contents I. Functions of several variables. II. R n. III. Derivatives. IV. Taylor s Theorem. V. Differential Geometry. VI. Applications. 1. Best affine approximations.
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationWeak convergence and large deviation theory
First Prev Next Go To Go Back Full Screen Close Quit 1 Weak convergence and large deviation theory Large deviation principle Convergence in distribution The Bryc-Varadhan theorem Tightness and Prohorov
More informationLagrangian discretization of incompressibility, congestion and linear diffusion
Lagrangian discretization of incompressibility, congestion and linear diffusion Quentin Mérigot Joint works with (1) Thomas Gallouët (2) Filippo Santambrogio Federico Stra (3) Hugo Leclerc Seminaire du
More informationSINC PACK, and Separation of Variables
SINC PACK, and Separation of Variables Frank Stenger Abstract This talk consists of a proof of part of Stenger s SINC-PACK computer package (an approx. 400-page tutorial + about 250 Matlab programs) that
More informationApproximation of fluid-structure interaction problems with Lagrange multiplier
Approximation of fluid-structure interaction problems with Lagrange multiplier Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 30, 2016 Outline
More informationIntroduction to the Numerical Solution of IVP for ODE
Introduction to the Numerical Solution of IVP for ODE 45 Introduction to the Numerical Solution of IVP for ODE Consider the IVP: DE x = f(t, x), IC x(a) = x a. For simplicity, we will assume here that
More informationPart I. Discrete Models. Part I: Discrete Models. Scientific Computing I. Motivation: Heat Transfer. A Wiremesh Model (2) A Wiremesh Model
Part I: iscrete Models Scientific Computing I Module 5: Heat Transfer iscrete and Continuous Models Tobias Neckel Winter 04/05 Motivation: Heat Transfer Wiremesh Model A Finite Volume Model Time ependent
More informationOptimal stopping time formulation of adaptive image filtering
Optimal stopping time formulation of adaptive image filtering I. Capuzzo Dolcetta, R. Ferretti 19.04.2000 Abstract This paper presents an approach to image filtering based on an optimal stopping time problem
More informationPoisson Jumps in Credit Risk Modeling: a Partial Integro-differential Equation Formulation
Poisson Jumps in Credit Risk Modeling: a Partial Integro-differential Equation Formulation Jingyi Zhu Department of Mathematics University of Utah zhu@math.utah.edu Collaborator: Marco Avellaneda (Courant
More informationSingular Perturbations of Stochastic Control Problems with Unbounded Fast Variables
Singular Perturbations of Stochastic Control Problems with Unbounded Fast Variables Joao Meireles joint work with Martino Bardi and Guy Barles University of Padua, Italy Workshop "New Perspectives in Optimal
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 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 informationOptimization. Escuela de Ingeniería Informática de Oviedo. (Dpto. de Matemáticas-UniOvi) Numerical Computation Optimization 1 / 30
Optimization Escuela de Ingeniería Informática de Oviedo (Dpto. de Matemáticas-UniOvi) Numerical Computation Optimization 1 / 30 Unconstrained optimization Outline 1 Unconstrained optimization 2 Constrained
More informationδ Substituting into the differential equation gives: x i+1 x i δ f(t i,x i ) (3)
Solving Differential Equations Numerically Differential equations are ubiquitous in Physics since the laws of nature often take on a simple form when expressed in terms of infinitesimal changes of the
More informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
More informationExit times of diffusions with incompressible drift
Exit times of diffusions with incompressible drift Gautam Iyer, Carnegie Mellon University gautam@math.cmu.edu Collaborators: Alexei Novikov, Penn. State Lenya Ryzhik, Stanford University Andrej Zlatoš,
More informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
More informationarxiv: v1 [math.ap] 15 Jul 2016
Mean Field Games models of segregation Yves Achdou, Martino Bardi and Marco Cirant September 4, 8 Abstract arxiv:7.445v [math.ap] 5 Jul This paper introduces and analyses some models in the framework of
More informationNumerical explorations of a forward-backward diffusion equation
Numerical explorations of a forward-backward diffusion equation Newton Institute KIT Programme Pauline Lafitte 1 C. Mascia 2 1 SIMPAF - INRIA & U. Lille 1, France 2 Univ. La Sapienza, Roma, Italy September
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
More informationNotes on Mean Field Games
Notes on Mean Field Games (from P.-L. Lions lectures at Collège de France) Pierre Cardaliaguet September 27, 203 Contents Introduction 2 2 Nash equilibria in games with a large number of players 4 2. Symmetric
More informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
More information