par Ahusborde, Etienne;Azaïez, Mejdi;Deville, Michel;Mund, Ernest
Référence Computer methods in applied mechanics and engineering, 196, 45-48, page (4538-4547)
Publication Publié, 2007-09
Article révisé par les pairs
Résumé : This paper describes two Legendre spectral methods for the -grad(div) eigenvalue problem in R2. The first method uses a single grid resulting from the PN ⊗ PN discretization in primal and dual variational formulations. As is well-known, this method is unstable and exhibits spectral 'pollution' effects: increased number of singular eigenvalues, and increased multiplicity of some eigenvalues belonging to the regular spectrum. Our study aims at the understanding of these effects. The second spectral method is based on a staggered grid of the PN ⊗ PN - 1 discretization. This discretization leads to a stable algorithm, free of spurious eigenmodes and with spectral convergence of the regular eigenvalues/eigenvectors towards their analytical values. In addition, divergence-free vector fields with sufficient regularity properties are spectrally projected onto the discrete kernel of -grad(div), a clear indication of the robustness of this algorithm. © 2007 Elsevier B.V. All rights reserved.