The Form of Swarms: A Primer of Swarm Equilibria
|
|
- David York
- 5 years ago
- Views:
Transcription
1 The Form of Swarms: A Primer of Swarm Equilibria Andrew J. Bernoff Harvey Mudd College Chad M. Topaz Macalester College Supported by: National Science Foundation DMS
2 Andrew J. Bernoff Harvey Mudd College Chad M. Topaz Macalester College Abstract: Biological aggregations (swarms) ehibit morphologies governed by social interactions and responses to environment. Starting from a particle model we derive a nonlocal PDE, known as the aggregation equation, which describes evolving population density. The solutions to the aggregation equation can ehibit a variety of behaviors including spreading without bound, concentrating into δ- functions and formation of compactly supported equilibria. We describe some tools for investigating the asymptotic behavior of solutions. We also study equilibria and their stability via the calculus of variations which yields analytical solutions. Finally we present a case study about how these methods can be used to construct a model of locust swarms. Useful References: A. J. Bernoff & C. M. Topaz, Nonlocal aggregation models: A primer of swarm equilibria. SIAM Review 55 (2013) Note the introduction and epilogue to this paper contain a potentially useful literature review. A. J. Leverentz, C. M. Topaz & A. J. Bernoff. Asymptotic Dynamics of Attractive-Repulsive Swarms, SIAM J. Appl. Dyn. Sys. 8 (2009) C. M. Topaz, A. J. Bernoff, S. Logan & W. Toolson, Aggregations, Interactions, and Boundaries: A Minimal Model for Rolling Swarms of Locusts, Eur. Phys. J. - Spec. Top. 157 (2008)
3 Many aggregations have sharp boundaries. Sinclair, 1977
4 Biological question How are individual behavior and group behavior connected? Mathematical question How do different social interaction kernels affect the (asymptotic) macroscopic behavior of a solution?
5 A Continuum Model of Swarming Three Key Ideas: Continuum model of swarm density,. Individuals interact pairwise and social forces are additive. Repel at short distances (collision avoidance) Attract at longer distances (swarming behavior) Kinematic Model (velocity social forces) Previous Work: Bodnar & Velazquez (2006) considered well-posedness and blow-up in 1D. Bertozzi & Laurent (2007) proved aisymmetric blow-up in N-dimensions. Bertozzi & collaborators ( ) working on general blow-up in N-dimensions.
6 Computing the Velocity (y) Velocity for a test mass at v() = Z 1 y= 1 dy from the density (y) f s ( y) (y) dy f s Where f s () is the social force induced by the mass distribution with density (y). f s ()
7 Conservation of Mass J J The flu, J, is given by: J = v So the equation of conservation of @ =0 ) t +(v ) =0 And the velocity is given by v = Z 1 1 f s ( y) (y) dy f s f s ()
8 Conservation of Mass (y) The flu, J, is given by: J = v So the equation of conservation of @ =0 ) t +(v ) =0 And the velocity is given by v = Z 1 1 f s ( y) (y) dy f s f s ()
9 Continuum vs. Discrete Continuum Model: Difficult to simulate, particularly, in higher dimensions Simulation difficulty independent of number of particles. Relationship to biology not always clear. Much better tools for analysis Discrete Model: Easy to simulate a small number of particles. Easy to generalize to higher dimensions. Relationship to biology well understood. Difficult to analyze or find eact solutions. z (b)
10 Meet Andy. SIAM J. APPLIED DYNAMICAL SYSTEMS Vol. 8, No Society for Industrial and Applied Mathematics Asymptotic Dynamics of Attractive-Repulsive Swarms Andrew J. Leverentz, Chad M. Topaz, and Andrew J. Bernoff Andy Leverentz HMC 2008
11 The Continuum Model (1D) Conservation of Mass Odd Finite first moment Cont. & piecewise diff. ( 0) Jump discontinuity ( = 0) {Crosses zero for eactly one Morse Social Force f s () (0, 1 F ) Piecewise Linear Social Force f s () (0, F ) ( L, 1) fs() fs() ( 1, 0) (1, 0) (0, 1 +F ) ( L, 1) (0, F )
12 A Menagerie of Behaviors 0.2 Spreading 0.4 Steady state !() 0.1!() !20! !4! (c) 1.5 (d) 1.5 Blow-up 1 Blow-up!() 1!() !2! !4!
13 Morse Social Force: A Phase Diagram F L
14 Cumulative Mass Function Define the cumulative mass function: Where satisfies a wave equation And the velocity can be approimated as Where. Note that as the density is positive,, is increasing.
15 Burger s Equation The dynamics are dominated by - If social forces are repulsive at small distances and density is bounded for all t. - If social forces are attractive at small distances and will form a shock. But as, Shocks = Density δ-functions!! See Bodnar & Velazquez (2006)
16 Spreading Swarms Epand velocity at long lengthscales: Where is the first moment Substitute into C of M,, to yield the porous media equation, Barenblatt Solution which has spreading solutions for κ< 0.
17 Barenblatt Solution Rescaled Density Approaches Barenblatt Solution RMS width has proper time dependence
18 Finding Steady States!() (b) Steady state When β > 0, κ< 0 we epect steady solutions !4! β α Energy first variation Fredholm Integral Equation Details Fredholm integral equation V ( z z )ρ(z ) dz = λ f(z), Coming Later Where α < z < β is the support of the solution Mass Constraint β α ρ(z ) dz = M
19 Phase Diagram Revisited Morse Social Force: F L
20 Some conclusions Repulsion (attraction) at the origin suppresses (causes) blow-up and density concentrations. First moment governs spreading at long lengthscales. Spreading governed by porous media equation. Steady states can be found via integral equations. Open Questions Necessary and sufficient conditions for spreading and convergence to Barenblatt solution.
21 A fly in the ointment... This potential has β > 0, κ> 0 which suggests spreading Social Force Density Plot RMS Width Final Steady State
22 Let s focus on locusts...
23 Locusts form giant, destructive swarms. National Geographic, A Perfect Swarm $10 billion/yr for pesticides (all insects) $100 million/yr for control (Africa) $100 million/yr in crop loss (Africa) (EPA, UNFAO) locusts km/day days or weeks 100 km 2
24 With gravity, swarm forms a bubble. Catastrophic Gravity G = 1 Sticky boundary
25 With wind, swarms roll Catastrophic Wind U = 1 Takeoff Resting Landing Uvarov, Grasshoppers & Locusts (1977)
26 Can we eplain the swarm morphology? Steady Catastrophic Locust Swarm 1.2 (b) z Empty Gap Mass concentration on ground
27 A Continuum Swarming Model Velocity Conservation of Mass q() Eternal Force (0, 1 F ) (0, 1 +F ) Morse Social Force Finite Mass Compact Support Pairwise, additive Antisymmetric Finite first moment Jump discontinuity ( = 0) Bodnar & Velasquez (2005,2006)
28 Energy for the Continuum Model Energy Eternal Potential Energy Dissipation Morse Potential Energy dissipated unless equilibrium Connection to optimal transport Pairwise, additive [ Q ()= - q(), F () = - f() ] Repulsive at small and attractive at large Symmetric Pointy at ( = 0)
29 First Variation and Steady Solutions Let where is a steady solution and is a zero mass perturbation. Then Where the first variation is given by: For the first variation to vanish, in the support of the solution: Z Q( y) (y) dy + F () = for 2. which is a Fredholm Integral Equation. Constant energy per unit mass
30 Second Variation and Stability Suppose Z is an equilibrium solution satisfying: Q( y) (y) dy + F () = for 2. Then for stability two conditions are sufficient: I) The energy per unit mass minimizes the induced potential, (): Z Q( y) (y) dy + F () = () II) Second Variation is positive definite: W 2 [, ] = Z for zero mass perturbations,. Z Q( y) () (y) d dy > 0 () () =
31 Finding Equilibrium States Equilibrium solutions satisfy an integral equation 0.4 Equilibrium (b) state 0.3 Energy first variation Fredholm integral equation!() !4! Fredholm Integral Equation Mass Constraint Where α < z < β is the support of the solution
32 Finding Equilibrium States!() Equilibrium (b) state 0.1 0!4! Energy first variation Fredholm integral equation invert integral operator ODE If then So applying to yields
33 Finding Equilibrium States!() Equilibrium (b) state 0.1 0!4! Energy first variation Fredholm integral equation invert integral operator ODE Linear ODE Solving Plug in and satisfy mass constraint Solution Yields ρ(z) = A + B cos (μz) in the support of the solution.
34 Continuum solutions agree with discrete (numerical) ones. 30 ρ vs. z Steady Solutions (N = 100) F = 0.5, L = ρ(z) = A + B cos (μz)
35 A Primer of Swarm Equilibria
36 Repulsive Forces on a Finite Interval Laplace Repulsive Potential: Finite Interval: -d < < d Equipartition of mass between concentrations at the endpoints and classical solution in interior. Mass Concentrations can t be found by classical integral equation methods Numerics (a) 0.3 ρ
37 Repulsive Forces in a Quadratic Potential Laplace Repulsive Potential: Quadratic eternal potential: Solution is compactly supported on. Jump discontinuities at edge of support Global energy minimizer. 1.5 Numerics (c) 1.4 ρ = = =
38 Repulsive forces in a gravitational potential Laplace Repulsive Potential: Half line: If M < g all mass concentrates at origin. If M > g some mass concentrates at origin and the remainder forms a compact swarm. Global minimizer and probably global attractor. 0.6 Numerics (b) 0.4 ρ As a swarm model: - Concentration on ground (good) but no gap (bad).
39 Catastrophic Morse Potential in Free Space Morse Potential in catastrophic regime. Solution is compactly supported on. Width of support independent of mass!! Jump discontinuities at edge of support. Local energy minimizer. Probably a global energy minimizer. Numerics (a) Morse Potential Analytic Solution 0.3 ρ
40 H-Stable Morse Potential on bounded interval Morse Potential in H-stable regime. No solution in free space (swarm spreads forever). Mass concentrations at endpoints of interval. Global energy minimizer. Probably global attractor. Morse Potential Numerics (b) ρ 0.26 = = =
41 Quasi 2D-Potential Suppose that we look at a quasi-2d mass distribution This has an associated potential, : Q 2D (z) = Z 1 = 1 Q( p 2 + z 2 ) d = Z 1 e p 2 +z 2 = 1 This allows steady swarms with a gap at zero: d Airborne swarm z Repulsion from mass concentrated on ground ~ dz 2 ln(z) Gravity + gz = Stable gap (for some parameters) Empty Gap Mass concentration on ground
42 The quasi 2D potential yields dynamically stable bubbles. Organism column Density profile ρ 1D z ρ Quasi 2D 35 insects z
43 Quasi 2D Locust Swarms Quasi 2D Laplace Repulsion Potential. Gravity & Ground: - Linear potential on a half space. Mass concentration on ground Gap above ground. Continuous family of solutions Solution is only a local swarm minimizer. Quasi-2D Numerics (d) 8 6 ρ() (b) (c) 12.8 z 0.8 Airborne swarm Λ() Empty Gap Mass concentration on ground
44 2D Morse Equilibria 2D Catastrophic Swarm (d) 0.6 Why are there mass concentrations at the edges? 8 6 2D Morse Ribbons (d) 0 ρ() 4 2!0.6! Mass concentration on ground z Airborne swarm Louis Ryan HMC 2012 Empty Gap
45 2D Morse - Ribbon Morse Potential Quasi-2D Morse Potential Q 2D (z) = z 1 Quasi-2D Morse Potential ( ) Q 2D (z) Z L L Z 1 = 1 z 2 ln z Q( p 2 + z 2 ) d Integral Equation (Carleman) y 2 ln y (z y) dz = z <L Density z L (z) -L Numerical Quasi 2D Swarm Carleman (z) = C L 2 z 2 z <L Position
46 2D Morse - Circular Swarm Density [ ρ(r) ] density D Catastrophic Swarm (Continuum Model) Plot of Swarm Density for C = e+000, L = 2, n = radius Radius [r] Found as an Energy Minimizer Singular mass concentration at the edges (r) C p R r 2D Catastrophic Swarm (Particle Simulation) 0.6 0!0.6 (d)! Dimensionality of local minimizers of the interaction energy D Balagué, JA Carrillo, T Laurent, G Raoul Archive for Rational Mechanics and Analysis (2013)
47 Conclusions Energy methods allow us to analyze equlibria and their stability. Equilibria satisfy a Fredholm Integral Equation. This yields some analytical solutions. Equlibria often contain concentrations on boundaries. For locust model, quasi-2d repulsive potential and gravity yield a ground concentration together with a gap. 2D Morse yields inverse square-root mass singularities on boundaries. Inverse square-root singularities appear generic in 2D for pointy repulsive potentials. z 1D Quasi 2D Empty Gap Mass concentration on ground Bernoff & Topaz Nonlinear Aggregation Equations: A Primer of Swarm Equilibria. SIREV (2013).
An Integrodifferential Equation Modeling 1-D Swarming Behavior
An Integrodifferential Equation Modeling 1-D Swarming Behavior Andrew Leverentz Andrew Bernoff, Advisor Chad Topaz, Reader May, 28 Department of Mathematics Copyright c 28 Andrew Leverentz. The author
More informationAsymptotic Dynamics of Attractive-Repulsive Swarms
Claremont Colleges Scholarship @ Claremont All HMC Faculty Publications and Research HMC Faculty Scholarship 7-15-2009 Asymptotic Dynamics of Attractive-Repulsive Swarms Andrew J. Leverentz '08 Harvey
More informationSynchronization and Swarming: Clocks and Flocks
Synchronization and Swarming: Clocks and Flocks Andrew J. Bernoff, Harvey Mudd College Thanks to Chad Topaz, Macalester College Andrew J. Bernoff Harvey Mudd College ajb@hmc.edu Chad M. Topaz Macalester
More informationSwarming by Nature and by Design
Swarming by Nature and by Design Andrea Bertozzi Sonia Kovalevsky Lecture Department of Mathematics, UCLA Thanks to support from ARO, NSF, and ONR and to contributions from numerous collaborators. Students,
More informationROBERT SIMIONE, DEJAN SLEPČEV, AND IHSAN TOPALOGLU
EXISTENCE OF MINIMIZERS OF NONLOCAL INTERACTION ENERGIES ROBERT SIMIONE, DEJAN SLEPČEV, AND IHSAN TOPALOGLU Abstract. We investigate nonlocal-interaction energies on the space of probability measures.
More informationA model for rolling swarms of locusts
Macalester College From the SelectedWorks of Chad M. Topaz 2008 A model for rolling swarms of locusts Chad M. Topaz, Macalester College Andrew J. Bernoff Sheldon Logan, Harvey Mudd College Wyatt Toolson,
More informationFinite-time blow-up of solutions of an aggregation equation in R n
Finite-time blow-up of solutions of an aggregation equation in R n Andrea L. Bertozzi a,b and Thomas Laurent b a Department of Mathematics, University of California, Los Angeles 90095 b Department of Mathematics,
More informationMultiscale Analysis of Many Particle Systems with Dynamical Control
Michael Herrmann Multiscale Analysis of Many Particle Systems with Dynamical Control joint work with Barbara Niethammer and Juan J.L. Velázquez Kinetic description of multiscale phenomena Archimedes Center
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 informationSample Final Questions: Solutions Math 21B, Winter y ( y 1)(1 + y)) = A y + B
Sample Final Questions: Solutions Math 2B, Winter 23. Evaluate the following integrals: tan a) y y dy; b) x dx; c) 3 x 2 + x dx. a) We use partial fractions: y y 3 = y y ) + y)) = A y + B y + C y +. Putting
More informationThe Global Stability of the Solution to the Morse Potential in a Catastrophic Regime
Claremont Colleges Scholarship @ Claremont HMC Senior Theses HMC Student Scholarship 2016 The Global Stability of the Solution to the Morse Potential in a Catastrophic Regime Weerapat Pittayakanchit Harvey
More informationPropagation and quenching in a reactive Burgers-Boussinesq system.
Propagation and quenching in a reactive system., University of New Meico, Los Alamos National Laboratory, University of Chicago October 28 Between Math and Astrophysics Motivation The model Reaction Compression
More informationNotes for Introduction to Partial Di erential Equations, Fall 2018
Notes for Introduction to Partial Di erential Equations, Fall 2018 Esteban G. Tabak 1 Tra c Flow We eplore a simple mathematical model for tra c flow. This will be our introduction to systems of conservation
More informationEINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science. CASA-Report June 2014
EINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science CASA-Report 14-20 June 2014 Anisotropic interactions in a first-order aggregation model: a proof of concept by J.H.M. Evers,
More informationCrime Modeling with Lévy Flights
Crime Modeling with Lévy Flights Jonah Breslau, 1 Sorathan (Tum) Chaturapruek, 2 Daniel Yazdi 3 Mentors: Professor Theodore Kolokolnikov, 4 Professor Scott McCalla 3 1 Pomona College 2 Harvey Mudd College
More informationPiecewise Smooth Solutions to the Burgers-Hilbert Equation
Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang
More informationMath 256: Applied Differential Equations: Final Review
Math 256: Applied Differential Equations: Final Review Chapter 1: Introduction, Sec 1.1, 1.2, 1.3 (a) Differential Equation, Mathematical Model (b) Direction (Slope) Field, Equilibrium Solution (c) Rate
More informationy x is symmetric with respect to which of the following?
AP Calculus Summer Assignment Name: Note: Unless otherwise specified, the domain of a function f is assumed to be the set of all real numbers for which f () is a real number. Part : Multiple Choice Solve
More informationCalculus Math 21B, Winter 2009 Final Exam: Solutions
Calculus Math B, Winter 9 Final Exam: Solutions. (a) Express the area of the region enclosed between the x-axis and the curve y = x 4 x for x as a definite integral. (b) Find the area by evaluating the
More informationA high order adaptive finite element method for solving nonlinear hyperbolic conservation laws
A high order adaptive finite element method for solving nonlinear hyperbolic conservation laws Zhengfu Xu, Jinchao Xu and Chi-Wang Shu 0th April 010 Abstract In this note, we apply the h-adaptive streamline
More informationPredator-swarm interactions
Predator-swarm interactions Hypothesis: swarming behaviour is an evolutionary adaptation that confers certain benefits on the individuals or group as a whole [Parrish,Edelstein-Keshet 1999; Sumpter 2010,
More informationEquilibria for a two-dimensional aggregation model with two species
Equilibria for a two-dimensional aggregation model with two species Joep H.M. Evers Razvan C. Fetecau Theodore Kolokolnikov August 1, 17 Abstract We consider an aggregation model for two interacting species.
More information3 Hydrostatic Equilibrium
3 Hydrostatic Equilibrium Reading: Shu, ch 5, ch 8 31 Timescales and Quasi-Hydrostatic Equilibrium Consider a gas obeying the Euler equations: Dρ Dt = ρ u, D u Dt = g 1 ρ P, Dɛ Dt = P ρ u + Γ Λ ρ Suppose
More informationChapter 3 Burgers Equation
Chapter 3 Burgers Equation One of the major challenges in the field of comple systems is a thorough understanding of the phenomenon of turbulence. Direct numerical simulations (DNS) have substantially
More informationCHAPTER 1 Introduction to Differential Equations 1 CHAPTER 2 First-Order Equations 29
Contents PREFACE xiii CHAPTER 1 Introduction to Differential Equations 1 1.1 Introduction to Differential Equations: Vocabulary... 2 Exercises 1.1 10 1.2 A Graphical Approach to Solutions: Slope Fields
More informationChapter 1. Introduction to Nonlinear Space Plasma Physics
Chapter 1. Introduction to Nonlinear Space Plasma Physics The goal of this course, Nonlinear Space Plasma Physics, is to explore the formation, evolution, propagation, and characteristics of the large
More informationIndex. C 2 ( ), 447 C k [a,b], 37 C0 ( ), 618 ( ), 447 CD 2 CN 2
Index advection equation, 29 in three dimensions, 446 advection-diffusion equation, 31 aluminum, 200 angle between two vectors, 58 area integral, 439 automatic step control, 119 back substitution, 604
More informationCalculus I Practice Test Problems for Chapter 2 Page 1 of 7
Calculus I Practice Test Problems for Chapter Page of 7 This is a set of practice test problems for Chapter This is in no way an inclusive set of problems there can be other types of problems on the actual
More informationElliptical-like orbits on a warped spandex fabric
Elliptical-like orbits on a warped spandex fabric Prof. Chad A. Middleton CMU Physics Seminar September 3, 2015 Submitted to the American Journal of Physics Kepler s 3 Laws of planetary motion Kepler s
More informationDyson series for the PDEs arising in Mathematical Finance I
for the PDEs arising in Mathematical Finance I 1 1 Penn State University Mathematical Finance and Probability Seminar, Rutgers, April 12, 2011 www.math.psu.edu/nistor/ This work was supported in part by
More informationQuantum Evaporation of 2-d Black Holes
p. Quantum Evaporation of 2-d Black Holes Abhay Ashtekar Institute for Gravitation and the Cosmos, Penn State CGHS Black holes: A conceptually new angle and high precision numerics. This talk will put
More informationProblems in diffusion and absorption: How fast can you hit a target with a random walk?
Problems in diffusion and absorption: How fast can you hit a target with a random walk? Andrew J. Bernoff Harvey Mudd College In collaboration with Alan Lindsay (Notre Dame) Thanks to Alan Lindsay, Michael
More informationSpotlight on Laplace s Equation
16 Spotlight on Laplace s Equation Reference: Sections 1.1,1.2, and 1.5. Laplace s equation is the undriven, linear, second-order PDE 2 u = (1) We defined diffusivity on page 587. where 2 is the Laplacian
More informationP = 1 3 (σ xx + σ yy + σ zz ) = F A. It is created by the bombardment of the surface by molecules of fluid.
CEE 3310 Thermodynamic Properties, Aug. 27, 2010 11 1.4 Review A fluid is a substance that can not support a shear stress. Liquids differ from gasses in that liquids that do not completely fill a container
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 information1. Find A and B so that f x Axe Bx. has a local minimum of 6 when. x 2.
. Find A and B so that f Ae B has a local minimum of 6 when.. The graph below is the graph of f, the derivative of f; The domain of the derivative is 5 6. Note there is a cusp when =, a horizontal tangent
More informationResearch interests: Nonlinear partial differential equations and their applications
Research interests: Nonlinear partial differential equations and their applications Degenerate diffusion, higher order equations, reaction diffusion equations, heat flow of harmonic maps, pseudo-parabolic
More informationForaging Fruit Flies: Lagrangian and Eulerian Descriptions of Insect Swarming
Foraging Fruit Flies: Lagrangian and Eulerian Descriptions of Insect Swarming Joseph D. Majkut Andrew Bernoff, Advisor Chad Topaz, UCLA, Reader May, 26 Department of Mathematics Copyright c 26 Joseph D.
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 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 informationStep Bunching in Epitaxial Growth with Elasticity Effects
Step Bunching in Epitaxial Growth with Elasticity Effects Tao Luo Department of Mathematics The Hong Kong University of Science and Technology joint work with Yang Xiang, Aaron Yip 05 Jan 2017 Tao Luo
More informationCalculus - Chapter 2 Solutions
Calculus - Chapter Solutions. a. See graph at right. b. The velocity is decreasing over the entire interval. It is changing fastest at the beginning and slowest at the end. c. A = (95 + 85)(5) = 450 feet
More informationLocust Dynamics: Behavioral Phase Change and Swarming
Locust Dynamics: Behavioral Phase Change and Swarming Chad M. Topaz 1 *, Maria R. D Orsogna 2, Leah Edelstein-Keshet 3, Andrew J. Bernoff 4 1 Department of Mathematics, Statistics, and Computer Science,
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 informationPopulation Games and Evolutionary Dynamics
Population Games and Evolutionary Dynamics William H. Sandholm The MIT Press Cambridge, Massachusetts London, England in Brief Series Foreword Preface xvii xix 1 Introduction 1 1 Population Games 2 Population
More informationEffective dynamics of many-particle systems with dynamical constraint
Michael Herrmann Effective dnamics of man-particle sstems with dnamical constraint joint work with Barbara Niethammer and Juan J.L. Velázquez Workshop From particle sstems to differential equations WIAS
More informationProbing Holographic Superfluids with Solitons
Probing Holographic Superfluids with Solitons Sean Nowling Nordita GGI Workshop on AdS4/CFT3 and the Holographic States of Matter work in collaboration with: V. Keränen, E. Keski-Vakkuri, and K.P. Yogendran
More informationSolutions of Burgers Equation
ISSN 749-3889 (print, 749-3897 (online International Journal of Nonlinear Science Vol.9( No.3,pp.9-95 Solutions of Burgers Equation Ch. Srinivasa Rao, Engu Satyanarayana Department of Mathematics, Indian
More informationChapter 6 Laminar External Flow
Chapter 6 aminar Eternal Flow Contents 1 Thermal Boundary ayer 1 2 Scale analysis 2 2.1 Case 1: δ t > δ (Thermal B.. is larger than the velocity B..) 3 2.2 Case 2: δ t < δ (Thermal B.. is smaller than
More informationThreshold of singularity formation in the semilinear wave equation
PHYSICAL REVIEW D 71, 044019 (2005) Threshold of singularity formation in the semilinear wave equation Steven L. Liebling Department of Physics, Long Island University-C.W. Post Campus, Brookville, New
More informationRates of Convergence to Self-Similar Solutions of Burgers Equation
Rates of Convergence to Self-Similar Solutions of Burgers Equation by Joel Miller Andrew Bernoff, Advisor Advisor: Committee Member: May 2 Department of Mathematics Abstract Rates of Convergence to Self-Similar
More informationPhase-field systems with nonlinear coupling and dynamic boundary conditions
1 / 46 Phase-field systems with nonlinear coupling and dynamic boundary conditions Cecilia Cavaterra Dipartimento di Matematica F. Enriques Università degli Studi di Milano cecilia.cavaterra@unimi.it VIII
More informationPattern formation is one of the most studied aspects of animal
Comple spatial group patterns result from different animal communication mechanisms R. Eftimie, G. de Vries, and M. A. Lewis Center for Mathematical Biology, Department of Mathematical and Statistical
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 informationAim: Understand equilibrium of galaxies
8. Galactic Dynamics Aim: Understand equilibrium of galaxies 1. What are the dominant forces? 2. Can we define some kind of equilibrium? 3. What are the relevant timescales? 4. Do galaxies evolve along
More informationLyapunov stability ORDINARY DIFFERENTIAL EQUATIONS
Lyapunov stability ORDINARY DIFFERENTIAL EQUATIONS An ordinary differential equation is a mathematical model of a continuous state continuous time system: X = < n state space f: < n! < n vector field (assigns
More informationTraffic flow on networks: conservation laws models. Benedetto Piccoli I.A.C. C.N.R. Rome
Traffic flow on networks: conservation laws models Benedetto Piccoli I.A.C. C.N.R. Rome Conservation laws on networks u t + f(u) x=0 Dynamics at nodes? 1. The only conservation at nodes does not determine
More informationLie Symmetry Analysis and Approximate Solutions for Non-linear Radial Oscillations of an Incompressible Mooney Rivlin Cylindrical Tube
Ž Journal of Mathematical Analysis and Applications 45, 34639 doi:116jmaa6748, available online at http:wwwidealibrarycom on Lie Symmetry Analysis and Approimate Solutions for Non-linear Radial Oscillations
More informationAsymptotic and Perturbation Methods: Lecture Notes.
Asymptotic and Perturbation Methods: Lecture Notes. Blerta Shtylla University of Utah Mathematics Department Fall 2008 1 Lecture 1 and Lecture 2 Week of August 26, 2008 Introduction: For many problems
More informationMULTIVARIABLE INTEGRATION
MULTIVARIABLE INTEGRATION (SPHERICAL POLAR COORDINATES) Question 1 a) Determine with the aid of a diagram an expression for the volume element in r, θ, ϕ. spherical polar coordinates, ( ) [You may not
More informationarxiv: v1 [math.ap] 8 Feb 2018
The dynamics of disappearing pulses in a singularly perturbed reaction-diffusion system with parameters that vary in time and space Robbin Bastiaansen and Arjen Doelman July 7, 18 arxiv:18.737v1 [math.ap]
More informationAncient solutions to Geometric Flows Lecture No 2
Ancient solutions to Geometric Flows Lecture No 2 Panagiota Daskalopoulos Columbia University Frontiers of Mathematics and Applications IV UIMP 2015 July 20-24, 2015 Topics to be discussed In this lecture
More information1993 AP Calculus AB: Section I
99 AP Calculus AB: Section I 9 Minutes Scientific Calculator Notes: () The eact numerical value of the correct answer does not always appear among the choices given. When this happens, select from among
More informationF n = F n 1 + F n 2. F(z) = n= z z 2. (A.3)
Appendix A MATTERS OF TECHNIQUE A.1 Transform Methods Laplace transforms for continuum systems Generating function for discrete systems This method is demonstrated for the Fibonacci sequence F n = F n
More informationThe Multi-domain Hybrid Method for Calculating the Power Law Decay of Gravitational Waveforms with the Analytic Radiation Boundary Conditions
The Multi-domain Hybrid Method for Calculating the Power Law Decay of Gravitational Waveforms with the Analytic Radiation Boundary Conditions Joshua Buli (UCR) in collaboration with Dr. Jae-Hun Jung (SUNY
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 informationIntegrodifference Models with Temporally Varying Environments
Integrodifference Models with Temporally Varying Environments Jon Jacobsen Harvey Mudd College Claremont, CA 91711 Integrodifference Models with Temporally Varying Environments Jon Jacobsen Harvey Mudd
More informationHamiltonian partial differential equations and Painlevé transcendents
The 6th TIMS-OCAMI-WASEDA Joint International Workshop on Integrable Systems and Mathematical Physics March 22-26, 2014 Hamiltonian partial differential equations and Painlevé transcendents Boris DUBROVIN
More informationNumerical simulation of nonlinear continuity equations by evolving diffeomorphisms
www.oeaw.ac.at Numerical simulation of nonlinear continuity equations by evolving diffeomorphisms J.A. Carrillo, H. Ranetbauer, M-T. Wolfram RICAM-Report 26-9 www.ricam.oeaw.ac.at NUMERICAL SIMULATION
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 informationNonlinear Integral Equation Formulation of Orthogonal Polynomials
Nonlinear Integral Equation Formulation of Orthogonal Polynomials Eli Ben-Naim Theory Division, Los Alamos National Laboratory with: Carl Bender (Washington University, St. Louis) C.M. Bender and E. Ben-Naim,
More informationSwarm shape and its dynamics in a predator-swarm model
Swarm shape and its dynamics in a predator-swarm model Hayley Tomkins Advised by Theodore Kolokolnikov June 5, 015 Abstract We consider a particle predator-swarm model introduced in [1]. In the continuum
More informationUpscaling Wave Computations
Upscaling Wave Computations Xin Wang 2010 TRIP Annual Meeting 2 Outline 1 Motivation of Numerical Upscaling 2 Overview of Upscaling Methods 3 Future Plan 3 Wave Equations scalar variable density acoustic
More informationSelf-Similar Gas Dynamics of Voids and Supernova Ejecta
Self-Similar Gas Dynamics of Voids and Supernova Ejecta arxiv:1109.2682v2 [astro-ph.sr] Li-Le Wang Director: Prof. Yu-Qing Lou Department of Physics and Tsinghua Center for Astrophysics (THCA) Tsinghua
More informationStudy Guide and Intervention. The Quadratic Formula and the Discriminant. Quadratic Formula. Replace a with 1, b with -5, and c with -14.
Study Guide and Intervention Quadratic Formula The Quadratic Formula can be used to solve any quadratic equation once it is written in the form a 2 + b + c = 0. Quadratic Formula The solutions of a 2 +
More information1.2. Introduction to Modeling. P (t) = r P (t) (b) When r > 0 this is the exponential growth equation.
G. NAGY ODE January 9, 2018 1 1.2. Introduction to Modeling Section Objective(s): Review of Exponential Growth. The Logistic Population Model. Competing Species Model. Overview of Mathematical Models.
More informationMath 75B Practice Problems for Midterm II Solutions Ch. 16, 17, 12 (E), , 2.8 (S)
Math 75B Practice Problems for Midterm II Solutions Ch. 6, 7, 2 (E),.-.5, 2.8 (S) DISCLAIMER. This collection of practice problems is not guaranteed to be identical, in length or content, to the actual
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 informationUBC-SFU-UVic-UNBC Calculus Exam Solutions 7 June 2007
This eamination has 15 pages including this cover. UBC-SFU-UVic-UNBC Calculus Eam Solutions 7 June 007 Name: School: Signature: Candidate Number: Rules and Instructions 1. Show all your work! Full marks
More informationCS 542G: The Poisson Problem, Finite Differences
CS 542G: The Poisson Problem, Finite Differences Robert Bridson November 10, 2008 1 The Poisson Problem At the end last time, we noticed that the gravitational potential has a zero Laplacian except at
More informationA mechanism for the invasion of a sub-tropical high by a hurricane: Part 1, static analysis
A mechanism for the invasion of a sub-tropical high by a hurricane: Part 1, static analysis Les Hatton Computing Laboratory, University of Kent October 5, 003 Abstract The major point of this paper is
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 informationMath 111 Final Exam Review KEY
Math 111 Final Eam Review KEY 1. Use the graph of y = f in Figure 1 to answer the following. Approimate where necessary. a Evaluate f 1. f 1 = 0 b Evaluate f0. f0 = 6 c Solve f = 0. =, = 1, =,or = 3 Solution
More informationBoost-invariant dynamics near and far from equilibrium physics and AdS/CFT.
Boost-invariant dynamics near and far from equilibrium physics and AdS/CFT. Micha l P. Heller michal.heller@uj.edu.pl Department of Theory of Complex Systems Institute of Physics, Jagiellonian University
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 informationPlanet Formation. XIII Ciclo de Cursos Especiais
Planet Formation Outline 1. Observations of planetary systems 2. Protoplanetary disks 3. Formation of planetesimals (km-scale bodies) 4. Formation of terrestrial and giant planets 5. Evolution and stability
More informationSolve Wave Equation from Scratch [2013 HSSP]
1 Solve Wave Equation from Scratch [2013 HSSP] Yuqi Zhu MIT Department of Physics, 77 Massachusetts Ave., Cambridge, MA 02139 (Dated: August 18, 2013) I. COURSE INFO Topics Date 07/07 Comple number, Cauchy-Riemann
More informationCausality, hyperbolicity, and shock formation in Lovelock theories
Causality, hyperbolicity, and shock formation in Lovelock theories Harvey Reall DAMTP, Cambridge University HSR, N. Tanahashi and B. Way, arxiv:1406.3379, 1409.3874 G. Papallo, HSR arxiv:1508.05303 Lovelock
More informationTHE CAHN-HILLIARD EQUATION WITH A LOGARITHMIC POTENTIAL AND DYNAMIC BOUNDARY CONDITIONS
THE CAHN-HILLIARD EQUATION WITH A LOGARITHMIC POTENTIAL AND DYNAMIC BOUNDARY CONDITIONS Alain Miranville Université de Poitiers, France Collaborators : L. Cherfils, G. Gilardi, G.R. Goldstein, G. Schimperna,
More informationSection 7.4: Inverse Laplace Transform
Section 74: Inverse Laplace Transform A natural question to ask about any function is whether it has an inverse function We now ask this question about the Laplace transform: given a function F (s), will
More informationF. MENZAQUE, R. R. ROSALES, E. G. TABAK, AND C. V. TURNER and dissipation are many decades apart, the intermediate scales span the so{called turbulent
The Forced Inviscid Burgers Equation as a Model for Nonlinear Interactions among Dispersive Waves Fernando Menzaque, Rodolfo R. Rosales, Esteban G. Tabak, and Cristina V. Turner Abstract. The forced inviscid
More informationMULTIVARIABLE CALCULUS
MULTIVARIABLE CALCULUS Summer Assignment Welcome to Multivariable Calculus, Multivariable Calculus is a course commonly taken by second and third year college students. The general concept is to take the
More informationChapter 1 Direct Modeling for Computational Fluid Dynamics
Chapter 1 Direct Modeling for Computational Fluid Dynamics Computational fluid dynamics (CFD) is a scientific discipline, which aims to capture fluid motion in a discretized space. The description of the
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationBlow up of solutions for a 1D transport equation with nonlocal velocity and supercritical dissipation
Blow up of solutions for a 1D transport equation with nonlocal velocity and supercritical dissipation Dong Li a,1 a School of Mathematics, Institute for Advanced Study, Einstein Drive, Princeton, NJ 854,
More informationModeling group formation and activity. patterns in self-organizing collectives of. individuals
Modeling group formation and activity patterns in self-organizing collectives of individuals R. Eftimie a,, G. de Vries a, M. A. Lewis a,b, F. Lutscher a,c a Centre for Mathematical Biology, Department
More information( ) ( ) = (2) u u u v v v w w w x y z x y z x y z. Exercise [17.09]
Eercise [17.09] Suppose is a region of (Newtonian) space, defined at a point in time t 0. Let T(t,r) be the position at time t of a test particle (of insignificantly tiny mass) that begins at rest at point
More informationSolutions in the sense of distributions. Solutions in the sense of distributions Definition, non uniqueness
Solutions in the sense of distributions Definition, non uniqueness 1. Notion of distributions In order to build weak solutions to the Hopf equation, we need to define derivatives of non smooth functions,
More informationMath 231 Final Exam Review
Math Final Eam Review Find the equation of the line tangent to the curve 4y y at the point (, ) Find the slope of the normal line to y ) ( e at the point (,) dy Find d if cos( y) y 4 y 4 Find the eact
More informationWest Essex Regional School District. AP Calculus AB. Summer Packet
West Esse Regional School District AP Calculus AB Summer Packet 05-06 Calculus AB Calculus AB covers the equivalent of a one semester college calculus course. Our focus will be on differential and integral
More information