An adaptive ANOVA-based data-driven stochastic method for elliptic PDEs with random coefficient

In this paper, we present an adaptive, analysis of variance (ANOVA)-based data-driven stochastic method (ANOVA-DSM) to study the stochastic partial differential equations (SPDEs) in the multi-query setting. Our new method integrates the advantages of both the adaptive ANOVA decomposition technique and the data-driven stochastic method. To handle high-dimensional stochastic problems, we investigate the use of adaptive ANOVA decomposition in the stochastic space as an effective dimension-reduction technique. To improve the slow convergence of the generalized polynomial chaos (gPC) method or stochastic collocation (SC) method, we adopt the data-driven stochastic method (DSM) for speed up. An essential ingredient of the DSM is to construct a set of stochastic basis under which the stochastic solutions enjoy a compact representation for a broad range of forcing functions and/or boundary conditions. Our ANOVA-DSM consists of offline and online stages. In the offline stage, the original high-dimensional stochastic problem is decomposed into a series of low-dimensional stochastic subproblems, according to the ANOVA decomposition technique. Then, for each subproblem, a data-driven stochastic basis is computed using the Karhunen-Loeve expansion (KLE) and a two-level preconditioning optimization approach. Multiple trial functions are used to enrich the stochastic basis and improve the accuracy. In the online stage, we solve each stochastic subproblem for any given forcing function by projecting the stochastic solution into the data-driven stochastic basis constructed offline. In our ANOVA-DSM framework, solving the original high-dimensional stochastic problem is reduced to solving a series of ANOVA-decomposed stochastic subproblems using the DSM. An adaptive ANOVA strategy is also provided to further reduce the number of the stochastic subproblems and speed up our method. To demonstrate the accuracy and efficiency of our method, numerical examples are presented for one- and two-dimensional elliptic PDEs with random coefficients.

[1]  Habib N. Najm,et al.  Uncertainty Quantification and Polynomial Chaos Techniques in Computational Fluid Dynamics , 2009 .

[2]  H. Rabitz,et al.  General foundations of high‐dimensional model representations , 1999 .

[3]  Thomas Y. Hou,et al.  A Multiscale Finite Element Method for Elliptic Problems in Composite Materials and Porous Media , 1997 .

[4]  R. Ghanem,et al.  Stochastic Finite Elements: A Spectral Approach , 1990 .

[5]  Yanzhao Cao,et al.  International Journal of C 2009 Institute for Scientific Numerical Analysis and Modeling Computing and Information Anova Expansions and Efficient Sampling Methods for Parameter Dependent Nonlinear Pdes , 2022 .

[6]  W. Hoeffding.  A Class of Statistics with Asymptotically Normal Distribution , 1948 .

[7]  K. Ritter,et al.  High dimensional integration of smooth functions over cubes , 1996 .

[8]  Xiu Yang,et al.  Adaptive ANOVA decomposition of stochastic incompressible and compressible flows , 2012, J. Comput. Phys..

[9]  G. Karniadakis,et al.  An adaptive multi-element generalized polynomial chaos method for stochastic differential equations , 2005 .

[10]  Thomas Y. Hou,et al.  Wiener Chaos expansions and numerical solutions of randomly forced equations of fluid mechanics , 2006, J. Comput. Phys..

[11]  Martin Berggren,et al.  Optimization of an acoustic horn with respect to efficiency and directivity , 2008 .

[12]  George Em Karniadakis,et al.  Random roughness enhances lift in supersonic flow. , 2007, Physical review letters.

[13]  George Em Karniadakis,et al.  Sensitivity analysis and stochastic simulations of non‐equilibrium plasma flow , 2009 .

[14]  M. Kendall.  Statistical Methods for Research Workers , 1937, Nature.

[15]  Christoph Schwab,et al.  Karhunen-Loève approximation of random fields by generalized fast multipole methods , 2006, J. Comput. Phys..

[16]  George Em Karniadakis,et al.  Stochastic modeling of random roughness in shock scattering problems: Theory and simulations , 2008 .

[17]  Xiang Ma,et al.  An adaptive high-dimensional stochastic model representation technique for the solution of stochastic partial differential equations , 2010, J. Comput. Phys..

[18]  D. Xiu,et al.  Modeling uncertainty in flow simulations via generalized polynomial chaos , 2003 .

[19]  R. Cools,et al.  Monomial cubature rules since “Stroud”: a compilation , 1993 .

[20]  Baskar Ganapathysubramanian,et al.  Sparse grid collocation schemes for stochastic natural convection problems , 2007, J. Comput. Phys..

[21]  Guang Lin,et al.  An efficient, high-order probabilistic collocation method on sparse grids for three-dimensional flow and solute transport in randomly heterogeneous porous media , 2009 .

[22]  L. Sirovich.  Turbulence and the dynamics of coherent structures. I. Coherent structures , 1987 .

[23]  Daniele Venturi,et al.  Stochastic low-dimensional modelling of a random laminar wake past a circular cylinder , 2008, Journal of Fluid Mechanics.

[24]  BabuskaIvo,et al.  A Stochastic Collocation Method for Elliptic Partial Differential Equations with Random Input Data , 2007 .

[25]  Erich Novak,et al.  High dimensional polynomial interpolation on sparse grids , 2000, Adv. Comput. Math..

[26]  R. Ghanem,et al.  Uncertainty propagation using Wiener-Haar expansions , 2004 .

[27]  K. Ritter,et al.  Simple Cubature Formulas with High Polynomial Exactness , 1999 .

[28]  C. Schwab,et al.  Sparse high order FEM for elliptic sPDEs , 2009 .

[29]  Li Jun Jiang,et al.  An Adaptive Hierarchical Sparse Grid Collocation Method , 2014 .

[30]  Yalchin Efendiev,et al.  Coarse-gradient Langevin algorithms for dynamic data integration and uncertainty quantification , 2006, J. Comput. Phys..

[31]  R. Cools.  Monomial cubature rules since “Stroud”: a compilation—part 2 , 1999 .

[32]  H. Rabitz,et al.  Efficient input-output model representations , 1999 .

[33]  Alexandre M. Tartakovsky,et al.  Numerical Studies of Three-dimensional Stochastic Darcy’s Equation and Stochastic Advection-Diffusion-Dispersion Equation , 2010, J. Sci. Comput..

[34]  Daniel M. Tartakovsky,et al.  Multivariate sensitivity analysis of saturated flow through simulated highly heterogeneous groundwater aquifers , 2005, J. Comput. Phys..

[35]  Nathan Halko,et al.  Finding Structure with Randomness: Probabilistic Algorithms for Constructing Approximate Matrix Decompositions , 2009, SIAM Rev..

[36]  Dongbin Xiu,et al.  High-Order Collocation Methods for Differential Equations with Random Inputs , 2005, SIAM J. Sci. Comput..

[37]  B. Øksendal.  Stochastic differential equations : an introduction with applications , 1987 .

[38]  S. Rahman,et al.  A generalized dimension‐reduction method for multidimensional integration in stochastic mechanics , 2004 .

[39]  Michael Griebel,et al.  Adaptive sparse grid multilevel methods for elliptic PDEs based on finite differences , 1998, Computing.

[40]  N. Nguyen,et al.  An ‘empirical interpolation’ method: application to efficient reduced-basis discretization of partial differential equations , 2004 .

[41]  Dongbin Xiu,et al.  The Wiener-Askey Polynomial Chaos for Stochastic Differential Equations , 2002, SIAM J. Sci. Comput..

[42]  Thomas Gerstner,et al.  Dimension–Adaptive Tensor–Product Quadrature , 2003, Computing.

[43]  Guang Lin,et al.  Predicting shock dynamics in the presence of uncertainties , 2006, J. Comput. Phys..

[44]  Kari Karhunen,et al.  Über lineare Methoden in der Wahrscheinlichkeitsrechnung , 1947 .

[45]  A. Owen,et al.  Valuation of mortgage-backed securities using Brownian bridges to reduce effective dimension , 1997 .

[46]  George E. Karniadakis,et al.  Multi-element probabilistic collocation method in high dimensions , 2010, J. Comput. Phys..

[47]  Jasmine Yen-teng Foo,et al.  Multi-element probabilistic collocation in high dimensions: Applications to systems biology and physical systems , 2008 .

[48]  Zhiwen Zhang,et al.  A Data-Driven Stochastic Method for Elliptic PDEs with Random Coefficients , 2013, SIAM/ASA J. Uncertain. Quantification.

[49]  Xiang Ma,et al.  An adaptive hierarchical sparse grid collocation algorithm for the solution of stochastic differential equations , 2009, J. Comput. Phys..

[50]  George Em Karniadakis,et al.  Anchor Points Matter in ANOVA Decomposition , 2011 .

[51]  M. Griebel.  Sparse Grids and Related Approximation Schemes for Higher Dimensional Problems , 2006 .

[52]  Fabio Nobile,et al.  A Stochastic Collocation Method for Elliptic Partial Differential Equations with Random Input Data , 2007, SIAM Rev..

[53]  Jan S. Hesthaven,et al.  On ANOVA expansions and strategies for choosing the anchor point , 2010, Appl. Math. Comput..