### 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| <2^{m} 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 language | English (US) |
---|---|

Pages (from-to) | 82-111 |

Number of pages | 30 |

Journal | Journal of Computational Physics |

Volume | 113 |

Issue number | 1 |

DOIs | |

State | Published - Jul 1994 |

### Fingerprint

### ASJC Scopus subject areas

- Physics and Astronomy(all)
- Computer Science Applications

### Cite this

*Journal of Computational Physics*,

*113*(1), 82-111. https://doi.org/10.1006/jcph.1994.1120

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

Research output: Contribution to journal › Article

*Journal of Computational Physics*, vol. 113, no. 1, pp. 82-111. https://doi.org/10.1006/jcph.1994.1120

}

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 -