An Analysis of New Mixed Finite Elements for the Approximation of Wave Propagation Problems

International audienceWe construct and analyze a new family of rectangular (two-dimensional) or cubic (three-dimensional) mixed finite elements for the approximation of the acoustic wave equations. The main advantage of this element is that it permits us to obtain through mass lumping an explicit scheme even in an anisotropic medium. Nonclassical error estimates are given for this new element

Mixed Finite Elements, Strong Symmetry and Mass Lumping for Elastic Waves

We present here the continuation of our work on mixed finite elements for wave propagation problems. In a previous report, we constructed and analysed a new family of quadrangular (2D) or cubic (3D) mixed finite elements, for the approximation of the scalar anisotropic wave equation. This work is extended here to the elastic wave equation, including in the case of an anisotropic medium. These new elements present the specificity to enforce the symmetry of the stress tensor in a str ong way and lead to explicit schemes (via mass lumping), after time discretization. The convergence analysis of these mixed finite elements is not straightforward: neither the standard abstract theory nor the theory we developed for the scalar case can be applied. That is why we introduce a new abstract theory which allows to get error estimates

A novel 2.5D spectral approach for studying thin-walled waveguides with fluid-acoustic interaction

This paper presents a novel formulation of two spectral elements to study guided waves in coupled problems involving thin-walled structures and fluid-acoustic enclosures. The aim of the proposed work is the development of a new efficient computational method to study problems where geometry and properties are invariant in one direction, commonly found in the analysis of guided waves. This assumption allows using a two-and-a-half dimensional (2.5D) spectral formulation in the wavenumber-frequency domain. The novelty of the proposed work is the formulation of spectral plate and fluid elements with an arbitrary order in 2.5D. A plate element based on a Reissner-Mindlin/Kirchhoff-Love mixed formulation is proposed to represent the thin-walled structure. This element uses approximation functions to overcome the difficulties to formulate elements with an arbitrary order from functions. The proposed element uses a substitute transverse shear strain field to avoid shear locking effects. Three benchmark problems are studied to check the convergence and the computational effort for different strategies. Accurate results are found with an appropriate combination of element size and order of the approximation functions allowing at least six nodes per wavelength. The effectiveness of the proposed elements is demonstrated studying the wave propagation in a water duct with a flexible side and an acoustic cavity coupled to a Helmholtz resonator.Ministerio de Economía y Competitividad BIA2013-43085-P y BIA2016-75042-C2-1-RCentro Informático Científico de Andalucía (CICA

A stabilized finite element method for the mixed wave equation in an ALE framework with application to diphthong production

The archived file is not the final published version of the article. © (2016) S. Hirzel Verlag/European Acoustics Association The definitive publisher-authenticated version is available online at http://www.ingentaconnect.com/contentone/dav/aaua/2016/00000102/00000001/art00012 Readers must contact the publisher for reprint or permission to use the material in any form.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.Peer ReviewedPostprint (author's final draft

Numerical wave propagation for the triangular P1DGP1_{DG}-P2P2 finite element pair

Inertia-gravity mode and Rossby mode dispersion properties are examined for discretisations of the linearized rotating shallow-water equations using the $P1_{DG}$-$P2$ finite element pair on arbitrary triangulations in planar geometry. A discrete Helmholtz decomposition of the functions in the velocity space based on potentials taken from the pressure space is used to provide a complete description of the numerical wave propagation for the discretised equations. In the $f$-plane case, this decomposition is used to obtain decoupled equations for the geostrophic modes, the inertia-gravity modes, and the inertial oscillations. As has been noticed previously, the geostrophic modes are steady. The Helmholtz decomposition is used to show that the resulting inertia-gravity wave equation is third-order accurate in space. In general the \pdgp finite element pair is second-order accurate, so this leads to very accurate wave propagation. It is further shown that the only spurious modes supported by this discretisation are spurious inertial oscillations which have frequency $f$, and which do not propagate. The Helmholtz decomposition also allows a simple derivation of the quasi-geostrophic limit of the discretised $P1_{DG}$-$P2$ equations in the $\beta$-plane case, resulting in a Rossby wave equation which is also third-order accurate.Comment: Revised version prior to final journal submissio

Residual-based stabilization of the finite element approximation to the acoustic perturbation equations for low Mach number aeroacoustics

This is the peer reviewed version of the following article: [Guasch, O., Sánchez-Martín, P., Pont, A., Baiges, J., and Codina, R. (2016) Residual-based stabilization of the finite element approximation to the acoustic perturbation equations for low Mach number aeroacoustics. Int. J. Numer. Meth. Fluids, 82: 839–857. doi: 10.1002/fld.4243], which has been published in final form at http://onlinelibrary.wiley.com/doi/10.1002/fld.4243/abstract. This article may be used for non-commercial purposes in accordance with Wiley Terms and Conditions for Self-Archiving.The acoustic perturbation equations (APE) are suitable to predict aerodynamic noise in the presence of a non-uniform mean flow. As for any hybrid computational aeroacoustics approach, a first computational fluid dynamics simulation is carried out from which the mean flow characteristics and acoustic sources are obtained. In a second step, the APE are solved to get the acoustic pressure and particle velocity fields. However, resorting to the finite element method (FEM) for that purpose is not straightforward. Whereas mixed finite elements satisfying an appropriate inf–sup compatibility condition can be built in the case of no mean flow, that is, for the standard wave equation in mixed form, these are difficult to implement and their good performance is yet to be checked for more complex wave operators. As a consequence, strong simplifying assumptions are usually considered when solving the APE with FEM. It is possible to avoid them by resorting to stabilized formulations. In this work, a residual-based stabilized FEM is presented for the APE at low Mach numbers, which allows one to deal with the APE convective and reaction terms in its full extent. The key of the approach resides in the design of the matrix of stabilization parameters. The performance of the formulation and the contributions of the different terms in the equations are tested for an acoustic pulse propagating in sheared-solenoidal mean flow, and for the aeolian tone generated by flow past a two-dimensional cylinder.Peer ReviewedPostprint (author's final draft

A discontinuous Galerkin method for a new class of Green-Naghdi equations on simplicial unstructured meshes

In this paper, we introduce a discontinuous Finite Element formulation on simplicial unstructured meshes for the study of free surface flows based on the fully nonlinear and weakly dispersive Green-Naghdi equations. Working with a new class of asymptotically equivalent equations, which have a simplified analytical structure, we consider a decoupling strategy: we approximate the solutions of the classical shallow water equations supplemented with a source term globally accounting for the non-hydrostatic effects and we show that this source term can be computed through the resolution of scalar elliptic second-order sub-problems. The assets of the proposed discrete formulation are: (i) the handling of arbitrary unstructured simplicial meshes, (ii) an arbitrary order of approximation in space, (iii) the exact preservation of the motionless steady states, (iv) the preservation of the water height positivity, (v) a simple way to enhance any numerical code based on the nonlinear shallow water equations. The resulting numerical model is validated through several benchmarks involving nonlinear wave transformations and run-up over complex topographies