arxiv: v2 [cond-mat.stat-mech] 27 Sep 2011

Transcription

1 The sensiiviy of he vorex filamen mehod o differen reconnecion models A. W. Baggaley School of Mahemaics and Saisics, Universiy of Newcasle, Newcasle upon Tyne, NE1 7RU, UK arxiv: v2 [cond-ma.sa-mech] 27 Sep 2011 We presen a deailed analysis on he effec of using differen algorihms o model he reconnecion of vorices in quanum urbulence, using he hin-filamen approach. We examine differences beween four main algorihms for he case of urbulence driven by a counerflow. In calculaing he velociy field we use boh he local inducion approximaion (LIA) and he full Bio-Savar inegral. We show ha resuls of Bio-Savar simulaions are no sensiive o he paricular reconnecion mehod used, bu LIA resuls are. I. INTRODUCTION Turbulence in he quanum, low emperaure phase of liquid helium ( 4 He), also known as quanum urbulence [1 3], consiss of reconnecing quanized vorex filamens, arranged in a random, disordered angle. Due o quanum mechanical consrains, each vorex carries he same fixed circulaion h/m, where h is Plack s consan and m is he mass of one aom. This quaniy is called he quanum of circulaion κ. A number of experimenal mehods have been used o creae his form of urbulence, for example agiaing superfluid liquid helium wih propellers [4, 5], forks[6], or grids [7]; hese echniques are also used o creae urbulence in ordinary fluids. A finie emperaures superfluid helium is a wo fluid sysem: a viscous normal fluid componen coexising wih an inviscid superfluid componen. The superfluid vorices inerac wih he hermal exciaions which make up he normal fluid, hus inroducing a muual fricion force beween he wo fluid componens. This means ha urbulence in he quanum fluid can be driven by he flow of he normal fluid or vice-versa. One paricular example is superfluid urbulence driven by a hea flow [8 10]. Noe ha his form of urbulence has no classical analogy. Recenly new flow visualizaion echniques, such as racer paricles[11, 12], Andreev scaering[13] and laser-induced fluorescence[14], have added subsanially o our knowledge of he naure of quanum urbulence. Even wih hese advances in flow visualizaion numerical simulaions will coninue o play a crucial role in furhering our knowledge. Indeed, building on he pioneering work of Schwarz [15], numerical simulaions have always been imporan in he field [16 26], and he recen experimenal progress has highlighed heir imporance in inerpreing experimenal daa [27, 28]. This aricle is concerned wih verifying ha he vorex reconnecion procedure used in he popular vorex filamen mehod (VFM) is robus. Before we discuss he modelling of reconnecions in he VFM, we shall give a brief ouline of he VFM mehod. In superfluid helium, he vorex core radius (a cm) is many orders of magniude smaller han he average separaionbeweenvorexlines (ypicallyfrom10 2 o10 4 cm) oranyoherrelevanlengh scalein he flow. Saring from his key observaion, Schwarz [15] modelled vorex lines as spaces curves s = s(ξ, ) of infiniesimal hickness, where is he ime and ξ is arc lengh, using he classical heory of vorex filamens [29]. In he VFM hese space curves are numerically discreized by a large, variable number of poins s i (i = 1, N), which hereafer we refer o as vorex poins. In all implemenaions of he VFM o superfluids he number of vorex poins varies in ime, as i is desirable o mainain a relaively consan resoluion along he filamens by adding or removing poins as he lengh of he filamens changes. In a recen paper [25] we have discussed he various numerical deails a lengh, here i is suffice o say ha he resoluion along he filamens lies beween an upper, δ, and lower, δ/2, bound. The moion of a vorex filamen is deermined by he normal fluid, hrough muual fricion, and by he induced velociy of he oher vorices. The governing equaion of moion of he superfluid vorex lines, a poin s is given by he Schwarz equaion [15] ds d = v s +αs (v n v s ) α s [s (v n v s )], (1) where α, α are emperaure dependen fricion coefficiens [30, 31], v n is he normal fluid s velociy, and he velociy v s is he velociy field he quanised vorices alone induce. The prime denoes derivaive wih respec o arclengh, e.g. s = ds/dξ. Noe ha, sricly speaking, one should also model he back reacion of he superfluid on he normal Elecronic address: a.w.baggaley@ncl.ac.uk

2 FIG. 1: Schemaic reconnecion procedure for Type I and II reconnecions described in secion II. If wo vorex poins are closer han a criical disance = δ/2 hen a flag swapping operaion changes he opology of he filamens and performs he reconnecion. From (1 2 3) and (4 5 6) o ( ) and (1 4). componen; v n should be he soluion of he Navier-Sokes equaion modified by muual fricion. However, for he sake of simpliciy, here we shall simply prescribe he normal flow. In our sysem he quanised vorices define he voriciy field; we recover he velociy field v s by numerically solving he Bio-Savar (BS) inegral, v s (s) = κ (s r) dr, (2) 4π L s r 3 where he line inegral exends over he enire vorex configuraion L. We de-singularize he inegral in a sandard way [15]. Schwarz [15] proposed he use of he Local Inducion Approximaion (LIA) [32, 33] as an alernaive o he Bio- Savar mehod, which is compuaionally very expensive. The LIA ignores he non-local conribuion o he moion of a secion of he vorex filamen. Insead a vorex line a he poin s i moves along he binormal o s i, wih he following equaion of moion ( cr ds i d = κ 4π ln a 0 ) s i s i, (3) where c is a consan of order uniy and R is he local radius of curvaure s i 1 of he vorex filamen. We shall es various reconnecion algorihms using boh he full BS law, and he LIA in a series of numerical experimens. Again, for furher deails abou spaial derivaives and mehods o ime-sep he vorices, we refer he reader o [25]. We know from experimens [34] and from more microscopic models[35 38] ha superfluid vorex lines can reconnec wih each oher when hey come sufficienly close, as envisaged by Feynman [39]. Superfluid vorex reconnecions do no violae Kelvin s heorem as near he axis of he vorex core, where densiy and pressure vanish and velociy diverges, he governing Gross-Piaevski equaion (GPE) differs from he classical Euler equaion. Whils vorex reconnecions are naural soluions of he GPE, wihin he VFM reconnecions mus be modelled by supplemening Eq. (2) wih an algorihmical reconnecion procedure. This was originally proposed by Schwarz [15], and since hen a number of alernaive algorihms have been proposed. Unil now no deailed es of he effecs of varying his procedure has been performed. Whils a number of sudies have shown good agreemen beween resuls using he VFM and experimenal resuls [26, 40, 41], a deailed sudy of he reconnecion algorihm is imely; ha is he purpose of his sudy. II. RECONNECTION ALGORITHMS In his pioneering paper [15], Schwarz s suggesed ha vorex lines reconnec whenever he disance beween a pair of vorices is less han = 2R/[cln(R/a 0 )], where R is he radius of curvaure a he reconnecion poin, and c is a consan of order uniy. However, his approach can lead o non-physical reconnecions. Consider, for example, wo

3 FIG. 2: Schemaic reconnecion procedure for Type III reconnecion algorihm described in secion II. The segmens (1 2 3) and (4 5 6) evolve ino (6 3) and (1 4). Poins 2 and 5 are eliminaed. FIG. 3: Schemaic represenaion of he moion of vorex filamens under he Type IV algorihm. If he segmens will collide (as shown) hen he reconnecion is performed as in Fig. 1. almos sraigh vorices; under Schwarz s crierion hese vorices mus reconnec even if hey are very far apar, as a large radius of curvaure R, resuls in a large value for. In his work we herefore avoid he use of his reconnecion crierion, insead focusing on mehods which have been used in more recen sudies. Mos recen sudies have relaed he reconnecion disance o he space resoluion along he filamens [42, 43]. We consider hree differing reconnecion algorihms which ake he criical reconnecion disance o be relaed o he maximum spaial resoluion, = δ/2. We define he firs reconnecion algorihm as Type I. In a Type I reconnecion vorices are simply reconneced if heir separaion if less han. Moivaed by he fac ha a reconnecion is dissipaive even, leading o phonon emission [44], we consider an improved algorihm denoed Type II. As vorex line lengh is aproxyfor he kineic energy, foratypeii reconnecionno only mus he disancebeween he reconnecing filamen be less ha δ/2, bu also he line lengh mus be reduced by he change of opology. A schemaic for he Type I and II algorihms can be seen in Fig. 1. The hird reconnecion model, Type III, can be considered an ulra dissipaive algorihm. Type III differs from Type I because he poins which riggered he reconnecion are eliminaed leading o a greaer loss of line lengh, and hence energy. A schemaic of his algorihm is displayed in Fig. 2. In a recen work Kondaurova and Nemirovskii [21] inroduced a new model (Type IV) of reconnecions wihin he VFM, which ess wheher line segmens will mee during he ime-sep. Candidae pairs of vorex poins (i & j), poins which are close enough o reconnec, are idenified. The Type IV algorihm hen assumes ha he line

4 segmens beween each of he pair will move wih a consan velociy (v(s i ) & v(s j )) during he ime-sep. Under his assumpion a se of simulaneous equaions can be consruced, x i +v x (s i ) +(x i+1 x i )ψ = x j +v x (s j ) +(x j+1 x j )φ 4 y i +v y (s i ) +(y i+1 y i )ψ = y j +v y (s j ) +(y j+1 y j )φ z i +v z (s i ) +(z i+1 z i )ψ = z j +v z (s j ) +(z j+1 z j )φ 0 φ 1,0 ψ 1 such ha if a soluion for φ and ψ can be found, hen he filamens will collide during he ime-sep. Here s i = (x i,y i,z i ), s i+1 = (x i+1,y i+1,z i+1 ) s j = (x j,y j,z j ) and s j+1 = (x j+1,y j+1,z j+1 ) are he coordinaes of he pair of poins and heir neighbours along he filamen, and is he numerical ime-sep. If he line segmens will mee, hen he filamens are reconneced in he same manner as in Fig. 1. We solve for ψ and φ algebraically using he equaions for he x and y componens. The values of ψ and φ obained are hen esed in he equaion for he z componen by verifying ha, z i +v z (s i ) +(z i+1 z i )ψ z j +v z (s j ) +(z j+1 z j )φ < ǫ, (4) where ǫ is he imposed he olerance of he mehod. We define candidae poins (i & j) as vorex poins wihin disance = δ. A schemaic diagram of he Type IV algorihm is shown in Fig. 3. Before performing any reconnecion we es he disance from a poin i o all oher poins, which are no he neares neighbours along he filamen. We hen begin he es for a reconnecion based on he closes reconnecion poin, assuming he disance is smaller han. This means self-reconnecions (which can arise if a vorex filamen has wised by a large amoun) are reaed in he same manner as reconnecions beween differen filamens. Finally as reconnecions mus preserve he orienaion of he voriciy, we check ha he wo reconnecing filamens are no parallel. To do his we form local (uni) angen vecors, for example if he reconnecing poins are i and j we can readily calculae ŝ i and ŝ j. We hen es ha ŝ i ŝ j < 0.965, which ensures ha he minimum angle beween reconnecing filamens is approximaely 15 degrees. We have esed ha alering his hreshold o smaller angles does no lead o any discernible difference in he resuls, hence for breviy all resuls here are subjec o his crierion. III. NUMERICAL SIMULATIONS We now he numerical simulaions used o es he various reconnecion algorihms described in he previous secion. As a benchmark we choose counerflow urbulence, he relaive moion of he normal fluid and superfluid componens susained by an applied hea flow in he direcion of v n. The superfluid flows in he opposie direcion so ha ρ n v n + ρ s v s = 0, where ρ n and ρ s are respecively he normal fluid and superfluid densiies. We choose his form of urbulence as here has been a wealh of experimenal sudies [8, 10, 45], as well as some recen deailed numerical simulaions [40, 46]. Counerflow urbulence was also used in a recen numerical sudy by Kondaurova and Nemirovskii [47], where hey repored ha he use of a reconnecion mehod, similar o he Type IV algorihm, gave a seady sae soluion for LIA simulaions. We shall also es hese claims in his sudy. Our calculaions are performed in a periodic cube wih sides of lengh D = 0.1 cm. Superfluid and normal fluid velociies v n and v s are imposed in he posiive and negaive x direcions respecively, where v ns = v n v s is proporionalo he applied hea flux. Simulaions are performed wih he same numerical resoluionδ = cm and ime-sep = 10 4 s. The iniial condiion consiss of a se of vorex loops se a random locaions, which are displayed in Fig. 4. During he evoluion we monior he vorex line densiy calculaed as, L = Λ V, (5) where V = D 3 is he volume of he compuaional domain, L is he enire vorex configuraion and Λ = L dξ is he oal vorex line lengh. We also monior he number of reconnecions per uni ime (he reconnecion rae), ζ. We run simulaions, wih he differen reconnecion algorihms, using boh he LIA and he full BS law wih v ns = 0.55cm/s and emperaure T = 1.6K (α = 0.098, α = 0.016). We hen run furher simulaions, using only he

5 5 FIG. 4: The iniial condiions used in every simualion; a se of random loops wih radius cm. L L ag replacemens PSfrag replacemens L FIG. 5: The vorex line densiy L (cm) ploed as a funcion of ime (cm) for simulaions wih he LIA (lef) and he full BS law (righ) wih differen reconnecion algorihms, v ns = 0.55cm/s and T = 1.6K. The symbols are as follows: (magena) diamonds: Type I reconnecion; (yellow) riangles: Type II; (cyan) penagons: Type III; (black) circles: Type IV (ǫ = 10 4 ); (blue) squares: Type IV (ǫ = 10 3 ); (red) crosses: Type IV (ǫ = 10 2 ) and (green) aserisks: Type IV (ǫ = 10 1 ). full BS law, for hree furher values of v ns. This allows us o es he heoreically prediced [8] and experimenally and numerically [40] verified law for he seady sae line lengh, L = γ 2 v 2 ns, (6) where γ is a emperaure-dependen parameer. Our aim is o deermine wheher he differen reconnecion algorihms yield differen values of γ. We sress ha γ is a macroscopic quaniy which is known from experimens.

6 ζ/λ ζ/λ ag replacemens 10 0 PSfrag replacemens ζ/λ FIG. 6: The reconnecion rae, ζ (s 1 ), scaled by he vorex line lengh, Λ (cm), ploed as a funcion of ime, (s), for he LIA (lef) and BS (righ) simulaions; ploing symbols and colors is as in Fig. 5. Noe ha in he case of he LIA simulaions a number of he simulaions reach a degenerae seady sae, as is visible in Fig. 5 (lef). In his case he reconnecion rae drops o zero and hence he resuls are no visible on his semi-log scale plo. This does no happen in he full BS simulaions. IV. RESULTS The resuls for he iniial simulaions using boh he LIA and he full BS law are presened in Figs Figure 5 shows he vorex line densiy ploed as a funcion of ime for he LIA (lef) and he BS (righ) simulaion, using various reconnecion algorihms. The resuls for he Bio-Savar law are very encouraging showing lile difference beween he simulaions wih he differen reconnecion algorihms. In all simulaions here is an iniial rapid growh in he vorex line densiy, followed by a more gradual increase unil evenual sauraion o a flucuaing seady sae. As we discussed earlier, he reconnecion mehod is he one ad-hoc aspec of he VFM, and i seems ha, a leas for quanum urbulence driven by counerflow, he resuls are very robus. Noe infac ha all calculaions performed wih he BS law (Fig. 5 righ) converge o approximaely he same value of L. We believe his is a very useful resul as his means ha one can confidenly draw conclusions from simulaions which use he BS mehod. We noe, however, ha one should probably always err on he side of cauion wih VFM simulaions and check he effec of differen reconnecion algorihms wih a differen numerical experimen. Figure 6 (righ) shows he number of reconnecions per uni ime per uni lengh, ζ/λ, for he BS simulaions, and he picure is he same. For each of he differen algorihms he reconnecion rae per uni line lengh is approximaely he same. The only excepion is he simulaion wih he Type I reconnecion algorihm, which shows a slighly higher reconnecion rae. This is undersandable as his is he algorihm which has he leas resricions and simply reconnecs filamens if he disance beween hem is less han δ/2; however here is no discernible impac on he evoluion of he vorex line densiy. Snapshos of he sysem a he end of simulaions are ploed in Fig. 8. We shall reurn o discuss he furher simulaions using he BS mehod a he end of his secion. The picure for he LIA simulaions is very differen. We find ha even wih he same iniial condiions, sligh differences in he reconnecion algorihms lead o large changes in he overall evoluion of he sysem. Kondaurova and Nemirovskii [47] found ha using a reconnecion scheme similar o he Type IV algorihm used here hey reached a seady sae. We find ha cerain simulaions do reach a seady sae, alhough he seady saes we find are degenerae. By degenerae we mean ha he sysem evolves o a se of very sraigh vorices, arranged in planes parallel o he counerflow direcion. These planes are hen simply adveced in he direcion of he counerflow wih a negligible change in vorex line densiy. In his configuraion he reconnecion rae drops o zero. This is he evenual fae of he simulaions wih Type I, II and IV (ǫ = 10 3,10 2 ) algorihms. In he oher simulaions a seady sae is no found and large flucuaions in he line densiy are seen. We find he same bundles of vorex lines, sraified in layers, observed by Schwarz [15] and recenly reproduced by Adachi [40]. Schwarz recognised ha his is a spurious effec and inroduced his arificial mixing procedure o avoid i [15], see Fig. 7 (righ). The reconnecion rae per uni line lengh, Fig. 6 (lef), also show ha LIA simulaions are highly sensiive o he reconnecion algorihm used. We hus agree wih oher auhors ha, due o he absence of he non-local componen of he velociy field, LIA is unsuiable for simulaions of urbulence[40]. We now reurn o he BS simulaions, focusing on he affec of varying he counerflow velociy v ns. We are

7 7 FIG. 7: Snapshos of he vorex filamens from he simulaions using he LIA, as ploed in Fig. 5 (lef). The lef panel shows he ypical sae of each of he simulaions a he peak of vorex line densiy, a 6(s); his paricular plo is aken from he simulaion using he Type III reconnecion algorihm. The middle panel shows he degenerae soluion aken a = 80s for he simulaion wih Type IV reconnecion algorihm (ǫ = 10 3 ). The righ panel shows a plo in he plane perpendicular o he counerflow for he simulaion wih he Type IV (ǫ = 10 1 ) algorihm a = 60s: his layered vorex srucure was recognised by Schwarz [15] and Adachi e al. [40]. FIG. 8: Snapshos of he vorex filamens in he plane perpendicular o he counerflow for BS simulaions, from Fig. 5 (righ), each a = 10s. The lef panel is for he simulaion wih he Type II reconnecion algorihm, Type III (middle) and Type IV (ǫ = 1) 4 (righ). ineresed in he relaionship beween he seady sae line densiy and he counerflow velociy. As above, he emperaure is fixed a 1.6K Figure 5 (righ) is very represenaive of he ime series plo for he line densiy, a differen emperaures, and so we do no reproduce furher plos for each of he differen values of v ns. Insead in Fig. 9 we display he linear relaionship beween L and v ns, according o Eq. 5. Again, we find ha he resul is no sensiive o which reconnecion algorihm is used. Finally i is insrucive o compare he effec of he reconnecion algorihm used on he scaling parameer, γ, and compare wih experimenal resuls. Our aim here is no o claim ha he algorihm which yields he value closes o ha experimenal value should be used in he VFM. Insead we wish o add furher weigh o our claims ha resuls from he VFM can be informaive, so long as a physically jusifiable reconnecion algorihm is used. Table 1. shows he values of γ for he differen reconnecion algorihms. The experimenal value of γ a his emperaure is 93 [10, 45] ; Adachi e al. [40] gave numerical resuls of γ =

8 L PSfrag replacemens v ns FIG. 9: A plo of he roo of he vorex line densiy, L, and he magniude of he counerflow velociy, v ns. The linear relaionship beween L and v ns for each of he reconnecion schemes is apparen, ploing symbols and colors is as in Fig. 5. reconnecion algorihm γ Type I Type II Type III Type IV (ǫ = 10 4 ) Type IV (ǫ = 10 3 ) Type IV (ǫ = 10 2 ) Type IV (ǫ = 10 1 ) TABLE I: A able of values of γ for he differen reconnecion algorihms in simulaions using he Bio-Savar law a T=1.6K V. CONCLUSIONS In conclusion we have shown ha he vorex filamen mehod is very robus o he reconnecion algorihm used, provided one uses he full Bio-Savar inegral o deermine he velociy field. In conras resuls from simulaions using he local inducion approximaion are highly dependen on he algorihm used, lending furher weigh o he criicism of his mehod. VI. ACKNOWLEDGEMENTS We hank C. F. Barneghi for useful discussions. This work was suppored by he Leverhulme Trus. [1] R.J. Donnelly, Quanized Vorices In Helium II, Cambridge Universiy Press, Cambridge (1991). [2] C.F. Barenghi & Y.A. Sergeev (eds.), Vorices and Turbulence a Very Low Temperaures, CISM Courses and Lecure Noes, Springer (2008). [3] W. P. Halperin & M. Tsuboa (eds.), Progress in Low Temperaure Physics: Quanum Turbulence, vol. XVI, Elsevier (2008). [4] J. Maurer & P. Tabeling, Europhys. Le. 43, 29 (1998). [5] P.-E. Roche, P. Diribarne, T. Didelo, O. Français, L. Rousseau, & H. Willaime, Europhys. Le. 77, (2007). [6] M. Bla zková, D. Schmoranzer, L. Skrbek, & W. F. Vinen Phys. Rev. B 79, (2009). [7] M.R. Smih, R.J. Donnelly, N. Goldenfeld, & W.F. Vinen, Phys. Rev. Le. 71, 2583 (1993).

9 [8] W.F. Vinen, Proc. Roy. Soc. A240, 114 (1957); W.F. Vinen, Proc. Roy. Soc. A240, 128 (1957); W.F. Vinen, Proc. Roy. Soc. A242, 494 (1957); W.F. Vinen, Proc. Roy. Soc. A243, 400 (1957). [9] M.S. Paolei, M.E. Fisher, K.R. Sreenivasan, and D.P. Lahrop, Phys. Rev. Le (2008). [10] J.T. Tough, Progress in Low Temperaure Physics, vol. VIII, Elsevier (2008). [11] G.P. Bewley, D.P. Lahrop, & K.R. Sreenivasan, Naure 441, 588 (2006). [12] T. Zhang & S.W. Van Sciver, Naure Physics 1, 36 (2005). [13] D.I. Bradley, S.N. Fisher, A.M. Guénaul, M.R. Lowe, G.R. Picke, A. Rahm, & R.C.V. Whiehead, Phys. Rev. Le. 93, (2004). [14] W. Guo, S.B. Cahn, J.A. Nikkel, W.F. Vinen, & D.N. McKinsey Phys. Rev. Le. 105, (2010). [15] K.W. Schwarz, Phys. Rev. B 38, 2398 (1988). [16] D.C. Samuels, Phys. Rev. B 46, (1992). [17] A.T.A.M. de Waele and R.G.K.M. Aars, Phys. Rev. Le. 72, 482 (1994). [18] C.F. Barenghi, D.C. Samuels, G.H. Bauer, & R.J. Donnelly, Phys. Fluids 9, 2631 (1997). [19] M. Tsuboa, T. Araki, & C.F. Barenghi, Phys. Rev. Le. 90, (2003). [20] V.B. Elsov, A.I. Golov, R. de Graaf, R. Hänninen, M. Krusius, V.S. L vov, & R. E. Solnsev, Phys. Rev. Le. 99, (2007). [21] L. Kondaurova & S.K. Nemirovskii, J. Low Temp. Phys. 138, 555 (2005). [22] D. Kivoides, J.C. Vassilicos, D.C. Samuels, & C.F. Barenghi, Phys. Rev. Le. 86, 3080 (2001). [23] K. Morris, J. Koplik, & D.W.I. Rouson, Phys. Rev. Le. 101, (2008). [24] D. Kivoides, Phys. Rev. Le. 96, (2006). [25] A.W. Baggaley & C.F. Barenghi, Phys. Rev. B 83, (2011). [26] A.W. Baggaley & C.F. Barenghi, Phys. Rev. B 84, (2011). [27] D. Kivoides, C.F. Barenghi, & Y.A. Sergeev, Phys. Rev. B (2008). [28] A.P. Finne, T. Araki, R. Blaauwgeers, V.B. Elsov, N.B. Kopnin, M. Krusius, L. Skrbek, M. Tsuboa, G.E. Volovik, Naure 424, 1022 (2003). [29] P.G. Saffman, Vorex dynamics, Cambridge Universiy Press (1995). [30] C.F. Barenghi, R.J. Donnelly, & W.F. Vinen, J. Low Temp. Phys. 52, 189 (1983). [31] R.J. Donnelly & C.F. Barenghi, J. Phys. Chem. Reference Daa 27, 1217 (1998). [32] L.S. Da Rios, Rend. Circ. Ma. Palermo 22, 117 (1905) [33] R.J. Arms and F.R. Hama, Phys. Fluids 8, 553 (1965). [34] M.S. Paolei, M.E. Fisher, & D.P. Lahrop, Physica D (2010). [35] J. Koplik and H. Levine, Phys. Rev. Le (1993). [36] R. Tebbs, A.J. Youd and C.F. Barenghi, J. Low Temp. Physics 162, 314 (2011) [37] R.M. Kerr Phys. Rev. Le. 106, (2011) [38] M. Kursa, K. Bajer, & T. Lipniacki, Phys. Rev. B 83, (2011) [39] R.P. Feynman, Progress in Low Temperaure Physics 1 17 (1955). [40] H. Adachi, S. Fujiyama, & M. Tsuboa, Phys. Rev. B, 81, (2010). [41] S. Fujiyama, A. Miani, M. Tsuboa, D.I. Bradley, S.N. Fisher, A. M. Guénaul, R. P. Haley, G. R. Picke, & V. Tsepelin Phys. Rev. B, 81, (2010). [42] M. Tsuboa, T. Araki, & S.K. Nemirovskii, Phys. Rev. B 62, (2000). [43] M. Tsuboa & H. Adachi, J. Low Temp. Phys. 162, 367 (2011). [44] M. Leadbeaer, T. Winiecki, D.C. Samuels, C.F. Barenghi & C.S. Adams, Phys. Rev. Leers 86, 1410 (2001). [45] K. Childers & J. T. Tough, Phys. Rev. B 13, 1040, (1976). [46] H. Adachi, & M. Tsuboa, J. Low Temp. Phys 158, 422 (2010). [47] L. Kondaurova & S.K. Nemirovskii, J. Low Temp. Phys. 150, 415 (2008). 9