A statistically accurate modified quasilinear Gaussian closure for uncertainty quantification in turbulent dynamical systems

Themistoklis P. Sapsis, Andrew J. Majda

Research output: Contribution to journalArticle

Abstract

We develop a novel second-order closure methodology for uncertainty quantification in damped forced nonlinear systems with high dimensional phase-space that possess a high-dimensional chaotic attractor. We focus on turbulent systems with quadratic nonlinearities where the finite size of the attractor is caused exclusively by the synergistic activity of persistent, linearly unstable directions and a nonlinear energy transfer mechanism. We first illustrate how existing UQ schemes that rely on the Gaussian assumption will fail to perform reliable UQ in the presence of unstable dynamics. To overcome these difficulties, a modified quasilinear Gaussian (MQG) closure is developed in two stages. First we exploit exact statistical relations between second order correlations and third order moments in statistical equilibrium in order to decompose the energy flux at equilibrium into precise additional damping and enhanced noise on suitable modes, while preserving statistical symmetries; in the second stage, we develop a nonlinear MQG dynamical closure which has this statistical equilibrium behavior as a stable fixed point of the dynamics. Our analysis, UQ schemes, and conclusions are illustrated through a specific toy-model, the forty-modes Lorenz 96 system, which despite its simple formulation, presents strongly turbulent behavior with a large number of unstable dynamical components in a variety of chaotic regimes. A suitable version of MQG successfully captures the mean and variance, in transient dynamics with initial data far from equilibrium and with large random fluctuations in forcing, very cheaply at the cost of roughly two ensemble members in a Monte-Carlo simulation.

Original languageEnglish (US)
Pages (from-to)34-45
Number of pages12
JournalPhysica D: Nonlinear Phenomena
Volume252
DOIs
StatePublished - Jun 1 2013

Fingerprint

dynamical systems
closures
nonlinear systems
preserving
damping
energy transfer
nonlinearity
methodology
moments
formulations
symmetry
simulation
energy

Keywords

  • Modeling of nonlinear energy fluxes
  • Quasilinear Gaussian closure
  • Turbulent systems
  • Uncertainty quantification

ASJC Scopus subject areas

  • Condensed Matter Physics
  • Statistical and Nonlinear Physics

Cite this

A statistically accurate modified quasilinear Gaussian closure for uncertainty quantification in turbulent dynamical systems. / Sapsis, Themistoklis P.; Majda, Andrew J.

In: Physica D: Nonlinear Phenomena, Vol. 252, 01.06.2013, p. 34-45.

Research output: Contribution to journalArticle

@article{830f1af3023a4fe2a2aed8e82220acab,
title = "A statistically accurate modified quasilinear Gaussian closure for uncertainty quantification in turbulent dynamical systems",
abstract = "We develop a novel second-order closure methodology for uncertainty quantification in damped forced nonlinear systems with high dimensional phase-space that possess a high-dimensional chaotic attractor. We focus on turbulent systems with quadratic nonlinearities where the finite size of the attractor is caused exclusively by the synergistic activity of persistent, linearly unstable directions and a nonlinear energy transfer mechanism. We first illustrate how existing UQ schemes that rely on the Gaussian assumption will fail to perform reliable UQ in the presence of unstable dynamics. To overcome these difficulties, a modified quasilinear Gaussian (MQG) closure is developed in two stages. First we exploit exact statistical relations between second order correlations and third order moments in statistical equilibrium in order to decompose the energy flux at equilibrium into precise additional damping and enhanced noise on suitable modes, while preserving statistical symmetries; in the second stage, we develop a nonlinear MQG dynamical closure which has this statistical equilibrium behavior as a stable fixed point of the dynamics. Our analysis, UQ schemes, and conclusions are illustrated through a specific toy-model, the forty-modes Lorenz 96 system, which despite its simple formulation, presents strongly turbulent behavior with a large number of unstable dynamical components in a variety of chaotic regimes. A suitable version of MQG successfully captures the mean and variance, in transient dynamics with initial data far from equilibrium and with large random fluctuations in forcing, very cheaply at the cost of roughly two ensemble members in a Monte-Carlo simulation.",
keywords = "Modeling of nonlinear energy fluxes, Quasilinear Gaussian closure, Turbulent systems, Uncertainty quantification",
author = "Sapsis, {Themistoklis P.} and Majda, {Andrew J.}",
year = "2013",
month = "6",
day = "1",
doi = "10.1016/j.physd.2013.02.009",
language = "English (US)",
volume = "252",
pages = "34--45",
journal = "Physica D: Nonlinear Phenomena",
issn = "0167-2789",
publisher = "Elsevier",

}

TY - JOUR

T1 - A statistically accurate modified quasilinear Gaussian closure for uncertainty quantification in turbulent dynamical systems

AU - Sapsis, Themistoklis P.

AU - Majda, Andrew J.

PY - 2013/6/1

Y1 - 2013/6/1

N2 - We develop a novel second-order closure methodology for uncertainty quantification in damped forced nonlinear systems with high dimensional phase-space that possess a high-dimensional chaotic attractor. We focus on turbulent systems with quadratic nonlinearities where the finite size of the attractor is caused exclusively by the synergistic activity of persistent, linearly unstable directions and a nonlinear energy transfer mechanism. We first illustrate how existing UQ schemes that rely on the Gaussian assumption will fail to perform reliable UQ in the presence of unstable dynamics. To overcome these difficulties, a modified quasilinear Gaussian (MQG) closure is developed in two stages. First we exploit exact statistical relations between second order correlations and third order moments in statistical equilibrium in order to decompose the energy flux at equilibrium into precise additional damping and enhanced noise on suitable modes, while preserving statistical symmetries; in the second stage, we develop a nonlinear MQG dynamical closure which has this statistical equilibrium behavior as a stable fixed point of the dynamics. Our analysis, UQ schemes, and conclusions are illustrated through a specific toy-model, the forty-modes Lorenz 96 system, which despite its simple formulation, presents strongly turbulent behavior with a large number of unstable dynamical components in a variety of chaotic regimes. A suitable version of MQG successfully captures the mean and variance, in transient dynamics with initial data far from equilibrium and with large random fluctuations in forcing, very cheaply at the cost of roughly two ensemble members in a Monte-Carlo simulation.

AB - We develop a novel second-order closure methodology for uncertainty quantification in damped forced nonlinear systems with high dimensional phase-space that possess a high-dimensional chaotic attractor. We focus on turbulent systems with quadratic nonlinearities where the finite size of the attractor is caused exclusively by the synergistic activity of persistent, linearly unstable directions and a nonlinear energy transfer mechanism. We first illustrate how existing UQ schemes that rely on the Gaussian assumption will fail to perform reliable UQ in the presence of unstable dynamics. To overcome these difficulties, a modified quasilinear Gaussian (MQG) closure is developed in two stages. First we exploit exact statistical relations between second order correlations and third order moments in statistical equilibrium in order to decompose the energy flux at equilibrium into precise additional damping and enhanced noise on suitable modes, while preserving statistical symmetries; in the second stage, we develop a nonlinear MQG dynamical closure which has this statistical equilibrium behavior as a stable fixed point of the dynamics. Our analysis, UQ schemes, and conclusions are illustrated through a specific toy-model, the forty-modes Lorenz 96 system, which despite its simple formulation, presents strongly turbulent behavior with a large number of unstable dynamical components in a variety of chaotic regimes. A suitable version of MQG successfully captures the mean and variance, in transient dynamics with initial data far from equilibrium and with large random fluctuations in forcing, very cheaply at the cost of roughly two ensemble members in a Monte-Carlo simulation.

KW - Modeling of nonlinear energy fluxes

KW - Quasilinear Gaussian closure

KW - Turbulent systems

KW - Uncertainty quantification

UR - http://www.scopus.com/inward/record.url?scp=84875830966&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84875830966&partnerID=8YFLogxK

U2 - 10.1016/j.physd.2013.02.009

DO - 10.1016/j.physd.2013.02.009

M3 - Article

AN - SCOPUS:84875830966

VL - 252

SP - 34

EP - 45

JO - Physica D: Nonlinear Phenomena

JF - Physica D: Nonlinear Phenomena

SN - 0167-2789

ER -