Simulation of Incompressible Flows in Two-Sided Lid-Driven Square Cavities. Part II - LBM
|
|
- Jocelyn Bryant
- 6 years ago
- Views:
Transcription
1 Perumal & Dass CFD Letters Vol. 2(1) Vol. 2 (1) March 2010 Smulaton of Incompressble Flows n Two-Sded Ld-Drven Square Cavtes. Part II - LBM D. Arumuga Perumal 1c and Anoop K. Dass 1 1 Department of Mechancal Engneerng, Indan Insttute of Technology Guwahat, Guwahat , INDIA Receved: 15/10/2009 Revsed 15/01/2010 Accepted 04/03/2010 Abstract The present paper computes the flow n a two-sded ld-drven square cavty by the Lattce Boltzmann Method (LBM). For some aspect ratos there exsts a multplcty of steady solutons, but the square cavty problem gves only a sngle steady soluton for both the parallel and antparallel moton of the walls. It s found that for parallel moton of the walls, there appears a par of counter-rotatng secondary vortces of equal sze near the centre of a wall. Because of symmetry, ths par of counter-rotatng vortces has smlar shapes and ther detaled study as to how they grow wth ncreasng Reynolds number has not yet been made by lattce Boltzmann Method. Such a study s attempted n ths paper through the LBM, as the results of the problem have the potental of beng used for testng varous soluton methods for ncompressble vscous flows. The results for the antparallel moton of the walls are also presented n some detal. To lend credblty to the LBM results they are further compared wth those obtaned from a fnte dfference method (FDM). Keywords: two-sded ld-drven square cavty; lattce Boltzmann method; D2Q9 model; bounce-back boundary condton; fnte dfference method. 1. Introducton Smulaton of ncompressble flows n two-sded ld-drven square cavty flow by FDM s dscussed n some detal n the frst part of our paper [1]. The present work focuses on the computaton and valdaton of two-sded ld-drven square cavty flows by the Lattce Boltzmann Method. The LBM s a relatvely novel technque that has become an alternatve to the tradtonal numercal methods for computng flud flow problems. The method s dscussed n suffcent detals n the books by Wolf-Gladrow [2] and Succ [3]. Chen and Doolen [4] have wrtten an excellent revew paper on the subject. Hstorcally, LBM orgnated from the method of Lattce Gas Cellular automata (LGCA), whch was frst ntroduced n 1973 by Hardy, Pomeau and de Pazzs (HPP) [5]. In 1986, Frsch, Hasslacher and Pomeau (FHP) obtaned the correct Naver-Stokes equatons usng a hexagonal lattce. Lattce Boltzmann Equatons has been used at the cradle of Lattce Gas Automata (LGA) by Frsch et al. [6] to calculate vscosty. To elmnate statstcal nose n 1988 c Correspondng Author: D. Arumuga Perumal Emal: d.perumal@tg.ernet.n Telephone: Fax: All rghts reserved. ISSR Journals 25
2 Smulaton of Incompressble Flows n Two-Sded Ld-Drven Square Cavtes. Part II - LBM McNamara and Zanett [7] dd away wth the Boolean operaton of LGA nvolvng the partcle occupaton varables by neglectng partcle correlatons and ntroducng averaged dstrbuton functons gvng rse to the LBM. Hguera and Jmenez [8] brought about an mportant smplfcaton n LBM by presentng an lattce Boltzmann Equaton (LBE) wth a lnearzed collson operator that assumes that the dstrbuton s close to the local equlbrum state. A partcularly smple verson of lnearzed collson operator based on the Bhatnagar-Gross-Krook (BGK) [9] collson model was ndependently ntroduced by several authors ncludng Koelman [10] and Chen et al. [11]. The lattce BGK (LBGK) model [12, 13] utlzes the local equlbrum dstrbuton functon to recover the macroscopc Naver-Stokes equatons. A revew of computatonal and expermental studes on ld-drven cavty flow can be found n Shankar & Deshpande [14]. They have studed and analyzed corner eddes, nonunqueness, transton and turbulence n the ld-drven cavty. The two-sded ld-drven rectangular cavty problem has been nvestgated n some detals by Kuhlmann and coworkers [15]. Many researchers [17-20] carred out smulatons of one-sded ld-drven cavty flow by the lattce Boltzmann method. Yong G La et al. [19] compared the lattce Boltzmann method and the fnte volume Naver-Stokes solver and concluded that bounce-back boundary condton has better than frst order accuracy. The present work uses Lattce Boltzmann BGK model (LBGK) wth sngle tme relaxaton and bounceback boundary condton to nvestgate the flow drven by parallel and antparallel moton of two facng walls n a square cavty for Reynolds number up to A nne-velocty ncompressble LB model n 2D space has been used n the present work snce t s known to gve more accurate results compared to seven-velocty ncompressble LB model. Ths paper s organzed n four sectons. In Secton 2 numercal methods ncludng LBGK wth sngle tme relaxaton scheme and two-dmensonal nne-velocty lattce model s descrbed. In Secton 3 the two-sded ld-drven cavty problem s descrbed and the results wth parallel and antparallel moton of the walls are presented and valdated. Concludng remarks are made n Secton Numercal Methods 2.1. Lattce Boltzmann Method As has already been mentoned the Lattce Boltzmann method (LBM) represents an alternatve possblty for the drect smulaton of the ncompressble flow. It s seen that the accuracy of the Lattce Boltzmann method s of second order both n space and tme [19]. The Lattce Boltzmann equaton whch can be lnked to the Boltzmann equaton n knetc theory s formulated as [17] f ( x c t, t t) f ( x, t) = (1) where f s the partcle dstrbuton functon, c s the partcle velocty along the th drecton and s the collson term. The so-called LBGK model wth sngle tme relaxaton, whch s a commonly used lattce Boltzmann method s gven by [17] 1 (, ) (, ) = - (, ) (0) f x c t t t f x t f x t f ( x, t) (2) where f (0) ( x, t) s the equlbrum dstrbuton functon at x,t and s the tme relaxaton parameter Two-Dmensonal Nne-Velocty Square Lattce Model The D2Q9 square lattce used here has nne dscrete veloctes. A square lattce s used, each node of whch has eght neghbours connected by eght lnks as shown n Fgure 1. Partcles resdng on a node move to ther nearest neghbours along these lnks n unt tme 26
3 Perumal & Dass CFD Letters Vol. 2(1) 2010 step. The occupaton of the rest partcle s desgnated as f 0. The occupaton of the partcles movng along the x- and y-axes are desgnated as dagonally movng partcles are desgnated as f, 5 f, 1 f, 6 f, 2 f, 7 f, 3 defned as c = 0, = 0 c = (cos( / 4( 1)),sn( / 4( 1))), = 1, 2, 3, 4 c = 2(cos( / 4( 1)),sn( / 4( 1))), = 5, 6, 7, 8. f, whle the occupaton of 4 f. The partcle veloctes are 8 (3) Fgure 1. Two-Dmensonal Nne-Velocty Square Lattce (D2Q9) Model. The macroscopc quanttes such as densty and momentum densty u are defned as velocty moments of the dstrbuton functon f as follows: N ρ= f, =0 (4) N ρ u= f =0 c (5) The densty (whch s drectly related to the pressure) s determned from the partcle dstrbuton functon. The densty and the veloctes satsfy the Naver-Stokes equatons n the low-mach number lmt. Ths can be demonstrated by usng the Chapman-Enskog expanson. In the nne-speed square lattce, a sutable equlbrum dstrbuton functon that has been proposed s [17] (0) 3 2 f = w 1 u, = 0 2 f (0) = 1 3( ) ( ) w c.u c.u u, = 1, 2, 3, 4 (6) f (0) = w 1 3( ) ( ) c.u c.u u, = 5, 6, 7, 8 where the lattce weghts are gven by w = 4/9, w = w = w = w = 1/9 and w = w = w = w = 1/36. The relaxaton tme s related to the vscosty by [19] 6 1 = (7) 2 where s the knematc vscosty measured n lattce unts. It s seen that = 0.5 s the crtcal value for ensurng a non-negatve knematc vscosty. Numercal nstablty can occur for a close to ths crtcal value. Ths stuaton takes place at hgh Reynolds numbers. In ths work Reynolds numbers up to 2000 n a lattce sze of have been nvestgated. 4 27
4 Smulaton of Incompressble Flows n Two-Sded Ld-Drven Square Cavtes. Part II - LBM Boundary Condtons In LBM several boundary condtons have been proposed [17-20]. The bounce-back scheme was used n these smulatons to copy the velocty no-slp condton on walls. In ths scheme, the partcle dstrbuton functon at the wall lattce node s assgned to be the partcle dstrbuton functon of ts opposte drecton. The basc argument for the use of ongrd bounce-back model' s that t s both mathematcally applcable and qute relevant for LBE smulatons of flud flows n smple bounded domans. For ths reason, ths boundary condton has been employed here on the two statonary walls. However for the movng walls, the equlbrum boundary condton s appled [17]. At the lattce nodes on the movng walls, flow-varables are re-set to ther pre-assumed values at the end of every streamngstep. A ld-velocty of U = 0.1 has been consdered n ths work. Snce Mach number s U/c s, where c s equals 1/ 3 a Mach number of , well wthn the ncompressble lmt, s obtaned Numercal Procedure The veloctes are assumed to be zero at the tme of startng the smulatons for all nodes. Intally, the equlbrum dstrbuton functon that corresponds to the flow-varables s assumed as the unknown dstrbuton functon for all node at t = 0. Also a unform flud densty =1.0 s mposed ntally. The boundary condtons for the parallel and antparallel wall moton are shown n Fgures 2 (a) and 2 (b). The soluton procedure of the above LBM at each tme step comprse the streamng and collson step, applcaton of boundary condtons, calculaton of partcle dstrbuton functon followed by calculaton of macroscopc varables. The LBE s solved n the soluton doman subjected to the above ntal and boundary condtons on a unform 2D square lattce structure. It s seen that the numercal algorthm of the lattce Boltzmann method s relatvely smpler compared wth conventonal Naver-Stokes methods that use technques lke fnte dfference, fnte volume or fnte element to dscretze the equatons at the macroscopc level. Another beneft of the present approach s the easness of programmng. Fgure 2. Two-Sded Ld-Drven Cavty for (a) parallel wall moton (b) antparallel wall moton wth LBM boundary condtons for the movng walls Fnte Dfference Results for Comparson As the LBM method s ntended to be used to compute the flow n a relatvely unexplored problem, Ref [1] that uses FDM to compute the flow n the same geometry serves as an effectve bass for comparson of most of the results presented n ths paper. The FDM code 28
5 Perumal & Dass CFD Letters Vol. 2(1) 2010 numercally solves the 2D Naver-Stokes equatons n the stream functon-vortcty form gven by = x y u + v = +. (9) t x y Re x y Ths form of the Naver-Stokes equaton s known to be amenable for hghly accurate numercal soluton. The code, whch s valdated aganst establshed results, uses ADI technque and second-order accurate central dfferencng for space dscretzaton. Consequently these results, whch are the only ones avalable for comparson, are hghly relable. Thus, favourable comparson of the present LBM results wth these accurate FDM results wll grant legtmacy to them. Ths n turn would lend added credblty to the results for ths relatvely unexplored problem, so that the flow confguraton could be used as a good test case for algorthm valdaton Valdaton of the LBM Code The developed LBM code s used to compute the sngle ld-drven square cavty flow for Re = 1000 on a lattce structure. Well-establshed results computed by Gha et al. [16] exst for the same problem whch s used for the present code-valdaton exercse. Fgures 3(a) and 3(b) shows the steady-state x-component of the velocty along the vertcal centrelne and the y-component of the velocty along the horzontal centrelne of the cavty at Re = Here the top ld moves from left to rght and t s observed that the agreement between our LBM results and those of Gha et al. [16] s excellent. The close agreement gves credblty to the result of our LBM code and t stands valdated. The same fgure also dsplays the FDM results [1] of the present authors, whch agan are n excellent agreement wth these results. (8) Fgure 3. Code valdaton: (a) u-velocty along vertcal centrelne and (b) v-velocty along horzontal centrelne for sngle ld-drven square-cavty (Re = 1000). 3. Two-Sded Ld-Drven Cavty Flow An ncompressble vscous flow n a square cavty whose top and bottom walls move n the same (parallel moton) or opposte (antparallel moton) drecton wth a unform velocty s the problem nvestgated n the present work. In the case of parallel wall moton, a free shear layer exsts mdway between the top and bottom walls apart from the wall bounded shear layers, whereas n the case of antparallel wall moton, only wall bounded shear layers exst. 29
6 Smulaton of Incompressble Flows n Two-Sded Ld-Drven Square Cavtes. Part II - LBM 3.1. Parallel Wall Moton In the parallel moton we consder, both the upper and lower plates move from left to rght n the x drecton wth the same velocty. Fgure 4 shows the streamlne patterns for varous Reynolds numbers on a lattce structure. Expectedly, the streamlnes are found to be symmetrcal wth respect to the horzontal centerlne. Fgure 4(a) shows the streamlne pattern for Re = 100. Two counter-rotatng prmary vortces symmetrcal to each other are seen to form wth a `free shear layer n between. At ths Reynolds number the prmary vortex cores are seen to be somewhat away from the centres of the top and the bottom halves of the cavty towards the rghthand top and rghthand bottom corners respectvely. Fgure 4(b) shows the streamlne pattern for Re = 400. At ths Reynolds number s seen also a par of counter-rotatng secondary vortces symmetrcally placed about the horzontal centrelne near the centre of the rght wall. Fgures 4(c) and 4(d) show the streamlne patterns for Re = 1000 and Re = 2000 respectvely. It s seen that wth the ncrease n Reynolds number the prmary vortex cores move towards the centres of the top and bottom halves of the cavty and the secondary vortex par grow n sze. At all the Reynolds numbers the counter-rotatng pars of prmary and secondary vortces mantan ther symmetry about the horzontal centerlne. (a) Re = 100 (b) Re = 400 (c) Re = 1000 (d) Re = 2000 Fgure 4. Streamlne pattern for parallel wall moton at (a) Re = 100 (b) Re = 400 (c) Re = 1000 and (d) Re = 2000 on a lattce. 30
7 Perumal & Dass CFD Letters Vol. 2(1) 2010 (a) Re = 100 (b) Re = 400 (c) Re = 1000 (d) Re = 2000 Fgure 5. Vortcty contours for parallel wall moton at (a) Re = 100 (b) Re = 400 (c) Re = 1000 and (d) Re = 2000 on a lattce. Fgure 6. A magnfed vew of secondary vortces for parallel wall moton at Re = 1000 (a) LBM, (b) FDM [1]. 31
8 Smulaton of Incompressble Flows n Two-Sded Ld-Drven Square Cavtes. Part II - LBM Fgure 5 shows the vortcty contours. In Fgure 6 a magnfed vew of the secondary vortces at Re = 1000 s compared wth those obtaned by FDM [1] to show that the vortces have been captured accurately. The other results are now authentcated by comparson wth the same reference. Fgures 7-10 show the comparson for horzontal velocty profles along vertcal lnes and vertcal velocty profles along horzontal lnes passng through dfferent ponts of the cavty for varous Reynolds numbers. Agreement of the velocty profles gven by LBM and FDM s once agan excellent. Table 1 gves the locatons of the vortces gven by the LBM and FDM for Re = 100, 400, 1000, 1500 and All these results show that the agreement s very good, whch further substantates the accuracy of the present LBM computatons. We thus accurately lst varous flow detals for ths relatvely unexplored confguraton. Fgure 7. Parallel wall moton, Re = 100: (a) horzontal velocty u along vertcal lnes (x=0.25, 0.50, 0.75), (b) vertcal velocty v along horzontal lnes (y=0.25, 0.50, 0.75). Fgure 8. Parallel wall moton, Re = 400: (a) horzontal velocty u along vertcal lnes (x=0.25, 0.50, 0.75), (b) vertcal velocty v along horzontal lnes (y=0.25, 0.50, 0.75). Fgure 9. Parallel wall moton, Re = 1000: (a) horzontal velocty u along vertcal lnes (x=0.25, 0.50, 0.75), (b) vertcal velocty v along horzontal lnes (y=0.25, 0.50, 0.75). 32
9 Perumal & Dass CFD Letters Vol. 2(1) 2010 Fgure 10. Parallel wall moton, Re = 2000: (a) horzontal velocty u along vertcal lnes (x=0.25, 0.50, 0.75), (b) vertcal velocty v along horzontal lnes (y=0.25, 0.50, 0.75). TABLE 1: LOCATIONS OF THE VORTICES FOR PARALLEL WALL MOTION: a. FDM [1], b. LBM. Prmary vortex centres Secondary vortex centres Re 100 Bottom Top Bottom Top x y x y x y x y a b a b a b a b a b Antparallel Wall Moton In the antparallel wall moton the upper and lower walls move n opposte drectons along the x-axs wth the same velocty. Fgure 11 shows the streamlne patterns on a lattce sze of for the antparallel wall moton. A sngle prmary vortex centred at the geometrc centre of the cavty s formed at low Reynolds numbers as shown n Fgure 11(a) for Re = 100 and Fgure 11(b) for Re = 400. The streamlne patterns for Re = 1000 and 2000 are shown n Fgures 11(c) and 11(d) respectvely. The ncreased Reynolds number results n the appearance of two secondary vortces near the top left and the bottom rght corners of the cavty and a very small shft of the prmary vortex centre from the geometrc centre of the cavty. It wll not be out of place to menton here that for the much-examned sngle ld-drven cavty flow the secondary vortex near the tralng edge of the movng wall does not appear at a Reynolds number as low as 1000 but much beyond that (at some value hgher than 2000). It s also seen that wth the ncrease n Reynolds number the shft of the prmary vortex centre from the geometrc centre of the cavty s small so that t remans very close to the geometrc centre even 33
10 Smulaton of Incompressble Flows n Two-Sded Ld-Drven Square Cavtes. Part II - LBM for these hgher values of Re = 1000 and However, n lne wth parallel wall moton, between Re = 1000 and 2000 the secondary vortces are seen to grow n sze. The vortcty contours for varous Reynolds numbers are shown n Fgure 12. A magnfed vew of the secondary vortces (the ones at the top left corner) at Re = 1000 gven by the present LBM and FDM [1] s shown n Fgure 13. The agreement s found to be very good. Fgures show the comparsons for horzontal velocty profles along vertcal lnes and vertcal velocty profles along horzontal lnes passng through dfferent ponts of the cavty for the varous Reynolds numbers and the agreement s excellent once agan. Table 2 gves the locatons of the vortces gven by the present LBM and FDM [1] for Re = 100, 400, 1000, 1500 and Clearly the results gven by the LBM agree very well wth those of the FDM. We thus see that the LBM results gven by all the fgures and tables are n excellent agreement wth those gven by the FDM [1]. Ths lends credblty to the current LBM results for ths problem. (a) Re = 100 (b) Re = 400 (c) Re = 1000 (d) Re = 2000 Fgure 11. Streamlne pattern for antparallel wall moton at (a) Re = 100 (b) Re = 400 (c) Re = 1000 and (d) Re = 2000 on a lattce. 34
11 Perumal & Dass CFD Letters Vol. 2(1) 2010 (a) Re = 100 (b) Re = 400 (c) Re = 1000 (d) Re = 2000 Fgure 12. Vortcty contours for antparallel wall moton at (a) Re = 100 (b) Re = 400 (c) Re = 1000 and (d) Re = 2000 on a lattce. Fgure 13. A magnfed vew of secondary vortces for antparallel wall moton at Re = 1000 (a) LBM, (b) FDM. 35
12 Smulaton of Incompressble Flows n Two-Sded Ld-Drven Square Cavtes. Part II - LBM Fgure 14. Antparallel wall moton, Re = 100: (a) horzontal velocty u along vertcal lnes (x=0.25, 0.50 and 0.75), (b) vertcal velocty v along horzontal lnes (y=0.25, 0.50 and 0.75). Fgure 15. Antparallel wall moton, Re = 400: (a) horzontal velocty u along vertcal lnes (x=0.25, 0.50 and 0.75), (b) vertcal velocty v along horzontal lnes (y=0.25, 0.50 and 0.75). Fgure 16. Antparallel wall moton, Re = 1000: (a) horzontal velocty u along vertcal lnes (x=0.25, 0.50 and 0.75), (b) vertcal velocty v along horzontal lnes (y=0.25, 0.50 and 0.75). Fgure 17. Antparallel wall moton, Re = 2000: (a) horzontal velocty u along vertcal lnes (x=0.25, 0.50 and 0.75), (b) vertcal velocty v along horzontal lnes (y=0.25, 0.50 and 0.75). 36
13 Perumal & Dass CFD Letters Vol. 2(1) 2010 TABLE 2: LOCATIONS OF THE VORTICES FOR ANTIPARALLEL WALL MOTION: (a) FDM [1], (b) LBM. Re Prmary Vortex (PV) x 100 a b a b a b a b a b y Secondary Vortces (SV) Bottom Rght Top Left x y x y Concluson In ths work a relatvely unexplored flow confguraton n a two-sded ld-drven square cavty s computed wth the LBM. The flow s nvestgated for both the parallel and antparallel moton of two of the facng walls. In the case of parallel wall moton, besdes two prmary vortces, there also appears a par of counter-rotatng secondary vortces symmetrcally placed about the centerlne parallel to the moton of the walls. About ths centerlne also appears a free shear layer wth the ncrease n Reynolds number. In the case of ant-parallel wall moton, besdes a sngle prmary vortex, there appears two secondary vortces near the tralng edges of the movng walls. Each of these near-tarlng-edge vortces are seen to appear at a much lower Reynolds number compared wth that n the sngle ld-drven cavty flow. The results of nvestgaton of ths relatvely unexplored problem comprse fgures and tables detalng streamlnes and vortcty contours, varous velocty profles and features of vortces lke ther centres and shapes. All these LBM results compare very well wth the only exstng (and accurate) set of FDM results. Ths not only lends credblty to, but also benchmarks these results for the present flow confguraton. Consequently these results, lke those of the sngle ld-drven cavty flow, may be used for valdatng the algorthms for computng steady flows governed by the 2D ncompressble Naver- Stokes equatons. References [1] Perumal, D.A, and A.K. Dass, Smulaton of Incompressble flows n two-sded ld-drven square cavtes. Part I FDM, CFD Letters, 2(1), pp [2] Wolf-Gladrow, D.A., Lattce-Gas Cellular Automata and Lattce Boltzmann Models: An Introducton. 2000: Sprnger-Verlag Berln-Hedelberg. [3] Succ, S., The Lattce Boltzmann Method for Flud Dynamcs and Beyond. 2000: Oxford Unversty Press. [4] Chen, S., and G.D. Doolen, Lattce Boltzmann method for flud flows. Annual Revew of Flud Mechancs, : p [5] Hardy, J., Y. Pomeau, and O. de Pazzs, Tme Evoluton of a Two-Dmensonal classcal lattce system. Physcal Revew Letters, : p [6] Frsch, U., B. Hasslacher, and Y. Pomeau, Lattce-Gas automata for the Naver-Stokes equatons. Physcal Revew Letters, : p [7] McNamara, G.R., and G. Zanett, Use of the Boltzmann Equaton to Smulate Lattce-Gas Automata. Physcal Revew Letters, : p
14 Smulaton of Incompressble Flows n Two-Sded Ld-Drven Square Cavtes. Part II - LBM [8] Hguera., and Jmenez, Boltzmann approach to Lattce Gas Smulatons. Europhyscs Letters, (7): p [9] Bhatnagar P.L., E.P. Gross and M. Krook, A Model for Collson Processes n Gases. I. Small Ampltude Processes n Charged and Neutral One-Component Systems. Physcal Revew, : p [10] Koelman J.M.V.A., A Smple Lattce Boltzmann Scheme for Naver-Stokes Flud Flow. Europhyscs Letters, (6): p [11] Chen S., H. Chen, Martnez and Matthaeus, Lattce Boltzmann Model for Smulaton of Magnetohydrodynamcs. Physcal Revew Letters, (27): p [12] Qan Y.H., D. d'humeres and P. Lallamand, Lattce BGK Models for Naver-Stokes Equaton. Europhyscs Letters, (6): p [13] Chen H., S. Chen, and Matthaeus, Recovery of the Naver-Stokes Equaton usng a lattcegas Boltzmann method. Physcal Revew A, (8): p. R [14] Shankar, P.N, and M.D. Deshpande, Flud Mechancs n the Drven Cavty. Annual Revew of Flud Mechancs, : p [15] Kuhlmann H.C., M. Wanschura, and H.J. Rath, Flow n two-sded ld-drven cavtes: nonunqueness, nstablty, and cellular structures. Journal of Computatonal Physcs, : p [16] Gha U., K.N. Gha, and C.T. Shn, Hgh-Re solutons for ncompressble flow usng Naver- Stokes equatons and a multgrd method. Journal of Computatonal Physcs, : p [17] Hou, S., Q. Zou, S. Chen, G. Doolen, and A.C. Cogley, Smulaton of Cavty Flow by the Lattce Boltzmann Method. Journal of Computatonal Physcs, : p [18] Perumal, D.A, and A.K. Dass, Lattce Boltzmann Smulaton of Two and Three Dmensonal Flow n a Ld-Drven Cavty. n 2008 Proceedngs of the Int. Conference on Advances n Mechancal Engg, IC-ICAME IISc Bangalore, INDIA. [19] La Y.G., C.L. Ln, and J. Huang, Accuracy and Effcency study of a Lattce Boltzmann method for steady-state flow smulatons. Numercal Heat Transfer, Part B, : p [20] Yu D., R. Me, L.S. Luo, and W. Shyy, Vscous flow computatons wth the method of lattce Boltzmann equaton. Progress n Aerospace Scences, : p
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 informationVarious Speed Ratios of Two-Sided Lid-Driven Cavity Flow using Lattice Boltzmann Method
Penerbt and Thermal Scences ISSN (onlne): 2289-7879 Varous Speed Ratos of Two-Sded Ld-Drven Cavty Flow usng Lattce Boltzmann Method Nor Azwad Che Sdk *,a and St Asyah Razal b Department of Thermo-Fluds,
More information(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 informationHYBRID LBM-FVM AND LBM-MCM METHODS FOR FLUID FLOW AND HEAT TRANSFER SIMULATION
HYBRID LBM-FVM AND LBM-MCM METHODS FOR FLUID FLOW AND HEAT TRANSFER SIMULATION Zheng L a,b, Mo Yang b and Yuwen Zhang a* a Department of Mechancal and Aerospace Engneerng, Unversty of Mssour, Columba,
More informationSIMULATION OF SOUND WAVE PROPAGATION IN TURBULENT FLOWS USING A LATTICE-BOLTZMANN SCHEME. Abstract
SIMULATION OF SOUND WAVE PROPAGATION IN TURBULENT FLOWS USING A LATTICE-BOLTZMANN SCHEME PACS REFERENCE: 43.20.Mv Andreas Wlde Fraunhofer Insttut für Integrerte Schaltungen, Außenstelle EAS Zeunerstr.
More informationLATTICE BOLTZMANN SIMULATION OF FLOW OVER A CIRCULAR CYLINDER AT MODERATE REYNOLDS NUMBERS
THERMAL SCIENCE: Year 014, Vol. 18, No. 4, pp. 135-146 135 LATTICE BOLTZMANN SIMULATION OF FLOW OVER A CIRCULAR CYLINDER AT MODERATE REYNOLDS NUMBERS by Dharmaraj ARUMUGA PERUMAL a*, Gundavarapu V. S.
More informationOn Higher Order Dynamics in Lattice-based Models Using Chapman-Enskog Method
NASA/CR-1999-209346 ICASE Report No. 99-21 On Hgher Order Dynamcs n Lattce-based Models Usng Chapman-Enskog Method Yue-Hong Qan Columba Unversty, New York, New York Ye Zhou IBM, Yorktown Heghts, New York
More informationResearch Article A Multilevel Finite Difference Scheme for One-Dimensional Burgers Equation Derived from the Lattice Boltzmann Method
Appled Mathematcs Volume 01, Artcle ID 9590, 13 pages do:10.1155/01/9590 Research Artcle A Multlevel Fnte Dfference Scheme for One-Dmensonal Burgers Equaton Derved from the Lattce Boltzmann Method Qaoe
More informationA Solution of the Harry-Dym Equation Using Lattice-Boltzmannn and a Solitary Wave Methods
Appled Mathematcal Scences, Vol. 11, 2017, no. 52, 2579-2586 HIKARI Ltd, www.m-hkar.com https://do.org/10.12988/ams.2017.79280 A Soluton of the Harry-Dym Equaton Usng Lattce-Boltzmannn and a Soltary Wave
More informationNumerical study of the effect of wedge angle of a wedge-shaped body in a channel using lattice Boltzmann method
Avalable onlne at www.scencedrect.com Proceda Engneerng 56 (2013 ) 232 238 5 th BSME Internatonal Conference on Thermal Engneerng Numercal study of the effect of wedge angle of a wedge-shaped body n a
More informationApplication of the lattice Boltzmann method for solving conduction problems with heat flux boundary condition.
November 5-7, 9 - Sousse Tunsa Applcaton of the lattce Boltzmann method for solvng conducton problems wth heat flux boundary condton. Raoudha CHAABANE, Faouz ASKRI, Sass Ben NASRALLAH Laboratore d Etudes
More informationWeek3, 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 informationOne-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 informationNumerical Study of Lid-Driven Square Cavity with Heat Generation Using LBM
Amercan Journal of Flud Dynamcs 23, 3(2): 4-47 DOI:.5923/j.ajfd.2332.4 Numercal Study of Ld-Drven Square Cavty wth Heat Generaton Usng LBM M. A. Tahe r, S. C. Saha 2, Y. W. Lee 3, H. D. Km,* Dept of Mechancal
More informationHigher Order Wall Boundary Conditions for Incompressible Flow Simulations
THE 5 TH ASIAN COMPUTAITIONAL FLUID DYNAMICS BUSAN KOREA OCTOBER 7-30 003 Hgher Order Wall Boundary Condtons for Incompressble Flow Smulatons Hdetosh Nshda. Department of Mechancal and System Engneerng
More informationA Solution of Porous Media Equation
Internatonal Mathematcal Forum, Vol. 11, 016, no. 15, 71-733 HIKARI Ltd, www.m-hkar.com http://dx.do.org/10.1988/mf.016.6669 A Soluton of Porous Meda Equaton F. Fonseca Unversdad Naconal de Colomba Grupo
More informationFinite Wings Steady, incompressible flow
Steady, ncompressble flow Geometrc propertes of a wng - Fnte thckness much smaller than the span and the chord - Defnton of wng geometry: a) Planform (varaton of chord and sweep angle) b) Secton/Arfol
More informationPhysics 207: Lecture 20. Today s Agenda Homework for Monday
Physcs 207: Lecture 20 Today s Agenda Homework for Monday Recap: Systems of Partcles Center of mass Velocty and acceleraton of the center of mass Dynamcs of the center of mass Lnear Momentum Example problems
More informationIntroduction to the lattice Boltzmann method
Introducton to LB Introducton to the lattce Boltzmann method Burkhard Dünweg Max Planck Insttute for Polymer Research Ackermannweg 10, D-55128 Manz, Germany duenweg@mpp-manz.mpg.de Introducton Naver-Stokes
More informationA large scale tsunami run-up simulation and numerical evaluation of fluid force during tsunami by using a particle method
A large scale tsunam run-up smulaton and numercal evaluaton of flud force durng tsunam by usng a partcle method *Mtsuteru Asa 1), Shoch Tanabe 2) and Masaharu Isshk 3) 1), 2) Department of Cvl Engneerng,
More informationNumerical 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 informationSTUDY 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 informationSimulated 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 informationSimulation for Pedestrian Dynamics by Real-Coded Cellular Automata (RCA)
Smulaton for Pedestran Dynamcs by Real-Coded Cellular Automata (RCA) Kazuhro Yamamoto 1*, Satosh Kokubo 1, Katsuhro Nshnar 2 1 Dep. Mechancal Scence and Engneerng, Nagoya Unversty, Japan * kazuhro@mech.nagoya-u.ac.jp
More informationAN INTRODUCTION TO ENTROPIC LATTICE BOLTZMANN SCHEME
AN INTRODUCTION TO ENTROPIC LATTICE BOLTZMANN SCHEME Francesca Tos 1 and Sauro Succ 2 1 Dp. d Matematca U.Dn Unverstá degl Stud d Frenze, Vale Morgagn, 67-a, 50134 Frenze, Italy tos@math.unf.t 2 Isttuto
More informationNON-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 informationCHAPTER 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 informationInvestigation 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 informationPHYS 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 informationA Hybrid Variational Iteration Method for Blasius Equation
Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method
More informationConsideration of 2D Unsteady Boundary Layer Over Oscillating Flat Plate
Proceedngs of the th WSEAS Internatonal Conference on Flud Mechancs and Aerodynamcs, Elounda, Greece, August -, (pp-) Consderaton of D Unsteady Boundary Layer Over Oscllatng Flat Plate N.M. NOURI, H.R.
More informationStatistical Energy Analysis for High Frequency Acoustic Analysis with LS-DYNA
14 th Internatonal Users Conference Sesson: ALE-FSI Statstcal Energy Analyss for Hgh Frequency Acoustc Analyss wth Zhe Cu 1, Yun Huang 1, Mhamed Soul 2, Tayeb Zeguar 3 1 Lvermore Software Technology Corporaton
More informationModule 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 informationA 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 informationTensor 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 informationFORCED CONVECTION HEAT TRANSFER FROM A RECTANGULAR CYLINDER: EFFECT OF ASPECT RATIO
ISTP-,, PRAGUE TH INTERNATIONAL SYMPOSIUM ON TRANSPORT PHENOMENA FORCED CONVECTION HEAT TRANSFER FROM A RECTANGULAR CYLINDER: EFFECT OF ASPECT RATIO Mohammad Rahnama*, Seyed-Mad Hasheman*, Mousa Farhad**
More informationIndeterminate 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 informationCHAPTER 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 informationElectrical 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 informationTurbulent 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 informationTools 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 informationPerfect 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 informationLattice Boltzmann Method in Theory and in Application to Coupled Problems
Lattce Boltzmann Method n Theory and n Applcaton to Coupled Problems Master Thess Danel Heubes supervsng Prof. Dr. Mchael Günther Dr. Andreas Bartel Unversty of Wuppertal Faculty of Mathematcs and Natural
More informationLattice Boltzmann Method and its Application to Flow Analysis in Porous Media
Specal Issue Multscale Smulatons for Materals 7 Research Report Lattce Boltzmann Method and ts Applcaton to Flow Analyss n Porous Meda Hdemtsu Hayash Abstract Under the exstence of an external force, a
More informationPhysics 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 informationChapter 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 informationSurvey of applications of discrete vortex method in civil engineering
Budownctwo Archtektura 5 (2009) 29-38 Survey of applcatons of dscrete vortex method n cvl engneerng Tomasz Nowck Lubln Unversty of Technology, Faculty of Cvl Engneerng and Archtecture, Department of Structural
More informationCOMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD
COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,
More informationcoordinates. Then, the position vectors are described by
Revewng, what we have dscussed so far: Generalzed coordnates Any number of varables (say, n) suffcent to specfy the confguraton of the system at each nstant to tme (need not be the mnmum number). In general,
More informationCHAPTER 6. LAGRANGE S EQUATIONS (Analytical Mechanics)
CHAPTER 6 LAGRANGE S EQUATIONS (Analytcal Mechancs) 1 Ex. 1: Consder a partcle movng on a fxed horzontal surface. r P Let, be the poston and F be the total force on the partcle. The FBD s: -mgk F 1 x O
More informationComputational 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 informationInductance 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 informationSimulation of Turbulent Flow Using FEM
Internatonal Journal of Engneerng and Technology Volume 2 No. 8, August, 2012 Smulaton of Turbulent Flow Usng FEM Sabah Tamm College of Computng, AlGhurar Unversty, Duba, Unted Arab Emrates. ABSTRACT An
More informationA Cartesian-grid integrated-rbf method for viscoelastic flows
Home Search Collectons Journals About Contact us My IOPscence A Cartesan-grd ntegrated-rbf method for vscoelastc flows Ths artcle has been downloaded from IOPscence. Please scroll down to see the full
More informationSimulation of Gas Exothermic Chemical Reaction in Porous Media Reactor with Lattice Boltzmann Method
Journal of Thermal Scence Vol.22, No.1 (2013) 42 47 DOI: 10.1007/s11630-013-0590-5 Artcle ID: 1003-2169(2013)01-0042-06 Smulaton of Gas Exothermc Chemcal Reacton n Porous Meda Reactor wth Lattce Boltzmann
More informationCOMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
More informationA new Approach for Solving Linear Ordinary Differential Equations
, ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of
More informationJournal 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 information3D Lattice Boltzmann Simulation of Droplet Formation in a Cross-Junction Microchannel
3D Lattce Boltzmann Smulaton of Droplet Formaton n a Cross-Juncton Mcrochannel SURESH ALAPATI, SANGMO KANG AND YONG KWEON SUH * Department of Mechancal Engneerng Dong-A Unversty 840 Hadan-dong, Saha-gu,
More informationPhysics 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 informationModule 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 informationWeek 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 informationPhysics 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 informationSimulation of Fluid Flow Inside a Back-ward-facing Step by MRT- LBM
2012 Internatonal Conference on Flud Dynacs and Therodynacs Technologes (FDTT 2012) IPCSIT vol.33(2012) (2012) IACSIT Press, Sngapore Sulaton of Flud Flow Insde a Back-ward-facng Step by MRT- LBM Mohaad
More informationProf. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model
EXACT OE-DIMESIOAL ISIG MODEL The one-dmensonal Isng model conssts of a chan of spns, each spn nteractng only wth ts two nearest neghbors. The smple Isng problem n one dmenson can be solved drectly n several
More informationHigh 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 informationAmplification 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 informationPhysics 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 informationSupplemental 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 informationThe 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 informationIntroduction to Computational Fluid Dynamics
Introducton to Computatonal Flud Dynamcs M. Zanub 1, T. Mahalakshm 2 1 (PG MATHS), Department of Mathematcs, St. Josephs College of Arts and Scence for Women-Hosur, Peryar Unversty 2 Assstance professor,
More informationLattice Boltzmann simulation of turbulent natural convection in tall enclosures.
Lattce Boltzmann smulaton of turbulent natural convecton n tall enclosures Hasan SAJJADI 1*, Reza KEFAYATI 1* Department of Engneerng, Shahd Bahonar Unversty of Kerman, Kerman, I. R. Iran School of Computer
More informationOn the Interval Zoro Symmetric Single-step Procedure for Simultaneous Finding of Polynomial Zeros
Appled Mathematcal Scences, Vol. 5, 2011, no. 75, 3693-3706 On the Interval Zoro Symmetrc Sngle-step Procedure for Smultaneous Fndng of Polynomal Zeros S. F. M. Rusl, M. Mons, M. A. Hassan and W. J. Leong
More informationWorkshop: 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 informationMagnetic Diffusion using Lattice-Boltzmann
Revsta Mexcana de Físca S 58 2 188 14 DICIEMBRE 212 Magnetc Dffuson usng Lattce-Boltzmann F. Fonseca Physcs Department, Unversdad Naconal de Colomba, Bogotá Colomba. e-mal: frfonsecaf@unal.edu.co Recbdo
More informationIC 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 informationThe equation of motion of a dynamical system is given by a set of differential equations. That is (1)
Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence
More informationTurbulent 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 informationKernel 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 informationFirst 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 informationMathematical 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 informationLecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES
COMPUTATIONAL FLUID DYNAMICS: FDM: Appromaton of Second Order Dervatves Lecture APPROXIMATION OF SECOMD ORDER DERIVATIVES. APPROXIMATION OF SECOND ORDER DERIVATIVES Second order dervatves appear n dffusve
More informationAppendix 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 informationSTATISTICAL MECHANICS
STATISTICAL MECHANICS Thermal Energy Recall that KE can always be separated nto 2 terms: KE system = 1 2 M 2 total v CM KE nternal Rgd-body rotaton and elastc / sound waves Use smplfyng assumptons KE of
More informationDurban Watson for Testing the Lack-of-Fit of Polynomial Regression Models without Replications
Durban Watson for Testng the Lack-of-Ft of Polynomal Regresson Models wthout Replcatons Ruba A. Alyaf, Maha A. Omar, Abdullah A. Al-Shha ralyaf@ksu.edu.sa, maomar@ksu.edu.sa, aalshha@ksu.edu.sa Department
More informationGlobal 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 informationIrregular 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 informationThe Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method
Journal of Electromagnetc Analyss and Applcatons, 04, 6, 0-08 Publshed Onlne September 04 n ScRes. http://www.scrp.org/journal/jemaa http://dx.do.org/0.46/jemaa.04.6000 The Exact Formulaton of the Inverse
More informationPerformance analysis of a suction muffler in a hermetic reciprocating compressor using CAA techniques based on Lattice Boltzmann Method
Performance analyss of a sucton muffler n a hermetc recprocatng compressor usng CAA technques based on Lattce Boltzmann Method Songjune LEE 1 ; Cheolung CHEONG 2 Pusan Natonal Unversty, Republc of Korea
More informationThe Solution of the Two-Dimensional Gross-Pitaevskii Equation Using Lattice-Boltzmann and He s Semi-Inverse Method
Internatonal Journal of Mathematcal Analyss Vol., 7, no., 69-77 HIKARI Ltd, www.m-hkar.com https://do.org/.988/jma.7.634 The Soluton of the Two-Dmensonal Gross-Ptaevsk Equaton Usng Lattce-Boltzmann and
More informationAvailable online at ScienceDirect. Procedia IUTAM 15 (2015 ) 34 40
Avalable onlne at www.scencedrect.com ScenceDrect Proceda IUTAM 5 (5 ) 34 4 IUTAM Symposum on Multphase flows wth phase change: challenges and opportuntes, Hyderabad, Inda (December 8 December, 4) Turbulent
More informationSupplementary 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 informationAerodynamic analysis involving moving parts with XFlow
Aerodynamc analyss nvolvng movng parts wth XFlow Ths report presents some of the possbltes of XFlow 2011, specfcally those related to the aerodynamc analyss of movng parts. Usng tradtonal Computatonal
More informationNUMERICAL MODEL FOR NON-DARCY FLOW THROUGH COARSE POROUS MEDIA USING THE MOVING PARTICLE SIMULATION METHOD
THERMAL SCIENCE: Year 2018, Vol. 22, No. 5, pp. 1955-1962 1955 NUMERICAL MODEL FOR NON-DARCY FLOW THROUGH COARSE POROUS MEDIA USING THE MOVING PARTICLE SIMULATION METHOD Introducton by Tomok IZUMI a* and
More informationTHE 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 informationarxiv: v1 [physics.flu-dyn] 16 Sep 2013
Three-Dmensonal Smoothed Partcle Hydrodynamcs Method for Smulatng Free Surface Flows Rzal Dw Prayogo a,b, Chrstan Fredy Naa a a Faculty of Mathematcs and Natural Scences, Insttut Teknolog Bandung, Jl.
More informationarxiv: v1 [physics.comp-ph] 20 Nov 2018
arxv:8.08530v [physcs.comp-ph] 20 Nov 208 Internatonal Journal of Modern Physcs C c World Scentfc Publshng Company Valdty of the Molecular-Dynamcs-Lattce-Gas Global Equlbrum Dstrbuton Functon M. Reza Parsa
More informationChapter 02: Numerical methods for microfluidics. Xiangyu Hu Technical University of Munich
Chapter 02: Numercal methods for mcrofludcs Xangyu Hu Techncal Unversty of Munch Possble numercal approaches Macroscopc approaches Fnte volume/element method Thn flm method Mcroscopc approaches Molecular
More informationNew Method for Solving Poisson Equation. on Irregular Domains
Appled Mathematcal Scences Vol. 6 01 no. 8 369 380 New Method for Solvng Posson Equaton on Irregular Domans J. Izadan and N. Karamooz Department of Mathematcs Facult of Scences Mashhad BranchIslamc Azad
More informationAP Physics 1 & 2 Summer Assignment
AP Physcs 1 & 2 Summer Assgnment AP Physcs 1 requres an exceptonal profcency n algebra, trgonometry, and geometry. It was desgned by a select group of college professors and hgh school scence teachers
More information