A new flocking model through body attitude coordination
|
|
- Constance Shields
- 5 years ago
- Views:
Transcription
1 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 Transport phenomena in collective dynamics: from micro to social hydrodynamics
2 Section 1 Introduction
3 Emergent phenomena: auto-organisation Many agents Local interactions
4 Figure: Scale bar, 10 µm. Source: A study of synchronisation between the flagella of bull spermatozoa, with related observations, David M. Woolley, Rachel F. Crockett, William D. I. Groom, Stuart G. Revell Journal of Experimental Biology :
5 Description of the system to model Vicsek model
6 Description of the system to model Vicsek model Body attitude coordination
7 Previous results for the Vicsek model Equations at different scales Vicsek model [Degond, Motsch, 08] (X i (t), ω i (t)) R n S n, i = 1,..., N dx i = cω i dt, ( c > 0 dω i = P ω ν ω i + ) 2DB i i t ν, D > 0 ω i := J i J i, J i := N K( X i (t) X j (t) )ω j Fokker-Planck type equation Relations Macroscopic eq. SOH
8 Previous results for the Vicsek model Equations at different scales Vicsek model [Degond, Motsch, 08] (X i (t), ω i (t)) R n S n, i = 1,..., N dx i = cω i dt, ( c > 0 dω i = P ω ν ω i + ) 2DB i i t ν, D > 0 ω i := J i J i, J i := N K( X i (t) X j (t) )ω j Fokker-Planck type equation f = f (t, x, ω) ( t f ε + cω x f ε = 1 ε f ω M Ω[f ] ε ω M )+O(ε) Ω[f ] M Ω (ω) exp ( ν D ω Ω) (von-mises equilibria) Ω ε [f ] S ω f ε (t, x, ω ) dω n Macroscopic eq. SOH Relations ft N = 1 N N i=1 δ x i (t),ω i (t) (empirical distribution) N [Bolley, Carrillo, Cañizo]
9 Previous results for the Vicsek model Equations at different scales Vicsek model [Degond, Motsch, 08] (X i (t), ω i (t)) R n S n, i = 1,..., N dx i = cω i dt, ( c > 0 dω i = P ω ν ω i + ) 2DB i i t ν, D > 0 ω i := J i J i, J i := N K( X i (t) X j (t) )ω j Fokker-Planck type equation f = f (t, x, ω) ( t f ε + cω x f ε = 1 ε f ω M Ω[f ] ε ω M )+O(ε) Ω[f ] M Ω (ω) exp ( ν D ω Ω) (von-mises equilibria) Ω ε [f ] S ω f ε (t, x, ω ) dω n Macroscopic eq. SOH Relations ft N = 1 N N i=1 δ x i (t),ω i (t) (empirical distribution) N [Bolley, Carrillo, Cañizo] f ε = f ε (t, x, v) [Figalli, Kang, Morales] NON conserved quant. ρ ε = f ε dv ρ ω f ε dω Ω ε Ω S n ε 0
10 Previous results for the Vicsek model Equations at different scales Vicsek model [Degond, Motsch, 08] (X i (t), ω i (t)) R n S n, i = 1,..., N dx i = cω i dt, ( c > 0 dω i = P ω ν ω i + ) 2DB i i t ν, D > 0 ω i := J i J i, J i := N K( X i (t) X j (t) )ω j Fokker-Planck type equation f = f (t, x, ω) ( t f ε + cω x f ε = 1 ε f ω M Ω[f ] ε ω M )+O(ε) Ω[f ] M Ω (ω) exp ( ν D ω Ω) (von-mises equilibria) Ω ε [f ] S ω f ε (t, x, ω ) dω n Macroscopic eq. SOH ρ = ρ(t, x), Ω = Ω(t, x) { t ρ + x (c 1 ρω) = 0, ρ ( t Ω + c 2 (Ω x )Ω) + c 3 P Ω ( x ρ) = 0 Relations ft N = 1 N N i=1 δ x i (t),ω i (t) (empirical distribution) N [Bolley, Carrillo, Cañizo] f ε = f ε (t, x, v) [Figalli, Kang, Morales] NON conserved quant. ρ ε = f ε dv ρ ω f ε dω Ω ε Ω S n ε 0 f ε f = ρ(t, x)m Ω (ω) (local equilibria) [Jiang, Xiong, Zhang] Generalised Collision Invariant Simulation by Sébastien Motsch (ASU)
11 Section 2 New results: body attitude coordination
12 Derivation of the Individual Based Model (X i, A i ) i {1,...,N}, X i R 3, A i. Averaged body orientation: M k (t) := 1 N N K( X i (t) X k (t) )A i (t). (1) i=1
13 Derivation of the Individual Based Model (X i, A i ) i {1,...,N}, X i R 3, A i. Averaged body orientation: M k (t) := 1 N N K( X i (t) X k (t) )A i (t). (1) i=1 Lemma (Polar decomposition of a square matrix, Golub 12) Given a matrix M M, if det(m) 0 then there exists a unique orthogonal matrix A (given by A = M( M T M) 1 ) and a unique symmetric positive definite matrix S such that M = AS. The matrix A minimizes the quantity 1 N N i=1 K( X i(t) X k (t) ) A i (t) A 2 among the elements of.
14 Derivation of the Individual Based Model (X i, A i ) i {1,...,N}, X i R 3, A i. Averaged body orientation: M k (t) := 1 N N K( X i (t) X k (t) )A i (t). (1) i=1 Lemma (Polar decomposition of a square matrix, Golub 12) Given a matrix M M, if det(m) 0 then there exists a unique orthogonal matrix A (given by A = M( M T M) 1 ) and a unique symmetric positive definite matrix S such that M = AS. The matrix A minimizes the quantity 1 N N i=1 K( X i(t) X k (t) ) A i (t) A 2 among the elements of. The equations: dx k (t) = v 0 A k (t)e 1 dt (2) [ da k (t) = P TAk νpd(m k )dt + 2 ] DdWt k (3)
15 Mean-field limit Proposition (Formal) When the number of agents in (2)-(3) N, its corresponding empirical distribution f N (x, A, t) = 1 N N k=1 δ (Xk (t),a k (t)) converges weakly to f = f (x, A, t), (x, A, t) R 3 [0, ) satisfying t f + v 0 Ae 1 x f = D A f A (ff [f ]) (4) F [f ] := νp TA ( M[f ]) M[f ] = PD(M[f ]), M[f ](x, t) := R 3 K ( x x ) f (x, A, t)a da dx where PD(M[f ]) corresponds to the orthogonal matrix obtained on the Polar Decomposition of M[f ].
16 Macroscopic limit Adimensional analysis and rescaling: t f ε + Ae 1 x f ε = 1 ε Q(f ε ) + O(ε) (5) F 0 [f ] := ν P TA (Λ[f ]) Λ[f ] = PD(λ[f ]), λ[f ](x, t) := Q(f ) := D A f A (ff 0 [f ]) f (x, A, t)a da
17 Macroscopic limit Adimensional analysis and rescaling: t f ε + Ae 1 x f ε = 1 ε Q(f ε ) + O(ε) (5) F 0 [f ] := ν P TA (Λ[f ]) Λ[f ] = PD(λ[f ]), λ[f ](x, t) := Q(f ) := D A f A (ff 0 [f ]) = D A with equilibria given by ρ(t, x)m Λ (A): M Λ (A) = 1 ( ) ν(a Λ) Z exp, D f (x, A, t)a da [ M Λ[f ] A ( Consistency relation: λ[m Λ0 ] = c 1 Λ 0 where c 1 (0, 1). f M Λ[f ] )]. M Λ (A) da = 1, Λ, (6)
18 Macroscopic limit Adimensional analysis and rescaling: t f ε + Ae 1 x f ε = 1 ε Q(f ε ) + O(ε) (5) F 0 [f ] := ν P TA (Λ[f ]) Λ[f ] = PD(λ[f ]), λ[f ](x, t) := Q(f ) := D A f A (ff 0 [f ]) = D A with equilibria given by ρ(t, x)m Λ (A): M Λ (A) = 1 ( ) ν(a Λ) Z exp, D f (x, A, t)a da [ M Λ[f ] A ( Consistency relation: λ[m Λ0 ] = c 1 Λ 0 where c 1 (0, 1). As ε 0, f ε f = f (t, x, A) = ρm Λ (A), Λ = Λ(t, x), ρ = ρ(t, x) 0. f M Λ[f ] )]. M Λ (A) da = 1, Λ, (6)
19 Key steps of the proof: the continuity equation and the non-conservation of momentum t f ε + Ae 1 x f ε = 1 ε Q(f ε ) + O(ε) By conservation of mass: t ρ ε + x Ae 1 f ε da = O(ε), ρ ε (t, x) := } {{ } =:j[f ε] in the limit (formally) ρ ε ρ (remember that f ε ρm Λ ) j[f ε ] ρ Ae 1 M Λ (A) da =: j by the consistency relation, we have that Conclude: j = ρc 1 Λe 1. t ρ(t, x) + x (ρ(t, x) c 1 Λ(t, x)e 1 ) = 0. f ε (x, A, t) da
20 Key steps of the proof: the continuity equation and the non-conservation of momentum t f ε + Ae 1 x f ε = 1 ε Q(f ε ) + O(ε) Equation for Λ? By the consistency relation: AM Λ (A) da = c 1 Λ but the momentum is not conserved! 1 AQ(f ε )(A) da 0 ε
21 The Generalised Collision Invariant Define the operator notice in particular that Q(f, Λ 0 ) := A ( ( )) f M Λ0 A M Λ0 Q(f ) = Q(f, Λ[f ]).
22 The Generalised Collision Invariant Define the operator notice in particular that Q(f, Λ 0 ) := A Using this operator we define: ( ( )) f M Λ0 A M Λ0 Q(f ) = Q(f, Λ[f ]). Definition (Generalised Collision Invariant) For a given Λ 0 we say that a real-valued function ψ : R is a Generalised Collision Invariant associated to Λ 0, ψ GCI (Λ 0 ) for short, if Q(f, Λ 0 )ψ da = 0 for all f s.t P TΛ0 (λ[f ]) = 0.
23 Why this helps? Consider ψ f GCI (Λ f ), then [ t f ε + Ae 1 x f ε] ψ f da = 1 ε = 1 ε = O(ε) Q(f ε )ψ f da + O(ε) Q(f ε, Λ f )ψ f da + O(ε)
24 Why this helps? Consider ψ f GCI (Λ f ), then [ t f ε + Ae 1 x f ε] ψ f da = 1 ε = 1 ε = O(ε) Q(f ε )ψ f da + O(ε) Q(f ε, Λ f )ψ f da + O(ε) Reference: Ning Jiang, Linjie Xiong, Teng-Fei Zhang, Hydrodynamic limits of the kinetic self-organized models; Existence of the GCI: Lax-Milgram; Computation of the GCI: differential geometry in
25 Computation of the GCI Proposition (Equation for the GCI) We have that ψ GCI (Λ 0 ) if and only if exists C T Λ0 such that A (M Λ0 A ψ) = C AM Λ0. (7)
26 Computation of the GCI Proposition (Equation for the GCI) We have that ψ GCI (Λ 0 ) if and only if exists C T Λ0 such that A (M Λ0 A ψ) = C AM Λ0. (7) Proposition (Existence of the GCI) For a given B T Λ fixed, there exists a unique ψ B H 1 () up to a constant, satisfying the relation (7).
27 Computation of the GCI Proposition (Equation for the GCI) We have that ψ GCI (Λ 0 ) if and only if exists C T Λ0 such that A (M Λ0 A ψ) = C AM Λ0. (7) Proposition (Existence of the GCI) For a given B T Λ fixed, there exists a unique ψ B H 1 () up to a constant, satisfying the relation (7). Proposition (Change of variables) Let ψ(b) := ψ(λb). Then ψ satisfies B (M Id B ψ ) = P BM Id, P A. (8)
28 Proposition (Non-constant GCI) Let P A, then the unique solution ψ H0 1 () of (8) is given by Going back to the GCI ψ(a), we can write it as ψ(b) = P B ψ 0 ( 1 2tr(B)). (9) ψ(a) = P (Λ T A) ψ 0 (Λ A) (10) where ψ 0 is constructed as follows: Let ψ 0 : R R be the unique solution to 1 )) sin 2 (θ/2) θ (sin 2 m(θ) sin θ (θ/2)m(θ) θ (sin θ ψ 0 ψ 2 sin 2 0 = sin θm(θ) θ/2 (11) where m(θ) = M Id (B) = exp(d 1 σ( cos θ))/z. Then ψ 0 (θ) = ψ 0 ( 1 2 tr(b)) by the relation 1 2 tr(b) = cos θ. ψ 0 is 2π-periodic, even and negative.
29 Back to the limit t f ε + Ae 1 x f ε = 1 ε Q(f ε ) + O(ε). With GCI ψλ P(A) = (P ΛT A)ψ 0 (Λ A) for any P antisymmetric: [ ] [ t (ρm Λ ) + Ae 1 x (ρm Λ )] P Λ T A ψ 0 (Λ A) da = 0 for all P A therefore X := since [ t (ρm Λ ) + Ae 1 x (ρm Λ )] (Λ T A A T Λ)ψ 0 (Λ A) da = 0. M = P S P: antisymmetric matrices, S, symmetric matrices. Compute the value of each term.
30 The macroscopic equations Theorem ((Formal) macroscopic limit) When ε 0 in the kinetic equation it holds that f ε f = f (x, A, t) = ρm Λ (A), Λ = Λ(t, x), ρ = ρ(t, x) 0. Consider the orthogonal basis {Λe 1 =: Ω, Λe 2 =: u, Λe 3 =: v}, where {e 1, e 2, e 3 } is the canonical basis. The equation for Λ expressed as the evolution of this basis is ρd t Ω + P Ω (c 3 x ρ + c 4 ρr) = 0, (12) ρd t u u (c 3 x ρ + c 4 ρr) Ω + c 4 ρδv = 0, (13) ρd t v v (c 3 x ρ + c 4 ρr) Ω c 4 ρδu = 0, (14) where D t := t + c 2 (Ω x ), (15) δ := [(Ω x )u] v + [(u x )v] Ω + [(v x )Ω] u, (16) r := ( x Ω)Ω + ( x u)u + ( x v)v. (17)
31 Proposition (Properties of δ and r) Let Λ = Λ(x) : R 3 be a function taking values on. It holds that δ Λ = δ ΛA and r Λ = r ΛA for any A. (18) where δ Λ and r Λ correspond to the expressions in (16) and (17), respectively. Suppose now that Λ = Λ(x) = exp ([b(x)] ) Λ(x 0 ) where b = b(x) : R 3 R 3 is a smooth function around a fixed x 0 R 3 where b(x 0 ) = 0, then δ Λ (x 0 ) = ( x b) (x 0 ), r Λ (x 0 ) = ( x b) (x 0 ).
32 Future work Simulations; new modelling using quaternions: joint with Pierre Degond (Imperial College London), Amic Frouvelle (Université Paris Dauphine) and Ariane Trescases (University of Cambridge). Study emergent structures. Compare them with the ones appearing in the Vicsek model (sperm dynamics). Some preliminary results: Maciek Biskupiak s video. Simulation by Maciek Biskupiak
33 Special thanks to: THE ORGANISERS Gianluca Crippa (Universität Basel) Siddhartha Mishra (ETH Zrich) Eitan Tadmor (University of Maryland / ETH-ITS Zrich) AND Andrea Bianda and Anil Zenginoglu AND Sponsors: EPSRC (Engineering and Physical Sciences Research Council, UK) & KI-NET & FIM at ETH Zürich
34 Summary: body attitude coordination Equations at different scales Individual Based Model (X i (t), A i (t)) R 3, i = 1,..., N dx i = v 0 A i e 1 dt, ( v 0 > 0 da i = P TAi νpd(m k )dt + 2 ) DdWt i M i (t) := 1 N N k=1 K( X k(t) X i (t) )A k (t). Fokker-Planck type equation f = f (t, x, A) ( t f ε + cω x f ε = 1 ε f A M Λ[f ] ε A M )+O(ε) Λ[f ] M Λ (A) exp ( ν D A Λ) (von-mises eq. ) Λ ε [f ] A f ε (t, x, A ) da Macroscopic eq. ρ = ρ(t, x), Λ = Λ(t, x) t ρ + x (c 1 ρλe 1 ) = 0, ρd t Ω + P Ω (c 3 x ρ + c 4 ρr) = 0, ρd t u u (c 3 x ρ + c 4 ρr) Ω + c 4 ρδv = 0, ρd t v v (c 3 x ρ + c 4 ρr) Ω c 4 ρδu = 0, Relations ft N = 1 N N i=1 δ X i (t),a i (t) (empirical distribution) N f ε = f ε (t, x, v) ρ ε = f ε dv ρ NON conserved quant. Consistency relation: A f ε da = c 1 Λ ε Λ ε 0 f ε f = ρ(t, x)m Λ (A)
A new flocking model through body attitude coordination
A new flocking model through body attitude coordination arxiv:1605.03509v1 [math-ph] 11 May 016 Pierre Degond 1, Amic Frouvelle, and Sara Merino-Aceituno 3 13 Department of Mathematics, Imperial College
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 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 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 informationA model of alignment interaction for oriented particles with phase transition
A model of alignment interaction for oriented particles with phase transition Amic Frouvelle Institut de Mathématiques de Toulouse Joint work with Jian-Guo Liu (Duke Univ.) and Pierre Degond (IMT) Amic
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 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 informationarxiv: v1 [math.ap] 16 Oct 2018
Alignment of self-propelled rigid bodies: from particle systems to macroscopic equations arxiv:1810.06903v1 [math.ap] 16 Oct 2018 Pierre Degond 1, Amic Frouvelle 2, Sara Merino-Aceituno 3, and Ariane Trescases
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 informationPhase transitions, hysteresis, and hyperbolicity for self-organized alignment dynamics
Phase transitions, hysteresis, and hyperbolicity for self-organized alignment dynamics Pierre Degond (,), Amic Frouvelle (3), Jian-Guo Liu (4) - Université de Toulouse; UPS, INSA, UT, UT; Institut de athématiques
More informationAnomalous energy transport in FPU-β chain
Anomalous energy transport in FPU-β chain Sara Merino Aceituno Joint work with Antoine Mellet (University of Maryland) http://arxiv.org/abs/1411.5246 Imperial College London 9th November 2015. Kinetic
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 informationA note on phase transitions for the Smoluchowski equation with dipolar potential
A note on phase transitions for the Smoluchowski equation with dipolar potential Pierre Degond (,2), Amic Frouvelle (3), Jian-Guo Liu (4) arxiv:22.3920v [math-ph] 7 Dec 202 -Université de Toulouse; UPS,
More informationarxiv: v1 [math.na] 14 Jan 2019
Noname manuscript No. will be inserted by the editor) Kinetic equations and self-organized band formations Quentin Griette Sebastien Motsch arxiv:191.43v1 [math.na] 14 Jan 19 the date of receipt and acceptance
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 informationPhase Transitions, Hysteresis, and Hyperbolicity for Self-Organized Alignment Dynamics
Arch. Rational ech. Anal. 6 (05) 63 5 Digital Object Identifier (DOI) 0.007/s0005-04-0800-7 Phase Transitions, Hysteresis, and Hyperbolicity for Self-Organized Alignment Dynamics Pierre Degond, Amic Frouvelle
More informationMicro-macro methods for Boltzmann-BGK-like equations in the diffusion scaling
Micro-macro methods for Boltzmann-BGK-like equations in the diffusion scaling Anaïs Crestetto 1, Nicolas Crouseilles 2, Giacomo Dimarco 3 et Mohammed Lemou 4 Saint-Malo, 14 décembre 2017 1 Université de
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 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 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 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 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 informationHamilton-Jacobi theory on Lie algebroids: Applications to nonholonomic mechanics. Manuel de León Institute of Mathematical Sciences CSIC, Spain
Hamilton-Jacobi theory on Lie algebroids: Applications to nonholonomic mechanics Manuel de León Institute of Mathematical Sciences CSIC, Spain joint work with J.C. Marrero (University of La Laguna) D.
More informationAn asymptotic-preserving micro-macro scheme for Vlasov-BGK-like equations in the diffusion scaling
An asymptotic-preserving micro-macro scheme for Vlasov-BGK-like equations in the diffusion scaling Anaïs Crestetto 1, Nicolas Crouseilles 2 and Mohammed Lemou 3 Saint-Malo 13 December 2016 1 Université
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 informationNumerical Simulation of Herding and Flocking Models
ERASMUS MUNDUS MASTER OF INDUSTRIAL MATHEMATICS Numerical Simulation of Herding and Flocking Models Robert Setekera Supervisor: Prof. Dr. Axel Klar A Thesis submitted for the degree of Master of Science
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 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 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 informationBoltzmann Equation and Hydrodynamics beyond Navier-Stokes
Boltzmann Equation and Hydrodynamics beyond Navier-Stokes Alexander Bobylev Keldysh Institute for Applied Mathematics, RAS, Moscow Chapman-Enskog method and Burnett equations Notation: f (x, v, t) - distribution
More informationNumerical methods for kinetic equations
Numerical methods for kinetic equations Lecture 6: fluid-kinetic coupling and hybrid methods Lorenzo Pareschi Department of Mathematics and Computer Science University of Ferrara, Italy http://www.lorenzopareschi.com
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 informationMicrolocal analysis and inverse problems Lecture 4 : Uniqueness results in admissible geometries
Microlocal analysis and inverse problems Lecture 4 : Uniqueness results in admissible geometries David Dos Santos Ferreira LAGA Université de Paris 13 Wednesday May 18 Instituto de Ciencias Matemáticas,
More informationFOUR-WAY COUPLED SIMULATIONS OF TURBULENT
FOUR-WAY COUPLED SIMULATIONS OF TURBULENT FLOWS WITH NON-SPHERICAL PARTICLES Berend van Wachem Thermofluids Division, Department of Mechanical Engineering Imperial College London Exhibition Road, London,
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 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 informationFourier Law and Non-Isothermal Boundary in the Boltzmann Theory
in the Boltzmann Theory Joint work with Raffaele Esposito, Yan Guo, Rossana Marra DPMMS, University of Cambridge ICERM November 8, 2011 Steady Boltzmann Equation Steady Boltzmann Equation v x F = Q(F,
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 informationOn the interplay between kinetic theory and game theory
1 On the interplay between kinetic theory and game theory Pierre Degond Department of mathematics, Imperial College London pdegond@imperial.ac.uk http://sites.google.com/site/degond/ Joint work with J.
More informationDecompositions of Modules and Comodules
Decompositions of Modules and Comodules Robert Wisbauer University of Düsseldorf, Germany Abstract It is well-known that any semiperfect A ring has a decomposition as a direct sum (product) of indecomposable
More informationEffective Theories and Minimal Energy Configurations for Heterogeneous Multilayers
Effective Theories and Minimal Energy Configurations for Universität Augsburg, Germany Minneapolis, May 16 th, 2011 1 Overview 1 Motivation 2 Overview 1 Motivation 2 Effective Theories 2 Overview 1 Motivation
More informationKinetic/Fluid micro-macro numerical scheme for Vlasov-Poisson-BGK equation using particles
Kinetic/Fluid micro-macro numerical scheme for Vlasov-Poisson-BGK equation using particles Anaïs Crestetto 1, Nicolas Crouseilles 2 and Mohammed Lemou 3. Workshop Asymptotic-Preserving schemes, Porquerolles.
More informationarxiv: v1 [math.ap] 1 Feb 2015
GLOBAL WEAK SOLUTIONS FOR KOLMOGOROV-VICSEK TYPE EQUATIONS WITH ORIENTATIONAL INTERACTION IRENE M. GAMBA AN MOON-JIN KANG arxiv:1502.00293v1 [math.ap] 1 Feb 2015 Abstract. We rove the global existence
More informationPhenomena in high dimensions in geometric analysis, random matrices, and computational geometry Roscoff, France, June 25-29, 2012
Phenomena in high dimensions in geometric analysis, random matrices, and computational geometry Roscoff, France, June 25-29, 202 BOUNDS AND ASYMPTOTICS FOR FISHER INFORMATION IN THE CENTRAL LIMIT THEOREM
More informationEDP with strong anisotropy : transport, heat, waves equations
EDP with strong anisotropy : transport, heat, waves equations Mihaï BOSTAN University of Aix-Marseille, FRANCE mihai.bostan@univ-amu.fr Nachos team INRIA Sophia Antipolis, 3/07/2017 Main goals Effective
More informationEmergence of collective dynamics from a purely stochastic orogin
Emergence of collective dynamics from a purely stochastic orogin Yann Brenier CNRS-Centre de Mathématiques Laurent SCHWARTZ Ecole Polytechnique FR 91128 Palaiseau. Transport phenomena in collective dynamics,
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 informationNon-equilibrium mixtures of gases: modeling and computation
Non-equilibrium mixtures of gases: modeling and computation Lecture 1: and kinetic theory of gases Srboljub Simić University of Novi Sad, Serbia Aim and outline of the course Aim of the course To present
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 informationTD 1: Hilbert Spaces and Applications
Université Paris-Dauphine Functional Analysis and PDEs Master MMD-MA 2017/2018 Generalities TD 1: Hilbert Spaces and Applications Exercise 1 (Generalized Parallelogram law). Let (H,, ) be a Hilbert space.
More informationFast-slow systems with chaotic noise
Fast-slow systems with chaotic noise David Kelly Ian Melbourne Courant Institute New York University New York NY www.dtbkelly.com May 1, 216 Statistical properties of dynamical systems, ESI Vienna. David
More informationUne décomposition micro-macro particulaire pour des équations de type Boltzmann-BGK en régime de diffusion
Une décomposition micro-macro particulaire pour des équations de type Boltzmann-BGK en régime de diffusion Anaïs Crestetto 1, Nicolas Crouseilles 2 et Mohammed Lemou 3 La Tremblade, Congrès SMAI 2017 5
More informationEXTREMELY CHARGED STATIC DUST DISTRIBUTIONS IN GENERAL RELATIVITY
arxiv:gr-qc/9806038v1 8 Jun 1998 EXTREMELY CHARGED STATIC DUST DISTRIBUTIONS IN GENERAL RELATIVITY METÍN GÜRSES Mathematics department, Bilkent University, 06533 Ankara-TURKEY E-mail: gurses@fen.bilkent.edu.tr
More informationDynamical Systems & Lyapunov Stability
Dynamical Systems & Lyapunov Stability Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University Outline Ordinary Differential Equations Existence & uniqueness Continuous dependence
More informationLarge-scale atmospheric circulation, semi-geostrophic motion and Lagrangian particle methods
Large-scale atmospheric circulation, semi-geostrophic motion and Lagrangian particle methods Colin Cotter (Imperial College London) & Sebastian Reich (Universität Potsdam) Outline 1. Hydrostatic and semi-geostrophic
More informationQuantum Hydrodynamics models derived from the entropy principle
1 Quantum Hydrodynamics models derived from the entropy principle P. Degond (1), Ch. Ringhofer (2) (1) MIP, CNRS and Université Paul Sabatier, 118 route de Narbonne, 31062 Toulouse cedex, France degond@mip.ups-tlse.fr
More informationHonors Differential Equations
MIT OpenCourseWare http://ocw.mit.edu 18.034 Honors Differential Equations Spring 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. UNIT III: HIGHER-ORDER
More informationLecture 41: Highlights
Lecture 41: Highlights The goal of this lecture is to remind you of some of the key points that we ve covered this semester Note that this is not the complete set of topics that may appear on the final
More informationLecture Models for heavy-ion collisions (Part III): transport models. SS2016: Dynamical models for relativistic heavy-ion collisions
Lecture Models for heavy-ion collisions (Part III: transport models SS06: Dynamical models for relativistic heavy-ion collisions Quantum mechanical description of the many-body system Dynamics of heavy-ion
More informationFrom the Newton equation to the wave equation in some simple cases
From the ewton equation to the wave equation in some simple cases Xavier Blanc joint work with C. Le Bris (EPC) and P.-L. Lions (Collège de France) Université Paris Diderot, FRACE http://www.ann.jussieu.fr/
More informationWaves in plasma. Denis Gialis
Waves in plasma Denis Gialis This is a short introduction on waves in a non-relativistic plasma. We will consider a plasma of electrons and protons which is fully ionized, nonrelativistic and homogeneous.
More information(TRAVELLING) 1D WAVES. 1. Transversal & Longitudinal Waves
(TRAVELLING) 1D WAVES 1. Transversal & Longitudinal Waves Objectives After studying this chapter you should be able to: Derive 1D wave equation for transversal and longitudinal Relate propagation speed
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 informationSRB measures for non-uniformly hyperbolic systems
SRB measures for non-uniformly hyperbolic systems Vaughn Climenhaga University of Maryland October 21, 2010 Joint work with Dmitry Dolgopyat and Yakov Pesin 1 and classical results Definition of SRB measure
More informationBrownian motion and the Central Limit Theorem
Brownian motion and the Central Limit Theorem Amir Bar January 4, 3 Based on Shang-Keng Ma, Statistical Mechanics, sections.,.7 and the course s notes section 6. Introduction In this tutorial we shall
More informationTransport Continuity Property
On Riemannian manifolds satisfying the Transport Continuity Property Université de Nice - Sophia Antipolis (Joint work with A. Figalli and C. Villani) I. Statement of the problem Optimal transport on Riemannian
More informationON A NEW CLASS OF WEAK SOLUTIONS TO THE SPATIALLY HOMOGENEOUS BOLTZMANN AND LANDAU EQUATIONS
ON A NEW CLASS OF WEAK SOLUTIONS TO THE SPATIALLY HOMOGENEOUS BOLTZMANN AND LANDAU EQUATIONS C. VILLANI Abstract. This paper deals with the spatially homogeneous Boltzmann equation when grazing collisions
More informationKinetic 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 informationVariational approach to mean field games with density constraints
1 / 18 Variational approach to mean field games with density constraints Alpár Richárd Mészáros LMO, Université Paris-Sud (based on ongoing joint works with F. Santambrogio, P. Cardaliaguet and F. J. Silva)
More informationEffective dynamics for the (overdamped) Langevin equation
Effective dynamics for the (overdamped) Langevin equation Frédéric Legoll ENPC and INRIA joint work with T. Lelièvre (ENPC and INRIA) Enumath conference, MS Numerical methods for molecular dynamics EnuMath
More 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 informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
More informationGlobal Weak Solutions for Kolmogorov Vicsek Type Equations with Orientational Interactions
Arch. Rational Mech. Anal. 222 (2016) 317 342 igital Object Identifier (OI) 10.1007/s00205-016-1002-2 Global Weak Solutions for Kolmogorov Vicsek Tye Equations with Orientational Interactions Irene M.
More informationInformation and Credit Risk
Information and Credit Risk M. L. Bedini Université de Bretagne Occidentale, Brest - Friedrich Schiller Universität, Jena Jena, March 2011 M. L. Bedini (Université de Bretagne Occidentale, Brest Information
More informationCaltech Ph106 Fall 2001
Caltech h106 Fall 2001 ath for physicists: differential forms Disclaimer: this is a first draft, so a few signs might be off. 1 Basic properties Differential forms come up in various parts of theoretical
More informationLecture 5: Kinetic theory of fluids
Lecture 5: Kinetic theory of fluids September 21, 2015 1 Goal 2 From atoms to probabilities Fluid dynamics descrines fluids as continnum media (fields); however under conditions of strong inhomogeneities
More informationAPPENDIX Z. USEFUL FORMULAS 1. Appendix Z. Useful Formulas. DRAFT 13:41 June 30, 2006 c J.D Callen, Fundamentals of Plasma Physics
APPENDIX Z. USEFUL FORMULAS 1 Appendix Z Useful Formulas APPENDIX Z. USEFUL FORMULAS 2 Key Vector Relations A B = B A, A B = B A, A A = 0, A B C) = A B) C A B C) = B A C) C A B), bac-cab rule A B) C D)
More informationUne décomposition micro-macro particulaire pour des équations de type Boltzmann-BGK en régime de diffusion
Une décomposition micro-macro particulaire pour des équations de type Boltzmann-BGK en régime de diffusion Anaïs Crestetto 1, Nicolas Crouseilles 2 et Mohammed Lemou 3 Rennes, 14ème Journée de l équipe
More informationElectromagnetic waves in magnetized plasma The dispersion relation
Electromagnetic waves in magnetized plasma The dispersion relation Bruno Després (LJLL-UPMC) Electromagnetic waves in magnetized plasma The dispersion relation p. 1 / 32 Vlasov-Maxwell Vectors are in bold.
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 informationHydrodynamic Limit with Geometric Correction in Kinetic Equations
Hydrodynamic Limit with Geometric Correction in Kinetic Equations Lei Wu and Yan Guo KI-Net Workshop, CSCAMM University of Maryland, College Park 2015-11-10 1 Simple Model - Neutron Transport Equation
More informationTutorial on Control and State Constrained Optimal Control Pro. Control Problems and Applications Part 3 : Pure State Constraints
Tutorial on Control and State Constrained Optimal Control Problems and Applications Part 3 : Pure State Constraints University of Münster Institute of Computational and Applied Mathematics SADCO Summer
More informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 23 March 2001 with E. A. Spiegel
More informationTraveling waves dans le modèle de Kuramoto quenched
Traveling waves dans le modèle de Kuramoto quenched Eric Luçon MAP5 - Université Paris Descartes Grenoble - Journées MAS 2016 31 août 2016 Travail en commun avec Christophe Poquet (Lyon 1) [Giacomin, L.,
More informationp 1 ( Y p dp) 1/p ( X p dp) 1 1 p
Doob s inequality Let X(t) be a right continuous submartingale with respect to F(t), t 1 P(sup s t X(s) λ) 1 λ {sup s t X(s) λ} X + (t)dp 2 For 1 < p
More 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 informationMacroscopic limits and phase transition in a system of self-propelled particles
Macroscopic limits and phase transition in a system of self-propelled particles Pierre Degond (,), Amic Frouvelle (,), Jian-Guo Liu (3) arxiv:09.404v [math-ph] ep 0 -Université de Toulouse; UP, INA, UT,
More informationSet 3: Galaxy Evolution
Set 3: Galaxy Evolution Environment. Galaxies are clustered, found in groups like the local group up to large clusters of galaxies like the Coma cluster Small satellite galaxies like the LMC and SMC are
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 information13 The martingale problem
19-3-2012 Notations Ω complete metric space of all continuous functions from [0, + ) to R d endowed with the distance d(ω 1, ω 2 ) = k=1 ω 1 ω 2 C([0,k];H) 2 k (1 + ω 1 ω 2 C([0,k];H) ), ω 1, ω 2 Ω. F
More informationMATH & MATH FUNCTIONS OF A REAL VARIABLE EXERCISES FALL 2015 & SPRING Scientia Imperii Decus et Tutamen 1
MATH 5310.001 & MATH 5320.001 FUNCTIONS OF A REAL VARIABLE EXERCISES FALL 2015 & SPRING 2016 Scientia Imperii Decus et Tutamen 1 Robert R. Kallman University of North Texas Department of Mathematics 1155
More informationSurfaces of Arbitrary Constant Negative Gaussian Curvature and Related Sine-Gordon Equations
Mathematica Aeterna, Vol.1, 011, no. 01, 1-11 Surfaces of Arbitrary Constant Negative Gaussian Curvature and Related Sine-Gordon Equations Paul Bracken Department of Mathematics, University of Texas, Edinburg,
More informationarxiv: v7 [quant-ph] 22 Aug 2017
Quantum Mechanics with a non-zero quantum correlation time Jean-Philippe Bouchaud 1 1 Capital Fund Management, rue de l Université, 75007 Paris, France. (Dated: October 8, 018) arxiv:170.00771v7 [quant-ph]
More informationHydrodynamics, Thermodynamics, and Mathematics
Hydrodynamics, Thermodynamics, and Mathematics Hans Christian Öttinger Department of Mat..., ETH Zürich, Switzerland Thermodynamic admissibility and mathematical well-posedness 1. structure of equations
More informationHomogenization of micro-resonances and localization of waves.
Homogenization of micro-resonances and localization of waves. Valery Smyshlyaev University College London, UK July 13, 2012 (joint work with Ilia Kamotski UCL, and Shane Cooper Bath/ Cardiff) Valery Smyshlyaev
More informationApplications of the periodic unfolding method to multi-scale problems
Applications of the periodic unfolding method to multi-scale problems Doina Cioranescu Université Paris VI Santiago de Compostela December 14, 2009 Periodic unfolding method and multi-scale problems 1/56
More informationHybrid and Moment Guided Monte Carlo Methods for Kinetic Equations
Hybrid and Moment Guided Monte Carlo Methods for Kinetic Equations Giacomo Dimarco Institut des Mathématiques de Toulouse Université de Toulouse France http://perso.math.univ-toulouse.fr/dimarco giacomo.dimarco@math.univ-toulouse.fr
More informationA diamagnetic inequality for semigroup differences
A diamagnetic inequality for semigroup differences (Irvine, November 10 11, 2001) Barry Simon and 100DM 1 The integrated density of states (IDS) Schrödinger operator: H := H(V ) := 1 2 + V ω =: H(0, V
More informationEigenvalues and Eigenvectors
/88 Chia-Ping Chen Department of Computer Science and Engineering National Sun Yat-sen University Linear Algebra Eigenvalue Problem /88 Eigenvalue Equation By definition, the eigenvalue equation for matrix
More information