ADAPTIVE NONINTRUSIVE RECONSTRUCTION OF SOLUTIONS TO HIGH-DIMENSIONAL PARAMETRIC PDEs.
Saved in:
| Title: | ADAPTIVE NONINTRUSIVE RECONSTRUCTION OF SOLUTIONS TO HIGH-DIMENSIONAL PARAMETRIC PDEs. |
|---|---|
| Authors: | EIGEL, MARTIN1 martin.eigel@wias-berlin.de, FARCHMIN, NANDO2 sebastian.heidenreich@ptb.de, HEIDENREICH, SEBASTIAN2 sebastian.heidenreich@ptb.de, TRUNSCHKE, PHILIPP3 philipp.trunschke@tu-berlin.de |
| Source: | SIAM Journal on Scientific Computing. 2023, Vol. 45 Issue 2, pA457-A479. 23p. |
| Subject Terms: | *STOCHASTIC analysis, *GALERKIN methods, *FINITE element method, *MONTE Carlo method, *COLLOCATION methods |
| Abstract: | Numerical methods for random parametric PDEs can greatly benefit from adaptive refinement schemes, in particular when functional approximations are computed as in stochastic Galerkin and stochastic collocations methods. This work is concerned with a nonintrusive generalization of the adaptive Galerkin finite element method with residual-based error estimation. It combines the nonintrusive character of a randomized least squares method with the a posteriori error analysis of stochastic Galerkin methods. The proposed approach uses the variational Monte Carlo method to obtain a quasi-optimal low-rank approximation of the Galerkin projection in a highly efficient hierarchical tensor format. We derive an adaptive refinement algorithm which is steered by a reliable error estimator. Opposite to stochastic Galerkin methods, the approach is easily applicable to a wide range of problems, enabling a fully automated adjustment of all discretization parameters. Benchmark examples with affine and (unbounded) lognormal coefficient fields illustrate the performance of the nonintrusive adaptive algorithm, showing the expected convergence rates of single-level strategies. [ABSTRACT FROM AUTHOR] |
| Database: | Academic Search Index |
Nájsť tento článok vo Web of Science Be the first to leave a comment!