APhO 2017, Теория — Задача 1. Vortices in superfluid

Theory APhO YAKUTSK 2017 RUSSIA Q1 English (Official)

Vortices in Superfluid

MODD-Problems

May 5, 2017

Условие

Introduction

Superfluidity is a property of flowing without friction. Everyday experience tells us that motion of an ordinary fluid (say, water at room temperature) is always accompanied by viscous dissipation of energy, so that the flow gradually becomes slower unless it is maintained by external forces. In contrast, superfluid exhibits no loss of kinetic energy: once excited the motion of superfluid can continue indefinitely. Superfluidity was originally discovered experimentally in liquid helium.

We study properties of superfluid helium at zero temperature. It will be treated as an incompressible fluid with density ρ\rho. Flow continuity (the fact that the mass flowing into and the mass flowing out of a given infinitesimal volume are equal) implies that the flux of helium velocity v⃗\vec{v} through a closed surface is always zero. Superfluid velocity in this aspect is analogous to the magnetic field intensity. By analogy with the magnetic field lines, “streamlines” are tangential to the fluid velocity at each point and their density is proportional to its magnitude.

True superflow has an important property of being irrotational: circulation of superfluid velocity v⃗\vec{v} along any closed path within helium is zero

∫Lv⃗⋅dl⃗=0(1)\int\limits_L \vec{v}\cdot d\vec{l} = 0 \qquad (1)

This statement must be amended if superfluidity is absent along a thin “vortex filament”. The thickness of the filament itself is of approximately atomic dimensions aa, but the vortex induces long range velocity field in surrounding superfluid: velocity circulation around such filament is the circulation quantum1^1

∣∫Lv⃗⋅dl⃗∣=2πκ,(2)\left|\int\limits_L \vec{v}\cdot d\vec{l}\right| = 2\pi\kappa, \qquad (2)

and zero if the path can be contracted to a single point without crossing the vortex, see Fig. 1. This supports the analogy between superflow and the magnetic field created by wires carrying current: superposition of two valid velocity distributions is a valid velocity distribution and the velocity at any point is equal (up to a dimensional factor) to the magnetic field produced by the unit currents running through a system of wires representing vortex filaments.

1^1Circulation quantization is a macroscopic quantum effect and corresponds to the angular momentum quantization in Bohr model. The circulation quantum can be expressed as κ=ℏ/mHe\kappa = \hbar/m_{\mathrm{He}}, where mHem_{\mathrm{He}} is the mass of helium atom.

A light-blue rectangle of superfluid crossed from top to bottom by a curved red vortex filament of thickness a (marked by a dimension at the top). Six closed paths with arrows are drawn: L1 and L3 loop near the filament, L2 is a small loop at the upper right, L4 is a large loop, L5 is a loop around part of the filament, L6 loops across the filament near the bottom.
Fig. 1: Vortex filament (red) in superfluid (light blue). Velocity circulations along paths L1L_1, L2L_2, L5L_5, and L6L_6 are all zero, but those for L3L_3 and L4L_4 are equal to ±2πκ\pm 2\pi\kappa. Note that circulations along L3L_3 and L4L_4 have opposite signs.

Part A. Steady filament (0.75 points) (0,75 т.)

Consider a cylindrical beaker (radius R0≫aR_0 \gg a) of superfluid helium and a straight vertical vortex filament in its center Fig. 2.

A cylindrical beaker filled with light-blue helium; a red vertical vortex line along its axis with a circular arrow around it; a radius arrow labelled R0 at the bottom; a downward arrow labelled g on the right; coordinate axes x, y, z at the lower left.
Fig. 2: Straight vortex along the axis of a beaker.

A.1 Plot the streamlines. Find out the velocity vv at a point r⃗\vec{r}. [0,25 т.]

A.2 Work out the free surface shape (height as a function of coordinate z(r⃗)z(\vec{r})) around the vortex. Free fall acceleration is gg. Surface tension can be neglected. [0,5 т.]

Part B. Vortex motion (1.4 points) (1,4 т.)

Free vortices move about in space with the flow2^2. In other words each element of the filament moves with the velocity v⃗\vec{v} of the fluid at the position of that element.

2^2This is a consequence of momentum conservation, see next section.

As an example, consider a pair of counter-rotating straight vortices placed initially at distance r0r_0 from each other, see Fig. 3. Each vortex produces velocity v0=κ/r0v_0 = \kappa/r_0 at the axis of another. As a result, the vortex pair moves rectilinearly with constant speed v0=κ/r0v_0 = \kappa/r_0 so that the distance between them remains unchanged.

Two red dots (vortex cross-sections) a distance r0 apart, marked by a double arrow between dashed lines. Each has an upward velocity arrow labelled v0; the left one has a counter-clockwise curved arrow, the right one a clockwise curved arrow.
Fig. 3: Parallel vortex filaments with opposite circulations.
Two red dots a distance r0 apart, marked by a double arrow between dashed lines; each has a clockwise curved arrow beside it.
Fig. 4: Parallel vortex filaments with equal circulations.
Top view of a circular beaker of diameter 2R0 filled with a regular array of small vortex symbols; arrows mark vortex A near the wall, vortex B midway and vortex C at the center. A green square region is magnified at the right, showing the triangular lattice with spacing u marked along three directions by dash-dotted lines.
Fig. 5: Triangular lattice of vortices in a beaker. The view from above.

B.1 Consider two identical straight vortices initially placed at distance r0r_0 from each other as shown in Fig. 4. Find initial velocities of the vortices and draw their trajectories. [0,25 т.]

A beaker of helium (see Part A) is filled with triangular lattice (u≪R0u \ll R_0) of identical vertical vortices, see Fig. 5.

B.2 Draw the trajectories of vortices A, B, and C (located in the center). [0,15 т.]

B.3 Find velocity v(r⃗)v(\vec{r}) of a vortex positioned at r⃗\vec{r}. [0,4 т.]

B.4 Find the distance AB(t)(t) between the vortices A and B at time tt. Treat AB(0)(0) as given. [0,35 т.]

B.5 Work out the "smoothed out" (omitting the lattice structure) free helium surface shape z(r⃗)z(\vec{r}). [0,25 т.]

Part C. Momentum and energy (1.75 points) (1,75 т.)

The long range velocity field is the major contribution to the energy of a system of vortices, it is insensitive to exact structure of the filament. The filament itself can not be properly described by the macroscopic theory and apparent singularities (infinities) are insignificant. Real physical quantities, such as the energy, of the region inside a thin tube of radius aa around the filament should be neglected. Outside of this tube the density of superflow kinetic energy ρv2/2\rho v^2/2 (where ρ=const\rho = \mathrm{const}) is analogous to the energy density of the magnetic field B2/(2μ0)B^2/(2\mu_0) — they are both quadratic in respective variables. This analogy together with the correspondence between magnetic field and superfluid velocity generated by vortices (currents) facilitates calculation of the flow energy for a given system. For instance, given the inductance of a circular wire loop L≈μ0Rlog⁡(R/a)L \approx \mu_0 R \log(R/a), where RR is the loop radius and aa is wire radius, we get the superfluid vortex loop energy3^3

U≈2Rρπ2κ2log⁡(R/a)(3)U \approx 2R\rho\pi^2\kappa^2 \log(R/a) \qquad (3)

Total fluid momentum is also determined by the long range velocity distribution. It is obtained by integration of the momentum density ρv⃗\rho\vec{v}. Again, consider a flow generated by a circular vortex loop placed in xyxy plane. It is obvious from the symmetry considerations, that total momentum has only zz component:

P=∫ρvzdV=ρ∫∫(∫vzdz)⏟q(x,y)dxdy(4)P = \int \rho v_z dV = \rho \int\int \underbrace{\left(\int v_z dz\right)}_{q(x,y)} dxdy \qquad (4)

The innermost integration is in fact an integration along appropriate paths parallel to zz-axis, see Fig. 6. From the circulation identity (2) it follows that

q(x,y)=∫L(x,y)v⃗⋅dl⃗(5)q(x,y) = \int\limits_{L(x,y)} \vec{v}\cdot d\vec{l} \qquad (5)

is piecewise constant. Particularly, it is zero for paths passing outside the ring and 2πκ2\pi\kappa for paths inside it. Total momentum is therefore

P=ρ⋅πR2⋅2πκ=2π2ρR2κ(6)P = \rho\cdot\pi R^2\cdot 2\pi\kappa = 2\pi^2\rho R^2\kappa \qquad (6)

3^3This expression is also valid only if log⁡R/a≫1\log R/a \gg 1.

Left: a red circular vortex loop in perspective with black streamlines looping around it and passing upward through its interior. Right: cross-section in the x–z plane with the two red cross-section points of the loop, black closed and open streamlines around them, and seven vertical green lines with upward arrows marking the integration paths; axes z and x.
Fig. 6: Velocity field of a circular vortex loop and integration paths (green) for q(x,y)q(x,y) calculation.
A red elongated rounded rectangle of length d (dimension line above) and height b (dimension at the right), with small curved arrows at the top and bottom sides showing the circulation.
Fig. 7: A nearly rectangular vortex loop, b≪db \ll d.
Axes x, y, z in perspective; a red straight vortex filament along z with a circulation arrow, and a faded copy shifted by a distance b in the x direction (arrow labelled b).
Fig. 8: Momentum changes whenever the vortex shifts with respect to the fluid.

C.1 Consider a nearly rectangular vortex loop b×db \times d, b≪db \ll d, Fig. 7. Indicate the direction of its momentum P⃗\vec{P}. Find out the momentum magnitude. [0,3 т.]

C.2 Calculate its energy UU. [0,7 т.]

C.3 Suppose we shift a long straight vortex filament by a distance bb in xx direction, see Fig. 8. How much does the fluid momentum change? Indicate the momentum change direction. The filament length (constrained by the vessel walls) is dd. [0,75 т.]

Part D. Trapped charges (2.85 points) (2,85 т.)

Electrons, if injected in helium, get trapped in the vortex filaments. Here and below polarizability of helium can be neglected (ϵ=1\epsilon = 1).

Axes x and y; a red dot (vortex cross-section) with a curved circulation arrow beneath it, and to the right a horizontal arrow labelled E pointing in the +x direction.
Fig. 9: Straight vortex in a uniform electric field.
Left: a red vortex ring seen obliquely with its center C and a downward arrow labelled E. Right: axes x and z with the two red cross-section points of the ring, each with a circulation arrow, the center C between them, and a downward arrow labelled E.
Figure 10: (left) Vortex ring in a uniform electric field. (right) Cross section of the ring.

D.1 Consider a straight vortex charged with uniform linear density λ<0\lambda < 0 in a uniform electric field E⃗\vec{E}. Draw the vortex trajectory. Find its velocity as a function of time. [0,5 т.]

A circular vortex loop of radius R0R_0 initially charged with uniform linear density λ<0\lambda < 0 is placed in a uniform electric field E⃗\vec{E} perpendicular to its plane, opposite to its momentum P⃗0\vec{P}_0.

D.2 Draw the trajectory of the loop center CC. Find the radius of the loop as a function of time. [0,6 т.]

D.3 Find its velocity v(t)v(t) as a function of time. [1,5 т.]

D.4 The field is switched off at a time t∗t^* when the velocity reaches the value v∗=v(t∗)v^* = v(t^*). Find the loop velocity v(t)v(t) at a later time t>t∗t > t^*. [0,25 т.]

Part E. Influence of the boundaries (3.25 points) (3,25 т.)

Solid walls alter the velocity field created by a vortex filament, because the fluid cannot flow through them. Mathematically this means that the wall-normal velocity component vanishes at the wall surface.

A red dot (vortex cross-section) with a clockwise curved arrow above it, at height h0 above a flat grey wall, the distance marked by a vertical double arrow labelled h0.
Fig. 11: Straight vortex filament near a flat wall.
A right-angle corner formed by a grey vertical wall on the left and a grey horizontal wall at the bottom; a red dot (vortex) with a curved circulation arrow, at distance h0 from the vertical wall and h0 from the bottom wall, both marked by double arrows labelled h0.
Fig. 12: Straight vortex filament in a corner.

E.1 Draw the trajectory of a straight vortex, initially placed at a distance h0h_0 from a flat wall. Find its velocity as a function of time. [0,5 т.]

Consider a straight vortex placed in a corner at a distance h0h_0 from both walls.

E.2 What is the initial velocity v0v_0 of the vortex? [0,75 т.]

E.3 Draw the trajectory of the vortex. [0,5 т.]

E.4 What is the velocity of the vortex v∞v_\infty after very long time? [1,5 т.]

Отговори

Решение


Оригинал в Архива: 1.pdf · официални решения: 1_sol.pdf