Open Access

ARTICLE

Approximation of the grad div Operator in Nonconvex Domains

D. Boffi¹, C. Chinosi², L. Gastaldi³

1 Dipartimento di Matematica “F. Casorati”, Università di Pavia, Italy and Department of Mathematics, Pennsylvania State University, USA
2 Dipartimento di Scienze e Tecnologie Avanzate, Università del Piemonte Orientale, Alessandria, Italy
3 Dipartimento di Matematica, Università di Brescia, Italy

Computer Modeling in Engineering & Sciences 2000, 1(2), 31-44. https://doi.org/10.3970/cmes.2000.001.191

Abstract

In this paper we are dealing with the approximation of the grad-div operator in nonconvex polygonal domains. A penalization strategy is considered in order to obtain a formulation of the original eigenproblem which is associated with an elliptic operator. However the presence of singular eigensolutions, in the case of nonconvex domains, is the origin of major troubles in the numerical approximation of the problem. A mixed-type approximation, based on a projection procedure, is introduced and analyzed from the theoretical and numerical point of view. Several numerical experiments confirm that in presence of singularities the projection is needed in order to reproduce the features of the continuous problem.

---