SciCADE 2013
International Conference on Scientific Computation and Differential Equations
September 16-20, 2013, Valladolid (Spain)

Contributed Talk

Show full list of Contributed Talks Show talk context (CS05)

Stochastic correction of kinetic energy spectra in fluids

J. Frank, J. Bajars, K. Myerscough and B. Leimkuhler

Both conservative and diffusive truncations of inviscid fluid models fail to capture the correct power law scaling of the kinetic energy spectrum at small scales. This scaling is important both for downscale transport of vorticity and energy, as for small upscale `backscatter' that influences variance. In this talk we discuss a simple correction using a stochastic thermostat approach from molecular dynamics. Thermostats are used in MD to perturb dynamics such that trajectories are ergodic in the canonical Gibbs measure (constant temperature). To apply these methods to discretized fluids, several challenges must be met: (i) we perturb only the smallest scales, hence controllability must be established; (ii) we are given expectations (kinetic energy spectrum) instead of invariant measure; (iii) we have to deal with forcing at low wave numbers; (iv) experience from heat conduction problems suggests that artifacts may occur. We report recent progress on these fronts [1, 2].

[1] J. Bajars, J. Frank and B. Leimkuhler, Weakly coupled heat bath models for Gibbs-like invariant states in nonlinear wave equations, Nonlinearity (2013), in press.
[2] S. Dubinkina, J. Frank and B. Leimkuhler, Simplified Modelling of a Thermal Bath, with Application to a Fluid Vortex System, SIAM Multiscale Model. Simul. 8 (2010) 1882-1901.

Organized by         Universidad de Valladolid     IMUVA