A Multi-Field Space-Time Finite Element Method for Structural Acoustics

Author(s):  
Lonny L. Thompson

Abstract A Computational Structural Acoustics (CSA) capability for solving scattering, radiation, and other problems related to the acoustics of submerged structures has been developed by employing some of the recent algorithmic trends in Computational Fluid Dynamics (CFD), namely time-discontinuous Galerkin Least-Squares finite element methods. Traditional computational methods toward simulation of acoustic radiation and scattering from submerged elastic bodies have been primarily based on frequency domain formulations. These classical time-harmonic approaches (including boundary element, finite element, and finite difference methods) have been successful for problems involving a limited range of frequencies (narrow band response) and scales (wavelengths) that are large compared to the characteristic dimensions of the elastic structure. Attempts at solving large-scale structural acoustic systems with dimensions that are much larger than the operating wavelengths and which are complex, consisting of many different components with different scales and broadband frequencies, has revealed limitations of many of the classical methods. As a result, there has been renewed interest in new innovative approaches, including time-domain approaches. This paper describes recent advances in the development of a new class of high-order accurate and unconditionally stable space-time methods for structural acoustics which employ finite element discretization of the time domain as well as the usual discretization of the spatial domain. The formulation is based on a space-time variational equation for both the acoustic fluid and elastic structure together with their interaction. Topics to be discussed include the development and implementation of higher-order accurate non-reflecting boundary conditions based on the exact impedance relation through the. Dirichlet-to-Neumann (DtN) map, and a multi-field representation for the acoustic fluid based on independent pressure and velocity potential variables. Numerical examples involving radiation and scattering of acoustic waves are presented to illustrate the high-order accuracy achieved by the new methodology for CSA.

1995 ◽  
Vol 117 (1) ◽  
pp. 145-151 ◽  
Author(s):  
L. F. Kallivokas ◽  
J. Bielak

A finite element (FE)-based procedure is presented for the solution directly in the time-domain of transient problems involving axisymmetric three-dimensional structures submerged in an infinite acoustic fluid. The central component of the procedure is a novel impedance element that is used to render the computational domain finite. The element is local in both time and space, and is completely defined by a pair of stiffness and damping matrices. It is shown that the exterior structural acoustics problem can be solved accurately and efficiently by using the same tools as those used for interior problems. The method is illustrated with numerical results for a submerged shell.


2014 ◽  
Vol 24 (08) ◽  
pp. 1495-1539 ◽  
Author(s):  
Francesco Bassi ◽  
Lorenzo Botti ◽  
Alessandro Colombo

In this work we consider agglomeration-based physical frame discontinuous Galerkin (dG) discretization as an effective way to increase the flexibility of high-order finite element methods. The mesh free concept is pursued in the following (broad) sense: the computational domain is still discretized using a mesh but the computational grid should not be a constraint for the finite element discretization. In particular the discrete space choice, its convergence properties, and even the complexity of solving the global system of equations resulting from the dG discretization should not be influenced by the grid choice. Physical frame dG discretization allows to obtain mesh-independent h-convergence rates. Thanks to mesh agglomeration, high-order accurate discretizations can be performed on arbitrarily coarse grids, without resorting to very high-order approximations of domain boundaries. Agglomeration-based h-multigrid techniques are the obvious choice to obtain fast and grid-independent solvers. These features (attractive for any mesh free discretization) are demonstrated in practice with numerical test cases.


2018 ◽  
Vol 52 (5) ◽  
pp. 2065-2082 ◽  
Author(s):  
Erik Burman ◽  
Jonathan Ish-Horowicz ◽  
Lauri Oksanen

We consider a finite element discretization for the reconstruction of the final state of the heat equation, when the initial data is unknown, but additional data is given in a sub domain in the space time. For the discretization in space we consider standard continuous affine finite element approximation, and the time derivative is discretized using a backward differentiation. We regularize the discrete system by adding a penalty on the H2-semi-norm of the initial data, scaled with the mesh-parameter. The analysis of the method uses techniques developed in E. Burman and L. Oksanen [Numer. Math. 139 (2018) 505–528], combining discrete stability of the numerical method with sharp Carleman estimates for the physical problem, to derive optimal error estimates for the approximate solution. For the natural space time energy norm, away from t = 0, the convergence is the same as for the classical problem with known initial data, but contrary to the classical case, we do not obtain faster convergence for the L2-norm at the final time.


Sign in / Sign up

Export Citation Format

Share Document