Non-Markovian Closure Models for Large Eddy Simulations using the Mori-Zwanzig Formalism

Abstract

This work uses the Mori-Zwanzig (M-Z) formalism, a concept originating from non-equilibrium statistical mechanics, as a basis for the development of coarse-grained models of turbulence. The mechanics of the generalized Langevin equation (GLE) are considered and insight gained from the orthogonal dynamics equation is used as a starting point for model development. A class of sub-grid models is considered which represent non-local behavior via a finite memory approximation (Stinis, P., "Mori-Zwanzig reduced models for uncertainty quantification I: Parametric uncertainty," \textit{arXiv:1211.4285}, 2012.), the length of which is determined using a heuristic that is related to the spectral radius of the Jacobian of the resolved variables. The resulting models are intimately tied to the underlying numerical resolution and are capable of approximating non-Markovian effects. Numerical experiments on the Burgers equation demonstrate that the M-Z-based models can accurately predict the temporal evolution of the total kinetic energy and the total dissipation rate at varying mesh resolutions. The trajectory of each resolved mode in phase-space is accurately predicted for cases where the coarse-graining is moderate. LES of homogeneous isotropic turbulence and the Taylor Green Vortex show that the M-Z-based models are able to provide excellent predictions, accurately capturing the sub-grid contribution to energy transfer. Lastly, LES of fully developed channel flow demonstrate the applicability of M-Z-based models to non-decaying problems. It is notable that the form of the closure is not imposed by the modeler, but is rather derived from the mathematics of the coarse-graining, highlighting the potential of M-Z-based techniques to define LES closures.

arXiv comment32 pages, this is a pre-print of an article submitted to physical review fluids

References in corpus (2)

• Application of the t-model of optimal prediction to the estimation of the rate of decay of solutions of the Euler equations in two and three dimensions
• Renormalized reduced models for singular PDEs

Cited by in corpus (19)

• Perspectives on Machine Learning-augmented Reynolds-averaged and Large Eddy Simulation Models of Turbulence
• Physics-Informed Probabilistic Learning of Linear Embeddings of Non-linear Dynamics With Guaranteed Stability
• Data-driven model reduction, Wiener projections, and the Koopman-Mori-Zwanzig formalism
• A Priori Estimation Of Memory Effects In Coarse-Grained Nonlinear Systems Using The Mori-Zwanzig Formalism
• A Dynamic Subgrid Scale Model for Large Eddy Simulations Based on the Mori-Zwanzig Formalism
• Reduced order models for Lagrangian hydrodynamics
• Non-modal analysis of spectral element methods: Towards accurate and robust large-eddy simulations
• Theoretical tools for understanding the climate crisis from Hasselmann's program and beyond
• Predictive Reduced Order Modeling of Chaotic Multi-scale Problems Using Adaptively Sampled Projections
• Learning Closed-form Equations for Subgrid-scale Closures from High-fidelity Data: Promises and Challenges
• Data Driven Learning of Mori-Zwanzig Operators for Isotropic Turbulence
• Generalized Langevin equations for systems with local interactions
• Understanding probability and irreversibility in the Mori-Zwanzig projection operator formalism
• Nonlocality of Mean Scalar Transport in Two-Dimensional Rayleigh-Taylor Instability Using the Macroscopic Forcing Method
• Variational Multiscale Closures for Finite Element Discretizations Using the Mori-Zwanzig Approach
• Optimal renormalization of multi-scale systems
• Hypoellipticity and the Mori-Zwanzig formulation of stochastic differential equations
• Efficient Adjoint Petrov-Galerkin Reduced Order Models for fluid flows governed by the incompressible Navier-Stokes equations
• Data-driven Mori-Zwanzig modeling of Lagrangian particle dynamics in turbulent flows