×

Analytic regularity and stochastic collocation of high-dimensional Newton iterates. (English) Zbl 1462.65054

Summary: In this paper, we introduce concepts from uncertainty quantification (UQ) and numerical analysis for the efficient evaluation of stochastic high-dimensional Newton iterates. In particular, we develop complex analytic regularity theory of the solution with respect to the random variables. This justifies the application of sparse grids for the computation of statistical measures. Convergence rates are derived and are shown to be subexponential or algebraic with respect to the number of realizations of random perturbations. Due to the accuracy of the method, sparse grids are well suited for computing low-probability events with high confidence. We apply our method to the power flow problem. Numerical experiments on the non-trivial, 39-bus New England power system model with large stochastic loads are consistent with the theoretical convergence rates. Moreover, compared with the Monte Carlo method, our approach is at least \(10^{11}\) times faster for the same accuracy.

MSC:

65H10 Numerical computation of solutions to systems of equations
60H35 Computational methods for stochastic equations (aspects of stochastic analysis)
PDFBibTeX XMLCite
Full Text: DOI arXiv

References:

[1] Ablowitz, MJ; Fokas, AS, Complex Variables : Introduction and Applications (2003), Cambridge: Cambridge University Press, Cambridge · Zbl 1088.30001
[2] Argyros, I.K.: Convergence and Applications of Newton-Type Iterations. Springer, Berlin (2008) · Zbl 1153.65057
[3] Babuska, I.; Nobile, F.; Tempone, R., A stochastic collocation method for elliptic partial differential equations with random input data, SIAM Rev., 52, 2, 317-355 (2010) · Zbl 1226.65004 · doi:10.1137/100786356
[4] Bäck, J, Nobile, F., Tamellini, L., Tempone, R.: Stochastic Spectral Galerkin and Collocation Methods for PDEs with Random Coefficients: a Numerical Comparison. In: Hesthaven, J.S., Rønquist, E.M. (eds.) Spectral and High Order Methods for Partial Differential Equations, Lecture Notes in Computational Science and Engineering, vol. 76, pp 43-62. Springer, Berlin (2011) · Zbl 1216.65004
[5] Bank, RE; Rose, DJ, Analysis of a multilevel iterative method for nonlinear finite element equations, Math. Comput., 39, 160, 453-465 (1982) · Zbl 0496.65060 · doi:10.1090/S0025-5718-1982-0669639-X
[6] Barthelmann, V.; Novak, E.; Ritter, K., High dimensional polynomial interpolation on sparse grids, Adv. Comput. Math., 12, 273-288 (2000) · Zbl 0944.41001 · doi:10.1023/A:1018977404843
[7] Beck, J.; Nobile, F.; Tamellini, L.; Tempone, R., Convergence of quasi-optimal stochastic Galerkin methods for a class of PDEs with random coefficients, Comput. Math. Appl., 67, 4, 732-751 (2014) · Zbl 1350.65003 · doi:10.1016/j.camwa.2013.03.004
[8] Bergen, A.R., Vittal, V.: Power systems analysis, 2nd edn. Pearson/Prentice Hall (2000)
[9] Castrillón-Candás, JE; Nobile, F.; Tempone, R., Analytic regularity and collocation approximation for PDEs with random domain deformations, Comput. Math. Appl., 71, 6, 1173-1197 (2016) · Zbl 1443.65316 · doi:10.1016/j.camwa.2016.01.005
[10] Fiandrino, C.: How can I do a power electric system in circuitikz? tex.stackexchange.com/questions/145197/how-can-i-do-a-power-electric-system-in-circuitikz (2013)
[11] Fishman, G.S.: Monte Carlo : concepts, algorithms, and applications. Springer series in operations research. Springer, Berlin. http://opac.inria.fr/record=b1079070, with 98 illustrations (p. de titre) (1996) · Zbl 0859.65001
[12] Gerstner, T.; Griebel, M., Dimension-adaptive tensor-product quadrature, Computing, 71, 1, 65-87 (2003) · Zbl 1030.65015 · doi:10.1007/s00607-003-0015-5
[13] Golub, GH; Van Loan, CF, Matrix Computations (1996), Baltimore: Johns Hopkins University Press, Baltimore · Zbl 0865.65009
[14] Gunning, R., Rossi, H.: Analytic Functions of Several Complex Variables. American Mathematical Society (1965) · Zbl 0141.08601
[15] Hockenberry, J.R., Lesieutre, B.C.: Evaluation of uncertainty in dynamic simulations of power system models: The probabilistic collocation method. IEEE Transactions On Power Systems 19(3) (2004)
[16] Holst, M.: The Poisson-Boltzmann Equation: Analysis and Multilevel Numerical Solution. 1St Edn, Applied Mathematics and CRPC. California Institute of Technology (1994)
[17] J Tang, F.N., Ponci, F., Monti, A.: Dimension-adaptive sparse grid interpolation for uncertainty quantification in modern power systems: Probabilistic power flow. IEEE Transactions On Power Systems 19 (2015)
[18] Klimke, A.: Sparse Grid Interpolation Toolbox – user’s guide. Tech. Rep. IANS report 2007/017, University of Stuttgart (2007)
[19] Klimke, A., Wohlmuth, B.: Algorithm 847: spinterp: Piecewise multilinear hierarchical sparse grid interpolation in MATLAB. ACM Trans. Math. Softw. 31(4) (2005) · Zbl 1136.65308
[20] Krantz, S.G.: Function Theory of Several Complex Variables. AMS Chelsea Publishing, Providence (1992) · Zbl 0776.32001
[21] Murillo-Sánchez, CE; Zimmerman, RD; Anderson, CL; Thomas, RJ, Secure planning and operations of systems with stochastic sources, energy storage, and active demand, IEEE Trans. Smart Grid, 4, 4, 2220-2229 (2013) · doi:10.1109/TSG.2013.2281001
[22] National Academies of Sciences, Engineering, and Medicine. Analytic Research Foundations for the Next-Generation Electric Grid. The National Academies Press, Washington (2016). 10.17226/21919
[23] Nobile, F.; Tempone, R., Analysis and implementation issues for the numerical approximation of parabolic equations with random coefficients, Int. J. Numer. Methods Eng., 80, 6-7, 979-1006 (2009) · Zbl 1176.76110 · doi:10.1002/nme.2656
[24] Nobile, F., Tempone, R., Webster, C.: An anisotropic sparse grid stochastic collocation method for partial differential equations with random input data. SIAM J. Numer. Anal. 46(5), 2411-2442 (2008a). 10.1137/070680540 · Zbl 1176.65007
[25] Nobile, F., Tempone, R., Webster, C.: A sparse grid stochastic collocation method for partial differential equations with random input data. SIAM J. Numer. Anal. 46(5), 2309-2345 (2008b). 10.1137/060663660 · Zbl 1176.65137
[26] Nobile, F.; Tamellini, L.; Tempone, R., Convergence of quasi-optimal sparse-grid approximation of Hilbert-space-valued functions: application to random elliptic PDEs, Numer. Math., 134, 2, 343-388 (2016) · Zbl 1351.41004 · doi:10.1007/s00211-015-0773-y
[27] Prempraneerach, P.; Hover, F.; Triantafyllou, M.; Karniadakis, G., Uncertainty quantification in simulations of power systems: Multi-element polynomial chaos methods, Reliab. Eng. Syst. Safety, 95, 6, 632-646 (2010) · doi:10.1016/j.ress.2010.01.012
[28] Roald, L., Oldewurtel, F., Krause, T., Andersson, G.: Analytical reformulation of security constrained optimal power flow with probabilistic constraints. In: 2013 IEEE Grenoble Conference, pp 1-6. 10.1109/PTC.2013.6652224(2013)
[29] Smolyak, S., Quadrature and interpolation formulas for tensor products of certain classes of functions, Soviet Math. Doklady, 4, 240-243 (1963) · Zbl 0202.39901
[30] Tamellini, L., Nobile, F.: Sparse grids matlab kit. http://csqi.epfl.ch/page-107231-en.html (Unknown Month 2009)
[31] Trefethen, L.N.: Approximation Theory and Approximation Practice (Other Titles in Applied Mathematics). Society for Industrial and Applied Mathematics, Philadelphia (2012) · Zbl 1264.41001
[32] Zimmerman, RD; Murillo-Sánchez, CE; Thomas, RJ, Matpower: Steady-state operations, planning, and analysis tools for power systems research and education, IEEE Trans. Power Syst., 26, 1, 12-19 (2011) · doi:10.1109/TPWRS.2010.2051168
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.