A dynamical systems approach to simulating macroscale spatial dynamics in multiple dimensions

Abstract

Developments in dynamical systems theory provide new support for the macroscale modelling of pdes and other microscale systems such as lattice Boltzmann, Monte Carlo or molecular dynamics simulators. By systematically resolving subgrid microscale dynamics the dynamical systems approach constructs accurate closures of macroscale discretisations of the microscale system. Here we specifically explore reaction–diffusion problems in two spatial dimensions as a prototype of generic systems in multiple dimensions. Our approach unifies into one the discrete modelling of systems governed by known pdes and the ‘equation-free’ macroscale modelling of microscale simulators efficiently executing only on small patches of the spatial domain. Centre manifold theory ensures that a closed model exists on the macroscale grid, is emergent, and is systematically approximated. Dividing space into either overlapping finite elements or spatially separated small patches, the specially crafted inter-element/patch coupling also ensures that the constructed discretisations are consistent with the microscale system/pde to as high an order as desired. Computer algebra handles the considerable algebraic details, as seen in the specific application to the Ginzburg–Landau pde. However, higher-order models in multiple dimensions require a mixed numerical and algebraic approach that is also developed. The modelling here may be straightforwardly adapted to a wide class of reaction–diffusion pdes and lattice equations in multiple space dimensions. When applied to patches of microscopic simulations our coupling conditions promise efficient macroscale simulation.