Spatial localization for nonlinear dynamical stochastic models for excitable media
arXiv:1901.07318
Abstract
Nonlinear dynamical stochastic models are ubiquitous in different areas. Excitable media models are typical examples with large state dimensions. Their statistical properties are often of great interest but are also very challenging to compute. In this article, a theoretical framework to understand the spatial localization for a large class of stochastically coupled nonlinear systems in high dimensions is developed. Rigorous mathematical theories show the covariance decay behavior due to both local and nonlocal effects, which result from the diffusion and the mean field interaction, respectively. The analysis is based on a comparison with an appropriate linear surrogate model, of which the covariance propagation can be computed explicitly. Two important applications of these theoretical results are discussed. They are the spatial averaging strategy for efficiently sampling the covariance matrix and the localization technique in data assimilation. Test examples of a surrogate linear model and a stochastically coupled FitzHugh-Nagumo model for excitable media are adopted to validate the theoretical results. The latter is also used for a systematical study of the spatial averaging strategy in efficiently sampling the covariance matrix in different dynamical regimes.
arXiv comment34 pages, 9 figures
References in corpus (9)
- Data Assimilation: A Mathematical Introduction
- Dimension-independent likelihood-informed MCMC
- Likelihood-informed dimension reduction for nonlinear inverse problems
- Nonlinear stability and ergodicity of ensemble based Kalman filters
- Efficient Statistically Accurate Algorithms for the Fokker-Planck Equation in Large Dimensions
- Localization of response functions of spiral waves in the FitzHugh-Nagumo system
- A computational framework for infinite-dimensional Bayesian inverse problems. Part I: The linearized case, with application to global seismic inversion
- MCMC Methods for Functions: Modifying Old Algorithms to Make Them Faster
- A multilevel Monte Carlo method for a class of McKean-Vlasov processes