Reduced-basis output bounds for approximately parametrized elliptic coercive partial differential equations

We present a technique for the rapid and reliable prediction of linear-functional outputs of elliptic coercive partial differential equations with (approximately) affine parameter dependence. The essential components are (i) (provably) rapidly convergent global reduced-basis approximations – Galerkin projection onto a space WN spanned by solutions of the governing partial differential equation at N selected points in parameter space; (ii) a posteriori error estimation – relaxations of the error-residual equation that provide inexpensive yet sharp bounds for the error in the outputs of interest; and (iii) off-line/on-line computational procedures – methods which decouple the generation and projection stages of the approximation process. The operation count for the on-line stage – in which, given a new parameter value, we calculate the output of interest and associated error bound – depends only on N, typically very small, and the (approximate) parametric complexity of the problem; the method is thus ideally suited for the repeated and rapid evaluations required in the context of parameter estimation, design, optimization, and real-time control.In our earlier work, we develop a rigorous a posteriori error bound framework for the case in which the parametrization of the partial differential equation is exact; in this paper, we address the situation in which our mathematical model is not “complete.” In particular, we permit error in the data that define our partial differential operator: this error may be introduced, for example, by imperfect specification, measurement, calculation, or parametric expansion of a coefficient function. We develop both accurate predictions for the outputs of interest and associated rigorous a posteriori error bounds; and the latter incorporate both numerical discretization and “model truncation” effects. Numerical results are presented for a particular instantiation in which the model error originates in the (approximately) prescribed velocity field associated with a three-dimensional convection-diffusion problem.

[1]  W. Wendland.  Elliptic systems in the plane , 1979 .

[2]  Jacques-Louis Lions,et al.  Nonlinear partial differential equations and their applications: Collège de France seminar. Volume XIV , 2002 .

[3]  E. Allgower,et al.  Simplicial and Continuation Methods for Approximating Fixed Points and Solutions to Systems of Equations , 1980 .

[4]  Janet S. Peterson,et al.  The Reduced Basis Method for Incompressible Viscous Flow Calculations , 1989 .

[5]  Ahmed K. Noor,et al.  Reduced Basis Technique for Nonlinear Analysis of Structures , 1980 .

[6]  D. Rovas,et al.  Reliable Real-Time Solution of Parametrized Partial Differential Equations: Reduced-Basis Output Bound Methods , 2002 .

[7]  D. Rovas,et al.  A blackbox reduced-basis output bound method for noncoercive linear problems , 2002 .

[8]  G. Reddien,et al.  On the reduced basis method , 1995 .

[9]  W. Rheinboldt.  Numerical analysis of continuation methods for nonlinear structural problems , 1981 .

[10]  T. A. Porsching,et al.  Estimation of the error in the reduced basis method solution of nonlinear equations , 1985 .

[11]  Anthony T. Patera,et al.  A Priori Convergence Theory for Reduced-Basis Approximations of Single-Parameter Elliptic Partial Differential Equations , 2002, J. Sci. Comput..

[12]  V. Lakshmikantham,et al.  Nonlinear Analysis: Theory, Methods and Applications , 1978 .

[13]  Anthony T. Patera,et al.  Global a priori convergence theory for reduced-basis approximations of single-parameter symmetric coercive elliptic partial differential equations , 2002 .

[14]  G. Strang,et al.  An Analysis of the Finite Element Method , 1974 .

[15]  Karen Paula L. Veroy,et al.  Reduced-basis methods applied to problems in elasticity : analysis and applications , 2003 .

[16]  C. Farhat,et al.  Extending substructure based iterative solvers to multiple load and repeated analyses , 1994 .

[17]  Tony F. Chan,et al.  Analysis of Projection Methods for Solving Linear Systems with Multiple Right-Hand Sides , 1997, SIAM J. Sci. Comput..

[18]  Werner C. Rheinboldt,et al.  On the theory and error estimation of the reduced basis method for multi-parameter problems , 1993 .

[19]  E. Yip.  A Note on the Stability of Solving a Rank-p Modification of a Linear System by the Sherman–Morrison–Woodbury Formula , 1986 .

[20]  P. Stern,et al.  Automatic choice of global shape functions in structural analysis , 1978 .

[21]  D. Rovas,et al.  Output bounds for reduced-basis approximations of symmetric positive definite eigenvalue problems , 2000 .

[22]  Etienne Balmes,et al.  PARAMETRIC FAMILIES OF REDUCED FINITE ELEMENT MODELS. THEORY AND APPLICATIONS , 1996 .

[23]  S. Fischler.  Formes linéaires en polyzêtas et intégrales multiples , 2002 .

[24]  Werner C. Rheinboldt,et al.  On the Error Behavior of the Reduced Basis Technique for Nonlinear Finite Element Approximations , 1983 .

[25]  Raphael T. Haftka,et al.  Fast exact linear and non‐linear structural reanalysis and the Sherman–Morrison–Woodbury formulas , 2001 .

[26]  Anthony T. Patera,et al.  A posteriori error estimation for reduced-basis approximation of parametrized elliptic coercive partial differential equations : “convex inverse” bound conditioners , 2002 .