Lecture 1: Main Models & Basics of Wasserstein Distance
|
|
- Cassandra Holland
- 5 years ago
- Views:
Transcription
1 Lecture 1: Main Models & Basics of Wasserstein Distance J. A. Carrillo ICREA - Universitat Autònoma de Barcelona Methods and Models of Kinetic Theory
2 Outline 1 Presentation of models Nonlinear diffusions Kinetic models for granular gases 2 Wasserstein Distance: Basics Definition Properties 3 Contractivity in 1D Diffusive Models Dissipative Models Towards any dimension
3 Leading examples Nonlinear Diffusions.- u t = um, (x R d, t > 0) Inelastic Dissipative Models: Nonlinear friction equations.- 1 : f t + v f x =»Z (v w) v w γ( v w )f (x, w, t) dw f (x, v, t) v R Inelastic Dissipative Models: Nonlinear Boltzmann-type kinetic equations.- f + (v x)f =Qe(f, f ) t Conservation of mass and center of mass/mean velocity. Spreading (diffusions) versus Concentration (dissipative models). 1 Benedetto, D., Caglioti, E., Pulvirenti, M., Math. Mod. and Num. An. (1997); Benedetto, D., Caglioti, E., Carrillo, J. A., Pulvirenti, M., J. Stat. Phys. (1998).
4 Leading examples Nonlinear Diffusions.- u t = um, (x R d, t > 0) Inelastic Dissipative Models: Nonlinear friction equations.- 1 : f t + v f x =»Z (v w) v w γ( v w )f (x, w, t) dw f (x, v, t) v R Inelastic Dissipative Models: Nonlinear Boltzmann-type kinetic equations.- f + (v x)f =Qe(f, f ) t Conservation of mass and center of mass/mean velocity. Spreading (diffusions) versus Concentration (dissipative models). 1 Benedetto, D., Caglioti, E., Pulvirenti, M., Math. Mod. and Num. An. (1997); Benedetto, D., Caglioti, E., Carrillo, J. A., Pulvirenti, M., J. Stat. Phys. (1998).
5 Leading examples Nonlinear Diffusions.- u t = um, (x R d, t > 0) Inelastic Dissipative Models: Nonlinear friction equations.- 1 : f t + v f x =»Z (v w) v w γ( v w )f (x, w, t) dw f (x, v, t) v R Inelastic Dissipative Models: Nonlinear Boltzmann-type kinetic equations.- f + (v x)f =Qe(f, f ) t Conservation of mass and center of mass/mean velocity. Spreading (diffusions) versus Concentration (dissipative models). 1 Benedetto, D., Caglioti, E., Pulvirenti, M., Math. Mod. and Num. An. (1997); Benedetto, D., Caglioti, E., Carrillo, J. A., Pulvirenti, M., J. Stat. Phys. (1998).
6 Leading examples Nonlinear Diffusions.- u t = um, (x R d, t > 0) Inelastic Dissipative Models: Nonlinear friction equations.- 1 : f t + v f x =»Z (v w) v w γ( v w )f (x, w, t) dw f (x, v, t) v R Inelastic Dissipative Models: Nonlinear Boltzmann-type kinetic equations.- f + (v x)f =Qe(f, f ) t Conservation of mass and center of mass/mean velocity. Spreading (diffusions) versus Concentration (dissipative models). 1 Benedetto, D., Caglioti, E., Pulvirenti, M., Math. Mod. and Num. An. (1997); Benedetto, D., Caglioti, E., Carrillo, J. A., Pulvirenti, M., J. Stat. Phys. (1998).
7 Leading examples Nonlinear Diffusions.- u t = um, (x R d, t > 0) Inelastic Dissipative Models: Nonlinear friction equations.- 1 : f t + v f x =»Z (v w) v w γ( v w )f (x, w, t) dw f (x, v, t) v R Inelastic Dissipative Models: Nonlinear Boltzmann-type kinetic equations.- f + (v x)f =Qe(f, f ) t Conservation of mass and center of mass/mean velocity. Spreading (diffusions) versus Concentration (dissipative models). 1 Benedetto, D., Caglioti, E., Pulvirenti, M., Math. Mod. and Num. An. (1997); Benedetto, D., Caglioti, E., Carrillo, J. A., Pulvirenti, M., J. Stat. Phys. (1998).
8 Nonlinear diffusions Outline 1 Presentation of models Nonlinear diffusions Kinetic models for granular gases 2 Wasserstein Distance: Basics Definition Properties 3 Contractivity in 1D Diffusive Models Dissipative Models Towards any dimension
9 Nonlinear diffusions Nonlinear diffusions Self-similar solution: Barenblatt profile 2.- An explicit self-similar solution that is integrable for m > (d 2)/d: B m(x, t) = t d/λ C 2m x 2 1 m t 2/λ «1/(1 m) for m 1 where λ = d(m 1) + 2 and C > 0 is determined to have unit mass. It verifies that B m(x, t) converges weakly-* as measures towards δ 0 as t Zeldovich, Ya. B., Barenblatt, G. I. Doklady, USSR Academy of Sciences (1958).
10 Nonlinear diffusions Nonlinear diffusions Asymptotic behaviour: Comparison methods.- Friedman & Kamin (1980) (completed by Vázquez (1997)). Given an initial data in the class then, for any (d 2)/d < m X 0 = {u 0 L 1 (R d ) : u 0 0}, lim t u(, t) Bm(, t) L 1 = 0, and lim t td/λ u(, t) B m(, t) L = 0.
11 Nonlinear diffusions Nonlinear diffusions ρ t = div(xρ + ρm ), (x R d, t > 0) ρ(x, t = 0) = ρ 0(x) 0, (x R d ) ρ(x, t) = e dt u(e t x, 1 λ (eλt 1)) λ = d(m 1) + 2 u t = um, (x R d, t > 0) u(x, t = 0) = u 0(x) 0, (x R d ) 3 Exponential decay to equilibria translates into algebraic decay to Barenblatt profiles. 3 J. A. Carrillo, G. Toscani, Indiana Math. Univ. J. (2000); F. Otto, Comm. PDE (2001); J. Dolbeault, M. del Pino, J. Math. Pures Appl. (2002); J.L. Vázquez, J. Evol. Eq
12 Nonlinear diffusions Nonlinear diffusions ρ t = div(xρ + ρm ), (x R d, t > 0) ρ(x, t = 0) = ρ 0(x) 0, (x R d ) ρ(x, t) = e dt u(e t x, 1 λ (eλt 1)) λ = d(m 1) + 2 u t = um, (x R d, t > 0) u(x, t = 0) = u 0(x) 0, (x R d ) 3 Exponential decay to equilibria translates into algebraic decay to Barenblatt profiles. 3 J. A. Carrillo, G. Toscani, Indiana Math. Univ. J. (2000); F. Otto, Comm. PDE (2001); J. Dolbeault, M. del Pino, J. Math. Pures Appl. (2002); J.L. Vázquez, J. Evol. Eq
13 Nonlinear diffusions Nonlinear diffusions Generic finite speed of propagation/moving free boundary for degenerate diffusion equations and creation of thick tails for fast diffusion equations. No better rate of convergence can be established under the generality u 0 a probability density. At which rate does this self-similarity take over? General nonlinear diffusion equation: 8 >< >: No explicit source-type solutions, so... What is the typical asymptotic profile? u t = P(u), (x Rd, t > 0), u(x, t = 0) = u 0(x) 0 (x R d )
14 Nonlinear diffusions Nonlinear diffusions Generic finite speed of propagation/moving free boundary for degenerate diffusion equations and creation of thick tails for fast diffusion equations. No better rate of convergence can be established under the generality u 0 a probability density. At which rate does this self-similarity take over? General nonlinear diffusion equation: 8 >< >: No explicit source-type solutions, so... What is the typical asymptotic profile? u t = P(u), (x Rd, t > 0), u(x, t = 0) = u 0(x) 0 (x R d )
15 Nonlinear diffusions Nonlinear diffusions Generic finite speed of propagation/moving free boundary for degenerate diffusion equations and creation of thick tails for fast diffusion equations. No better rate of convergence can be established under the generality u 0 a probability density. At which rate does this self-similarity take over? General nonlinear diffusion equation: 8 >< >: No explicit source-type solutions, so... What is the typical asymptotic profile? u t = P(u), (x Rd, t > 0), u(x, t = 0) = u 0(x) 0 (x R d )
16 Kinetic models for granular gases Outline 1 Presentation of models Nonlinear diffusions Kinetic models for granular gases 2 Wasserstein Distance: Basics Definition Properties 3 Contractivity in 1D Diffusive Models Dissipative Models Towards any dimension
17 Kinetic models for granular gases Rapid Granular Flows Pattern formation in a vertically oscillated granular layer. 4 4 Bizon, C., Shattuck, M. D., Swift, J.B., Swinney, H.L., Phys. Rev. E (1999); Carrillo, J.A., Poschel, T., Salueña, C., in preparation (2006).
18 Kinetic models for granular gases Rapid Granular Flows Schock waves in supersonic sand. 5 5 Rericha, E., Bizon, C., Shattuck, M. D., Swinney, H.L., Phys. Rev. Letters (2002).
19 Kinetic models for granular gases Inelastic Collisions Binary Collisions: Spheres of diameter r > 0. Given (x, v) and (x rn, w), where n S 2 is the unit vector along the impact direction, the post-collisional velocities are found assuming conservation of momentum and a loose of normal relative velocity after the collision: (v w ) n = e((v w) n) where 0 < e 1 is called the restitution coefficient. Postcollisional velocities: v = 1 u (v + w) w = 1 u (v + w) 2 2 where u = u (1 + e)(u n)n, u = v w and u = v w.
20 Kinetic models for granular gases Inelastic Collisions Binary Collisions: Spheres of diameter r > 0. Given (x, v) and (x rn, w), where n S 2 is the unit vector along the impact direction, the post-collisional velocities are found assuming conservation of momentum and a loose of normal relative velocity after the collision: (v w ) n = e((v w) n) where 0 < e 1 is called the restitution coefficient. Postcollisional velocities: v = 1 u (v + w) w = 1 u (v + w) 2 2 where u = u (1 + e)(u n)n, u = v w and u = v w.
21 Kinetic models for granular gases Inelastic Collisions Postcollisional velocities: New Parameterization where v = 1 u (v + w) w = 1 u (v + w) 2 2 u = 1 e 4 u e 4 u σ, u = v w and u = v w.
22 Kinetic models for granular gases Boltzmann equation for granular gases Hypothesis: Binary, localized in t and x and inelastic collisions. Molecular chaos. Boltzmann equation for inelastic particles: f t Z Z» (v x)f =Qe(f, f )= ((v w) n) 4π R 3 S+ 2 e f 2 (v )f (w ) f (v)f (w) dndw. a Weak form of the Boltzmann equation: < ϕ, Q e(f, f ) >= 1 Z Z Z h i v w f (v)f (w) ϕ(v ) ϕ(v) dσ dv dw 4π S 2 R 3 R 3 a Jenkins, J. T., Richman, M. W., Arch. Rat. Mech. Anal. (1985).
23 Kinetic models for granular gases Boltzmann equation for granular gases Hypothesis: Binary, localized in t and x and inelastic collisions. Molecular chaos. Boltzmann equation for inelastic particles: f t Z Z» (v x)f =Qe(f, f )= ((v w) n) 4π R 3 S+ 2 e f 2 (v )f (w ) f (v)f (w) dndw. a Weak form of the Boltzmann equation: < ϕ, Q e(f, f ) >= 1 Z Z Z h i v w f (v)f (w) ϕ(v ) ϕ(v) dσ dv dw 4π S 2 R 3 R 3 a Jenkins, J. T., Richman, M. W., Arch. Rat. Mech. Anal. (1985).
24 Kinetic models for granular gases Boltzmann equation for granular gases Hypothesis: Binary, localized in t and x and inelastic collisions. Molecular chaos. Boltzmann equation for inelastic particles: f t Z Z» (v x)f =Qe(f, f )= ((v w) n) 4π R 3 S+ 2 e f 2 (v )f (w ) f (v)f (w) dndw. a Weak form of the Boltzmann equation: < ϕ, Q e(f, f ) >= 1 Z Z Z h i v w f (v)f (w) ϕ(v ) ϕ(v) dσ dv dw 4π S 2 R 3 R 3 a Jenkins, J. T., Richman, M. W., Arch. Rat. Mech. Anal. (1985).
25 Kinetic models for granular gases Boltzmann equation for granular gases Hypothesis: Binary, localized in t and x and inelastic collisions. Molecular chaos. Boltzmann equation for inelastic particles: f t Z Z» (v x)f =Qe(f, f )= ((v w) n) 4π R 3 S+ 2 e f 2 (v )f (w ) f (v)f (w) dndw. a Weak form of the Boltzmann equation: < ϕ, Q e(f, f ) >= 1 Z Z Z h i v w f (v)f (w) ϕ(v ) ϕ(v) dσ dv dw 4π S 2 R 3 R 3 a Jenkins, J. T., Richman, M. W., Arch. Rat. Mech. Anal. (1985).
26 Kinetic models for granular gases Basic Properties 6 f t = Qe(f, f ) Mass and momentum are preserved while energy decreases: v 2 + w 2 v 2 w 2 = 1 e2 u 2 2 Let us fix unit mass and zero mean velocity for the rest. Temperature cools off down to 0 as t (Haff s law): dθ dt 1 e2 4 What are the typical asymptotic profiles? Do self-similar solutions exist? If so At which rate does this self-similarity take over? θ 3 2, 6 J.A. Carrillo, C. Cercignani, I.M. Gamba, Phys. Rev. E (2000); C. Cercignani, R. Illner, C. Stoica, J. Stat. Phys. (2001); A. Bobylev, C. Cercignani, J. Stat. Phys. (2002); I.M. Gamba, V. Panferov, C. Villani, Comm. Math. Phys. (2004); C. Mouhot, S. Mischer, Martínez-Ricard, J. Stat. Phys. (2006).
27 Kinetic models for granular gases Basic Properties 6 f t = Qe(f, f ) Mass and momentum are preserved while energy decreases: v 2 + w 2 v 2 w 2 = 1 e2 u 2 2 Let us fix unit mass and zero mean velocity for the rest. Temperature cools off down to 0 as t (Haff s law): dθ dt 1 e2 4 What are the typical asymptotic profiles? Do self-similar solutions exist? If so At which rate does this self-similarity take over? θ 3 2, 6 J.A. Carrillo, C. Cercignani, I.M. Gamba, Phys. Rev. E (2000); C. Cercignani, R. Illner, C. Stoica, J. Stat. Phys. (2001); A. Bobylev, C. Cercignani, J. Stat. Phys. (2002); I.M. Gamba, V. Panferov, C. Villani, Comm. Math. Phys. (2004); C. Mouhot, S. Mischer, Martínez-Ricard, J. Stat. Phys. (2006).
28 Kinetic models for granular gases Simplified Granular media model Another dissipative models (toy models) keeping the same properties: conservation of mass, mean velocity and dissipation of energy are of the form: where f + (v x)f = I(f, f ) t I(f, f ) = div v (v w) v w f (x, w, t) dw f (x, v, t).»zr 3 This dissipative collision in the one dimensional case simplifies to 7 : I(f, f ) =»Z (v w) v w f (x, w, t) dw f (x, v, t) v R 7 Benedetto, D., Caglioti, E., Pulvirenti, M., Math. Mod. and Num. An. (1997).
29 Kinetic models for granular gases Simplified Granular media model Another dissipative models (toy models) keeping the same properties: conservation of mass, mean velocity and dissipation of energy are of the form: where f + (v x)f = I(f, f ) t I(f, f ) = div v (v w) v w f (x, w, t) dw f (x, v, t).»zr 3 This dissipative collision in the one dimensional case simplifies to 7 : I(f, f ) =»Z (v w) v w f (x, w, t) dw f (x, v, t) v R 7 Benedetto, D., Caglioti, E., Pulvirenti, M., Math. Mod. and Num. An. (1997).
30 Kinetic models for granular gases Simplified Granular media model Another dissipative models (toy models) keeping the same properties: conservation of mass, mean velocity and dissipation of energy are of the form: where f + (v x)f = I(f, f ) t I(f, f ) = div v (v w) v w f (x, w, t) dw f (x, v, t).»zr 3 This dissipative collision in the one dimensional case simplifies to 7 : I(f, f ) =»Z (v w) v w f (x, w, t) dw f (x, v, t) v R 7 Benedetto, D., Caglioti, E., Pulvirenti, M., Math. Mod. and Num. An. (1997).
31 Kinetic models for granular gases Summary: Spreading versus Concentration Diffusive models - Spreading: u t = P(u), (x Rd, t > 0) Typical asymptotic profiles for P(u) = u m (Barenblatt): with θ m(t). B m(x, t) = θ m(t) d/2 B m(θ m(t) 1/2 x, t o,m) Homogeneous Kinetic Dissipative models - Concentration: 8 >< Q e(f, f ) f t = Z >: div v»f (v, t) W(v w)f (w, t) dw R d Typical asymptotic profiles (Homogeneous Cooling States), if they exist: with θ hc(t) 0. f hc(v, t) = ρ θ 2 3 hc (t) g ((v u) θ 2 1 hc (t))
32 Kinetic models for granular gases Summary: Spreading versus Concentration Diffusive models - Spreading: u t = P(u), (x Rd, t > 0) Typical asymptotic profiles for P(u) = u m (Barenblatt): with θ m(t). B m(x, t) = θ m(t) d/2 B m(θ m(t) 1/2 x, t o,m) Homogeneous Kinetic Dissipative models - Concentration: 8 >< Q e(f, f ) f t = Z >: div v»f (v, t) W(v w)f (w, t) dw R d Typical asymptotic profiles (Homogeneous Cooling States), if they exist: with θ hc(t) 0. f hc(v, t) = ρ θ 2 3 hc (t) g ((v u) θ 2 1 hc (t))
33 Definition Outline 1 Presentation of models Nonlinear diffusions Kinetic models for granular gases 2 Wasserstein Distance: Basics Definition Properties 3 Contractivity in 1D Diffusive Models Dissipative Models Towards any dimension
34 Definition Definition of the distance 8 Transporting measures: Given T : R d R d mesurable, we say that T transports µ P onto ν P or that ν is the push-forward of µ through T, ν = T#µ, if ν[k] := µ[t 1 (K)] for all mesurable sets K R d, equivalently Z Z ϕ dν = (ϕ T) dµ R d R d for all ϕ C o(r d ). Random variables: Say that X is a random variable with law given by µ, is to say X : (Ω, A, P) (R d, B d) is a mesurable map such that X#P = µ, i.e., Z Z ϕ(x) dµ = (ϕ X) dp = E [ϕ(x)]. R d Ω Coming back to the transport of measures: if X and Y are random variables with law µ and ν respectively, then ν = T#µ is equivalent to Y = T(X). 8 C. Villani, AMS Graduate Texts (2003).
35 Definition Definition of the distance 8 Transporting measures: Given T : R d R d mesurable, we say that T transports µ P onto ν P or that ν is the push-forward of µ through T, ν = T#µ, if ν[k] := µ[t 1 (K)] for all mesurable sets K R d, equivalently Z Z ϕ dν = (ϕ T) dµ R d R d for all ϕ C o(r d ). Random variables: Say that X is a random variable with law given by µ, is to say X : (Ω, A, P) (R d, B d) is a mesurable map such that X#P = µ, i.e., Z Z ϕ(x) dµ = (ϕ X) dp = E [ϕ(x)]. R d Ω Coming back to the transport of measures: if X and Y are random variables with law µ and ν respectively, then ν = T#µ is equivalent to Y = T(X). 8 C. Villani, AMS Graduate Texts (2003).
36 Definition Definition of the distance 8 Transporting measures: Given T : R d R d mesurable, we say that T transports µ P onto ν P or that ν is the push-forward of µ through T, ν = T#µ, if ν[k] := µ[t 1 (K)] for all mesurable sets K R d, equivalently Z Z ϕ dν = (ϕ T) dµ R d R d for all ϕ C o(r d ). Random variables: Say that X is a random variable with law given by µ, is to say X : (Ω, A, P) (R d, B d) is a mesurable map such that X#P = µ, i.e., Z Z ϕ(x) dµ = (ϕ X) dp = E [ϕ(x)]. R d Ω Coming back to the transport of measures: if X and Y are random variables with law µ and ν respectively, then ν = T#µ is equivalent to Y = T(X). 8 C. Villani, AMS Graduate Texts (2003).
37 Definition Definition of the distance 8 Transporting measures: Given T : R d R d mesurable, we say that T transports µ P onto ν P or that ν is the push-forward of µ through T, ν = T#µ, if ν[k] := µ[t 1 (K)] for all mesurable sets K R d, equivalently Z Z ϕ dν = (ϕ T) dµ R d R d for all ϕ C o(r d ). Random variables: Say that X is a random variable with law given by µ, is to say X : (Ω, A, P) (R d, B d) is a mesurable map such that X#P = µ, i.e., Z Z ϕ(x) dµ = (ϕ X) dp = E [ϕ(x)]. R d Ω Coming back to the transport of measures: if X and Y are random variables with law µ and ν respectively, then ν = T#µ is equivalent to Y = T(X). 8 C. Villani, AMS Graduate Texts (2003).
38 Definition Definition of the distance Euclidean Wasserstein jzz Distance: ff 1/2 n W 2(µ, ν)= inf x y 2 dπ(x, y) = inf E h X Y 2io 1/2 π R d R d (X,Y) where the transference plan π runs over the set of joint probability measures on R d R d with marginals f and g P 2(R d ) and (X, Y) are all possible couples of random variables with µ and ν as respective laws. Monge s optimal mass transport problem: Find jz ff 1/2 I := inf x T(x) 2 dµ(x); ν = T#µ. T R d Take γ T = (1 R d T)#µ as transference plan π.
39 Definition Definition of the distance Euclidean Wasserstein jzz Distance: ff 1/2 n W 2(µ, ν)= inf x y 2 dπ(x, y) = inf E h X Y 2io 1/2 π R d R d (X,Y) where the transference plan π runs over the set of joint probability measures on R d R d with marginals f and g P 2(R d ) and (X, Y) are all possible couples of random variables with µ and ν as respective laws. Monge s optimal mass transport problem: Find jz ff 1/2 I := inf x T(x) 2 dµ(x); ν = T#µ. T R d Take γ T = (1 R d T)#µ as transference plan π.
40 Definition Definition of the distance Euclidean Wasserstein jzz Distance: ff 1/2 n W 2(µ, ν)= inf x y 2 dπ(x, y) = inf E h X Y 2io 1/2 π R d R d (X,Y) where the transference plan π runs over the set of joint probability measures on R d R d with marginals f and g P 2(R d ) and (X, Y) are all possible couples of random variables with µ and ν as respective laws. Monge s optimal mass transport problem: Find jz ff 1/2 I := inf x T(x) 2 dµ(x); ν = T#µ. T R d Take γ T = (1 R d T)#µ as transference plan π.
41 Definition Definition of the distance Wasserstein Distances: Given 1 p <, we define jzz ff 1/p W p(µ, ν)= inf x y p dπ(x, y) = inf {E [ X π R d R d (X,Y) Y p ]} 1/p where the transference plan π runs over the set of joint probability measures on R d R d with marginals f and g P p(r d ) and (X, Y) are all possible couples of random variables with µ and ν as respective laws. We define the -Wasserstein distance as W (µ, ν) = lim Wp(µ, ν) = inf {esssup { x T(x), x supp(µ)} with ν = T#µ} p T for measures with all moments bounded.
42 Definition Definition of the distance Wasserstein Distances: Given 1 p <, we define jzz ff 1/p W p(µ, ν)= inf x y p dπ(x, y) = inf {E [ X π R d R d (X,Y) Y p ]} 1/p where the transference plan π runs over the set of joint probability measures on R d R d with marginals f and g P p(r d ) and (X, Y) are all possible couples of random variables with µ and ν as respective laws. We define the -Wasserstein distance as W (µ, ν) = lim Wp(µ, ν) = inf {esssup { x T(x), x supp(µ)} with ν = T#µ} p T for measures with all moments bounded.
43 Definition Definition of the distance Wasserstein Distances: Given 1 p <, we define jzz ff 1/p W p(µ, ν)= inf x y p dπ(x, y) = inf {E [ X π R d R d (X,Y) Y p ]} 1/p where the transference plan π runs over the set of joint probability measures on R d R d with marginals f and g P p(r d ) and (X, Y) are all possible couples of random variables with µ and ν as respective laws. We define the -Wasserstein distance as W (µ, ν) = lim Wp(µ, ν) = inf {esssup { x T(x), x supp(µ)} with ν = T#µ} p T for measures with all moments bounded.
44 Properties Outline 1 Presentation of models Nonlinear diffusions Kinetic models for granular gases 2 Wasserstein Distance: Basics Definition Properties 3 Contractivity in 1D Diffusive Models Dissipative Models Towards any dimension
45 Properties Euclidean Wasserstein Distance Basic Properties 1 Convexity: f 1, f 2, g 1, g 2 P 2(R d ) and α [0, 1], then W 2 2 (αf 1 + (1 α)f 2, αg 1 + (1 α)g 2) αw 2 2 (f 1, g 1) + (1 α)w 2 2 (f 2, g 2). 2 Relation to Temperature: f P 2(R d ) and a R 3, then Z W2 2 (f, δ a) = v a 2 df (v). R d 3 Scaling: Given f in P 2(R d ) and θ > 0, let us define S θ [f ] = θ d/2 f (θ 1/2 v) for absolutely continuous measures with respect to Lebesgue or its corresponding definition by duality for general measures; then for any f and g in P 2(R d ), we have W 2(S θ [f ], S θ [g]) = θ 1/2 W 2(f, g).
46 Properties Euclidean Wasserstein Distance Basic Properties 1 Convexity: f 1, f 2, g 1, g 2 P 2(R d ) and α [0, 1], then W 2 2 (αf 1 + (1 α)f 2, αg 1 + (1 α)g 2) αw 2 2 (f 1, g 1) + (1 α)w 2 2 (f 2, g 2). 2 Relation to Temperature: f P 2(R d ) and a R 3, then Z W2 2 (f, δ a) = v a 2 df (v). R d 3 Scaling: Given f in P 2(R d ) and θ > 0, let us define S θ [f ] = θ d/2 f (θ 1/2 v) for absolutely continuous measures with respect to Lebesgue or its corresponding definition by duality for general measures; then for any f and g in P 2(R d ), we have W 2(S θ [f ], S θ [g]) = θ 1/2 W 2(f, g).
47 Properties Euclidean Wasserstein Distance Basic Properties 1 Convexity: f 1, f 2, g 1, g 2 P 2(R d ) and α [0, 1], then W 2 2 (αf 1 + (1 α)f 2, αg 1 + (1 α)g 2) αw 2 2 (f 1, g 1) + (1 α)w 2 2 (f 2, g 2). 2 Relation to Temperature: f P 2(R d ) and a R 3, then Z W2 2 (f, δ a) = v a 2 df (v). R d 3 Scaling: Given f in P 2(R d ) and θ > 0, let us define S θ [f ] = θ d/2 f (θ 1/2 v) for absolutely continuous measures with respect to Lebesgue or its corresponding definition by duality for general measures; then for any f and g in P 2(R d ), we have W 2(S θ [f ], S θ [g]) = θ 1/2 W 2(f, g).
48 Properties Euclidean Wasserstein Distance Convergence Properties 1 Convergence of measures: W 2(f n, f ) 0 is equivalent to f n f weakly-* as measures and convergence of second moments. 2 Weak lower semicontinuity: Given f n f and g n g weakly-* as measures, then W 2(f, g) lim inf W2(fn, gn). n 3 Completeness: The space P 2(R d ) endowed with the distance W 2 is a complete metric space. Moreover, the set j Z ff M θ = µ P 2(R d ) such that v 2 df (v) = θ R 3, i.e. the sphere of radius θ in P 2(R d ), is a complete metric space.
49 Properties Euclidean Wasserstein Distance Convergence Properties 1 Convergence of measures: W 2(f n, f ) 0 is equivalent to f n f weakly-* as measures and convergence of second moments. 2 Weak lower semicontinuity: Given f n f and g n g weakly-* as measures, then W 2(f, g) lim inf W2(fn, gn). n 3 Completeness: The space P 2(R d ) endowed with the distance W 2 is a complete metric space. Moreover, the set j Z ff M θ = µ P 2(R d ) such that v 2 df (v) = θ R 3, i.e. the sphere of radius θ in P 2(R d ), is a complete metric space.
50 Properties Euclidean Wasserstein Distance Convergence Properties 1 Convergence of measures: W 2(f n, f ) 0 is equivalent to f n f weakly-* as measures and convergence of second moments. 2 Weak lower semicontinuity: Given f n f and g n g weakly-* as measures, then W 2(f, g) lim inf W2(fn, gn). n 3 Completeness: The space P 2(R d ) endowed with the distance W 2 is a complete metric space. Moreover, the set j Z ff M θ = µ P 2(R d ) such that v 2 df (v) = θ R 3, i.e. the sphere of radius θ in P 2(R d ), is a complete metric space.
51 Properties Wasserstein Distances Properties 1 Convexity of W p p 2 Scaling: for any f and g in P p(r d ), we have W p(s θ [f ], S θ [g]) = θ 1/2 W p(f, g). 3 Completeness and weak lower semicontinuity. 4 Convergence weak-* convergence as measures and convergence of p-moments.
52 Properties Wasserstein Distances Properties 1 Convexity of W p p 2 Scaling: for any f and g in P p(r d ), we have W p(s θ [f ], S θ [g]) = θ 1/2 W p(f, g). 3 Completeness and weak lower semicontinuity. 4 Convergence weak-* convergence as measures and convergence of p-moments.
53 Properties Wasserstein Distances Properties 1 Convexity of W p p 2 Scaling: for any f and g in P p(r d ), we have W p(s θ [f ], S θ [g]) = θ 1/2 W p(f, g). 3 Completeness and weak lower semicontinuity. 4 Convergence weak-* convergence as measures and convergence of p-moments.
54 Properties Wasserstein Distances Properties 1 Convexity of W p p 2 Scaling: for any f and g in P p(r d ), we have W p(s θ [f ], S θ [g]) = θ 1/2 W p(f, g). 3 Completeness and weak lower semicontinuity. 4 Convergence weak-* convergence as measures and convergence of p-moments.
55 Properties One dimensional Case Distribution functions: In one dimension, denoting by F(x) the distribution function of µ, F(x) = we can define its pseudo inverse: F 1 (η) = inf{x : F(x) > η} for η (0, 1), we have F 1 : ((0, 1), B 1), dη) (R, B 1) is a random variable with law µ, i.e., F 1 #dη = µ Z Z 1 ϕ(x) dµ = ϕ(f 1 (η)) dη = E [ϕ(x)]. R 0 Z x dµ,
56 Properties One dimensional Case Distribution functions: In one dimension, denoting by F(x) the distribution function of µ, F(x) = we can define its pseudo inverse: F 1 (η) = inf{x : F(x) > η} for η (0, 1), we have F 1 : ((0, 1), B 1), dη) (R, B 1) is a random variable with law µ, i.e., F 1 #dη = µ Z Z 1 ϕ(x) dµ = ϕ(f 1 (η)) dη = E [ϕ(x)]. R 0 Z x dµ,
57 Properties One dimensional Case Wasserstein distance: Given µ and ν P 2(R d ) with distribution functions F(x) and G(y) then, ZZ W2 2 (µ, ν) = x y 2 dh(x, y) R d R d where H(x, y) = min(f(x), G(y)) is a joint distribution function and the above integral is defined in the Riemann-Stieljes sense. a In one dimension, it is easy to check that given two measures µ and ν with distribution functions F(x) and G(y) then, (F 1 G 1 )#dη has joint distribution function H(x, y) = min(f(x), G(y)). Therefore, in one dimension, the optimal plan is given by π opt(x, y) = (F 1 G 1 )#dη, and thus Z 1 W 2(µ, ν) = [F 1 (η) G 1 (η)] 2 dη 0 «1/2 a W. Hoeffding (1940); M. Fréchet (1951); A. Pulvirenti, G. Toscani, Annali Mat. Pura Appl. (1996).
58 Properties One dimensional Case Wasserstein distance: Given µ and ν P 2(R d ) with distribution functions F(x) and G(y) then, ZZ W2 2 (µ, ν) = x y 2 dh(x, y) R d R d where H(x, y) = min(f(x), G(y)) is a joint distribution function and the above integral is defined in the Riemann-Stieljes sense. a In one dimension, it is easy to check that given two measures µ and ν with distribution functions F(x) and G(y) then, (F 1 G 1 )#dη has joint distribution function H(x, y) = min(f(x), G(y)). Therefore, in one dimension, the optimal plan is given by π opt(x, y) = (F 1 G 1 )#dη, and thus Z 1 W 2(µ, ν) = [F 1 (η) G 1 (η)] 2 dη 0 «1/2 a W. Hoeffding (1940); M. Fréchet (1951); A. Pulvirenti, G. Toscani, Annali Mat. Pura Appl. (1996).
59 Properties One dimensional Case Wasserstein distance: Given µ and ν P p(r), 1 p < then and Z 1 1/p W p(µ, ν) = [F 1 (η) G 1 (η)] dη«p = F 1 G 1 L p (R) 0 W (µ, ν) = F 1 G 1 L (R).
60 Properties Relation 1D-anyD Euclidean Wasserstein distance: If the probability measures µ i j on R are the successive one-dimensional marginals of the measure µ i on R d, for i = 1, 2 and j = 1,..., d, then dx W2 2 (µ 1 j, µ 2 j ) W2 2 (µ 1, µ 2 ). j=1
61 Properties Relation 1D-anyD Proof. Let π be a measure on R d x R d y with marginals µ 1 and µ 2, optimal in the sense that ZZ W2 2 (µ 1, µ 2 ) = x y 2 dπ(x, y). R d R d Then its marginal π j on R vj R wj has itself marginals µ 1 j and µ 2 j, so ZZ W2 2 (µ 1 j, µ 2 j ) v j w j 2 dπ j(v j, w j). R R The inequality follows by noting that dx v j w j 2 = v w 2. j=1
62 Properties Relation 1D-anyD Proof. Let π be a measure on R d x R d y with marginals µ 1 and µ 2, optimal in the sense that ZZ W2 2 (µ 1, µ 2 ) = x y 2 dπ(x, y). R d R d Then its marginal π j on R vj R wj has itself marginals µ 1 j and µ 2 j, so ZZ W2 2 (µ 1 j, µ 2 j ) v j w j 2 dπ j(v j, w j). R R The inequality follows by noting that dx v j w j 2 = v w 2. j=1
63 Properties Relation 1D-anyD Proof. Let π be a measure on R d x R d y with marginals µ 1 and µ 2, optimal in the sense that ZZ W2 2 (µ 1, µ 2 ) = x y 2 dπ(x, y). R d R d Then its marginal π j on R vj R wj has itself marginals µ 1 j and µ 2 j, so ZZ W2 2 (µ 1 j, µ 2 j ) v j w j 2 dπ j(v j, w j). R R The inequality follows by noting that dx v j w j 2 = v w 2. j=1
64 Diffusive Models Outline 1 Presentation of models Nonlinear diffusions Kinetic models for granular gases 2 Wasserstein Distance: Basics Definition Properties 3 Contractivity in 1D Diffusive Models Dissipative Models Towards any dimension
65 Diffusive Models 1D Nonlinear Diffusions 9 The flows of nonlinear diffusions are contractions in Wasserstein distance. u t = (P(u)) xx u(x, 0) = u 0(x) 0 L 1 (R) with P : R + 0 R, continuous, nondecreasing with P(0) = 0. Diccionary for F 1 (η): F 1 (η, t) = η F 1 t» 1 u (η, t) = (F 1 (η, t), t),» 1 (F 1 (η, t), t) F u t (F 1 (η, t), t) 2 F 1 h ux i (η, t) = (F 1 (η, t), t) η 2 u 3 Equation for F 1 : F 1 (η, t) =» P t η» F 1 η «1, 9 H. Li, G. Toscani, ARMA (2004);J.A. Carrillo, G. Toscani, proceedings conference in honor of S. Rionero (2005).
66 Diffusive Models 1D Nonlinear Diffusions 9 The flows of nonlinear diffusions are contractions in Wasserstein distance. u t = (P(u)) xx u(x, 0) = u 0(x) 0 L 1 (R) with P : R + 0 R, continuous, nondecreasing with P(0) = 0. Diccionary for F 1 (η): F 1 (η, t) = η F 1 t» 1 u (η, t) = (F 1 (η, t), t),» 1 (F 1 (η, t), t) F u t (F 1 (η, t), t) 2 F 1 h ux i (η, t) = (F 1 (η, t), t) η 2 u 3 Equation for F 1 : F 1 (η, t) =» P t η» F 1 η «1, 9 H. Li, G. Toscani, ARMA (2004);J.A. Carrillo, G. Toscani, proceedings conference in honor of S. Rionero (2005).
67 Diffusive Models 1D Nonlinear Diffusions 9 The flows of nonlinear diffusions are contractions in Wasserstein distance. u t = (P(u)) xx u(x, 0) = u 0(x) 0 L 1 (R) with P : R + 0 R, continuous, nondecreasing with P(0) = 0. Diccionary for F 1 (η): F 1 (η, t) = η F 1 t» 1 u (η, t) = (F 1 (η, t), t),» 1 (F 1 (η, t), t) F u t (F 1 (η, t), t) 2 F 1 h ux i (η, t) = (F 1 (η, t), t) η 2 u 3 Equation for F 1 : F 1 (η, t) =» P t η» F 1 η «1, 9 H. Li, G. Toscani, ARMA (2004);J.A. Carrillo, G. Toscani, proceedings conference in honor of S. Rionero (2005).
68 Diffusive Models 1D Nonlinear Diffusions 10 Contraction for nonlinear diffusions: Given any two solutions of the nonlinear diffusion equation W p(u 1(t), u 2(t)) W p(u 1(0), u 2(0)) for all 1 p. Proof. Assume 1 < p <, the other ones by approximation in p, then 1 d p(p 1) dt Z 1 0 Z 1 0 F 1 G 1 p 2 (F 1 η F 1 G 1 p dη = j G 1 η ) P» 1 Fη 1 P» ff 1 G 1 η dη 0 10 H. Li, G. Toscani, ARMA (2004);J.A. Carrillo, G. Toscani, proceedings conference in honor of S. Rionero (2005).
69 Diffusive Models 1D Nonlinear Diffusions 10 Contraction for nonlinear diffusions: Given any two solutions of the nonlinear diffusion equation W p(u 1(t), u 2(t)) W p(u 1(0), u 2(0)) for all 1 p. Proof. Assume 1 < p <, the other ones by approximation in p, then 1 d p(p 1) dt Z 1 0 Z 1 0 F 1 G 1 p 2 (F 1 η F 1 G 1 p dη = j G 1 η ) P» 1 Fη 1 P» ff 1 G 1 η dη 0 10 H. Li, G. Toscani, ARMA (2004);J.A. Carrillo, G. Toscani, proceedings conference in honor of S. Rionero (2005).
70 Diffusive Models 1D Nonlinear Diffusions 11 Strict Contraction for nonlinear diffusions with confining potential: Given any two probability density solutions of the nonlinear diffusion equation «dv ρ t = (t, x)ρ + (P(ρ)) dx xx where V(t, x) is a smooth strictly uniform convex potential ( d2 V dx 2 (t, x) α > 0), then for all 1 p. W p(ρ 1(t), ρ 2(t)) W p(ρ 1(0), ρ 2(0)) e αt x Proof. Exercise; F 1 (η, t) = dv t dx (t, F 1 )» P η» F 1 η «1, The case of V(x) = x 2 /2 and P(ρ) = ρ m comes by the scaling property of the Wasserstein distances and the previous non strict contraction. 11 H. Li, G. Toscani, ARMA (2004);J.A. Carrillo, G. Toscani, proceedings conference in honor of S. Rionero (2005).
71 Diffusive Models 1D Nonlinear Diffusions 11 Strict Contraction for nonlinear diffusions with confining potential: Given any two probability density solutions of the nonlinear diffusion equation «dv ρ t = (t, x)ρ + (P(ρ)) dx xx where V(t, x) is a smooth strictly uniform convex potential ( d2 V dx 2 (t, x) α > 0), then for all 1 p. W p(ρ 1(t), ρ 2(t)) W p(ρ 1(0), ρ 2(0)) e αt x Proof. Exercise; F 1 (η, t) = dv t dx (t, F 1 )» P η» F 1 η «1, The case of V(x) = x 2 /2 and P(ρ) = ρ m comes by the scaling property of the Wasserstein distances and the previous non strict contraction. 11 H. Li, G. Toscani, ARMA (2004);J.A. Carrillo, G. Toscani, proceedings conference in honor of S. Rionero (2005).
72 Diffusive Models 1D Nonlinear Diffusions Convergence towards equilibrium: All probability solutions of the equation «dv ρ t = (t, x)ρ + (P(ρ)) dx xx converges exponentially fast towards the unique equilibrium ρ characterized by x V(x) + h(ρ (x)) = C ρh (ρ) = P (ρ). Here, P is strictly increasing for being h invertible.
73 Diffusive Models 1D Nonlinear Diffusions 12 Uniform Finite Speed of Propagation: Given any two probability density solutions of the degenerate nonlinear diffusion equation ρ t = (ρ m ) xx, we have the following control of the "relative speed" of the support spreading for any t > 0. inf{supp{ρ 1(t)}} inf{supp{ρ 2(t)}} W (ρ 1(0), ρ 2(0)) sup{supp{ρ 1(t)}} sup{supp{ρ 2(t)}} W (ρ 1(0), ρ 2(0)) 12 J.A. Carrillo, M.P. Gualdani, G. Toscani, CRAS (2004).
74 Dissipative Models Outline 1 Presentation of models Nonlinear diffusions Kinetic models for granular gases 2 Wasserstein Distance: Basics Definition Properties 3 Contractivity in 1D Diffusive Models Dissipative Models Towards any dimension
75 Dissipative Models 1D Dissipative Kinetic Equations 13 Given the 1D granular model f t =»Z (v w) v w γ f (w, t) dw f (v, t) v R Equation for F 1 : F 1 (η, t) = t Z 1 After some manipulations (good exercise) 0 F 1 (η, t) F 1 (p, t) γ (F 1 (η, t) F 1 (p, t)) dp. implying d 1 W2(f (t), g(t)) W2(f (t), g(t))1+γ dt 2γ 1 W 2(f (t), g(t)) 2 γ 1 «1/γ. γt + 2 γ 1 W 2(f (0), g(0)) γ 13 H. Li, G. Toscani, ARMA (2004);J.A. Carrillo, G. Toscani, proceedings conference in honor of S. Rionero (2005).
76 Dissipative Models 1D Dissipative Kinetic Equations 13 Given the 1D granular model f t =»Z (v w) v w γ f (w, t) dw f (v, t) v R Equation for F 1 : F 1 (η, t) = t Z 1 After some manipulations (good exercise) 0 F 1 (η, t) F 1 (p, t) γ (F 1 (η, t) F 1 (p, t)) dp. implying d 1 W2(f (t), g(t)) W2(f (t), g(t))1+γ dt 2γ 1 W 2(f (t), g(t)) 2 γ 1 «1/γ. γt + 2 γ 1 W 2(f (0), g(0)) γ 13 H. Li, G. Toscani, ARMA (2004);J.A. Carrillo, G. Toscani, proceedings conference in honor of S. Rionero (2005).
77 Dissipative Models 1D Dissipative Kinetic Equations 13 Given the 1D granular model f t =»Z (v w) v w γ f (w, t) dw f (v, t) v R Equation for F 1 : F 1 (η, t) = t Z 1 After some manipulations (good exercise) 0 F 1 (η, t) F 1 (p, t) γ (F 1 (η, t) F 1 (p, t)) dp. implying d 1 W2(f (t), g(t)) W2(f (t), g(t))1+γ dt 2γ 1 W 2(f (t), g(t)) 2 γ 1 «1/γ. γt + 2 γ 1 W 2(f (0), g(0)) γ 13 H. Li, G. Toscani, ARMA (2004);J.A. Carrillo, G. Toscani, proceedings conference in honor of S. Rionero (2005).
78 Dissipative Models 1D Dissipative Kinetic Equations 13 Given the 1D granular model f t =»Z (v w) v w γ f (w, t) dw f (v, t) v R Equation for F 1 : F 1 (η, t) = t Z 1 After some manipulations (good exercise) 0 F 1 (η, t) F 1 (p, t) γ (F 1 (η, t) F 1 (p, t)) dp. implying d 1 W2(f (t), g(t)) W2(f (t), g(t))1+γ dt 2γ 1 W 2(f (t), g(t)) 2 γ 1 «1/γ. γt + 2 γ 1 W 2(f (0), g(0)) γ 13 H. Li, G. Toscani, ARMA (2004);J.A. Carrillo, G. Toscani, proceedings conference in honor of S. Rionero (2005).
79 Towards any dimension Outline 1 Presentation of models Nonlinear diffusions Kinetic models for granular gases 2 Wasserstein Distance: Basics Definition Properties 3 Contractivity in 1D Diffusive Models Dissipative Models Towards any dimension
80 Towards any dimension Heat Equation 14 Contraction for the Heat Equation: Given any two probability density solutions u 1(x, t), u 2(x, t) of the heat equation, then W p(u 1(t), u 2(t)) W p(u 1(0), u 2(0)) for t 0, 1 p. Proof. Let γ o the optimal transference plan between u 1(0) y u 2(0) for W p, 1 p <. The solutions are u 1(t) = K(t, x) u 1(0) and u 2(t) = K(t, x) u 2(0). Let uszdefine γ t by duality: Z Z ϕ(x, y) dγ t(x, y) = ϕ(x + z, y + z)k(t, z) dz dγ o(x, y) R d R d R d R d R d Check that γ t is a transference plan between u 1(t) and u 2(t), and thus Wp p (u 1(t), u 2(t)) Z x y p dγ t(x, y)= R d R d Z x y p dγ o(x, y)=wp p (u 1(0), u 2(0)). R d R d 14 J.A. Carrillo, Boletín SEMA (2004).
81 Towards any dimension Heat Equation 14 Contraction for the Heat Equation: Given any two probability density solutions u 1(x, t), u 2(x, t) of the heat equation, then W p(u 1(t), u 2(t)) W p(u 1(0), u 2(0)) for t 0, 1 p. Proof. Let γ o the optimal transference plan between u 1(0) y u 2(0) for W p, 1 p <. The solutions are u 1(t) = K(t, x) u 1(0) and u 2(t) = K(t, x) u 2(0). Let uszdefine γ t by duality: Z Z ϕ(x, y) dγ t(x, y) = ϕ(x + z, y + z)k(t, z) dz dγ o(x, y) R d R d R d R d R d Check that γ t is a transference plan between u 1(t) and u 2(t), and thus Wp p (u 1(t), u 2(t)) Z x y p dγ t(x, y)= R d R d Z x y p dγ o(x, y)=wp p (u 1(0), u 2(0)). R d R d 14 J.A. Carrillo, Boletín SEMA (2004).
82 Towards any dimension Heat Equation 14 Contraction for the Heat Equation: Given any two probability density solutions u 1(x, t), u 2(x, t) of the heat equation, then W p(u 1(t), u 2(t)) W p(u 1(0), u 2(0)) for t 0, 1 p. Proof. Let γ o the optimal transference plan between u 1(0) y u 2(0) for W p, 1 p <. The solutions are u 1(t) = K(t, x) u 1(0) and u 2(t) = K(t, x) u 2(0). Let uszdefine γ t by duality: Z Z ϕ(x, y) dγ t(x, y) = ϕ(x + z, y + z)k(t, z) dz dγ o(x, y) R d R d R d R d R d Check that γ t is a transference plan between u 1(t) and u 2(t), and thus Wp p (u 1(t), u 2(t)) Z x y p dγ t(x, y)= R d R d Z x y p dγ o(x, y)=wp p (u 1(0), u 2(0)). R d R d 14 J.A. Carrillo, Boletín SEMA (2004).
83 Towards any dimension Heat Equation 14 Contraction for the Heat Equation: Given any two probability density solutions u 1(x, t), u 2(x, t) of the heat equation, then W p(u 1(t), u 2(t)) W p(u 1(0), u 2(0)) for t 0, 1 p. Proof. Let γ o the optimal transference plan between u 1(0) y u 2(0) for W p, 1 p <. The solutions are u 1(t) = K(t, x) u 1(0) and u 2(t) = K(t, x) u 2(0). Let uszdefine γ t by duality: Z Z ϕ(x, y) dγ t(x, y) = ϕ(x + z, y + z)k(t, z) dz dγ o(x, y) R d R d R d R d R d Check that γ t is a transference plan between u 1(t) and u 2(t), and thus Wp p (u 1(t), u 2(t)) Z x y p dγ t(x, y)= R d R d Z x y p dγ o(x, y)=wp p (u 1(0), u 2(0)). R d R d 14 J.A. Carrillo, Boletín SEMA (2004).
84 Towards any dimension Heat Equation 14 Contraction for the Heat Equation: Given any two probability density solutions u 1(x, t), u 2(x, t) of the heat equation, then W p(u 1(t), u 2(t)) W p(u 1(0), u 2(0)) for t 0, 1 p. Proof. Let γ o the optimal transference plan between u 1(0) y u 2(0) for W p, 1 p <. The solutions are u 1(t) = K(t, x) u 1(0) and u 2(t) = K(t, x) u 2(0). Let uszdefine γ t by duality: Z Z ϕ(x, y) dγ t(x, y) = ϕ(x + z, y + z)k(t, z) dz dγ o(x, y) R d R d R d R d R d Check that γ t is a transference plan between u 1(t) and u 2(t), and thus Wp p (u 1(t), u 2(t)) Z x y p dγ t(x, y)= R d R d Z x y p dγ o(x, y)=wp p (u 1(0), u 2(0)). R d R d 14 J.A. Carrillo, Boletín SEMA (2004).
85 Towards any dimension Homogeneous Boltzmann Equation 15 Convergence Towards Monokinetic State: Given f (v, t) P 2(R 3 ) with zero mean velocity then para t 0. Proof. Evolution of temperature: < v 2, Q e(f, f ) >= 1 e2 4 W 2(f (t), δ 0) C(1 + t) 2 Z R 3 Z R 3 v w 3 f (v)f (w) dv dw and thus using Jensen s and Hölder s inequality, we get the Haff s law: dθ dt 1 e2 4 Using the relation to temperature of the Euclidean Wasserstein distance, we conclude. θ J.A. Carrillo, Boletín SEMA (2004).
86 Towards any dimension Homogeneous Boltzmann Equation 15 Convergence Towards Monokinetic State: Given f (v, t) P 2(R 3 ) with zero mean velocity then para t 0. Proof. Evolution of temperature: < v 2, Q e(f, f ) >= 1 e2 4 W 2(f (t), δ 0) C(1 + t) 2 Z R 3 Z R 3 v w 3 f (v)f (w) dv dw and thus using Jensen s and Hölder s inequality, we get the Haff s law: dθ dt 1 e2 4 Using the relation to temperature of the Euclidean Wasserstein distance, we conclude. θ J.A. Carrillo, Boletín SEMA (2004).
87 Towards any dimension Homogeneous Boltzmann Equation 15 Convergence Towards Monokinetic State: Given f (v, t) P 2(R 3 ) with zero mean velocity then para t 0. Proof. Evolution of temperature: < v 2, Q e(f, f ) >= 1 e2 4 W 2(f (t), δ 0) C(1 + t) 2 Z R 3 Z R 3 v w 3 f (v)f (w) dv dw and thus using Jensen s and Hölder s inequality, we get the Haff s law: dθ dt 1 e2 4 Using the relation to temperature of the Euclidean Wasserstein distance, we conclude. θ J.A. Carrillo, Boletín SEMA (2004).
Kinetic models of Maxwell type. A brief history Part I
. A brief history Part I Department of Mathematics University of Pavia, Italy Porto Ercole, June 8-10 2008 Summer School METHODS AND MODELS OF KINETIC THEORY Outline 1 Introduction Wild result The central
More informationTanaka Theorem for Inelastic Maxwell Models
Tanaka Theorem for Inelastic Maxwell Models François Bolley, José Antonio Carrillo To cite this version: François Bolley, José Antonio Carrillo. Tanaka Theorem for Inelastic Maxwell Models. Communications
More informationAnomalous transport of particles in Plasma physics
Anomalous transport of particles in Plasma physics L. Cesbron a, A. Mellet b,1, K. Trivisa b, a École Normale Supérieure de Cachan Campus de Ker Lann 35170 Bruz rance. b Department of Mathematics, University
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 informationExistence of global solutions to the Cauchy problem for the inelastic Boltzmann equation with near-vacuum data
Existence of global solutions to the Cauchy problem for the inelastic Boltzmann equation with near-vacuum data Ricardo J. Alonso July 12, 28 Abstract The Cauchy problem for the inelastic Boltzmann equation
More informationLecture No 2 Degenerate Diffusion Free boundary problems
Lecture No 2 Degenerate Diffusion Free boundary problems Columbia University IAS summer program June, 2009 Outline We will discuss non-linear parabolic equations of slow diffusion. Our model is the porous
More informationEntropic structure of the Landau equation. Coulomb interaction
with Coulomb interaction Laurent Desvillettes IMJ-PRG, Université Paris Diderot May 15, 2017 Use of the entropy principle for specific equations Spatially Homogeneous Kinetic equations: 1 Fokker-Planck:
More informationEntropy-dissipation methods I: Fokker-Planck equations
1 Entropy-dissipation methods I: Fokker-Planck equations Ansgar Jüngel Vienna University of Technology, Austria www.jungel.at.vu Introduction Boltzmann equation Fokker-Planck equations Degenerate parabolic
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 informationLong time behavior of non autonomous Fokker Planck equations and the cooling of granular gases
Long time behavior of non autonomous Fokker Planck equations and the cooling of granular gases Bertrand Lods Politecnico di Torino, Dipartimento di Matematica, Corso Duca degli Abruzzi, 1129 Torino, Italy.
More informationExponential tail behavior for solutions to the homogeneous Boltzmann equation
Exponential tail behavior for solutions to the homogeneous Boltzmann equation Maja Tasković The University of Texas at Austin Young Researchers Workshop: Kinetic theory with applications in physical sciences
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 informationFluid-Particles Interaction Models Asymptotics, Theory and Numerics I
Fluid-Particles Interaction Models Asymptotics, Theory and Numerics I J. A. Carrillo collaborators: T. Goudon (Lille), P. Lafitte (Lille) and F. Vecil (UAB) (CPDE 2005), (JCP, 2008), (JSC, 2008) ICREA
More informationGlobal Weak Solution to the Boltzmann-Enskog equation
Global Weak Solution to the Boltzmann-Enskog equation Seung-Yeal Ha 1 and Se Eun Noh 2 1) Department of Mathematical Science, Seoul National University, Seoul 151-742, KOREA 2) Department of Mathematical
More informationModelling and numerical methods for the diffusion of impurities in a gas
INERNAIONAL JOURNAL FOR NUMERICAL MEHODS IN FLUIDS Int. J. Numer. Meth. Fluids 6; : 6 [Version: /9/8 v.] Modelling and numerical methods for the diffusion of impurities in a gas E. Ferrari, L. Pareschi
More informationAnalisis para la ecuacion de Boltzmann Soluciones y Approximaciones
Analisis para la ecuacion de Boltzmann Soluciones y Approximaciones Irene M. Gamba Department of Mathematics and ICES The University of Texas at Austin Mar del PLata, Julio 2012 Collaborators: R. Alonso,
More informationin Bounded Domains Ariane Trescases CMLA, ENS Cachan
CMLA, ENS Cachan Joint work with Yan GUO, Chanwoo KIM and Daniela TONON International Conference on Nonlinear Analysis: Boundary Phenomena for Evolutionnary PDE Academia Sinica December 21, 214 Outline
More informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
More informationApplications of the time derivative of the L 2 -Wasserstein distance and the free entropy dissipation
Applications of the time derivative of the L 2 -Wasserstein distance and the free entropy dissipation Hiroaki YOSHIDA Ochanomizu University Tokyo, Japan at Fields Institute 23 July, 2013 Plan of talk 1.
More informationHypocoercivity for kinetic equations with linear relaxation terms
Hypocoercivity for kinetic equations with linear relaxation terms Jean Dolbeault dolbeaul@ceremade.dauphine.fr CEREMADE CNRS & Université Paris-Dauphine http://www.ceremade.dauphine.fr/ dolbeaul (A JOINT
More informationAlignment processes on the sphere
Alignment processes on the sphere Amic Frouvelle CEREMADE Université Paris Dauphine Joint works with : Pierre Degond (Imperial College London) and Gaël Raoul (École Polytechnique) Jian-Guo Liu (Duke University)
More informationNonlinear stabilization via a linear observability
via a linear observability Kaïs Ammari Department of Mathematics University of Monastir Joint work with Fathia Alabau-Boussouira Collocated feedback stabilization Outline 1 Introduction and main result
More informationExponential moments for the homogeneous Kac equation
Exponential moments for the homogeneous Kac equation Maja Tasković University of Pennsylvania Young Researchers Workshop: Stochastic and deterministic methods in kinetic theory, Duke University, November
More informationAsymptotic stability of an evolutionary nonlinear Boltzmann-type equation
Acta Polytechnica Hungarica Vol. 14, No. 5, 217 Asymptotic stability of an evolutionary nonlinear Boltzmann-type equation Roksana Brodnicka, Henryk Gacki Institute of Mathematics, University of Silesia
More informationGlobal Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations
Global Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations C. David Levermore Department of Mathematics and Institute for Physical Science and Technology
More informationSemigroup factorization and relaxation rates of kinetic equations
Semigroup factorization and relaxation rates of kinetic equations Clément Mouhot, University of Cambridge Analysis and Partial Differential Equations seminar University of Sussex 24th of february 2014
More informationQUANTUM MODELS FOR SEMICONDUCTORS AND NONLINEAR DIFFUSION EQUATIONS OF FOURTH ORDER
QUANTUM MODELS FOR SEMICONDUCTORS AND NONLINEAR DIFFUSION EQUATIONS OF FOURTH ORDER MARIA PIA GUALDANI The modern computer and telecommunication industry relies heavily on the use of semiconductor devices.
More informationConvergence to self-similarity for the Boltzmann equation for strongly inelastic Maxwell molecules
Ann. I. H. Poincaré A 7 00 79 737 www.elsevier.com/locate/anihpc Convergence to self-similarity for the Boltzmann equation for strongly inelastic Maxwell molecules G. Furioli a,, A. Pulvirenti b, E. Terraneo
More informationKinetic theory of gases
Kinetic theory of gases Toan T. Nguyen Penn State University http://toannguyen.org http://blog.toannguyen.org Graduate Student seminar, PSU Jan 19th, 2017 Fall 2017, I teach a graduate topics course: same
More informationUne approche hypocoercive L 2 pour l équation de Vlasov-Fokker-Planck
Une approche hypocoercive L 2 pour l équation de Vlasov-Fokker-Planck Jean Dolbeault dolbeaul@ceremade.dauphine.fr CEREMADE CNRS & Université Paris-Dauphine http://www.ceremade.dauphine.fr/ dolbeaul (EN
More informationAnalysis of a Mollified Kinetic Equation for Granular Media. William Thompson B.Sc., University of Victoria, 2014
Analysis of a Mollified Kinetic Equation for Granular Media by William Thompson B.Sc., University of Victoria, 2014 A Thesis Submitted in Partial Fulfillment of the Requirements for the Degree of MASTER
More informationScore functions, generalized relative Fisher information and applications
Score functions, generalized relative Fisher information and applications Giuseppe Toscani January 19, 2016 Abstract Generalizations of the linear score function, a well-known concept in theoretical statistics,
More informationALMOST EXPONENTIAL DECAY NEAR MAXWELLIAN
ALMOST EXPONENTIAL DECAY NEAR MAXWELLIAN ROBERT M STRAIN AND YAN GUO Abstract By direct interpolation of a family of smooth energy estimates for solutions near Maxwellian equilibrium and in a periodic
More informationEXPONENTIAL AND ALGEBRAIC RELAXATION IN KINETIC MODELS FOR WEALTH DISTRIBUTION
1 EXPONENTIAL AND ALGEBAIC ELAXATION IN KINETIC MODELS FO WEALTH DISTIBUTION B. DÜING Institut für Analysis und Scientific Computing, Technische Universität Wien, 1040 Wien, Austria. D. MATTHES and G.
More informationClassical inequalities for the Boltzmann collision operator.
Classical inequalities for the Boltzmann collision operator. Ricardo Alonso The University of Texas at Austin IPAM, April 2009 In collaboration with E. Carneiro and I. Gamba Outline Boltzmann equation:
More informationContractive metrics for scalar conservation laws
Contractive metrics for scalar conservation laws François Bolley 1, Yann Brenier 2, Grégoire Loeper 34 Abstract We consider nondecreasing entropy solutions to 1-d scalar conservation laws and show that
More informationFrom a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction
From a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction Carnegie Mellon University Center for Nonlinear Analysis Working Group, October 2016 Outline 1 Physical Framework 2 3 Free Energy
More informationVALIDITY OF THE BOLTZMANN EQUATION
VALIDITY OF THE BOLTZMANN EQUATION BEYOND HARD SPHERES based on joint work with M. Pulvirenti and C. Saffirio Sergio Simonella Technische Universität München Sergio Simonella - TU München Academia Sinica
More informationEXAMPLE OF A FIRST ORDER DISPLACEMENT CONVEX FUNCTIONAL
EXAMPLE OF A FIRST ORDER DISPLACEMENT CONVEX FUNCTIONAL JOSÉ A. CARRILLO AND DEJAN SLEPČEV Abstract. We present a family of first-order functionals which are displacement convex, that is convex along the
More informationPart 1 Introduction Degenerate Diffusion and Free-boundaries
Part 1 Introduction Degenerate Diffusion and Free-boundaries Columbia University De Giorgi Center - Pisa June 2012 Introduction We will discuss, in these lectures, certain geometric and analytical aspects
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 informationPATH FUNCTIONALS OVER WASSERSTEIN SPACES. Giuseppe Buttazzo. Dipartimento di Matematica Università di Pisa.
PATH FUNCTIONALS OVER WASSERSTEIN SPACES Giuseppe Buttazzo Dipartimento di Matematica Università di Pisa buttazzo@dm.unipi.it http://cvgmt.sns.it ENS Ker-Lann October 21-23, 2004 Several natural structures
More informationRecent result on porous medium equations with nonlocal pressure
Recent result on porous medium equations with nonlocal pressure Diana Stan Basque Center of Applied Mathematics joint work with Félix del Teso and Juan Luis Vázquez November 2016 4 th workshop on Fractional
More informationL 1 stability of conservation laws for a traffic flow model
Electronic Journal of Differential Equations, Vol. 2001(2001), No. 14, pp. 1 18. ISSN: 1072-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu ftp ejde.math.unt.edu (login:
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 informationRelative entropies, suitable weak solutions, and weak-strong uniqueness for the compressible Navier-Stokes system
Relative entropies, suitable weak solutions, and weak-strong uniqueness for the compressible Navier-Stokes system Institute of Mathematics, Academy of Sciences of the Czech Republic, Prague joint work
More informationThe logarithmic Sobolev inequality in W 1,p and related questions
The logarithmic Sobolev inequality in W 1,p and related questions Jean DOLBEAULT Ceremade, Université Paris IX-Dauphine, Place de Lattre de Tassigny, 75775 Paris Cédex 16, France E-mail: dolbeaul@ceremade.dauphine.fr
More informationFrom Hamiltonian particle systems to Kinetic equations
From Hamiltonian particle systems to Kinetic equations Università di Roma, La Sapienza WIAS, Berlin, February 2012 Plan of the lectures 1 Particle systems and BBKGY hierarchy (the paradigm of the Kinetic
More informationOPTIMAL CONVERGENCE RATES FOR THE COMPRESSIBLE NAVIER-STOKES EQUATIONS WITH POTENTIAL FORCES
OPTIMAL CONVERGENCE RATES FOR THE COMPRESSIBLE NAVIER-STOKES EQUATIONS WITH POTENTIAL FORCES RENJUN DUAN Department of Mathematics, City University of Hong Kong 83 Tat Chee Avenue, Kowloon, Hong Kong,
More informationOBSERVABILITY INEQUALITY AND DECAY RATE FOR WAVE EQUATIONS WITH NONLINEAR BOUNDARY CONDITIONS
Electronic Journal of Differential Equations, Vol. 27 (27, No. 6, pp. 2. ISSN: 72-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu OBSERVABILITY INEQUALITY AND DECAY RATE FOR WAVE EQUATIONS
More informationRegularity for the optimal transportation problem with Euclidean distance squared cost on the embedded sphere
Regularity for the optimal transportation problem with Euclidean distance squared cost on the embedded sphere Jun Kitagawa and Micah Warren January 6, 011 Abstract We give a sufficient condition on initial
More informationMATHS 730 FC Lecture Notes March 5, Introduction
1 INTRODUCTION MATHS 730 FC Lecture Notes March 5, 2014 1 Introduction Definition. If A, B are sets and there exists a bijection A B, they have the same cardinality, which we write as A, #A. If there exists
More informationMoment Inequalities and High-Energy Tails for Boltzmann Equations with Inelastic Interactions
Journal of Statistical Physics, Vol. 6, No. 5/6, September 004 ( 004) Moment Inequalities and High-Energy Tails for Boltzmann Equations with Inelastic Interactions A. V. Bobylev, I. M. Gamba, and V. A.
More informationVelocity Distribution in a Viscous Granular Gas
Velocity Distribution in a Viscous Granular Gas Alexandre Rosas,, Daniel ben-avraham 3, and Katja Lindenberg,4 Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA
More informationMathematical modelling of collective behavior
Mathematical modelling of collective behavior Young-Pil Choi Fakultät für Mathematik Technische Universität München This talk is based on joint works with José A. Carrillo, Maxime Hauray, and Samir Salem
More informationNonlinear elliptic systems with exponential nonlinearities
22-Fez conference on Partial Differential Equations, Electronic Journal of Differential Equations, Conference 9, 22, pp 139 147. http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu
More informationA model of alignment interaction for oriented particles with phase transition
A model of alignment interaction for oriented particles with phase transition Amic Frouvelle ACMAC Joint work with Jian-Guo Liu (Duke University, USA) and Pierre Degond (Institut de Mathématiques de Toulouse,
More informationHeat Flows, Geometric and Functional Inequalities
Heat Flows, Geometric and Functional Inequalities M. Ledoux Institut de Mathématiques de Toulouse, France heat flow and semigroup interpolations Duhamel formula (19th century) pde, probability, dynamics
More informationPreconditioning techniques in frame theory and probabilistic frames
Preconditioning techniques in frame theory and probabilistic frames Department of Mathematics & Norbert Wiener Center University of Maryland, College Park AMS Short Course on Finite Frame Theory: A Complete
More informationThe Exponential Formula for the Wasserstein Metric
The Exponential Formula for the Wasserstein Metric A dissertation submitted to Rutgers, The State University of New Jersey, in partial fulfillment of the requirements for the degree of Doctor of Philosophy
More informationDifferent types of phase transitions for a simple model of alignment of oriented particles
Different types of phase transitions for a simple model of alignment of oriented particles Amic Frouvelle CEREMADE Université Paris Dauphine Joint work with Jian-Guo Liu (Duke University, USA) and Pierre
More informationStability in geometric & functional inequalities
Stability in geometric & functional inequalities Alessio Figalli Abstract. The aim of this note is to review recent stability results for some geometric and functional inequalities, and to describe applications
More informationA Selective Modern History of the Boltzmann and Related Eq
A Selective Modern History of the Boltzmann and Related Equations October 2014, Fields Institute Synopsis 1. I have an ambivalent relation to surveys! 2. Key Words, Tools, People 3. Powerful Tools, I:
More informationLecture 21: Convergence of transformations and generating a random variable
Lecture 21: Convergence of transformations and generating a random variable If Z n converges to Z in some sense, we often need to check whether h(z n ) converges to h(z ) in the same sense. Continuous
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 informationA model of alignment interaction for oriented particles with phase transition
A model of alignment interaction for oriented particles with phase transition Amic Frouvelle Archimedes Center for Modeling, Analysis & Computation (ACMAC) University of Crete, Heraklion, Crete, Greece
More informationOn the Interior Boundary-Value Problem for the Stationary Povzner Equation with Hard and Soft Interactions
On the Interior Boundary-Value Problem for the Stationary Povzner Equation with Hard and Soft Interactions Vladislav A. Panferov Department of Mathematics, Chalmers University of Technology and Göteborg
More informationMATH MEASURE THEORY AND FOURIER ANALYSIS. Contents
MATH 3969 - MEASURE THEORY AND FOURIER ANALYSIS ANDREW TULLOCH Contents 1. Measure Theory 2 1.1. Properties of Measures 3 1.2. Constructing σ-algebras and measures 3 1.3. Properties of the Lebesgue measure
More information1/12/05: sec 3.1 and my article: How good is the Lebesgue measure?, Math. Intelligencer 11(2) (1989),
Real Analysis 2, Math 651, Spring 2005 April 26, 2005 1 Real Analysis 2, Math 651, Spring 2005 Krzysztof Chris Ciesielski 1/12/05: sec 3.1 and my article: How good is the Lebesgue measure?, Math. Intelligencer
More informationProbability and Measure
Part II Year 2018 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2018 84 Paper 4, Section II 26J Let (X, A) be a measurable space. Let T : X X be a measurable map, and µ a probability
More informationComplex systems: Self-organization vs chaos assumption
1 Complex systems: Self-organization vs chaos assumption P. Degond Institut de Mathématiques de Toulouse CNRS and Université Paul Sabatier pierre.degond@math.univ-toulouse.fr (see http://sites.google.com/site/degond/)
More informationBoltzmann Equations For Mixtures of Maxwell Gases: Exact Solutions and Power Like Tails
Journal of Statistical Physics ( C 6 ) DOI: 1.17/s1955-6-944-8 Boltzmann Equations For Mixtures of Maxwell Gases: Exact Solutions and Power Like Tails A. V. Bobylev 1 and I. M. Gamba Received July 1, 5;
More informationStability results for Logarithmic Sobolev inequality
Stability results for Logarithmic Sobolev inequality Daesung Kim (joint work with Emanuel Indrei) Department of Mathematics Purdue University September 20, 2017 Daesung Kim (Purdue) Stability for LSI Probability
More informationCOMPLETE METRIC SPACES AND THE CONTRACTION MAPPING THEOREM
COMPLETE METRIC SPACES AND THE CONTRACTION MAPPING THEOREM A metric space (M, d) is a set M with a metric d(x, y), x, y M that has the properties d(x, y) = d(y, x), x, y M d(x, y) d(x, z) + d(z, y), x,
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 informationFractal Conservation Laws: Global Smooth Solutions and Vanishing Regularization
Progress in Nonlinear Differential Equations and Their Applications, Vol. 63, 217 224 c 2005 Birkhäuser Verlag Basel/Switzerland Fractal Conservation Laws: Global Smooth Solutions and Vanishing Regularization
More informationNodal solutions of a NLS equation concentrating on lower dimensional spheres
Presentation Nodal solutions of a NLS equation concentrating on lower dimensional spheres Marcos T.O. Pimenta and Giovany J. M. Figueiredo Unesp / UFPA Brussels, September 7th, 2015 * Supported by FAPESP
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationThe inviscid limit to a contact discontinuity for the compressible Navier-Stokes-Fourier system using the relative entropy method
The inviscid limit to a contact discontinuity for the compressible Navier-Stokes-Fourier system using the relative entropy method Alexis Vasseur, and Yi Wang Department of Mathematics, University of Texas
More informationFirst order differential equations
First order differential equations Samy Tindel Purdue University Differential equations and linear algebra - MA 262 Taken from Differential equations and linear algebra by Goode and Annin Samy T. First
More information3 (Due ). Let A X consist of points (x, y) such that either x or y is a rational number. Is A measurable? What is its Lebesgue measure?
MA 645-4A (Real Analysis), Dr. Chernov Homework assignment 1 (Due ). Show that the open disk x 2 + y 2 < 1 is a countable union of planar elementary sets. Show that the closed disk x 2 + y 2 1 is a countable
More informationdegenerate Fokker-Planck equations with linear drift
TECHNISCHE UNIVERSITÄT WIEN Entropy method for hypocoercive & non-symmetric Fokker-Planck equations with linear drift Anton ARNOLD with Jan Erb; Franz Achleitner Paris, April 2017 Anton ARNOLD (TU Vienna)
More informationSMOOTH OPTIMAL TRANSPORTATION ON HYPERBOLIC SPACE
SMOOTH OPTIMAL TRANSPORTATION ON HYPERBOLIC SPACE JIAYONG LI A thesis completed in partial fulfilment of the requirement of Master of Science in Mathematics at University of Toronto. Copyright 2009 by
More informationNonlinear Diffusion. The Porous Medium Equation. From Analysis to Physics and Geometry
Nonlinear Diffusion. The Porous Medium Equation. From Analysis to Physics and Geometry Juan Luis Vázquez Departamento de Matemáticas Universidad Autónoma de Madrid http://www.uam.es/personal pdi/ciencias/jvazquez/
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 informationA TWO PARAMETERS AMBROSETTI PRODI PROBLEM*
PORTUGALIAE MATHEMATICA Vol. 53 Fasc. 3 1996 A TWO PARAMETERS AMBROSETTI PRODI PROBLEM* C. De Coster** and P. Habets 1 Introduction The study of the Ambrosetti Prodi problem has started with the paper
More informationHeat equation and the sharp Young s inequality
Noname manuscript No. will be inserted by the editor) Heat equation and the sharp Young s inequality Giuseppe Toscani the date of receipt and acceptance should be inserted later Abstract We show that the
More informationA Relaxed Explicit Extragradient-Like Method for Solving Generalized Mixed Equilibria, Variational Inequalities and Constrained Convex Minimization
, March 16-18, 2016, Hong Kong A Relaxed Explicit Extragradient-Like Method for Solving Generalized Mixed Equilibria, Variational Inequalities and Constrained Convex Minimization Yung-Yih Lur, Lu-Chuan
More informationStability for a GNS inequality and the Log-HLS inequality, with application to the critical mass Keller-Segel equation
Stability for a GNS inequality and the Log-HLS inequality, with application to the critical mass Keller-Segel equation Eric A. Carlen and Alessio Figalli 2. Department of Mathematics, Hill Center, Rutgers
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 informationNonlinear Diffusion. Porous Medium and Fast Diffusion. From Analysis to Physics and Geometry
Juan L. Vázquez - Nonlinear Diffusion. Porous Medium and Fast Diffusion Equations p. 1/77 Nonlinear Diffusion. Porous Medium and Fast Diffusion. From Analysis to Physics and Geometry Juan Luis Vázquez
More informationOPTIMAL TRANSPORTATION PLANS AND CONVERGENCE IN DISTRIBUTION
OPTIMAL TRANSPORTATION PLANS AND CONVERGENCE IN DISTRIBUTION J.A. Cuesta-Albertos 1, C. Matrán 2 and A. Tuero-Díaz 1 1 Departamento de Matemáticas, Estadística y Computación. Universidad de Cantabria.
More informationAbout the method of characteristics
About the method of characteristics Francis Nier, IRMAR, Univ. Rennes 1 and INRIA project-team MICMAC. Joint work with Z. Ammari about bosonic mean-field dynamics June 3, 2013 Outline The problem An example
More informationLARGE TIME BEHAVIOR OF THE RELATIVISTIC VLASOV MAXWELL SYSTEM IN LOW SPACE DIMENSION
Differential Integral Equations Volume 3, Numbers 1- (1), 61 77 LARGE TIME BEHAVIOR OF THE RELATIVISTIC VLASOV MAXWELL SYSTEM IN LOW SPACE DIMENSION Robert Glassey Department of Mathematics, Indiana University
More informationCopyright by Ricardo Jose Alonso 2008
Copyright by Ricardo Jose Alonso 2008 The Dissertation Committee for Ricardo Jose Alonso certifies that this is the approved version of the following dissertation: The Boltzmann equation: Sharp Povzner
More informationHOMEOMORPHISMS OF BOUNDED VARIATION
HOMEOMORPHISMS OF BOUNDED VARIATION STANISLAV HENCL, PEKKA KOSKELA AND JANI ONNINEN Abstract. We show that the inverse of a planar homeomorphism of bounded variation is also of bounded variation. In higher
More informationarxiv: v2 [math.fa] 17 May 2016
ESTIMATES ON SINGULAR VALUES OF FUNCTIONS OF PERTURBED OPERATORS arxiv:1605.03931v2 [math.fa] 17 May 2016 QINBO LIU DEPARTMENT OF MATHEMATICS, MICHIGAN STATE UNIVERSITY EAST LANSING, MI 48824, USA Abstract.
More informationAN ASYMPTOTIC MEAN VALUE CHARACTERIZATION FOR p-harmonic FUNCTIONS. To the memory of our friend and colleague Fuensanta Andreu
AN ASYMPTOTIC MEAN VALUE CHARACTERIZATION FOR p-harmonic FUNCTIONS JUAN J. MANFREDI, MIKKO PARVIAINEN, AND JULIO D. ROSSI Abstract. We characterize p-harmonic functions in terms of an asymptotic mean value
More informationDifferent types of phase transitions for a simple model of alignment of oriented particles
Different types of phase transitions for a simple model of alignment of oriented particles Amic Frouvelle Université Paris Dauphine Joint work with Jian-Guo Liu (Duke University, USA) and Pierre Degond
More informationStatistical mechanics of random billiard systems
Statistical mechanics of random billiard systems Renato Feres Washington University, St. Louis Banff, August 2014 1 / 39 Acknowledgements Collaborators: Timothy Chumley, U. of Iowa Scott Cook, Swarthmore
More information