Classical (viscous) turbulence In a 3D classical turbulent t + ( - - PowerPoint PPT Presentation
Classical (viscous) turbulence In a 3D classical turbulent t + ( - - PowerPoint PPT Presentation
Classical (viscous) turbulence In a 3D classical turbulent t + ( v r ) v = 1 v flow, large scale eddies r p + r 2 v break up into smaller eddies, these into smaller ones and so on...(Richardson Cascade) v = sin( x )
Classical (viscous) turbulence
- In a 3D classical turbulent
flow, large scale eddies break up into smaller eddies, these into smaller ones and so on...(Richardson Cascade)
- If there is a large inertial
range between the forcing and dissipation scale (i.e. high Re) then the flow of energy through scales is characterized by a constant energy flux .
- Dimensional analysis leads
to a power-law scaling for the energy spectrum,
E(k) = C✏2/3k−5/3
∂v ∂t + (v · r)v = 1 ρrp + νr2v
v = sin(x) ⇒ v ∂v ∂x ∼ sin(2x)
Classical Vorticity
ω = r ⇥ u
Quantum Fluids
Γ = ∮$𝐰 ⋅ 𝑒𝐦 ∈ ℝ Γ = ∮$𝐰 ⋅ 𝑒𝐦 = 2𝜌ℏ 𝑛 𝑜
Kuchemann:
“vortices are the sinews and muscles of fluid motions”
If this is true then Quantum Turbulence represents the ‘skeleton’
Yet we still see ‘classical’ behaviour
- Salort et al., 2011
10−1 100 101 10−7 10−6 10−5 10−4
- 5/3
kL0/(2π) φ(k)
Probe cut-off Vortex shedding
−4 −2 2 4 10≠5 10≠4 10≠3
v≠ÈvÍ σp(v) 10−1 5 · 10−1 10−2 10−1 100 r/L0 − 5
4
+ ”v3, /(‘r)
10−1 100 ≠0.2 ≠0.1 r/L0 Èδv3Í Èδv2Í
3/2T = 2.2K, ρs/ρn = 0
<latexit sha1_base64="ERTrTxM/p8YBZYG6rgJVkexURhQ=">ACF3icbVDLSgMxFM34tr6qLt0Ei+BCxpkq6KYguhHcVGhV6JSydzaYCYzJnfEMvQv3Pgrblwo4lZ3/o3p2IWvAyGHc+69uTlhKoVBz/twxsYnJqemZ2ZLc/MLi0vl5ZUzk2SaQ5MnMtEXITMghYImCpRwkWpgcSjhPLw6GvrnN6CNSFQD+ym0Y3apRFdwhlbqlN0A4RaLObmGaJA3alW3erJFg+uMRTQvaRjtotL1bxBp1zxXK8A/Uv8EamQEeqd8nsQJTyLQSGXzJiW76XYzplGwSUMSkFmIGX8il1Cy1LFYjDtvNhnQDesEtFuou1RSAv1e0fOYmP6cWgrY4Y989sbiv95rQy7+1cqDRDUPzroW4mKSZ0GBKNhAaOsm8J41rYXSnvMc042ihLNgT/95f/krOq6+41dPdysHhKI4ZskbWySbxyR45IMekTpqEkzvyQJ7Is3PvPDovzutX6Zgz6lklP+C8fQI2LZ9M</latexit>T = 1.56K, ρs/ρn = 6
<latexit sha1_base64="jyBL1KM+gbkPaO/UeKaqXwkHDk=">ACGXicbVDLSgMxFM34tr6qLt0Ei+BC6ozvTaHoRnCjYG2hU0omvbXBTGZM7ohl6G+48VfcuFDEpa78G9NpF1o9EHI4596bmxPEUh03S9nbHxicmp6ZjY3N7+wuJRfXrkyUaI5VHgkI10LmAEpFRQoIRarIGFgYRqcHPS96t3oI2I1CV2Y2iE7FqJtuAMrdTMuz7CPWZz0kAm0EsvS15x/+Bsi/q3CWtRX3eiptnOLlU6DXzBbfoZqB/iTckBTLEeTP/4bcinoSgkEtmTN1zY2ykTKPgEno5PzEQM37DrqFuqWIhmEabdSjG1Zp0Xak7VFIM/VnR8pCY7phYCtDh0z6vXF/7x6gu2jRipUnCAoPnionUiKEe3HRFtCA0fZtYRxLeyulHeYZhxtmDkbgjf65b/kaqfo7RZ3LvYK5eNhHDNkjayTeKRQ1Imp+ScVAgnD+SJvJBX59F5dt6c90HpmDPsWSW/4Hx+A6WloAs=</latexit>Coherent structures
- In classical turbulence vorticity
is concentrated in vortical ‘worms’ (She & al, Nature, 1990 ; Goto, JFM, 2008)
- Are there vortex bundles in
quantum turbulence ?
- Would allow a mechanism for
vortex stretching, i.e. stretch the bundle.
Dω Dt = (ω · r)v + νr2ω
Mathematical approach
3 distinct scales/numerical approaches
Gross-Pitaevskii Point Vortex/VFM Course-Grained HVBK
Barenghi et al. (2014)
Vortex filament method
Biot-Savart Integral Model reconnections algorithmically ‘cut and paste’
Mutual friction
Normal viscous fluid coupled to inviscid superfluid via mutual friction. Superfluid component extracts energy from normal fluid component via Donelly- Glaberson instability, amplification of Kelvin waves. Kelvin wave grows with amplitude:
Counterflow Turbulence
vext
n (s, t) = (c, 0, 0)
Andronikashvili, 1946
Generation of bundles at finite temperatures
Vortex Locking - Morris, Koplik & Rouson, PRL, 2008 Gaussian normal fluid vortex – Samuels, PRB, 1993
Reconnections: Bundles remain intact
Alamri et al. 2008
Numerical simulations using both GPE and vortex filament method.
Decomposition of a tangle
0.2 0.4 0.6 0.8 1 AWB, Laurie & Barenghi, 2012
Motivation
Roussel, Schneider & Farge, 2005
Numerical results
10 10
1
10
2
10
−7
10
−6
10
−5
10
−4
10
−3
10
−2
10
−1
PSD
f f−5/3
10
2
10
3
10
−5
10
−4
10
−3
k E(k) k−5/3 k−1
a b
Left, frequency spectra (red polarised ; black total), right energy spectrum, upper random component, lower polarised component.
0.2 0.4 0.6 0.8 1Experimental detection
Presence of coherent structures inferred from intermittent pressure drops
Pumped He bath Pressurized HeI / He II (Ø 780 mm cell) Bottom propeller Heat exchanger 702 mm Parietal pressure probes ( Ø 1 mm tap holes 34 mm below equator) Transmission shaft Top propeller Mixing layer
20 40 60 80 100
Time [number of turns] P [arbitr. units and offset]
ρs/ρ = 0 % , Re = 6.6e7 [θ=0.12] ρs/ρ= 19 %, Re = 8.6e7 [θ=0.12] ρs/ρ= 83 %, Re = 8.9e7 [θ=0.11]
Rusaouen et al., 2017
- 10
- 5
10-4 10-3 10-2 10-1
P [standard deviation unit] Probability density
ρs/ρ= 0 %, Re=5.5e7 [θ=0.12] ρs/ρ= 0 %, Re=6.6e7 [θ=0.12] ρs/ρ= 19 % Re=5.9e7 [θ=0.20] ρs/ρ= 19 % Re=8.6e7 [θ=0.12] ρs/ρ= 19 % Re=1.1e8 [θ=0.12] ρs/ρ= 79 % Re=1.3e8 [θ=0.20] ρs/ρ= 79 % Re=1.3e8 [θ=0.11] ρs/ρ= 83 % Re=8.9e7 [θ=0.11] gaussian (standard deviation=1) ρs/ρ= 0 %, Re=5.5e7 [θ=0.20]∂vn ∂t + (vn · r) vn = 1 ρrP + µr2vn + ρs ρ F, r · vn = 0, ∂vs ∂t + (vs · r) vs = 1 ρrP ρn ρ F, r · vs = 0.
F ' αρsh|ωs|i(vs vn),
ρs ρn : r2P ⇠ ρs 2 (ω2
s σ2 s)
<latexit sha1_base64="KlG78PnM1cDk81IXtyIu4X3+Bxg=">ACOnicbZDNThsxFIU9QAsNpQRYsrGIKsGCaGZAnWF2g3LIDWAFIfRHczsfDPxPZUikZ5LjZ9CnYsumEBQmz7AHUmWfDTK1n6dO65ur4nLQS3LgzvgoXFpQ8fl1c+NVY/r31Zb25snltdGsq6VAtLlOwTHDFuo47wS4Lw0Cmgl2k1z+m/YtfzFiu1U83LlhfQq54xik4LyXNM2KGOrGY5DmuUX3DZDQqYCJglTAVYw7mFguMckM0Grmn1TxZJdoyXJI7FW87w25rHEvabCdlgXfg/RHFpoXp2keUsGmpaSKUcFWNuLwsL1KzCOU8EmDVJaVgC9hpz1PCqQzPar+vQJ/uqVAc608U85XKsvJyqQ1o5l6p0S3NC+7U3F/V6pcuO+xVXRemYorNFWSmw03iaIx5w6gTYw9ADfd/xXQIPiHn0274EK3J7+H87gdHbTjs8PWyfd5HCtoG+2gXRShI3SCTlEHdRFN+gPekCPwe/gPngKnmfWhWA+s4VeVfD3HxRFrQI=</latexit>- 10
- 5
10-4 10-3 10-2 10-1
P [standard deviation unit] Probability density
ρs/ρ= 0 %, Re=5.5e7 [θ=0.12] ρs/ρ= 0 %, Re=6.6e7 [θ=0.12] ρs/ρ= 19 % Re=5.9e7 [θ=0.20] ρs/ρ= 19 % Re=8.6e7 [θ=0.12] ρs/ρ= 19 % Re=1.1e8 [θ=0.12] ρs/ρ= 79 % Re=1.3e8 [θ=0.20] ρs/ρ= 79 % Re=1.3e8 [θ=0.11] ρs/ρ= 83 % Re=8.9e7 [θ=0.11] gaussian (standard deviation=1) ρs/ρ= 0 %, Re=5.5e7 [θ=0.20]Hall-Vinen-Bekarevich-Khalatnikov Equations
Course-grained, macroscopic model
A single bundle in isolation
∂vs ∂t + (vs · r) vs = 1 ρrP ρn ρ F,
vs = (vr, vθ, vz) = ✓ 0, NΓ 2⇡r, 0 ◆
P = P0 ⇢sN 2Γ2 8⇡2r2 ,
min
V
P(N) ⇠ N 2.
- 0.5
- 9
- 8
- 7
- 6
- 5
- 4
- 3
- 2
- 1
ˆ F(|k|) = exp |k|2 2k2
f
!
AWB & Laurie arxiv:1910.00276
Turbulent Tangle
1 2 3 4 5 10 -3 10 -2 10 -1 10 0- 5
- 4
- 3
- 2
- 1
ˆ F(|k|) = exp |k|2 2k2
f
!
r kf = 2⇡/lf tering process
x y −0.1 −0.05 0.05 −0.1 −0.08 −0.06 −0.04 −0.02 0.02 0.04 0.06 0.08 0.2 0.4 0.6 0.8 1 1.2vn
<latexit sha1_base64="AEiwS5mjYRGefKXhdlFkIRy+Lek=">AB83icbVDLSgMxFL2pr1pfVZdugkVwVWaqoMuiG5cV7AM6Q8mkmTY0kxmSTKEM/Q03LhRx68+482/MtLPQ1gOBwzn3ck9OkAiujeN8o9LG5tb2Tnm3srd/cHhUPT7p6DhVlLVpLGLVC4hmgkvWNtwI1ksUI1EgWDeY3Od+d8qU5rF8MrOE+REZSR5ySoyVPC8iZhyE2XQ+kINqzak7C+B14hakBgVag+qXN4xpGjFpqCBa910nMX5GlOFUsHnFSzVLCJ2QEetbKknEtJ8tMs/xhVWGOIyVfdLghfp7IyOR1rMosJN5Rr3q5eJ/Xj814a2fcZmkhkm6PBSmApsY5wXgIVeMGjGzhFDFbVZMx0QRamxNFVuCu/rldJp1N2reuPxuta8K+owxmcwyW4cANeIAWtIFCAs/wCm8oRS/oHX0sR0uo2DmFP0CfP4Ukf8=</latexit>- 5
- 4
- 3
- 2
- 1
ωs
<latexit sha1_base64="2enihajQwg4lcYFsMcsB5cRLIqs=">AB73icbVDLSgNBEOyNrxhfUY9eBoPgKexGQY9BLx4jmAckS5idzCZD5rHOzAphyU948aCIV3/Hm3/jJNmDJhY0FXdHdFCWfG+v63V1hb39jcKm6Xdnb39g/Kh0cto1JNaJMornQnwoZyJmnTMstpJ9EUi4jTdjS+nfntJ6oNU/LBThIaCjyULGYEWyd1ekrQIe6bfrniV/050CoJclKBHI1+as3UCQVFrCsTHdwE9smGFtGeF0WuqlhiaYjPGQdh2VWFATZvN7p+jMKQMUK+1KWjRXf09kWBgzEZHrFNiOzLI3E/zuqmNr8OMyS1VJLFojlyCo0ex4NmKbE8okjmGjmbkVkhDUm1kVUciEy+vklatGlxUa/eXlfpNHkcRTuAUziGAK6jDHTSgCQ4PMrvHmP3ov37n0sWgtePnMf+B9/gAe7JAF</latexit>P
<latexit sha1_base64="bzZbWfGTwOVHfOXRUHQsbNaA+g=">AB6HicbVDLSgMxFL1TX7W+qi7dBIvgqsxUQZdFNy5bsA9oB8mkd9rYTGZIMkIZ+gVuXCji1k9y59+YtrPQ1gOBwznkntPkAiujet+O4W19Y3NreJ2aWd3b/+gfHjU1nGqGLZYLGLVDahGwSW2DcCu4lCGgUCO8H4duZ3nlBpHst7M0nQj+hQ8pAzaqzUbDyUK27VnYOsEi8nFch81/9QczSCKVhgmrd89zE+BlVhjOB01I/1ZhQNqZD7FkqaYTaz+aLTsmZVQYkjJV90pC5+nsio5HWkyiwyYiakV72ZuJ/Xi814bWfcZmkBiVbfBSmgpiYzK4mA6QGTGxhDLF7a6EjaizNhuSrYEb/nkVdKuVb2Laq15Wanf5HU4QRO4Rw8uI63EDWsA4Rle4c15dF6cd+djES04+cwx/IHz+QOq54zY</latexit>Random ‘Vinen’ Tangle
Quasiclassical UltraquantumWalmsley et al. 2013
- 5
- 4
- 3
- 2
- 1
- 5
- 4
- 3
- 2
- 1
1 2 3 4 5 10 -3 10 -2 10 -1 10 0
F ' αρsh|ωs|i(vs vn),
- 5
- 4
- 3
- 2
- 1
Summary
- Coherent vortical structures are present in the quasi-
classical regime of Quantum Turbulence.
- Important (essential?) for K41 like statistical
properties of QT.
- Good agreement between macroscopic HVBK model
and mesoscale vortex approach.
- Interesting high pressure signal found in the Vinen
regime.