Méthode des éléments spectraux — Wikipédia

Dans la solution numérique des équations aux dérivées partielles, un sujet en mathématiques, la méthode des éléments spectraux est une formulation de la méthode des éléments finis qui utilise des degrés élevés de polynômes par morceaux comme fonctions de base.

La méthode spectrale étend la solutions de l'équation aux séries trigonométriques, permettant ainsi d'obtenir une solution d'un ordre bien plus élevé. Cette approche se fonde sur l'idée que les polynômes trigonométriques forment une base orthonormale dans [1]. La méthode des éléments spectraux quant à elle utilise des polynômes par morceaux à degré élevé comme fonctions de bases, atteignant ainsi une grande précision. De tels polynômes sont souvent des polynômes de Tchebychev orthogonaux ou des polynômes de Lagrange d'ordre très élevé sur des espaces de noeuds non uniformes. Avec cette méthode, l'erreur induite décroît de manière exponentielle lorsque le degré des polynômes approximant augmentent. Ainsi, une convergence rapide vers la solution exacte peut être obtenue avec quelques degrés de liberté sur la structure, à l'inverse de la méthode des éléments finis. En contrôle de santé intégré, la méthode des éléments finis peut être utilisée pour détecter les plus gros défauts dans une structure. En revanche, au fur et à mesure que la taille des défaults diminue, des fréquences plus élevées sont nécessaires. Pour simuler la propagation d'une onde à fréquence élevée, la méthode des éléments finis demande un maillage de l'espace très fin, ce qui accroît donc le temps de calcul. À l'inverse, la méthode des éléments spectraux permet d'obtenir une plus grande précision pour moins de degrés de liberté. La non-uniformité des nœuds contribue à rendre la matrice de masse diagonale, ce qui permet d'économiser du temps et de la mémoire et est également utile pour adopter une méthode de différence centrale. Les inconvénients de la méthode des éléments spectraux incluent la difficulté de modéliser une géométrie complexe, par rapport à la flexibilité de la méthode des éléments finis.

Bien que la méthode puisse être appliquée avec une base polynomiale orthogonale modale par morceaux, elle est le plus souvent implémentée avec une base de Lagrange de produit tensoriel nodal[2]. La méthode gagne son efficacité en plaçant les points nodaux aux points Legendre-Gauss-Lobatto et en effectuant les intégrations de la méthode Galerkin avec une quadrature de Gauss-Lobatto réduite en utilisant les mêmes nœuds. Avec cette combinaison, il en résulte des simplifications telles qu'un regroupement de masse se produit à tous les nœuds et qu'une procédure de collocation résulte aux points intérieurs.

Les applications les plus populaires de la méthode sont la dynamique des fluides computationnelle[2] et la modélisation de la propagation des ondes sismiques[3].

La méthode des éléments spectraux a été introduite dans un article de 1984 par A. T. Patera.

Notes et références

[modifier | modifier le code]
  1. Aliki D. Muradova, « The spectral method and numerical continuation algorithm for the von Kármán problem with postbuckling behaviour of solutions », Adv Comput Math, vol. 29, no 2,‎ , p. 179–206, 2008 (DOI 10.1007/s10444-007-9050-7, hdl 1885/56758 Accès libre, S2CID 46564029)
  2. a et b Karniadakis, G. and Sherwin, S.: Spectral/hp Element Methods for Computational Fluid Dynamics, Oxford Univ. Press, (2013), (ISBN 9780199671366)
  3. Komatitsch, D. and Villote, J.-P.: “The Spectral Element Method: An Efficient Tool to Simulate the Seismic Response of 2D and 3D Geologic Structures,” Bull. Seismological Soc. America, 88, 2, 368-392 (1998)

Bibliographie

[modifier | modifier le code]
  • A. T. Patera. A spectral element method for fluid dynamics - Laminar flow in a channel expansion. Journal of Computational Physics, 54:468--488, 1984.