We comment on the paper by Van Gorder [“Motion of a helical vortex filament in superfluid 4He under the extrinsic form of the local induction approximation,” Phys. Fluids 25, 085101 (2013)]. We point out that the flow of the normal fluid component parallel to the vortex will often lead into the Donnelly–Glaberson instability, which will cause the amplification of the Kelvin wave. We explain why the comparison to local nonlinear equation is unreasonable, and remark that neglecting the motion in the x-direction is not reasonable for a Kelvin wave with an arbitrary wavelength and amplitude. The correct equations in the general case are also derived.
Kelvin waves, which are helical perturbations of a vortex filament, are much studied in the context of superfluid turbulence due to the role they are believed to have in the dissipation of energy in the zero temperature limit. In his article,1 Van Gorder studies a helical vortex at finite temperature and the normal fluid velocity being aligned along the axis of the vortex. Van Gorder states that “for non-zero values of the superfluid friction parameters, there are no helical solutions studied for any of the models in the literature, full or approximate.” We would argue that this is not the case.
When the normal fluid flow is parallel to the vortex, solutions with a constant amplitude are only exceptions. In the presence of mutual friction helical solutions are usually decaying or unstable. It is well known that normal fluid velocity parallel to the vortex may lead to the Donnelly–Glaberson instability.2–4 This instability causes the amplification of Kelvin waves. A simple derivation to show the Donnelly–Glaberson instability (in the small amplitude limit and when α′ = 0) can be found, for example, in an article by Barenghi et al.5 In his paper, Van Gorder fails to mention the Donnelly–Glaberson instability completely.
Let us consider a straight vortex with a small amplitude Kelvin wave. If one is in the zero temperature limit, where there is no dissipation due to the mutual friction (the mutual friction parameters α and α′ are zero), then all helical solutions stay as they are. If α ≠ 0 and we have a flow of the normal fluid component parallel to the vortex, then there is only one critical wave number for which a helical vortex keeps its shape: k = U/γ, where k is the wave number, U is the velocity of the normal component, and γ = (Γ/4π)ln (c/⟨κ⟩a0) (Γ ≈ 9.97 × 10−4 cm2/s is the quantum of circulation, c is a scaling factor of the order of unity, ⟨κ⟩ is the average curvature, and a0 ≈ 1.3 × 10−8 cm is the effective vortex core radius).6 The condition for the Donnelly–Glaberson instability to occur is k < U/γ. In this wave number regime the Kelvin wave will grow exponentially. If k > U/γ, the wave will be damped leading eventually to a straight vortex.
Mathematically, a Kelvin wave with the critical wave number will remain as it is. Physically this kind of solution is not interesting, because small perturbations would cause it either to damp or grow. Also, for infinitely long vortex, with non-zero U, the instability always occurs at small enough k. The question of the stability of helical vortices even at the zero temperature is interesting. Already Betchov has shown that they are unstable.7 Umeki8 and Salman9 have shown additionally that under the local induction approximation (LIA) there is a breather solution of vortex motion, where the shape of the vortex filament is helical as t → ±∞. Salman9 also showed that the breather excitations of the vortex can lead to self-reconnections and generation of small vortex rings.
The effects of the mutual friction on the dispersion relation of a Kelvin wave are also known.10,11 Mutual friction will shift the resonance frequency by a factor (1 − α′) and cause the resonance frequency to broaden (induce dissipation) by a factor proportional to α. The presence of a normal fluid flow and a finite amplitude will further slightly modify the resonance frequency and damping, as will be shown below.
After this brief review of some known results, let us now have a look at Van Gorder's article. First of all, an important part of Van Gorder's article is to compare his results with approximations to the fully nonlinear form of the LIA. The approximations that were considered were a cubic model of Shivamoggi12 and local nonlinear equation (LNE) of Laurie et al.13 Comparison to Shivamoggi's approximation is reasonable, since Van Gorder's and Shivamoggi's equations are obtained in a similar way. On the other hand, comparison to LNE is unreasonable. Van Gorder's equation includes mutual friction and LNE does not. This is very clear from the original paper by Laurie et al.13 The local nonlinear equation was given in the context of studying the Kelvin wave cascade, which is believed to be essential for the energy dissipation in quantum turbulence at zero temperature. In the zero temperature limit there is no mutual friction.
When using the local induction approximation, the velocity of the vortex is
where κ is the local curvature (not average curvature, as one might conclude from Van Gorder's paper), t and n are the unit tangent and normal vectors to the vortex filament, and α and α′ are the dimensionless friction parameters. Van Gorder assumes the normal fluid velocity to be U = Uix and says that then the equation of motion reduces to
However, the second last term should actually be −α′(t · U)t. Van Gorder refers to a paper by Shivamoggi12 when he gives this equation. Shivamoggi12 assumed that the vortex is essentially aligned along the x-axis and that is why he apparently considered that t · U ≈ U. This assumption is not reasonable for a vortex with a Kelvin wave of arbitrary wavelength and amplitude. The dot product should give an extra factor of
One should be able to discard the tangential term, because movement in the direction of the tangent does not change the form of the vortex, hence it is not important. However, the discard of this term would change the dispersion relation. This is not the only peculiarity that catches reader's eye. In Van Gorder's Eq. (2) the term −αUΦx, corresponding to the component in the direction of t × ix, misses the factor
Let us also consider, what is Van Gorder's dispersion relation of a Kelvin wave in the absence of the mutual friction. We get
In the small amplitude limit, Ak ≪ 1 (where x ≈ s, or
In order to explain all these deviations from known results, let us see how Van Gorder has derived his equation. The very first assumption one must make is to assume the vortex to be single-valued in x-direction. Then one can give the y- and z-coordinates as functions of x and time: y ≡ y(x, t) and z ≡ z(x, t). For a Kelvin wave we have r = xix + Acos (kx − ωt + x0)iy + Asin (kx − ωt + x0)iz. Then from the fact ds2 = dx2 + dy2 + dz2 it also follows that
For a Kelvin wave this is
Because Van Gorder chooses to represent the vortex as r = xix + y(x, t)iy + z(x, t)iz, where x ≡ x(s) is not a function of time, he has to assume that the velocity has no component in x-direction. This assumption is unreasonable. The motion of a helical vortex is purely rotational only when Ak → 0. If Ak → ∞, then a segment of the vortex is almost like a segment of a vortex ring, and the motion of the vortex consists of translation to x-direction. Usually the motion is a combination of both rotation and translation.
So, if Van Gorder assumes that v = ytiy + ztiz, he should assume that Ak ≪ 1. This should be also clear from the equations he gets, when setting
where
Equation (8) can be combined defining variable Φ = y + iz, the derivative (5) is given by
which gives the incorrect dispersion relation (3). This is the equation which Van Gorder obtained and used in his earlier works.14,18,19 In most of his papers Van Gorder does not mention the constraint resulting from xt = 0 (Eq. (7)), nor does Shivamoggi and van Heijst12,17 or Umeki20 ever consider this. Only in his first article18 concerning LIA in extrinsic coordinate space Van Gorder does mention this.
The condition yxzxx − zxyxx = 0 can be given with the help of Φ in the following way:
In the paper under consideration, it seems that in making the change of variable for the velocity term in the direction of the normal vector, he exploits this requirement. However, this is not totally clear, since if the mutual friction terms are included, then the expression for xt is more complicated. Typos in the equation giving the normal vector κn may cause even more confusion (same typos appear also in Refs. 14,15). The normal vector should be
If Van Gorder would simply neglect all terms with A2k2, he would recover all the known results obtained in the small amplitude limit. The amplitude of a helical wave would not be anymore related to the wave number. One would obtain the same critical wave number k = U/γ as with Shivamoggi's model17 and which was shown also in the article by Barenghi et al.5
Still, one would be left wondering why the term coming from the tangential velocity shows in the dispersion relation. One should be able to discard the tangential velocity components. This problem is solved if one takes into account the movement in the x-direction. Having x ≡ x(s, t), and then
Using x ≡ x(s, t), Eq. (2) can be expressed in the potential form as
The resulting dispersion relation in zero temperature is the dispersion relation for a large amplitude Kelvin wave (4). The dispersion relation with mutual friction is
where ω0 = γk2(1 + A2k2)−1/2 is the zero temperature dispersion relation.
Requiring Im(ω) = 0 will give us the critical wave number for the Donnelly–Glaberson instability, i.e., the wave number for which the amplitude stays constant in time. We will notice that the critical k depends actually on the amplitude. We get either k = 0 (a straight vortex) or
In the small amplitude limit this is the familiar k = U/γ. From the denominator we get an additional condition
This is a reasonable condition. If one had γ < UA, then a vortex ring with radius A moving in the direction of the flow of the normal component would grow.
In the zero temperature limit Eq. (11) reduces to
An easier way to derive Eq. (16) would have been using the Hamiltonian formulation. The Hamiltonian corresponding to LIA21 and the equation of the motion are given by
Equation (16) is the same as, for example, Eq. (4) in article by Boffetta et al.22 The advantage of deriving the equation via formulating LIA in the Cartesian coordinates is that the inclusion of the mutual friction is straightforward. In their paper Boffetta et al.22 have also an equation of motion with a mutual friction term (Eq. (27)), where this term is proportional to Φ, when it evidently should be proportional to Φxx.
The authors would like to thank H. Salman for useful discussions on this comment. This work is supported by the Academy of Finland (Grant Nos. 218211 and 263798) and in part by the European Union 7th Framework Programme (FP7/2007-2013, Grant No. 228464 Microkelvin) and by the Academy of Finland through its LTQ CoE grant (Project No. 250280).