|Auteur||Bouillard, Philippe (Philippe.Bouillard@ulb.ac.be)|
|Titre||Méthodes de contrôle de la qualité de solutions éléments finis (application à l'acoustique)|
|Département||F509 - Faculté des sciences appliquées - Construction|
|Intitulé du diplôme||Doctorat en sciences appliquées|
|Date de défense||1997-12-05|
Coyette, Jean-Pierre (Membre du jury/Committee Member)
De Wilde, Patrick W. (Membre du jury/Committee Member)
Ihlenburg, Frank (Membre du jury/Committee Member)
Ladeveze, Pierre (Membre du jury/Committee Member)
Preumont, André (Membre du jury/Committee Member)
Beaufays, Oscar (Président du jury/Committee Chair)
Warzee, Guy (Promoteur/Director)
|Mots-clés||mécanique des milieux continus, maillages adaptatifs, estimation d'erreur à postériori, acoustique, analyse numérique, méthodes de lissage, méthode des éléments finis|
|Résumé||This work is dedicated to the control of the accuracy of computational simulations of sound propagation and scattering. Assuming time-harmonic behaviour, the mathematical models are given as boundary value problems for the Helmholtz equation Delta u+k2u=0 in Oméga. A distinction is made between interior, exterior and coupled problems and this work focuses mainly on interior uncoupled problems for which the Helmholtz equation becomes singular at eigenfrequencies.
As in other application fields, error control is an important issue in acoustic computations. It is clear that the numerical parameters (mesh size h and degree of approximation p) must be adapted to the physical parameter k. The well known ‘rule of the thumb’ for the h version with linear elements is to resolve the wavelength lambda=2 pi k-1 by six elements characterising the approximability of the finite element mesh. If the numerical model is stable, the quality of the numerical solution is entirely controlled by the approximability of the finite element mesh. The situation is quite different in the presence of singularities. In that case, stability (or the lack thereof) is equally (sometimes more) important. In our application, the solutions are ‘rough’, i.e., highly oscillatory if the wavenumber is large. This is a singularity inherent to the differential operator rather than to the domain or the boundary conditions. This effect is called the k-singularity. Similarly, the discrete operator (“stiffness” matrix) becomes singular at eigenvalues of the discretised interior problem (or nearly singular at damped eigenvalues in solid-fluid interaction). This type of singularities is called the lambda-singularities. Both singularities are of global character. Without adaptive correction, their destabilizing effect generally leads to large error of the finite element results, even if the finite element mesh satisfies the ‘rule of the thumb’.
The k- and lambda-singularities are first extensively demonstrated by numerical examples. Then, two a posteriori error estimators are developed and the numerical tests show that, due to these specific phenomena of dynamo-acoustic computations, error control cannot, in general, be accomplished by just ‘transplanting’ methods that worked well in static computations. However, for low wavenumbers, it is necessary to also control the influence of the geometric (reentrants corners) or physical (discontinuities of the boundary conditions) singularities. An h-adaptive version with refinements has been implemented. These tools have been applied to two industrial examples : the GLT, a bi-mode bus from Bombardier Eurorail, and the Vertigo, a sport car from Gillet Automobiles.
As a conclusion, it is recommanded to replace the rule of the thumb by a criterion based on the control of the influence of the specific singularities of the Helmholtz operator. As this aim cannot be achieved by the a posteriori error estimators, it is suggested to minimize the influence of the singularities by modifying the formulation of the finite element method or by formulating a “meshless” method.