Search

We revisit the problem of solving the one-dimensional wave equation on a domain with moving boundary. In J. Math. Phys. 11, 2679 (1970), Moore introduced an interesting method to do so. As only in rare cases, a closed analytical solution is possible, one must turn to perturbative expansions of Moore's method. We investigate the then made minimal assumption for convergence of the perturbation series, namely that the boundary position should be an analytic function of time. Though, we prove here that the latter requirement is not a sufficient condition for Moore's method to converge. We then introduce a novel numerical approach based on interpolation which also works for fast boundary dynamics. In comparison with other state-of-the-art numerical methods, our method offers greater speed if the wave solution needs to be evaluated at many points in time or space, whilst preserving accuracy. We discuss two variants of our method, either based on a conformal coordinate transformation or on the method of characteristics, together with interpolation., Comment: 18 pages

A new solution strategy for quadratic eigenvalue problems, and the derivatives of the eigenvalues, is proposed, by combining the generalized reduction method with dual numbers. To demonstrate the method, we use the quadratic eigenvalue problem encountered in the semi-analytical finite element method (SAFE) as a guiding example. The SAFE method is designed to calculate the spectrum of Lamb wave phase, group and energy velocities in (visco)elastic orthotropic media, over a wide frequency range. It was found that the new approach essentially doubles the computational speed and efficiency, without sacrificing accuracy., Comment: 5 pages, 1 figure

Matched layers are commonly used in numerical simulations of wave propagation to model (semi-)infinite domains. Attenuation functions describe the damping in layers, and provide a matching of the wave impedance at the interface between the domain of interest and the absorbing region. Selecting parameters in the attenuation functions is non-trivial. In this work, an optimisation procedure for automatically calibrating matched layers is presented. The procedure is based on solving optimisation problems constrained by partial differential equations with polynomial and piecewise-constant attenuation functions. We show experimentally that, for finite element time domain simulations, piecewise-constant attenuation function are at least as efficient as quadratic attenuation functions. This observation leads us to introduce consecutive matched layers as an alternative to perfectly matched layers, which can easily be employed for problems with arbitrary geometries. Moreover, the use of consecutive matched layers leads to a reduction in computational cost compared to perfectly matched layers. Examples are presented for acoustic, elastodynamic and electromagnetic problems. Numerical simulations are performed with the libraries FEniCS/DOLFIN and dolfin-adjoint, and the computer code to reproduce all numerical examples is made freely available.