Search

A resolution-independent data-driven stochastic parametrization method for subgrid-scale processes in coarsened fluid descriptions is proposed. The method enables the inclusion of high-fidelity data into the coarsened flow model, thereby enabling accurate simulations also with the coarser representation. The small-scale parametrization is introduced at the level of the Fourier coefficients of the coarsened numerical solution. It is designed to reproduce the kinetic energy spectra observed in high-fidelity data of the same system. The approach is based on a control feedback term reminiscent of continuous data assimilation. The method relies solely on the availability of high-fidelity data from a statistically steady state. No assumptions are made regarding the adopted discretization method or the selected coarser resolution. The performance of the method is assessed for the two-dimensional Euler equations on the sphere. Applying the method at two significantly coarser resolutions yields good results for the mean and variance of the Fourier coefficients. Stable and accurate large-scale dynamics can be simulated over long integration times.

A data-driven turbulence model for coarse-grained numerical simulations of two-dimensional Rayleigh-Bénard convection is proposed. The model starts from high-fidelity data and is based on adjusting the Fourier coefficients of the numerical solution, with the aim of accurately reproducing the kinetic energy spectra as seen in the high-fidelity reference findings. No assumptions about the underlying PDE or numerical discretization are used in the formulation of the model. We also develop a constraint on the heat flux to guarantee accurate Nusselt number estimates on coarse computational grids and high Rayleigh numbers. Model performance is assessed in coarse numerical simulations at $Ra=10^{10}$. We focus on key features including kinetic energy spectra, wall-normal flow statistics, and global flow statistics. The method of data-driven modeling of small scale flow dynamics is found to reproduce the reference kinetic energy spectra well across all scales and yields good results for flow statistics and average heat transfer., 23 pages, 10 figures

Casimir preserving integrators for stochastic Lie-Poisson equations with Stratonovich noise are developed extending Runge-Kutta Munthe-Kaas methods. The underlying Lie-Poisson structure is preserved along stochastic trajectories. A related stochastic differential equation on the Lie algebra is derived. The solution of this differential equation updates the evolution of the Lie-Poisson dynamics by means of the exponential map. The constructed numerical method conserves Casimir-invariants exactly, which is important for long time integration. This is illustrated numerically for the case of the stochastic heavy top and the stochastic sine-Euler equations., Comment: 27 pages, 9 figures, fifth version, all comments are welcome!