Based on the thermodynamic concept of a reservoir, we investigate a computational model for interaction with unresolved degrees of freedom (a thermal bath). We assume that a finite restricted system can be modelled by a generalized canonical ensemble, described by a density which is a smooth function of the energy of the restricted system. A thermostat is constructed to continuously perturb the resolved dynamics, while leaving the desired equilibrium distribution invariant. We build on a thermostatting framework developed and tested in the setting of molecular dynamics, using stochastic perturbations to control (and stabilize) the invariant measure. We also apply these techniques in the setting of a simplified point vortex flow on a disc, in which a modified Gibbs distribution (modelling a finite, rather than infinite, bath of weak vortices) provides a regularizing formulation for restricted system dynamics. Numerical experiments, effectively replacing many vortices by a few artificial degrees of freedom, are in excellent agreement with the two-scale simulations of Bühler [Phys. Fluids, 14 (2002), pp. 2139–2149].
, , , , ,
S.I.A.M.
Multiscale Modeling and Simulation
Symplectic Integration of Atmospheric Dynamics: Long-term Statistical Accuracy for Ensemble Climate Simulations
Computational Dynamics

Dubinkina, S., Frank, J., & Leimkuhler, B. (2010). Simplified modelling of a thermal bath, with application to a fluid vortex system. Multiscale Modeling and Simulation, 8(2010), 1882–1901.