Gregor Tanner

David Chappell, Gregor Tanner, Niels Sondergaard et al.

Energy distributions of high frequency linear wave fields are often modelled in terms of flow or transport equations with ray dynamics given by a Hamiltonian vector field in phase space. Applications arise in underwater and room acoustics, vibro-acoustics, seismology, electromagnetics, and quantum mechanics. Related flow problems based on general conservation laws are used, for example, in weather forecasting or molecular dynamics simulations. Solutions to these flow equations are often large scale, complex and high-dimensional, leading to formidable challenges for numerical approximation methods. This paper presents an efficient and widely applicable method, called discrete flow mapping, for solving such problems on triangulated surfaces. An application in structural dynamics - determining the vibro-acoustic response of a cast aluminium car body component - is presented.

David J. Chappell, Gregor Tanner

Intensity distributions of linear wave fields are, in the high frequency limit, often approximated in terms of flow or transport equations in phase space. Common techniques for solving the flow equations for both time dependent and stationary problems are ray tracing or level set methods. In the context of predicting the vibro-acoustic response of complex engineering structures, reduced ray tracing methods such as Statistical Energy Analysis or variants thereof have found widespread applications. Starting directly from the stationary Liouville equation, we develop a boundary element method for solving the transport equations for complex multi-component structures. The method, which is an improved version of the Dynamical Energy Analysis technique introduced recently by the authors, interpolates between standard statistical energy analysis and full ray tracing, containing both of these methods as limiting cases. We demonstrate that the method can be used to efficiently deal with complex large scale problems giving good approximations of the energy distribution when compared to exact solutions of the underlying wave equation.

David J. Chappell, Gregor Tanner, Stefano Giani

Dynamical energy analysis was recently introduced as a new method for determining the distribution of mechanical and acoustic wave energy in complex built up structures. The technique interpolates between standard statistical energy analysis and full ray tracing, containing both of these methods as limiting cases. As such the applicability of the method is wide ranging and additionally includes the numerical modelling of problems in optics and more generally of linear wave problems in electromagnetics. In this work we consider a new approach to the method with enhanced versatility, enabling three-dimensional problems to be handled in a straightforward manner. The main challenge is the high dimensionality of the problem: we determine the wave energy density both as a function of the spatial coordinate and momentum (or direction) space. The momentum variables are expressed in separable (polar) coordinates facilitating the use of products of univariate basis expansions. However this is not the case for the spatial argument and so we propose to make use of automated mesh generating routines to both localise the approximation, allowing quadrature costs to be kept moderate, and give versatility in the code for different geometric configurations.

Janis Bajars, David Chappell, Niels Sondergaard et al.

Discrete flow mapping was recently introduced as an efficient ray based method determining wave energy distributions in complex built up structures. Wave energy densities are transported along ray trajectories through polygonal mesh elements using a finite dimensional approximation of a ray transfer operator. In this way the method can be viewed as a smoothed ray tracing method defined over meshed surfaces. Many applications require the resolution of wave energy distributions in three-dimensional domains, such as in room acoustics, underwater acoustics and for electromagnetic cavity problems. In this work we extend discrete flow mapping to three-dimensional domains by propagating wave energy densities through tetrahedral meshes. The geometric simplicity of the tetrahedral mesh elements is utilised to efficiently compute the ray transfer operator using a mixture of analytic and spectrally accurate numerical integration. The important issue of how to choose a suitable basis approximation in phase space whilst maintaining a reasonable computational cost is addressed via low order local approximations on tetrahedral faces in the position coordinate and high order orthogonal polynomial expansions in momentum space.