arxiv: v1 [physics.flu-dyn] 31 May 2018

Size: px
Start display at page:

Download "arxiv: v1 [physics.flu-dyn] 31 May 2018"

Transcription

1 Self-consstent feedback mechansm for the sudden vscous dsspaton of fnte-mach-number compressng turbulence arxv: v1 [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

2 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

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

4 ρ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

5 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

6 <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

7 D(κ)/u 3 η κη 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 κ = 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 = e10 7 (cgs unts), and the molar mass used s that of Deuterum,.e. M = 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

8 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) = 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) = 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 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

9 k (s) /k (s) S 0 = 0.50 S 0 = 5.0 S 0 = 50 S 0 = 500 RDT 1/ k (d) /k (d) (a) (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 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

10 k (s) /k (s) S 0 = 0.50 S 0 = 5.0 S 0 = 50 S 0 = Ss k (d) /k (d) 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

11 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 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

12 8 103 ρ ɛ (s) P D 6 4 P (s) IA (s) 6 4 ρ ɛ (d) P (d) IA (d) 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/ 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

13 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 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

14 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 e e e+03 IA (d) -.3e e e e+03 P (s).89e e e e+07 P (d).80e e e e+05 ρ ɛ (s) 1.46e e e e+07 ρ ɛ (d) 1.87e e e e+06 P D 4.89e e+0.64e e+06 M W 3.34e+05.69e e+09.14e+10 E (s) (κ) E (d) (κ) κ (a) κ (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 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

15 T / T S 0 = 500 S 0 = 50 S 0 = 5.0 S 0 = / Fgure 10: Evoluton of temperature as a functon of the sze of the doman. The ntal temperature s denoted by T 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

16 10 1 P D ρ ɛ (d) 10 9 ρ ɛ (s) MW 10 6 RDT 1/ Fgure 11: Mean nternal energy budget n log scale on the left and lnear scale on the rght, for S = 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 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

17 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 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

18 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 , , , 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 = 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 = 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 = Usng eq. (9), ths gves T (a) f / T f = 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 The slghtly lower rato for the fastest compresson s most lkely due to 18

19 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 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 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 = 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

20 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 = 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

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:

More information

Turbulent Flow. Turbulent Flow

Turbulent Flow. Turbulent Flow http://www.youtube.com/watch?v=xoll2kedog&feature=related http://br.youtube.com/watch?v=7kkftgx2any http://br.youtube.com/watch?v=vqhxihpvcvu 1. Caothc fluctuatons wth a wde range of frequences and

More information

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look

More information

Handout: Large Eddy Simulation I. Introduction to Subgrid-Scale (SGS) Models

Handout: Large Eddy Simulation I. Introduction to Subgrid-Scale (SGS) Models Handout: Large Eddy mulaton I 058:68 Turbulent flows G. Constantnescu Introducton to ubgrd-cale (G) Models G tresses should depend on: Local large-scale feld or Past hstory of local flud (va PDE) Not all

More information

Flow equations To simulate the flow, the Navier-Stokes system that includes continuity and momentum equations is solved

Flow equations To simulate the flow, the Navier-Stokes system that includes continuity and momentum equations is solved Smulaton of nose generaton and propagaton caused by the turbulent flow around bluff bodes Zamotn Krll e-mal: krart@gmal.com, cq: 958886 Summary Accurate predctons of nose generaton and spread n turbulent

More information

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle

More information

Lecture 5.8 Flux Vector Splitting

Lecture 5.8 Flux Vector Splitting Lecture 5.8 Flux Vector Splttng 1 Flux Vector Splttng The vector E n (5.7.) can be rewrtten as E = AU (5.8.1) (wth A as gven n (5.7.4) or (5.7.6) ) whenever, the equaton of state s of the separable form

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc

More information

Lecture 12. Modeling of Turbulent Combustion

Lecture 12. Modeling of Turbulent Combustion Lecture 12. Modelng of Turbulent Combuston X.S. Ba Modelng of TC Content drect numercal smulaton (DNS) Statstcal approach (RANS) Modelng of turbulent non-premxed flames Modelng of turbulent premxed flames

More information

STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS

STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS Blucher Mechancal Engneerng Proceedngs May 0, vol., num. www.proceedngs.blucher.com.br/evento/0wccm STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS Takahko Kurahash,

More information

Turbulence and its Modelling

Turbulence and its Modelling School of Mechancal Aerospace and Cvl Engneerng 3rd Year Flud Mechancs Introducton In earler lectures we have consdered how flow nstabltes develop, and noted that above some crtcal Reynolds number flows

More information

Publication 2006/01. Transport Equations in Incompressible. Lars Davidson

Publication 2006/01. Transport Equations in Incompressible. Lars Davidson Publcaton 2006/01 Transport Equatons n Incompressble URANS and LES Lars Davdson Dvson of Flud Dynamcs Department of Appled Mechancs Chalmers Unversty of Technology Göteborg, Sweden, May 2006 Transport

More information

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate Internatonal Journal of Mathematcs and Systems Scence (018) Volume 1 do:10.494/jmss.v1.815 (Onlne Frst)A Lattce Boltzmann Scheme for Dffuson Equaton n Sphercal Coordnate Debabrata Datta 1 *, T K Pal 1

More information

One-sided finite-difference approximations suitable for use with Richardson extrapolation

One-sided finite-difference approximations suitable for use with Richardson extrapolation Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,

More information

A Numerical Study of Heat Transfer and Fluid Flow past Single Tube

A Numerical Study of Heat Transfer and Fluid Flow past Single Tube A Numercal Study of Heat ransfer and Flud Flow past Sngle ube ZEINAB SAYED ABDEL-REHIM Mechancal Engneerng Natonal Research Center El-Bohos Street, Dokk, Gza EGYP abdelrehmz@yahoo.com Abstract: - A numercal

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

Tools for large-eddy simulation

Tools for large-eddy simulation Center for Turbulence Research Proceedngs of the Summer Program 00 117 Tools for large-eddy smulaton By Davd A. Caughey AND Grdhar Jothprasad A computer code has been developed for solvng the ncompressble

More information

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Structure and Drive Paul A. Jensen Copyright July 20, 2003 Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.

More information

Investigation of a New Monte Carlo Method for the Transitional Gas Flow

Investigation of a New Monte Carlo Method for the Transitional Gas Flow Investgaton of a New Monte Carlo Method for the Transtonal Gas Flow X. Luo and Chr. Day Karlsruhe Insttute of Technology(KIT) Insttute for Techncal Physcs 7602 Karlsruhe Germany Abstract. The Drect Smulaton

More information

Physics 181. Particle Systems

Physics 181. Particle Systems Physcs 181 Partcle Systems Overvew In these notes we dscuss the varables approprate to the descrpton of systems of partcles, ther defntons, ther relatons, and ther conservatons laws. We consder a system

More information

High resolution entropy stable scheme for shallow water equations

High resolution entropy stable scheme for shallow water equations Internatonal Symposum on Computers & Informatcs (ISCI 05) Hgh resoluton entropy stable scheme for shallow water equatons Xaohan Cheng,a, Yufeng Ne,b, Department of Appled Mathematcs, Northwestern Polytechncal

More information

FTCS Solution to the Heat Equation

FTCS Solution to the Heat Equation FTCS Soluton to the Heat Equaton ME 448/548 Notes Gerald Recktenwald Portland State Unversty Department of Mechancal Engneerng gerry@pdx.edu ME 448/548: FTCS Soluton to the Heat Equaton Overvew 1. Use

More information

( ) = ( ) + ( 0) ) ( )

( ) = ( ) + ( 0) ) ( ) EETOMAGNETI OMPATIBIITY HANDBOOK 1 hapter 9: Transent Behavor n the Tme Doman 9.1 Desgn a crcut usng reasonable values for the components that s capable of provdng a tme delay of 100 ms to a dgtal sgnal.

More information

Turbulence. Lecture 21. Non-linear Dynamics. 30 s & 40 s Taylor s work on homogeneous turbulence Kolmogorov.

Turbulence. Lecture 21. Non-linear Dynamics. 30 s & 40 s Taylor s work on homogeneous turbulence Kolmogorov. Turbulence Lecture 1 Non-lnear Dynamcs Strong non-lnearty s a key feature of turbulence. 1. Unstable, chaotc behavor.. Strongly vortcal (vortex stretchng) 3 s & 4 s Taylor s work on homogeneous turbulence

More information

Lecture Note 3. Eshelby s Inclusion II

Lecture Note 3. Eshelby s Inclusion II ME340B Elastcty of Mcroscopc Structures Stanford Unversty Wnter 004 Lecture Note 3. Eshelby s Incluson II Chrs Wenberger and We Ca c All rghts reserved January 6, 004 Contents 1 Incluson energy n an nfnte

More information

Inductance Calculation for Conductors of Arbitrary Shape

Inductance Calculation for Conductors of Arbitrary Shape CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors

More information

Appendix B. The Finite Difference Scheme

Appendix B. The Finite Difference Scheme 140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton

More information

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the

More information

Second Order Analysis

Second Order Analysis Second Order Analyss In the prevous classes we looked at a method that determnes the load correspondng to a state of bfurcaton equlbrum of a perfect frame by egenvalye analyss The system was assumed to

More information

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction ECONOMICS 5* -- NOTE (Summary) ECON 5* -- NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also

More information

Time-Varying Systems and Computations Lecture 6

Time-Varying Systems and Computations Lecture 6 Tme-Varyng Systems and Computatons Lecture 6 Klaus Depold 14. Januar 2014 The Kalman Flter The Kalman estmaton flter attempts to estmate the actual state of an unknown dscrete dynamcal system, gven nosy

More information

A PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.

A PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY. Proceedngs of the th Brazlan Congress of Thermal Scences and Engneerng -- ENCIT 006 Braz. Soc. of Mechancal Scences and Engneerng -- ABCM, Curtba, Brazl,- Dec. 5-8, 006 A PROCEDURE FOR SIMULATING THE NONLINEAR

More information

Turbulence classification of load data by the frequency and severity of wind gusts. Oscar Moñux, DEWI GmbH Kevin Bleibler, DEWI GmbH

Turbulence classification of load data by the frequency and severity of wind gusts. Oscar Moñux, DEWI GmbH Kevin Bleibler, DEWI GmbH Turbulence classfcaton of load data by the frequency and severty of wnd gusts Introducton Oscar Moñux, DEWI GmbH Kevn Blebler, DEWI GmbH Durng the wnd turbne developng process, one of the most mportant

More information

Week 9 Chapter 10 Section 1-5

Week 9 Chapter 10 Section 1-5 Week 9 Chapter 10 Secton 1-5 Rotaton Rgd Object A rgd object s one that s nondeformable The relatve locatons of all partcles makng up the object reman constant All real objects are deformable to some extent,

More information

Composite Hypotheses testing

Composite Hypotheses testing Composte ypotheses testng In many hypothess testng problems there are many possble dstrbutons that can occur under each of the hypotheses. The output of the source s a set of parameters (ponts n a parameter

More information

The Finite Element Method

The Finite Element Method The Fnte Element Method GENERAL INTRODUCTION Read: Chapters 1 and 2 CONTENTS Engneerng and analyss Smulaton of a physcal process Examples mathematcal model development Approxmate solutons and methods of

More information

THE COUPLED LES - SUBGRID STOCHASTIC ACCELERATION MODEL (LES-SSAM) OF A HIGH REYNOLDS NUMBER FLOWS

THE COUPLED LES - SUBGRID STOCHASTIC ACCELERATION MODEL (LES-SSAM) OF A HIGH REYNOLDS NUMBER FLOWS /2 THE COUPLED LES - SUBGRID STOCHASTIC ACCELERATION MODEL LES-SSAM OF A HIGH REYNOLDS NUMBER FLOWS Vladmr Sabel nov DEFA/EFCA ONERA, France In collaboraton wth: Anna Chtab CORIA, Unversté de Rouen, France

More information

9 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations

9 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations Physcs 171/271 - Chapter 9R -Davd Klenfeld - Fall 2005 9 Dervaton of Rate Equatons from Sngle-Cell Conductance (Hodgkn-Huxley-lke) Equatons We consder a network of many neurons, each of whch obeys a set

More information

2 Finite difference basics

2 Finite difference basics Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T

More information

Linear Approximation with Regularization and Moving Least Squares

Linear Approximation with Regularization and Moving Least Squares Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...

More information

The Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites

The Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites 7 Asa-Pacfc Engneerng Technology Conference (APETC 7) ISBN: 978--6595-443- The Two-scale Fnte Element Errors Analyss for One Class of Thermoelastc Problem n Perodc Compostes Xaoun Deng Mngxang Deng ABSTRACT

More information

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng

More information

2016 Wiley. Study Session 2: Ethical and Professional Standards Application

2016 Wiley. Study Session 2: Ethical and Professional Standards Application 6 Wley Study Sesson : Ethcal and Professonal Standards Applcaton LESSON : CORRECTION ANALYSIS Readng 9: Correlaton and Regresson LOS 9a: Calculate and nterpret a sample covarance and a sample correlaton

More information

Uncertainty in measurements of power and energy on power networks

Uncertainty in measurements of power and energy on power networks Uncertanty n measurements of power and energy on power networks E. Manov, N. Kolev Department of Measurement and Instrumentaton, Techncal Unversty Sofa, bul. Klment Ohrdsk No8, bl., 000 Sofa, Bulgara Tel./fax:

More information

Negative Binomial Regression

Negative Binomial Regression STATGRAPHICS Rev. 9/16/2013 Negatve Bnomal Regresson Summary... 1 Data Input... 3 Statstcal Model... 3 Analyss Summary... 4 Analyss Optons... 7 Plot of Ftted Model... 8 Observed Versus Predcted... 10 Predctons...

More information

IC Engine Flow Simulation using KIVA code and A Modified Reynolds Stress Turbulence Model

IC Engine Flow Simulation using KIVA code and A Modified Reynolds Stress Turbulence Model IC Engne Flow Smulaton usng KIVA code and A Modfed Reynolds Stress Turbulence Model Satpreet Nanda and S.L. Yang Mechancal Engneerng-Engneerng Mechancs Department Mchgan Technologcal Unversty Houghton,

More information

Mathematical Preparations

Mathematical Preparations 1 Introducton Mathematcal Preparatons The theory of relatvty was developed to explan experments whch studed the propagaton of electromagnetc radaton n movng coordnate systems. Wthn expermental error the

More information

More metrics on cartesian products

More metrics on cartesian products More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of

More information

Kernel Methods and SVMs Extension

Kernel Methods and SVMs Extension Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general

More information

1 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations

1 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations Physcs 171/271 -Davd Klenfeld - Fall 2005 (revsed Wnter 2011) 1 Dervaton of Rate Equatons from Sngle-Cell Conductance (Hodgkn-Huxley-lke) Equatons We consder a network of many neurons, each of whch obeys

More information

Chapter 13: Multiple Regression

Chapter 13: Multiple Regression Chapter 13: Multple Regresson 13.1 Developng the multple-regresson Model The general model can be descrbed as: It smplfes for two ndependent varables: The sample ft parameter b 0, b 1, and b are used to

More information

829. An adaptive method for inertia force identification in cantilever under moving mass

829. An adaptive method for inertia force identification in cantilever under moving mass 89. An adaptve method for nerta force dentfcaton n cantlever under movng mass Qang Chen 1, Mnzhuo Wang, Hao Yan 3, Haonan Ye 4, Guola Yang 5 1,, 3, 4 Department of Control and System Engneerng, Nanng Unversty,

More information

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM Ganj, Z. Z., et al.: Determnaton of Temperature Dstrbuton for S111 DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM by Davood Domr GANJI

More information

In this section is given an overview of the common elasticity models.

In this section is given an overview of the common elasticity models. Secton 4.1 4.1 Elastc Solds In ths secton s gven an overvew of the common elastcty models. 4.1.1 The Lnear Elastc Sold The classcal Lnear Elastc model, or Hooean model, has the followng lnear relatonshp

More information

Lecture Notes on Linear Regression

Lecture Notes on Linear Regression Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume

More information

Note 10. Modeling and Simulation of Dynamic Systems

Note 10. Modeling and Simulation of Dynamic Systems Lecture Notes of ME 475: Introducton to Mechatroncs Note 0 Modelng and Smulaton of Dynamc Systems Department of Mechancal Engneerng, Unversty Of Saskatchewan, 57 Campus Drve, Saskatoon, SK S7N 5A9, Canada

More information

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential Open Systems: Chemcal Potental and Partal Molar Quanttes Chemcal Potental For closed systems, we have derved the followng relatonshps: du = TdS pdv dh = TdS + Vdp da = SdT pdv dg = VdP SdT For open systems,

More information

CHARACTERIZATION OF INERTIAL AND PRESSURE EFFECTS IN HOMOGENEOUS TURBULENCE

CHARACTERIZATION OF INERTIAL AND PRESSURE EFFECTS IN HOMOGENEOUS TURBULENCE CHARACTERIZATION OF INERTIAL AND PRESSURE EFFECTS IN HOMOGENEOUS TURBULENCE A Thess by RAVI KIRAN BIKKANI Submtted to the Offce of Graduate Studes of Texas A&M Unversty n partal fulfllment of the requrements

More information

A Hybrid Variational Iteration Method for Blasius Equation

A Hybrid Variational Iteration Method for Blasius Equation Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method

More information

This column is a continuation of our previous column

This column is a continuation of our previous column Comparson of Goodness of Ft Statstcs for Lnear Regresson, Part II The authors contnue ther dscusson of the correlaton coeffcent n developng a calbraton for quanttatve analyss. Jerome Workman Jr. and Howard

More information

See Book Chapter 11 2 nd Edition (Chapter 10 1 st Edition)

See Book Chapter 11 2 nd Edition (Chapter 10 1 st Edition) Count Data Models See Book Chapter 11 2 nd Edton (Chapter 10 1 st Edton) Count data consst of non-negatve nteger values Examples: number of drver route changes per week, the number of trp departure changes

More information

NON LINEAR ANALYSIS OF STRUCTURES ACCORDING TO NEW EUROPEAN DESIGN CODE

NON LINEAR ANALYSIS OF STRUCTURES ACCORDING TO NEW EUROPEAN DESIGN CODE October 1-17, 008, Bejng, Chna NON LINEAR ANALYSIS OF SRUCURES ACCORDING O NEW EUROPEAN DESIGN CODE D. Mestrovc 1, D. Czmar and M. Pende 3 1 Professor, Dept. of Structural Engneerng, Faculty of Cvl Engneerng,

More information

A large scale tsunami run-up simulation and numerical evaluation of fluid force during tsunami by using a particle method

A large scale tsunami run-up simulation and numerical evaluation of fluid force during tsunami by using a particle method A large scale tsunam run-up smulaton and numercal evaluaton of flud force durng tsunam by usng a partcle method *Mtsuteru Asa 1), Shoch Tanabe 2) and Masaharu Isshk 3) 1), 2) Department of Cvl Engneerng,

More information

Tensor Smooth Length for SPH Modelling of High Speed Impact

Tensor Smooth Length for SPH Modelling of High Speed Impact Tensor Smooth Length for SPH Modellng of Hgh Speed Impact Roman Cherepanov and Alexander Gerasmov Insttute of Appled mathematcs and mechancs, Tomsk State Unversty 634050, Lenna av. 36, Tomsk, Russa RCherepanov82@gmal.com,Ger@npmm.tsu.ru

More information

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,

More information

Temperature. Chapter Heat Engine

Temperature. Chapter Heat Engine Chapter 3 Temperature In prevous chapters of these notes we ntroduced the Prncple of Maxmum ntropy as a technque for estmatng probablty dstrbutons consstent wth constrants. In Chapter 9 we dscussed the

More information

12. The Hamilton-Jacobi Equation Michael Fowler

12. The Hamilton-Jacobi Equation Michael Fowler 1. The Hamlton-Jacob Equaton Mchael Fowler Back to Confguraton Space We ve establshed that the acton, regarded as a functon of ts coordnate endponts and tme, satsfes ( ) ( ) S q, t / t+ H qpt,, = 0, and

More information

Resource Allocation with a Budget Constraint for Computing Independent Tasks in the Cloud

Resource Allocation with a Budget Constraint for Computing Independent Tasks in the Cloud Resource Allocaton wth a Budget Constrant for Computng Independent Tasks n the Cloud Wemng Sh and Bo Hong School of Electrcal and Computer Engneerng Georga Insttute of Technology, USA 2nd IEEE Internatonal

More information

Physics 53. Rotational Motion 3. Sir, I have found you an argument, but I am not obliged to find you an understanding.

Physics 53. Rotational Motion 3. Sir, I have found you an argument, but I am not obliged to find you an understanding. Physcs 53 Rotatonal Moton 3 Sr, I have found you an argument, but I am not oblged to fnd you an understandng. Samuel Johnson Angular momentum Wth respect to rotatonal moton of a body, moment of nerta plays

More information

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1 P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the

More information

Introduction to Turbulence Modeling

Introduction to Turbulence Modeling Introducton to Turbulence Modelng Professor Ismal B. Celk West Vrgna nversty Ismal.Celk@mal.wvu.edu CFD Lab. - West Vrgna nversty I-1 Introducton to Turbulence CFD Lab. - West Vrgna nversty I-2 Introducton

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

4DVAR, according to the name, is a four-dimensional variational method.

4DVAR, according to the name, is a four-dimensional variational method. 4D-Varatonal Data Assmlaton (4D-Var) 4DVAR, accordng to the name, s a four-dmensonal varatonal method. 4D-Var s actually a drect generalzaton of 3D-Var to handle observatons that are dstrbuted n tme. The

More information

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests Smulated of the Cramér-von Mses Goodness-of-Ft Tests Steele, M., Chaselng, J. and 3 Hurst, C. School of Mathematcal and Physcal Scences, James Cook Unversty, Australan School of Envronmental Studes, Grffth

More information

Implicit Integration Henyey Method

Implicit Integration Henyey Method Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure

More information

Amplification and Relaxation of Electron Spin Polarization in Semiconductor Devices

Amplification and Relaxation of Electron Spin Polarization in Semiconductor Devices Amplfcaton and Relaxaton of Electron Spn Polarzaton n Semconductor Devces Yury V. Pershn and Vladmr Prvman Center for Quantum Devce Technology, Clarkson Unversty, Potsdam, New York 13699-570, USA Spn Relaxaton

More information

Uncertainty as the Overlap of Alternate Conditional Distributions

Uncertainty as the Overlap of Alternate Conditional Distributions Uncertanty as the Overlap of Alternate Condtonal Dstrbutons Olena Babak and Clayton V. Deutsch Centre for Computatonal Geostatstcs Department of Cvl & Envronmental Engneerng Unversty of Alberta An mportant

More information

Optimal Control of Temperature in Fluid Flow

Optimal Control of Temperature in Fluid Flow Kawahara Lab. 5 March. 27 Optmal Control of Temperature n Flud Flow Dasuke YAMAZAKI Department of Cvl Engneerng, Chuo Unversty Kasuga -3-27, Bunkyou-ku, Tokyo 2-855, Japan E-mal : d33422@educ.kc.chuo-u.ac.jp

More information

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 30 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 2 Remedes for multcollnearty Varous technques have

More information

Global Sensitivity. Tuesday 20 th February, 2018

Global Sensitivity. Tuesday 20 th February, 2018 Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values

More information

Computational Fluid Dynamics. Smoothed Particle Hydrodynamics. Simulations. Smoothing Kernels and Basis of SPH

Computational Fluid Dynamics. Smoothed Particle Hydrodynamics. Simulations. Smoothing Kernels and Basis of SPH Computatonal Flud Dynamcs If you want to learn a bt more of the math behnd flud dynamcs, read my prevous post about the Naver- Stokes equatons and Newtonan fluds. The equatons derved n the post are the

More information

TURBULENT FLOW A BEGINNER S APPROACH. Tony Saad March

TURBULENT FLOW A BEGINNER S APPROACH. Tony Saad March TURBULENT FLOW A BEGINNER S APPROACH Tony Saad March 2004 http://tsaad.uts.edu - tsaad@uts.edu CONTENTS Introducton Random processes The energy cascade mechansm The Kolmogorov hypotheses The closure problem

More information

8 Derivation of Network Rate Equations from Single- Cell Conductance Equations

8 Derivation of Network Rate Equations from Single- Cell Conductance Equations Physcs 178/278 - Davd Klenfeld - Wnter 2019 8 Dervaton of Network Rate Equatons from Sngle- Cell Conductance Equatons Our goal to derve the form of the abstract quanttes n rate equatons, such as synaptc

More information

The Order Relation and Trace Inequalities for. Hermitian Operators

The Order Relation and Trace Inequalities for. Hermitian Operators Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence

More information

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

The equation of motion of a dynamical system is given by a set of differential equations. That is (1) Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence

More information

1. Governing Equations

1. Governing Equations 1. Governng Equatons 1a. Governng Equatons for Mean Varables The governng equatons descrbe the varaton n space and tme of the zonal, merdonal and vertcal wnd components, densty, temperature, specfc humdty

More information

Supplementary Notes for Chapter 9 Mixture Thermodynamics

Supplementary Notes for Chapter 9 Mixture Thermodynamics Supplementary Notes for Chapter 9 Mxture Thermodynamcs Key ponts Nne major topcs of Chapter 9 are revewed below: 1. Notaton and operatonal equatons for mxtures 2. PVTN EOSs for mxtures 3. General effects

More information

CFD VALIDATION OF STRATIFIED TWO-PHASE FLOWS IN A HORIZONTAL CHANNEL

CFD VALIDATION OF STRATIFIED TWO-PHASE FLOWS IN A HORIZONTAL CHANNEL CFD VALIDATION OF STRATIFIED TWO-PHASE FLOWS IN A HORIZONTAL CHANNEL 1. Introducton Chrstophe Vallée and Thomas Höhne In dfferent scenaros of small break Loss of Coolant Accdent (SB-LOCA), stratfed twophase

More information

Thermodynamics General

Thermodynamics General Thermodynamcs General Lecture 1 Lecture 1 s devoted to establshng buldng blocks for dscussng thermodynamcs. In addton, the equaton of state wll be establshed. I. Buldng blocks for thermodynamcs A. Dmensons,

More information

Finite Element Modelling of truss/cable structures

Finite Element Modelling of truss/cable structures Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures

More information

Pulse Coded Modulation

Pulse Coded Modulation Pulse Coded Modulaton PCM (Pulse Coded Modulaton) s a voce codng technque defned by the ITU-T G.711 standard and t s used n dgtal telephony to encode the voce sgnal. The frst step n the analog to dgtal

More information

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity 1 Module 1 : The equaton of contnuty Lecture 1: Equaton of Contnuty 2 Advanced Heat and Mass Transfer: Modules 1. THE EQUATION OF CONTINUITY : Lectures 1-6 () () () (v) (v) Overall Mass Balance Momentum

More information

Electrical double layer: revisit based on boundary conditions

Electrical double layer: revisit based on boundary conditions Electrcal double layer: revst based on boundary condtons Jong U. Km Department of Electrcal and Computer Engneerng, Texas A&M Unversty College Staton, TX 77843-318, USA Abstract The electrcal double layer

More information

Normally, in one phase reservoir simulation we would deal with one of the following fluid systems:

Normally, in one phase reservoir simulation we would deal with one of the following fluid systems: TPG4160 Reservor Smulaton 2017 page 1 of 9 ONE-DIMENSIONAL, ONE-PHASE RESERVOIR SIMULATION Flud systems The term sngle phase apples to any system wth only one phase present n the reservor In some cases

More information

EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES

EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES Manuel J. C. Mnhoto Polytechnc Insttute of Bragança, Bragança, Portugal E-mal: mnhoto@pb.pt Paulo A. A. Perera and Jorge

More information

FREQUENCY DISTRIBUTIONS Page 1 of The idea of a frequency distribution for sets of observations will be introduced,

FREQUENCY DISTRIBUTIONS Page 1 of The idea of a frequency distribution for sets of observations will be introduced, FREQUENCY DISTRIBUTIONS Page 1 of 6 I. Introducton 1. The dea of a frequency dstrbuton for sets of observatons wll be ntroduced, together wth some of the mechancs for constructng dstrbutons of data. Then

More information

Modeling of Dynamic Systems

Modeling of Dynamic Systems Modelng of Dynamc Systems Ref: Control System Engneerng Norman Nse : Chapters & 3 Chapter objectves : Revew the Laplace transform Learn how to fnd a mathematcal model, called a transfer functon Learn how

More information