Lagrangian discretization of incompressibility, congestion and linear diffusion
|
|
- Leslie Wade
- 5 years ago
- Views:
Transcription
1 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 LJLL 19 Octobre 2018
2 Outline 1. General formalism 2. Incompressible Euler equations 3. Crowd motion 4. Linear diffusion
3 Lagrangian formulation of (some) PDEs We track the evolution of a population of particles (= probability distribution) ρ 0 P 2 (R d ), whose displacement field s : [0, T ] L 2 (ρ 0, R d ) satisfies ṡ E(s) (gradient flow) or s E(s) (hamiltonian flow) The distribution of particles at time t is then ρ(t) = s(t) # ρ 0. Main assumption: The energy/entropy E only depends on the distribution of particles, i.e. E(s) = E(s # ρ 0 ), with E : P 2 (R d ) R {+ }. { ρ log ρ if ρ P ac 2 (Ω) E ent (ρ) = + if not { 0 if ρ Leb Ω E cong (ρ) = + if not E inc (ρ) = E pot (ρ) = V d ρ { 0 if ρ = Leb Ω + if not NB: E is typically non-convex, with values in R {+ }.
4 Lagrangian-Eulerian dictionnary ṡ E(s) s E(s) E inc incompressible Euler equation E cong + E pot crowd motion pressureless Euler equation E ent + E pot linear Fokker-Planck equation isentropic Euler equation Wasserstein gradient-flows: [Otto 99], [Jordan-Kinderlehrer-Otto 98], [Ambrosio-Gigli-Savaré 08], etc. [Evans-Gangbo-Savin 08] incompressible Euler as geodesics: [Arnold 66], [Ebin Marsden 70], [Brenier], Numerical advantages of using a Lagrangian formulation: 1) handling many phases 2) tracking of individual particles / interfaces 3) allowing singular solutions (e.g. Dirac masses) Moreau-Yosida regularization: E ε (s) := inf s ṡ ε = E ε (s ε ) or s ε = E ε (s ε ) 1 2ε s s 2 + E(s )
5 2. Incompressible Euler equations Joint work with Thomas Gallouët
6 Geodesics on a submanifold of R d K = submanifold of R d, d K = distance to K, Π K = projection on K { { E(s) = 0 if s K Nor s K if s K E(s) = + if not if not E ε (s) = 1 2ε d2 K(m) E ε (s) = s Π K(s) ε a.e. Equation of geodesics: s E(s) ( s Nor s K = (T s K) and s K) approximated by s ε = E ε (s ε ) s ε + s ε Π K (s ε ) ε = 0 Nor K (Π K s ε ) Example: K = R {0} R 2, s(0) = 0, s ε (0) = (0, h ε ), ṡ(0) = ṡ ε (0) = (1, 0) h ε with s ε = (x ε, y ε ) we have { ẍ ε = 0 ÿ ε + 1 ε y ε = 0 i.e. { x ε (t) = t y ε (t) = h ε cos(t/ ε)
7 Incompressible fluids as geodesics Setting: Ω R d bounded connected, Ω = 1, ρ 0 = Leb Ω [Arnold 1966] Measure-preserving maps: K = {s L 2 (ρ 0, R d ) s # ρ 0 = ρ 0 }, E = i K s E(s) ( s Nor K (s) and s K) s 0 (x) s t (x) ( s = p s and s K) [ ] Lemma: a) Norid K = (T id K) { p p H 1 (Ω)} b) Nor s K = {w s w Nor id K} With u t := s t s 1 t (= velocity in Eulerian coordinates), one recovers Euler s eqs: t u + div(u u) = p in Ω divu = 0 in Ω u n = 0 on Ω Can this formulation be used for numerical computations (Brenier)? Minimizing geodesics (with JM Mirebeau, 15) / Cauchy pb (with T. Gallouet, 16).
8 Distance to measure-preserving maps I Measure-preserving maps: K = {s L 2 (ρ 0, R d ) s # ρ 0 = ρ 0 }, ρ 0 = Leb Ω Wasserstein distance: W 2 2(ρ, µ) = min{ x T (x) 2 d ρ(x) T # ρ = µ} Polar factorization (Brenier): (i) d 2 K(s) = W 2 2(ρ 0, s # ρ 0 ) (ii) If µ := s # ρ 0 is absolutely continuous,! optimal transport map S between µ and ρ 0 and then Π K (s) = {S s}. Moreover, S = φ where φ : Ω R is a convex function. [Brenier 92] With E = i K, E ε (s) = 1 2ε d2 K(s) E ε (s) = W 2 2(s # ρ 0, ρ 0 ) E ε (s) = s Π K(s) ε E ε (s) = id S ε s S = OT between s # ρ and ρ s ε = p ε s ε p ε = 1 ε ( 1 2 x 2 φ ε ) Formally, s ε = E ε (s ε ) [assuming µ ε Leb Ω ] µ ε = s ε# ρ 0 φ ε# µ ε = ρ 0, φ ε convex
9 Distance to measure-preserving maps II Definition: Given y 1,..., y N R d, let D(y 1,..., y N ) := N W 2 2(ρ 0, 1 N i δ y i ). Example: N = 900, ρ 0 = Leb [0,1] 2. y 1,..., y N 1 N D 0, 031 (y ( yi D) 1 i N i yi D) 1 i N
10 Distance to measure-preserving maps II Definition: Given y 1,..., y N R d, let D(y 1,..., y N ) := N W 2 2(ρ 0, 1 N i δ y i ). Example: N = 740, ρ 0 = Leb [0,1] 2. y 1,..., y N 1 N D 0, 14 ( y i D) 1 i N (y i yi D) 1 i N yi D(y 1,..., y N ) = B i (y 1,..., y N ) where a) T = is the optimal transport map from ρ 0 to 1 N b) V i := T 1 ({y i }) i δ y i c) B i (y 1,..., y N ) = N V i x d ρ 0 (x) = barycenter of V i
11 Space discretization of Euler s equations [Brenier, CMP 2000] [Gallouët M., 2016] Setting: let M N := {s = i y i1 ωi y 1,..., y N R d } L 2 (ρ 0, R d ) where ρ 0 (ω i ) = 1 N, ρ 0(ω i ω j ) = 0 and diam(ω i ) N 1 d. Given s = i y i1 ωi M N, (1) s + E ε MN (s) = 0 (2) ÿ i + 1 ε N yi D(y 1,..., y N ) = 0 (3) ÿ i + 1 ε N (y i B i (y 1,..., y N )) = 0 y i (t) is attached by a spring to the time-dependent barycenter B i (y 1 (t),..., y N (t)). Theorem: Let (u, p) be a regular (e.g. C 1,1 ) solution to Euler s equations and let s satisfy (1) with initial condition (s(0), ṡ(0)) = (Π MN (id), Π MN (u 0 )). Then, t [0, T ], ṡ t u t s t 2 L 2 (ρ 0,R d ) h2 N ε + ε + h N w. h N = N 1/d Proof: Gronwall on modulated energy E u (t) = 1 2 ṁ t u t m t ε 2 d 2 S(m t )
12 Numerical computation of D Computation of the Wasserstein distance: D(y 1,..., y N ) := N W 2 2(ρ 0, 1 N i δ y i ) We seek φ convex s.t. φ # ρ = 1 N i δ y i = φ = max i y i ψ i for some ψ R N. = φ # ρ = i ρ(v i(ψ))δ yi, where V i (ψ) = {x R d j, x y i ψ i x y j ψ j }. The partition (V i (ψ)) 1 i N can be computed in near-linear time in 2D/3D Corollary: Computing D(y 1,..., y N ) Finding φ convex such that φ # ρ = ν, Finding ψ R N s.t. i, ρ(v i (ψ)) = 1 N. Numerics: [Cullen Purser 84], [Oliker Prussner 88], [Aurenhammer Aronov Hoffman 98] [M. 11], [de Goes et al 12], [Lévy 15], [Kitagawa M. Thibert 16] [Lévy 15] + [KMT 16] handles points in R 3, within 20min on a laptop.
13 Numerical computation of D Computation of the Wasserstein distance: D(y 1,..., y N ) := N W 2 2(ρ 0, 1 N We seek φ convex s.t. φ # ρ = 1 N V i (ψ) = {x R d j, x y i ψ i x y j ψ j }. i δ y i = φ = max i y i ψ i for some ψ R N. = φ # ρ = i ρ(v i(ψ))δ yi, where The partition (V i (ψ)) 1 i N can be computed in near-linear time in 2D/3D i δ y i ) Corollary: Computing D(y 1,..., y N ) Finding φ convex such that φ # ρ = ν, Finding ψ R N s.t. i, ρ(v i (ψ)) = 1 N. For such a ψ R N, the Wasserstein distance and its gradient are given by D(y 1,..., y N ) = N i V i (ψ) x y i 2 d ρ 0 (x) 1 2 y i D(y 1,..., y N ) = y i bary(v i (ψ)) =: B i (y 1,..., y N )
14 Numerical result: Stationary flow on [0, 1] 2 Stationary flow on [0, 1] 2 : speed: u(x) = (cos(πx 1 ) sin(πx 2 ), sin(πx 1 ) cos(πx 2 )) pressure: p(x) = 1 4 (sin2 (πx 1 ) + sin 2 (πx 2 )) X = [0, 1] 2
15 Numerical result: Instabilities Objectives: Larger computations, with more complex behaviour. Preservation of the Hamiltonian by the discrete scheme. A. Discontinuous initial velocity u 0 = (.5, 0) u 0 = (1, 0) Hamiltonian preservation X = [0, 2] [.5,.5]/(x = 0 x = 2) 200k particles, 2000 timesteps, t max = 8
16 Numerical result: Instabilities Objectives: Larger computations, with more complex behaviour. Preservation of the Hamiltonian by the discrete scheme. B. Rayleigh-Taylor instability (Inhomogeneous fluid) ρ = 3 gravity velocity Verlet symplectic Euler ρ = 1 Hamiltonian preservation X = [ 1, 1] [ 3, 3] 50k particles, 2000 timesteps, t max = 2
17 3. Crowd motion Joint work with F. Santambrogio and F. Stra
18 Constrained gradient flows in R d K = compact convex subset of R d V C 1 (R d ), V 0 in all cases, y(0) = y 0 R d Constrained gradient flow: ẏ ( V (y) + i K (y)) ẏ = V (y) w y K w N y K ( ) ε-constrained gradient flow: ẏ ε = (V + 1 2ε d2 K)(y ε ) ẏ ε = V (y ε ) 1 ε (yε Π K (y ε )) Nor ΠK (y ε )K As ε 0, y ε converges locally uniformly, up to subsequence, to a weak solution of ( ). K Π K (y ε ) y ε
19 Crowd motion model [Maury, Roudneff-Chupin, Santambrogio 09] Ω R d compact, Ω > 1, V C 1 (R d ), V 0 K = {s L 2 (Ω, R d ) s # ρ 0 Leb Ω } where ρ 0 P ac (Ω) E pot (s) = V s d ρ 0 E cong (s) = i K (s) Lagrangian formulation: ṡ (E pot + E cong )(s) If s K and ρ = s # ρ 0 is the density of people, one expects T s K = {v s div(v) 0 on {ρ = 1} + no flux b.c.} i.e. E cong (s) = N s K = { p s p 0, p(1 ρ) = 0}. ṡ = V s p s ρ = s # ρ 0 Leb Ω p 0, p(1 ρ) = 0 Eulerian formulation: t ρ + div(ρv) = 0, v = V p p 0, p(1 ρ) = 0 ρ Leb Ω Gradient flow in (P(Ω), W 2 ) of the energy F(ρ) := V d ρ + i LebΩ (ρ)
20 Moreau-Yosida regularisation of E cong Admissible displacements: K = {s L 2 (ρ 0, R d ) s # ρ 0 Leb Ω }, E cong = i K E cong,ε (s) = 1 2ε d2 K(s), Formally, ṡ ε = ( E pot + E cong,ε )(s ε ) E cong,ε (s) = s Π K(s) ε ṡ ε = V (s ε ) p ε Π K (s ε ) σ ε = Π K (s ε ) # ρ 0 p ε 0, p ε (1 σ ε ) = 0 The computation of E cong,ε is justified using optimal transport techniques. Proposition: (i) d 2 K(s) = min σ P(Ω), σ LebΩ W 2 2(σ, s # ρ 0 ) (ii)!σ realizing the minimum (iii) If µ := s # ρ 0 is a.c. then Π K (s) = {S s} where S = OT between µ and σ Moreover, S 1 = id p with p 1-concave, p 0, p(1 σ) = 0. [M, RC, S 09] i.e. ε E cong,ε (s) = s Π K (s) = s S s = (S 1 id) S s = p Π K (s)
21 Discretization of the crowd motion model Eulerian schemes: [Maury, Roudneff-Chupin, Santambrogio 10] [Benamou, Carlier, Laborde 15] catching up ALG2-JKO Lagrangian setting: Let M N := {s = i y i1 ωi y 1,..., y N R d } L 2 (ρ 0, R d ) [with s = 1 N i δ y i ] We let µ N (t) = 1 N where ρ 0 (ω i ) = 1 N and ρ 0(ω i ω j ) = 0, { ṡ = ( E pot MN + E cong,εn MN )(s) s(0) = Π MN (id) { y i (t) = V (y i (t)) 1 ε N yi D (y 1 (t),..., y N (t)) y i (0) = N ω i x d ρ 0 where D (y 1,..., y N ) := N 2 inf σ 1 W 2 2(σ, 1 N i δ y i ) i δ y i (t) P(Ω) the corresponding distribution of particles.
22 Interpretation of the discretization Given y 1,..., y N R Nd : (example: y 1,..., y N iid in [0,.8] 2, Ω = [0, 2] 2 ) y 1,..., y N σ 1)!σ arg min ρ 1 W 2 2(ρ, 1 N i δ y i ) [Maury, Roudneff-Chupin, Santambrogio] 2)!T OT map from σ to 1 N [Brenier] i δ y i T 1 ({y i }) B i (y 1,..., y N ) 3) With V i := T 1 ({y i }) spt(σ), D (y 1,..., y N ) = N i V i x y i 2 d x 1 2 y i D (y 1,..., y N ) = N V i y i x d x =: y i B i (y 1,..., y N ) y i = V (y i ) 1 ε N (B i (y 1,..., y N ) y i )
23 Convergence theorem Let M N := {s = i y i1 ωi y 1,..., y N R d } L 2 (ρ 0, R d ), { ṡ N = ( E pot MN + E cong,εn MN )(s) s N (0) = Π MN (id) We denote µ N (t) = 1 N i δ y i (t) P(Ω) the corresponding distribution of particles. Theorem: If W 2 2(ρ 0, µ N (0)) ε N and ε N 0, then µ N ρ C 0 ([0, T ], P(Ω)), where ρ satisfies in the sense of distribution, for some p L 2 ([0, T ], H 1 (Ω)), t ρ + div(ρv) = 0, ρ K, v = V p, p 0, p(1 ρ) = 0, [M., Stra, Santambrogio 18]
24 Numerical examples 1) Square + gravity 1) Nonconvex domain + gravity 3) Gravity + aggregation 4) Two populations Fast solver for computing E cong,ε by Hugo Leclerc (LMO):.1s for N = 10 5!
25 4. Linear diffusion
26 Linear diffusion Setting: E ent (s) = E ent (s # ρ 0 ) with E ent (ρ) = { ρ log ρ if ρ P ac 2 (R d ) + if not Moreau-Yosida regularization in L 2 M-Y regularization in (P(R d ), W 2 ): i.e. E ent,ε (s) = E ent,ε (s # ρ 0 ) where E ent,ε (s) = min s E ent,ε (µ) = min σ P ac (R d ) 1 2ε s s 2 + E ent (s) 1 2ε W2 2(σ, µ) + E ent (σ) Letting D ent,ε (y 1,..., y N ) = E ent,ε ( 1 N i δ y i ), we can compute the ε-entropy of a point cloud in near-linear time in 2D/3D. Particle discretization of the heat equation: ẏ i = yi D ent,εn (y 1,..., y N ) yields a time dependent measure µ N = 1 N i δ y i. { t ρ + div(ρv) = 0 Under assumptions, µ N converges to a weak solution of v = log ρ
27 Summary Optimal transport can be used to construct particle discretization of (some) evolution PDEs with congestion/incompressibility constraints or diffusion terms Some perspectives: more general models: non-linear diffusion, reaction terms (e.g. tumor growth), multiple population, surface tension, interaction forces (e.g. Keller-Segel)? convergence towards Lagrangian solutions (à la [Evans-Gangbo-Savin 05])? getting rid of ε for 2nd order models (cf [de Goes et al, SIGGRAPH 2015])? Thank you for your attention!
Convex discretization of functionals involving the Monge-Ampère operator
Convex discretization of functionals involving the Monge-Ampère operator Quentin Mérigot CNRS / Université Paris-Dauphine Joint work with J.D. Benamou, G. Carlier and É. Oudet GeMeCod Conference October
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 informationConvex discretization of functionals involving the Monge-Ampère operator
Convex discretization of functionals involving the Monge-Ampère operator Quentin Mérigot CNRS / Université Paris-Dauphine Joint work with J.D. Benamou, G. Carlier and É. Oudet Workshop on Optimal Transport
More informationFrom Euler to Monge and vice versa. Y.B. Lecture 2: From Euler to Vlasov through Monge-Ampère and Kantorovich.
From Euler to Monge and vice versa. Y.B. Lecture 2: From Euler to Vlasov through Monge-Ampère and Kantorovich. Yann Brenier, Mikaela Iacobelli, Filippo Santambrogio, Paris, Durham, Lyon. MFO SEMINAR 1842,
More informationarxiv: v1 [math.na] 13 May 2015
Minimal geodesics along volume preserving maps, through semi-discrete optimal transport Quentin Mérigot Jean-Marie Mirebeau arxiv:1505.03306v1 [math.na] 13 May 2015 August 18, 2018 Abstract We introduce
More informationThe semi-geostrophic equations - a model for large-scale atmospheric flows
The semi-geostrophic equations - a model for large-scale atmospheric flows Beatrice Pelloni, University of Reading with M. Cullen (Met Office), D. Gilbert, T. Kuna INI - MFE Dec 2013 Introduction - Motivation
More informationBranched transport limit of the Ginzburg-Landau functional
Branched transport limit of the Ginzburg-Landau functional Michael Goldman CNRS, LJLL, Paris 7 Joint work with S. Conti, F. Otto and S. Serfaty Introduction Superconductivity was first observed by Onnes
More informationDiscrete Ricci curvature via convexity of the entropy
Discrete Ricci curvature via convexity of the entropy Jan Maas University of Bonn Joint work with Matthias Erbar Simons Institute for the Theory of Computing UC Berkeley 2 October 2013 Starting point McCann
More informationMinimal time mean field games
based on joint works with Samer Dweik and Filippo Santambrogio PGMO Days 2017 Session Mean Field Games and applications EDF Lab Paris-Saclay November 14th, 2017 LMO, Université Paris-Sud Université Paris-Saclay
More informationNumerical Optimal Transport and Applications
Numerical Optimal Transport and Applications Gabriel Peyré Joint works with: Jean-David Benamou, Guillaume Carlier, Marco Cuturi, Luca Nenna, Justin Solomon www.numerical-tours.com Histograms in Imaging
More informationCONGESTED AGGREGATION VIA NEWTONIAN INTERACTION
CONGESTED AGGREGATION VIA NEWTONIAN INTERACTION KATY CRAIG, INWON KIM, AND YAO YAO Abstract. We consider a congested aggregation model that describes the evolution of a density through the competing effects
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 informationA Modest Proposal for MFG with Density Constraints
A Modest Proposal for MFG with Density Constraints Mean Field Games and related topics, Rome, May 011 Filippo Santambrogio October 31, 011 Abstract We consider a typical problem in Mean Field Games: the
More informationA DIFFUSIVE MODEL FOR MACROSCOPIC CROWD MOTION WITH DENSITY CONSTRAINTS
A DIFFUSIVE MODEL FOR MACROSCOPIC CROWD MOTION WITH DENSITY CONSTRAINTS Abstract. In the spirit of the macroscopic crowd motion models with hard congestion (i.e. a strong density constraint ρ ) introduced
More informationUNIQUENESS ISSUES FOR EVOLUTIVE EQUATIONS WITH DENSITY CONSTRAINTS
UNIQUENESS ISSUES FOR EVOLUTIVE EQUATIONS WITH DENSITY CONSTRAINTS SIMONE DI MARINO AND ALPÁR RICHÁRD MÉSZÁROS Abstract. In this paper we present some basic uniqueness results for evolutive equations under
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 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 informationMean field games via probability manifold II
Mean field games via probability manifold II Wuchen Li IPAM summer school, 2018. Introduction Consider a nonlinear Schrödinger equation hi t Ψ = h2 1 2 Ψ + ΨV (x) + Ψ R d W (x, y) Ψ(y) 2 dy. The unknown
More informationMean Field Limits for Ginzburg-Landau Vortices
Mean Field Limits for Ginzburg-Landau Vortices Sylvia Serfaty Courant Institute, NYU & LJLL, Université Pierre et Marie Curie Conference in honor of Yann Brenier, January 13, 2017 Fields Institute, 2003
More informationFree Energy, Fokker-Planck Equations, and Random walks on a Graph with Finite Vertices
Free Energy, Fokker-Planck Equations, and Random walks on a Graph with Finite Vertices Haomin Zhou Georgia Institute of Technology Jointly with S.-N. Chow (Georgia Tech) Wen Huang (USTC) Yao Li (NYU) Research
More informationHamiltonian aspects of fluid dynamics
Hamiltonian aspects of fluid dynamics CDS 140b Joris Vankerschaver jv@caltech.edu CDS 01/29/08, 01/31/08 Joris Vankerschaver (CDS) Hamiltonian aspects of fluid dynamics 01/29/08, 01/31/08 1 / 34 Outline
More informationOn Lagrangian solutions for the semi-geostrophic system with singular initial data
On Lagrangian solutions for the semi-geostrophic system with singular initial data Mikhail Feldman Department of Mathematics University of Wisconsin-Madison Madison, WI 5376, USA feldman@math.wisc.edu
More informationLEAST ACTION PRINCIPLES FOR INCOMPRESSIBLE FLOWS AND OPTIMAL TRANSPORT BETWEEN SHAPES JIAN-GUO LIU
LEAST ACTION PRINCIPLES FOR INCOMPRESSIBLE FLOWS AND OPTIMAL TRANSPORT BETWEEN SHAPES JIAN-GUO LIU Department of Physics and Department of Mathematics Duke University, Durham, NC 27708, USA ROBERT L. PEGO
More informationMeasure-valued - strong uniqueness for hyperbolic systems
Measure-valued - strong uniqueness for hyperbolic systems joint work with Tomasz Debiec, Eduard Feireisl, Ondřej Kreml, Agnieszka Świerczewska-Gwiazda and Emil Wiedemann Institute of Mathematics Polish
More informationA quantum heat equation 5th Spring School on Evolution Equations, TU Berlin
A quantum heat equation 5th Spring School on Evolution Equations, TU Berlin Mario Bukal A. Jüngel and D. Matthes ACROSS - Centre for Advanced Cooperative Systems Faculty of Electrical Engineering and Computing
More informationarxiv: v1 [math.ap] 10 Apr 2013
QUASI-STATIC EVOLUTION AND CONGESTED CROWD TRANSPORT DAMON ALEXANDER, INWON KIM, AND YAO YAO arxiv:1304.3072v1 [math.ap] 10 Apr 2013 Abstract. We consider the relationship between Hele-Shaw evolution with
More informationSpaces with Ricci curvature bounded from below
Spaces with Ricci curvature bounded from below Nicola Gigli February 23, 2015 Topics 1) On the definition of spaces with Ricci curvature bounded from below 2) Analytic properties of RCD(K, N) spaces 3)
More informationSobolev regularity for the Monge-Ampère equation, with application to the semigeostrophic equations
Sobolev regularity for the Monge-Ampère equation, with application to the semigeostrophic equations Alessio Figalli Abstract In this note we review some recent results on the Sobolev regularity of solutions
More informationAcoustic Energy Estimates in Inhomogeneous Moving Media
Subject: FORUM ACUSTICUM 1999 Abstract Acoustic Energy Estimates in Inhomogeneous Moving Media, NASA Langley Research Center, Hampton, Virginia Mark Farris, Midwestern State University, Wichita Falls,
More informationOn Lagrangian solutions for the semi-geostrophic system with singular initial data
On Lagrangian solutions for the semi-geostrophic system with singular initial data Mikhail Feldman Department of Mathematics University of Wisconsin-Madison Madison, WI 5376, USA feldman@math.wisc.edu
More informationMean Field Limits for Ginzburg-Landau Vortices
Mean Field Limits for Ginzburg-Landau Vortices Sylvia Serfaty Université P. et M. Curie Paris 6, Laboratoire Jacques-Louis Lions & Institut Universitaire de France Jean-Michel Coron s 60th birthday, June
More informationLogarithmic Sobolev Inequalities
Logarithmic Sobolev Inequalities M. Ledoux Institut de Mathématiques de Toulouse, France logarithmic Sobolev inequalities what they are, some history analytic, geometric, optimal transportation proofs
More information{Euclidean, metric, and Wasserstein} gradient flows: an overview
Bull. Math. Sci. (2017) 7:87 154 DOI 10.1007/s13373-017-0101-1 {Euclidean, metric, and Wasserstein} gradient flows: an overview Filippo Santambrogio 1 Received: 11 September 2016 / Revised: 9 February
More informationStationary states for the aggregation equation with power law attractive-repulsive potentials
Stationary states for the aggregation equation with power law attractive-repulsive potentials Daniel Balagué Joint work with J.A. Carrillo, T. Laurent and G. Raoul Universitat Autònoma de Barcelona BIRS
More informationStochastic Homogenization for Reaction-Diffusion Equations
Stochastic Homogenization for Reaction-Diffusion Equations Jessica Lin McGill University Joint Work with Andrej Zlatoš June 18, 2018 Motivation: Forest Fires ç ç ç ç ç ç ç ç ç ç Motivation: Forest Fires
More informationA shor t stor y on optimal transpor t and its many applications
Snapshots of modern mathematics from Oberwolfach 13/2018 A shor t stor y on optimal transpor t and its many applications Filippo Santambrogio We present some examples of optimal transport problems and
More informationarxiv: v2 [math.ap] 28 Nov 2016
ONE-DIMENSIONAL SAIONARY MEAN-FIELD GAMES WIH LOCAL COUPLING DIOGO A. GOMES, LEVON NURBEKYAN, AND MARIANA PRAZERES arxiv:1611.8161v [math.ap] 8 Nov 16 Abstract. A standard assumption in mean-field game
More informationLarge-scale atmospheric circulation, semi-geostrophic motion and Lagrangian particle methods
Large-scale atmospheric circulation, semi-geostrophic motion and Lagrangian particle methods Colin Cotter (Imperial College London) & Sebastian Reich (Universität Potsdam) Outline 1. Hydrostatic and semi-geostrophic
More 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 informationUniqueness of the solution to the Vlasov-Poisson system with bounded density
Uniqueness of the solution to the Vlasov-Poisson system with bounded density Grégoire Loeper December 16, 2005 Abstract In this note, we show uniqueness of weak solutions to the Vlasov- Poisson system
More informationÉquation de Burgers avec particule ponctuelle
Équation de Burgers avec particule ponctuelle Nicolas Seguin Laboratoire J.-L. Lions, UPMC Paris 6, France 7 juin 2010 En collaboration avec B. Andreianov, F. Lagoutière et T. Takahashi Nicolas Seguin
More informationKramers formula for chemical reactions in the context of Wasserstein gradient flows. Michael Herrmann. Mathematical Institute, University of Oxford
eport no. OxPDE-/8 Kramers formula for chemical reactions in the context of Wasserstein gradient flows by Michael Herrmann Mathematical Institute, University of Oxford & Barbara Niethammer Mathematical
More informationSEPARABILITY AND COMPLETENESS FOR THE WASSERSTEIN DISTANCE
SEPARABILITY AND COMPLETENESS FOR THE WASSERSTEIN DISTANCE FRANÇOIS BOLLEY Abstract. In this note we prove in an elementary way that the Wasserstein distances, which play a basic role in optimal transportation
More informationOptimal Transportation. Nonlinear Partial Differential Equations
Optimal Transportation and Nonlinear Partial Differential Equations Neil S. Trudinger Centre of Mathematics and its Applications Australian National University 26th Brazilian Mathematical Colloquium 2007
More informationFrom slow diffusion to a hard height constraint: characterizing congested aggregation
3 2 1 Time Position -0.5 0.0 0.5 0 From slow iffusion to a har height constraint: characterizing congeste aggregation Katy Craig University of California, Santa Barbara BIRS April 12, 2017 2 collective
More informationJean-David Benamou 1, Guillaume Carlier 2 and Maxime Laborde 3
ESAIM: PROCEEDINGS AND SURVEYS, June 6, Vol. 54, p. -7 B. Düring, C.-B. Schönlieb and M.-T. Wolfram, Editors AN AUGMENTED LAGRANGIAN APPROACH TO WASSERSTEIN GRADIENT FLOWS AND APPLICATIONS Jean-David Benamou,
More informationDistance-Divergence Inequalities
Distance-Divergence Inequalities Katalin Marton Alfréd Rényi Institute of Mathematics of the Hungarian Academy of Sciences Motivation To find a simple proof of the Blowing-up Lemma, proved by Ahlswede,
More informationLocal semiconvexity of Kantorovich potentials on non-compact manifolds
Local semiconvexity of Kantorovich potentials on non-compact manifolds Alessio Figalli, Nicola Gigli Abstract We prove that any Kantorovich potential for the cost function c = d / on a Riemannian manifold
More informationWeak and Measure-Valued Solutions of the Incompressible Euler Equations
Weak and Measure-Valued Solutions for Euler 1 / 12 Weak and Measure-Valued Solutions of the Incompressible Euler Equations (joint work with László Székelyhidi Jr.) October 14, 2011 at Carnegie Mellon University
More informationIntroduction to optimal transport
Introduction to optimal transport Nicola Gigli May 20, 2011 Content Formulation of the transport problem The notions of c-convexity and c-cyclical monotonicity The dual problem Optimal maps: Brenier s
More informationNumerical Approximation of L 1 Monge-Kantorovich Problems
Rome March 30, 2007 p. 1 Numerical Approximation of L 1 Monge-Kantorovich Problems Leonid Prigozhin Ben-Gurion University, Israel joint work with J.W. Barrett Imperial College, UK Rome March 30, 2007 p.
More informationOptimal Transport: A Crash Course
Optimal Transport: A Crash Course Soheil Kolouri and Gustavo K. Rohde HRL Laboratories, University of Virginia Introduction What is Optimal Transport? The optimal transport problem seeks the most efficient
More informationOptimal transportation and action-minimizing measures
Scuola Normale Superiore of Pisa and École Normale Supérieure of Lyon Phd thesis 24th October 27 Optimal transportation and action-minimizing measures Alessio Figalli a.figalli@sns.it Advisor Prof. Luigi
More informationGradient Flow. Chang Liu. April 24, Tsinghua University. Chang Liu (THU) Gradient Flow April 24, / 91
Gradient Flow Chang Liu Tsinghua University April 24, 2017 Chang Liu (THU) Gradient Flow April 24, 2017 1 / 91 Contents 1 Introduction 2 Gradient flow in the Euclidean space Variants of Gradient Flow in
More informationMean Field Limits for Ginzburg-Landau Vortices
Mean Field Limits for Ginzburg-Landau Vortices Sylvia Serfaty Courant Institute, NYU Rivière-Fabes symposium, April 28-29, 2017 The Ginzburg-Landau equations u : Ω R 2 C u = u ε 2 (1 u 2 ) Ginzburg-Landau
More informationEULERIAN DERIVATIONS OF NON-INERTIAL NAVIER-STOKES EQUATIONS
EULERIAN DERIVATIONS OF NON-INERTIAL NAVIER-STOKES EQUATIONS ML Combrinck, LN Dala Flamengro, a div of Armscor SOC Ltd & University of Pretoria, Council of Scientific and Industrial Research & University
More informationG: Uniform Convergence of Fourier Series
G: Uniform Convergence of Fourier Series From previous work on the prototypical problem (and other problems) u t = Du xx 0 < x < l, t > 0 u(0, t) = 0 = u(l, t) t > 0 u(x, 0) = f(x) 0 < x < l () we developed
More information0.2. CONSERVATION LAW FOR FLUID 9
0.2. CONSERVATION LAW FOR FLUID 9 Consider x-component of Eq. (26), we have D(ρu) + ρu( v) dv t = ρg x dv t S pi ds, (27) where ρg x is the x-component of the bodily force, and the surface integral is
More informationGeodesic Density Tracking with Applications to Data Driven Modeling
Geodesic Density Tracking with Applications to Data Driven Modeling Abhishek Halder and Raktim Bhattacharya Department of Aerospace Engineering Texas A&M University College Station, TX 77843-3141 American
More informationVarious lecture notes for
Various lecture notes for 18311. R. R. Rosales (MIT, Math. Dept., 2-337) April 12, 2013 Abstract Notes, both complete and/or incomplete, for MIT s 18.311 (Principles of Applied Mathematics). These notes
More informationA variational proof of partial regularity for optimal transportation maps between sets
A variational proof of partial regularity for optimal transportation maps between sets M. Goldman F. Otto October 24, 27 Abstract We provide a new proof of the known partial regularity result for the optimal
More informationarxiv: v2 [math.na] 2 Aug 2013
A VISCOSITY SOLUTION APPROACH TO THE MONGE-AMPÈRE FORMULATION OF THE OPTIMAL TRANSPORTATION PROBLEM arxiv:1208.4873v2 [math.na] 2 Aug 2013 JEAN-DAVID BENAMOU, BRITTANY D. FROESE, AND ADAM M. OBERMAN Abstract.
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 informationarxiv: v1 [math.ap] 10 Oct 2013
The Exponential Formula for the Wasserstein Metric Katy Craig 0/0/3 arxiv:30.292v math.ap] 0 Oct 203 Abstract We adapt Crandall and Liggett s method from the Banach space case to give a new proof of the
More informationEULER SPRAYS AND WASSERSTEIN GEOMETRY OF THE SPACE OF SHAPES JIAN-GUO LIU
EULER SPRAYS AND WASSERSTEIN GEOMETRY OF THE SPACE OF SHAPES JIAN-GUO LIU Department of Physics and Department of Mathematics Duke University, Durham, NC 27708, USA ROBERT L. PEGO AND DEJAN SLEPČEV Department
More informationWeak Solutions to a Fractional Fokker-Planck Equation via Splitting and Wasserstein Gradient Flow
Weak Solutions to a Fractional Fokker-Planck Equation via Splitting and Wasserstein Gradient Flow by Malcolm Bowles B.Sc., University of Victoria, 212 A Thesis Submitted in Partial Fulfillment of the Requirements
More informationEntropy and Relative Entropy
Entropy and Relative Entropy Joshua Ballew University of Maryland October 24, 2012 Outline Hyperbolic PDEs Entropy/Entropy Flux Pairs Relative Entropy Weak-Strong Uniqueness Weak-Strong Uniqueness for
More informationOn the velocities of flows consisting of cyclically monotone maps
On the velocities of flows consisting of cyclically monotone maps A. Tudorascu Department of Mathematics West Virginia University Morgantown, WV 2656, USA adriant@math.wvu.edu August 25, 21 Abstract Motivated
More informationNumerical Schemes for Calculating the Discrete Wasserstein Distance
Numerical Schemes for Calculating the Discrete Wasserstein Distance Benjamin Cabrera Born 9th April 99 in Groß-Gerau, Germany February 6, 6 Master s Thesis Mathematics Advisor: Prof. Dr. Martin Rumpf Mathematical
More informationOn dissipation distances for reaction-diffusion equations the Hellinger-Kantorovich distance
Weierstrass Institute for! Applied Analysis and Stochastics On dissipation distances for reaction-diffusion equations the Matthias Liero, joint work with Alexander Mielke, Giuseppe Savaré Mohrenstr. 39
More informationAbsolutely continuous curves in Wasserstein spaces with applications to continuity equation and to nonlinear diffusion equations
Università degli Studi di Pavia Dipartimento di Matematica F. Casorati Dottorato di ricerca in Matematica e Statistica Absolutely continuous curves in Wasserstein spaces with applications to continuity
More informationThe quasi-periodic Frenkel-Kontorova model
The quasi-periodic Frenkel-Kontorova model Philippe Thieullen (Bordeaux) joint work with E. Garibaldi (Campinas) et S. Petite (Amiens) Korean-French Conference in Mathematics Pohang, August 24-28,2015
More informationA regularised particle method for linear and nonlinear diffusion
1/25 A regularised particle method for linear and nonlinear diffusion Francesco Patacchini Department of Mathematical Sciences, Carnegie Mellon University Joint work with J. A. Carrillo (Imperial College
More informationOptimal Transport for Applied Mathematicians
Optimal Transport for Applied Mathematicians Calculus of Variations, PDEs and Modelling Filippo Santambrogio 1 1 Laboratoire de Mathématiques d Orsay, Université Paris Sud, 91405 Orsay cedex, France filippo.santambrogio@math.u-psud.fr,
More informationDynamic and Stochastic Brenier Transport via Hopf-Lax formulae on Was
Dynamic and Stochastic Brenier Transport via Hopf-Lax formulae on Wasserstein Space With many discussions with Yann Brenier and Wilfrid Gangbo Brenierfest, IHP, January 9-13, 2017 ain points of the
More informationPhysics 451/551 Theoretical Mechanics. G. A. Krafft Old Dominion University Jefferson Lab Lecture 18
Physics 451/551 Theoretical Mechanics G. A. Krafft Old Dominion University Jefferson Lab Lecture 18 Theoretical Mechanics Fall 018 Properties of Sound Sound Waves Requires medium for propagation Mainly
More informationParallel-in-time integrators for Hamiltonian systems
Parallel-in-time integrators for Hamiltonian systems Claude Le Bris ENPC and INRIA Visiting Professor, The University of Chicago joint work with X. Dai (Paris 6 and Chinese Academy of Sciences), F. Legoll
More informationOptimal Transport for Data Analysis
Optimal Transport for Data Analysis Bernhard Schmitzer 2017-05-16 1 Introduction 1.1 Reminders on Measure Theory Reference: Ambrosio, Fusco, Pallara: Functions of Bounded Variation and Free Discontinuity
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 informationA Very Brief Introduction to Conservation Laws
A Very Brief Introduction to Wen Shen Department of Mathematics, Penn State University Summer REU Tutorial, May 2013 Summer REU Tutorial, May 2013 1 / The derivation of conservation laws A conservation
More informationScalar conservation laws with moving density constraints arising in traffic flow modeling
Scalar conservation laws with moving density constraints arising in traffic flow modeling Maria Laura Delle Monache Email: maria-laura.delle monache@inria.fr. Joint work with Paola Goatin 14th International
More informationSeparation for the stationary Prandtl equation
Separation for the stationary Prandtl equation Anne-Laure Dalibard (UPMC) with Nader Masmoudi (Courant Institute, NYU) February 13th-17th, 217 Dynamics of Small Scales in Fluids ICERM, Brown University
More informationExistence and uniqueness results for the MFG system
Existence and uniqueness results for the MFG system P. Cardaliaguet Paris-Dauphine Based on Lasry-Lions (2006) and Lions lectures at Collège de France IMA Special Short Course Mean Field Games", November
More informationLocal Density Approximation for the Almost-bosonic Anyon Gas. Michele Correggi
Local Density Approximation for the Almost-bosonic Anyon Gas Michele Correggi Università degli Studi Roma Tre www.cond-math.it QMATH13 Many-body Systems and Statistical Mechanics joint work with D. Lundholm
More informationGeometric Inference for Probability distributions
Geometric Inference for Probability distributions F. Chazal 1 D. Cohen-Steiner 2 Q. Mérigot 2 1 Geometrica, INRIA Saclay, 2 Geometrica, INRIA Sophia-Antipolis 2009 July 8 F. Chazal, D. Cohen-Steiner, Q.
More informationWeek 2 Notes, Math 865, Tanveer
Week 2 Notes, Math 865, Tanveer 1. Incompressible constant density equations in different forms Recall we derived the Navier-Stokes equation for incompressible constant density, i.e. homogeneous flows:
More informationEntropic curvature-dimension condition and Bochner s inequality
Entropic curvature-dimension condition and Bochner s inequality Kazumasa Kuwada (Ochanomizu University) joint work with M. Erbar and K.-Th. Sturm (Univ. Bonn) German-Japanese conference on stochastic analysis
More informationSUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS
SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS A. RÖSCH AND R. SIMON Abstract. An optimal control problem for an elliptic equation
More informationLagrangian Dynamics & Mixing
V Lagrangian Dynamics & Mixing See T&L, Ch. 7 The study of the Lagrangian dynamics of turbulence is, at once, very old and very new. Some of the earliest work on fluid turbulence in the 1920 s, 30 s and
More informationThe incompressible Navier-Stokes equations in vacuum
The incompressible, Université Paris-Est Créteil Piotr Bogus law Mucha, Warsaw University Journées Jeunes EDPistes 218, Institut Elie Cartan, Université de Lorraine March 23th, 218 Incompressible Navier-Stokes
More informationOn a Class of Multidimensional Optimal Transportation Problems
Journal of Convex Analysis Volume 10 (2003), No. 2, 517 529 On a Class of Multidimensional Optimal Transportation Problems G. Carlier Université Bordeaux 1, MAB, UMR CNRS 5466, France and Université Bordeaux
More informationA STRATEGY FOR NON-STRICTLY CONVEX TRANSPORT COSTS AND THE EXAMPLE OF x y p IN R 2
A STRATEGY FOR NON-STRICTLY CONVEX TRANSPORT COSTS AND THE EXAMPLE OF x y p IN R 2 GUILLAUME CARLIER, LUIGI DE PASCALE, AND FILIPPO SANTAMBROGIO Abstract. This paper deals with the existence of optimal
More informationAPPLICATIONS DU TRANSPORT OPTIMAL DANS LES JEUX À POTENTIEL
APPLICATIONS DU TRANSPORT OPTIMAL DANS LES JEUX À POTENTIEL Adrien Blanchet IAST/TSE/Université Toulouse 1 Capitole Journées SMAI-MODE 216 En collaboration avec G. Carlier & P. Mossay & F. Santambrogio
More informationThe 2-d isentropic compressible Euler equations may have infinitely many solutions which conserve energy
The -d isentropic compressible Euler equations may have infinitely many solutions which conserve energy Simon Markfelder Christian Klingenberg September 15, 017 Dept. of Mathematics, Würzburg University,
More informationarxiv: v1 [math.ap] 29 May 2018
Non-uniqueness of admissible weak solution to the Riemann problem for the full Euler system in D arxiv:805.354v [math.ap] 9 May 08 Hind Al Baba Christian Klingenberg Ondřej Kreml Václav Mácha Simon Markfelder
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 informationVARIATIONAL PARTICLE SCHEMES FOR THE POROUS MEDIUM EQUATION AND FOR THE SYSTEM OF ISENTROPIC EULER EQUATIONS
VARIATIONAL PARTICLE SCHEMES FOR THE POROUS MEDIUM EQUATION AND FOR THE SYSTEM OF ISENTROPIC EULER EQUATIONS MICHAEL WESTDICKENBERG AND JON WILKENING Abstract. Both the porous medium equation and the system
More informationBlowup dynamics of an unstable thin-film equation
Blowup dynamics of an unstable thin-film equation IMA Summer Program Geometrical Singularities and Singular Geometries July 15, 2008. (IMA, July 15, 2008.) 1 / 12 Thin-film equations (TFE) u t = (u n u
More informationNumerische Mathematik
Numer. Math. (2) 84: 375 393 Digital Object Identifier (DOI) 1.17/s21199117 Numerische Mathematik c Springer-Verlag 2 A computational fluid mechanics solution to the Monge-Kantorovich mass transfer problem
More informationFREE BOUNDARY PROBLEMS IN FLUID MECHANICS
FREE BOUNDARY PROBLEMS IN FLUID MECHANICS ANA MARIA SOANE AND ROUBEN ROSTAMIAN We consider a class of free boundary problems governed by the incompressible Navier-Stokes equations. Our objective is to
More information