The common methods of spectral analysis for multivariate (n-dimensional) time series, like discrete Frourier transform (FT) or Wavelet transform, are based on Fourier series to decompose discrete data into a set of trigonometric model components, e. g. amplitude and phase. Applied to discrete data with a finite range several limitations of (time discrete) FT can be observed which are caused by the orthogonality mismatch of the trigonometric basis functions on a finite interval. However, in the general situation of non-equidistant or fragmented sampling FT based methods will cause significant errors in the parameter estimation. Therefore, the classical Lomb-Scargle method (LSM), which is not based on Fourier series, was developed as a statistical tool for one dimensional data to circumvent the inconsistent and erroneous parameter estimation of FT. The present work deduces LSM for n-dimensional data sets by a redefinition of the shifting parameter \tau, to maintain orthogonality of the trigonometric basis. An analytical derivation shows, that n-D LSM extents the traditional 1D case preserving all the statistical benefits, such as the improved noise rejection. Here, we derive the parameter confidence intervals for LSM and compare it with FT. Applications with ideal test data and experimental data will illustrate and support the proposed method.

The electrical current through an incompressible, viscous and resistive liquid conductor produces an azimuthal magnetic field that becomes unstable when the corresponding Hartmann number exceeds a critical value in the order of 20. This Tayler instability, which is not only discussed as a key ingredient of a non-linear stellar dynamo model (Tayler-Spruit dynamo), but also as a limiting factor for the maximum size of large liquid metal batteries, was recently observed experimentally in a column of a liquid metal (Seilmayer et al., Phys. Rev. Lett. 108, 244501, 2012}. On the basis of an integro-differential equation approach, we have developed a fully three-dimensional numerical code, and have utilized it for the simulation of the Tayler instability at typical viscosities and resistivities of liquid metals. The resulting growth rates are in good agreement with the experimental data. We illustrate the capabilities of the code for the detailed simulation of liquid metal battery problems in realistic geometries.