Here we describe our approach to reduced basis approximation and rigorous a posteriori error estimation.
Methodology: Reduced Basis Methods and A Posteriori Error Estimators
Our methodology focuses on the rapid and reliable prediction of engineering outputs associated with parametrized partial differential equations.
In particular, we consider a (say, single) "output of interest" — related to energies or forces, stresses or strains, flowrates or pressure drops, temperatures or fluxes — as a function of an "input parameter" -vector — related to geometry, physical properties, boundary conditions, or loads. The input parameter domain — the set of possible inputs — is denoted , which is a subset of .
The output of interest is a (say, linear) functional of a field variable , . Here — say displacement, velocity, or temperature — satisfies a partial differential equation parametrized with respect to . We thus arrive at an input-output statement evaluation of which requires solution of a parametrized partial differential equation.
We consider partial differential equations for which the parametric dependence is strictly or approximately affine; "affine" dependence implies that the parametrized differential operator can be expressed as a sum of products of [parameter-dependent functions] x [parameter-independent operators].
As regards "rapid,'' our method minimizes the marginal cost associated with (approximate) input-output evaluation, and is thus most useful either (a) in the real-time or interactive context, or (b) in the limit of many queries.
Engineering situations which satisfy these criteria include in-the-field robust parameter estimation (or inverse problems, or nondestructive evaluation), design and optimization, and control.
Many educational situations also satisfy these criteria — from in-class demonstrations that require extensive parameter exploration and immediate gratification to homework assignments and projects that must be completed (by many parties) rapidly on modest platforms.
As regards "reliable,'' we provide certificates of fidelity with every prediction: an estimate that rigorously bounds the error in our (rapid, approximate) input-output evaluation or field variable relative to a highly accurate (and hence very expensive) "truth" finite element solution.
In many engineering situations, the certainty provided by these error bounds is crucial. For example, in the real-time context, critical decisions must be made in the field — quickly, without recourse to extensive Offline resources — that are at least feasible and safe if not optimal.
Educational situations also demand certainty: a demonstration or project founded upon a-physical numerical artifacts is obviously anathema to the development of sound engineering principles and practices.
back to top
The essential components of our approach are threefold.
(i) Rapidly convergent global Reduced-Basis (RB) approximations — (Galerkin) projection onto a space
spanned by solution of the governing partial differential equation at (optimally) selected points in the parameter set . Typically, will be small, as we focus attention on the (smooth) low-dimensional parametrically-induced manifold of interest. The RB approximations to the field variable and output are denoted and , respectively.
Our approach is premised upon a classical Finite Element (FE) method ''truth'' approximation space of (typically very large) dimension . It is the FE truth approximation — our introduced above — upon which we build our RB approximation, and with respect to which we measure the RB error (see (ii) below).
(ii) Rigorous a posteriori error estimation procedures — relaxations of the error-residual equation that provide inexpensive yet sharp bounds for the error in the RB field-variable approximation, , and output(s) approximation, . Our error indicators are rigorous upper bounds for the error (relative to the FE truth approximation) for all and for all ; furthermore, in many cases, we can prove that the effectivity of our error estimators — the ratio of the error bound to the true error — is O(1)–O(10).
Our inexpensive error estimators also serve to construct the optimal RB samples and spaces which ensure an efficient and well-conditioned RB approximation.
(iii) Offline/Online computational procedures — decomposition stratagems which decouple the generation and projection stages of the RB approximation: very extensive (parameter-independent) pre-processing performed Offline once that then prepares the way for subsequent very inexpensive calculations performed Online for each new input-output evaluation required.
The operation count for the Online stage — in which, given a new parameter value , the RB Online Evaluator calculates the RB output and associated error bound (relative to the expensive FE truth approximation) — depends only on and the parametric complexity of the problem. The Online computational complexity and mathematical stability does not depend on , the dimension of the underlying "truth" FE approximation space; we may thus consider a highly accurate truth approximation.
(Note that for visualization of the RB field variable approximation, the complexity does scale with — roughly as — since we must recreate and render the field over the entire physical domain.)
back to top