David J. Chappell

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.

David J. Chappell

We consider a class of potential problems on a periodic half-space for the modelling of electrified oil films, which are used in the development of novel switchable liquid optical devices (diffraction gratings). A boundary integral formulation which reduces the problem to the study of the oil-air interface alone is derived and solved in a highly efficient manner using the Nyström method. The oil films encountered experimentally are typically very thin and thus an interface-only integral representation is important for avoiding the near-singularity problems associated with boundary integral methods for long slender domains. The super-algebraic convergence of the proposed methods is discussed and demonstrated via appropriate numerical experiments.