Source author record

Rainer Hollerbach

Rainer Hollerbach appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2020arXiv

Powering Central Compact Objects with a Tangled Crustal Magnetic Field

Central Compact Objects (CCOs) are X-ray sources with luminosity ranging between $10^{32}$-$10^{34}$ erg~s$^{-1}$, located at the centres of supernova remnants. Some of them have been confirmed to be neutron stars. Timing observations have allowed the estimation of their dipole magnetic field, placing them in the range $\sim10^{10}$-$10^{11}$ G. The decay of their weak dipole fields, mediated by the Hall effect and Ohmic dissipation, cannot provide sufficient thermal energy to power their X-ray luminosity, as opposed to magnetars whose X-ray luminosities are comparable. Motivated by the question of producing high X-ray power through magnetic field decay while maintaining a weak dipole field, we explore the evolution of a crustal magnetic field that does not consist of an ordered axisymmetric structure, but rather comprises a tangled configuration. This can be the outcome of a non-self-excited dynamo, buried inside the crust by fallback material following the supernova explosion. We find that such initial conditions lead to the emergence of the magnetic field from the surface of the star and the formation of a dipolar magnetic field component. An internal tangled magnetic field of the order of $10^{14}$ G can provide sufficient Ohmic heating to the crust and power CCOs, while the dipole field it forms is approximately $10^{10}$ G, as observed in CCOs.

preprint2020arXiv

Time-dependent Probability Density Functions and Information Geometry in the Fusion L-H Transition

We report a first study of time-dependent Probability Density Functions (PDFs) in the Low-to- High confinement mode (L-H) transition by extending the previous prey-predator-type model (Kim & Diamond, Phys. Rev. Lett. 91, 185006, 2003) to a stochastic model. We highlight the limited utility of mean value and variance in understanding the L-H transition by showing strongly non-Gaussian PDFs, with the number of peaks changing in time. We also propose a new information geometric method by using information length, dynamical time scale, and information phase portrait, and show their utility in forecasting transitions and self-regulation between turbulence and zonal flows. In particular, we demonstrate the importance of intermittency (rare events of large amplitude) of zonal flows that can play an important role in promoting the L-H transition.

preprint2020arXiv

Wave propagation in semi-convective regions of giant planets

Recent observations of Jupiter and Saturn suggest that heavy elements may be diluted in the gaseous envelope, providing a compositional gradient that could stabilise ordinary convection and produce a stably-stratified layer near the core of these planets. This region could consist of semi-convective layers with a staircase-like density profile, which have multiple convective zones separated by thin stably-stratified interfaces, as a result of double-diffusive convection. These layers could have important effects on wave propagation and tidal dissipation that have not been fully explored. We analyse the effects of these layers on the propagation and transmission of internal waves within giant planets, extending prior work in a local Cartesian model. We adopt a simplified global Boussinesq planetary model in which we explore the internal waves in a non-rotating spherical body. We begin by studying the free modes of a region containing semi-convective layers. We then analyse the transmission of internal waves through such a region. The free modes depend strongly on the staircase properties, and consist of modes with both internal and interfacial gravity wave-like behaviour. We determine the frequency shifts of these waves as a function of the number of steps to explore their potential to probe planetary internal structures. We also find that wave transmission is strongly affected by the presence of a staircase. Very large-wavelength waves are transmitted efficiently, but small-scale waves are only transmitted if they are resonant with one of the free modes. The effective size of the core is therefore larger for non-resonant modes.

preprint2016arXiv

Dynamical model for spindown of solar-type stars

Since their formation, stars slow down their rotation rates by the removal of angular momentum from their surfaces, e.g. via stellar winds. Despite the complexity of the processes involved, a traditional model, where the removal of angular momentum loss by magnetic fields is prescribed, has provided a useful framework to understand observational relations between stellar rotation and age and magnetic field strength. Here, a spindown model is proposed where loss of angular momentum by magnetic fields is evolved dynamically, instead of being kinematically prescribed. To this end, we evolve the stellar rotation and magnetic field simultaneously over stellar evolution time by extending our previous work on a dynamo model which incorporates the nonlinear feedback mechanisms on rotation and magnetic fields. Our extended model reproduces key observations and explains the presence of the two branches of (fast and slow rotating) stars which have different relations between rotation rate $Ω$ vs. time (age), magnetic field strength $|B|$ vs. rotation rate, and frequency of magnetic field $ω_{cyc}$ vs. rotation rate. For fast rotating stars we find: (i) an exponential spindown $Ω\propto e^{-1.35t}$, with $t$ measured in Gyrs, (ii) $|B|$ saturates for higher rotation rate, (iii) $ω_{cyc} \propto Ω^{0.85}$. For slow rotating stars we obtain: (i) a power law spindown $Ω\propto t^{-0.52}$, (ii) $|B|$ scales almost linearly with rotation rate, (iii) $ω_{cyc} \propto Ω^{1.16}$. The results obtained are in good agreement with observations. The Vaughan-Preston gap is consistently explained in our model by the shortest spindown timescale in this transition from fast to slow rotators. Our results highlight the importance of self-regulation of magnetic fields and rotation by direct and indirect interactions involving nonlinear feedback in stellar evolution.

preprint2016arXiv

Generalised quasi-linear approximation of the HMRI

Motivated by recent advances in Direct Statistical Simulation (DSS) of astrophysical phenomena such as out of equilibrium jets, we perform a Direct Numerical Simulation (DNS) of the helical magnetorotational instability (HMRI) under the generalised quasilinear approximation (GQL). This approximation generalises the quasilinear approximation (QL) to include the self-consistent interaction of large-scale modes, interpolating between fully nonlinear DNS and QL DNS whilst still remaining formally linear in the small scales. In this paper we address whether GQL can more accurately describe low-order statistics of axisymmetric HMRI when compared with QL by performing DNS under various degrees of GQL approximation. We utilise various diagnostics, such as energy spectra in addition to first and second cumulants, for calculations performed for a range of Reynolds and Hartmann numbers (describing rotation and imposed magnetic field strength respectively). We find that GQL performs significantly better than QL in describing the statistics of the HMRI even when relatively few large-scale modes are kept in the formalism. We conclude that DSS based on GQL (GCE2) will be significantly more accurate than that based on QL (CE2).

preprint2016arXiv

Magnetic field evolution in magnetar crusts through three dimensional simulations

Current models of magnetars require extremely strong magnetic fields to explain their observed quiescent and bursting emission, implying that the field strength within the star's outer crust is orders of magnitude larger than the dipole component inferred from spin-down measurements. This presents a serious challenge to theories of magnetic field generation in a proto-neutron star. Here, we present detailed modelling of the evolution of the magnetic field in the crust of a neutron star through 3-D simulations. We find that, in the plausible scenario of equipartition of energy between global-scale poloidal and toroidal magnetic components, magnetic instabilities transfer energy to non-axisymmetric, kilometre-sized magnetic features, in which the local field strength can greatly exceed that of the global-scale field. These intense small-scale magnetic features can induce high energy bursts through local crust yielding, and the localised enhancement of Ohmic heating can power the star's persistent emission. Thus, the observed diversity in magnetar behaviour can be explained with mixed poloidal-toroidal fields of comparable energies.

preprint2016arXiv

Resistive Tearing Instability in Electron-MHD: Application to Neutron Star Crusts

We study a resistive tearing instability developing in a system evolving through the combined effect of Hall drift in the Electron-MHD limit and Ohmic dissipation. We explore first the exponential growth of the instability in the linear case and we find the fastest growing mode, the corresponding eigenvalues and dispersion relation. The instability growth rate scales as $γ\propto B^{2/3} σ^{-1/3}$ where $B$ is the magnetic field and $σ$ the electrical conductivity. We confirm the development of the tearing resistive instability in the fully non-linear case, in a plane parallel configuration where the magnetic field polarity reverses, through simulations of systems initiating in Hall equilibrium with some superimposed perturbation. Following a transient phase, during which there is some minor rearrangement of the magnetic field, the perturbation grows exponentially. Once the instability is fully developed the magnetic field forms the characteristic islands and X-type reconnection points, where Ohmic decay is enhanced. We discuss the implications of this instability for the local magnetic field evolution in neutron stars' crusts, proposing that it can contribute to heating near the surface of the star, as suggested by models of magnetar post-burst cooling. In particular, we find that a current sheet a few meters thick, covering as little as $1\%$ of the total surface can provide $10^{42}~$erg in thermal energy within a few days. We briefly discuss applications of this instability in other systems where the Hall effect operates such as protoplanetary discs and space plasmas.

preprint2016arXiv

Suppression of a kinematic dynamo by large shear

We numerically solve the magnetic induction equation in a spherical shell geometry, with a kinematically prescribed axisymmetric flow that consists of a superposition of a small-scale helical flow and a large-scale shear flow. The small-scale flow is chosen to be a local analog of the classical Roberts cells, consisting of strongly helical vortex rolls. The large-scale flow is a shearing motion in either the radial or the latitudinal directions. In the absence of large-scale shear, the small-scale flow is an efficient dynamo, in agreement with previous results. Adding increasingly large shear flows strongly suppresses the dynamo efficiency, indicating that shear is not always a favourable ingredient in dynamo action.

preprint2015arXiv

Magnetar Activity via the Density-Shear Instability in Hall-MHD

We investigate the density-shear instability in Hall-MHD via numerical simulation of the full non-linear problem, in the context of magnetar activity. We confirm the development of the instability of a plane-parallel magnetic field with an appropriate intensity and electron density profile, in accordance with analytic theory. We find that the instability also appears for a monotonically decreasing electron number density and magnetic field, a plane-parallel analogue of an azimuthal or meridional magnetic field in the crust of a magnetar. The growth rate of the instability depends on the Hall properties of the field (magnetic field intensity, electron number density and the corresponding scale-heights), while being insensitive to weak resistivity. Since the Hall effect is the driving process for the evolution of the crustal magnetic field of magnetars, we argue that this instability is critical for systems containing strong meridional or azimuthal fields. We find that this process mediates the formation of localised structures with much stronger magnetic field than the average, which can lead to magnetar activity and accelerate the dissipation of the field and consequently the production of Ohmic heating. Assuming a $5\times10^{14}$G magnetic field at the base of crust, we anticipate that magnetic field as strong as $10^{15}$G will easily develop in regions of typical size of a few $10^{2}$ meters, containing magnetic energy of $10^{43}$erg, sufficient to power magnetar bursts. These active regions are more likely to appear in the magnetic equator where the tangential magnetic field is stronger.

preprint2015arXiv

Nonaxisymmetric linear instability of cylindrical magnetohydrodynamic Taylor-Couette flow

We consider the nonaxisymmetric modes of instability present in Taylor-Couette flow under the application of helical magnetic fields, mainly for magnetic Prandtl numbers close to the inductionless limit, and conduct a full examination of marginal stability in the resulting parameter space. We allow for the azimuthal magnetic field to be generated by a combination of currents in the inner cylinder and fluid itself, and introduce a parameter governing the relation between the strength of these currents. A set of governing eigenvalue equations for the nonaxisymmetric modes of instability are derived and solved by spectral collocation with Chebyshev polynomials over the relevant parameter space, with the resulting instabilities examined in detail. We find that by altering the azimuthal magnetic field profiles the azimuthal magnetorotational instability, nonaxisymmetric helical magnetorotational instability, and Tayler instability yield interesting dynamics, such as different preferred mode types, and modes with azimuthal wave number $m>1$. Finally, a comparison is given to the recent WKB analysis performed by Kirillov $et$ $al$. [Kirillov, Stefani, and Fukumoto, J. Fluid Mech. 760, 591 (2014)] and its validity in the linear regime.

preprint2015arXiv

Three-dimensional simulations of the magnetic stress in a neutron star crust

We present the first fully self-consistent three-dimensional model of a neutron star's magnetic field, generated by electric currents in the star's crust via the Hall effect. We find that the global-scale field converges to a Hall-attractor state, as seen in recent axisymmetric models, but that small-scale features in the magnetic field survive even on much longer timescales. These small-scale features propagate toward the dipole equator, where the crustal electric currents organize themselves into a strong equatorial jet. By calculating the distribution of magnetic stresses in the crust, we predict that neutron stars with fields stronger than $10^{14}$G can still be subject to starquakes more than $10^5$yr after their formation.

preprint2014arXiv

Density-shear instability in electron MHD

We discuss a novel instability in inertia-less electron magneto-hydrodynamics (EMHD), which arises from a combination of electron velocity shear and electron density gradients. The unstable modes have a lengthscale longer than the transverse density scale, and a growth-rate of the order of the inverse Hall timescale. We suggest that this density-shear instability may be of importance in magnetic reconnection regions on scales smaller than the ion skin depth, and in neutron star crusts. We demonstrate that the so-called Hall drift instability, previously argued to be relevant in neutron star crusts, is a resistive tearing instability rather than an instability of the Hall term itself. We argue that the density-shear instability is of greater significance in neutron stars than the tearing instability, because it generally has a faster growth-rate and is less sensitive to geometry and boundary conditions. We prove that, for uniform electron density, EMHD is "at least as stable" as regular, incompressible MHD, in the sense that any field configuration that is stable in MHD is also stable in EMHD. We present a connection between the density-shear instability in EMHD and the magneto-buoyancy instability in anelastic MHD.

preprint2014arXiv

Experimental evidence for non-axisymmetric magnetorotational instability in a rotating liquid metal exposed to an azimuthal magnetic field

The azimuthal version of the magnetorotational instability (MRI) is a non-axisymmetric instability of a hydrodynamically stable differentially rotating flow under the influence of a purely or predominantly azimuthal magnetic field. It may be of considerable importance for destabilizing accretion disks, and plays a central role in the concept of the MRI dynamo. We report the results of a liquid metal Taylor-Couette experiment that shows the occurrence of an azimuthal MRI in the expected range of Hartmann numbers.

preprint2014arXiv

The effects of rotation on a double-diffusive layer in a rotating spherical shell

So far, numerical studies of double-diffusive layering in turbulent convective flows have neglected the effects of rotation. We undertake a first step into that direction by investigating how Coriolis forces affect a double-diffusive layer inside a rotating spherical shell. For this purpose we have run simulations in a parameter regime where these layers are expected to form and successively increased the rate of rotation with the result that fast rotation is found to have a similar stabilising effect on the overall convective flux as an increase of the stability ratio $R_ρ$ has in a non-rotating setup. We have also studied to what extent the regimes of rotational constraints suggested by King, Stellmach, and Buffett (2013) for rotation in the case of Rayleigh-Bénard convection are influenced by double-diffusive convection: their classification could also be applicable to the case of double-diffusive convection in a spherical shell if it is extended to be also a function of the stability ratio $R_ρ$. Furthermore, we examined the ratio of saline and thermal Nusselt numbers and compared our results with models of Spruit (2013), Rosenblum et al. (2011) and Wood, Garaud, and Stellmach (2013). We find our data to be fitted best by Spruit's model. Our result that fast rotation further decreases the convective transport, which is already lowered by double-diffusive convection, could play a major role for e.g. the modeling of the interior of some rapidly rotating giant planets, as gaseous giant planets have recently been proposed to be influenced by double-diffusive convection.

preprint2014arXiv

Tidally driven dynamos in a rotating sphere

Large-scale planetary or stellar magnetic fields generated by a dynamo effect are mostly attributed to flows forced by buoyancy forces in electrically conducting fluid layers. However, these large-scale fields may also be controlled by tides, as previously suggested for the star $τ$-boo, Mars or the Early Moon. By simulating a small local patch of a rotating fluid, \cite{Barker2014} have recently shown that tides can drive small-scale dynamos by exciting a hydrodynamic instability, the so-called elliptical (or tidal) instability. By performing global magnetohydrodynamic simulations of a rotating spherical fluid body, we investigate if this instability can also drive the observed large-scale magnetic fields. We are thus interested by the dynamo threshold and the generated magnetic field in order to test if such a mechanism is relevant for planets and stars. Rather than solving the problem in a geometry deformed by tides, we consider a spherical fluid body and add a body force to mimic the tidal deformation in the bulk of the fluid. This allows us to use an efficient spectral code to solve the magnetohydrodynamic problem. We first compare the hydrodynamic results with theoretical asymptotic results, and numerical results obtained in a truely deformed ellipsoid, which confirms the presence of the elliptical instability. We then perform magnetohydrodynamic simulations, and investigate the dynamo capability of the flow. Kinematic and self-consistent dynamos are finally simulated, showing that the elliptical instability is capable of generating dipole dominated large-scale magnetic field in global simulations of a fluid rotating sphere.

preprint2013arXiv

On the necessary conditions for bursts of convection within the rapidly rotating cylindrical annulus

Zonal flows are often found in rotating convective systems. Not only are these jet-flows driven by the convection, they can also have a profound effect on the nature of the convection. In this work the cylindrical annulus geometry is exploited in order to perform nonlinear simulations seeking to produce strong zonal flows and multiple jets. The parameter regime is extended to Prandtl numbers that are not unity. Multiple jets are found to be spaced according to a Rhines scaling based on the zonal flow speed, not the convective velocity speed. Under certain conditions the nonlinear convection appears in quasi-periodic bursts. A mean field stability analysis is performed around a basic state containing both the zonal flow and the mean temperature gradient found from the nonlinear simulations. The convective growth rates are found to fluctuate with both of these mean quantities suggesting that both are necessary in order for the bursting phenomenon to occur.

preprint2012arXiv

Kinematic dynamo in spherical Couette flow

We investigate numerically kinematic dynamos driven by flow of electrically conducting fluid in the shell between two concentric differentially rotating spheres, a configuration normally referred to as spherical Couette flow. We compare between axisymmetric (2D) and fully three dimensional flows, between low and high global rotation rates, between prograde and retrograde differential rotations, between weak and strong nonlinear inertial forces, between insulating and conducting boundaries, and between two aspect ratios. The main results are as follows. Azimuthally drifting Rossby waves arising from the destabilisation of the Stewartson shear layer are crucial to dynamo action. Differential rotation and helical Rossby waves combine to contribute to the spherical Couette dynamo. At a slow global rotation rate, the direction of differential rotation plays an important role in the dynamo because of different patterns of Rossby waves in prograde and retrograde flows. At a rapid global rotation rate, stronger flow supercriticality (namely the difference between the differential rotation rate of the flow and its critical value for the onset of nonaxisymmetric instability) facilitates the onset of dynamo action. A conducting magnetic boundary condition and a larger aspect ratio both favour dynamo action.

preprint2011arXiv

The shear-Hall instability in newborn neutron stars

Aims. In the first few minutes of a newborn neutron star's life the Hall effect and differential rotation may both be important. We demonstrate that these two ingredients are sufficient for generating a 'shear-Hall instability' and for studying its excitation conditions, growth rates, and characteristic magnetic field patterns. Methods. We numerically solve the induction equation in a spherical shell, with a kinematically prescribed differential rotation profile Ω(s), where s is the cylindrical radius. The Hall term is linearized about an imposed uniform axial field. The linear stability of individual azimuthal modes, both axisymmetric and non-axisymmetric, is then investigated. Results. For the shear-Hall instability to occur, the axial field must be parallel to the rotation axis if Ω(s) decreases outward, whereas if Ω(s) increases outward it must be anti-parallel. The instability draws its energy from the differential rotation, and occurs on the short rotational timescale rather than on the much longer Hall timescale. It operates most efficiently if the Hall time is comparable to the diffusion time. Depending on the precise field strengths B0, either axisymmetric or non-axisymmetric modes may be the most unstable. Conclusions. Even if the differential rotation in newborn neutron stars is quenched within minutes, the shear-Hall instability may nevertheless amplify any seed magnetic fields by many orders of magnitude.

preprint2010arXiv

Instabilities of Shercliff and Stewartson layers in spherical Couette flow

We explore numerically the flow induced in a spherical shell by differentially rotating the inner and outer spheres. The fluid is also taken to be electrically conducting (in the low magnetic Reynolds number limit), and a magnetic field is imposed parallel to the axis of rotation. If the outer sphere is stationary, the magnetic field induces a Shercliffe layer on the tangent cylinder, the cylinder just touching the inner sphere and parallel to the field. If the magnetic field is absent, but a strong overall rotation is present, Coriolis effects induce a Stewartson layer on the tangent cylinder. The non-axisymmetric instabilities of both types of layer separately have been studied before; here we consider the two cases side by side, as well as the mixed case, and investigate how magnetic and rotational effects interact. We find that if the differential rotation and the overall rotation are in the same direction, the overall rotation may have a destabilizing influence, whereas if the differential rotation and the overall rotation are in the opposite direction, the overall rotation always has a stabilizing influence.

preprint2010arXiv

Magnetic spherical Couette flow in linear combinations of axial and dipolar fields

We present axisymmetric numerical calculations of the fluid flow induced in a spherical shell with inner sphere rotating and outer sphere stationary. A magnetic field is also imposed, consisting of particular linear combinations of axial and dipolar fields, chosen to make $B_r=0$ at either the outer sphere, or the inner, or in between. This leads to the formation of Shercliff shear layers at these particular locations. We then consider the effect of increasingly large inertial effects, and show that an outer Shercliff layer is eventually de-stabilized, an inner Shercliff layer appears to remain stable, and an in-between Shercliff layer is almost completely disrupted even before the onset of time-dependence, which does eventually occur though.

preprint2009arXiv

Helical magnetorotational instability in a Taylor-Couette flow with strongly reduced Ekman pumping

The magnetorotational instability (MRI) is thought to play a key role in the formation of stars and black holes by sustaining the turbulence in hydrodynamically stable Keplerian accretion discs. In previous experiments the MRI was observed in a liquid metal Taylor-Couette flow at moderate Reynolds numbers by applying a helical magnetic field. The observation of this helical MRI (HMRI) was interfered with a significant Ekman pumping driven by solid end-caps that confined the instability only to a part of the Taylor-Couette cell. This paper describes the observation of the HMRI in an improved Taylor-Couette setup with the Ekman pumping significantly reduced by using split end-caps. The HMRI, which now spreads over the whole height of the cell, appears much sharper and in better agreement with numerical predictions. By analyzing various parameter dependencies we conclude that the observed HMRI represents a self-sustained global instability rather than a noise-sustained convective one.

preprint2009arXiv

Non-axisymmetric Magnetorotational Instabilities in Cylindrical Taylor-Couette Flow

We study the stability of cylindrical Taylor-Couette flow in the presence of azimuthal magnetic fields, and show that one obtains non-axisymmetric magnetorotational instabilities, having azimuthal wavenumber m=1. For Omega_o/Omega_i only slightly greater than the Rayleigh value (r_i/r_o)^2, the critical Reynolds and Hartmann numbers are Re_c ~ 10^3 and Ha_c ~ 10^2, independent of the magnetic Prandtl number Pm. These values are sufficiently small that it should be possible to obtain these instabilities in the PROMISE experimental facility.

preprint2006arXiv

The traveling wave MRI in cylindrical Taylor-Couette flow: comparing wavelengths and speeds in theory and experiment

We study experimentally the flow of a liquid metal confined between differentially rotating cylinders, in the presence of externally imposed axial and azimuthal magnetic fields. For increasingly large azimuthal fields a wave-like disturbance arises, traveling along the axis of the cylinders. The wavelengths and speeds of these structures, as well as the field strengths and rotation rates at which they arise, are broadly consistent with theoretical predictions of such a traveling wave magnetorotational instability.

preprint2005arXiv

End-effects in rapidly rotating cylindrical Taylor-Couette flow

We present numerical simulations of the flow in a rapidly rotating cylindrical annulus. We show that at the rotation rates relevant to the magneto-rotational instability, the flow is strongly constrained by the Taylor-Proudman theorem. As a result, it is controlled almost entirely by the end-plates. We then consider two possible options for minimizing these end-effects, namely (i) simply taking a very long cylinder, and (ii) splitting the end-plates into a series of differentially rotating rings. Regarding option (i), we show that the cylinder would have to be hundreds of times as long as it is wide before end-effects become unimportant in the interior. Since this is clearly not feasible, we turn to option (ii), and show that in order to obtain a smooth angular velocity profile, the end-plates would have to be split into around ten rings. If the end-plates are split into fewer rings, perhaps 3-5, the angular velocity profile will not be smooth, but will instead consist of a series of Stewartson layers at the boundaries from one ring to the next. We suggest therefore that the instabilities one obtains in this system will be the familiar Kelvin-Helmholtz instabilities of these Stewartson layers, rather than the magneto-rotational instability. At best, one might hope to obtain the MRI superimposed on these Kelvin-Helmholtz modes. Any subsequent interpretation of results is thus likely to be quite complicated.