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. Fluids25, 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

(1)

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

(2)

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 · UU. 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

$\frac{\text{d} x}{\text{d} s}$
dxds⁠.

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

$\frac{\text{d} x}{\text{d} s}$
dxds coming from the tangent vector. The same mistakes appear also in Van Gorder's other articles.14,15 It is also questionable, why ω is real in his Eqs. (3) and (4), when it should be complex. The imaginary part indicates the presence of the dissipation.

Let us also consider, what is Van Gorder's dispersion relation of a Kelvin wave in the absence of the mutual friction. We get

(3)

In the small amplitude limit, Ak ≪ 1 (where xs, or

$\frac{\text{d} x}{\text{d} s} \approx 1$
dxds1⁠), this correctly produces the known dispersion relation ω = γk2. However, this differs from the known results for a large amplitude Kelvin wave, within the local induction approximation, which is given by16 

(4)

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: yy(x, t) and zz(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

(5)

For a Kelvin wave this is

$\frac{\text{d}x}{\text{d}s} =(1+ A^2k^2)^{-1/2}$
dxds=(1+A2k2)1/2⁠. If Ak ≪ 1, this could be approximated to be
$\frac{\text{d}x}{\text{d}s} = 1- (1/2)(y_x^2+ z_x^2) = 1- (1/2)(Ak)^2$
dxds=1(1/2)(yx2+zx2)=1(1/2)(Ak)2
(as Shivamoggi and van Heijst17 did), or even
$\frac{\text{d}x}{\text{d}s}=1$
dxds=1
, but Van Gorder makes no such approximation.

Because Van Gorder chooses to represent the vortex as r = xix + y(x, t)iy + z(x, t)iz, where xx(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

(6)

where

$\mathbf {t} = \frac{\text{d} \mathbf {r}}{\text{d}s}$
t=drds and
$\kappa \mathbf {n} = \frac{\text{d} \mathbf {t}}{\text{d}s}$
κn=dtds
. For simplicity, we consider the zero temperature case (α = 0 and α = 0) here. With Van Gorder's assumption, vx = 0, this leads to three equations, namely,

(7)
(8)

Equation (8) can be combined defining variable Φ = y + iz, the derivative (5) is given by

$\frac{\text{d} x}{\text{d}s} = (1 + |\Phi _x|^2)^{-1/2}$
dxds=(1+|Φx|2)1/2⁠, and we get equation

(9)

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 yxzxxzxyxx = 0 can be given with the help of Φ in the following way:

$\Phi _x^* \Phi _{xx} - \Phi _x \Phi _{xx}^* =0$
Φx*ΦxxΦxΦxx*=0 (in other words the Wronskian of Φx and
$\Phi _x^*$
Φx*
is zero). For a Kelvin wave, Φ = Aexp [i(kx − ωt + x0)], Eq. (7) gives
$\gamma A^2k^3(\frac{\text{d}x}{\text{d}s})^3 = 0$
γA2k3(dxds)3=0
. Here we see that one has to assume that Ak ≪ 1 in order to neglect the term on the left hand side.

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

(10)

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 xx(s, t), and then

$\mathbf {v} = \frac{\text{d} \mathbf {r}}{\text{d} t} = x_t \mathbf {i}_x + (y_t + y_x x_t) \mathbf {i}_y + (z_t + z_x x_t) \mathbf {i}_z$
v=drdt=xtix+(yt+yxxt)iy+(zt+zxxt)iz⁠, will give Φt = yt + izt = 0, if one has
$\mathbf {v} = K \mathbf {t} = K \frac{\text{d} x}{\text{d} s}( \mathbf {i}_x +y_x \mathbf {i}_y + z_x \mathbf {i}_z)$
v=Kt=Kdxds(ix+yxiy+zxiz)
, where K is an arbitrary constant. In other words, terms parallel to t have no effect on Φt.

Using xx(s, t), Eq. (2) can be expressed in the potential form as

(11)

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

(12)
(13)

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

(14)

In the small amplitude limit this is the familiar k = U/γ. From the denominator we get an additional condition

(15)

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

(16)

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

(17)
(18)

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).

1.
R. A.
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
).
2.
D. K.
Cheng
,
M. W.
Cromar
, and
R. J.
Donnelly
, “
Influence of an axial heat current on negative-ion trapping in rotating helium II
,”
Phys. Rev. Lett.
31
,
433
(
1973
).
3.
W. I.
Glaberson
,
W. W.
Johnson
, and
R. M.
Ostermeier
, “
Instability of a vortex array in He II
,”
Phys. Rev. Lett.
33
,
1197
(
1974
).
4.
R. M.
Ostermeier
and
W. I.
Glaberson
, “
Instability of vortex lines in the presence of axial normal fluid flow
,”
J. Low Temp. Phys.
21
,
191
(
1975
).
5.
C. F.
Barenghi
,
M.
Tsubota
,
A.
Mitani
, and
T.
Araki
, “
Transient growth of Kelvin waves on quantized vortices
,”
J. Low Temp. Phys.
134
,
489
(
2004
).
6.
In our notation we have included the factor 1/4π in γ and kept the units, since Van Gorder's notation is only partly dimensionless due to the units in a0. This has no effect on any of the equations below.
7.
R.
Betchov
, “
On the curvature and torsion of an isolated vortex filament
,”
J. Fluid Mech.
22
,
471
(
1965
).
8.
M.
Umeki
, “
A locally induced homoclinic motion of a vortex filament
,”
Theor. Comput. Fluid Dyn.
24
,
383
(
2010
).
9.
H.
Salman
, “
Breathers on quantized superfluid vortices
,”
Phys. Rev. Lett.
111
,
165301
(
2013
).
10.
C. F.
Barenghi
,
R. J.
Donnelly
, and
W. F.
Vinen
, “
Thermal excitation of waves on quantized vortices
,”
Phys. Fluids
28
,
498
(
1985
).
11.
R. J.
Donnelly
,
Quantized Vortices in Helium II
(
Cambridge University Press
,
Cambridge
,
1991
), p.
214
.
12.
B. K.
Shivamoggi
, “
Vortex motion in superfluid 4He: Reformulation in the extrinsic vortex-filament coordinate space
,”
Phys. Rev. B
84
,
012506
(
2011
).
13.
J.
Laurie
,
V. S.
L'vov
,
S.
Nazarenko
, and
O.
Rudenko
, “
Interaction of Kelvin waves and nonlocality of energy transfer in superfluids
,”
Phys. Rev. B
81
,
104526
(
2010
).
14.
R. A.
Van Gorder
, “
Fully nonlinear local induction equation describing the motion of a vortex filament in superfluid 4He
,”
J. Fluid Mech.
707
,
585
(
2012
).
15.
R. A.
Van Gorder
, “
Self-similar vortex dynamics in superfluid 4He under the Cartesian representation of the Hall-Vinen model including superfluid friction
,”
Phys. Fluids
25
,
095105
(
2013
).
16.
E. B.
Sonin
, “
Dynamics of helical vortices and helical-vortex rings
,”
Europhys. Lett.
97
,
46002
(
2012
).
17.
B. K.
Shivamoggi
and
G. J. F.
van Heijst
, “
Motion of a vortex filament in the local induction approximation: Reformulation of the Da Rios–Betchov equations in the extrinsic filament coordinate space
,”
Phys. Lett. A
374
,
1742
(
2010
).
18.
R. A.
Van Gorder
, “
Motion of a vortex filament in the local induction approximation: a perturbative approach
,”
Theor. Comput. Fluid Dyn.
26
,
161
(
2012
).
19.
R. A.
Van Gorder
, “
Integrable stationary solution for the fully nonlinear local induction equation describing the motion of a vortex filament
,”
Theor. Comput. Fluid Dyn.
26
,
591
(
2012
).
20.
M.
Umeki
, “
A real-space representation for a locally induced motion of a vortex filament
,”
Theor. Appl. Mech. Jpn.
61
,
195
(
2013
).
21.
B. V.
Svistunov
, “
Superfluid turbulence in the low-temperature limit
,”
Phys. Rev. B
52
,
3647
(
1995
).
22.
G.
Boffetta
,
A.
Celani
,
D.
Dezzani
,
J.
Laurie
, and
S.
Nazarenko
, “
Modeling Kelvin wave cascades in superfluid helium
,”
J. Low Temp. Phys.
156
,
193
(
2009
).