Mathematical modelling of collective behavior
|
|
- Randall Fowler
- 5 years ago
- Views:
Transcription
1 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 Mathematical Topics in Kinetic Theory University of Cambridge May 13, / 29
2 Introduction Singular communication weights Normalized communication weights Sharp sensitivity regions Swarming in Nature 2 / 29
3 Individual Based Models Figure: 3 Zones 3 / 29
4 Cucker-Smale model Let (x i, v i ) be the position and velocity of i-th particle. Then Cucker-Smale model reads as dx i dt = v i, t > 0, i = 1,, N, dv i dt = 1 N 1 ψ(x j x i )(v j v i ), ψ(x) := N (1 + x 2, β 0. ) β/2 j=1 Definition: Let P := {(x i, v i )} N i=1 be an N-body interacting system. Then P exhibits a flocking if and only if the following two relations hold: lim v i (t) v j (t) = 0 and sup x i (t) x j (t) <, 1 i, j N. t 0 t< Property: Conservation of momentum. Flocking estimate: If ψ(r) dr = ( β 1), then we have sup max x i (t) x j (t) < and max v i (t) v j (t) Ce Ct for t 0. 0 t< 1 i,j N 1 i,j N Ref.- Cucker-Smale(2007), Ha-Tadmor(2008), Ha-Liu(2009), Carrillo-Fornasier-Rosado-Toscani(2010),... 4 / 29
5 Numerical Simulation by Sergio Pérez 5 / 29
6 Kinetic Cucker-Smale equation When N, the macroscopic observables for Cucker-Smale model can be calculated from the velocity moments of the density function f = (x, v, t) which is a solution to the following Vlasov-type equation: Flocking: tf + v x f + v [F (f )f ] = 0, (x, v) R d R d, F (f )(x, v, t) := ψ(x y)(w v)f (y, w, t) dydw. R d R d v v c 2 f dxdv 0 as t where v c = v f dxdv f dxdv. Several drawbacks: - Collision between individuals - Dynamics away from equilibrium - Vision geometrical constraints -... Ref.- Cucker-Dong(2010), Park-Kim-Ha(2010), Motsch-Tadmor(2011), Ahn-Choi-Ha-Lee(2012), Agueh-Illner-Richardson(2011),... 6 / 29
7 Collision avoidance between individuals 7 / 29
8 Previous work Cucker-Smale model with collision avoiding forces: dx i dt = v i, t > 0, i = 1,, N, dv i dt = 1 N ψ(x i x j )(v j v i ) 1 N 1 K(x i x j ), ψ(x) = N N (1 + x 2 ) β/2. j i j i }{{}}{{} Alignment Attraction/Repulsion Ref.- Cucker-Dong(2010), Ha-Park-Kim(2010), Yang-Chen(2014),... 8 / 29
9 Numerical simulation K(x) = kφ(x) α x 2 2 where φ = δ 0. 9 / 29
10 Numerical simulation K(x) = kφ(x) α x 2 2 where φ = δ 0. by Sergio Pérez 9 / 29
11 Singular communication weights Cucker-Smale model with singular communication weights: dx i dt = v i, t > 0, i = 1,, N, dv i dt = 1 N ψ(x j x i )(v j v i ), ψ(x) = ψ s (x) := 1 N x α, α > 0. j=1 (1) 10 / 29
12 Singular communication weights Cucker-Smale model with singular communication weights: dx i dt = v i, t > 0, i = 1,, N, dv i dt = 1 N ψ(x j x i )(v j v i ), ψ(x) = ψ s (x) := 1 N x α, α > 0. j=1 (1) Flocking estimate: Let α 1. Suppose that the initial data satisfies v 0 v c(0) l < 1 ψ(s) ds. 2 2 x 0 x c (0) l Then there exist positive constants c 1 > c 0 > 0 such that x(t) x c(0) l [c 0, c 1 ] and v(t) v c(0) l v 0 v c(0) l e ψ(2c1)t. Ref.- Ahn-Choi-Ha-Lee(2012) 10 / 29
13 Sharp condition to avoid collision Finite-time collision estimate(α < 1): Let d = 1. If we consider the two-particle system, then the following statements are equivalent: 1. There exists a time t 0 < such that φ(t 0 ) = φ(t 0 ) = Initial data satisfy φ(0) = 2Ψ(φ(0)), where Ψ(s) := 1 1 α s1 α is a primitive of φ and φ(t) = x 1 (t) x 2 (t). Non collision estimate(α 1): Suppose that the initial data (x 0, v 0 ) are non-collisional, i.e., they satisfy x i0 x j0 for 1 i j N. Then the system (1) admits a unique smooth solution. Moreover, the trajectories of this solution are also non-collisional, i.e., x i (t) x j (t) for 1 i j N, t 0. Ref.- Peszek(2014), Carrillo-C.-Mucha-Peszek(work in progress) 11 / 29
14 Cucker-Smale model with singular communication weights Main equation: tf + v x f + v (F (f )f ) = 0, (x, v) R d R d, t > 0, (2) where F denotes the alignment force between particles: F (f )(x, v, t) = ψ(x y)(w v)f (y, w) dydw, ψ(x) = 1 with α (0, d 1). x α Definition: For a given T (0, ), f is a weak solution of (3) on the time-interval [0, T ) if and only if the following condition are satisfied: (i) f L (0, T ; (L 1 + Lp )(R d R d )) C([0, T ]; P 1 (R d R d )), (ii) f satisfies the equation (3) in the weak sense. 12 / 29
15 Local existence and uniqueness of weak solutions Theorem(Carrillo-C.-Hauray, 2014): Suppose that f 0 is compactly supported in velocity and (α + 1)p < d, and the initial data f 0 satisfies f 0 (L 1 + L p )(R d R d ) P 1 (R d R d ). Then there exist T > 0 and unique weak solution f to the system (3) in the sense of Definition 1 on the time interval [0, T ]. Furthermore, if f i, i = 1, 2 are two such solutions to (3), then we have the following d 1 -stability estimate. d dt d 1(f 1 (t), f 2 (t)) Cd 1 (f 1 (t), f 2 (t)), for t [0, T ], where C is a positive constant. 13 / 29
16 Idea of Proof (Regularization) The approximate solutions are obtained by solving: tf ε + v x f ε + v [F ε ] (f ε)f ε = 0, (x, v) R d R d, t > 0, F ε (f ε)(x, v, t) := ψ ε(x y)(w v)f ε(y, w, t)dydw, R d R d f ε(x, v, 0) =: f 0 (x, v). where ψ ε := ψ θ ε. (Support estimate of f ε in velocity) Consider the forward bi-characteristics (X ε(s; 0, x, v), V ε(s; 0, x, v)) satisfying the following ODE system: dx ε(s) = V ε(s), ds dv ε(s) = ψ ε (X ε(s) y) (w V ε(s))f ε(y, w, s)dydw. ds R d R d Set Ω ε(t) and R v ε (t) the v-projection of compact suppf ε(, t) and maximum value of v in Ω ε(t), respectively: Ω ε(t) := {v R d : (x, v) R d R d such that f ε(x, v, t) 0}, Rε v (t) := max v. v Ω ε(t) Then we have R v ε (t) R v ε (0) = R v 0 := max v Ω(0) v. 14 / 29
17 Idea of Proof(Conti.) (Uniform bound estimate): Let f ε be the solution to regularized equation. Then there exists a T > 0 such that the uniform L 1 L p -estimate of f ε sup f ε L 1 L p C, t [0,T ] holds, where C is a positive constant independent of ε. (f ε is a Cauchy sequence): Let f ε and f ε be two solutions of the regularized equation. Then there exists C independent of ε and ε such that holds for all t 0. d dt d 1(f ε(t), f ε (t)) C(d 1 (f ε(t), f ε (t)) + ε + ε ), (Passing to the limit ε 0): The limit curve of measure f obtained from the above is a solution of the equation (3). Open question: Rigorous derivation of the mean-field limit 15 / 29
18 Further extensions: Flocking models with singular kernels Main equation: tf + v x f + v (F (f )f ) = 0, (x, v) R d R d, t > 0, (3) where F denotes the alignment force between particles: F (f )(x, v, t) = ψ(x y) v φ(v w)f (y, w, t)dydw. R d R d Here, the potential function φ(v) for the velocity coupling is given by Two different cases: φ(v) = 1 β v β, where β > 0. Singular communication weight and super-linear velocity coupling: ψ(x) = ψ 1 (x) := 1 x α with α (0, d 1) and β 2. Regular communication weight and sub-linear velocity coupling: ( ) 3 d ψ(x) = ψ 2 (x) 0 symmetric, with ψ 2 (L loc Lip loc ) (Rd ) and β, / 29
19 Normalized communication weights Cucker-Smale model: dx i dt = v i, t > 0, i = 1,, N, dv i dt = 1 N 1 ψ(x j x i )(v j v i ), ψ(x) = N (1 + x 2 ) β/2. j=1 Cucker-Smale model with normalized communication weights: dx i dt = v i, t > 0, i = 1,, N, dv i dt = 1 N k=1 ψ(x i x k ) N ψ(x j x i )(v j v i ). j=1 Motsch-Tadmor(2011) ψ 2 (r) dr = = Flocking Ref.- Tadmor-Tan (2014): ψ(r) dr = = Flocking 17 / 29
20 Numerical simulation by Sergio Pérez 18 / 29
21 Interacting diffusing particle systems and kinetic equation Interacting diffusing particle system: dx i = v i dt, t > 0, i = 1,, N, dv i = 1 N k=1 ψ(x i x k ) N ψ(x j x i )(v j v i )dt + 2dB i. j=1 19 / 29
22 Interacting diffusing particle systems and kinetic equation Interacting diffusing particle system: Kinetic model: dx i = v i dt, t > 0, i = 1,, N, dv i = where ũ F is given by with 1 N k=1 ψ(x i x k ) N ψ(x j x i )(v j v i )dt + 2dB i. j=1 tf + v x F + v ((ũ F v)f ) = v F, x R d, v R d, ũ F (x, t) := (ψ bf )(x, t) (ψ a F )(x, t), a F (x, t) := F (x, v, t) dv and b F (x, t) := vf (x, v, t) dv. R d R d 19 / 29
23 Interacting diffusing particle systems and kinetic equation Interacting diffusing particle system: Kinetic model: dx i = v i dt, t > 0, i = 1,, N, dv i = where ũ F is given by with 1 N k=1 ψ(x i x k ) N ψ(x j x i )(v j v i )dt + 2dB i. j=1 tf + v x F + v ((ũ F v)f ) = v F, x R d, v R d, ũ F (x, t) := (ψ bf )(x, t) (ψ a F )(x, t), a F (x, t) := F (x, v, t) dv and b F (x, t) := vf (x, v, t) dv. R d R d Consider the singular limit in which the communication weight ψ converges to a Dirac distribution. ũ F u F := bf a F = vf dv Rd R d F dv. Ref.- Karper-Mellet-Trivisa(2014) 19 / 29
24 Perturbation framework Nonlinear Fokker-Planck equation: Define the perturbation f = f (x, v, t) by tf + v x F = v ( v F + (v u F )F ). F = M + ) 1 Mf, M = M(v) = exp ( v 2 (2π) d/2 2 The equation for the perturbation f satisfies tf + v x f + u F v f = Lf + Γ(f, f ), (4) where the linear part Lf and the nonlinear part Γ(f, f ) are given by Lf := 1 ( ( )) ( ) f 1 v M v and Γ(f, f ) = u F M M 2 vf + v M, respectively. Γ(f, f ) is not bilinear! Ref.- Duan-Fornasier-Toscani(2010), Karper-Mellet-Trivisa(2013) 20 / 29
25 Global classical solutions and time-asymptotic behavior Theorem(C. 2015): Let d 3 and s 2[d/2] + 2. Suppose F 0 M + Mf 0 0 and f 0 H s ɛ 0 1. Then we have the global existence of the unique classical solution f (x, v, t) to the equation (4) satisfying f C([0, ); H s (R d R d )), F M + Mf 0, and t t f (t) 2 H s + c 1 x (a, b) 2 H 0 s 1 ds + c 1 k x l v {I P}f 2 µds c 2 f 0 2 H s, 0 0 k+l s for some positive constants c 1, c 2 > 0. Furthermore, if f 0 L 2 v (L 1 ) is bounded, we have ( ) f (t) H s C f 0 H s + f 0 L 2 v (L 1 ) (1 + t) d 4, t 0, where C is a positive constant independent of t. Idea of proof: Careful analysis of the nonlinear dissipation term + Coercivity of the linear operator L + Hypocoercivity for a modified linearized Cauchy problem with a non-homogeneous microscopic source. 21 / 29
26 Collective behavior models with sharp sensitivity regions 22 / 29
27 Collective behavior models at particle and kinetic levels Particle system: dx i dt = v i, t > 0, i = 1,, N, dv i dt = m j 1 K(vi )(x i x j )(v j v i ), j i where 1 K(vi ) is the indicator function on the set K(v i ). Here {x i } N i=1, {v i } N i=1, {m i } N i=1 are the position, velocity, and weight of i-th particle, respectively. Kinetic equation: As a mean-field limit, we consider tf + v x f + v [F (f )f ] = 0, (x, v) R d R d, t > 0, (5) where F (f ) is a alignment force of velocities given by F (f )(x, v, t) = 1 K(v) (x y)(w v)f (y, w) dydw. R d R d 23 / 29
28 Main assumptions (H1): K(v) is globally compact, i.e., K(v) is compact and there exists a compact set K such that K(v) K for all v R d. (H2): There exists a family of closed sets v Θ(v) and a constant C such that: - K(v) Θ(v), for all v R d, - Θ(v) ε,+ Cε, for all ε (0, 1), - K(v) K(w) Θ(v) C v w,+, for v, w R d, - Θ(w) Θ(v) C v w,+, for v, w R d. Figure: A sketch of the sets K, K ε,+, and K ε, 24 / 29
29 Examples A ball with a fixed radius: K(v) B(0, r) where B(0, r) := {x R d : x r} with r > 0. A ball with a radius depending on the speed: K(v) = B(0, r( v )) where r : R + R + is bounded and Lipschitz. Θ(v) = K(v) 25 / 29
30 Examples A vision cone: K(v) = C(r, v, θ( v )) := { ( ) } x v x : x r and θ( v ) cos 1 θ( v ), x v with 0 < θ(z) C (R) satisfying θ(z) = π for 0 z 1, θ(z) is decreasing for z 1, and θ(z) θ > 0 as z +. { C(r, v, θ( v )) R(v) if v (1/2, 1), Θ(v) := C(r, v, θ( v )) else, where R(v) = [a(v), b(v)] with a(v) = r v v, b(v) = 2r( v 1) v v. 26 / 29
31 Global weak solutions & stability Definition: For a given T (0, ), f is a weak solution of (5) on the time-interval [0, T ) if and only if the following condition are satisfied: (i) f L (0, T ; (L 1 + L )(R d R d )) C([0, T ]; P 1 (R d R d )), (ii) f satisfies the equation (5) in the weak sense. Theorem(Carrillo-C.-Hauray-Salem 2016): Given an initial data satisfying f 0 (L 1 + L )(R d R d ) P 1 (R d R d ), (6) and assume further that f 0 is compactly supported in velocity. Then there exists a positive time T > 0 such that the system (5) with the sensitivity region set-valued function K(v) satisfying (H1)-(H2) admits a unique weak solution f L (0, T : (L 1 + L )(R d R d )), which is also compactly supported in velocity. Furthermore if f i, i = 1, 2 are two such solutions to the system (5) with initial data f i (0) satisfying (6), we have d 1 (f 1 (t), f 2 (t)) d 1 (f 1 (0), f 2 (0))e Ct, for t [0, T ], where C is a positive constant that depends only on the L (R d R d (0, T )) norm of f 1. Weak-strong stability! 27 / 29
32 Further extensions Second-order model: tf + v x f + v (F (f )f ) = 0, (x, v) R d R d, t > 0, where the force term F (f ) can be chosen from two different types: ψ(x y)1 K(v) (x y)h(w v)f (y, w) dydw F (f )(x, v, t) := x ϕ(x y)1 K(v) (x y)f (y, w) dydw (Cucker-Smale) (Attractive-Repulsive). Here ψ, h, and ϕ denote the communication weight, velocity coupling, and interaction potential, respectively. First-order model: tρ + x (ρu) = 0, x R d, t > 0, where the vector field u is given by u(x, t) := 1 K(w(x)) (x y) x ϕ(x y)ρ(y) dy, R d and ϕ W 1, (R d ) and w is a given orientational field satisfying w W 1, (R d ) and w w 0 > / 29
33 Thank you for your attention 29 / 29
Kinetic swarming models and hydrodynamic limits
Kinetic swarming models and hydrodynamic limits Changhui Tan Department of Mathematics Rice University Young Researchers Workshop: Current trends in kinetic theory CSCAMM, University of Maryland October
More informationFractional parabolic models arising in flocking. dynamics and fluids. Roman Shvydkoy jointly with Eitan Tadmor. April 18, 2017
Fractional April 18, 2017 Cucker-Smale model (2007) We study macroscopic versions of systems modeling self-organized collective dynamics of agents : ẋ i = v i, v i = λ N φ( x i x j )(v j v i ), N j=1 (x
More informationModels of collective displacements: from microscopic to macroscopic description
Models of collective displacements: from microscopic to macroscopic description Sébastien Motsch CSCAMM, University of Maryland joint work with : P. Degond, L. Navoret (IMT, Toulouse) SIAM Analysis of
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 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 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 informationFirst-order aggregation models with alignment
First-order aggregation models with alignment R. C. Fetecau W. Sun Changhui Tan March 1, 216 Abstract We include alignment interactions in a well-studied first-order attractive-repulsive macroscopic model
More informationInvisible Control of Self-Organizing Agents Leaving Unknown Environments
Invisible Control of Self-Organizing Agents Leaving Unknown Environments (joint work with G. Albi, E. Cristiani, and D. Kalise) Mattia Bongini Technische Universität München, Department of Mathematics,
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 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 informationRegularity and compactness for the DiPerna Lions flow
Regularity and compactness for the DiPerna Lions flow Gianluca Crippa 1 and Camillo De Lellis 2 1 Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy. g.crippa@sns.it 2 Institut für Mathematik,
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 informationFrom Boltzmann Equations to Gas Dynamics: From DiPerna-Lions to Leray
From Boltzmann Equations to Gas Dynamics: From DiPerna-Lions to Leray C. David Levermore Department of Mathematics and Institute for Physical Science and Technology University of Maryland, College Park
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 informationfor all subintervals I J. If the same is true for the dyadic subintervals I D J only, we will write ϕ BMO d (J). In fact, the following is true
3 ohn Nirenberg inequality, Part I A function ϕ L () belongs to the space BMO() if sup ϕ(s) ϕ I I I < for all subintervals I If the same is true for the dyadic subintervals I D only, we will write ϕ BMO
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More 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 informationDecay rates for partially dissipative hyperbolic systems
Outline Decay rates for partially dissipative hyperbolic systems Basque Center for Applied Mathematics Bilbao, Basque Country, Spain zuazua@bcamath.org http://www.bcamath.org/zuazua/ Numerical Methods
More informationPARABOLIC LIMIT AND STABILITY OF THE VLASOV FOKKER PLANCK SYSTEM
Mathematical Models and Methods in Applied Sciences Vol. 10, No. 7 (2000 1027 1045 c World Scientific Publishing Company PARABOLIC LIMIT AND STABILITY OF THE VLASOV FOKKER PLANCK SYSTEM F. POUPAUD Laboratoire
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 informationanalysis for transport equations and applications
Multi-scale analysis for transport equations and applications Mihaï BOSTAN, Aurélie FINOT University of Aix-Marseille, FRANCE mihai.bostan@univ-amu.fr Numerical methods for kinetic equations Strasbourg
More informationGLOBAL CLASSICAL SOLUTIONS OF THE VLASOV-FOKKER-PLANCK EQUATION WITH LOCAL ALIGNMENT FORCES
GLOBAL CLASSICAL SOLUTIONS OF THE VLASOV-FOKKER-PLANCK EQUATION WITH LOCAL ALIGNMENT FORCES YOUNG-PIL CHOI arxiv:164.86v1 [math.ap] 7 Apr 16 Abstract. In this paper, we are concerne with the global well-poseness
More informationOn the Boltzmann equation: global solutions in one spatial dimension
On the Boltzmann equation: global solutions in one spatial dimension Department of Mathematics & Statistics Colloque de mathématiques de Montréal Centre de Recherches Mathématiques November 11, 2005 Collaborators
More informationWeak-Strong Uniqueness of the Navier-Stokes-Smoluchowski System
Weak-Strong Uniqueness of the Navier-Stokes-Smoluchowski System Joshua Ballew University of Maryland College Park Applied PDE RIT March 4, 2013 Outline Description of the Model Relative Entropy Weakly
More informationThe Dirichlet boundary problems for second order parabolic operators satisfying a Carleson condition
The Dirichlet boundary problems for second order parabolic operators satisfying a Carleson condition Sukjung Hwang CMAC, Yonsei University Collaboration with M. Dindos and M. Mitrea The 1st Meeting of
More informationA LOCAL SENSITIVITY ANALYSIS FOR THE KINETIC CUCKER-SMALE EQUATION WITH RANDOM INPUTS
A LOCAL SENSITIVITY ANALYSIS FOR THE KINETIC CUCKER-SMALE EQUATION WITH RANDOM INPUTS SEUNG-YEAL HA, SHI JIN, AND JINWOOK JUNG Abstract. We present a local sensitivity analysis for the kinetic Cucker-Smale
More informationTRANSPORT IN POROUS MEDIA
1 TRANSPORT IN POROUS MEDIA G. ALLAIRE CMAP, Ecole Polytechnique 1. Introduction 2. Main result in an unbounded domain 3. Asymptotic expansions with drift 4. Two-scale convergence with drift 5. The case
More informationSUMMER SCHOOL - KINETIC EQUATIONS LECTURE II: CONFINED CLASSICAL TRANSPORT Shanghai, 2011.
SUMMER SCHOOL - KINETIC EQUATIONS LECTURE II: CONFINED CLASSICAL TRANSPORT Shanghai, 2011. C. Ringhofer ringhofer@asu.edu, math.la.asu.edu/ chris OVERVIEW Quantum and classical description of many particle
More informationGlobal weak solutions for a Vlasov-Fokker-Planck/Navier-Stokes system of equations
Global weak solutions for a Vlasov-Fokker-Planck/Navier-Stokes system of equations A.Mellet and A.Vasseur Department of Mathematics University of Texas at Austin Abstract We establish the existence of
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 informationTOPOLOGICAL MODELS FOR EMERGENT DYNAMICS WITH SHORT-RANGE INTERACTIONS
TOPOLOGICAL MODELS FOR EMERGENT DYNAMICS WITH SHORT-RANGE INTERACTIONS ROMAN SHVYDKOY AND EITAN TADMOR Abstract. We introduce a new class of models for emergent dynamics. It is based on a new communication
More information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More informationResearch Statement Changhui Tan
Research Statement Changhui Tan (ctan@rice.edu) My research interests lie in the theoretical and numerical analysis of nonlinear partial differential equations. In particular, my current research focuses
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 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 informationMATH 425, FINAL EXAM SOLUTIONS
MATH 425, FINAL EXAM SOLUTIONS Each exercise is worth 50 points. Exercise. a The operator L is defined on smooth functions of (x, y by: Is the operator L linear? Prove your answer. L (u := arctan(xy u
More informationYAN GUO, JUHI JANG, AND NING JIANG
LOCAL HILBERT EXPANSION FOR THE BOLTZMANN EQUATION YAN GUO, JUHI JANG, AND NING JIANG Abstract. We revisit the classical ork of Caflisch [C] for compressible Euler limit of the Boltzmann equation. By using
More informationAn Introduction to Variational Inequalities
An Introduction to Variational Inequalities Stefan Rosenberger Supervisor: Prof. DI. Dr. techn. Karl Kunisch Institut für Mathematik und wissenschaftliches Rechnen Universität Graz January 26, 2012 Stefan
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 informationConservation law equations : problem set
Conservation law equations : problem set Luis Silvestre For Isaac Neal and Elia Portnoy in the 2018 summer bootcamp 1 Method of characteristics For the problems in this section, assume that the solutions
More informationElliptic Operators with Unbounded Coefficients
Elliptic Operators with Unbounded Coefficients Federica Gregorio Universitá degli Studi di Salerno 8th June 2018 joint work with S.E. Boutiah, A. Rhandi, C. Tacelli Motivation Consider the Stochastic Differential
More informationGLOBAL REGULARITY OF TWO-DIMENSIONAL FLOCKING HYDRODYNAMICS
GLOBAL REGULARITY OF TWO-DIMENSIONAL FLOCKING HYDRODYNAMICS SIMING HE AND EITAN TADMOR Abstract. We study the systems of Euler equations which arise from agent-based dynamics driven by velocity alignment.
More informationGlobal regularity for the fractional Euler alignment system
Global regularity for the fractional Euler alignment system am Do Alexander Kiselev Lenya Ryzhik Changhui an January 23, 2017 Abstract We study a pressureless Euler system with a non-linear density-dependent
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 informationKinetic relaxation models for reacting gas mixtures
Kinetic relaxation models for reacting gas mixtures M. Groppi Department of Mathematics and Computer Science University of Parma - ITALY Main collaborators: Giampiero Spiga, Giuseppe Stracquadanio, Univ.
More informationDynamical systems with Gaussian and Levy noise: analytical and stochastic approaches
Dynamical systems with Gaussian and Levy noise: analytical and stochastic approaches Noise is often considered as some disturbing component of the system. In particular physical situations, noise becomes
More informationTraveling waves of a kinetic transport model for the KPP-Fisher equation
Traveling waves of a kinetic transport model for the KPP-Fisher equation Christian Schmeiser Universität Wien and RICAM homepage.univie.ac.at/christian.schmeiser/ Joint work with C. Cuesta (Bilbao), S.
More informationInput to state Stability
Input to state Stability Mini course, Universität Stuttgart, November 2004 Lars Grüne, Mathematisches Institut, Universität Bayreuth Part IV: Applications ISS Consider with solutions ϕ(t, x, w) ẋ(t) =
More informationarxiv: v1 [math.ap] 28 Apr 2009
ACOUSTIC LIMIT OF THE BOLTZMANN EQUATION: CLASSICAL SOLUTIONS JUHI JANG AND NING JIANG arxiv:0904.4459v [math.ap] 28 Apr 2009 Abstract. We study the acoustic limit from the Boltzmann equation in the framework
More informationPresenter: Noriyoshi Fukaya
Y. Martel, F. Merle, and T.-P. Tsai, Stability and Asymptotic Stability in the Energy Space of the Sum of N Solitons for Subcritical gkdv Equations, Comm. Math. Phys. 31 (00), 347-373. Presenter: Noriyoshi
More informationAN EXACT RESCALING VELOCITY METHOD FOR SOME KINETIC FLOCKING MODELS
AN EXACT RESCALING VELOCITY METHOD FOR SOME KINETIC FLOCKING MODELS THOMAS REY AND CHANGHUI TAN Abstract. In this work, we discuss kinetic descriptions of flocking models, of the so-called Cucker- Smale
More informationNonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems. p. 1/1
Nonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems p. 1/1 p. 2/1 Converse Lyapunov Theorem Exponential Stability Let x = 0 be an exponentially stable equilibrium
More informationFluid Dynamics from Kinetic Equations
Fluid Dynamics from Kinetic Equations François Golse Université Paris 7 & IUF, Laboratoire J.-L. Lions golse@math.jussieu.fr & C. David Levermore University of Maryland, Dept. of Mathematics & IPST lvrmr@math.umd.edu
More informationHalf of Final Exam Name: Practice Problems October 28, 2014
Math 54. Treibergs Half of Final Exam Name: Practice Problems October 28, 24 Half of the final will be over material since the last midterm exam, such as the practice problems given here. The other half
More informationDynamic programming using radial basis functions and Shepard approximations
Dynamic programming using radial basis functions and Shepard approximations Oliver Junge, Alex Schreiber Fakultät für Mathematik Technische Universität München Workshop on algorithms for dynamical systems
More informationMacroscopic models for pedestrian flow and meshfree particle methods
Macroscopic models for pedestrian flow and meshfree particle methods Axel Klar TECHNISCHE UNIVERSITÄT KAISERSLAUTERN Fachbereich Mathematik Arbeitsgruppe Technomathematik In cooperation with R. Borsche,
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 informationThe first order quasi-linear PDEs
Chapter 2 The first order quasi-linear PDEs The first order quasi-linear PDEs have the following general form: F (x, u, Du) = 0, (2.1) where x = (x 1, x 2,, x 3 ) R n, u = u(x), Du is the gradient of u.
More informationStability estimates in stationary inverse transport
Stability estimates in stationary inverse transport Guillaume Bal and Alexandre Jollivet April 8, 28 Abstract We study the stability of the reconstruction of the scattering and absorption coefficients
More informationMathematical Aspects of Self-Organized Dynamics
October 30, 205 Mathematical Aspects of Self-Organized Dynamics Consensus, Emergence of Leaders, and Social Hydrodynamics By Eitan Tadmor Simulation of agent-based V icsek dynamics and the corresponding
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 informationNavier-Stokes equations in thin domains with Navier friction boundary conditions
Navier-Stokes equations in thin domains with Navier friction boundary conditions Luan Thach Hoang Department of Mathematics and Statistics, Texas Tech University www.math.umn.edu/ lhoang/ luan.hoang@ttu.edu
More informationMcGill University Department of Mathematics and Statistics. Ph.D. preliminary examination, PART A. PURE AND APPLIED MATHEMATICS Paper BETA
McGill University Department of Mathematics and Statistics Ph.D. preliminary examination, PART A PURE AND APPLIED MATHEMATICS Paper BETA 17 August, 2018 1:00 p.m. - 5:00 p.m. INSTRUCTIONS: (i) This paper
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 informationHydrodynamic limit and numerical solutions in a system of self-pr
Hydrodynamic limit and numerical solutions in a system of self-propelled particles Institut de Mathématiques de Toulouse Joint work with: P.Degond, G.Dimarco, N.Wang 1 2 The microscopic model The macroscopic
More informationMATH 819 FALL We considered solutions of this equation on the domain Ū, where
MATH 89 FALL. The D linear wave equation weak solutions We have considered the initial value problem for the wave equation in one space dimension: (a) (b) (c) u tt u xx = f(x, t) u(x, ) = g(x), u t (x,
More informationUniformly accurate averaging numerical schemes for oscillatory evolution equations
Uniformly accurate averaging numerical schemes for oscillatory evolution equations Philippe Chartier University of Rennes, INRIA Joint work with M. Lemou (University of Rennes-CNRS), F. Méhats (University
More informationThe Dirichlet boundary problems for second order parabolic operators satisfying a Carleson condition
The Dirichlet boundary problems for second order parabolic operators satisfying a Martin Dindos Sukjung Hwang University of Edinburgh Satellite Conference in Harmonic Analysis Chosun University, Gwangju,
More informationBose-Einstein Condensation and Global Dynamics of Solutions to a Hyperbolic Kompaneets Equation
Bose-Einstein Condensation and Global Dynamics of Solutions to a Hyperbolic Kompaneets Equation Joshua Ballew Abstract In this article, a simplified, hyperbolic model of the non-linear, degenerate parabolic
More informationHomogenization with stochastic differential equations
Homogenization with stochastic differential equations Scott Hottovy shottovy@math.arizona.edu University of Arizona Program in Applied Mathematics October 12, 2011 Modeling with SDE Use SDE to model system
More informationFree energy estimates for the two-dimensional Keller-Segel model
Free energy estimates for the two-dimensional Keller-Segel model dolbeaul@ceremade.dauphine.fr CEREMADE CNRS & Université Paris-Dauphine in collaboration with A. Blanchet (CERMICS, ENPC & Ceremade) & B.
More informationInput to state Stability
Input to state Stability Mini course, Universität Stuttgart, November 2004 Lars Grüne, Mathematisches Institut, Universität Bayreuth Part III: Lyapunov functions and quantitative aspects ISS Consider with
More informationBackward Stochastic Differential Equations with Infinite Time Horizon
Backward Stochastic Differential Equations with Infinite Time Horizon Holger Metzler PhD advisor: Prof. G. Tessitore Università di Milano-Bicocca Spring School Stochastic Control in Finance Roscoff, March
More informationActive Matter Lectures for the 2011 ICTP School on Mathematics and Physics of Soft and Biological Matter Lecture 3: Hydrodynamics of SP Hard Rods
Active Matter Lectures for the 2011 ICTP School on Mathematics and Physics of Soft and Biological Matter Lecture 3: of SP Hard Rods M. Cristina Marchetti Syracuse University Baskaran & MCM, PRE 77 (2008);
More informationPropagating terraces and the dynamics of front-like solutions of reaction-diffusion equations on R
Propagating terraces and the dynamics of front-like solutions of reaction-diffusion equations on R P. Poláčik School of Mathematics, University of Minnesota Minneapolis, MN 55455 Abstract We consider semilinear
More informationIntroduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems
p. 1/5 Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems p. 2/5 Time-varying Systems ẋ = f(t, x) f(t, x) is piecewise continuous in t and locally Lipschitz in x for all t
More informationOn Schrödinger equations with inverse-square singular potentials
On Schrödinger equations with inverse-square singular potentials Veronica Felli Dipartimento di Statistica University of Milano Bicocca veronica.felli@unimib.it joint work with Elsa M. Marchini and Susanna
More informationOn Fluid-Particle Interaction
Complex Fluids On Fluid-Particle Interaction Women in Applied Mathematics University of Crete, May 2-5, 2011 Konstantina Trivisa Complex Fluids 1 Model 1. On the Doi Model: Rod-like Molecules Colaborator:
More informationRegularity of local minimizers of the interaction energy via obstacle problems
Regularity of local minimizers of the interaction energy via obstacle problems J. A. Carrillo, M. G. Delgadino, A. Mellet September 22, 2014 Abstract The repulsion strength at the origin for repulsive/attractive
More informationParticle, Kinetic, and Hydrodynamic Models of Swarming
Particle, Kinetic, and Hydrodynamic Models of Swarming J. A. Carrillo 1, M. Fornasier, G. Toscani 3, and F. Vecil 1 ICREA - Departament de Matemàtiques, Universitat Autònoma de Barcelona, 8193 Bellaterra,
More informationParticles approximation for Vlasov equation with singular interaction
Particles approximation for Vlasov equation with singular interaction M. Hauray, in collaboration with P.-E. Jabin. Université d Aix-Marseille GRAVASCO Worshop, IHP, November 2013 M. Hauray (UAM) Particles
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 informationξ,i = x nx i x 3 + δ ni + x n x = 0. x Dξ = x i ξ,i = x nx i x i x 3 Du = λ x λ 2 xh + x λ h Dξ,
1 PDE, HW 3 solutions Problem 1. No. If a sequence of harmonic polynomials on [ 1,1] n converges uniformly to a limit f then f is harmonic. Problem 2. By definition U r U for every r >. Suppose w is a
More informationChap. 1. Some Differential Geometric Tools
Chap. 1. Some Differential Geometric Tools 1. Manifold, Diffeomorphism 1.1. The Implicit Function Theorem ϕ : U R n R n p (0 p < n), of class C k (k 1) x 0 U such that ϕ(x 0 ) = 0 rank Dϕ(x) = n p x U
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 informationExercises. T 2T. e ita φ(t)dt.
Exercises. Set #. Construct an example of a sequence of probability measures P n on R which converge weakly to a probability measure P but so that the first moments m,n = xdp n do not converge to m = xdp.
More informationKPP Pulsating Traveling Fronts within Large Drift
KPP Pulsating Traveling Fronts within Large Drift Mohammad El Smaily Joint work with Stéphane Kirsch University of British olumbia & Pacific Institute for the Mathematical Sciences September 17, 2009 PIMS
More informationA new flocking model through body attitude coordination
A new flocking model through body attitude coordination Sara Merino Aceituno (Imperial College London) Pierre Degond (Imperial College London) Amic Frouvelle (Paris Dauphine) ETH Zürich, November 2016
More informationQuantitative Homogenization of Elliptic Operators with Periodic Coefficients
Quantitative Homogenization of Elliptic Operators with Periodic Coefficients Zhongwei Shen Abstract. These lecture notes introduce the quantitative homogenization theory for elliptic partial differential
More informationParticles approximation for Vlasov equation with singular interaction
Particles approximation for Vlasov equation with singular interaction M. Hauray, in collaboration with P.-E. Jabin. Université d Aix-Marseille Oberwolfach Worshop, December 2013 M. Hauray (UAM) Particles
More informationTeoría de pesos para operadores multilineales, extensiones vectoriales y extrapolación
Teoría de pesos para operadores multilineales, extensiones vectoriales y extrapolación Encuentro Análisis Funcional Sheldy Ombrosi Universidad Nacional del Sur y CONICET 8 de marzo de 208 If f is a locally
More informationKinetic models for chemotaxis
Kinetic models for chemotaxis FABIO A. C. C. CHALUB Centro de Matemática e Aplicações Fundamentais Universidade de Lisboa Kinetic models for chemotaxis p. Main Goal To study two different levels of description
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationSobolev Spaces. Chapter 10
Chapter 1 Sobolev Spaces We now define spaces H 1,p (R n ), known as Sobolev spaces. For u to belong to H 1,p (R n ), we require that u L p (R n ) and that u have weak derivatives of first order in L p
More informationLiquid crystal flows in two dimensions
Liquid crystal flows in two dimensions Fanghua Lin Junyu Lin Changyou Wang Abstract The paper is concerned with a simplified hydrodynamic equation, proposed by Ericksen and Leslie, modeling the flow of
More informationBlow-up on manifolds with symmetry for the nonlinear Schröding
Blow-up on manifolds with symmetry for the nonlinear Schrödinger equation March, 27 2013 Université de Nice Euclidean L 2 -critical theory Consider the one dimensional equation i t u + u = u 4 u, t > 0,
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 informationInverse Transport Problems (Generalized) Stability Estimates. Guillaume Bal
Inverse Transport Problems (Generalized) Stability Estimates Guillaume Bal Department of Applied Physics & Applied Mathematics Columbia University Joint with Alexandre Jollivet Data Acquisition in CT-scan
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 informationTD M1 EDP 2018 no 2 Elliptic equations: regularity, maximum principle
TD M EDP 08 no Elliptic equations: regularity, maximum principle Estimates in the sup-norm I Let be an open bounded subset of R d of class C. Let A = (a ij ) be a symmetric matrix of functions of class
More information