Multilevel quasi Monte Carlo methods for elliptic PDEs with random field coefficients via fast white noise sampling
When solving partial differential equations (PDEs) with random fields as coefficients, the efficient sampling of random field realizations can be challenging. In this paper we focus on the fast sampling of Gaussian fields using quasi-random points in a finite element and multilevel quasi Monte Carlo...
| Main Authors: | , , |
|---|---|
| Format: | Journal article |
| Language: | English |
| Published: |
Society for Industrial and Applied Mathematics
2021
|
| _version_ | 1826260220202975232 |
|---|---|
| author | Croci, M Giles, M Farrell, PE |
| author_facet | Croci, M Giles, M Farrell, PE |
| author_sort | Croci, M |
| collection | OXFORD |
| description | When solving partial differential equations (PDEs) with random fields as coefficients, the efficient sampling of random field realizations can be challenging. In this paper we focus on the fast sampling of Gaussian fields using quasi-random points in a finite element and multilevel quasi Monte Carlo (MLQMC) setting. Our method uses the stochastic PDE (SPDE) approach of Lindgren et al. combined with a new fast algorithm for white noise sampling which is tailored to (ML)QMC. We express white noise as a wavelet series expansion that we divide into two parts. The first part is sampled using quasi-random points and contains a finite number of terms in order of decaying importance to ensure good quasi Monte Carlo (QMC) convergence. The second part is a correction term which is sampled using standard pseudo-random numbers. We show how the sampling of both terms can be performed in linear time and memory complexity in the number of mesh cells via a supermesh construction, yielding an overall linear cost. Furthermore, our technique can be used to enforce the MLQMC coupling even in the case of nonnested mesh hierarchies. We demonstrate the efficacy of our method with numerical experiments. |
| first_indexed | 2024-03-06T19:02:12Z |
| format | Journal article |
| id | oxford-uuid:13e9796c-c11d-409a-a1b0-9ca80a036b45 |
| institution | University of Oxford |
| language | English |
| last_indexed | 2024-03-06T19:02:12Z |
| publishDate | 2021 |
| publisher | Society for Industrial and Applied Mathematics |
| record_format | dspace |
| spelling | oxford-uuid:13e9796c-c11d-409a-a1b0-9ca80a036b452022-03-26T10:16:36ZMultilevel quasi Monte Carlo methods for elliptic PDEs with random field coefficients via fast white noise samplingJournal articlehttp://purl.org/coar/resource_type/c_dcae04bcuuid:13e9796c-c11d-409a-a1b0-9ca80a036b45EnglishSymplectic ElementsSociety for Industrial and Applied Mathematics2021Croci, MGiles, MFarrell, PEWhen solving partial differential equations (PDEs) with random fields as coefficients, the efficient sampling of random field realizations can be challenging. In this paper we focus on the fast sampling of Gaussian fields using quasi-random points in a finite element and multilevel quasi Monte Carlo (MLQMC) setting. Our method uses the stochastic PDE (SPDE) approach of Lindgren et al. combined with a new fast algorithm for white noise sampling which is tailored to (ML)QMC. We express white noise as a wavelet series expansion that we divide into two parts. The first part is sampled using quasi-random points and contains a finite number of terms in order of decaying importance to ensure good quasi Monte Carlo (QMC) convergence. The second part is a correction term which is sampled using standard pseudo-random numbers. We show how the sampling of both terms can be performed in linear time and memory complexity in the number of mesh cells via a supermesh construction, yielding an overall linear cost. Furthermore, our technique can be used to enforce the MLQMC coupling even in the case of nonnested mesh hierarchies. We demonstrate the efficacy of our method with numerical experiments. |
| spellingShingle | Croci, M Giles, M Farrell, PE Multilevel quasi Monte Carlo methods for elliptic PDEs with random field coefficients via fast white noise sampling |
| title | Multilevel quasi Monte Carlo methods for elliptic PDEs with random field coefficients via fast white noise sampling |
| title_full | Multilevel quasi Monte Carlo methods for elliptic PDEs with random field coefficients via fast white noise sampling |
| title_fullStr | Multilevel quasi Monte Carlo methods for elliptic PDEs with random field coefficients via fast white noise sampling |
| title_full_unstemmed | Multilevel quasi Monte Carlo methods for elliptic PDEs with random field coefficients via fast white noise sampling |
| title_short | Multilevel quasi Monte Carlo methods for elliptic PDEs with random field coefficients via fast white noise sampling |
| title_sort | multilevel quasi monte carlo methods for elliptic pdes with random field coefficients via fast white noise sampling |
| work_keys_str_mv | AT crocim multilevelquasimontecarlomethodsforellipticpdeswithrandomfieldcoefficientsviafastwhitenoisesampling AT gilesm multilevelquasimontecarlomethodsforellipticpdeswithrandomfieldcoefficientsviafastwhitenoisesampling AT farrellpe multilevelquasimontecarlomethodsforellipticpdeswithrandomfieldcoefficientsviafastwhitenoisesampling |