Linear stability of magnetohydrodynamic flow in a perfectly conducting rectangular duct

We analyse numerically the linear stability of a liquid metal flow in a rectangular duct with perfectly electrically conducting walls subject to a uniform transverse magnetic field. A non-standard three dimensional vector stream function/vorticity formulation is used with Chebyshev collocation method to solve the eigenvalue problem for small-amplitude perturbations. A relatively weak magnetic field is found to render the flow linearly unstable as two weak jets appear close to the centre of the duct at the Hartmann number Ha \approx 9.6. In a sufficiently strong magnetic field, the instability following the jets becomes confined in the layers of characteristic thickness \delta \sim Ha^{-1/2} located at the walls parallel to the magnetic field. In this case the instability is determined by \delta, which results in both the critical Reynolds and wavenumbers numbers scaling as \sim \delta^{-1}. Instability modes can have one of the four different symmetry combinations along and across the magnetic field. The most unstable is a pair of modes with an even distribution of vorticity along the magnetic field. These two modes represent strongly non-uniform vortices aligned with the magnetic field, which rotate either in the same or opposite senses across the magnetic field. The former enhance while the latter weaken one another provided that the magnetic field is not too strong or the walls parallel to the field are not too far apart. In a strong magnetic field, when the vortices at the opposite walls are well separated by the core flow, the critical Reynolds and wavenumbers for both of these instability modes are the same: Re_c \approx 642Ha^{1/2}+8.9x10^3Ha^{-1/2} and k_c \approx 0.477Ha^{1/2}. The other pair of modes, which differs from the previous one by an odd distribution of vorticity along the magnetic field, is more stable with approximately four times higher critical Reynolds number.Comment: 16 pages, 8 figures, revised version, to appear in J. Fluid Mec

Contactless Electromagnetic Phase-Shift Flowmeter for Liquid Metals

We present a concept and test results of an eddy-current flowmeter for liquid metals. The flow rate is determined by applying a weak ac magnetic field to a liquid metal flow and measuring the flow-induced phase disturbance in the external electromagnetic field. The phase disturbance is found to be more robust than that of the amplitude used in conventional eddy-current flowmeters. The basic characteristics of this type of flowmeter are analysed using simple theoretical models, where the flow is approximated by a solid body motion. Design of such a flowmeter is presented and its test results reported.Comment: 19 pages, 13 figures, to appear in Meas. Sci. Technol (final version

Elementary model of internal electromagnetic pinch-type instability

We analyse numerically a pinch-type instability in a semi-infinite planar layer of inviscid conducting liquid bounded by solid walls and carrying a uniform electric current. Our model is as simple as possible but still captures the salient features of the instability which otherwise may be obscured by the technical details of more comprehensive numerical models and laboratory experiments. Firstly, we show the instability in liquid metals, which are relatively poor conductors, differs significantly from the astrophysically-relevant Tayler instability. In liquid metals, the instability develops on the magnetic response time scale, which depends on the conductivity and is much longer than the Alfv\'en time scale, on which the Tayler instability develops in well conducting fluids. Secondly, we show that this instability is an edge effect caused by the curvature of the magnetic field, and its growth rate is determined by the linear current density and independent of the system size. Our results suggest that this instability may affect future liquid metal batteries when their size reaches a few meters.Comment: 14 pages, 5 figures (to appear in J Fluid Mech

Oscillations of weakly viscous conducting liquid drops in a strong magnetic field

We analyse small-amplitude oscillations of a weakly viscous electrically conducting liquid drop in a strong uniform DC magnetic field. An asymptotic solution is obtained showing that the magnetic field does not affect the shape eigenmodes, which remain the spherical harmonics as in the non-magnetic case. Strong magnetic field, however, constrains the liquid flow associated with the oscillations and, thus, reduces the oscillation frequencies by increasing effective inertia of the liquid. In such a field, liquid oscillates in a two-dimensional (2D) way as solid columns aligned with the field. Two types of oscillations are possible: longitudinal and transversal to the field. Such oscillations are weakly damped by a strong magnetic field - the stronger the field, the weaker the damping, except for the axisymmetric transversal and inherently 2D modes. The former are overdamped because of being incompatible with the incompressibility constraint, whereas the latter are not affected at all because of being naturally invariant along the field. Since the magnetic damping for all other modes decreases inversely with the square of the field strength, viscous damping may become important in a sufficiently strong magnetic field. The viscous damping is found analytically by a simple energy dissipation approach which is shown for the longitudinal modes to be equivalent to a much more complicated eigenvalue perturbation technique. This study provides a theoretical basis for the development of new measurement methods of surface tension, viscosity and the electrical conductivity of liquid metals using the oscillating drop technique in a strong superimposed DC magnetic field.Comment: 17 pages, 3 figures, substantially revised (to appear in J. Fluid Mech.

Edge pinch instability of oblate liquid metal drops in a transverse AC magnetic field

This paper considers the stability of liquid metal drops subject to a high-frequency AC magnetic field. An energy variation principle is derived in terms of the surface integral of the scalar magnetic potential. This principle is applied to a thin perfectly conducting liquid disk, which is used to model the drops constrained in a horizontal gap between two parallel insulating plates. Firstly, the stability of a circular disk is analysed with respect to small-amplitude harmonic edge perturbations. Analytical solution shows that the edge deformations with the azimuthal wavenumbers m=2,3,4... start to develop as the magnetic Bond number exceeds the critical threshold Bm_c=3pi(m+1)/2. The most unstable is m=2 mode, which corresponds to an elliptical deformation. Secondly, strongly deformed equilibrium shapes are modelled numerically by minimising the associated energy in combination with the solution of a surface integral equation for the scalar magnetic potential on an unstructured triangular mesh. The edge instability is found to result in the equilibrium shapes of either two- or threefold rotational symmetry depending on the magnetic field strength and the initial perturbation. The shapes of higher rotational symmetries are unstable and fall back to one of these two basic states. The developed method is both efficient and accurate enough for modelling of strongly deformed drop shapes.Comment: 18 pages, 11 figures, corrected final revision, to appear in J. Fluid Mec

Weakly nonlinear stability analysis of MHD channel flow using an efficient numerical approach

We analyze weakly nonlinear stability of a flow of viscous conducting liquid driven by pressure gradient in the channel between two parallel walls subject to a transverse magnetic field. Using a non-standard numerical approach, we compute the linear growth rate correction and the first Landau coefficient, which in a sufficiently strong magnetic field vary with the Hartmann number as $\mu_{1}\sim(0.814-\mathrm{i}19.8)\times10^{-3}\textit{Ha}$ and $\mu_{2}\sim(2.73-\mathrm{i}1.50)\times10^{-5}\textit{Ha}^{-4}$. These coefficients describe a subcritical transverse velocity perturbation with the equilibrium amplitude $|A|^{2}=\Re[\mu_{1}]/\Re[\mu_{2}](\textit{Re}_{c}-\textit{Re})\sim29.8\textit{Ha}^{5}(\textit{Re}_{c}-\textit{Re})$ which exists at Reynolds numbers below the linear stability threshold $\textit{Re}_{c}\sim 4.83\times10^{4}\textit{Ha}.$ We find that the flow remains subcritically unstable regardless of the magnetic field strength. Our method for computing Landau coefficients differs from the standard one by the application of the solvability condition to the discretized rather than continuous problem. This allows us to bypass both the solution of the adjoint problem and the subsequent evaluation of the integrals defining the inner products, which results in a significant simplification of the method.Comment: 16 pages, 10 figures, revised version (to appear in Phys Fluids

Pseudo–magnetorotational instability in a Taylor-Dean flow between electrically connected cylinders

We consider a Taylor-Dean-type flow of an electrically conducting liquid in an annulus between two infinitely long perfectly conducting cylinders subject to a generally helical magnetic field. The cylinders are electrically connected through a remote, perfectly conducting endcap, which allows a radial electric current to pass through the liquid. The radial current interacting with the axial component of magnetic field gives rise to the azimuthal electromagnetic force, which destabilizes the base flow by making its angular momentum decrease radially outwards. This instability, which we refer to as the pseudo--magnetorotational instability (MRI), looks like an MRI although its mechanism is basically centrifugal. In a helical magnetic field, the radial current interacting with the azimuthal component of the field gives rise to an axial electromagnetic force, which drives a longitudinal circulation. First, this circulation advects the Taylor vortices generated by the centrifugal instability, which results in a traveling wave as in the helical MRI (HMRI). However, the direction of travel of this wave is opposite to that of the true HMRI. Second, at sufficiently strong differential rotation, the longitudinal flow becomes hydrodynamically unstable itself. For electrically connected cylinders in a helical magnetic field, hydrodynamic instability is possible at any sufficiently strong differential rotation. In this case, there is no hydrodynamic stability limit defined in the terms of the critical ratio of rotation rates of inner and outer cylinders that would allow one to distinguish a hydrodynamic instability from the HMRI. These effects can critically interfere with experimental as well as numerical determination of MRI.Comment: 10 pages, 5 figures, minor revision, to appear in Phys. Rev.