Stochastic superparameterization in a one-dimensional model for wave turbulence

Ian G. Grooms, Andrew J. Majda

Research output: Contribution to journalArticle

Abstract

Superparameterization is a multiscale numerical method wherein solutions of prognostic equations for small scale processes on local domains embedded within the computational grid of a large scale model are computed and used to force the large scales. It was developed initially in the atmospheric sciences, but stands on its own as a nascent numerical method for the simulation of multiscale phenomena. Here we develop a stochastic version of superparameterization in a difficult one dimensional test problem involving self-similarly collapsing solitons, dispersive waves, and an intermittent inverse cascade of energy from small to large scales. We derive the nonlinear model equations by imposing a formal scale separation between resolved large scales and unresolved small scales; this allows the use of subdomains embedded within the large scale grid to describe the local small scale processes. To decrease the computational cost, we make a systematic quasi-linear stochastic approximation of the nonlinear small scale equations and use the statistical mean of the nonlinear small scale forcing (the covariance) in the large scale equations. The stochastic approximation allows the embedded domains to be formally infinite (unrealistically large scales are suppressed on the embedded domains). Further simplifications allow us to precompute the small scale forcing terms in the large scale equations as functions of the large scale variables only, which results in significant computational savings. The results are positive. The method increases the energy in overdamped simulations, decreases the energy in underdamped simulations, and improves the spatial distribution and frequency of collapsing solitons.

Original languageEnglish (US)
Pages (from-to)509-525
Number of pages17
JournalCommunications in Mathematical Sciences
Volume12
Issue number3
DOIs
StatePublished - 2014

Fingerprint

One-dimensional Model
Solitons
Turbulence
Numerical methods
Spatial distribution
Stochastic Approximation
Collapsing
Costs
Energy
Numerical Methods
Decrease
Simulation
Multiscale Methods
Forcing Term
Computational Grid
Linear Approximation
Spatial Distribution
Simplification
Forcing
Test Problems

Keywords

  • Intermittent inverse cascade
  • No scale separation
  • Stochastic modeling

ASJC Scopus subject areas

  • Mathematics(all)
  • Applied Mathematics

Cite this

Stochastic superparameterization in a one-dimensional model for wave turbulence. / Grooms, Ian G.; Majda, Andrew J.

In: Communications in Mathematical Sciences, Vol. 12, No. 3, 2014, p. 509-525.

Research output: Contribution to journalArticle

Grooms, Ian G. ; Majda, Andrew J. / Stochastic superparameterization in a one-dimensional model for wave turbulence. In: Communications in Mathematical Sciences. 2014 ; Vol. 12, No. 3. pp. 509-525.
@article{b5c15773992f450fbac8753f5fa3615d,
title = "Stochastic superparameterization in a one-dimensional model for wave turbulence",
abstract = "Superparameterization is a multiscale numerical method wherein solutions of prognostic equations for small scale processes on local domains embedded within the computational grid of a large scale model are computed and used to force the large scales. It was developed initially in the atmospheric sciences, but stands on its own as a nascent numerical method for the simulation of multiscale phenomena. Here we develop a stochastic version of superparameterization in a difficult one dimensional test problem involving self-similarly collapsing solitons, dispersive waves, and an intermittent inverse cascade of energy from small to large scales. We derive the nonlinear model equations by imposing a formal scale separation between resolved large scales and unresolved small scales; this allows the use of subdomains embedded within the large scale grid to describe the local small scale processes. To decrease the computational cost, we make a systematic quasi-linear stochastic approximation of the nonlinear small scale equations and use the statistical mean of the nonlinear small scale forcing (the covariance) in the large scale equations. The stochastic approximation allows the embedded domains to be formally infinite (unrealistically large scales are suppressed on the embedded domains). Further simplifications allow us to precompute the small scale forcing terms in the large scale equations as functions of the large scale variables only, which results in significant computational savings. The results are positive. The method increases the energy in overdamped simulations, decreases the energy in underdamped simulations, and improves the spatial distribution and frequency of collapsing solitons.",
keywords = "Intermittent inverse cascade, No scale separation, Stochastic modeling",
author = "Grooms, {Ian G.} and Majda, {Andrew J.}",
year = "2014",
doi = "10.4310/CMS.2014.v12.n3.a6",
language = "English (US)",
volume = "12",
pages = "509--525",
journal = "Communications in Mathematical Sciences",
issn = "1539-6746",
publisher = "International Press of Boston, Inc.",
number = "3",

}

TY - JOUR

T1 - Stochastic superparameterization in a one-dimensional model for wave turbulence

AU - Grooms, Ian G.

AU - Majda, Andrew J.

PY - 2014

Y1 - 2014

N2 - Superparameterization is a multiscale numerical method wherein solutions of prognostic equations for small scale processes on local domains embedded within the computational grid of a large scale model are computed and used to force the large scales. It was developed initially in the atmospheric sciences, but stands on its own as a nascent numerical method for the simulation of multiscale phenomena. Here we develop a stochastic version of superparameterization in a difficult one dimensional test problem involving self-similarly collapsing solitons, dispersive waves, and an intermittent inverse cascade of energy from small to large scales. We derive the nonlinear model equations by imposing a formal scale separation between resolved large scales and unresolved small scales; this allows the use of subdomains embedded within the large scale grid to describe the local small scale processes. To decrease the computational cost, we make a systematic quasi-linear stochastic approximation of the nonlinear small scale equations and use the statistical mean of the nonlinear small scale forcing (the covariance) in the large scale equations. The stochastic approximation allows the embedded domains to be formally infinite (unrealistically large scales are suppressed on the embedded domains). Further simplifications allow us to precompute the small scale forcing terms in the large scale equations as functions of the large scale variables only, which results in significant computational savings. The results are positive. The method increases the energy in overdamped simulations, decreases the energy in underdamped simulations, and improves the spatial distribution and frequency of collapsing solitons.

AB - Superparameterization is a multiscale numerical method wherein solutions of prognostic equations for small scale processes on local domains embedded within the computational grid of a large scale model are computed and used to force the large scales. It was developed initially in the atmospheric sciences, but stands on its own as a nascent numerical method for the simulation of multiscale phenomena. Here we develop a stochastic version of superparameterization in a difficult one dimensional test problem involving self-similarly collapsing solitons, dispersive waves, and an intermittent inverse cascade of energy from small to large scales. We derive the nonlinear model equations by imposing a formal scale separation between resolved large scales and unresolved small scales; this allows the use of subdomains embedded within the large scale grid to describe the local small scale processes. To decrease the computational cost, we make a systematic quasi-linear stochastic approximation of the nonlinear small scale equations and use the statistical mean of the nonlinear small scale forcing (the covariance) in the large scale equations. The stochastic approximation allows the embedded domains to be formally infinite (unrealistically large scales are suppressed on the embedded domains). Further simplifications allow us to precompute the small scale forcing terms in the large scale equations as functions of the large scale variables only, which results in significant computational savings. The results are positive. The method increases the energy in overdamped simulations, decreases the energy in underdamped simulations, and improves the spatial distribution and frequency of collapsing solitons.

KW - Intermittent inverse cascade

KW - No scale separation

KW - Stochastic modeling

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

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

U2 - 10.4310/CMS.2014.v12.n3.a6

DO - 10.4310/CMS.2014.v12.n3.a6

M3 - Article

VL - 12

SP - 509

EP - 525

JO - Communications in Mathematical Sciences

JF - Communications in Mathematical Sciences

SN - 1539-6746

IS - 3

ER -