A Wavelet Monte Carlo Method for Turbulent Diffusion with Many Spatial Scales

Frank W. Elliott, Andrew J. Majda

Research output: Contribution to journalArticle

Abstract

A new algorithm is developed here for Monte Carlo simulation of turbulent diffusion with random velocity fields having a power law spatial spectrum, arbitrarily many statistical scales, long range correlations, and infrared divergence. This algorithm uses an expansion of the moving average representation of a Gaussian random field via a specific orthonormal basis which exploits the scaling behavior of the random field and yields a compact representation of the field despite infrared divergence. The orthonormal basis involves a family of wavelets due to Alpert and Rokhlin and the vanishing of a large number of moments for this basis (order 3 or 4) is crucial for compact representation of the statistical velocity field. The authors also develop a rigorous practical a priori energy criterion for truncation of the wavelet expansion; furthermore, for moment cancellations with orders 3 or 4 and a fixed energy tolerance, the number of Gaussian random basis elements needed by the method is sublinear in the exponent m, where 1 <|x| <2m denotes the range of active velocity scales. The algorithm developed here has remarkably low variance as regards sampling errors. For example, for the infrared divergent velocity spectrum corresponding to the Kolmogoroff spectrum in an exactly solvable model, the velocity structure function is simulated accurately for separation distances over more than five decades with only 792 random basis elements and 100 ensemble realizations; with this same singular spectrum in the model problem, the pair dispersion statistic for the passive scalar can be determined within errors of 3% throughout over five decades of pair separation distance with only 792 random elements, 1000 realizations, and a few hours calculation on a workstation. All of the computational algorithms and results are presented for a nontrivial exactly solvable model and are described within a mathematical framework of error analysis where analytic, sampling, and discretization errors are treated separately.

Original languageEnglish (US)
Pages (from-to)82-111
Number of pages30
JournalJournal of Computational Physics
Volume113
Issue number1
DOIs
StatePublished - Jul 1994

Fingerprint

turbulent diffusion
Monte Carlo method
Monte Carlo methods
divergence
Infrared radiation
velocity distribution
sampling
moments
expansion
workstations
Sampling
error analysis
cancellation
Error analysis
statistics
exponents
scalars
scaling
energy
Statistics

ASJC Scopus subject areas

  • Physics and Astronomy(all)
  • Computer Science Applications

Cite this

A Wavelet Monte Carlo Method for Turbulent Diffusion with Many Spatial Scales. / Elliott, Frank W.; Majda, Andrew J.

In: Journal of Computational Physics, Vol. 113, No. 1, 07.1994, p. 82-111.

Research output: Contribution to journalArticle

Elliott, Frank W. ; Majda, Andrew J. / A Wavelet Monte Carlo Method for Turbulent Diffusion with Many Spatial Scales. In: Journal of Computational Physics. 1994 ; Vol. 113, No. 1. pp. 82-111.
@article{b975c49fe34343fe935d0af636829184,
title = "A Wavelet Monte Carlo Method for Turbulent Diffusion with Many Spatial Scales",
abstract = "A new algorithm is developed here for Monte Carlo simulation of turbulent diffusion with random velocity fields having a power law spatial spectrum, arbitrarily many statistical scales, long range correlations, and infrared divergence. This algorithm uses an expansion of the moving average representation of a Gaussian random field via a specific orthonormal basis which exploits the scaling behavior of the random field and yields a compact representation of the field despite infrared divergence. The orthonormal basis involves a family of wavelets due to Alpert and Rokhlin and the vanishing of a large number of moments for this basis (order 3 or 4) is crucial for compact representation of the statistical velocity field. The authors also develop a rigorous practical a priori energy criterion for truncation of the wavelet expansion; furthermore, for moment cancellations with orders 3 or 4 and a fixed energy tolerance, the number of Gaussian random basis elements needed by the method is sublinear in the exponent m, where 1 <|x| <2m denotes the range of active velocity scales. The algorithm developed here has remarkably low variance as regards sampling errors. For example, for the infrared divergent velocity spectrum corresponding to the Kolmogoroff spectrum in an exactly solvable model, the velocity structure function is simulated accurately for separation distances over more than five decades with only 792 random basis elements and 100 ensemble realizations; with this same singular spectrum in the model problem, the pair dispersion statistic for the passive scalar can be determined within errors of 3{\%} throughout over five decades of pair separation distance with only 792 random elements, 1000 realizations, and a few hours calculation on a workstation. All of the computational algorithms and results are presented for a nontrivial exactly solvable model and are described within a mathematical framework of error analysis where analytic, sampling, and discretization errors are treated separately.",
author = "Elliott, {Frank W.} and Majda, {Andrew J.}",
year = "1994",
month = "7",
doi = "10.1006/jcph.1994.1120",
language = "English (US)",
volume = "113",
pages = "82--111",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",
number = "1",

}

TY - JOUR

T1 - A Wavelet Monte Carlo Method for Turbulent Diffusion with Many Spatial Scales

AU - Elliott, Frank W.

AU - Majda, Andrew J.

PY - 1994/7

Y1 - 1994/7

N2 - A new algorithm is developed here for Monte Carlo simulation of turbulent diffusion with random velocity fields having a power law spatial spectrum, arbitrarily many statistical scales, long range correlations, and infrared divergence. This algorithm uses an expansion of the moving average representation of a Gaussian random field via a specific orthonormal basis which exploits the scaling behavior of the random field and yields a compact representation of the field despite infrared divergence. The orthonormal basis involves a family of wavelets due to Alpert and Rokhlin and the vanishing of a large number of moments for this basis (order 3 or 4) is crucial for compact representation of the statistical velocity field. The authors also develop a rigorous practical a priori energy criterion for truncation of the wavelet expansion; furthermore, for moment cancellations with orders 3 or 4 and a fixed energy tolerance, the number of Gaussian random basis elements needed by the method is sublinear in the exponent m, where 1 <|x| <2m denotes the range of active velocity scales. The algorithm developed here has remarkably low variance as regards sampling errors. For example, for the infrared divergent velocity spectrum corresponding to the Kolmogoroff spectrum in an exactly solvable model, the velocity structure function is simulated accurately for separation distances over more than five decades with only 792 random basis elements and 100 ensemble realizations; with this same singular spectrum in the model problem, the pair dispersion statistic for the passive scalar can be determined within errors of 3% throughout over five decades of pair separation distance with only 792 random elements, 1000 realizations, and a few hours calculation on a workstation. All of the computational algorithms and results are presented for a nontrivial exactly solvable model and are described within a mathematical framework of error analysis where analytic, sampling, and discretization errors are treated separately.

AB - A new algorithm is developed here for Monte Carlo simulation of turbulent diffusion with random velocity fields having a power law spatial spectrum, arbitrarily many statistical scales, long range correlations, and infrared divergence. This algorithm uses an expansion of the moving average representation of a Gaussian random field via a specific orthonormal basis which exploits the scaling behavior of the random field and yields a compact representation of the field despite infrared divergence. The orthonormal basis involves a family of wavelets due to Alpert and Rokhlin and the vanishing of a large number of moments for this basis (order 3 or 4) is crucial for compact representation of the statistical velocity field. The authors also develop a rigorous practical a priori energy criterion for truncation of the wavelet expansion; furthermore, for moment cancellations with orders 3 or 4 and a fixed energy tolerance, the number of Gaussian random basis elements needed by the method is sublinear in the exponent m, where 1 <|x| <2m denotes the range of active velocity scales. The algorithm developed here has remarkably low variance as regards sampling errors. For example, for the infrared divergent velocity spectrum corresponding to the Kolmogoroff spectrum in an exactly solvable model, the velocity structure function is simulated accurately for separation distances over more than five decades with only 792 random basis elements and 100 ensemble realizations; with this same singular spectrum in the model problem, the pair dispersion statistic for the passive scalar can be determined within errors of 3% throughout over five decades of pair separation distance with only 792 random elements, 1000 realizations, and a few hours calculation on a workstation. All of the computational algorithms and results are presented for a nontrivial exactly solvable model and are described within a mathematical framework of error analysis where analytic, sampling, and discretization errors are treated separately.

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

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

U2 - 10.1006/jcph.1994.1120

DO - 10.1006/jcph.1994.1120

M3 - Article

VL - 113

SP - 82

EP - 111

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

IS - 1

ER -