Crime Modeling with Lévy Flights
|
|
- Lillian Banks
- 6 years ago
- Views:
Transcription
1 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 3 University of California, Los Angeles 4 Dalhousie University 8/8/2012
2 Figure: Crime Hot-Spot Pattern, Long Beach, LA. Short et al [2]
3 Hotspot Modeling Our goal is to model crime, specifically the hotspot phenomenon.
4 Hotspot Modeling Our goal is to model crime, specifically the hotspot phenomenon. Crime hotspots are when several crimes occur in a short period of time and in a small area.
5 Hotspot Modeling Our goal is to model crime, specifically the hotspot phenomenon. Crime hotspots are when several crimes occur in a short period of time and in a small area. Considerable empirical evidence behind them.
6 Hotspot Modeling Our goal is to model crime, specifically the hotspot phenomenon. Crime hotspots are when several crimes occur in a short period of time and in a small area. Considerable empirical evidence behind them. We focus on burglaries for simplicity.
7 Hotspot Modeling: Theory What causes hotspots?
8 Hotspot Modeling: Theory What causes hotspots? Repeat/ Near Repeat Effects: increased knowledge of location after successful crime. Figure: Repeat Effects. (Short et al. 2009, [1])
9 Hotspot Modeling: Theory What causes hotspots? Repeat/ Near Repeat Effects: increased knowledge of location after successful crime. Figure: Repeat Effects. (Short et al. 2009, [1]) Broken Windows Theory: crime causes sense of lawlessness.
10 The Short et al Model (2008) Starts with discrete model
11 The Short et al Model (2008) Starts with discrete model Criminals are agents on a lattice that has attractiveness field
12 The Short et al Model (2008) Starts with discrete model Criminals are agents on a lattice that has attractiveness field Every time step, criminals move to neighboring lattice spaces based on attractiveness.
13 The Short et al Model (2008) Starts with discrete model Criminals are agents on a lattice that has attractiveness field Every time step, criminals move to neighboring lattice spaces based on attractiveness. Then decide to burgle or not based off attractiveness.
14 The Short et al Model (2008) Starts with discrete model Criminals are agents on a lattice that has attractiveness field Every time step, criminals move to neighboring lattice spaces based on attractiveness. Then decide to burgle or not based off attractiveness. Crimes are self-exciting; cause attractiveness at lattice point and neighbors to increase.
15 The Short et al Model (2008) (cont.) Then, taking limits as grid spacing and time steps go to zero, we get system of PDEs The Model ρ t = A t = η A A + ρa + A 0 ( ( ρ A A) ρ ) A A ρa + A A 0 (1a) (1b) A = attractiveness at a point ρ = criminal density at a location A 0 = a constant, background level of attractiveness A = the spatially homogeneous equilibrium solution for A.
16 Lévy Flights Short et al model assumes Brownian motion (and thus distance traveled normally distributed)
17 Lévy Flights Short et al model assumes Brownian motion (and thus distance traveled normally distributed) We use Lévy Flights
18 Lévy Flights Short et al model assumes Brownian motion (and thus distance traveled normally distributed) We use Lévy Flights Power law distribution of step sizes (P(k) k (2s+1), 0 < s 1)
19 Lévy Flights Short et al model assumes Brownian motion (and thus distance traveled normally distributed) We use Lévy Flights Power law distribution of step sizes (P(k) k (2s+1), 0 < s 1) Fractal-like motion; reflects many scales of human movement.
20 Implementation: Discrete We change the transition probability as follows: p i i+1 (t) = A i+1 (t) A i+1 (t) + A i 1 (t) (Brownian) p i j (t) = A j (t) i j (2s+1) Σ k Z,k i A k (t) i k (2s+1) (Lévy)
21 Implementation: Continuous From before we had: A t = η A A + ρa + A 0, ( ( ρ ρ t = A A) ρ ) A A ρa + A A 0.
22 Implementation: Continuous From before we had: A t = η A A + ρa + A 0, ( ( ρ ρ t = A A) ρ ) A A ρa + A A 0. Using our new transition probabilities and taking limits as before we get: A t = η A A + ρa + A 0 (No change) ( ( ρ t = A s ρ ) ρ ) A A s A ρa + A A 0,
23 Implementation: Continuous From before we had: A t = η A A + ρa + A 0, ( ( ρ ρ t = A A) ρ ) A A ρa + A A 0. Using our new transition probabilities and taking limits as before we get: where A t = η A A + ρa + A 0 (No change) ( ( ρ t = A s ρ ) ρ ) A A s A ρa + A A 0, s A = c s c s is a constant and 0 < s 1. A(y) A(x) dy, (2) y x 2s+1
24 Fractional Calculus The operator s is the Riesz Derivative,
25 Fractional Calculus The operator s is the Riesz Derivative, Generalization of the second derivative such that lim s 1 s A = 2 A x. 2
26 Fractional Calculus The operator s is the Riesz Derivative, Generalization of the second derivative such that lim s 1 s A = 2 A x. 2 But note that 1 2 A A x
27 Fractional Calculus The operator s is the Riesz Derivative, Generalization of the second derivative such that lim s 1 s A = 2 A x. 2 But note that 1 2 A A x Is a non-local operator. Leads to super-diffusion. Degree of non-locality controlled by s.
28 Fractional Calculus The operator s is the Riesz Derivative, Generalization of the second derivative such that lim s 1 s A = 2 A x. 2 But note that 1 2 A A x Is a non-local operator. Leads to super-diffusion. Degree of non-locality controlled by s. Fourier Transform has nice property: F x q { s A} = q 2s Â
29 Numerical Solutions: Spectral Method of Lines Discretize in space and the calculate derivatives in Fourier space. This turns into ODE in time.
30 Numerical Solutions: Spectral Method of Lines Discretize in space and the calculate derivatives in Fourier space. This turns into ODE in time. We make use of the fact that F x q { s A} = q 2s Â.
31 Numerical Solutions: Spectral Method of Lines Discretize in space and the calculate derivatives in Fourier space. This turns into ODE in time. We make use of the fact that F x q { s A} = q 2s Â. Then use MATLAB s stiff ODE solver.
32 Numerical Solutions: Spectral Method of Lines Discretize in space and the calculate derivatives in Fourier space. This turns into ODE in time. We make use of the fact that F x q { s A} = q 2s Â. Then use MATLAB s stiff ODE solver. We used to use a forward Euler spectral method, but that could not handle much of the parameter space.
33 Examples Figure: Single Hot-Spot.
34 Examples Figure: Four Hot-Spots, D = 1, s = 1, ε = 0.05.
35 Examples Figure: No Hot-Spot, D = 1, s = 0.5, ε = 0.05.
36 Examples Figure: Oscillating Hot-Spots, s = 0.7, η = 0.1, ρ = 0.4, A = 0.12.
37 Linear Stability of Homogeneous Equilibrium The homogeneous equilibrium does not change from Short et al. We get A = A 0 + B and ρ = B A 0 + B. (3)
38 Linear Stability of Homogeneous Equilibrium The homogeneous equilibrium does not change from Short et al. We get A = A 0 + B and ρ = B A 0 + B. (3) Perturb from the homogeneous equilibrium as follows and plug into linearized equations: A(x, t) = A + δ A e σt e ik x, (4a) ρ(x, t) = ρ + δ ρ e σt e ik x. (4b)
39 Linear Stability of Homogeneous Equilibrium The homogeneous equilibrium does not change from Short et al. We get A = A 0 + B and ρ = B A 0 + B. (3) Perturb from the homogeneous equilibrium as follows and plug into linearized equations: A(x, t) = A + δ A e σt e ik x, (4a) ρ(x, t) = ρ + δ ρ e σt e ik x. (4b) Solving the resulting eigenvalue problem results in this condition for instability: there exists k such that η k 2s+2 k 2s (3ρ 1) + ηa k 2 + A < 0. (5)
40 Stability (cont.) The first attempt to solve Eq. (5): the condition for instability is equivalent to, there exists k such that ρ > 1 ( 1 + η k 2 ) ( 1 + A ) 3 k 2s. (6)
41 Stability (cont.) The first attempt to solve Eq. (5): the condition for instability is equivalent to, 1 ( ρ > inf 1 + η k 2 ) ( 1 + A ) k 3 k 2s. (6)
42 Stability (cont.) The first attempt to solve Eq. (5): the condition for instability is equivalent to, ρ > where k is a root of 1 ( 1 + η k 2) ( 1 + A ) 3 k 2s, (6) k 2s+2 + A(1 s) k 2 As η = 0. (7)
43 Stability (cont.) The second attempt to solve Eq. (5): the condition for instability is equivalent to, there exists k such that A < 3ρ k 2s 1 + η k 2 k 2s. (8)
44 Stability (cont.) The second attempt to solve Eq. (5): the condition for instability is equivalent to, A < sup k 3ρ k 2s 1 + η k 2 k 2s. (8)
45 Stability (cont.) The second attempt to solve Eq. (5): the condition for instability is equivalent to, A < 3ρ k 2s 1 + η k 2 k 2s, (8) where k is a root of the equation η 2 s k 4 + η(3ρ(1 s) + 2s) k 2 + s(1 3ρ) = 0. (9)
46 Stability (cont.) The second attempt to solve Eq. (5): the condition for instability is equivalent to, A < 3ρ k 2s 1 + η k 2 k 2s, (8) where k is a root of the equation η 2 s k 4 + η(3ρ(1 s) + 2s) k 2 + s(1 3ρ) = 0. (9)
47 Stability (cont.) For ρ > 1 3, this generates condition for linear instability for the system: ( 3ρ(1 s) 2s + ) s ( W 3ρ(1 + s) W A < A (ρ, η, s) 2ηs 3ρ(1 s) + W (10) where W = 3ρ(3ρ(1 s) 2 + 4s). ),
48 Stability (cont.) For ρ > 1 3, this generates condition for linear instability for the system: ( 3ρ(1 s) 2s + ) s ( W 3ρ(1 + s) W A < A (ρ, η, s) 2ηs 3ρ(1 s) + W (10) where W = 3ρ(3ρ(1 s) 2 + 4s). When s = 1, the above inequality reduces to Aη + 1 < 3ρ, (11) which agrees with the result from Short et al. ),
49 Stability (cont.) For ρ > 1 3, this generates condition for linear instability for the system: ( 3ρ(1 s) 2s + ) s ( W 3ρ(1 + s) W A < A (ρ, η, s) 2ηs 3ρ(1 s) + W (10) where W = 3ρ(3ρ(1 s) 2 + 4s). When s = 1, the above inequality reduces to Aη + 1 < 3ρ, (11) which agrees with the result from Short et al. Has bifurcations in s, so changing degree of Lévy Flight alters stability. ),
50 Changing Stability with Varying Parameter (a) A stable regime Figure: Different possibilities of the effect of fractional diffusion.
51 Changing Stability with Varying Parameter (a) An unstable regime Figure: Different possibilities of the effect of fractional diffusion.
52 Changing Stability with Varying Parameter (a) Fractional diffusion leads to stability (b) Fractional diffusion leads to stability Figure: Different possibilities of the effect of fractional diffusion.
53 Changing Stability with Varying Parameter (a) Fractional diffusion leads to unstability (b) Fractional diffusion leads to unstability Figure: Different possibilities of the effect of fractional diffusion.
54 Changing Stability with Varying Parameter (a) Fractional diffusion leads (b) Fractional diffusion leads to to stability and then instability stability and then instability Figure: Different possibilities of the effect of fractional diffusion.
55 Changing Stability with Varying Parameter Remark. The fixed points k of σ( k ) with respect s are given by 6(ρ 1 3 k 1 = 1, k 2 = ) + 3A. (12) 2η
56 Changing Stability with Varying Parameter (a) A regime in which k 2 < k 1 = 1 Figure: Different possibilities of the effect of fractional diffusion.
57 Numerical Verfication We need to automate a hot-spot detection.
58 Numerical Verfication We need to automate a hot-spot detection. The variance and its derivative.
59 Numerical Verfication We need to automate a hot-spot detection. The variance and its derivative. The difference of the maximum of the solution in the final frame from A as a percentage of A.
60 Numerical Verfication We need to automate a hot-spot detection. The variance and its derivative. The difference of the maximum of the solution in the final frame from A as a percentage of A. The difference between the maximum and the minimum.
61 Numerical Verfication We need to automate a hot-spot detection. The variance and its derivative. The difference of the maximum of the solution in the final frame from A as a percentage of A. The difference between the maximum and the minimum. An ensemble method.
62 Numerical Verfication (cont.) Figure: A parameter analysis with a bifurcation curve. Fix ρ and η.
63 Numerical Verfication (cont.) Figure: A hot-spot formation when A = 13, s = 0.8.
64 Numerical Verfication (cont.) Figure: A parameter analysis with a bifurcation curve. Fix ρ and η.
65 Numerical Verfication (cont.) Figure: No hot-spot formation when A = 21, s = 0.6
66 Numerical Verfication (cont.)
67 Hot-spot Shape Are the hotspots any different under the two regimes?
68 Hot-spot Shape Are the hotspots any different under the two regimes? Yes and no. Attractiveness hotspots do not change, and first order approximations from Kolokolnikov et al. still work very well.
69 Hot-spot Shape Are the hotspots any different under the two regimes? Yes and no. Attractiveness hotspots do not change, and first order approximations from Kolokolnikov et al. still work very well. But distribution of criminals changes.
70 Hotspot Shape (cont.) Figure: The inner region and the outer region of a hotspot.
71 Hotspot Shape (cont.) Let x = εy, and v = ρ/a 2.
72 Hotspot Shape (cont.) Let x = εy, and v = ρ/a 2. In the inner region ( x < ε), A ε 1 v 1/2 0 w(y), (12) v v 0, (13) where v 0 and v 1 are constants, and w = 2 sech y (same as Kolokolnikov et al.).
73 Hotspot Shape (cont.) Let x = εy, and v = ρ/a 2. In the inner region ( x < ε), A ε 1 v 1/2 0 w(y), (12) v v 0, (13) where v 0 and v 1 are constants, and w = 2 sech y (same as Kolokolnikov et al.). In the outer region (ε x l), we have where A = α + o(1), (14a) v = h 0 (x) + o(1), (14b) s h 0 (x) = ζ = α γ D 0 α 2 < 0, 0 < x l, (h 0) x (±l) = 0, (15)
74 Future Work Analyze the perturbation near hot-spots.
75 Future Work Analyze the perturbation near hot-spots. Study the dynamics of K-hot-spots with Lévy Flights. Brownian based model. Lévy based model.
76 Future Work Analyze the perturbation near hot-spots. Study the dynamics of K-hot-spots with Lévy Flights. Brownian based model. Add the police. Lévy based model.
77 Future Work Analyze the perturbation near hot-spots. Study the dynamics of K-hot-spots with Lévy Flights. Brownian based model. Lévy based model. Add the police. Improve the numerical simulation, especially hotspot detector.
78 Future Work Analyze the perturbation near hot-spots. Study the dynamics of K-hot-spots with Lévy Flights. Brownian based model. Lévy based model. Add the police. Improve the numerical simulation, especially hotspot detector. Analyze weakly-nonlinear stability.
79 Acknowledgments Thanks... Theodore Kolokolnikov Scott McCalla UCLA and REU Program Organizers Harvey Mudd College Mathematics Department for funding Tum Nestor Guillen for creating the Nonlocal Equations Wiki (
80 References M. Short, M. DOrsogna, P. Brantingham, and G. Tita. Measuring and modeling repeat and near-repeat burglary effects. Journal of Quantitative Criminology, 25: , /s M.B. Short, P.J. Brantingham, A.L. Bertozzi, and G.E. Tita. Dissipation and displacement of hotspots in reaction-diffusion models of crime. Proceedings of the National Academy of Sciences of the United States of America, 107: , 2010.
CRIME MODELLING WITH LÉVY FLIGHTS AND THE STABILITY OF HOT-SPOT PATTERNS
CRIME MODELLING WITH LÉVY FLIGHTS AND THE STABILITY OF HOT-SPOT PATTERNS JONAH BRESLAU, TUM CHATURAPRUEK, THEODORE KOLOKOLNIKOV, SCOTT MCCALLA, DANIEL YAZDI ABSTRACT. Working off of the pioneering crime
More informationHot spots in crime model
Hot spots in crime model t=0.0 t=1.8 t=14.5 t=17.1 t=19.9 t=34505.5 0 1 2 3 4 5 Theodore Kolokolnikov Joint works with Jonah Breslau, Tum Chaturapruek, Daniel Yazdi, Scott McCalla, Michael Ward, Juncheng
More informationThe Mathematics of Criminal Behavior: From Agent-Based Models to PDEs. Martin B. Short GA Tech School of Math
The Mathematics of Criminal Behavior: From Agent-Based Models to PDEs Martin B. Short GA Tech School of Math First things first: crime exhibits spatial and/or temporal patterning. April-June 2001 July-September
More informationNONLINEAR PATTERNS IN URBAN CRIME - HOTSPOTS, BIFURCATIONS, AND SUPPRESSION
NONLINEAR PATTERNS IN URBAN CRIME - HOTSPOTS, BIFURCATIONS, AND SUPPRESSION M.B. SHORT, A.L. BERTOZZI, AND P.J. BRANTINGHAM Abstract. We present a weakly nonlinear analysis of our recently developed model
More informationGang rivalry dynamics via coupled point process networks
Gang rivalry dynamics via coupled point process networks George Mohler Department of Mathematics and Computer Science Santa Clara University georgemohler@gmail.com Abstract We introduce a point process
More informationNonlinear Patterns in Urban Crime: Hotspots, Bifurcations, and Suppression
SIAM J. APPLIED DYNAMICAL SYSTEMS Vol. 9, No. 2, pp. 462 483 c 2010 Society for Industrial and Applied Mathematics Nonlinear Patterns in Urban Crime: Hotspots, Bifurcations, and Suppression M. B. Short,
More informationButterflies, Bees & Burglars
Butterflies, Bees & Burglars The foraging behavior of contemporary criminal offenders P. Jeffrey Brantingham 1 Maria-Rita D OrsognaD 2 Jon Azose 3 George Tita 4 Andrea Bertozzi 2 Lincoln Chayes 2 1 UCLA
More informationarxiv: v3 [physics.soc-ph] 8 May 2018
CRIME MODELING WITH TRUNCATED LÉVY FLIGHTS FOR RESIDENTIAL BURGLARY MODELS CHAOHAO PAN arxiv:161.3415v3 physics.soc-ph] 8 May 218 BO LI CHUNTIAN WANG YUQI ZHANG NATHAN GELDNER LI WANG ANDREA L. BERTOZZI
More informationA SOCIAL INTERACTION MODEL FOR CRIME HOT SPOTS
A SOCIAL INTERACTION MODEL FOR CRIME HOT SPOTS Evan C Haskell Division of Math, Science, and Technology Farquhar College of Arts and Sciences Nova Southeastern University Fort Lauderdale, FL 33314 USA
More informationFirst online - August 25, Draft version - November 2, 2016
First online - August 25, 16. Draft version - November 2, 16 Novi Sad J. Math. Vol. XX, No. Y, ZZ,??-?? LÉVY PROCESSES, SUBORDINATORS AND CRIME MODELLING Tijana Levajković 12, Hermann Mena 1 and Martin
More informationCrime modeling with truncated Lévy flights for residential burglary models
Mathematical Models and Methods in Applied Sciences Vol. 28, No. 9 (2018) 1857 1880 c The Author(s) DOI: 10.1142/S0218202518400080 Crime modeling with truncated Lévy flights for residential burglary models
More informationCOPS ON THE DOTS IN A MATHEMATICAL MODEL OF URBAN CRIME AND POLICE RESPONSE. Joseph R. Zipkin. Martin B. Short. Andrea L. Bertozzi
Manuscript submitted to AIMS Journals Volume X, Number 0X, XX 200X doi:10.3934/xx.xx.xx.xx pp. X XX COPS ON THE DOTS IN A MATHEMATICAL MODEL OF URBAN CRIME AND POLICE RESPONSE Joseph R. Zipkin Department
More informationA Convection-Diffusion Model for Gang Territoriality. Abstract
A Convection-Diffusion Model for Gang Territoriality Abdulaziz Alsenafi 1, and Alethea B. T. Barbaro 2, 1 Kuwait University, Faculty of Science, Department of Mathematics, P.O. Box 5969, Safat -13060,
More informationNonlinear convective stability of travelling fronts near Turing and Hopf instabilities
Nonlinear convective stability of travelling fronts near Turing and Hopf instabilities Margaret Beck Joint work with Anna Ghazaryan, University of Kansas and Björn Sandstede, Brown University September
More informationMath 46, Applied Math (Spring 2008): Final
Math 46, Applied Math (Spring 2008): Final 3 hours, 80 points total, 9 questions, roughly in syllabus order (apart from short answers) 1. [16 points. Note part c, worth 7 points, is independent of the
More informationDiffuse interface methods on graphs: Data clustering and Gamma-limits
Diffuse interface methods on graphs: Data clustering and Gamma-limits Yves van Gennip joint work with Andrea Bertozzi, Jeff Brantingham, Blake Hunter Department of Mathematics, UCLA Research made possible
More informationThe Stability of Hot Spot Patterns for Reaction-Diffusion Models of Urban Crime
The Stability of Hot Spot Patterns for Reaction-Diffusion Models of Urban Crime Michael J. Ward (UBC) PIMS Hot Topics: Computational Criminology and UBC (IAM), September 1 Joint With: Theodore Kolokolonikov
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 informationPattern formation in a reaction-diffusion system with space-dependent feed rate
Pattern formation in a reaction-diffusion system with space-dependent feed rate Theodore Kolokolnkov and Juncheng Wei Abstract. We consider a standard reaction-diffusion system (the Schnakenberg model)
More informationHydrodynamics. Stefan Flörchinger (Heidelberg) Heidelberg, 3 May 2010
Hydrodynamics Stefan Flörchinger (Heidelberg) Heidelberg, 3 May 2010 What is Hydrodynamics? Describes the evolution of physical systems (classical or quantum particles, fluids or fields) close to thermal
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 informationIntroduction to multiscale modeling and simulation. Explicit methods for ODEs : forward Euler. y n+1 = y n + tf(y n ) dy dt = f(y), y(0) = y 0
Introduction to multiscale modeling and simulation Lecture 5 Numerical methods for ODEs, SDEs and PDEs The need for multiscale methods Two generic frameworks for multiscale computation Explicit methods
More informationThe Form of Swarms: A Primer of Swarm Equilibria
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-0740484 Andrew J. Bernoff Harvey Mudd
More informationChapter 9b: Numerical Methods for Calculus and Differential Equations. Initial-Value Problems Euler Method Time-Step Independence MATLAB ODE Solvers
Chapter 9b: Numerical Methods for Calculus and Differential Equations Initial-Value Problems Euler Method Time-Step Independence MATLAB ODE Solvers Acceleration Initial-Value Problems Consider a skydiver
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 informationDecember 2, :0 WSPC/INSTRUCTION FILE LocExtDraftFinal LOCAL EXISTENCE AND UNIQUENESS OF SOLUTIONS TO A PDE MODEL FOR CRIMINAL BEHAVIOR
December, 9 8: WSPC/INSTRUCTION FILE LocExtDraftFinal Mathematical Models and Methods in Applied Sciences c World Scientific Publishing Company LOCAL EXISTENCE AND UNIQUENESS OF SOLUTIONS TO A PDE MODEL
More informationLab 5: Nonlinear Systems
Lab 5: Nonlinear Systems Goals In this lab you will use the pplane6 program to study two nonlinear systems by direct numerical simulation. The first model, from population biology, displays interesting
More informationPattern Formation in a Reaction-Diffusion System with Space-Dependent Feed Rate
SIAM REVIEW Vol. 60, No. 3, pp. 000 000 c 08 Society for Industrial and Applied Mathematics Pattern Formation in a Reaction-Diffusion System with Space-Dependent Feed Rate Theodore Kolokolnikov Juncheng
More informationSynchronization Transitions in Complex Networks
Synchronization Transitions in Complex Networks Y. Moreno 1,2,3 1 Institute for Biocomputation and Physics of Complex Systems (BIFI) University of Zaragoza, Zaragoza 50018, Spain 2 Department of Theoretical
More informationPhase transition in a model for territorial development
Phase transition in a model for territorial development Alethea Barbaro Joint with Abdulaziz Alsenafi Supported by the NSF through Grant No. DMS-39462 Case Western Reserve University March 23, 27 A. Barbaro
More informationThe Stability of Steady-State Hot-Spot Patterns for a Reaction-Diffusion Model of Urban Crime
Preprint The Stability of Steady-State Hot-Spot Patterns for a Reaction-Diffusion Model of Urban Crime T. KOLOKOLNIKOV, M. J. WARD, and J. WEI Theodore Kolokolnikov; Department of Mathematics, Dalhousie
More informationObstacle problems for nonlocal operators
Obstacle problems for nonlocal operators Camelia Pop School of Mathematics, University of Minnesota Fractional PDEs: Theory, Algorithms and Applications ICERM June 19, 2018 Outline Motivation Optimal regularity
More informationUncommon (in)stability properties in the numerical method of characteristic applied to hyperbolic equations
University of Vermont ScholarWorks @ UVM Graduate College Dissertations and Theses Dissertations and Theses 2016 Uncommon (in)stability properties in the numerical method of characteristic applied to hyperbolic
More informationCOMPLETENESS THEOREM FOR THE DISSIPATIVE STURM-LIOUVILLE OPERATOR ON BOUNDED TIME SCALES. Hüseyin Tuna
Indian J. Pure Appl. Math., 47(3): 535-544, September 2016 c Indian National Science Academy DOI: 10.1007/s13226-016-0196-1 COMPLETENESS THEOREM FOR THE DISSIPATIVE STURM-LIOUVILLE OPERATOR ON BOUNDED
More informationNumerical Methods for Differential Equations Mathematical and Computational Tools
Numerical Methods for Differential Equations Mathematical and Computational Tools Gustaf Söderlind Numerical Analysis, Lund University Contents V4.16 Part 1. Vector norms, matrix norms and logarithmic
More informationPDEs, part 3: Hyperbolic PDEs
PDEs, part 3: Hyperbolic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Hyperbolic equations (Sections 6.4 and 6.5 of Strang). Consider the model problem (the
More information3 Mathematical modeling of the torsional dynamics of a drill string
3 Mathematical modeling of the torsional dynamics of a drill string 3.1 Introduction Many works about torsional vibrations on drilling systems [1, 12, 18, 24, 41] have been published using different numerical
More informationMonsuru Adepeju 1 and Andy Evans 2. School of Geography, University of Leeds, LS21 1HB 1
Investigating the impacts of training data set length (T) and the aggregation unit size (M) on the accuracy of the self-exciting point process (SEPP) hotspot method Monsuru Adepeju 1 and Andy Evans 2 1,
More informationAnalysis Preliminary Exam Workshop: Hilbert Spaces
Analysis Preliminary Exam Workshop: Hilbert Spaces 1. Hilbert spaces A Hilbert space H is a complete real or complex inner product space. Consider complex Hilbert spaces for definiteness. If (, ) : H H
More informationLecture Notes for Math 251: ODE and PDE. Lecture 30: 10.1 Two-Point Boundary Value Problems
Lecture Notes for Math 251: ODE and PDE. Lecture 30: 10.1 Two-Point Boundary Value Problems Shawn D. Ryan Spring 2012 Last Time: We finished Chapter 9: Nonlinear Differential Equations and Stability. Now
More informationExam in TMA4195 Mathematical Modeling Solutions
Norwegian University of Science and Technology Department of Mathematical Sciences Page of 9 Exam in TMA495 Mathematical Modeling 6..07 Solutions Problem a Here x, y are two populations varying with time
More informationA review of stability and dynamical behaviors of differential equations:
A review of stability and dynamical behaviors of differential equations: scalar ODE: u t = f(u), system of ODEs: u t = f(u, v), v t = g(u, v), reaction-diffusion equation: u t = D u + f(u), x Ω, with boundary
More informationFree Entropy for Free Gibbs Laws Given by Convex Potentials
Free Entropy for Free Gibbs Laws Given by Convex Potentials David A. Jekel University of California, Los Angeles Young Mathematicians in C -algebras, August 2018 David A. Jekel (UCLA) Free Entropy YMC
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 informationStability of Krylov Subspace Spectral Methods
Stability of Krylov Subspace Spectral Methods James V. Lambers Department of Energy Resources Engineering Stanford University includes joint work with Patrick Guidotti and Knut Sølna, UC Irvine Margot
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 informationMS: Nonlinear Wave Propagation in Singular Perturbed Systems
MS: Nonlinear Wave Propagation in Singular Perturbed Systems P. van Heijster: Existence & stability of 2D localized structures in a 3-component model. Y. Nishiura: Rotational motion of traveling spots
More informationSolitons in a macroscopic traffic model
Solitons in a macroscopic traffic model P. Saavedra R. M. Velasco Department of Mathematics, Universidad Autónoma Metropolitana, Iztapalapa, 093 México, (e-mail: psb@xanum.uam.mx). Department of Physics,
More informationM.ARCIERO, G.LADAS AND S.W.SCHULTZ
Georgian Mathematical Journal 1(1994), No. 3, 229-233 SOME OPEN PROBLEMS ABOUT THE SOLUTIONS OF THE DELAY DIFFERENCE EQUATION x n+1 = A/x 2 n + 1/x p n k M.ARCIERO, G.LADAS AND S.W.SCHULTZ Abstract. We
More informationThe Motion of a Hanging Chain
The free end chain with little kick and free falling Altrichter Stacey Qiang Chen Department of Mathematical Sciences University of Delaware December 3, 2007 Outline 1 Why this problem? The Mathematical
More informationFokker-Planck Equation on Graph with Finite Vertices
Fokker-Planck Equation on Graph with Finite Vertices January 13, 2011 Jointly with S-N Chow (Georgia Tech) Wen Huang (USTC) Hao-min Zhou(Georgia Tech) Functional Inequalities and Discrete Spaces Outline
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 information35. RESISTIVE INSTABILITIES: CLOSING REMARKS
35. RESISTIVE INSTABILITIES: CLOSING REMARKS In Section 34 we presented a heuristic discussion of the tearing mode. The tearing mode is centered about x = 0, where F x that B 0 ( ) = k! B = 0. Note that
More informationAdaptive evolution : a population approach Benoît Perthame
Adaptive evolution : a population approach Benoît Perthame 30 20 t 50 t 10 x 0 1 0.5 0 0.5 1 x 0 1 0.5 0 0.5 1 OUTLINE OF THE LECTURE DIRECT COMPETITION AND POLYMORPHIC CONCENTRATIONS I. Direct competition
More informationA Model of Evolutionary Dynamics with Quasiperiodic Forcing
paper presented at Society for Experimental Mechanics (SEM) IMAC XXXIII Conference on Structural Dynamics February 2-5, 205, Orlando FL A Model of Evolutionary Dynamics with Quasiperiodic Forcing Elizabeth
More informationMatrix stabilization using differential equations.
Matrix stabilization using differential equations. Nicola Guglielmi Universitá dell Aquila and Gran Sasso Science Institute, Italia NUMOC-2017 Roma, 19 23 June, 2017 Inspired by a joint work with Christian
More informationThe circular law. Lewis Memorial Lecture / DIMACS minicourse March 19, Terence Tao (UCLA)
The circular law Lewis Memorial Lecture / DIMACS minicourse March 19, 2008 Terence Tao (UCLA) 1 Eigenvalue distributions Let M = (a ij ) 1 i n;1 j n be a square matrix. Then one has n (generalised) eigenvalues
More informationFinite difference method for solving Advection-Diffusion Problem in 1D
Finite difference method for solving Advection-Diffusion Problem in 1D Author : Osei K. Tweneboah MATH 5370: Final Project Outline 1 Advection-Diffusion Problem Stationary Advection-Diffusion Problem in
More informationFractional Laplacian
Fractional Laplacian Grzegorz Karch 5ème Ecole de printemps EDP Non-linéaire Mathématiques et Interactions: modèles non locaux et applications Ecole Supérieure de Technologie d Essaouira, du 27 au 30 Avril
More informationUsing web-based Java pplane applet to graph solutions of systems of differential equations
Using web-based Java pplane applet to graph solutions of systems of differential equations Our class project for MA 341 involves using computer tools to analyse solutions of differential equations. This
More informationStability of the Parareal Algorithm
Stability of the Parareal Algorithm Gunnar Andreas Staff and Einar M. Rønquist Norwegian University of Science and Technology Department of Mathematical Sciences Summary. We discuss the stability of the
More informationSelf-exciting point process modeling of crime
Self-exciting point process modeling of crime G. O. Mohler M. B. Short P. J. Brantingham F. P. Schoenberg G. E. Tita Abstract Highly clustered event sequences are observed in certain types of crime data,
More informationFall Exam II. Wed. Nov. 9, 2005
Fall 005 10.34 Eam II. Wed. Nov. 9 005 (SOLUTION) Read through the entire eam before beginning work and budget your time. You are researching drug eluting stents and wish to understand better how the drug
More informationSelective Manipulation of Molecules by Electrostatic Force and Detection of Single Molecules in Aqueous Solution
Supporting Information Selective Manipulation of Molecules by Electrostatic Force and Detection of Single Molecules in Aqueous Solution Zhongbo Yan, Ming Xia, Pei Zhang, and Ya-Hong Xie* Department of
More informationMath 46, Applied Math (Spring 2009): Final
Math 46, Applied Math (Spring 2009): Final 3 hours, 80 points total, 9 questions worth varying numbers of points 1. [8 points] Find an approximate solution to the following initial-value problem which
More informationGeneral introduction to Hydrodynamic Instabilities
KTH ROYAL INSTITUTE OF TECHNOLOGY General introduction to Hydrodynamic Instabilities L. Brandt & J.-Ch. Loiseau KTH Mechanics, November 2015 Luca Brandt Professor at KTH Mechanics Email: luca@mech.kth.se
More informationLong-wave Instability in Anisotropic Double-Diffusion
Long-wave Instability in Anisotropic Double-Diffusion Jean-Luc Thiffeault Institute for Fusion Studies and Department of Physics University of Texas at Austin and Neil J. Balmforth Department of Theoretical
More informationMCE/EEC 647/747: Robot Dynamics and Control. Lecture 8: Basic Lyapunov Stability Theory
MCE/EEC 647/747: Robot Dynamics and Control Lecture 8: Basic Lyapunov Stability Theory Reading: SHV Appendix Mechanical Engineering Hanz Richter, PhD MCE503 p.1/17 Stability in the sense of Lyapunov A
More information5 Applying the Fokker-Planck equation
5 Applying the Fokker-Planck equation We begin with one-dimensional examples, keeping g = constant. Recall: the FPE for the Langevin equation with η(t 1 )η(t ) = κδ(t 1 t ) is = f(x) + g(x)η(t) t = x [f(x)p
More informationPhysically Based Simulations (on the GPU)
Physically Based Simulations (on the GPU) (some material from slides of Mark Harris) CS535 Fall 2014 Daniel G. Aliaga Department of Computer Science Purdue University Simulating the world Floating point
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 informationPartial differential equations
Partial differential equations Many problems in science involve the evolution of quantities not only in time but also in space (this is the most common situation)! We will call partial differential equation
More informationDimensional Analysis - Concepts
Dimensional Analysis - Concepts Physical quantities: R j = v(r j )[R j ] = value unit, j = 1,..., m. Units: Dimension matrix of R 1,, R m : A = Change of units change of values: [R j ] = F a 1j 1 F a nj
More informationIntroduction to Mathematical Physics
Introduction to Mathematical Physics Methods and Concepts Second Edition Chun Wa Wong Department of Physics and Astronomy University of California Los Angeles OXFORD UNIVERSITY PRESS Contents 1 Vectors
More informationFINAL EXAM MATH303 Theory of Ordinary Differential Equations. Spring dx dt = x + 3y dy dt = x y.
FINAL EXAM MATH0 Theory of Ordinary Differential Equations There are 5 problems on 2 pages. Spring 2009. 25 points Consider the linear plane autonomous system x + y x y. Find a fundamental matrix of the
More informationSaturation of Information Exchange in Locally Connected Pulse-Coupled Oscillators
Saturation of Information Exchange in Locally Connected Pulse-Coupled Oscillators Will Wagstaff School of Computer Science, Georgia Institute of Technology, Atlanta, Georgia 30332, USA (Dated: 13 December
More informationModal Decomposition Methods on Aerodynamic Flows
498 Modal Decomposition Methods on Aerodynamic Flows Mohammad Niaz Murshed University of Michigan Department of Aerospace Engineering February 5, 2018 Abstract- This paper discusses the significance of
More informationMath 322. Spring 2015 Review Problems for Midterm 2
Linear Algebra: Topic: Linear Independence of vectors. Question. Math 3. Spring Review Problems for Midterm Explain why if A is not square, then either the row vectors or the column vectors of A are linearly
More informationBasic Concepts and the Discovery of Solitons p. 1 A look at linear and nonlinear signatures p. 1 Discovery of the solitary wave p. 3 Discovery of the
Basic Concepts and the Discovery of Solitons p. 1 A look at linear and nonlinear signatures p. 1 Discovery of the solitary wave p. 3 Discovery of the soliton p. 7 The soliton concept in physics p. 11 Linear
More informationHardening nonlinearity effects on forced vibration of viscoelastic dynamical systems to external step perturbation field
Int. J. Adv. Appl. Math. and Mech. 3(2) (215) 16 32 (ISSN: 2347-2529) IJAAMM Journal homepage: www.ijaamm.com International Journal of Advances in Applied Mathematics and Mechanics Hardening nonlinearity
More informationMATH 100 Introduction to the Profession
MATH 100 Introduction to the Profession Differential Equations in MATLAB Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2012 fasshauer@iit.edu MATH 100 ITP 1 What
More informationDecentralized Stabilization of Heterogeneous Linear Multi-Agent Systems
1 Decentralized Stabilization of Heterogeneous Linear Multi-Agent Systems Mauro Franceschelli, Andrea Gasparri, Alessandro Giua, and Giovanni Ulivi Abstract In this paper the formation stabilization problem
More informationUniversity of Bristol - Explore Bristol Research. Link to publication record in Explore Bristol Research PDF-document.
Dobra, T., Lawrie, A., & Dalziel, S. B. (2016). Nonlinear Interactions of Two Incident Internal Waves. 1-8. Paper presented at VIIIth International Symposium on Stratified Flows, San Diego, United States.
More informationof the Schnakenberg model
Pulse motion in the semi-strong limit of the Schnakenberg model Newton Institute 2005 Jens Rademacher, Weierstraß Institut Berlin joint work with Michael Ward (UBC) Angelfish 2, 6, 12 months old [Kondo,
More informationNonlinear Dynamic Systems Homework 1
Nonlinear Dynamic Systems Homework 1 1. A particle of mass m is constrained to travel along the path shown in Figure 1, which is described by the following function yx = 5x + 1x 4, 1 where x is defined
More informationAPPLICATIONS OF FD APPROXIMATIONS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS
LECTURE 10 APPLICATIONS OF FD APPROXIMATIONS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS Ordinary Differential Equations Initial Value Problems For Initial Value problems (IVP s), conditions are specified
More informationOrdinary Differential Equations II
Ordinary Differential Equations II February 9 217 Linearization of an autonomous system We consider the system (1) x = f(x) near a fixed point x. As usual f C 1. Without loss of generality we assume x
More informationComputing Spectra of Linear Operators Using Finite Differences
Computing Spectra of Linear Operators Using Finite Differences J. Nathan Kutz Department of Applied Mathematics University of Washington Seattle, WA 98195-2420 Email: (kutz@amath.washington.edu) Stability
More information154 Chapter 9 Hints, Answers, and Solutions The particular trajectories are highlighted in the phase portraits below.
54 Chapter 9 Hints, Answers, and Solutions 9. The Phase Plane 9.. 4. The particular trajectories are highlighted in the phase portraits below... 3. 4. 9..5. Shown below is one possibility with x(t) and
More informationα 2 )(k x ũ(t, η) + k z W (t, η)
1 3D Poiseuille Flow Over the next two lectures we will be going over the stabilization of the 3-D Poiseuille flow. For those of you who havent had fluids, that means we have a channel of fluid that is
More informationNumerical Methods for Partial Differential Equations. CAAM 452 Spring 2005 Lecture 2 Instructor: Tim Warburton
Numerical Methods for Partial Differential Equations CAAM 452 Spring 2005 Lecture 2 Instructor: Tim Warburton Note on textbook for finite difference methods Due to the difficulty some students have experienced
More informationModeling Spatial Relationships Using Regression Analysis
Esri International User Conference San Diego, California Technical Workshops July 24, 2012 Modeling Spatial Relationships Using Regression Analysis Lauren M. Scott, PhD Lauren Rosenshein Bennett, MS Answering
More informationSection 3.9. Matrix Norm
3.9. Matrix Norm 1 Section 3.9. Matrix Norm Note. We define several matrix norms, some similar to vector norms and some reflecting how multiplication by a matrix affects the norm of a vector. We use matrix
More informationSelf-exciting point process modeling of crime
Self-exciting point process modeling of crime G. O. Mohler M. B. Short P. J. Brantingham F. P. Schoenberg G. E. Tita Abstract Highly clustered event sequences are observed in certain types of crime data,
More informationNonlinear stability of time-periodic viscous shocks. Margaret Beck Brown University
Nonlinear stability of time-periodic viscous shocks Margaret Beck Brown University Motivation Time-periodic patterns in reaction-diffusion systems: t x Experiment: chemical reaction chlorite-iodite-malonic-acid
More information1. (i) Determine how many periodic orbits and equilibria can be born at the bifurcations of the zero equilibrium of the following system:
1. (i) Determine how many periodic orbits and equilibria can be born at the bifurcations of the zero equilibrium of the following system: ẋ = y x 2, ẏ = z + xy, ż = y z + x 2 xy + y 2 + z 2 x 4. (ii) Determine
More informationCOMPARISON OF TWO METHODS TO SOLVE PRESSURES IN SMALL VOLUMES IN REAL-TIME SIMULATION OF A MOBILE DIRECTIONAL CONTROL VALVE
COMPARISON OF TWO METHODS TO SOLVE PRESSURES IN SMALL VOLUMES IN REAL-TIME SIMULATION OF A MOBILE DIRECTIONAL CONTROL VALVE Rafael ÅMAN*, Heikki HANDROOS*, Pasi KORKEALAAKSO** and Asko ROUVINEN** * Laboratory
More informationME8230 Nonlinear Dynamics
ME8230 Nonlinear Dynamics Lecture 1, part 1 Introduction, some basic math background, and some random examples Prof. Manoj Srinivasan Mechanical and Aerospace Engineering srinivasan.88@osu.edu Spring mass
More informationAbstracts. Furstenberg The Dynamics of Some Arithmetically Generated Sequences
CHAOS AND DISORDER IN MATHEMATICS AND PHYSICS Monday 10:00-11:00 Okounkov Algebraic geometry of random surfaces 11:30-12:30 Furstenberg Dynamics of Arithmetically Generated Sequences 12:30-14:30 lunch
More informationExponential integrators and functions of the matrix exponential
Exponential integrators and functions of the matrix exponential Paul Matthews, Stephen Cox, Hala Ashi and Linda Cummings School of Mathematical Sciences, University of Nottingham, UK Introduction to exponential
More information