Improved successive constraint method based a posteriori error estimate for reduced basis approximation of 2D Maxwell's problem

¹ Division of Applied Mathematics, Brown University, 182 George St, Providence, RI 02912, USA. Yanlai_Chen@brown.edu; Jan.Hesthaven@Brown.edu
² Université Pierre et Marie Curie-Paris 6, UMR 7598, Laboratoire J.-L. Lions, 75005 Paris, France. maday@ann.jussieu.fr
³ Departamento de Matemática Aplicada, Universidade de Santiago de Compostela, 15782 Santiago de Compostela, Spain. jeronimo.rodriguez@usc.es

Received: 28 August 2008
Revised: 28 March 2009

Abstract


In a posteriori error analysis of reduced basis approximations to affinely parametrized partial differential equations, the construction of lower bounds for the coercivity and inf-sup stability constants is essential. In [Huynh et al., C. R. Acad. Sci. Paris Ser. I Math. 345 (2007) 473–478], the authors presented an efficient method, compatible with an off-line/on-line strategy, where the on-line computation is reduced to minimizing a linear functional under a few linear constraints. These constraints depend on nested sets of parameters obtained iteratively using a greedy algorithm. We improve here this method so that it becomes more efficient and robust due to two related properties: (i) the lower bound is obtained by a monotonic process with respect to the size of the nested sets; (ii) less eigen-problems need to be solved. This improved evaluation of the inf-sup constant is then used to consider a reduced basis approximation of a parameter dependent electromagnetic cavity problem both for the greedy construction of the elements of the basis and the subsequent validation of the reduced basis approximation. The problem we consider has resonance features for some choices of the parameters that are well captured by the methodology.