Hostname: page-component-77c89778f8-sh8wx Total loading time: 0 Render date: 2024-07-17T02:08:13.157Z Has data issue: false hasContentIssue false
  • English
  • Français

An Adaptive ANOVA-Based Data-Driven Stochastic Method for Elliptic PDEs with Random Coefficient

Published online by Cambridge University Press:  03 June 2015

Zhiwen Zhang ,
Xin Hu ,
Thomas Y. Hou ,
Guang Lin  and
Mike Yan
Zhiwen Zhang
Affiliation:
Computing and Mathematical Sciences, California Institute of Technology, Pasadena, CA 91125, USA
Xin Hu
Affiliation:
Computing and Mathematical Sciences, California Institute of Technology, Pasadena, CA 91125, USA
Thomas Y. Hou*
Affiliation:
Computing and Mathematical Sciences, California Institute of Technology, Pasadena, CA 91125, USA
Guang Lin*
Affiliation:
Pacific Northwest National Laboratory, Richland, WA 99352, USA
Mike Yan
Affiliation:
Computing and Mathematical Sciences, California Institute of Technology, Pasadena, CA 91125, USA
*
Corresponding author.Email:hou@cms.caltech.edu
Corresponding author.Email:Guang.Lin@pnnl.gov
Get access

Abstract

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-Loève 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 highdimensional 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.

Keywords

60H3565C3035R6062H99Analysis of variancestochastic partial differential equationsdata-driven methodsKarhunen-Loeve expansionuncertainty quantificationmodel reduction.
Type
Research Article
Information
Communications in Computational Physics , Volume 16 , Issue 2 , August 2014 , pp. 571 - 598
Copyright
Copyright © Global Science Press Limited 2014

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1]Babuska, I., Nobile, F., and Tempone, R.A stochastic collocation method for elliptic partial differential equations with random input data. SIAM J. Numer. Anal., 45(3): 1005–1034, 2007.Google Scholar
[2]Barrault, M., Maday, Y., Nguyen, N.C., and Patera, A.T.An ‘empirical interpolation’ method: application to efficient reduced basis discretization of partial differential equations. C.R. Acad. Sci. Paris, Series I, 339: 667–672,2004.Google Scholar
[3]Barthelmann, V., Novak, E., and Ritter, K.High dimensional polynomial interpolation on sparse grids. Adv. Comput. Math., 12: 273–288, 2000.Google Scholar
[4]Bieri, M. and Schwab, C.Sparse high order FEM for elliptic SPDEs. Comput. Methods Appl. Mech. Engrg., 198: 1149–1170, 2009.Google Scholar
[5]Cao, Y., Chen, Z., and Gunzburger, M.ANOVA Expansions and Efficient Sampling Methods for Parameter Dependent Nonlinear PDEs. Int. J. Numer. Analysis and Modeling, 6(2):256–273, 2009.Google Scholar
[6]Caflisch, R.E., Morokoff, W., and Owen, A.Valuation of mortgage-backed securities using brownian bridges to reduce the effective dimension. J. Comput. Finance, 1: 27–46, 1997.Google Scholar
[7]Cheng, M.L., Hou, T.Y., Yan, M. and Zhang, Z.W.A Data-driven Stochastic Method for elliptic PDEs with random coefficients. SIAM/ASA J. UQ, 1-1,452–493,2013.Google Scholar
[8]Dostert, P., Efendiev, Y., Hou, T. Y., and Luo, W., Coarse gradient Langevin algorithms for dynamic data integration and uncertainty quantification. J. Comput. Phys., 217: 123–142, 2006.Google Scholar
[9]Fisher, R.Statistical Methods for Research Workers. Oliver and Boyd, 1925.Google Scholar
[10]Foo, J. Y. and Karniadakis, G. E.Multi-element probabilistic collocation in high dimensions. J. Comput. Phys., 229: 1536–1557, 2009.Google Scholar
[11]Ganapathysubramanian, B. and Zabaras, N.Sparse grid collocation schemes for stochastic natural convection problems. J. Comput. Phys., 225(1): 652–685, 2007.Google Scholar
[12]Gao, Z. and Hesthaven, J.S.On ANOVA Expansions and Stragies for Choosing the Anchor Point. App. Math. Comp., 217(7): 3274–3285, 2009.Google Scholar
[13]Gerstner, T. and Griebel, M.Dimension-adaptive tensor-product quadrature. Computing, 71(1): 65–87, 2003.CrossRefGoogle Scholar
[14]Ghanem, R.G. and Spanos, P.D.Stochastic Finite Elements: A Spectral Approach. New York: Springer-Verlag, 1991.Google Scholar
[15]Griebel, M.Adaptive sparse grid multilevel methods for elliptic pdes based on finite differences. Computing, 61(2):151C180,1998CrossRefGoogle Scholar
[16]Griebel, M.Sparse grids and related approximation schemes for higher dimensional problems. Proceedings of the conference on Foundations of Computational Mathematics, Santander, Spain, 2005.Google Scholar
[17]Hoeffding, W.A class of statistics with asymptotically normal distributions. Annals of Math. Statist., 19:293C325,1948.Google Scholar
[18]Halko, N., Martinsson, P.G. and Tropp, J.Finding structure with randomness: probabilistic algorithms for constructing approximate matrix decompositions. SIAM Review, 53(2):217–288, 2011Google Scholar
[19]Hou, T.Y., Luo, W., Rozovskii, B., and Zhou, H.Wiener Chaos expansions and numerical solutions of randomly forced equations of fluid mechanics. J. Comput. Phys., 216(2):687–706, 2006.Google Scholar
[20]Karhunen, K.Uber lineare Methoden in der Wahrscheinlichkeitsrechnung. Ann. Acad. Sci. Fennicae. Ser. A. I. Math.-Phys. 37: 1–79, 1947.Google Scholar
[21]Lin, G., Su, C.H., and Karniadakis, G.E.Predicting shock dynamics in the presence of uncertainties. J. Comp. Phys., 217:260C276,2006.Google Scholar
[22]Lin, G. and Karniadakis, G.E.Sensitivity analysis and stochastic simulations of nonequilibrium plasma flow. Int. J. Numer. Meth. Engng, 80(6-7):738C766,2009.Google Scholar
[23]Lin, G., Su, C.H., and Karniadakis, G. E.Random roughness enhances lift in supersonic flow. Physical Review Letters, 99:104501,2007.Google Scholar
[24]Lin, G., Su, C.H., and Karniadakis, G. E.Stochastic modeling of random roughness in shock scattering problems: Theory and simulations. Comput. Methods Appl. Mech. Engrg., 197:3420C3434,2008.Google Scholar
[25]Loève, M.Probability theory. Vol. II, 4th ed. GTM. 46. Springer-Verlag. ISBN 0-387-90262-7.Google Scholar
[26]Ma, X. and Zabaras, N.An adaptive hierarchical sparse grid collocation algorithm for the solution of stochastic differential equations. J. Comput. Phys., 228(8): 3084–3113, 2009.CrossRefGoogle Scholar
[27]Ma, X. and Zabaras, N.An adaptive high-dimensional stochastic model representation technique for the solution of stochastic partial differential equations. J. Comp. Phys., 229:3884C3915,2010.Google Scholar
[28]Le Maitre, O.Uncertainty propagation using Wiener-Haar expansions. J. Comput. Phys., 197(1): 28–57, 2004.Google Scholar
[29]Najm, H.N.Uncertainty quantification and polynomial chaos techniques in computational fluid dynamics. Ann. Rev. Fluid Mech., 41(1): 35–52, 2009.Google Scholar
[30]Novak, E. and Ritter, K.High dimensional integration of smooth funcions over cubes. Numer. Math., 75:79C97,1996.Google Scholar
[31]Novak, E. and Ritter, K.Simple cubature formulas with high polynomial exactness. Constr. Approx., 15:499C522,1999.Google Scholar
[32]Oksendal, B. K.Stochastic Differential Equations: An Introduction with Applications. Sixth edition. Berlin: Springer. 2003.Google Scholar
[33]Rabitz, H. and Alis, O. F.General foundations of high-dimensional model representations. J. Math. Chem., 25:197C233,1999.Google Scholar
[34]Tatang, M. and McRae, G.Direct treatment of uncertainty in models of reaction and transport. Technical report, MIT Tech. Rep., 1994.Google Scholar
[35]Wan, X. and Karniadakis, G.E.An adaptive multi-element generalized polynomial chaos method for stochastic differential equations. J. Comput. Phys., 209(2): 617–642, 2005.Google Scholar
[36]Winter, C.L., Guadagnini, A., Nychka, D., and Tartakovsky, D.M.Multivariate sensitivity analysis of saturated flow through simulated highly heterogeneous groundwater aquifers. J. Comput. Phys., 217:166C175,2009.Google Scholar
[37]Xiu, D. and Karniadakis, G.E.The Wiener-Askey polynomial chaos for stochastic differential equations. SIAM J. Sci. Comput., 24(2): 614–644, 2002.CrossRefGoogle Scholar
[38]Xiu, D. and Karniadakis, G.E.Modeling uncertainty in flow simulations via generalized polynomial chaos. J. Comput. Phys., 187(1): 137–167, 2003.Google Scholar
[39]Xiu, D. and Hesthaven, J.S.High-order collocation methods for differential equations with random inputs. SIAM J. Sci. Comput., 27(3): 1118–1139, 2005.Google Scholar
[40]Yang, Xiu, Choi, Minseok, Lin, Guang, and Em Karniadakis, George. Adaptive anova decomposition of stochastic incomressible and comressible flows. J. Comp. Phys., 231(4):1587–1614, 2012.Google Scholar
[41]Xu, H., Rahman, S., A generalized dimensionCreduction method for multidimensional integration in stochastic mechanics, Int. J. Numer. Methods Eng. 61 (12) (2004) 1992–2019.Google Scholar
[42]Zhang, Z., Choi, M., Karniadakis, G.E., Anchor points matter in anova decomposition, in: Spectral and High Order Methods for Partial Differential Equations Lecture Notes in Computational Science and Engineering, Springer, 76: 347–355, 2011.Google Scholar