RESEARCH IN RESIDENCE
Dispersion correction for the convected Helmholtz equation
Correction de dispersion pour l’équation de Helmholtz convectée
20 – 24 January 2025
Participants
Pierre-Henri Coquet (Université de Pau et des Pays de l’Adour)
Martin Gander (Université de Genève)
Antoine Tonnoir (INSA Rouen Normandie)
Solving time-harmonic wave equations with iterative methods is very challenging because they combine multiple difficulties like e.g. highly oscillatory solutions and complex non-Hermitian and indefinite discrete operators. In addition, the wave number associated to the numerical solution is different from the one of the exact solution thus leading to so-called numerical dispersion which requires very fine mesh to be controlled.
There exists numerical schemes designed to reduce the dispersion error where the most used technique to build them relies on stencils involving free parameters that are next obtained by numerically minimizing the dispersion error. In the past few years, we develop an alternative approach called asymptotic dispersion correction. The latter is based on the introduction of a perturbation of the wavenumber (the shift) in the numerical scheme. We then explicitly compute the shift that minimizes the dispersion error for small enough meshsize. In addition, we show that our method can be applied to any finite difference scheme and reduces the relative error.
We emphasize that most of works dealing with dispersion reducing methods consider the Helmholtz equation. The goal of the Research in Residence is to extend the asymptotic dispersion correction to the convected Helmholtz equation that can be used to model time-harmonic wave propagation in a moving flow.
La résolution numérique de problèmes de propagation d’ondes en régime harmonique est difficile du fait de la présence de solutions fortement oscillantes et d’opérateurs à la fois complexes, symétriques, mais non-Hermitiens, et non-coercifs. De plus, l’approximation numérique de ces problèmes souffre de l’erreur de dispersion qui provient du fait que le nombre d’onde associé à la solution numérique est différent du nombre d’onde continu. Cette dernière nécessite des maillages extrêmement fins pour être contrôlée.
Il existe plusieurs schémas numériques proposant une correction de dispersion. La plupart sont dérivés en considérant un schéma contenant des paramètres libres qui sont ensuite obtenus en minimisant numériquement l’erreur de dispersion. Récemment, nous avons développé une autre approche, appelée correction de dispersion asymptotique, qui est basée sur l’introduction d’une perturbation, appelé le shift, du nombre d’onde dans le schéma. Le shift est ensuite déterminé explicitement en minimisant l’erreur de dispersion pour des maillages suffisamment fins. Nous avons de plus montré que notre approche était applicable pour tout schémas différences finies et qu’elle permettait de réduire l’erreur relative.
Actuellement, la plupart des travaux proposant des schémas numériques avec correction de dispersion concernent l’équation de Helmholtz. L’objectif de ce Recherche en Résidence est d’étendre la correction de dispersion asymptotique au cas des équations de Helmholtz convectée qui permettent de modéliser la propagation d’ondes en régime harmonique en présence d’un écoulement support.