Extracting turbulent spectral transfer from under-resolved velocity fields

Size: px
Start display at page:

Download "Extracting turbulent spectral transfer from under-resolved velocity fields"

Transcription

1 Extractng turbulent spectral transfer from under-resolved velocty felds Ru N, Greg A. Voth, and Ncholas T. Ouellette Ctaton: Physcs of Fluds (1994-present) 26, 117 (214); do: 1.163/ Vew onlne: Vew Table of Contents: Publshed by the AIP Publshng Artcles you may be nterested n Ansotropc energy transfers n quas-statc magnetohydrodynamc turbulence Phys. Plasmas 21, 1231 (214); 1.163/ Effect of shocklets on the velocty gradents n hghly compressble sotropc turbulence Phys. Fluds 23, 1213 (211); 1.163/ Rotatng helcal turbulence. I. Global evoluton and spectral behavor Phys. Fluds 22, 31 (21); 1.163/ Effects of vortex flaments on the velocty of tracers and heavy partcles n turbulence Phys. Fluds 18, 8172 (26); 1.163/ Energy spectrum n the enstrophy transfer range of two-dmensonal forced turbulence Phys. Fluds 13, 44 (21); 1.163/

2 PHYSICS OF FLUIDS 26, 117 (214) Extractng turbulent spectral transfer from under-resolved velocty felds Ru N, 1 Greg A. Voth, 2 and Ncholas T. Ouellette 1,a) 1 Department of Mechancal Engneerng and Materals Scence, Yale Unversty, New Haven, Connectcut 62, USA 2 Department of Physcs, Wesleyan Unversty, Mddletown, Connectcut 649, USA (Receved 17 August 214; accepted 8 October 214; publshed onlne 27 October 214) The strong nonlneartes n turbulent flows drve the transfer of energy and other quanttes between dfferent scales of moton. In three-dmensonal (3D) turbulence, ths transfer organzes nto the classc Rchardson Kolmogorov cascade of energy to small scales; n two-dmensonal (2D) turbulence, t leads to an nverse cascade of energy to large scales and a forward cascade of enstrophy to small scales. Drectly measurng ths spectral transfer s dffcult, partcularly n experments. Recently developed flterng technques allow spectral fluxes to be measured locally, but have been assumed to requre fnely resolved velocty felds that are typcally not avalable n 3D experments. Here we show, usng expermental data n 2D and the results of a 3D smulaton, that poorly resolved velocty felds can stll be used to extract nformaton about spectral transfer processes. Our results suggest new useful ways of analyzng data from turbulence experments wth lmted spatal resoluton. C 214 AIP Publshng LLC. [ I. INTRODUCTION In turbulent flow at hgh Reynolds numbers, the nonlnear term n the Naver Stokes equatons plays a key role n the dynamcs. Ths term ntroduces nteractons between wavenumber trads, whch spontaneously self-organze nto the Rchardson Kolmogorov energy cascade n three-dmensonal (3D) turbulence, drvng knetc energy from the large scales where t s njected nto the turbulence to small scales where t can be dsspated by vscous acton. 1 In two dmensonal (2D) turbulence, the addtonal conservaton of vortcty n the nvscd lmt produces a cascade of enstrophy from large to small scales whle reversng the drecton of the scale-to-scale flow of knetc energy. 2 4 Gross nformaton about the transfer of energy between scales s often obtaned va spectra or the velocty structure functons. 1, Nether of these quanttes, however, can easly provde the local spatotemporal detals of the spectral transfer processes at work n turbulence, as they requre averagng over tme and space. A recently developed alternatve approach s the so-called flter-space technque (FST) The essental dea of an FST s straghtforward. When one removes some degrees of freedom from a nonlnear equaton of moton by flterng (as s done n, for example, large-eddy smulaton), new terms appear that express the couplng between the degrees of freedom that are retaned and those that are suppressed. 8 By measurng these new terms, whch are related to the dfference between the fltered feld and the full feld and have prescrbed forms snce the equaton of moton s known, ths couplng can be drectly obtaned n a spatotemporally localzed way. When ths process s appled to, for example, the turbulent knetc energy, the couplng term descrbes the transfer of energy between scales. FSTs therefore allow the measurement of local spectral transfer propertes, and n prncple can be appled to any quantty n turbulence for whch the equaton of moton can be wrtten. a) Electronc mal: ncholas.ouellette@yale.edu /214/26(1)/117/9/$3. 26, C 214 AIP Publshng LLC

3 117-2 N, Voth, and Ouellette Phys. Fluds 26, 117 (214) In prncple, the result of an FST gves the net couplng between all the retaned scales and all the suppressed scales; thus, t would seem that an accurate measurement of the spectral transfer would requre a very hghly resolved velocty feld as an nput. Such hgh resolutons are typcally unavalable n experments, partcularly for 3D, hghly turbulent flows. In a turbulent cascade, however, we expect that spectral dynamcs are approxmately local, 14, 1 so that only modes near a gven length scale contrbute to spectral transfer across that scale. If localty holds, then as long we lmt our queres to scales somewhat larger than the resoluton of the velocty feld, FSTs should produce accurate measurements of spectral transfer, just as other tools such as velocty structure functons are expected to do. The queston remans, however, as to exactly how much resoluton we need to measure the spectral transfer across a gven scale wth suffcent accuracy. Prevous studes usng FSTs have prmarly nvestgated the double cascades of energy and enstrophy n 2D turbulence, 1, 11, 16, 17 where hgh spatal resoluton s much easer to obtan n experments. Here, we show n both 2D (usng expermental data) and 3D (usng numercal smulaton) turbulence that ndeed the localty assumpton holds, and that FSTs can gve accurate measurements of spectral transfer even when the velocty feld s poorly resolved. We quantfy the accuracy of the spectral energy transfer rates as a functon of the spatal resoluton of the velocty feld, and show that accurate measurements can be made for scales larger than about 6 tmes the resoluton of the velocty feld. For 3D turbulence, we also compare our FST results wth measurements of the secondand thrd-order Euleran velocty structure functons, and show that the FST gves comparable f not better estmates of the energy transfer rate. Our results have useful mplcatons for 3D turbulence experments, where hghly resolved velocty felds are typcally not avalable; even n these cases, FSTs can be appled and can provde good measurements of spectral processes. We begn n Sec. II wth a bref descrpton of our data sets and of the mechancs of applyng an FST. In Sec. III, we present our results for both 2D and 3D turbulence, and compare them wth more tradtonal estmates of spectral transfer processes. Fnally, n Sec. IV, we summarze our results and brefly dscuss ther mplcatons. II. METHODS A. Data sets To generate (weak) quas-2d turbulence, we drve a thn layer of electrolytc flud wth electromagnetc Lorentz forces. The detals of our experments are reported elsewhere; 18 2 here, we only brefly descrbe the setup. The workng flud s a thn layer of salt water (16% NaCl by mass n deonzed water), measurng cm 3, that rests on a glass plate. Beneath the glass les an array of permanent neodymum-ron-boron magnets arranged n a checkerboard pattern of alternatng polarty. To remove the possblty of long-range couplng between the partcles we use to measure the flow (see below) due to surface-tenson effects, we float a second layer of pure water (also. cm deep) above the salt water. To force the flud, we run d.c. electrc current laterally through the salt-water layer between a par of copper bar electrodes. The combnaton of electrc current and magnetc feld produces a Lorentz force that drves flud moton and njects energy nto the system. Thus, the energy njecton scale s approxmately set by the spacng between magnets L m, equal to 2.4 cm n our apparatus. We defne a Reynolds number Re = UL m /ν based on ths length scale, the root-mean-square (rms) velocty U, and the knematc vscosty ν of the workng flud. For the data shown here, Re = 18, suffcent to keep the flow nearly two-dmensonal. 18 The flow feld was measured usng partcle trackng velocmetry (PTV) by seedng the flud wth around 3 fluorescent polystyrene partcles wth dameters of 1 μm. The partcles have a specfc gravty of 1., between that of the salt water and the fresh water; they are therefore gravtatonally pnned at the nterface between the two flud layers. Images of the partcles n the central cm 2 of the apparatus were captured at a rate of 6 frames/s usng an IDT MotonPro M camera wth a resoluton of pxels. The hgh seedng densty allows us to measure very hghly resolved velocty felds. To ensure that the measured felds are as 2D as possble and to remove nose, we projected the measured veloctes onto a bass of numercally constructed streamfuncton egenmodes. 18

4 117-3 N, Voth, and Ouellette Phys. Fluds 26, 117 (214) The 3D turbulence data set we used was obtaned from drect numercal smulaton (DNS) results hosted at the Johns Hopkns Turbulence Database. 21 Ths DNS of forced sotropc turbulence was performed on a grd n a 2π 2π 2π doman; the Reynolds number based on the Taylor mcroscale s R λ = 433. Detals of the smulaton can be found n Ref. 21. B. Flter-space technques Implementng an FST requres suppressng small-scale varaton whle retanng the large scales. We accomplsh ths task by convolvng the quantty of nterest wth a spectral low-pass flter. As a concrete example, we defne the th component of the fltered velocty feld as u (r) (x) G(x x ; r)u (x )dx. (1) G here s the flter kernel, and the superscrpt (r) denotes a quantty fltered at a length scale r (where spatal scales smaller than r have been removed). In prncple, any functon that suppresses hghfrequency varaton wll work as a kernel for an FST. We choose to base our kernel on a Gaussan wth a Fourer-space wdth of 2π/r, followng prevous studes, 9 11 snce t has reasonable propertes n both real and Fourer space. Wth ths defnton, the equaton of moton for the fltered knetc energy E (r) = (1/2)[u (r) ] 2 s gven by 8 E (r) = J (r) ν u(r) u (r) (r), (2) t x x j x j where summaton s mpled over repeated ndces. The frst two terms on the rght-hand sde of ths equaton are analogous to terms n the equaton of moton for the full knetc energy, and represent spatal transport of knetc energy (the current J (r) ncludes contrbutons from flud advecton, pressure, and vscous dffuson) and drect vscous dsspaton. But the thrd term, gven by (r) = 1 [ ] [ (u u j ) (r) u (r) u (r) u (r) j + u ] (r) j, (3) 2 x j x s new. Ths term acts as a source or snk of knetc energy for the fltered velocty feld, and s expressed as the nner product of the fltered rate of stran wth a quantty that plays the role of a stress and that accounts for the momentum couplng between the fltered velocty feld and the small-scale part that we have removed. Wth our sgn conventon, (r) > mples energy flux from scales larger than r to those that are smaller, whle (r) < ndcates the opposte. C. Spatal resoluton To study the dependence of FSTs on the spatal resoluton of the velocty feld, we must be able to change ths resoluton. But for both the 2D and 3D data sets, the resoluton s fxed and we cannot ncrease t past ts gven value. What we can do, however, s artfcally reduce the resoluton by down-samplng the data sets. Here, we descrbe how we acheve ths down-samplng for both the expermental and numercal data. In a PTV experment lke ours, we obtan at each tme step velocty vectors for each tracked partcle. These partcles, however, are randomly seeded n the measurement area, and so the velocty feld s typcally sampled nonunformly. Ths nonunform samplng s undesrable when applyng Fourer flters as we do n the FST; thus, before flterng, we frst nterpolate our measured velocty vectors onto a regular grd. The spacng d of ths grd wll then set the effectve spatal resoluton of the velocty feld. Because the seedng densty s nhomogeneous, the fnest resoluton of the feld cannot be unquely defned; t should be, however, on the order of the typcal separaton between pars of partcles. To estmate the smallest d, we therefore calculated the dstrbuton of the nearest-neghbor dstance between partcle pars. More than 99.7% of pars were separated by dstances smaller than.11l m ; thus, we take d mn =.11L m to be the fnest spatal resoluton of the

5 117-4 N, Voth, and Ouellette Phys. Fluds 26, 117 (214) velocty feld, and down-sample the feld by nterpolatng the measured velocty vectors onto grds wth d > d mn. For the 3D DNS results, the data are already grdded, and so nterpolaton s not necessary; down-samplng, however, s. The orgnal smulaton was performed wth grd ponts spaced by d = 2.1η, where η s the Kolmogorov length scale. To down-sample the resoluton of the DNS data, we smply removed every nth grd pont, coarsenng the grd by a factor of n. We chose ths method rather than an nterpolaton-based scheme because t mmcs to some degree how data would be obtaned from a poorly sampled 3D partcle-trackng experment, where veloctes are known only at dscrete ponts and the velocty at other ponts can typcally not be obtaned by nterpolaton. We dd, however, only choose values of n that were commensurate wth 124, so that the boundary grd ponts remaned. Had we not done so, removng a few grd ponts near the boundares of the perodc box would have the effect of applyng a hgh-pass flter at a scale close to the box sze, whch would ntroduce complcated alasng effects. Here, we present data for fve values of n:4,8, 16, 32, and 64. Note that the way we ncrease d n both data sets s smlar to applyng a low-pass box flter n real space, whch wll certanly nduce some rngng at other scales n Fourer space. Ths effect, however, usually affects hgh frequences (small scale) much more strongly than t does low frequences (large scales). In Sec. III below, we use relatve error to provde emprcal, a posteror evdence that ths effect s not sgnfcant for scales larger than a few multples of d. III. RESULTS AND DISCUSSION A. 2D weak turbulence In Fg. 1(a), we plot the spatally averaged spectral energy flux (r) for a sngle velocty feld as measured n our quas-2d experment, for a range of dfferent spatal resolutons d. The closed symbol at the begnnng of each curve shows the smallest flter scale, chosen here to be r = 4d. Although the curves vary somewhat n detal, the overall trend s smlar. In partcular, all the resolutons capture the strong negatve flux for r 1.8L m, ndcatve of the nverse energy transfer that s expected for 2D turbulence. 1. (a) 3 (b) Π (r) (mm 2 /s 3 ) d=.11l m d=.16l m d=.22l m d=.27l m d=.32l m d=.38l m d=.43l m d=.48l m ξ (%) r/l m r ref /d FIG. 1. (a) Spatally averaged energy flux (r) as a functon of length scale r (scaled by the magnet spacng L m ) for one nstantaneous velocty feld n the quas-2d experment. Each curve shows data for a dfferent spatal resoluton d. The frst data pont on each curve, marked wth a symbol, shows the result for a resoluton of d = r/4. (b) Relatve error ξ,defnedby Eq. (4), for the same data as n (a) as a functon of the relatve spatal resoluton r/d. Each data pont shows results averaged over the range 2 < r/l m < 3, gvng an estmate of the value of ξ at r ref = 2.L m (see the text), and the error bars show the standard error n ξ.

6 117- N, Voth, and Ouellette Phys. Fluds 26, 117 (214) We can quantfy the effects of the lmted spatal resoluton by defnng a relatve error ξ = (r) (d) (r) (d mn ), (4) (r) (d mn ) where (r) (d mn ) s the mean energy flux for the velocty feld wth the fnest spatal resoluton d mn =.11L m. As wrtten, ξ s a functon of the flter scale r; thus, we must choose a value or range of r over whch to consder ξ. If our 2D turbulence were fully developed and we observed an extended nertal range (whch would be ndcated by a plateau over an apprecable range of r n Fg. 1(a)), a better choce would be to average ξ over ths plateau, as we do for the 3D case below. But because we (unsurprsngly, gven our Reynolds number) do not observe an nertal range, we cannot take ths approach. If we consdered ξ over the full range of r, we would ntroduce large fluctuatons snce ξ wll tend to become very large when (r) (d mn ). Thus, for the 2D data we report values of ξ averaged over the range 2 < r/l m < 3, where we observe net nverse energy transfer and where we emprcally observe that ξ fluctuates wthout a systematc dependence on r; ths procedure gves us an estmate of ξ(r = 2.L m ). The result s shown n Fg. 1(b) as a functon of r/d, the rato of the flter scale to the underlyng resoluton of the feld. Although ξ fluctuates, the trend of ts dependence on resoluton s clear. Unsurprsngly, ξ s large for very poorly resolved felds; but once r/d 7, ξ saturates to a value of about %. Ths result suggests that an FST can accurately provde an estmate of spectral transfer even for a poorly resolved data set, as long as one asks about scales larger than about 7 tmes the spatal resoluton of the data. In addton to average values, FSTs allow the determnaton of the spectral energy flux locally n both tme and space. It s thus natural to ask how reducng the resoluton affects the local values of the spectral flux. Here, we choose two dfferent spatal resolutons (d =.22 L m and.32 L m )to llustrate the effect of the resoluton d on the local spectral flux. In Fgs. 2(a) and 2(b), we show the same velocty feld but sampled at these two dfferent resolutons. The correspondng maps of the spectral flux calculated for the same flter length (r = 2.2L m, where the mean nverse energy flux s the most ntense) are shown n Fgs. 2(c) and 2(d). Despte the larger and somewhat fuzzer patches n Fg. 2(d), the spatal dstrbuton and magntude of the energy flux n the two fgures are nearly the same. Smlar effects are seen for other flter scales. In Fgs. 2(e) and 2(f), for example, we show the spectral energy flux for the same velocty feld at r = 3.2L m, where the mean spectral energy flux s much smaller. Comparng the results for the two values of r, t s clear that because the typcal (a) (c) 2 (e) 1 (b) 1 (d) (f) FIG. 2. (a) and (b) A sngle velocty feld, shown wth arrows, n the quas-2d flow wth a spatal resoluton of (a) d =.22L m and (b) d =.32L m. (c) and (d) Spatally resolved spectral energy flux for a flter length of r = 2.2L m and spatal resolutons of (c) d =.22L m and (d) d =.32L m. (e) and (f) Spatally resolved spectral energy flux for a flter length of r = 3.2L m and spatal resolutons of (e) d =.22L m and (f) d =.32L m. The color bars n (c) (f) show the magntude of the energy flux n unts of mm 2 /s 3.

7 117-6 N, Voth, and Ouellette Phys. Fluds 26, 117 (214) sze of a patch of energy flux of the same sgn s larger for r = 3.2L m (snce more of the small scales have been removed), these patches are well resolved for both resolutons (panels (e) and (f)). For the smaller flter length and lower resoluton (panel (d)), though, some patches are smlar n sze to d, and thus are poorly resolved. For Fg. 2(d), r = 6.8d, close to the pont below whch the relatve error ξ grows sgnfcantly. Thus, Fg. 2(d) gves a vsual confrmaton of the result seen n Fg. 1(b). To quantfy the effect of resoluton further, we plotted the ponts n panels (c) and (d) and panels (e) and (f) aganst each other and measured the R 2 coeffcents, whch would be unty f the data were dentcal. For panels (c) and (d), we found R 2 =.91, and for panels (e) and (f), we found R 2 =.97. B. 3D turbulence In Fg. 3(a), we show the spatally averaged spectral energy flux computed by applyng an FST to a sngle snapshot from the 3D DNS data at fve reduced resolutons d. As wth the 2D data, the curves agree well wth each other for all resolutons, agan (unsurprsngly) confrmng the quas-localty of the energy transfer n turbulence. In 3D turbulence, the mean energy transfer rate n the nertal range s expected to be equal to the mean energy dsspaton rate ɛ. Thus, t s natural to compare the spatally averaged energy flux as measured usng an FST to ɛ. InFg.3(a), we therefore also show the results of three tradtonal ways used to determne the energy dsspaton rate 22 for the fully resolved DNS data. By defnton, ɛ 2ν S j S j, where ν s the knematc vscosty and S j s the rate of stran. Ths defnton, however, s dffcult to apply n experments, snce resolvng the velocty gradents accurately requres sub-kolmogorov-scale resoluton of the velocty feld. Kolmogorov s scalng theory 1 provdes an alternatve way of determnng ɛ by relatng t to the Euleran velocty structure functons n the nertal range (that s, the range of scales that are far from both the scale of energy dsspaton and the scale of energy njecton). The velocty structure functons are defned to be the statstcal moments of the dfference n velocty between two ponts separated by a length scale r. The second-order structure functon s gven by D 2 (r) = [u(x + r) u(x)] 2. Projectng t along the separaton drecton gves the longtudnal second-order velocty structure functon D LL. In the nertal range, 1 D LL = C 2 (ɛr) 2/3, where C s expected to be a unversal constant, although wth some uncertanty Thus, the dsspaton rate s often estmated n experments by measurng D LL and scalng t approprately. In Fg. 3(a),we.1 (a) (b) 3D, N=128 3, d=17.2η 3D, N=64 3, d=34.3η 3D, N=32 3, d=68.7η ε, Π (r) (a.u.) , d=8.6η 128 3, d=17.2η 64 3, d=34.3η 32 3, d=68.7η 16 3, d=137.4η [D LL /(C 2 r 2/3 )] 3/2 D LLL /4r 2ν S j S j ξ (%) r/η 1 1 r/d FIG. 3. (a) The energy dsspaton rate ɛ andenergyflux (r) calculated va dfferent methods as a functon of length scale r, normalzed by the Kolmogorov length scale η. The symbols show (r) calculated from the same velocty feld wth dfferent spatal resolutons d, as shown n the legend. The legend also gves the number of DNS grd ponts used. The dashed and sold lnes show the compensated second- and thrd-order longtudnal Euleran structure functons, and the dashed-dotted lne shows the energy dsspaton rate obtaned drectly from the velocty gradent. (b) Relatve error ξ for (r) as shown n (a) for three representatve spatal resolutons d.

8 117-7 N, Voth, and Ouellette Phys. Fluds 26, 117 (214) plot the compensated second-order structure functon (that s, (D LL /C 2 ) 3/2 /r) computed from the DNS data, whch should be equal to ɛ n the nertal range. And ndeed, there s a reasonably extended range (though less than one decade) over whch the compensated structure functon has a constant value that agrees well wth the exact computaton of ɛ from the velocty gradents. The thrd-order longtudnal structure D LLL s also often used to extract ɛ; t s appealng gven that the so-called 4/ law, that D LLL = 4/ɛr n the nertal range, s exact. We therefore also show the compensated thrd-order structure functon n Fg. 3(a). However, as s often the case at fnte Reynolds numbers, the plateau regon s short, makng the estmate of ɛ somewhat unrelable. We note that ths structure functon s computed from only a sngle snapshot of the full DNS data set; usng more tme steps would gve a longer plateau regon. In general, though, the thrd-order structure functon s often dffcult to converge n experments, as t requres very large numbers of samples, and ts scalng regon can be qute small at low Reynolds numbers. Comparng our FST values to the three more tradtonal measurements of ɛ n Fg. 3(a), we can see that the FST performs qute well, even for reduced resoluton. Ths result s partcularly appealng snce, unlke the structure functons, there are no free parameters and no scalng assumptons necessary for nterpretng the FST results. To make the comparson quanttatve, we can agan compute the error ξ defned n Eq. (4); here, we use d mn = 8.6η as our baselne resoluton. As compared wth the 2D experment, the 3D DNS contans enough data ponts that we do not need to average over multple flter values to obtan smooth curves; thus, we do not do any averagng n computng ξ for the 3D case, and we can vary both r and d ndependently. Our results are shown n Fg. 3(b). The trend we observe s very smlar to the 2D case: the relatve error decreases wth r/d, but saturates to a constant, small value for r/d 6. For that reason, the symbols n Fg. 3(a) for the FST results wth dfferent resolutons d are computed only for flter lengths r 6d. In ths range, curves for dfferent d collapse well both wth each other and wth D LLL. In partcular, the estmate of ɛ from the plateau values of the FST s very accurate for d = 17.2η, and s adequate even for d = 34.3η, resolutons that are accessble n experments. As wth the 2D data set, we close by consderng the local spectral flux nformaton provded by FSTs at low resoluton n the 3D DNS. Fgure 4 shows slces through a sngle 3D spectral flux FIG. 4. Slces from the 3D feld of spectral energy flux (r) taken from an nstantaneous snapshot of the 3D DNS data. The color bar ndcates the magntude of the nondmensonal energy flux. Data are shown for four combnatons of flter length r and spatal resoluton d: (a)r = 2.2η and d = 8.6η (26 3 ); (b) r = 2.2η and d = 34.2η (64 3 ); (c) r = 342.1η and d = 8.6η (26 3 ); and (d) r = 342.1η and d = 34.2η (64 3 ).

9 117-8 N, Voth, and Ouellette Phys. Fluds 26, 117 (214) feld calculated for two dfferent flter lengths r and spatal resolutons d. For both flter scales, the spectral flux felds agree reasonably well despte the reduced resoluton, just as they dd n 2D. As above, we also computed the R 2 coeffcents for these data; for panels (a) and (b), we found R 2 =.7, and for panels (c) and (d) we found R 2 =.93. Even though the agreement s somewhat weaker than n the 2D case, our results ndcate that FSTs can be usefully appled to 3D turbulence even wth sgnfcantly under-resolved velocty felds, and that they are therefore potentally a very useful tool for analyzng partcle-trackng data. IV. CONCLUSIONS We have demonstrated that FSTs can be appled to study the average and local dynamcs of scale-to-scale energy transfer even when the velocty feld s not well resolved n both 2D and 3D turbulence. Our results show that as long as one restrcts the analyss to scales larger than about 6 tmes the spatal resoluton of the velocty feld, the spectral fluxes provded by an FST wll be accurate to wthn about %. Ths result confrms the expectaton that the energy cascade n turbulence s approxmately scale-local, snce the scales that are not present due to under-resoluton of the velocty feld do not contrbute strongly to the energy transfer processes at dstant scales. These results have useful mplcatons for partcle-trackng experments, where hghly resolved velocty felds are typcally not avalable. In partcular, snce they do not requre any assumptons about extended nertal ranges or scalng n order to estmate spectral transfer, FSTs may be a better way to estmate the energy dsspaton rate n turbulence experments than the more tradtonal velocty structure functons. FSTs are also lkely to be valuable for studyng spectral transfer n systems such as turbulence n non-newtonan fluds 26, 27 where hghly resolved numercal smulatons are very dffcult and expensve. ACKNOWLEDGMENTS Ths work was supported by the U.S. Natonal Scence Foundaton Grant Nos. DMR and DMS to Yale Unversty and DMR to Wesleyan Unversty. 1 A. N. Kolmogorov, The local structure of turbulence n ncompressble vscous flud for very large Reynolds numbers, Dokl. Akad. Nauk SSSR 3, (1941). 2 R. H. Krachnan, Inertal ranges n two-dmensonal turbulence, Phys. Fluds 1, (1967). 3 C. E. Leth, Dffuson approxmaton for two-dmensonal turbulence, Phys. Fluds 11, (1968). 4 G. K. Batchelor, Computaton of the energy spectrum n homegeneous two-dmensonal turbulence, Phys. Fluds 12, II-233 II-239 (1969). S. B. Pope, Turbulent Flows (Cambrdge Unversty Press, Cambrdge, England, 2). 6 M. Germano, Turbulence: The flterng approach, J. Flud Mech. 238, (1992). 7 S. Lu, C. Meneveau, and J. Katz, On the propertes of smlarty subgrd-scale models as deduced from measurements n a turbulent jet, J. Flud Mech. 27, (1994). 8 G. L. Eynk, Local energy flux and the refned smlarty hypothess, J. Stat. Phys. 78, (199). 9 M. K. Rvera, W. B. Danel, S. Y. Chen, and R. E. Ecke, Energy and enstrophy transfer n decayng two-dmensonal turbulence, Phys. Rev. Lett. 9, 142 (23). 1 S. Chen, R. E. Ecke, G. L. Eynk, X. Wang, and Z. Xao, Physcal mechansm of the two-dmensonal enstrophy cascade, Phys. Rev. Lett. 91, 2141 (23). 11 S. Chen, R. E. Ecke, G. L. Eynk, M. Rvera, M. Wan, and Z. Xao, Physcal mechansm of the two-dmensonal nverse energy cascade, Phys. Rev. Lett. 96, 842 (26). 12 M. Wan, Z. Xao, C. Meneveau, G. L. Eynk, and S. Chen, Dsspaton-energy flux correlatons as evdence for the Lagrangan energy cascade n turbulence, Phys. Fluds 22, 6172 (21). 13 Y. Lao and N. T. Ouellette, Spatal structure of spectral transport n two-dmensonal flow, J. Flud Mech. 72, (213). 14 G. Eynk, Localty of turbulent cascades, Physca D 27, (2). 1 J. A. Domaradzk and D. Carat, An analyss of the energy transfer and the localty of nonlnear nteractons n turbulence, Phys. Fluds 19, 8112 (27). 16 Z. Xao, M. Wan, S. Chen, and G. L. Eynk, Physcal mechansm of the nverse energy cascade of two-dmensonal turbulence: A numercal nvestgaton, J. Flud Mech. 619, 1 44 (29). 17 Y. Lao and N. T. Ouellette, Geometry of scale-to-scale energy and enstrophy transport n two-dmensonal flow, Phys. Fluds 26, 413 (214). 18 D. H. Kelley and N. T. Ouellette, Onset of three-dmensonalty n electromagnetcally drven thn-layer flows, Phys. Fluds 23, 413 (211).

10 117-9 N, Voth, and Ouellette Phys. Fluds 26, 117 (214) 19 D. H. Kelley and N. T. Ouellette, Spatotemporal persstence of spectral fluxes n two-dmensonal weak turbulence, Phys. Fluds 23, 1111 (211). 2 Y. Lao, D. H. Kelley, and N. T. Ouellette, Effects of forcng geometry on two-dmensonal weak turbulence, Phys. Rev. E 86, 3636 (212). 21 Y. L, E. Perlman, M. Wan, Y. Yang, C. Meneveau, R. Burns, S. Chen, A. Szalay, and G. Eynk, A publc turbulence database cluster and applcatons to study Lagrangan evoluton of velocty ncrements n turbulence, J. Turbul. 9, N31 (28). 22 M. Gbert, H. Xu, and E. Bodenschatz, Inertal effects on two-partcle relatve dsperson n turbulent flows, EPL 9, 64 (21). 23 K. R. Sreenvasan, On the unversalty of the Kolmogorov constant, Phys. Fluds 7, (199). 24 R. N and K.-Q. Xa, Kolmogorov constants for the second-order structure functon and the energy spectrum, Phys. Rev. E 87, 232 (213). 2 C.-C. Chen, D. B. Blum, and G. A. Voth, Effects of fluctuatng energy nput on the small scales n turbulence, J. Flud Mech. 737, 27 1 (213). 26 N. T. Ouellette, H. Xu, and E. Bodenschatz, Bulk turbulence n dlute polymer solutons, J. Flud Mech. 629, (29). 27 H.-D. X, E. Bodenschatz, and H. Xu, Elastc energy flux by flexble polymers n flud turbulence, Phys. Rev. Lett. 111, 241 (213).

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

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

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

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

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

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

A particle in a state of uniform motion remain in that state of motion unless acted upon by external force.

A particle in a state of uniform motion remain in that state of motion unless acted upon by external force. The fundamental prncples of classcal mechancs were lad down by Galleo and Newton n the 16th and 17th centures. In 1686, Newton wrote the Prncpa where he gave us three laws of moton, one law of gravty,

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

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

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

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

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

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

4 Analysis of Variance (ANOVA) 5 ANOVA. 5.1 Introduction. 5.2 Fixed Effects ANOVA

4 Analysis of Variance (ANOVA) 5 ANOVA. 5.1 Introduction. 5.2 Fixed Effects ANOVA 4 Analyss of Varance (ANOVA) 5 ANOVA 51 Introducton ANOVA ANOVA s a way to estmate and test the means of multple populatons We wll start wth one-way ANOVA If the populatons ncluded n the study are selected

More information

Turbulent Transport in Single-Phase Flow. Peter Bernard, University of Maryland

Turbulent Transport in Single-Phase Flow. Peter Bernard, University of Maryland Turbulent Transport n Sngle-Phase Flow Peter Bernard, Unversty of Maryland Assume that our goal s to compute mean flow statstcs such as U and One can ether: 1 u where U Pursue DNS (.e. the "honest" approach)

More information

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

Physics 5153 Classical Mechanics. Principle of Virtual Work-1 P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal

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

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

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

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

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

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

Department of Statistics University of Toronto STA305H1S / 1004 HS Design and Analysis of Experiments Term Test - Winter Solution

Department of Statistics University of Toronto STA305H1S / 1004 HS Design and Analysis of Experiments Term Test - Winter Solution Department of Statstcs Unversty of Toronto STA35HS / HS Desgn and Analyss of Experments Term Test - Wnter - Soluton February, Last Name: Frst Name: Student Number: Instructons: Tme: hours. Ads: a non-programmable

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

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

Chapter 11: Simple Linear Regression and Correlation

Chapter 11: Simple Linear Regression and Correlation Chapter 11: Smple Lnear Regresson and Correlaton 11-1 Emprcal Models 11-2 Smple Lnear Regresson 11-3 Propertes of the Least Squares Estmators 11-4 Hypothess Test n Smple Lnear Regresson 11-4.1 Use of t-tests

More information

Lecture 4. Macrostates and Microstates (Ch. 2 )

Lecture 4. Macrostates and Microstates (Ch. 2 ) Lecture 4. Macrostates and Mcrostates (Ch. ) The past three lectures: we have learned about thermal energy, how t s stored at the mcroscopc level, and how t can be transferred from one system to another.

More information

EXPERIMENTAL STUDY OF NEAR WALL TURBULENCE USING PIV

EXPERIMENTAL STUDY OF NEAR WALL TURBULENCE USING PIV EUROMECH 411 Rouen, 9-31 May EXPERIMENTAL STUDY OF NEAR WALL TURBULENCE USING PIV J. Carler, J. M. Foucaut and M. Stanslas LML URA 1441, Bv Paul Langevn, Cté Scentfque, 59655 Vlleneuve d'ascq Cedex, France

More information

PHYS 705: Classical Mechanics. Newtonian Mechanics

PHYS 705: Classical Mechanics. Newtonian Mechanics 1 PHYS 705: Classcal Mechancs Newtonan Mechancs Quck Revew of Newtonan Mechancs Basc Descrpton: -An dealzed pont partcle or a system of pont partcles n an nertal reference frame [Rgd bodes (ch. 5 later)]

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

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

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

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

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

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

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

arxiv:cs.cv/ Jun 2000

arxiv:cs.cv/ Jun 2000 Correlaton over Decomposed Sgnals: A Non-Lnear Approach to Fast and Effectve Sequences Comparson Lucano da Fontoura Costa arxv:cs.cv/0006040 28 Jun 2000 Cybernetc Vson Research Group IFSC Unversty of São

More information

Supplemental Material: Causal Entropic Forces

Supplemental Material: Causal Entropic Forces Supplemental Materal: Causal Entropc Forces A. D. Wssner-Gross 1, 2, and C. E. Freer 3 1 Insttute for Appled Computatonal Scence, Harvard Unversty, Cambrdge, Massachusetts 02138, USA 2 The Meda Laboratory,

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

Numerical Transient Heat Conduction Experiment

Numerical Transient Heat Conduction Experiment Numercal ransent Heat Conducton Experment OBJECIVE 1. o demonstrate the basc prncples of conducton heat transfer.. o show how the thermal conductvty of a sold can be measured. 3. o demonstrate the use

More information

First Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force.

First Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force. Secton 1. Dynamcs (Newton s Laws of Moton) Two approaches: 1) Gven all the forces actng on a body, predct the subsequent (changes n) moton. 2) Gven the (changes n) moton of a body, nfer what forces act

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

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

SIO 224. m(r) =(ρ(r),k s (r),µ(r))

SIO 224. m(r) =(ρ(r),k s (r),µ(r)) SIO 224 1. A bref look at resoluton analyss Here s some background for the Masters and Gubbns resoluton paper. Global Earth models are usually found teratvely by assumng a startng model and fndng small

More information

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown

More information

Robert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations

Robert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations Quantum Physcs 量 理 Robert Esberg Second edton CH 09 Multelectron atoms ground states and x-ray exctatons 9-01 By gong through the procedure ndcated n the text, develop the tme-ndependent Schroednger equaton

More information

Indeterminate pin-jointed frames (trusses)

Indeterminate pin-jointed frames (trusses) Indetermnate pn-jonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all

More information

Chapter 5 Multilevel Models

Chapter 5 Multilevel Models Chapter 5 Multlevel Models 5.1 Cross-sectonal multlevel models 5.1.1 Two-level models 5.1.2 Multple level models 5.1.3 Multple level modelng n other felds 5.2 Longtudnal multlevel models 5.2.1 Two-level

More information

Lecture 14: Forces and Stresses

Lecture 14: Forces and Stresses The Nuts and Bolts of Frst-Prncples Smulaton Lecture 14: Forces and Stresses Durham, 6th-13th December 2001 CASTEP Developers Group wth support from the ESF ψ k Network Overvew of Lecture Why bother? Theoretcal

More information

The Feynman path integral

The Feynman path integral The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space

More information

Workshop: Approximating energies and wave functions Quantum aspects of physical chemistry

Workshop: Approximating energies and wave functions Quantum aspects of physical chemistry Workshop: Approxmatng energes and wave functons Quantum aspects of physcal chemstry http://quantum.bu.edu/pltl/6/6.pdf Last updated Thursday, November 7, 25 7:9:5-5: Copyrght 25 Dan Dll (dan@bu.edu) Department

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

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

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

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

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

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

/ n ) are compared. The logic is: if the two

/ n ) are compared. The logic is: if the two STAT C141, Sprng 2005 Lecture 13 Two sample tests One sample tests: examples of goodness of ft tests, where we are testng whether our data supports predctons. Two sample tests: called as tests of ndependence

More information

2010 Black Engineering Building, Department of Mechanical Engineering. Iowa State University, Ames, IA, 50011

2010 Black Engineering Building, Department of Mechanical Engineering. Iowa State University, Ames, IA, 50011 Interface Energy Couplng between -tungsten Nanoflm and Few-layered Graphene Meng Han a, Pengyu Yuan a, Jng Lu a, Shuyao S b, Xaolong Zhao b, Yanan Yue c, Xnwe Wang a,*, Xangheng Xao b,* a 2010 Black Engneerng

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

Probabilistic method to determine electron correlation energy

Probabilistic method to determine electron correlation energy Probablstc method to determne electron elaton energy T.R.S. Prasanna Department of Metallurgcal Engneerng and Materals Scence Indan Insttute of Technology, Bombay Mumba 400076 Inda A new method to determne

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

EEE 241: Linear Systems

EEE 241: Linear Systems EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they

More information

Solution Thermodynamics

Solution Thermodynamics Soluton hermodynamcs usng Wagner Notaton by Stanley. Howard Department of aterals and etallurgcal Engneerng South Dakota School of nes and echnology Rapd Cty, SD 57701 January 7, 001 Soluton hermodynamcs

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

PHYS 705: Classical Mechanics. Calculus of Variations II

PHYS 705: Classical Mechanics. Calculus of Variations II 1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary

More information

Experiment 1 Mass, volume and density

Experiment 1 Mass, volume and density Experment 1 Mass, volume and densty Purpose 1. Famlarze wth basc measurement tools such as verner calper, mcrometer, and laboratory balance. 2. Learn how to use the concepts of sgnfcant fgures, expermental

More information

Professor Terje Haukaas University of British Columbia, Vancouver The Q4 Element

Professor Terje Haukaas University of British Columbia, Vancouver  The Q4 Element Professor Terje Haukaas Unversty of Brtsh Columba, ancouver www.nrsk.ubc.ca The Q Element Ths document consders fnte elements that carry load only n ther plane. These elements are sometmes referred to

More information

Georgia Tech PHYS 6124 Mathematical Methods of Physics I

Georgia Tech PHYS 6124 Mathematical Methods of Physics I Georga Tech PHYS 624 Mathematcal Methods of Physcs I Instructor: Predrag Cvtanovć Fall semester 202 Homework Set #7 due October 30 202 == show all your work for maxmum credt == put labels ttle legends

More information

Supporting Information

Supporting Information Supportng Informaton The neural network f n Eq. 1 s gven by: f x l = ReLU W atom x l + b atom, 2 where ReLU s the element-wse rectfed lnear unt, 21.e., ReLUx = max0, x, W atom R d d s the weght matrx to

More information

THE CURRENT BALANCE Physics 258/259

THE CURRENT BALANCE Physics 258/259 DSH 1988, 005 THE CURRENT BALANCE Physcs 58/59 The tme average force between two parallel conductors carryng an alternatng current s measured by balancng ths force aganst the gravtatonal force on a set

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

A NEW FILTERED DYNAMIC SUBGRID-SCALE MODEL FOR LARGE EDDY SIMULATION OF INDOOR AIRFLOW

A NEW FILTERED DYNAMIC SUBGRID-SCALE MODEL FOR LARGE EDDY SIMULATION OF INDOOR AIRFLOW A NEW FILTERED DYNAMIC SUBGRID-SCALE MODEL FOR LARGE EDDY SIMULATION OF INDOOR AIRFLOW We Zhang and Qngyan Chen Buldng Technology Program Massachusetts Insttute of Technology 77 Mass. Ave., Cambrdge, MA

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

Lecture 10 Support Vector Machines II

Lecture 10 Support Vector Machines II Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed

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

A Robust Method for Calculating the Correlation Coefficient

A Robust Method for Calculating the Correlation Coefficient A Robust Method for Calculatng the Correlaton Coeffcent E.B. Nven and C. V. Deutsch Relatonshps between prmary and secondary data are frequently quantfed usng the correlaton coeffcent; however, the tradtonal

More information

Lagrangian Field Theory

Lagrangian Field Theory Lagrangan Feld Theory Adam Lott PHY 391 Aprl 6, 017 1 Introducton Ths paper s a summary of Chapter of Mandl and Shaw s Quantum Feld Theory [1]. The frst thng to do s to fx the notaton. For the most part,

More information

Brownian-Dynamics Simulation of Colloidal Suspensions with Kob-Andersen Type Lennard-Jones Potentials 1

Brownian-Dynamics Simulation of Colloidal Suspensions with Kob-Andersen Type Lennard-Jones Potentials 1 Brownan-Dynamcs Smulaton of Collodal Suspensons wth Kob-Andersen Type Lennard-Jones Potentals 1 Yuto KIMURA 2 and Mcho TOKUYAMA 3 Summary Extensve Brownan-dynamcs smulatons of bnary collodal suspenton

More information

Quantitative Discrimination of Effective Porosity Using Digital Image Analysis - Implications for Porosity-Permeability Transforms

Quantitative Discrimination of Effective Porosity Using Digital Image Analysis - Implications for Porosity-Permeability Transforms 2004, 66th EAGE Conference, Pars Quanttatve Dscrmnaton of Effectve Porosty Usng Dgtal Image Analyss - Implcatons for Porosty-Permeablty Transforms Gregor P. Eberl 1, Gregor T. Baechle 1, Ralf Weger 1,

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

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

Psychology 282 Lecture #24 Outline Regression Diagnostics: Outliers

Psychology 282 Lecture #24 Outline Regression Diagnostics: Outliers Psychology 282 Lecture #24 Outlne Regresson Dagnostcs: Outlers In an earler lecture we studed the statstcal assumptons underlyng the regresson model, ncludng the followng ponts: Formal statement of assumptons.

More information

Generalized Linear Methods

Generalized Linear Methods Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set

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

Statistics Chapter 4

Statistics Chapter 4 Statstcs Chapter 4 "There are three knds of les: les, damned les, and statstcs." Benjamn Dsrael, 1895 (Brtsh statesman) Gaussan Dstrbuton, 4-1 If a measurement s repeated many tmes a statstcal treatment

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

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

(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

CONTRAST ENHANCEMENT FOR MIMIMUM MEAN BRIGHTNESS ERROR FROM HISTOGRAM PARTITIONING INTRODUCTION

CONTRAST ENHANCEMENT FOR MIMIMUM MEAN BRIGHTNESS ERROR FROM HISTOGRAM PARTITIONING INTRODUCTION CONTRAST ENHANCEMENT FOR MIMIMUM MEAN BRIGHTNESS ERROR FROM HISTOGRAM PARTITIONING N. Phanthuna 1,2, F. Cheevasuvt 2 and S. Chtwong 2 1 Department of Electrcal Engneerng, Faculty of Engneerng Rajamangala

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

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

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

Frequency dependence of the permittivity

Frequency dependence of the permittivity Frequency dependence of the permttvty February 7, 016 In materals, the delectrc constant and permeablty are actually frequency dependent. Ths does not affect our results for sngle frequency modes, but

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

PY2101 Classical Mechanics Dr. Síle Nic Chormaic, Room 215 D Kane Bldg

PY2101 Classical Mechanics Dr. Síle Nic Chormaic, Room 215 D Kane Bldg PY2101 Classcal Mechancs Dr. Síle Nc Chormac, Room 215 D Kane Bldg s.ncchormac@ucc.e Lectures stll some ssues to resolve. Slots shared between PY2101 and PY2104. Hope to have t fnalsed by tomorrow. Mondays

More information

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4) I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes

More information

Topic 23 - Randomized Complete Block Designs (RCBD)

Topic 23 - Randomized Complete Block Designs (RCBD) Topc 3 ANOVA (III) 3-1 Topc 3 - Randomzed Complete Block Desgns (RCBD) Defn: A Randomzed Complete Block Desgn s a varant of the completely randomzed desgn (CRD) that we recently learned. In ths desgn,

More information