Blending Modified Gaussian Closure and Non-Gaussian Reduced Subspace Methods for Turbulent Dynamical Systems

Turbulent dynamical systems are characterized by persistent instabilities which are balanced by nonlinear dynamics that continuously transfer energy to the stable modes. To model this complex statistical equilibrium in the context of uncertainty quantification all dynamical components (unstable mode...

Full description

Bibliographic Details
Main Authors: Sapsis, Themistoklis Panagiotis, Majda, Andrew J.
Other Authors: Massachusetts Institute of Technology. Department of Mechanical Engineering
Format: Article
Language:English
Published: Springer US 2016
Online Access:http://hdl.handle.net/1721.1/103349
https://orcid.org/0000-0003-0302-0691
Description
Summary:Turbulent dynamical systems are characterized by persistent instabilities which are balanced by nonlinear dynamics that continuously transfer energy to the stable modes. To model this complex statistical equilibrium in the context of uncertainty quantification all dynamical components (unstable modes, nonlinear energy transfers, and stable modes) are equally crucial. Thus, order-reduction methods present important limitations. On the other hand uncertainty quantification methods based on the tuning of the non-linear energy fluxes using steady-state information (such as the modified quasilinear Gaussian (MQG) closure) may present discrepancies in extreme excitation scenarios. In this paper we derive a blended framework that links inexpensive second-order uncertainty quantification schemes that model the full space (such as MQG) with high order statistical models in specific reduced-order subspaces. The coupling occurs in the energy transfer level by (i) correcting the nonlinear energy fluxes in the full space using reduced subspace statistics, and (ii) by modifying the reduced-order equations in the subspace using information from the full space model. The results are illustrated in two strongly unstable systems under extreme excitations. The blended method allows for the correct prediction of the second-order statistics in the full space and also the correct modeling of the higher-order statistics in reduced-order subspaces.