### Abstract

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.

Original language | English (US) |
---|---|

Pages (from-to) | 384-408 |

Number of pages | 25 |

Journal | Journal of Computational Physics |

Volume | 132 |

Issue number | 2 |

DOIs | |

State | Published - Apr 1997 |

### Fingerprint

### ASJC Scopus subject areas

- Computer Science Applications
- Physics and Astronomy(all)

### Cite this

*Journal of Computational Physics*,

*132*(2), 384-408. https://doi.org/10.1006/jcph.1996.5647

**A fourier-wavelet Monte Carlo method for fractal random fields.** / Elliott, Frank W.; Horntrop, David J.; Majda, Andrew J.

Research output: Contribution to journal › Article

*Journal of Computational Physics*, vol. 132, no. 2, pp. 384-408. https://doi.org/10.1006/jcph.1996.5647

}

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.

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

VL - 132

SP - 384

EP - 408

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

IS - 2

ER -