A finite element framework for simulating residential burglary in realistic urban geometries

We consider a partial differential equation (PDE) model to predict residential burglary derived from a probabilistic agent-based model through a mean-field limit operation. The PDE model is a nonlinear, coupled system of two equations in two variables (attractiveness of residential sites and density of criminals), similar to the Keller–Segel model for aggregation based on chemotaxis. Unlike previous works, which applied periodic boundary conditions, we enforce boundary conditions that arise natu