Abstract

An algebraic subgrid scale finite element method formally equivalent to the Galerkin Least-Squares method is presented to improve the accuracy of the Galerkin finite element solution to the two-dimensional convected Helmholtz equation. A stabilizing term has been added to the discrete weak formulation containing a stabilization parameter whose value turns to be the key for the good performance of the method. An appropriate value for this parameter has been obtained by means of a dispersion analysis. As an application, we have considered the case of aerodynamic sound radiated by incompressible flow past a two-dimensional cylinder. Following Lighthill’s acoustic analogy, we have used the time Fourier transform of the double divergence of the Reynolds stress tensor as a source term for the Helmholtz and convected Helmholtz equations and showed the benefits of using the subgrid scale stabilization.

Abstract

Working with the wave equation in mixed rather than irreducible form allows one to directly account for both, the acoustic pressure field and the acoustic particle velocity field. Indeed, this becomes the natural option in many problems, such as those involving waves propagating in moving domains, because the equations can easily be set in an arbitrary Lagrangian-Eulerian (ALE) frame of reference. Yet, when attempting a standard Galerkin finite element solution (FEM) for them, it turns out that an inf-sup compatibility constraint has to be satisfied, which prevents from using equal interpolations for the approximated acoustic pressure and velocity fields. In this work it is proposed to resort to a subgrid scale stabilization strategy to circumvent this condition and thus facilitate code implementation. As a possible application, we address the generation of diphthongs in voice production.