TY - JOUR

T1 - A fourier-wavelet Monte Carlo method for fractal random fields

AU - Elliott, Frank W.

AU - Horntrop, David J.

AU - Majda, Andrew J.

N1 - Funding Information:
Andrew J. Majda’s research is partially supported by the following grants: National Science Foundation DMS-9596102, Army Research Office DAAH04-95-1-0345, and Office of Naval Research N00014-96-1-0043. David J. Horntrop and Frank W. Elliott Jr. were supported as postdoctoral researchers by grants from the Army Research Office (DAAH04-95-1-0345) and the National Science Foundation (DMS-9596102), respectively.

PY - 1997/4

Y1 - 1997/4

N2 - A new hierarchical method for the Monte Carlo simulation of random fields called the Fourier-wavelet method is developed and applied to isotropic Gaussian random fields with power law spectral density functions. This technique is based upon the orthogonal decomposition of the Fourier stochastic integral representation of the field using wavelets. The Meyer wavelet is used here because its rapid decay properties allow for a very compact representation of the field. The Fourier-wavelet method is shown to be straightforward to implement, given the nature of the necessary precomputations and the run-time calculations, and yields comparable results with scaling behavior over as many decades as the physical space multiwavelet methods developed recently by two of the authors. However, the Fourier-wavelet method developed here is more flexible and, in particular, applies to anisotropic spectra generated through solutions of differential equations. Simulation results using this new technique and the well-known nonhierarchical simulation technique, the randomization method, are given and compared for both a simple shear layer model problem as well as a two-dimensional isotropic Gaussian random field. The Fourier-wavelet method results are more accurate for statistical quantities depending on moments higher than order 2, in addition to showing a quite smooth decay to zero on the scales smaller than the scaling regime when compared with the randomization method results. The only situation in which the nonhierarchical randomization method is more computationally efficient occurs when no more than four decades of scaling behavior are needed and the statistical quantities of interest depend only on second moments.

AB - A new hierarchical method for the Monte Carlo simulation of random fields called the Fourier-wavelet method is developed and applied to isotropic Gaussian random fields with power law spectral density functions. This technique is based upon the orthogonal decomposition of the Fourier stochastic integral representation of the field using wavelets. The Meyer wavelet is used here because its rapid decay properties allow for a very compact representation of the field. The Fourier-wavelet method is shown to be straightforward to implement, given the nature of the necessary precomputations and the run-time calculations, and yields comparable results with scaling behavior over as many decades as the physical space multiwavelet methods developed recently by two of the authors. However, the Fourier-wavelet method developed here is more flexible and, in particular, applies to anisotropic spectra generated through solutions of differential equations. Simulation results using this new technique and the well-known nonhierarchical simulation technique, the randomization method, are given and compared for both a simple shear layer model problem as well as a two-dimensional isotropic Gaussian random field. The Fourier-wavelet method results are more accurate for statistical quantities depending on moments higher than order 2, in addition to showing a quite smooth decay to zero on the scales smaller than the scaling regime when compared with the randomization method results. The only situation in which the nonhierarchical randomization method is more computationally efficient occurs when no more than four decades of scaling behavior are needed and the statistical quantities of interest depend only on second moments.

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

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

U2 - 10.1006/jcph.1996.5647

DO - 10.1006/jcph.1996.5647

M3 - Article

AN - SCOPUS:0031116394

VL - 132

SP - 384

EP - 408

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

IS - 2

ER -