GS and SGS Coherent Structures in Homogeneous Isotropic Turbulence

Size: px
Start display at page:

Download "GS and SGS Coherent Structures in Homogeneous Isotropic Turbulence"

Transcription

1 15 th CFD Symposum <C04-> GS and SGS Coherent Structures n Homogeneous Isotropc Turbulence Ο M. Ashraf Uddn Intellgent Modelng Laboratory, The Unersty of Tokyo, Yayo, Bunkyo-ku, Tokyo , Japan. E-mal: auddn@cebeer.s.u-tokyo.ac.jp Nobuyuk Tanguch Insttute of Industral Scence, The Unersty of Tokyo, Komaba, Meguro-ku, Tokyo , Japan. E-mal: ntan@s.u-tokyo.ac.jp Mamoru Tanahash Dept. of Mechancal and Aerospace Engneerng, Tokyo Insttute of Technology, -1-1 Ookayama, Meguro-ku, Tokyo , Japan. E-mal: mtanahas@mes.ttech.ac.jp Tosho Myauch Dept. of Mechancal and Aerospace Engneerng, Tokyo Insttute of Technology, -1-1 Ookayama, Meguro-ku, Tokyo , Japan. E-mal: tmyauch@mes.ttech.ac.jp Tosho Kobayash Insttute of Industral Scence, The Unersty of Tokyo, Komaba, Meguro-ku, Tokyo , Japan. E-mal: kobaya@s.u-tokyo.ac.jp Abstract: GS (grd-scale) and SGS (subgrd-scale) coherent structures n homogeneous sotropc turbulence are dentfed by a pror test. GS and SGS elocty felds are obtaned by flterng the elocty feld for Re λ =64.9 usng two classcal flters for LES: Gaussan flter and sharp cutoff flter, and the most mportant flter wdth s consdered as the length of Kolmogoro mcroscale n the feld wth a constant multplcaton. Coherent structures n the GS feld,.e., structures larger than flter wdth are consdered as grd-scale coherent structures and others are subgrd-scale coherent structures. Second narant Q of the elocty gradent tensor s used for dentfcaton of these structures n GS and SGS felds. By sualzng the contour surfaces of second narant Q, t s shown that GS and SGS felds tself contan lots of dstnct tube-lke coherent structures n homogeneous sotropc turbulence. The characterstc of GS and SGS coherent structures s somewhat smlar to structures, whch ndcates that the feld may contan mult-scale structures n turbulent flow. 1 INTRODUCTION Drect numercal and large-eddy smulatons ( and LES) hae been wdely used to study the physcs of turbulence. Howeer, drect numercal smulaton for hgh Reynolds number flow requres formdable computng power, and s only possble for low Reynolds numbers. Generally, ndustral, natural or expermental confguratons nole Reynolds numbers that are far too large to allow drect numercal smulaton, and the only possble method s large-eddy smulaton. Turbulent flow s completely nonlnear and complex and t has been shown n recent years that turbulent flows contan arous types of ortcal structures, more strctly coherent structures, but the lack of the proper knowledge about these structures preents the deelopment of turbulence theory and turbulent model. Durng the past decades, the studes on the coherent structures n turbulence hae been the subjects of consderable nterest among turbulence researchers. The study on these coherent structures s promsng not only for understandng turbulence phenomena such as entranment and mxng, scalar dsspaton, heat and mass transfer, chemcal reacton and combuston, drag and aerodynamc nose generaton, but also for modelng of turbulence. In the theoretcal study, t s beleed that tube-lke structure s a type of eddy or ortex, whch s the canddate of fne scale structure, partcularly, n the small-scale motons n turbulence. (1)()(3)(4)(5) Nowadays, from drect numercal smulaton of turbulence, (6)(7)(8)(9)(10)(11)(1)(13) fne scale tube-lke coherent structures n homogeneous turbulence are obsered, and the sualzaton of these small-scale structures n the turbulent flow becomes possble. In the recent studes, by drect use of local flow pattern, (10)(11) the cross-sectons of tube-lke coherent fne scale structures are nestgated from database of homogeneous sotropc turbulence n whch the crosssectons are selected to nclude the local maxmum of second narant of the elocty gradent tensor on the axs of the fne scale tube-lke coherent structures. In these studes, they hae shown that mean dameter of the coherent fne scale structures s about 10 tmes of Kolmogoro mcroscale (η) and the maxmum of mean azmuthal elocty s about a half of r.m.s of elocty fluctuatons (u rms ) and that the Reynolds numbers dependence of these characters s ery weak. The same analyses hae been appled to turbulent mxng layer (14) and showed that the characterstcs of tube-lke coherent structures n homogeneous sotropc turbulence and fully deeloped turbulent mxng layer obey the same scalng law. Snce the educed fne scale structures n ther study hae smlar mean azmuthal elocty profles and dstnct axes, they descrbed these structures as coherent fne scale structures n turbulence. The characterstcs of ortcal structures n turbulent channel flows (15) and MHD turbulence (16) also show the smlar behaor of tube-lke structures n homogeneous sotropc turbulence. These results suggest that the exstence of coherent fne scale structure n turbulence s unersal. 1

2 Although s the most exact approach to turbulence smulaton but too expense and s possble relately n smple flow felds whle large-eddy smulaton (LES) s less expense and can smulate ery complex flow felds n turbulence. Wth LES method, large-scale moton s drectly calculated but small-scale needs to be modeled by subgrdscale (SGS) model. Concernng the subgrd-scale model, t seems qute mportant to know what happened n the fltered feld for LES from the actual turbulence ncludng appearance or dsappearance of the tube-lke ortcal structures n resoled or unresoled feld. Therefore, the purpose of ths study s to flter elocty feld for obtanng the GS (grd-scale) and SGS (subgrd-scale) elocty felds n turbulent flow usng classcal flters for LES. Then we dentfy the grd-scale and subgrd-scale coherent structures and dscuss the statstcs of GS and SGS coherent structures n homogeneous sotropc turbulence. LARGE-EDDY SIMULATION.1 Data Base In ths study, data of decayng homogeneous sotropc turbulence has been used, whch s conducted by Tanahash et al. (10) and s calculated by usng 18 3 grd ponts. Reynolds number based on u rms and Taylor mcroscale, λ of the data s Re λ = GS and SGS Velocty Felds To obtan the grd-scale (GS) and gubgrd-scale (SGS) elocty felds, we hae drectly fltered the aboe elocty feld usng two classcal flters: Gaussan flter and sharp cutoff flter for LES. In LES, a elocty component u can be decomposed nto two components, one component s n the range of low wae-number of energy spectrum, large scale of moton, called GS component and s denoted by u, and the other component s n the range of hgh waenumber of energy spectrum, small scale of moton, called SGS component and s denoted by u. Ther relaton can be expressed as: u = u + u (1) The flterng s represented mathematcally n physcal space as a conoluton product. (17) The fltered part u of the arable u s defned formally by the relaton: u ( x) G( x x ) u( x ) dx =,, () D n whch G s flter kernel functon and s the flter wdth n -drecton. The dual defnton n the Fourer space û k of u(x) can be obtaned by multplyng the spectrum ( ) by the spectrum Ĝ ( k) of the kernel G(x) such that, ˆ ( k) = Gˆ ( k) uˆ ( k), k = 0, ± 1, ±,... u (3) The functon Ĝ s the transfer functon assocated wth the kernel G..3 Classcal Flters for LES In ths study, two classcal flters are used for performng the spatal scale separaton. For a flter wdth n -drecton, these flters n physcal space (PS) and Fourer space (FS) are wrtten as follows: (I) Gaussan flter: G 6 6x ( ) = x exp π ( ) ( ) k k = exp (PS) (4a) G ˆ 4 (FS) (4b) (II) Sharp cutoff flter: πx sn G( x ) = (PS) (5a) πx π 1, k Gˆ ( k ) = (FS) (5b) π 0, k > Usng the aboe two flters for LES n the Fourer space, one set of data s fltered and the exact GS elocty feld, u are obtaned. After generatng u, we subtract the GS elocty feld from elocty feld and obtan the SGS elocty feld. That s, the SGS elocty feld can be obtaned by the relaton (from Eq.1): u = u u (6) Flter wdth plays ery mportant role wth flter functons n ths process. The characterstc flter wdth s commonly used as the length, approxmately proportonal to the grd nteral n the preous researches. (18)(19)(0) The structures represented by the GS and SGS eloctes consequently depend both on the grd nteral and on the type of flter employed. In the preous studes, (10)(11) t s shown that the mean dameter of the coherent fne scale eddy s about 10 tmes of Kolmogoro mcroscale (η) n turbulent flows. Therefore, n ths study, the most mportant flter wdth,, s consdered as the length of Kolmogoro mcroscale n the feld wth a constant multplcaton. Snce we are dealng wth homogeneous sotropc turbulence, the flter wdth s same n each drecton and hereafter t s denoted by..4 Profles of the, GS and SGS Velocty Felds In the preous study, (1) usng seeral flter wdths and dfferent Reynolds numbers, t s shown that flter wdth depends on the Reynolds number of the flow, regardless of the flter employed. In ths study, to obtan the GS and SGS elocty felds from feld for Re λ =64.9, we hae consdered flter wdth =10η, where η s obtaned from the elocty feld. In order to understand the GS elocty feld from the elocty feld, we compared elocty wth fltered elocty n Fg.1. For ths purpose we hae randomly chosen one-dmensonal elocty profle n x 1 drecton and hae plotted u 1( x 1 ) and ( x ) 1 1 u for =10η n Fg.1. In all cases we

3 .0 GS-Gaussan u Fg. 1 Sample of fltered and unfltered elocty felds Slope -5/3 GS-Gaussan E(k) E(k) kη 10 1 Slope -5/3 SGS-Gaussan S 10 - kη 10 1 Fg. Three-dmensonal energy spectra of elocty fluctuatons. and GS felds, and SGS felds. obtaned these one-dmensonal profles for the aboe two flter functons: Gaussan and sharp cutoff flter. Although =10η s a small alue, but the profles n Fg.1 suggest that the generaton of GS elocty feld as well as SGS elocty fled from ths Re λ case s well usng both flter functons. Three-dmensonal energy spectra of fltered and unfltered elocty fluctuatons for ths flter wdth, =10η are presented n Fg., whch s calculated by the defnton wrtten as follows: E ( k ) = uˆ ( ) uˆ ( k) 1 1 k < k k+ 1 k (7) In each case, the GS and SGS spectra are obtaned by usng Gaussan and sharp cutoff flter, and then compared wth the spectrum. The spectrum shows the power decay close to k -5/3. In the preous study, (10) usng database, t s shown that the energy dsspaton rate s domnated by the fne scale eddes n homogeneous sotropc turbulence, whch s beyond the dscusson of ths paper. Wth the sharp cutoff flter, the SGS felds contan the elocty due to all the structures wth wae number k > π. On the other hand, full range of wae numbers contrbutes to the SGS elocty n the case of Gaussan flter. Fg. reeals ths dfference n behaor for GS and 3 Fg. 3 Contour surfaces of the second narant of the elocty gradent tensor (Q=3) n the feld. Vsualzed regon s the whole calculaton doman. Second narant s normalzed by Kolmogoro mcroscale and r.m.s of elocty fluctuatons. SGS spectra usng dfferent flters for LES. Ths fgure clearly ndcates that, when sharp cutoff flter s used, the GS spectrum exactly collapsed wth spectrum n the range of low wae numbers ( k π ), whle SGS spectrum s restrcted n the hgh wae number range, and the contrbuton of subgrd scale s entrely due to the hgh wae numbers. The behaor of these spectra also confrms the accuracy of the flterng process. When Gaussan flter s used, the subgrd scales account for a fracton of the total () knetc energy. Moreoer, SGS energy has a sgnfcant contrbuton from the low wae numbers (.e., the large scales) as we can see n Fg.. Howeer, wth ths small, t seems that GS felds contrbute to the large part of energy dsspaton rate for all flter functons. 3 GS AND SGS COHERENT STRUCTURES In order to dscuss about GS and SGS coherent structures, we notce the tube-lke coherent fne scale eddy by sualzaton of flows n the feld. The concept usually assocated wth an eddy s that of a regon n the flow where the flud elements are rotatng around a set of ponts. Identfcaton of the eddy or ortex from /LES database s a ery dffcult and complex task, requrng consderable computatonal efforts wth proper dentfcaton method. There are seeral methods for dentfcaton of the ortcal structures n turbulence wth sgnfcant dfferences () and most of them show threshold dependence. As we dscussed n the ntroducton, drect local flow pattern (3) can educe coherent structures n (10)(11)(1)(14) seeral flow felds, whch shows unersal characterstcs n turbulence. In our preous study, (13) usng the aboe method we hae dentfed the coherent fne scale eddes and ts axes wthout usng any thresholds and then dscussed the spatal dstrbuton of coherent fne scale eddes by sualzaton of axes n homogeneous sotropc turbulence. In ths study, we use the same method to dscuss about GS and SGS coherent structures n homogeneous sotropc turbulence. Fg.3 shows the contour surfaces of normalzed second narant of the elocty gradent tensor Q n the feld for Re λ =64.9. The second narant of the elocty gradent

4 I I I(c) II II II(c) Fg. 4 Contour surfaces of the second narant of the elocty gradent tensor n GS and SGS felds obtaned from feld usng (I) Gaussan flter and (II) sharp cutoff flter. Vsualzed regon s same as n Fg.3. GS feld (Q=), SGS feld (Q=1) (c) GS (Q=) & SGS (Q=1) felds tensor s defned as: Q 1 = ( S S W W ) j j j j, (8) where S j and W j are the symmetrc and asymmetrc part of the elocty gradent tensor A j. In Fg.3, the sualzed regon s whole calculaton doman and the leel of the sosurface s selected to be Q =3. Hereafter, denotes the normalzaton by Kolmogoro mcroscale η and root mean square of elocty fluctuatons, u rms, and for all cases η and u rms are obtaned from the feld. The normalzaton of Q by η and u rms s due to the fact that the dameter and the maxmum azmuthal elocty of tube-lke fne scale structures can be scaled by η and u rms (10) Fg.3 shows that lots of coherent tube-lke structures are randomly orented n homogeneous sotropc turbulence. Howeer, f we ncrease or decrease the alue of Q, we can also show dstnct tube-lke structures n turbulence, lttle bt dfferent from Fg.3, whch means the sualzaton of fne scale structures sgnfcantly depends on the alue of threshold of Q. (13) Fg.4 shows the contour surfaces of second narant of the elocty gradent tensor Q n GS and SGS felds obtaned by usng Gaussan and sharp cutoff flters wth flter wdth =10η. The sualzed regon and ewpont n Fg.4 s same as n Fg.3 for all cases. The leel of the 4 sosurface for all cases n Fg.4 s selected to be Q = for GS felds and Q =1 for SGS felds. As we dscussed aboe, the sualzaton of coherent structures depends on the threshold alue of Q and we do not concern the strength of the structures, therefore, n ths sualzaton we consdered these dfferent alues for Q for dfferent felds only to show the ortcal structures n GS and SGS felds by sualzaton. Fg.4 clearly ndcates that GS and SGS felds contan lots of dstnct tube-lke structures somewhat smlar to felds, whch can best be defned as coherent structures or eddes n turbulence. (13) It s also clear that the sze or length of GS structures n all cases seems to be larger than SGS structures. It s known from the classcal dea of flud dynamcs that seeral small-scale structure together form a large-scale structure,.e., seeral small (SGS) structures entrely le nsde a large (GS) structure n the order of ts sze. Howeer, Fg.4(c) clearly ndcates that GS (green) and SGS (red) structures are qute dstnct and unque n turbulence. 4. CHARACTERISTICS OF GS AND SGS COHERENT STRUCTURES 4.1 Identfcaton Scheme Snce for sualzaton of flows n Fgs.3-4, we use poste second narant Q and we can see the exstence of

5 many tube-lke coherent structures n, GS and SGS felds n homogeneous sotropc turbulence. Obously these tube-lke structures contan at least one local maxmum of Q on ts axs. Wth ths local Q maxmal, Tanahash et al (10) hae clarfed the statstcs of coherent fne scale structures n homogeneous sotropc turbulence. Ths method can best be used to study on the GS and SGS coherent structures n turbulence. Usng ths method, we can obtan a pont on the cross-secton on the axs of coherent structures wth local Q maxmum. The dentfcaton method conssts of the followng steps: Step : Ealuaton of Q at each collocaton pont from the results of, GS and SGS felds. Step : Probablty of exstence of poste local maxmum of Q near the collocaton ponts s ealuated at each collocaton pont from Q dstrbuton. Because the case that a local maxmum of Q concdes wth a collocaton pont s ery rare, t s necessary to defne probablty on collocaton ponts. Step (c): Collocaton ponts wth non-zero probablty are selected to surey actual maxma of Q. Locatons of maxmal Q are determned wthn the accuracy of 10-6 n terms of relate error of Q by applyng a 4 th order Lagrange-nterpolaton polynomal to, GS and SGS data. Step (d): A cylndrcal coordnate system (r,, z) s consdered by settng the maxmal pont as the orgn. The z drecton s selected to be parallel to the ortcty ector at the maxmal pont. The elocty ectors are projected on ths coordnate and azmuthal elocty u s calculated. Step (e): Pont that has small arance n azmuthal elocty compared wth the surroundngs s determned. If the azmuthal eloctes at r =1/5 computatonal grd space show same sgn for all, that pont s dentfed as the center of the swrlng moton. Step (f): Statstcal propertes are calculated around the ponts. 4. Statstcs of Coherent Structures Fgs.5-7 show the mean azmuthal elocty profle of coherent structures n, GS and SGS felds for two flter functons. In these fgures, r represents the radus of the coherent structures, whch s determned by the dstance between the center and the locaton where the mean azmuthal elocty reaches the maxmum alue. Mean azmuthal elocty profles n all cases are normalzed by u rms and η, whch are obtaned from the feld. In Fgs.5-7, symbols represent an azmuthal elocty of a Burgers ortex, whch s wrtten as follows: Γ αr = 1 exp πr 4ν z αz, (9) =, (10) where Γ s the crculaton of the Burgers ortex tube and α s a stretchng parameter. Usng database, (10)(14) t was shown n the preous study that the mean azmuthal elocty profle of coherent fne scale structures n homogeneous sotropc turbulence and turbulent mxng layer could be approxmated by a Burgers ortex. They (10) Burgers' ortex () r Fg. 5 Mean azmuthal elocty profle of the coherent structures n the feld, normalzed by Kolmogoro mcroscale and r.m.s of elocty fluctuatons. Symbols represent elocty profle of a Burgers ortex and error bars denote arances of azmuthal elocty. also hae shown that the mean dameter of the coherent fne scale structures n homogeneous sotropc turbulence s about 10η and the maxmum of mean azmuthal elocty s about half of u rms, and these characterstcs of coherent structures are ndependent on Reynolds numbers of the flow. Fg.5 confrms these behaors of the coherent structures for ths Reynolds number case. Our nterest s to dscuss the characterstcs of GS and SGS coherent structures comparng wth the results, hence, we hae shown approxmaton of coherent structures n the GS and SGS felds wth that of the Burgers ortex n Fgs.6-7. The mean azmuthal elocty profle of GS-Gaussan coherent structures n Fg.6 shows a good agreement wth that of Burgers ortex n the whole range as t as n feld. In Fg.7, the agreement of mean azmuthal elocty profle for GS-sharp cutoff coherent structures wth Burgers ortex s good for relately small dstance (r<1η), that s, the Burgers ortex profle does not collapse wth mean azmuthal elocty profle for the large dstance. For Gaussan flter, the contrbuton of GS and SGS elocty felds are obtaned from the whole regon of elocty feld, whle sharp cutoff flter separates the elocty feld n GS and SGS felds at the cutoff wae number. Maybe, that s one reason of these dfferences n Fg. 6 and 7. Moreoer, flter sze =10η s relately small. For large flter wdth, the large-scale structures,.e., most of the coherent structures wth large dameter accumulate n the GS feld. Usng large flter wdth, we hae also seen (not shown) that the approxmaton of mean azmuthal elocty profle by Burgers ortex becomes qute good for relately large dameter tube lke coherent structures n GS sharp cutoff feld as well as n GS-gaussan fled. On the other hand, although not n hole range, but Fgs. 6 and 7 also show that the approxmaton of mean azmuthal elocty profle of SGS coherent structures n both SGS-Gaussan fled and SGS-sharp cutoff fled by Burgers, ortex s well n a certan range, and shows almost smlar behaor for both flter functons. These results suggest that mean azmuthal elocty profle of coherent structures n LES can be approxmate by that of a Burgers ortex as well as n. The comparsons of normalzed mean azmuthal elocty profle of coherent structures n, GS and SGS

6 Burgers' ortex (GS-Gaussan) GS-Gaussan r r 0.3 Burgers' ortex (SGS-Gaussan) SGS-Gaussan S r Fg. 6 Same as Fg.5, but GS feld and SGS feld obtaned by usng Gaussan flter. - - Burgers' ortex () r Burgers' ortex (S) r Fg. 7 Same as Fg.5, but GS feld and SGS feld obtaned by usng sharp cutoff flter. felds for dfferent flter functons are shown n Fg.8. Mean azmuthal elocty profle of GS structures for dfferent flter functons wth ths flter wdth, collapse wth profle at r=0η (Fg.8 ). On the other hand, mean azmuthal elocty profle of SGS structures for dfferent flter functons collapse each other at r=10η (Fg.8 ), but not wth profle. In all cases, the maxmum of mean azmuthal elocty of GS structures s hgher than structures, and of SGS structures s lower than structures. The maxmum of mean azmuthal elocty and dameter of coherent structures are about 0.6u rms and 15η n r Fg. 8 Comparson of mean azmuthal elocty profle ( ) of the, GS and SGS coherent structures. and GS felds, and SGS felds. In all cases, mean azmuthal elocty profle s normalzed by Kolmogoro mcroscale and r.m.s of elocty fluctuatons obtaned from feld. GS-Gaussan feld, and 0.8u rms and 0η n the GS-sharp cutoff feld. On the other hand, the maxmum of mean azmuthal elocty and dameter of coherent structures are about 0.u rms and 5η n SGS-Gaussan feld, and 03u rms and 7η n the SGS-sharp cutoff feld. The maxmum of mean azmuthal elocty for GS structures n the Gaussan feld s close to the feld, but the maxmum of mean azmuthal elocty for GS structures n the sharp cutoff feld s larger than the feld. In actual LES, GS structures can be dentfed but SGS structures need to be model by SGS model. Coherent structures hae hgh and thn energy dsspaton regons around them and these regons contrbute total energy dsspaton n turbulence. (4) Snce n SGS felds we can see the exstence of coherent structures somewhat smlar to feld, of course ery small n sze, ths result s ery worthy to deelop a structures base SGS model for LES. The probablty densty functons (pdf) of dameter of tube-lke coherent structures n, GS and SGS felds, whch are normalzed by η obtaned from feld are shown n Fg.9 for dfferent flter functons. It s clear that the peak of pdf of GS-structures n all cases do not concde wth each other or wth profle. Moreoer, the peak of pdf of GS structure n the GS-sharp cutoff feld shows hgher alue than n the GS-Gaussan feld for ths flter wdth. It s also reealed that the small dameter tube-lke coherent structure n GS feld s rare n all cases and large dameter of SGS structures reaches about 5η. Howeer, the collapse of pdf of dameter n the SGS felds s good, and the peak of pdf for SGS profle ncreases from the profle for all flter functons. Fg.10 shows the pdf of maxmum of mean azmuthal elocty of coherent structures n the same, GS and

7 p p GS-Gaussan D SGS-Gaussan S do not concde each other or wth profle. The maxmum alue of mean azmuthal elocty profle s lower than u rms,.e., about half of alue. Usng large and dfferent flter wdth, we hae seen (not shown) that the pdf of SGS profle becomes ery close wth profle and the maxmum alue n the SGS fled neer exceed the maxmum alue n the feld. That s, the characterstcs of GS and SGS coherent structures sgnfcantly depend on the flter wdths. Ths behaor s nterestng to study on the coherent structure n actual LES or to deelop some LES modelng based on the coherent structures n turbulence. In the preous study, (10)(14) t was shown that the coherent fne scale structures could be scaled by Kolmogoro mcroscale and r. m. s of elocty fluctuatons. The aboe results n the present study also suggest that the GS and SGS coherent structures may possble to scale by η and u rms as well as n D Fg. 9 Probablty densty functon of dameter of coherent structures normalzed by Kolmogoro mcroscale. and GS felds, and SGS felds. p p GS-Gaussan.0.5,max.0.5,max SGS-Gaussan S Fg. 10 Probablty densty functon of the maxmum of mean azmuthal elocty of coherent structures normalzed by r.m.s of elocty fluctuatons. and GS felds, and SGS felds. SGS felds usng dfferent flter functons, whch s normalzed by u rms. It s reealed that, GS-Gaussan profle collapse wth profle well n the whole range. The profle n the GS-sharp cutoff feld s smaller than profle at small azmuthal elocty and larger than that at large azmuthal elocty of the tube-lke coherent structures. The maxmum alue of mean azmuthal elocty n GS and felds seems to be same. On the other hand, SGS profles of mean azmuthal elocty for both flter functons 7 5 CONCLUSIONS In ths study, GS and SGS coherent structures n homogeneous sotropc turbulence are dentfed. The GS and SGS elocty felds are obtaned by flterng the elocty feld usng classcal flters for LES. By sualzng the contour surfaces of second narant Q n GS and SGS felds as well as n feld, t s shown that GS and SGS felds tself contan lots of dstnct tube-lke coherent structures n homogeneous sotropc turbulence, whch ndcates that the feld contan mult-scale structures n turbulent flow. The characterstcs of these GS and SGS coherent structures are somewhat smlar to coherent structures. These GS and SGS coherent structures can be scaled by η and u rms as well as n. 6 REFERENCES (1) Townsend, A. A., On the fne-scale structure of turbulence, Proc. R. Soc. Lond., A08(1951), pp () Tennekes, H., Smple model for the small-scale structure of turbulence, Phys. Fluds, 11(1968), No. 3, pp (3) Lundgren, T. S., Straned Spral Vortex Model for Turbulent Fne Structure, Phys. Fluds, 5(198), pp (4) Pulln, D. I. and Saffman, P. G., On the Lundgren- Townsend Model of Turbulent Fne Scales, Phys. Fluds, A5(1993), pp (5) Saffman, P. G. & Pulln, D. I., Ansotropy of the Lundgren-Townsend Model of Fne-scale Turbulence, Phys. Fluds, A6(1994), pp (6) She, Z. S., Jackson, E. and Orszag, S. A., Intermttent Vortex Structures n Homogeneous Isotropc Turbulence, Nature, 344(1990), pp (7) Vncent, A. and Meneguzz, M., The spatal structure and statstcal propertes of homogeneous turbulence, J. Flud Mech., 5(1991), pp (8) Jmenez, J., Wray, A. A., Saffman, P. G. and Rogallo, R. S., The Structure of Intense Vortcty n Isotropc Turbulence, J. Flud Mech., 55(1993), pp (9) Vncent, A. & Meneguzz, M., The dynamcs of

8 ortcty tubes n homogeneous turbulence, J. Flud Mech., 58(1994), pp (10) Tanahash, M., Myauch, T. and Ikeda, J., Scalng law of coherent fne scale structure n homogeneous sotropc turbulence, Proc. of 11 th Symp. on Turbulent Shear Flows. 1(1997), pp (11) Tanahash, M., Iwase, S., Uddn, M. A. and Myauch, T., Three dmensonal features of coherent fne scale eddes n turbulence, Turbulence and Shear Flow Phenomena-1 1 st Int. Symp., Santa Barbara, Calforna, Sep. 1-15(1999)., pp (1) Tanahash, M., Uddn, M. A., Iwase, S. & Myauch, T., Three dmensonal feature of coherent fne scale eddes n homogeneous sotropc turbulence, Trans. JSME, 65-B(1999), No. 638, pp (13) Uddn, M.A., Tanahash, M., Iwase, S. and Myauch, T., Vsualzaton of Axes of Coherent Fne Scale Eddes n Homogeneous Isotropc Turbulence, the 3 rd Pac. Symp. on Flow Vsualzaton and Image Processng (PSFVIP-3), Mau, Hawa, March.18-1, CD-ROM Proc., C3-(001), No. F304. (14) Tanahash, M., Myauch, T. and Matsuoka, K., Coherent fne scale structures n temporally deelopng turbulent mxng layers, Proc. nd Int. Symp. on Turbulence, Heat and Mass Transfer, Vol. (1997), pp (15) Tanahash, M., Das, S. K., Shoj, K. and Myauch, T., Coherent fne scale structure n turbulent channel flows. Trans. JSME, Vol. 65-B(1999), No. 638, pp (16) Tanahash, M., Tsujmoto, T, Karm, M. F., Fujmura, D. and Myauch, T., Ansotropy of MHD homogeneous turbulence. Trans. JSME, Vol. 65- B(1999), No. 640, pp (17) Leonard, A. Energy Cascade n Large-eddy Smulatons of Turbulent Flud Flows, Ad. n Geophys., 18A(1974), pp (18) Pomell, U., Yu, Y. and Adran, R. J., Subgrd-scale Energy Transfer and Near-wall Turbulence Structure, Phys. Fluds, 8(1) (1996), pp (19) Pomell, U., Mon, P. and Ferzger, J. H., Model Consstency n Large Eddy Smulaton of Turbulent Channel Flows, Phys. Fluds, 31(7)(1988), pp (0) Horut, K., A Proper Velocty Scale for Modelng Subgrd-scale Eddy Vscostes n Large Eddy Smulaton, Phys. Fluds A5(1)(199), pp (1) Uddn, M.A., Tanahash, M., Nobuyuk, T., Iwase, S. and Myauch, T., Grd-scale and Subgrd-scale Eddes n Turbulent Flows ~ a Pror Test, the 4 th Int. Conf. on Mech. Eng. (ICME001), Dhaka, Bangladesh, Dec. 6-8(001), n press. () Jeong, J. and Hussan, F., On the Identfcaton of a Vortex, J. Flud Mech., 85(1995), pp (3) Chong, M.S., Perry, A.E. and Cantwell, B.J., A General Classfcaton of Three-dmensonal Flow Felds, Phys. Fluds A(1990), pp (4) Tanahash, M., Myauch, T. & Yoshda, T., Characterstcs of small scale ortces related to turbulent energy dsspaton, Transport Phenomena n Thermal Flud Engneerng, Vol. (1996), pp

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

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

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

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

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

FFT Based Spectrum Analysis of Three Phase Signals in Park (d-q) Plane

FFT Based Spectrum Analysis of Three Phase Signals in Park (d-q) Plane Proceedngs of the 00 Internatonal Conference on Industral Engneerng and Operatons Management Dhaka, Bangladesh, January 9 0, 00 FFT Based Spectrum Analyss of Three Phase Sgnals n Park (d-q) Plane Anuradha

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

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

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

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

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

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

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

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

Snce h( q^; q) = hq ~ and h( p^ ; p) = hp, one can wrte ~ h hq hp = hq ~hp ~ (7) the uncertanty relaton for an arbtrary state. The states that mnmze t

Snce h( q^; q) = hq ~ and h( p^ ; p) = hp, one can wrte ~ h hq hp = hq ~hp ~ (7) the uncertanty relaton for an arbtrary state. The states that mnmze t 8.5: Many-body phenomena n condensed matter and atomc physcs Last moded: September, 003 Lecture. Squeezed States In ths lecture we shall contnue the dscusson of coherent states, focusng on ther propertes

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

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

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

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY INFORMATION do: 0.08/nature09 I. Resonant absorpton of XUV pulses n Kr + usng the reduced densty matrx approach The quantum beats nvestgated n ths paper are the result of nterference between two exctaton paths of

More information

A Comparative Investigation into Aerodynamic Performances of Two Set Finned Bodies with Circular and Non Circular Cross Sections

A Comparative Investigation into Aerodynamic Performances of Two Set Finned Bodies with Circular and Non Circular Cross Sections A Comparatve Investgaton nto Aerodynamc Performances of Two Set Fnned Bodes wth Crcular and Non Crcular Cross Sectons MAHMOUD MANI and SHADI MAHJOOB Aerospace Engneerng Department Amrkabr Unversty of Technology

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

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

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD Journal of Appled Mathematcs and Computatonal Mechancs 7, 6(3), 7- www.amcm.pcz.pl p-issn 99-9965 DOI:.75/jamcm.7.3. e-issn 353-588 THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS

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

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

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

Perfect Fluid Cosmological Model in the Frame Work Lyra s Manifold

Perfect Fluid Cosmological Model in the Frame Work Lyra s Manifold Prespacetme Journal December 06 Volume 7 Issue 6 pp. 095-099 Pund, A. M. & Avachar, G.., Perfect Flud Cosmologcal Model n the Frame Work Lyra s Manfold Perfect Flud Cosmologcal Model n the Frame Work Lyra

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

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

(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

Thermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850)

Thermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850) hermal-fluds I Chapter 18 ransent heat conducton Dr. Prmal Fernando prmal@eng.fsu.edu Ph: (850) 410-6323 1 ransent heat conducton In general, he temperature of a body vares wth tme as well as poston. In

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

The Synchronous 8th-Order Differential Attack on 12 Rounds of the Block Cipher HyRAL

The Synchronous 8th-Order Differential Attack on 12 Rounds of the Block Cipher HyRAL The Synchronous 8th-Order Dfferental Attack on 12 Rounds of the Block Cpher HyRAL Yasutaka Igarash, Sej Fukushma, and Tomohro Hachno Kagoshma Unversty, Kagoshma, Japan Emal: {garash, fukushma, hachno}@eee.kagoshma-u.ac.jp

More information

Basic concept of reactive flows. Basic concept of reactive flows Combustion Mixing and reaction in high viscous fluid Application of Chaos

Basic concept of reactive flows. Basic concept of reactive flows Combustion Mixing and reaction in high viscous fluid Application of Chaos Introducton to Toshhsa Ueda School of Scence for Open and Envronmental Systems Keo Unversty, Japan Combuston Mxng and reacton n hgh vscous flud Applcaton of Chaos Keo Unversty 1 Keo Unversty 2 What s reactve

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

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

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

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

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

Application of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems

Application of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems Mathematca Aeterna, Vol. 1, 011, no. 06, 405 415 Applcaton of B-Splne to Numercal Soluton of a System of Sngularly Perturbed Problems Yogesh Gupta Department of Mathematcs Unted College of Engneerng &

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 3: Element Properties Lecture 1: Natural Coordinates

Module 3: Element Properties Lecture 1: Natural Coordinates Module 3: Element Propertes Lecture : Natural Coordnates Natural coordnate system s bascally a local coordnate system whch allows the specfcaton of a pont wthn the element by a set of dmensonless numbers

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

Linear Momentum. Equation 1

Linear Momentum. Equation 1 Lnear Momentum OBJECTIVE Obsere collsons between two carts, testng or the conseraton o momentum. Measure energy changes durng derent types o collsons. Classy collsons as elastc, nelastc, or completely

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

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

REAL-TIME DETERMINATION OF INDOOR CONTAMINANT SOURCE LOCATION AND STRENGTH, PART II: WITH TWO SENSORS. Beijing , China,

REAL-TIME DETERMINATION OF INDOOR CONTAMINANT SOURCE LOCATION AND STRENGTH, PART II: WITH TWO SENSORS. Beijing , China, REAL-TIME DETERMIATIO OF IDOOR COTAMIAT SOURCE LOCATIO AD STREGTH, PART II: WITH TWO SESORS Hao Ca,, Xantng L, Wedng Long 3 Department of Buldng Scence, School of Archtecture, Tsnghua Unversty Bejng 84,

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

8.592J: Solutions for Assignment 7 Spring 2005

8.592J: Solutions for Assignment 7 Spring 2005 8.59J: Solutons for Assgnment 7 Sprng 5 Problem 1 (a) A flament of length l can be created by addton of a monomer to one of length l 1 (at rate a) or removal of a monomer from a flament of length l + 1

More information

1. Why turbulence occur? Hydrodynamic Instability. Hydrodynamic Instability. Centrifugal Instability: Rayleigh-Benard Instability:

1. Why turbulence occur? Hydrodynamic Instability. Hydrodynamic Instability. Centrifugal Instability: Rayleigh-Benard Instability: . Why turbulence occur? Hydrodynamc Instablty Hydrodynamc Instablty T Centrfugal Instablty: Ω Raylegh-Benard Instablty: Drvng force: centrfugal force Drvng force: buoyancy flud Dampng force: vscous dsspaton

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

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

CHAPTER 14 GENERAL PERTURBATION THEORY

CHAPTER 14 GENERAL PERTURBATION THEORY CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves

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

Numerical Simulation of Lid-Driven Cavity Flow Using the Lattice Boltzmann Method

Numerical Simulation of Lid-Driven Cavity Flow Using the Lattice Boltzmann Method Proceedngs of the 3th WSEAS Internatonal Conference on APPLIED MATHEMATICS (MATH'8) Numercal Smulaton of Ld-Drven Cavty Flow Usng the Lattce Boltzmann Method M.A. MUSSA, S. ABDULLAH *, C.S. NOR AZWADI

More information

Using Large Eddy Simulation to Study Airflows in and around Buildings

Using Large Eddy Simulation to Study Airflows in and around Buildings Jang, Y., Su., M., and Chen, Q. 003. Usng large eddy smulaton to study arflows n and around buldngs, ASHRAE Transactons, 109(), 517-56. Usng Large Eddy Smulaton to Study Arflows n and around Buldngs Y

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

FEATURES OF TURBULENT TRANSPORT OF MOMENTUM AND HEAT IN STABLY STRATIFIED BOUNDARY LAYERS AND THEIR REPRODUCTION IN ATMOSPHERIC MESOSCALE MODELS

FEATURES OF TURBULENT TRANSPORT OF MOMENTUM AND HEAT IN STABLY STRATIFIED BOUNDARY LAYERS AND THEIR REPRODUCTION IN ATMOSPHERIC MESOSCALE MODELS C I T E S 009_ Krasnoyarsk 009 FEATURES OF TURBULENT TRANSPORT OF MOMENTUM AND HEAT IN STABLY STRATIFIED BOUNDARY LAYERS AND THEIR REPRODUCTION IN ATMOSPHERIC MESOSCALE MODELS A. F. Kurbatsky Insttute

More information

Numerical Modelling and Experimental Validation of a Turbulent Separated Reattached Flow

Numerical Modelling and Experimental Validation of a Turbulent Separated Reattached Flow Numercal Modellng and Expermental Valdaton of a Turbulent Separated Reattached Flow Florn Popescu, Tănase Panat Abstract An expermental study was conducted to analyse the feld velocty of a fully developed

More information

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed (2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected

More information

2.29 Numerical Fluid Mechanics Fall 2011 Lecture 12

2.29 Numerical Fluid Mechanics Fall 2011 Lecture 12 REVIEW Lecture 11: 2.29 Numercal Flud Mechancs Fall 2011 Lecture 12 End of (Lnear) Algebrac Systems Gradent Methods Krylov Subspace Methods Precondtonng of Ax=b FINITE DIFFERENCES Classfcaton of Partal

More information

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve

More information

Power law and dimension of the maximum value for belief distribution with the max Deng entropy

Power law and dimension of the maximum value for belief distribution with the max Deng entropy Power law and dmenson of the maxmum value for belef dstrbuton wth the max Deng entropy Bngy Kang a, a College of Informaton Engneerng, Northwest A&F Unversty, Yanglng, Shaanx, 712100, Chna. Abstract Deng

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

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

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 3. Ax x i a i. i i

Lecture 3. Ax x i a i. i i 18.409 The Behavor of Algorthms n Practce 2/14/2 Lecturer: Dan Spelman Lecture 3 Scrbe: Arvnd Sankar 1 Largest sngular value In order to bound the condton number, we need an upper bound on the largest

More information

Some Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)

Some Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS) Some Comments on Acceleratng Convergence of Iteratve Sequences Usng Drect Inverson of the Iteratve Subspace (DIIS) C. Davd Sherrll School of Chemstry and Bochemstry Georga Insttute of Technology May 1998

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

Canonical transformations

Canonical transformations Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,

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

Energy Storage Elements: Capacitors and Inductors

Energy Storage Elements: Capacitors and Inductors CHAPTER 6 Energy Storage Elements: Capactors and Inductors To ths pont n our study of electronc crcuts, tme has not been mportant. The analyss and desgns we hae performed so far hae been statc, and all

More information

Energy configuration optimization of submerged propeller in oxidation ditch based on CFD

Energy configuration optimization of submerged propeller in oxidation ditch based on CFD IOP Conference Seres: Earth and Envronmental Scence Energy confguraton optmzaton of submerged propeller n oxdaton dtch based on CFD To cte ths artcle: S Y Wu et al 01 IOP Conf. Ser.: Earth Envron. Sc.

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

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

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

Notes on Frequency Estimation in Data Streams

Notes on Frequency Estimation in Data Streams Notes on Frequency Estmaton n Data Streams In (one of) the data streamng model(s), the data s a sequence of arrvals a 1, a 2,..., a m of the form a j = (, v) where s the dentty of the tem and belongs to

More information

NUMERICAL SIMULATION OF FLOW AROUND A NACA 0012 AIRFOIL IN TRANSIENT PITCHING MOTION USING IMMERSED BOUNDARY METHOD WITH VIRTUAL PHYSICAL MODEL

NUMERICAL SIMULATION OF FLOW AROUND A NACA 0012 AIRFOIL IN TRANSIENT PITCHING MOTION USING IMMERSED BOUNDARY METHOD WITH VIRTUAL PHYSICAL MODEL NUMERICAL SIMULATION OF FLOW AROUND A NACA 0012 AIRFOIL IN TRANSIENT PITCHING MOTION USING IMMERSED BOUNDARY METHOD WITH VIRTUAL PHYSICAL MODEL José Eduardo Santos Olvera Laboratory of Heat and Mass Transfer

More information

CME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 13

CME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 13 CME 30: NUMERICAL LINEAR ALGEBRA FALL 005/06 LECTURE 13 GENE H GOLUB 1 Iteratve Methods Very large problems (naturally sparse, from applcatons): teratve methods Structured matrces (even sometmes dense,

More information

Analysis of Unsteady Aerodynamics of a Car Model with Radiator in Dynamic Pitching Motion using LS-DYNA

Analysis of Unsteady Aerodynamics of a Car Model with Radiator in Dynamic Pitching Motion using LS-DYNA Analyss of Unsteady Aerodynamcs of a Car Model wth Radator n Dynamc Ptchng Moton usng LS-DYNA Yusuke Nakae 1, Jro Takamtsu 1, Hrosh Tanaka 1, Tsuyosh Yasuk 1 1 Toyota Motor Corporaton 1 Introducton Recently,

More information

Journal of Fluid Science and Technology

Journal of Fluid Science and Technology Journal of Flud Scence and Technology Numercal Smulaton of Incompressble Flows around a Fsh Model at Low Reynolds Number Usng Seamless Vrtual Boundary Method * Hdetosh NISHIDA ** and Kyohe TAJIRI ** **Department

More information

modeling of equilibrium and dynamic multi-component adsorption in a two-layered fixed bed for purification of hydrogen from methane reforming products

modeling of equilibrium and dynamic multi-component adsorption in a two-layered fixed bed for purification of hydrogen from methane reforming products modelng of equlbrum and dynamc mult-component adsorpton n a two-layered fxed bed for purfcaton of hydrogen from methane reformng products Mohammad A. Ebrahm, Mahmood R. G. Arsalan, Shohreh Fatem * Laboratory

More information

Direct Numerical Simulation of Particles-Bubbles Collisions Kernel in Homogeneous Isotropic Turbulence

Direct Numerical Simulation of Particles-Bubbles Collisions Kernel in Homogeneous Isotropic Turbulence Drect Numercal Smulaton of Partcles-Bubbles Collsons Kernel n Hassan E. Fayed and Saad A. Ragab* Department of Engneerng Scence and Mechancs Vrgna Tech, Blacksburg, VA Receved: January, Accepted: June

More information

STATISTICAL MECHANICAL ENSEMBLES 1 MICROSCOPIC AND MACROSCOPIC VARIABLES PHASE SPACE ENSEMBLES. CHE 524 A. Panagiotopoulos 1

STATISTICAL MECHANICAL ENSEMBLES 1 MICROSCOPIC AND MACROSCOPIC VARIABLES PHASE SPACE ENSEMBLES. CHE 524 A. Panagiotopoulos 1 CHE 54 A. Panagotopoulos STATSTCAL MECHACAL ESEMBLES MCROSCOPC AD MACROSCOPC ARABLES The central queston n Statstcal Mechancs can be phrased as follows: f partcles (atoms, molecules, electrons, nucle,

More information

University of Washington Department of Chemistry Chemistry 453 Winter Quarter 2015

University of Washington Department of Chemistry Chemistry 453 Winter Quarter 2015 Lecture 2. 1/07/15-1/09/15 Unversty of Washngton Department of Chemstry Chemstry 453 Wnter Quarter 2015 We are not talkng about truth. We are talkng about somethng that seems lke truth. The truth we want

More information

Circuit Variables. Unit: volt (V = J/C)

Circuit Variables. Unit: volt (V = J/C) Crcut Varables Scentfc nestgaton of statc electrcty was done n late 700 s and Coulomb s credted wth most of the dscoeres. He found that electrc charges hae two attrbutes: amount and polarty. There are

More information

Grid Generation around a Cylinder by Complex Potential Functions

Grid Generation around a Cylinder by Complex Potential Functions Research Journal of Appled Scences, Engneerng and Technolog 4(): 53-535, 0 ISSN: 040-7467 Mawell Scentfc Organzaton, 0 Submtted: December 0, 0 Accepted: Januar, 0 Publshed: June 0, 0 Grd Generaton around

More information

Chapter 2. Pythagorean Theorem. Right Hand Rule. Position. Distance Formula

Chapter 2. Pythagorean Theorem. Right Hand Rule. Position. Distance Formula Chapter Moton n One Dmenson Cartesan Coordnate System The most common coordnate system or representng postons n space s one based on three perpendcular spatal axes generally desgnated x, y, and z. Any

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

Thermodynamics and Kinetics of Solids 33. III. Statistical Thermodynamics. Â N i = N (5.3) N i. i =0. Â e i = E (5.4) has a maximum.

Thermodynamics and Kinetics of Solids 33. III. Statistical Thermodynamics. Â N i = N (5.3) N i. i =0. Â e i = E (5.4) has a maximum. hermodynamcs and Knetcs of Solds 33 III. Statstcal hermodynamcs 5. Statstcal reatment of hermodynamcs 5.1. Statstcs and Phenomenologcal hermodynamcs. Calculaton of the energetc state of each atomc or molecular

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

Problem Points Score Total 100

Problem Points Score Total 100 Physcs 450 Solutons of Sample Exam I Problem Ponts Score 1 8 15 3 17 4 0 5 0 Total 100 All wor must be shown n order to receve full credt. Wor must be legble and comprehensble wth answers clearly ndcated.

More information

A new turbulence model for Large Eddy Simulation

A new turbulence model for Large Eddy Simulation Adv. Studes Theor. Phys., Vol. 1, 2007, no. 6, 247-270 A new turbulence model for Large Eddy Smulaton Francesco Gallerano Dpartmento d Idraulca Trasport e Strade Unversta degl Stud d Roma La Sapenza Va

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

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

Turbulent Flow in Curved Square Duct: Prediction of Fluid flow and Heat transfer Characteristics

Turbulent Flow in Curved Square Duct: Prediction of Fluid flow and Heat transfer Characteristics Internatonal Research Journal of Engneerng and Technology (IRJET) e-issn: 2395-56 Volume: 4 Issue: 7 July -217 www.ret.net p-issn: 2395-72 Turbulent Flow n Curved Square Duct: Predcton of Flud flow and

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

Numerical analysis on the agglomeration behavior of fine particles in plane jet

Numerical analysis on the agglomeration behavior of fine particles in plane jet Numercal analyss on the agglomeraton behavor of fne partcles n plane et Mn Guo 1, a, Ja L 2, Xn Su 1 and Gurong Yang 1 1 Tann Academy of Envronmental Scences, Tann, 300191, Chna; 2 The Fourth Research

More information

Design optimization of a high specific speed Francis turbine runner

Design optimization of a high specific speed Francis turbine runner IOP Conference Seres: Earth and Envronmental Scence Desgn optmzaton of a hgh specfc speed Francs turbne runner To cte ths artcle: Y Enomoto et al 2012 IOP Conf. Ser.: Earth Envron. Sc. 15 032010 Vew the

More information

A strong wind prediction model for the railway operation based on the onsite measurement and the numerical simulation

A strong wind prediction model for the railway operation based on the onsite measurement and the numerical simulation A strong wnd predcton for the ralway operaton based on the onste measurement and the numercal smulaton Yayo Msu a, Atsush Yamaguch b, and Takesh Ishhara c a Dsaster Preventon Research Laboratory, East

More information