A successive constraint linear optimization method for lower bounds of parametric coercivity and inf–sup stability constants

Abstract

We present an approach to the construction of lower bounds for the coercivity and inf–sup stability constants required in a posteriori error analysis of reduced basis approximations to affinely parametrized partial differential equations. The method, based on an Offline–Online strategy relevant in the reduced basis many-query and real-time context, reduces the Online calculation to a small Linear Program: the objective is a parametric expansion of the underlying Rayleigh quotient; the constraints reflect stability information at optimally selected parameter points. Numerical results are presented for coercive elasticity and non-coercive acoustics Helmholtz problems. To cite this article: D.B.P. Huynh et al., C. R. Acad. Sci. Paris, Ser. I 345 (2007). Nous présentons une méthode pour le calcul de bornes inférieures pour les constantes de stabilité (de coercivité ou d'inf–sup) nécessaires pour les estimateurs d'erreur a posteriori, associées à l'approximation par base réduite d'équations aux dérivées partielles ayant une dépendance affine en les paramétres. La méthode—basée sur une stratégie hors-ligne/en-ligne intéressante pour le calcul temps réel et les cas d'évaluations multiples—réduit le calcul en-ligne à un problème d'optimisation linéaire peu coûteux. La fonction objectif est un développement paramétrique du quotient de Rayleigh. Les contraintes traduisent la stabilité pour un ensemble optimal de paramétres. Nous présentons des résultats numériques pour un problème d'élasticité (coercif) ainsi que pour un problème d'acoustique de type Helmholtz (non-coercif). Pour citer cet article : D.B.P. Huynh et al., C. R. Acad. Sci. Paris, Ser. I 345 (2007).