An alternating low-rank projection approach for partial differential equations with random inputs
Fuente:
arXiv
Gespeichert in:
| Hauptverfasser: | , |
|---|---|
| Format: | Preprint |
| Veröffentlicht: |
2024
|
| Schlagworte: | |
| Online-Zugang: | |
| Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
| _version_ | 1866914004673232896 |
|---|---|
| author | Wang, Guanjie Liao, Qifeng |
| author_facet | Wang, Guanjie Liao, Qifeng |
| contents | It is known that standard stochastic Galerkin methods face challenges when solving partial differential equations (PDEs) with random inputs. These challenges are typically attributed to the large number of required physical basis functions and stochastic basis functions. Therefore, it becomes crucial to select effective basis functions to properly reduce the dimensionality of both the physical and stochastic approximation spaces. In this study, our focus is on the stochastic Galerkin approximation associated with generalized polynomial chaos (gPC). We delve into the low-rank approximation of the quasimatrix, whose columns represent the coefficients in the gPC expansions of the solution. We conduct an investigation into the singular value decomposition (SVD) of this quasimatrix, proposing a strategy to identify the rank required for a desired accuracy. Subsequently, we introduce both a simultaneous low-rank projection approach and an alternating low-rank projection approach to compute the low-rank approximation of the solution for PDEs with random inputs. Numerical results demonstrate the efficiency of our proposed methods for both diffusion and Helmholtz problems. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2410_22183 |
| institution | arXiv |
| publishDate | 2024 |
| record_format | arxiv |
| spellingShingle | An alternating low-rank projection approach for partial differential equations with random inputs Wang, Guanjie Liao, Qifeng Numerical Analysis It is known that standard stochastic Galerkin methods face challenges when solving partial differential equations (PDEs) with random inputs. These challenges are typically attributed to the large number of required physical basis functions and stochastic basis functions. Therefore, it becomes crucial to select effective basis functions to properly reduce the dimensionality of both the physical and stochastic approximation spaces. In this study, our focus is on the stochastic Galerkin approximation associated with generalized polynomial chaos (gPC). We delve into the low-rank approximation of the quasimatrix, whose columns represent the coefficients in the gPC expansions of the solution. We conduct an investigation into the singular value decomposition (SVD) of this quasimatrix, proposing a strategy to identify the rank required for a desired accuracy. Subsequently, we introduce both a simultaneous low-rank projection approach and an alternating low-rank projection approach to compute the low-rank approximation of the solution for PDEs with random inputs. Numerical results demonstrate the efficiency of our proposed methods for both diffusion and Helmholtz problems. |
| title | An alternating low-rank projection approach for partial differential equations with random inputs |
| topic | Numerical Analysis |
| url | https://arxiv.org/abs/2410.22183 |