Self-consstent feedback mechansm for the sudden vscous dsspaton of fnte-mach-number compressng turbulence arxv:1805.1336v1 [physcs.flu-dyn] 31 May 018 A. Campos, B. Morgan Prevous work (S. Davdovts and N. J. Fsch, Sudden vscous dsspaton of compressng turbulence, Phys. Rev. ett., 116(105004), 016) demonstrated that the compresson of a turbulent feld can lead to a sudden vscous dsspaton of turbulent knetc energy (TKE), and suggested ths mechansm could potentally be used to desgn new fast-gnton schemes for nertal confnement fuson. We expand on prevous work by accountng for fnte Mach numbers, rather than relyng on a zero-mach-lmt assumpton as prevously done. The fnte-mach-number formulaton s necessary to capture a self-consstent feedback mechansm n whch dsspated TKE ncreases the temperature of the system, whch n turn modfes the vscosty and thus the TKE dsspaton tself. Drect numercal smulatons wth a tenth-order accurate Padé scheme were carred out to analyze ths selfconsstent feedback loop for compressng turbulence. Results show that, for fnte Mach numbers, the sudden vscous dsspaton of TKE stll occurs, both for the solenodal and dlatatonal turbulent felds. As the doman s compressed, oscllatons n dlatatonal TKE are encountered due to the hghly-oscllatory nature of the pressure dlataton. An analyss of the source terms for the nternal energy shows that the mechancal work term domnates the vscous turbulent dsspaton. As a result, the effect of the suddenly dsspated TKE on temperature s mnmal for the Mach numbers tested. Moreover, an analytcal expresson s derved that confrms the dsspated TKE does not sgnfcantly alter the temperature evoluton for low Mach numbers, regardless of compresson speed. The self-consstent feedback mechansm s thus qute weak for subsonc turbulence, whch could lmt ts applcablty for nertal confnement fuson. 1 Introducton The compresson of a turbulent flow occurs n a broad array of applcatons. Examples nclude onedmensonal compressons n nternal combuston engnes [1] or across shock waves [], axsymmetrc compressons n Z-pnches [3], sphercally-symmetrc compressons n nertal confnement fuson [4, 5], and three-dmensonal complex contractons n the nterstellar medum [6]. Moreover, the compresson mechansm often leads to complex turbulence dynamcs, and the resultng evoluton of turbulence can have a strong effect on the overall behavor of the applcaton under consderaton. Thus, ncreased levels of understandng and mproved modelng capabltes for ths phenomenon are essental. Numerous drect numercal smulatons of compressng turbulence have been prevously carred out wth the am of mprovng engneerng turbulence models; see for example [7, 8, 9, 10, 11]. These studes treated the flud as a tradtonal gas, for whch the dependence of vscosty µ on temperature T s gven by µ T n, wth n havng a value of, or close to, 3/4. On the other hand, [1] demonstrated, through computatonal smulatons, that when a power law exponent 1
characterstc of weakly-coupled plasmas s used,.e. n = 5/, a sudden vscous mechansm occurs whch dsspates the turbulent energy. Ther results showed that a turbulent feld subjected to a contnuous sotropc compresson ntally creates an amplfcaton of turbulent knetc energy (TKE), untl vscous scalng domnates and TKE s rapdly dsspated nto heat. It was thus proposed n [1] that the resultng ncreases of temperature could be used to mprove the gnton condtons for nertal confnement fuson. Subsequent work has expanded on the smulatons of [1]. The effect of onzaton on the scalng of vscosty was accounted for n [13]. For that study, the onzaton state Z was assumed to depend solely on temperature, and thus the plasma vscosty µ T 5/ /Z 4 was smplfed to the form µ T β. Ther analyss of the evoluton of the energy spectrum showed that the sudden dsspaton of TKE occurs for β > 1 only. A TKE model that accounts for the vscous dsspatve mechansm for sotropc compressons s presented n [14]. Ths model was valdated aganst drect numercal smulatons, and showed excellent agreement for vscosty-power-law exponents greater than one. The model was then used to estmate the partton of energy between the turbulence and heat, as the compresson proceeds n tme. A two-pont spectral model based on the EDQNM formulaton was used by [15], along wth drect numercal smulatons, to reproduce the sudden vscous dsspaton mechansm. The lower computatonal cost of the EDQNM model allowed for the analyss of hgh- Reynolds-number effects and thus the dentfcaton of three dstnct regmes: turbulent producton, non-lnear energy transfer, and vscous dsspaton. Moreover, the assumpton of homogeneous turbulence was relaxed and a sphercal nhomogeneous turbulent layer under compresson was smulated wth both DNS and EDQNM closures. The sudden dsspaton of TKE was also observed for ths new case. Fnally, n [16], a stablty boundary for hot spot turbulence was derved to demarcate states of the compresson for whch a decrease of TKE s guaranteed. Moreover, an upper lmt for the amount of TKE that can be generated durng a compresson was proposed. Ths upper lmt was then compared to the nternal thermal energy of the system. The smulatons of the sudden vscous dsspaton mechansm prevously conducted have reled on the zero-mach-lmt assumpton. Gven a decomposton of the velocty U = U + u, where U s the Reynolds-averaged mean flow and u the fluctuatng velocty, the governng equatons for the fluctuatons n the zero-mach lmt take the form [7] u = 0, x (1) ( ) u ρ t + u u j = P + µ u + f. x () In the above, ρ s the Reynolds-averaged densty, P the pressure, µ the vscosty, and f a forcng functon that accounts for the effect of the compresson. The vscosty depends on temperature and thus an a-pror tme evoluton for temperature needs to be provded. For an adabatc sotropc compresson, ths s T = T 0, (3) where s a characterstc length of the doman beng compressed and T 0 the ntal temperature. The approach descrbed above s sutable for demonstratng the sudden dsspaton of TKE, but does not capture the self-consstent feedback mechansm mentoned n [1]. The feedback mechansm begns wth an energy transfer from the TKE towards the nternal energy, as a result of the sudden vscous dsspaton. Ths, n turn, causes ncreased temperatures that amplfy the
vscosty of the system. The stronger values of vscosty then precptate the vscous dsspaton of TKE, thus completng the self-consstent feedback loop. In the zero-mach lmt, an evoluton equaton for the nternal energy s not solved, and thus the effect of the dsspated TKE on the nternal energy and the vscosty cannot be reproduced. It s expected that accountng for the feedback mechansm would lead to vscous dsspatons that are more sudden, or of ncreased ntensty [1]. An alternatve to the assumpton of the zero-mach lmt s turbulence belongng to the fnte-mach number regme. For ths case, fully coupled governng equatons for densty, velocty, and energy are solved, whch allows for an explct accountng of the forward transfer of dsspated TKE nto heat, and the subsequent effect of ncreased temperature and vscosty on the dsspaton. The focus of ths study s the smulaton of turbulence n the fnte-mach number regme to nvestgate the complex self-consstent feedback mechansm, and thus further assess the benefts of vscous dsspaton for nertal-confnement fuson and other hgh-energy densty applcatons. The outlne of the paper s as follows. Secton ncludes a descrpton of the governng equatons for turbulence n the fnte-mach number regme. The mechansms that account for the energy transfer between the TKE and the nternal energy are also dscussed. In secton 3, detals of the numercal smulatons, such as the dscretzaton scheme and the creaton of realstc ntal condtons, are ncluded. The results of the smulatons are then provded n secton 4, whch s dvded nto two subsectons. Secton 4.1 focuses on the component of the feedback mechansm related to the TKE. Thus, the evoluton of the TKE, ts budget, and spectra are analyzed n ths subsecton. The component of the feedback loop assocated wth the nternal energy s then nvestgated n secton 4., where an analyss of the temperature evoluton and sources for the nternal energy are ncluded. Fnally, the paper ends wth secton 5, where concludng remarks and a dscusson of future work s provded. Governng equatons.1 Naver-Stokes equatons for sotropc compressons We denote Ũ and u as the Favre-averaged and Favre-fluctuatng veloctes, respectvely, so that U = Ũ + u [17]. In analogy to the zero-mach-lmt formulaton of [1], we analyze the effect of a compresson on a statstcally homogeneous turbulent feld u, where the compresson s acheved through a specfed Favre-averaged mean flow Ũ. The Favre-averaged velocty for homogeneous compressble turbulence needs to be restrcted to the form Ũ = G j x j [9]. The deformaton tensor G j correspondng to an sotropc compresson s gven n [7, 9, 1], and can be expressed as G j = δ j, (4) where s a tme-dependent characterstc length of the compressed doman, and s the constant tme rate of change of. Gven ths formalsm, one can derve, as detaled n appendx A, a set of Naver-Stokes equatons for the fluctuatng velocty undergong mean-flow compresson. These equatons, whch are summarzed below, consttute the fnte-mach-number analog of the low-machnumber eqs. (1) to (3), ρ t + ρu = f (ρ), (5) x ρu t + ρu u j 3 = σ j + f (u), (6)
ρe t t + ρe tu x = u σ j + ( κ T ) + f (e). (7) In the above ρ s the densty, u the Favre-fluctuatng velocty, and T the temperature. E t s the total energy, whch s gven by E t = U + K, where U = C v T s the nternal energy and K t = 1 u u s the knetc energy assocated wth the turbulent fluctuatons. C v s the specfc heat at constant volume. The stress tensor s gven by [ ( ) 1 u σ j = P δ j + µ + u j 1 u ] k δ j, (8) x 3 x k where P s the pressure and µ the vscosty. A power law of the form µ = µ 0 (T/T 0 ) n s used, where µ 0 and T 0 represent reference vscosty and temperature values, and n s the power-law exponent. The thermal conductvty κ s computed accordng to κ = µc p /P r, where C p s the specfc heat at constant pressure and P r the Prandtl number. An deal equaton of state P = ρrt s used, where R s the deal gas constant. The forcng functons f (ρ), f (u), and f (e) account for the effect of the mean compresson on the densty, velocty, and total energy, respectvely, and are defned as f (ρ) = ρ, (9) f (u) = 3 ρu, (10) f (e) = [ ρe t + ρu u + 3P ]. (11) The equatons above are sutable for numercal smulatons, snce the compressve effect of the mean flow Ũ has been abstracted nto the three forcng functons above. These are the equatons solved for the current study.. Energy exchange for compressble turbulence For turbulence n the fnte-mach-number regme, the Helmholtz decomposton s often employed to express the fluctuatng velocty as u = u (s) + u (d), where u (s) and u (d) are the solenodal and dlatatonal veloctes, respectvely. The solenodal component satsfes u (s) = w and u (s) = 0, where w = u s the vortcty vector, and the dlatatonal component satsfes u (d) = 0 and u (d) = d, where d = u s the dlataton. Gven ths decomposton, two TKEs can be defned. These are the solenodal TKE the dlatatonal TKE k (s) = 1 u(s) u (s), (1) k (d) = 1 u(d). (13) u (d) There are two addtonal energes n the system, namely, the mean knetc energy K = 1 ŨŨ, (14) and the mean nternal energy Ũ = C v T. (15) 4
Table 1: Sources n the evoluton equatons for the solenodal, dlatatonal, mean, and nternal energes. Upper scrpts α stand for ether s or d. τ j represents the Favre-averaged Reynolds stresses (τ j = ũ u j ). Name Symbol Defnton Intermode advecton IA (α) ρu u j (α) ρu + Producton P (α) 3 ρ k(α) G Solenodal dsspaton ρ ɛ (s) µw w Dlatatonal dsspaton ρ ɛ (d) 4 3 µd Pressure dlataton P D P d Mean knetc energy advecton AD ( K) K ρ Ũj Mean knetc energy transport T ( K) (Ũ ρ τ j + Ũj P Mechancal work MW P G ρu u(α) d ) A dervaton of the governng equatons for the solenodal and dlatatonal TKEs s gven n appendx B. Along wth the evoluton equatons for K and Ũ, one can summarze the governng dynamcs of the four relevant energes for homogeneous turbulence as follows ρ dk(s) dt ρ dk(d) dt = IA (s) + P (s) ρ ɛ (s), (16) = IA (d) + P (d) ρ ɛ (d) + P D, (17) ρ K t = AD( K) T ( K) MW P (s) P (d), (18) ρ dũ dt = MW + ρ ɛ(s) + ρ ɛ (d) P D. (19) Each of the sources n the evoluton equatons above s defned n table 1. We note that the dervaton of the evoluton equatons for the four energes assumed a generc yet sotropc deformaton tensor G j, and neglected the averaged heat flux snce for homogeneous turbulence the averaged temperature s unform n space [9]. The ntermode advecton represents a transfer of energy from the solenodal and dlatatonal modes, and thus satsfes IA (s) = IA (d). The producton terms transfer the compresson energy stored n the mean flow to the solenodal and dlatatonal TKEs. The solenodal and dlatatonal dsspatons then transfer energy stored n the solenodal and dlatatonal felds nto heat. The pressure dlataton represents a two-way energy transfer between the mean nternal energy and the dlatatonal TKE only, and the mechancal work transfers energy of the compresson drectly nto heat. The mean-knetc-energy advecton and transport are not dentcally zero, unlke the case for the three other energes. All of these energy transfer mechansms are depcted n fg. 1. We note that each energy component has a drect nteracton wth each of the other three energes. We also note that the drver for the nteractons s the mean knetc energy, snce t has a predetermned tme evoluton that emulates the compresson of the system. The other three energy components then respond n a self-consstent fashon to the tme-evoluton of 5
<latext sha1_base64="g/ji+xzxcjx31o94xgg83a/8y=">aaacrhcbzddshtbfmdnu1uj9pas+gbkexenzt5g7arefkwhxm0ahnnjthkdmazmrwwze/7n3vn/azfavxfvp7q0vxtg4m//fm35hqlnrjutdn4sfty1xjzzxxt9zu371rrgz+0tbufn0ou1sgkgjgt4btmoiwsbsqooqzdxdcypzwhpzkupyzybktuwaro8rya9pqbxseacxepp8groajo8gwoma7wtq86yytjpu160cpxdetqojuppurmwzcrnyzuzcq1hmtm8mjsom5fktbqehdehmmzskbj0jk+okfcwdwy4kso+yxdl/turk1jra5tzuzmmx6ak83/5capxynornjakdqh0vryrgrucsdz0wbntyzgldfsgj0jchcz7hw0ipf4ael4yg8gcqz6zx1niuyjy7wdftzd7n/e9j59pf5f6knbpr3ajhw/udze91b1zvpdfzao5ntina6apyub/toun0jhxe0u/0c1gs+fkuxfunbuh0ozt97xhj8k5/w0cabdp</latext> Mean Knetc Energy e U e K = 1 U <latext sha1_base64="zmsotbqkycss4m9jnpkkmcfdk=">aaackxcbvdtgjbejz1fhej14ookmhoza7krpegrjygg0mwovtbhzmczmtcnvsvxvup/bpp6tufcccnlpjj5xq7nr1esgjsjvoqzuxotu9m5uzy84vc4tr+rw5wijie6euzikw8rydsauqaawvuoaxopqcprhw/7jtuqorgqg9cahpcdahpcdzgut6s3ctbafzvf3jo0vnbpsvcvosrylqtzlrvturjoiqamkwuk3xcxu7xljtwdjtifisz93iwmoqhmonrxyhfbxmly/tcmgq0pvj/bssykzxgnpnkwpfu795q/k/xjr/i5peeyaavj5yi+yryu9fn/uuales4ehmehqvnqkhyumoq0dsutoq+xp8wnhfbnvcotmr94iikpj85smknud4oupu7ynp/u6tyldes0nf/yqvxpwu5zkv/ds3pmohw0gbar8qok5rddurqrzdowf0ad1zz9a9fy5omglotodnb7bkepzm=</latext> <latext sha1_base64="rlcmkgg0eqro53ztzh7wo55bc=">aaacvhcbvddstxagj3evpxu9r3gxdbc9ksbbqdu/ei30xlzxhcyzey+rmnozs8cchkpxya3ty3kh0yozm1tnxwmcz8/0fwghubbt9dmklz8vpx6ysrq1vvnx81dnavjdkaqyjpotsv5qyefzcyhir4krqqhiq4jopztxyxvqht5bsscjjmzsz5xrqyxpp1+wkba0fzoqjmlqcqeeohpnabo/yxoevzz97ac4rmjelbw00436ku4h9j3jiuane63qrwk1qxl/t+tbbjxnfuelfto8qof5nigcsdmzwdhtsxiwkpxxi13vjtgn//d49e+kutgldn1mtmx1+zprbh/fxs7m7fvfwwzojkd4myj7bvuhekp1wds60hddngwfyndgeebset6fkzshxl+squbdwxz5suvdurq68fjuorbod78cg7xr74jtdj8f+s4bxb/tg/d6wf3/pd4/wx9x0bv0fuga3qefqetteimxsvqhv6c74edyhs+hyqoytdwv0soem78atdy0+g==</latext> P <latext sha1_base64="ukqg8w7prmnzd9/zx+qewj5he=">aaacwxcbvfnbxmxehwwrztli6vhhypeocqk0aq5ifrxa4liqsytlo8jzrzwvouvpvsuwf4p/tw9wg9a6o/bsvdac0+y9obnjgbmmddkwozjn53ozt179x90n3qbdx89ftfevrv6syisivwpxwzkhjclkuqocknsbkrucy96t8sfnykzu1rdc1jarwklcykybmneh+9knbn30dm3ncsmey7xburp9k3mgfl4fi/l7fjnxl/ea/jne/jgnjgq4mm91nrnc6aecov1k6tumazywaluob7wwohzmbtmeaamvksdo3vtdtf0hjaafnebxstfp3holtcush8qs4zm9nvuj/8tngywozk5wdynqetbramoarqy+bsgec1diqji8ouvjyx4bmgu9m4vovkhebsmhcn+xxtkvs64a8o7z+7fexbvj17vjqz0f7ev+udyjgo+3zj8t8ddsfd5nnocjf+tglz8hz8pikzewoyqdyrfiyaw5jn/jj85vfexdqhddgnxanm1ya9hy+ltv4=</latext> Solenodal Turbulent Knetc Energy P (s) (d) k (s) = <latext sha1_base64="uvfvvdcgaecy4c7fv+vbohoby=">aaackxcbvdssnafj3u9/u1dbosqgupa/anehcl1wsvtpajpobonqmeymqgn5crf6b36nk3xrjznqa164mh3hu55x4vylrpx3m1cojy+mtk1ptm7nz8wus8uxssssqj0ijmtdwwoydaguqwbqcrg7jg48nrhg/7vhuhfrxh+xg0oe6gnkaeayndr9dukpk/ma53fotoclg/yvutooor6zzm0fefdqemdcvvdj1itxmsnsum0ulwrcdcpie70dq0xbxuo8kcp/aguxw7enjuqo1m/bmryk5un3tmkmn9q373buj/vwasg4nqsmo1hcsz0nbzgwt7mh7tk8lem36hmaqffqk1ssmdempkernha9jt1lpvehmpfmjhpe91gmaxj+cpqmzlzlf8vd0lnfo53alwdym5oqu5xfpxtcrxsnm0xd3fzhcfrklpdjescjpep6g6oggju7ra3q0nqxn68v6+xwtwpnochqc9f4bsf6lha==</latext> 1 ] (s) (s) u u AD(K ) + T (K ) <latext sha1_base64="j65k6ulhhpiwexjkwlevzwhpi=">aaacqncbvdahsxfnwkzan5osyg1e74bbjzvjyvusfkeiabgtgo0ysxphedambkltmgk+zpukz/or/qburdzlhzgvi7zqhb0tn3cu89nbvcg9//6c9edm/sj0anlldw19o7t5+lztdfomymkuqreg+ajta03a5tbssma7o8opyv/ggsnoztmwohv5mbgmpocpgsf0srrz/dgwrxuqupc13cryw9p/sr9u9uv+bph/ehskjaqc9te9lw4owrzdypggwnccpzu9s5thtec+3m00piqnyqa6jykbtzk1nyvoueedsuzcypfgnoyyjtr7f1fxgxfzrp95yfm7rzcz617m8stmdcxsyfgucg4nhuecqkbgjbwhthg3kbxrbfmxhozu6uq0oodpeeencj0vaeerfbrn9eviht36tcvqdvdzf1/w7wmq+ogim3gmb7xr1zd77sly8phyx0bv0flvrgbrogh1gz6ngpb39etuvn+e+8396fh9i5r+jzqjpw7v8czpkueg==</latext> MW IA <latext sha1_base64="la+ch3c4mw0xz0sz+x4eabwgznc=">aaacjhcbvbns8naen34wvs/wj16cvbbqyljrdbeh70ilsxvqlbatdukmg3y3qgn5c171h/hrpikh/4wtlqqz4yeyzyd48jbukst6n+bmfxaxljmrdzqvpgvrb5i3kocqiz1zcolgcoz60ffumbgmbhmytj3r6atfvgchkfev1tarochpnupwwo7v60d3v5olweph/zsikrr7bwm3ffk9adxxggpezyvqc6erakegzx96uegaywgpoaopjdq3sszg5p5w+qbhs5fmyn6cypcnprjz9gthlzd+bs3ef/rdulhncj6gehap9md7khmxu3j5+bfsqakdbwbbnbtvetdharol8zq44ni8udqt+pa+utjyxgxljh4uqr1dnj3fklwphjfvwinbxd+p1tvq1ztu7a/4wpvyvwxfvoqnappjbmjhbspbfrddxsomqfcbqb/sinoxn48v4nd6mo3ngurofzmb8fakkg6ob</latext> Mean Internal Energy e U = Cv Te <latext sha1_base64="cvasv3m8kyvqtqrx3zt+r71c4=">aaacrxcbvdshxbfk0ro3gv7j0uzgiwcjq7wsyo4kedsgoofgyhurzexmoqqpqo60dt9afkat8kf+av+rfbvqvhxvebglpnvk+ucm6s7197yx/gp05mtnanvk8oze/spjll1gzzthhsh9fofbwsvccz1knkeqct6phthu//y3accvpbj54gb5dfnyj10vtaiguqmstbcyag6x71kybognkfevug1w3/rhog9jujmgqxf4vujnhh3fsss1zgkm6qz+ansfamuzwhi6zaymwiywwk6jehi0vwj0tulxnnknsdusuth6vokahjj8rymqnyc/m6vonvxbqzjb/3c7tzkjkd4ptfcragun7xonzircewcau10puwanzdnyckqk1nbcznwlfyyd9amlysprizlcbz7oyz81dbwarc5xtkxpnwdmht1kqdpnrxwwum8hrwmn7o/zxyrzfknbkrftskeosqdwsgfckn+kn/elfffu/fuh1hvrmk3kb7+4eauyxaq==</latext> <latext sha1_base64="q5qcu1h8ojhiemrjwsatrtbwks=">aaacjhcbvbnt8jaen36eah6nfi5p4ikrffmhhtqbglesigs7tghdttvsbk1i07/gvf+bv8at8edf3+icjyr6kkle3sxk3jwnyfqqy3o35uyxfpewuyvpzora+kyt3kresginahnxdqdifrhxqkkgbnqadhaz3zvbs3+7byep9/ukicoh/s+dsnbaztxp7sdrn5qhnyp4ldkykeg9mzmy7r4noqe+igxbktqhubqldoj0o5qqydefwhp6mmpzceami3npa30tjcxb4yj+rpjqh7uo48r096wa3k795y/k/xcpv73imoh4qkfdi95ibmvnwcf7qacegttatvxk0ywaitpfozuejwpltykfqthrg6ynjybs5iq4uj4/jsorudkqichsy7p/k6tqlepj6rqf8xxkbwrrfuqlk+vkxxtabvtoh1kowqqoqturw1e0aa9oef0zdwb8ar8tydntosns00a+pje1zo4e=</latext> h (s) <latext sha1_base64="rbw5rkcbqb+apn3jjwsfmuyrdu=">aaacr3cbvbtsmwehxkr/wnlyfcsqujwut+muggusavfaakvlujpwqmnhtonurtkbpen+ag3iivyolbiv5psvxmzyxm5vkrz9q47rotgxufmjzkt8/mzs0vbawl+1jbwfopvcqmufaobmqn0ww+e6ukbcn8ov3z8a5q9uqwkmxyuzrnaksvewgffrnqubkw3orfddrpetj+wqcrzlykmrtb6xr7urbkj4/ey0grzthtzmzzy6kcqjcue60bnhuzfojuyzrdulmm9yqedonxwhykkgiupwm7knxhlu6ojdkpmhwsp3ekzbq60ho8qqmj7+nruk/+uasqkowgktuwxa0i9bqcyxkxhodu4wbdtwgswekmz3xbrhfkhgwvhj9l3xbf0s6efjzr0sm4cdnmatj+ffhmrjblf1tb8ntfz5x65vsq7gal6n/bvy6vqytsrfu7my95tirw0cbyuaxv0ak6rxve0r6rw/o0xlyxpxx5+jnodkpsvob3oo3bssiw=</latext> h (d) <latext sha1_base64="pxbqnjykqsjhvchr55s5c1+pm+m=">aaacr3cbvbtsmwehxk/19gycaiigeqqr8sncifraeragpzxjtfqrjh3zdliv5qsch3cgbtwc1aijw6j+d/j1ps3m5qz58ecaeo6z05hzhrsfgjyanpmdm5+ob4dkllojuqersxfteacc6oyzdtexah5hk783tegf3usjmpkw/hlzeookfjbjjpxzxfa3jehwwe3vlfba4g141cpbvbzrbwpbcsjse/ku8njrqjtpojptdcrnihcgcqj1w3nj00qjmoxyykabyay0b7pqmnsqsqrxr4t4bxrrgucr7hmfd9xthskt+5fvkynuvp3bd+lskjtxvpuzeqfbpwafccdg4oe5ogakqof9swhvzo6kazcoqo18mcux8qeib6lwqqwvohs6zr3+cjzon58wgwulvvvs1vdzuz9nlfrtusqjs5qxmf9tur5vrzo6uudnzyhyfrclpfg8hdvxsattapqok7ta9ekcpzppz4rw6bx+lbsfvwuy/uhdeavyssh0=</latext> <latext sha1_base64="t9wzjsmmdxy5pbkbvxy481m6go4=">aaacrxcbzdpbtnaembxauqbqp/qyy8reqouws3hzadutueyjffnymurnf6pu5wwxut3xfrzpkbebqu7rv0ffosnfcvdjnakacftnknby1m78glckg6945lsdpn0839yq1l683n7z3xv1xahmc/c5kkr3azagr8fchnpgcschf8zofvxefwgjvhkj8xrgmzskihkcoy3gu/xg8ksiayumifc+pgeja/osnzzngyx3kxov8arzr1uqo73nnqw1dxiyeuwtgddw3xvhonaouoagomwmp4zmgyg1cyvbjplnqv9y5oqrkrblybdpqs/chybm48dxkznjq/a4vwf7vbhtghus6sneni+oogkjmufvgoahqwfhorwfcc7sr5vomguccg1konqmwwdsjsfefv1yytyebzkdiv98flrkbrp9vum9oyospu8pty5v+7g0he83pv+w1wl5f4f1k+os4yy5ik/jw+krnjkhn0x+istb+q7usy3zq3zw/np3d+vpzyzz5zk/prafxsno=</latext> PD Dlatatonal Turbulent Knetc Energy <latext sha1_base64="kq05oku0zy43rjobegvmes4=">aaacvncbvdbtnafb0bskvi4ulmxereosqsttayk7izdyfntqsneu3zlcp6omz6x5ifmwf4svyqt/0k+bswsbby400te5d+7jsfik65kmopv3ny6tb1ze3fvzt179wf7dz5z7q3hgddsm3mgfqvqohpcstwvdubjj6xzavwp/umxgqttl1v4qkatrk54occtbymm47korfqxb8welyng6sn3javg0ezj4p3z3g9i1cs66a5wcyjjio9f+s9mriepws96o9bkw50i/shaezqublgdcv0lnruzt1gc38aa54eqknau6u68hj4jyormoqx9unup3/uufhbfsxufuau7n9ek/7pm3uxv1juqpxeoejxg3ivqdo0zyquhehuzbuiccpabpgfgoahsottmnybb8ygs1ayh7s6jeuytjju398+7kpk4pj84p0vebx/o7twnaznytaforvtnhadpkpxwoj8d9jjvkexlmnpkutmgxeudoyixw8ov8jd/i9+gy+hfvxdtxpxhu/3liree/aqj37xa</latext> k (d) = <latext sha1_base64="uccdow8wogmzmeje7ycphjcfio=">aaacjhcbvdssraejz4xn+uhr0ev8gdimvdw+gh5xmsqsqsazhxfymuyymqgh5be86h/4nz7egxe/xukmvgsaupuurqcmfglhefvghoegr0br01mtk3pzm7n1xfolugkay8ijurlgbuwgognqwzwguvapgbweqwov7fuhfrxsm0xh8jm8glkcdsgtda5xrucbttmpyf8lbkuaqenum5nxfuesthemjcsvnd1yu1nwgpkgostv4mcgjmbvoguorhmopysnjvbq0bpagqpjtl+r3jqxzpviememodv/97hxf71uosndp6nrngiymehmggfnbxudjeohmqsgysgq8qspjsba5ppjsdeqonamu96ejpks5mztwowqj6dhu3nmbpell3tzdze5/7lvrbovvkbb7gz+30ww33zonxt5wlwmnafltizc1ej76bh1kici6qm7di8erefryxqxj5gh6xqzxh9gpxdm4+o4s=</latext> <latext sha1_base64="xwere7o+5eyox7juyjrcjl6wc=">aaadf3cbvnbhmxepzu+snhpk5crfiuqsq0w5oscmevhcay0gevmqgypbojla83sxfkxwvghvwf 1 ] 00(d) 00(d) u u Fgure 1: Schematc of energy transfer between the four energy components: mean knetc energy, mean nternal energy, solenodal TKE, and dlatatonal TKE. the mean knetc energy. The self-consstent feedback mechansm for the sudden vscous dsspaton reles on these complex nteractons, and thus can only be represented usng the fnte-mach-number formulaton and not the low-mach-number assumpton. 3 Computatonal detals Drect numercal smulatons of eqs. (5) to (7) are carred out wth the Mranda code developed at the awrence vermore Natonal aboratory. Ths solver employs a tenth-order accurate Pade scheme [18] for the dscretzaton of the spatal dervatves, and a fourth-order, low-storage, fvestep Runge-Kutta solver [19] for the temporal dervatves. An eghth-order compact flter s appled to the conserved varables ρ, ρu, and Et after each substep of the Runge-Kutta scheme, for the purposes of stablty. Mranda reles on the artfcal-flud-property approach to stablze shock waves and contact dscontnutes. Thus, an artfcal bulk vscosty β s ntroduced n the defnton of the vscous stress tensor, and an artfcal thermal conductvty κ s added to the thermal conductvty κ of 6
D(κ)/u 3 η 3.0.5.0 1.5 1.0 0.5 0.0 0.0 0.5 1.0 1.5.0 κη Fgure : The dsspaton spectrum, normalzed by the Kolmogorov velocty u η smulatons of lnearly-forced compressble turbulence. = (ɛν) 1/4, for the flud. The artfcal bulk vscosty and artfcal thermal conductvty are computed as β = C β ρd(d), (0) κ = C κ ρ c v D(T ). (1) T t In the above, t s the tme step, the overbar denotes a truncated-gaussan flter, and D( ) s an eghth-order dervatve operator defned as ( 8 D( ) = max x 8 x10, 8 y 8 y10, 8 ) z 8 z10. () Ths operator strongly bases the artfcal propertes towards hgh wave numbers. The coeffcents C β = 0.07 and C κ = 0.001 have been calbrated for smulatons relevant to nertal confnement fuson, see for example [0, 4, 1]. For further detals or capabltes of the code, the reader s referred to [, 3, 4]. The computatonal doman conssts of a cube of length π, wth a unform dstrbuton of 56 3 grd ponts. Perodc boundary condtons are appled on all sdes of the cube. The rato of specfc heats has a value of γ = 5/3, and the Prandtl number s set to P r = 1. The gas constant s computed as R = R u /M, where the unversal gas constant s R u = 8.314474e10 7 (cgs unts), and the molar mass used s that of Deuterum,.e. M =.01410. Statstcal quanttes are obtaned by averagng over all nodes of the mesh. The ntal flow feld s extracted from a smulaton of lnearly-forced compressble turbulence [5, 6]. Ths prelmnary smulaton s carred out for a duraton of 18 ntal eddy-turn-over tmes. The forcng coeffcents ntroduced n [5] requre the specfcaton of a pror values for the solenodal and dlatatonal dsspatons. These two quanttes were obtaned from specfyng a total dsspaton ɛ = ɛ (s) + ɛ (d) and a dsspaton rato ɛ (d) /ɛ (s). The value of the total dsspaton s chosen a pror so that the correspondng Kolmogorov scale η = (ν 3 /ɛ) 1/4 s suffcently large compared to the grd spacng. Fgure plots the dsspaton spectrum at the end of the prelmnary forced-turbulence smulaton, and llustrates the range of scales resolved on the mesh,.e. 0 κη. Ths fgure 7
shows that we capture the long tal for hgh wavenumbers as t smoothly approaches a value of zero (a model spectrum for comparson s shown n [7]), and thus do not predct a fcttous energy pleup or unphyscal rapd decay at the hghest wave modes. Ths thus serves as further evdence that the combnaton of forcng coeffcents and mesh resoluton chosen approprately capture all of the dsspatve scales, as should be the case for any drect numercal smulaton. The rato of dsspatons was set to ɛ (d) /ɛ (s) = 0.01. Smulatons of compressng turbulence that were ntalzed from a lnearly-forced case wth ɛ (d) /ɛ (s) = 1.0 were also carred out. Results for ths hgher ntal dsspaton rato are qualtatvely smlar to those of the lower dsspaton rato and are therefore not ncluded n ths paper. The turbulent Mach number and Taylor-scale Reynolds number are defned as ũ M t = u, (3) c ( ) 0k 1/ Re λ =, (4) 3ɛν where c = γrt s the speed of sound, and ν = µ / ρ s the averaged knematc vscosty. The extracted turbulent feld at the end of the lnear forcng has M t 0.65 and Re λ 70. The correspondng rato of dlatatonal to solenodal TKE s k (d) /k (s) = 0.033. For the lnearly-forced smulatons, the power law exponent s set to the tradtonal flud value of n = 3/4. However, once the sotropc compresson s appled to the ntal flow feld, the power law exponent s swtched to the value used n [1], namely n = 5/, so as to reproduce the sudden vscous dsspaton mechansm. The ntal rato of the artfcal dsspaton [6] to physcal dsspaton s 0.016, whch decays rapdly as the compresson s ntated. Thus, the smulatons are mostly affected by physcal rather than artfcal dsspatve mechansms, as should be the case for a properly-refned drect numercal smulaton [4]. 4 Results The analyss of the self-consstent feedback mechansm s dvded nto two subsectons. The frst focuses on the behavor of the TKEs and the varous mechansms depcted n fg. 1 that modulate ther temporal evoluton. The second half of the analyss s centered around the resultng evoluton of the nternal energy, and the amplfcaton of the temperature due to the vscous dsspaton. 4.1 Turbulent knetc energes 4.1.1 Profle hstores The tme evoluton of the solenodal and dlatatonal TKEs are shown n fgs. 3a and 3b, respectvely. As done n [1], rather than plottng the TKE evoluton aganst tme, a parameterzaton n terms of the length of the doman s used, and thus tme progresses from rght to left. Also equvalent to the results n [1], the TKE evolutons for dfferent compresson speeds are shown. These dfferent cases are labelled by the ntal value of the RDT parameter S = Sk/ɛ, where S = / s the nverse tme scale of the compresson, and k/ɛ s the tme scale of the turbulence. For suffcently large values of S, the compresson s rapd enough that the non-lnear turbulenceturbulence nteractons are neglgble, and the evoluton of the turbulence s descrbed exactly by rapd-dstorton theory (RDT) [7, 11]. 8
k (s) /k (s) 0 10 3 10 10 1 S 0 = 0.50 S 0 = 5.0 S 0 = 50 S 0 = 500 RDT 1/ k (d) /k (d) 0 10 3 10 10 1 10 0 10 0 10 1 0.0 0. 0.4 0.6 0.8 1.0 (a) 10 1 0.0 0. 0.4 0.6 0.8 1.0 (b) Fgure 3: Evoluton of (a) solenodal and (b) dlatatonal TKE as a functon of the sze of the doman. The solenodal and dlatatonal TKEs are normalzed by ther ntal values k (s) 0 and k (d) 0, respectvely. The ntal length of the doman s = 1, whch decreases as tme progresses. The colors correspond to dfferent values of S0, whch s the shear parameter S at tme t = 0. As fg. 3a shows, the evoluton of solenodal TKE exhbts the sudden vscous dsspaton mechansm of [1]. Even though the compresson speeds used n ths study are dfferent from those of [1], there s strong qualtatve agreement between the results shown n fg. 3a and those n fg. 1 of [1]. Fgure 3b shows that the dlatatonal TKE also exhbts the sudden vscous dsspaton mechansm. Addtonally, strong agreement wth RDT [11] s observed for both the solenodal and dlatatonal modes gven the fastest compresson rates. A few dfferences between the behavor of solenodal and dlatatonal TKE are noted. Frst, ncreasngly strong oscllatons of dlatatonal TKE are observed as S0 s decreased. Ths hghly oscllatory behavor s dscussed further n secton 4.1.. Second, the dlatatonal energy lags the solenodal energy. At the last recorded nstance n tme, the solenodal TKE has decayed by more than two orders of magntude for case S0 = 50, whereas the dlatatonal TKE has decayed by about two orders of magntude only. For the S0 = 500 case, the solenodal TKE has decayed by more than three orders of magntude and the dlatatonal TKE has decreased by less than two orders of magntude. astly, for the slowest compresson rate, the solenodal and dlatatonal TKE dverge n ther ntal behavor: whereas the dlatatonal TKE slghtly ncreases untl t suddenly dsspates, the solenodal TKE decreases from the start. An alternate representaton of the evoluton of TKEs s shown n fgs. 4a and 4b. In fg. 4a the evoluton of the solenodal TKE s parameterzed by the solenodal shear parameter Ss = Sk (s) /ɛ (s), and n fg. 4b the evoluton of the dlatatonal TKE s parameterzed by the dlatatonal shear parameter Sd = Sk(d) /ɛ (d). The dashed vertcal lnes correspond to the pont n tme at whch producton s equal to dsspaton, that s, the pont at whch P (s) = ρ ɛ (s) for fg. 4a and P (d) = ρ ɛ (d) for fg. 4b. The dashed dagonal lnes correspond the assumpton of RDT scalng, for whch ɛ (α) ( k (α)) 3 for α = s, d. That s, assumng ths relatonshp between TKE and dsspaton, 9
k (s) /k (s) 0 10 3 10 10 1 10 0 10 1 10 S 0 = 0.50 S 0 = 5.0 S 0 = 50 S 0 = 500 10 3 10 10 1 10 0 10 1 10 10 3 Ss k (d) /k (d) 0 10 3 10 10 1 10 0 10 1 10 10 3 10 10 1 10 0 10 1 10 10 3 Sd (a) (b) Fgure 4: Profles of (a) solenodal and (b) dlatatonal TKE plotted aganst the shear parameter Sα = Sk (α) /ɛ (α), where α = s for fg. (a) and α = d for fg. (b). The solenodal and dlatatonal TKEs are normalzed by ther ntal values k (s) 0 and k (d) 0, respectvely. The vertcal dashed lne n the left plot corresponds to the tme at whch P (s) = ρ ɛ (s), and the vertcal dashed lne on the rght plot corresponds to the tme at whch P (d) = ρ ɛ (d). The dashed dagonal lnes correspond to eq. (5). the TKE can be expressed as k (α) = ( k (α) ) 3/ ( k (α) ) 1/ ( ɛ (α) ) 1/ ( k (α) ) 1/ S1/ S α 1/ (5) for α = s, d. It s mportant to note that k (α) does not scale smply as Sα 1/ snce S also depends on tme. However, the dependence of S on tme s gven by the predetermned and known compresson hstory of the doman. Fgure 4a shows that for compressons S0 = 5.0, 50, and 500, the ncrease n solenodal TKE s n close agreement wth RDT. That s, the RDT assumpton of ɛ (α) ( k (α)) 3 holds qute well durng the ntal phase of the compresson. The agreement wth the scalng of eq. (5) could be benefcal for modelng purposes. Sgnfcant dvergence from the RDT scalng occurs once the vertcal lne at whch solenodal producton equals solenodal dsspaton s reached. After ths pont, the solenodal dsspaton overtakes the solenodal producton, and the turbulence decays. For the S0 = 0.50 case, the compresson s slow enough that the solenodal producton s never larger than the solenodal dsspaton, and thus the entre curve s located to the left of the vertcal dashed lne. Fgure 4b shows a smlar trend. We frst note that the rapd oscllatons n the dlatatonal profles correspondng to slow compresson speeds are also evdent n ths fgure. The agreement wth the scalng ɛ (d) ( k (d)) 3 stll holds for the S 0 = 5.0, 50, and 500 cases, although, for the S0 = 5.0 case, ths agreement s not as strong as that of the correspondng solenodal feld. More mportantly, the vertcal lne at whch dlatatonal producton equals dlatatonal dsspaton no 10
longer demarcates the domans of ncreasng and decreasng turbulence for all four cases, snce for S0 = 50 the dlatatonal TKE keeps on ncreasng after ths vertcal lne s reached. astly, fg. 4b shows that the decrease n energy s slower than that observed n fg. 4a for cases S0 = 50 and 500. Ths suggests that the dlatatonal dsspaton s actng aganst an addtonal source, whch, as wll be shown n secton 4.1., turns out to be the pressure dlataton. 4.1. Budgets Fgures 5 to 8 contan the TKE budget for the solenodal and dlatatonal felds. For the S0 = 0.5 case shown n fg. 5, oscllatons n the pressure dlataton and dlatatonal dsspaton are observed. The magntude of the oscllatons n ρ ɛ (d) are sgnfcantly smaller than those of PD. We also note that the oscllatons of the pressure dlataton and dlatatonal dsspaton are correlated, wth the dlatatonal dsspaton slghtly laggng the pressure dlataton. Moreover, the oscllatons n k (d) shown n fg. 3b are also correlated wth P D, wth k (d) laggng behnd PD. Ths serves as evdence that pressure dlataton s responsble for the oscllatory behavor of the dlatatonal TKE. The strong oscllatory nature of P D has been observed elsewhere, see for example [8, 9] for the case of forced turbulence and [9, 30] for sheared turbulence. For the S0 = 5.0 case shown n fg. 6 the oscllatons n P D and ρ ɛ (d) have been attenuated. Fgures 7 and 8 show that as the compresson speed s ncreased to S0 = 50 and 500, P D and ρ ɛ(d) do not exhbt oscllatons up to the last smulated nstance n tme. Some trends shown by fgs. 5 to 8 are mportant to be noted. As the compresson speed ncreases, the peak values of the profles occur at smaller and smaller doman lengths, as hghlghted by the dfferent ranges used for the x axs. Ths s n agreement wth fg. 3 and the results n [1]. However, more relevant to the current study s that, as the compresson speed ncreases, peak values for the dlatatonal TKE sources occur at even smaller doman lengths than those of the solenodal energy. For example, for compresson speeds S0 = 50 and 500, peak values for the solenodal dsspaton and producton are reached before those of dlatatonal dsspaton and pressure dlataton. These last two terms begn to ncrease rapdly only after the solenodal dsspaton and producton are already decayng from ther maxmum values. Ths lag n the dlatatonal modes wth respect to the solenodal ones was also observed n fg. 3 as prevously noted, but s now more apparent due to the smaller scale of the x axs n fgs. 5 to 8. The second trend to hghlght s that, as the compresson speed s ncreased, the solenodal sources do not show as dstnct a change n behavor as the dlatatonal sources. The solenodal dsspaton and producton profles retan the same shape, though wth dfferent magntudes, for the dfferent compresson rates. On the other hand, the pressure dlataton and dlatatonal dsspaton go from a hghly-oscllatory pattern to smoothly-varyng non-oscllatory profles for the fastest compresson speed. Ths poses formdable challenges for modelng purposes. The last trend to hghlght s that the pressure dlataton s ether skewed towards postve values, as s the case for S0 = 0.5, or s postve throughout the entre compresson. Ths s further exemplfed by lookng at table, whch shows the ntegrated values of the energy transfer mechansms, from the ntal to the last avalable smulated tme. All ntegrated values for the pressure dlataton are postve. Thus, P D behaves more as a source rather than a snk or a neutral term n the balance of dlatatonal TKE. Thus, the dlatatonal dsspaton needs to counteract the effect of both the dlatatonal producton and pressure dlataton for the sudden vscous dsspaton to occur n the dlatatonal feld. Gven that for the two fastest compressons the ntegrated contrbuton of P D s almost as large as that of ρ ɛ (d), t s thus not unexpected that the dlatatonal TKE decays at a slower rate than the solenodal TKE, as shown n fg. 4. 11
8 103 ρ ɛ (s) 8 103 P D 6 4 P (s) IA (s) 6 4 ρ ɛ (d) P (d) IA (d) 0 0 0.1 0. 0.3 0.4 0.5 0.6 0.1 0. 0.3 0.4 0.5 0.6 Fgure 5: TKE budget for the solenodal mode on the left and the dlatatonal mode on the rght, for S0 = 0.5. The ntal length of the doman s = 1, whch decreases as tme progresses. All terms have been normalzed by ρ 0 U0 3/ 0. 8 106 8 106 6 6 4 4 0 0.05 0.10 0.15 0.0 0.5 0 0.05 0.10 0.15 0.0 0.5 Fgure 6: TKE budget for the solenodal mode on the left and the dlatatonal mode on the rght, for S0 = 5.0. The ntal length of the doman s = 1, whch decreases as tme progresses. All terms have been normalzed by ρ 0 U0 3/ 0. The same legend as that of fg. 5 apples to the plots above. 1
8 109 8 109 6 6 4 4 0 0.0 0.04 0.06 0.08 0.10 0 0.0 0.04 0.06 0.08 0.10 Fgure 7: TKE budget for the solenodal mode on the left and the dlatatonal mode on the rght, for S0 = 50. The ntal length of the doman s = 1, whch decreases as tme progresses. All terms have been normalzed by ρ 0 U0 3/ 0. The same legend as that of fg. 5 apples to the plots above. 8 101 8 101 6 6 4 4 0 0.01 0.0 0.03 0.04 0.05 0 0.01 0.0 0.03 0.04 0.05 Fgure 8: TKE budget for the solenodal mode on the left and the dlatatonal mode on the rght, for S0 = 500. The ntal length of the doman s = 1, whch decreases as tme progresses. All terms have been normalzed by ρ 0 U0 3/ 0. The same legend as that of fg. 5 apples to the plots above. 13
Table : Integrated energy sources gven dfferent compresson speeds. All values are normalzed by ρ 0 U 0. S 0 = 0.5 S 0 = 5.0 S 0 = 50 S 0 = 500 IA (s).3e-01 4.80e+01 4.45e+03 4.90e+03 IA (d) -.3e-01-4.80e+01-4.45e+03-4.90e+03 P (s).89e+01 3.08e+03 1.97e+05 1.03e+07 P (d).80e+00 1.57e+0 1.39e+04 3.3e+05 ρ ɛ (s) 1.46e+0 1.54e+04 9.8e+05 5.07e+07 ρ ɛ (d) 1.87e+01 1.07e+03 3.4e+05 6.59e+06 P D 4.89e+00 3.34e+0.64e+05 5.15e+06 M W 3.34e+05.69e+07 1.35e+09.14e+10 E (s) (κ) 10 14 10 11 10 8 0.7 0.38 0.0 0.10 0.04 E (d) (κ) 10 14 10 11 10 8 10 5 10 5 10 10 0 10 1 10 10 3 κ (a) 10 10 0 10 1 10 10 3 κ (b) Fgure 9: Energy spectra for the (a) solenodal and (b) dlatatonal TKE, at dfferent tmes (or doman lengths) throughout the compresson. The spectra correspond to the S0 = 5.0 case. 4.1.3 Spectra The energy spectra for the solenodal and dlatatonal felds are shown n fg. 9, for the compresson speed of S0 = 5.0. The profle obtaned at 0.10 corresponds to a pont n tme for whch the sudden vscous dsspaton mechansm s takng place, and the profle at 0.04 to a tme for whch most of the turbulence has already been dsspated. The shapes and trends exhbted by the solenodal and dlatatonal spectra are equal to each other. Addtonally, these profles are n qualtatve agreement wth results shown n [1]. As the compresson proceeds, the energy n the hgher modes decreases whereas the energy n the lower modes ncreases. The set of modes for whch the energy decreases expands as the compresson progresses, and eventually even the lower modes are dsspated, as shown by the profle at 0.04. 14
T / T0 10 6 10 5 10 4 10 3 10 10 1 10 0 10 1 S 0 = 500 S 0 = 50 S 0 = 5.0 S 0 = 0.50 1/ 10 0.0 0. 0.4 0.6 0.8 1.0 Fgure 10: Evoluton of temperature as a functon of the sze of the doman. The ntal temperature s denoted by T 0. 4. Internal energy The temperature evolutons as a functon of the doman length are shown n fg. 10 for all the compresson speeds. These are also compared aganst the 1/ temperature scalng correspondng to an adabatc sentropc process wth γ = 5/3, as assumed n [1]. As the fgure shows, the temperature evolutons are n very close agreement wth the adabatc scalng. Ths ndcates that the terms n the mean nternal energy equaton neglected under the assumpton of adabatc compresson, namely the solenodal dsspaton, dlatatonal dsspaton, and pressure dlataton, do not provde a strong contrbuton towards the ncrease of temperature for the current smulatons. The neglgble effect of the dsspatons and the pressure dlataton s confrmed by comparng the source terms of the mean nternal energy, as s done n fgs. 11 to 14. These fgures show that, throughout the compresson, the domnant source for the mean nternal energy equaton s the mechancal work, whch takes the form of MW = 3 P / for the gven sotropc compresson of eq. (4). For all compresson speeds tested, the solenodal dsspaton, dlatatonal dsspaton, and pressure dlataton are eclpsed by the mechancal work at all tmes durng the compresson. However, for the two fastest compresson rates, the peak values of the dlatatonal dsspaton and pressure dlataton are not acheved by the last-avalable smulated nstance n tme. Nonetheless, as shown n fg. 3, by ths last smulated nstance n tme the dlatatonal TKE has already surpassed ts peak value and has dsspated by more than an order of magntude, and t s thus unlkely that the dlatatonal dsspaton and pressure dlataton wll ever overtake the mechancal work. For the specfc-heat rato of γ = 5/3 and the assumpton of an adabatc compresson, the mechancal work scales as 1/ 6, whch s shown as black dots n fgs. 11 to 14. Snce the mechancal work overpowers the other sources of mean nternal energy, t s expected that the assumpton of an adabatc compresson would hold as well as shown n fgs. 11 to 14. The domnance of the mechancal work can be further exemplfed by consderng the ntegrated values of the mean nternal energy sources, shown n Table. The tme-ntegrated contrbuton towards the ncrease of temperature due to mechancal work s at least three orders of magn- 15
10 1 P D ρ ɛ (d) 10 9 ρ ɛ (s) MW 10 6 RDT 1/ 6 10 3 1.0 104 0.8 0.6 0.4 0. 10 0 0.0 10 3 0. 0.4 0.6 0.8 1.0 0. 0.1 0. 0.3 0.4 0.5 0.6 Fgure 11: Mean nternal energy budget n log scale on the left and lnear scale on the rght, for S = 0.50. The ntal length of the doman s = 1, whch decreases as tme progresses. All terms have been normalzed by ρ 0 U 3 0 / 0. The absolute value of P D s shown for the log-scale plot. 10 1 10 9 10 6 10 3 1.0 107 0.8 0.6 0.4 0. 10 0 10 3 0.0 0. 0.4 0.6 0.8 1.0 0.0 0. 0.05 0.10 0.15 0.0 0.5 Fgure 1: Mean nternal energy budget n log scale on the left and lnear scale on the rght, for S = 5.0. The ntal length of the doman s = 1, whch decreases as tme progresses. All terms have been normalzed by ρ 0 U 3 0 / 0. The absolute value of P D s shown for the log-scale plot. The same legend as that of fg. 11 apples to the plots above. 16
10 1 10 9 10 6 10 3 1.0 1010 0.8 0.6 0.4 0. 10 0 10 3 0.0 0. 0.4 0.6 0.8 1.0 0.0 0. 0.0 0.04 0.06 0.08 0.10 Fgure 13: Mean nternal energy budget n log scale on the left and lnear scale on the rght, for S = 50. The ntal length of the doman s = 1, whch decreases as tme progresses. All terms have been normalzed by ρ 0 U 3 0 / 0. The absolute value of P D s shown for the log-scale plot. The same legend as that of fg. 11 apples to the plots above. 10 1 10 9 10 6 10 3 1.0 1013 0.8 0.6 0.4 0. 10 0 0.0 10 3 0.0 0. 0.4 0.6 0.8 1.0 0. 0.01 0.0 0.03 0.04 0.05 Fgure 14: Mean nternal energy budget n log scale on the left and lnear scale on the rght, for S = 500. The ntal length of the doman s = 1, whch decreases as tme progresses. All terms have been normalzed by ρ 0 U 3 0 / 0. The absolute value of P D s shown for the log-scale plot. The same legend as that of fg. 11 apples to the plots above. 17
tude larger than the second most sgnfcant tme-ntegrated source term, namely, the solenodal dsspaton. An alternate metrc for hghlghtng the domnance of the mechancal work s the comparson of the tme-ntegrated contrbuton from the TKEs to mean nternal energy aganst the tme-ntegrated contrbuton from the mean knetc energy to the mean nternal energy. The rato of these two factors for the four cases S0 = 0.50, 5.0, 50, and, 500 s 0.0005, 0.0006, 0.0008, and 0.00, respectvely. Gven that, for the parameters used n these smulatons, the dsspated turbulent knetc energy does not sgnfcantly ncrease the temperature of the system above the adabatc predcton, t s crucal to determne under whch condtons would the dsspated TKE actually lead to meanngful ncreases n temperature. To do ths, we make use of the relaton ) d (Ũ + k = 0, (6) dt Ũ (a) whch was derved n appendx C. Ũ (a) s the mean nternal energy of the system gven the dealzed adabatc compresson, and s thus gven by Ũ (a) = Ũ0 where Ũ0 s the ntal value of Ũ. Integratng from the ntal tme t 0 to a fnal tme t f, one obtans Ũ + k Ũ (a) tf = Ũ + k Ũ (a) = 1 + k 0 = 1 + 5 t0 Ũ 0 9 M u,0. (7) In the above we have made use of the defnton of the fluctuatng Mach number [9] ũ M u = u c( T ), (8) whose ntal value s denoted by M u,0. We ntroduce T (a) = T 0 as the temperature correspondng to an adabatc compresson. If we defne t f as the tme by whch all of the turbulent knetc energy has dsspated, and T f and then eq. (7) can be expressed as T (a) f as the temperatures T and T (a) at tmes t > t f, respectvely, T f T (a) f = 1 + 5 9 M u,0. (9) The above relaton hghlghts a few notable aspects of the compresson mechansm. Gven T 0 and k 0, M u,0 s known, whch, along wth values of smaller than those correspondng to the tme t f, can be used n eq. (9) to obtan the temperature acheved after the TKE has been fully dsspated. It s therefore unnecessary to carry out expensve numercal smulatons to predct the temperature of the system post TKE depleton. The second aspect to hghlght s that the temperature rato (a) T f / T f s ndependent of the compresson speed. However, the relatonshp n eq. (9) holds for tmes t > t f, where t f s dfferent for the varous compresson speeds. For the smulatons descrbed n ths paper, the ntal fluctuatng Mach number mmedately precedng the start of the compresson s M u,0 = 0.651. Usng eq. (9), ths gves T (a) f / T f = 1.35. Table 3 lsts ths rato computed from smulaton data avalable at the last smulated nstance n tme, for the four compresson speeds. As the table shows, there s strong agreement wth the analytcal value of 1.35. The slghtly lower rato for the fastest compresson s most lkely due to 18
Table 3: Rato T (a) f / T f obtaned at the last smulated nstance n tme. S 0 = 0.50 S 0 = 5.0 S 0 = 50 S 0 = 500 T f T (a) f 1.3 1.3 1.3 1.30 the fact that all of the TKE, specally the dlatatonal TKE, has not yet fully dsspated nto heat. Equaton (9) can now thus be used to predct under whch condtons the dsspated TKE would lead to meanngful ncreases n temperature. For subsonc ntal fluctuatng Mach numbers, the temperature post TKE depleton can be up to about 1.5 tmes larger than that obtaned wth an adabatc compresson. If supersonc Mach numbers are used, such as M u,0 = and 5, then the temperature post TKE depleton would be about 3 and 15 tmes larger, respectvely, than for an adabatc compresson. For hghly supersonc turbulence such as that encountered n the nterstellar medum [31, 3], a Mach number of M u,0 = 17 would lead to fnal temperatures about 160 tmes hgher than those predcted assumng an adabatc scalng. As stated n [16], the hot spot of an nertal-confnement-fuson capsule can be characterzed by a turbulent Mach number M t 0.4. Usng ths value n eq. (9) leads to T (a) f / T f 1.09. Ths ncrease of temperature s mnmal, and s eclpsed by the effect of the mechancal work. For example, f we assume that the sudden vscous dsspaton of TKE occurs at = 0.1, a small reducton of the doman sze to = 0.0958 would already allow the mechancal work to generate an equvalent ncrease n temperature. It s thus expected that only for flow felds wth large ntal Mach numbers would the self-consstent feedback mechansm lead to sudden dsspatons wth sgnfcant effects. 5 Concludng remarks A sudden vscous dsspaton of plasma turbulence under compresson was demonstrated n [1]. We expand on ths prevous work by accountng for the self-consstent feedback loop assocated wth ths vscous mechansm. The feedback loop entals a transfer of energy from the turbulence towards the nternal energy, and the subsequent ncreased temperatures and vscostes that n turn accelerate the orgnal dsspaton of TKE. Although prevous efforts have reproduced the sudden dsspaton of TKE, these do not capture the subsequent effect of the dsspated energy on the temperature, and the consequences thereof. Ths lmtaton s due to the use of the zero-machlmt assumpton. To capture the ncrease of nternal energy resultng from the dsspated TKE, and thus account for the entre self-consstent feedback loop, drect numercal smulatons have been carred out usng a fnte-mach number formulaton that solves transport equatons for the densty, fluctuatng velocty, and total energy. The analyss of the self-consstent feedback loop was dvded nto two steps: the frst focused on the evoluton of the solenodal and dlatatonal TKEs, and the second on the evoluton of the mean nternal energy as t absorbs the dsspated TKE. Results show that both the solenodal and dlatatonal TKE exhbt the sudden vscous dsspaton mechansm, wth the dsspaton of dlatatonal TKE slghtly laggng that of solenodal TKE. Moreover, large oscllatons n the temporal evoluton of dlatatonal TKE for slow compresson rates are observed, whch are correlated wth the hghly-oscllatory nature of the pressure dlataton. The pressure dlataton consttutes a two-way energy transfer between the dlatatonal TKE and the mean nternal energy of the system. However, a detaled analyss of the dlatatonal TKE 19
budget shows that the tme-ntegrated effect of the pressure dlataton s to transfer energy from heat towards dlatatonal TKE, even for cases when the pressure dlataton transfers energy n both drectons on short tme scales. Thus, the dlatatonal dsspaton needs to overcome both the dlatatonal producton and pressure dlataton for the sudden vscous dsspaton of dlatatonal TKE to take place, whch explans the delayed dsspaton of dlatatonal TKE wth respect to solenodal TKE. An analyss of the sources for the mean nternal energy shows that mechancal work, whch transforms energy from the mean flow to ncrease heat, domnates all other sources of mean nternal energy for the turbulent Mach numbers chosen n ths study. For all nstances n tme, the mechancal work term s larger, often by multple orders of magntude, than the solenodal and dlatatonal dsspaton and the pressure dlataton. As a result, the contrbuton of the dsspated TKE towards the ncrease of temperature s mnmal, and the temperature evoluton closely follows an adabatc scalng. Ths valdates prevous efforts [1, 13, 14, 15] that reled on a fxed adabatc scalng for the temperature evoluton. An analytcal expresson was also derved for the rato of the temperature post TKE depleton to the temperature obtaned at the same nstance of the compresson assumng an adabatc scalng. Ths rato depends on the ntal fluctuatng Mach number only, ndcatng that the rate of compresson does not affect the magntudes of the temperature post TKE depleton. The derved analytcal expresson also shows that for subsonc ntal fluctuatng Mach numbers, the true temperature of the system s not substantally larger than that obtaned f the TKE contrbuton s neglected. To provde a pont of reference, t was shown that for an adabatc compresson where only the mechancal work s actve, reducng the doman from = 0.1 to = 0.0958 would have an equvalent effect on temperature as that of the suddenly dsspated TKE. Ths thus ndcates that the potental of the sudden vscous dsspaton mechansm to sgnfcantly enhance the heatng of the plasma by dsspatng the nherent turbulence could be lmted. Nonetheless, ths mechansm could serve as an effectve tool to prevent the occurrence of turbulent fluctuatons, whch are responsble for detrmental mxng n nertal-confnement-fuson capsules. It s crucal to hghlght that the fnte-mach-number framework chosen here, although more general than the zero-mach-number formalsm, s stll mssng physcs relevant to nertal-confnement fuson such as non-deal equatons of state, radaton transport, multple speces, plasma vscosty models, separate on and electron temperatures, and alpha heatng. Thus, these factors need to be explored to provde a defnte assessment on the ablty of the sudden vscous dsspaton mechansm to mprove the performance of nertal confnement fuson. 6 Acknowledgments The authors wsh to thank Dr. O. Schllng for helpful comments on drafts of the manuscrpt. Ths work was performed under the auspces of the U.S. Department of Energy by awrence vermore Natonal aboratory under Contract DE-AC5-07NA7344. A Dervaton of the fnte-mach-number Naver-Stokes equatons for sotropc mean compresson The dervaton of the governng equatons used for the computatonal smulatons n ths study s detaled below. Ths dervaton s dvded nto fve dstnct steps, each descrbed n the fve 0