Equation-free modeling of evolving diseases: Coarse-grained computations with individual-based models

We demonstrate how direct simulation of stochastic, individual-based models can be combined with continuum numerical analysis techniques to study the dynamics of evolving diseases. % Sidestepping the necessity of obtaining explicit population-level models, the approach analyzes the (unavailable in closed form) `coarse' macroscopic equations, estimating the necessary quantities through appropriately initialized, short `bursts' of individual-based dynamic simulation. % We illustrate this approach by analyzing a stochastic and discrete model for the evolution of disease agents caused by point mutations within individual hosts. % Building up from classical SIR and SIRS models, our example uses a one-dimensional lattice for variant space, and assumes a finite number of individuals. % Macroscopic computational tasks enabled through this approach include stationary state computation, coarse projective integration, parametric continuation and stability analysis.Comment: 16 pages, 8 figure

iMapD: intrinsic Map Dynamics exploration for uncharted effective free energy landscapes

We describe and implement iMapD, a computer-assisted approach for accelerating the exploration of uncharted effective Free Energy Surfaces (FES), and more generally for the extraction of coarse-grained, macroscopic information from atomistic or stochastic (here Molecular Dynamics, MD) simulations. The approach functionally links the MD simulator with nonlinear manifold learning techniques. The added value comes from biasing the simulator towards new, unexplored phase space regions by exploiting the smoothness of the (gradually, as the exploration progresses) revealed intrinsic low-dimensional geometry of the FES

A manifold learning approach to model reduction in combustion

We use a relatively recent nonlinear manifold learning technique (diffusion maps) to parameterize low dimensional attracting manifolds arising in the description of detailed chemical kinetics mechanisms. With no a priori knowledge about the shape and dimension of the manifold, such an approach provides a way of solving a reduced (and less stiff) set of equations in terms of automatically detected slow variables. Advantages as well as disadvantages of the approach are discussed