Multiple-relaxation-time lattice Boltzmann modeling of. incompressible flows in porous media
|
|
- Alban Parrish
- 5 years ago
- Views:
Transcription
1 Muliple-relaaion-ime laice Bolzmann modeling of incompressible flows in porous media Qing Liu, Ya-Ling He* Ke Laboraor of Thermo-Fluid Science and Engineering of MOE, School of Energ and Power Engineering, Xi an Jiaoong Universi, Xi an, Shaani, 749, P.R.China (*Corresponding auhor. Absrac In his paper, a wo-dimensional eigh-veloci (DQ8) muliple-relaaion-ime (MRT) laice Bolzmann (LB) model is proposed for incompressible porous flows a he represenaive elemenar volume scale based on he Brinkman-Forchheimer-eended Darc formulaion. In he model, he porosi is included ino he pressure-based equilibrium momens, and he linear and nonlinear drag forces of he porous media are incorporaed ino he model b adding a forcing erm o he MRT-LB equaion in he momen space. Through he Chapman-Enskog analsis, he generalized Navier-Sokes equaions can be recovered eacl wihou arificial compressible errors. Numerical simulaions of several pical wo-dimensional porous flows are carried ou o validae he presen MRT-LB model. The numerical resuls of he presen MRT-LB model are in good agreemen wih he analical soluions and/or oher numerical soluions repored in he lieraure. Kewords: Laice Bolzmann model; Muliple-relaaion-ime; Porous media; incompressible flows; REV scale.. Inroducion
2 Fluid flow in porous media has gained significan research ineres due o is imporance of relaed scienific and indusrial applicaions, which include conaminan ranspor in groundwaer, crude oil eploraion and eracion, radioacive wase managemen, hdrogeolog and so on. Comprehensive lieraure surves of hese applicaions can be found in Refs. [-3]. For incompressible flows in porous media a he represenaive elemenar volume (REV) scale, he Darc model, he Brinkman-eended Darc model and he Forchheimer-eended Darc model have been widel emploed in man sudies. However, he Darc model and he wo eended (Brinkman and Forchheimer) models have some inrinsic limiaions in simulaing porous flows [4]. In order o overcome he shorcomings of he above menioned models, he Brinkman-Forchheimer-eended Darc model (also called he generalized model) has been developed b several research groups [4-7]. In he generalized model, he viscous and inerial forces are incorporaed ino he momenum equaion b using he local volume-averaging echnique. The Darc model and he wo eended models can be regarded as he limiing cases of he generalized model. Moreover, he single phase flow and he ransien flow in porous media can be solved b he generalized model [7]. As poined ou b Vafai and Kim [8], numerical resuls based on he Brinkman-Forchheimer-eended Darc formulaion have been shown o be in good agreemen wih he eperimenal predicions, and a porous medium/free-fluid inerface can be bes deal wih b he Brinkman-Forchheimer-eended Darc formulaion and he coninui of sresses and velociies a he inerface. In he pas several decades, various radiional numerical mehods, such as he finie volume (FV) mehod, he finie difference (FD) mehod, and he finie elemen (FE) mehod, have been emploed o sud porous flows based on he generalized model. The laice Bolzmann (LB) mehod, as a mesoscopic numerical scheme originaes from he laice-gas auomaa (LGA) mehod [9], has achieved significan success in modeling comple fluid
3 flows and relaed ranspor phenomena due o is kineic characerisics [-6]. Owing o he kineic background, he LB mehod has some disincive advanages over he radiional numerical mehods (e.g., see Ref. [7]). Recenl, he LB mehod has been successfull applied o simulae fluid flows in porous media a he REV scale [8-5]. In he REV scale mehod, he deailed geomeric srucure of he media is ignored, and he sandard LB equaion is modified b adding an addiional erm o accoun for he presence of he porous media. Therefore, LB mehod a he REV scale can be used for ssems wih large compuaional domain. I is worh menioning ha he REV LB mehod is ver effecive for simulaing fluid flows in he region which is pariall filled wih a porous medium. As repored in Ref. [], he disconinui of he veloci-gradien a he porous medium/free-fluid inerface can be well capured b he LB mehod wihou including he sress boundar condiion ino he simulaion model. However, o he bes of our knowledge, mos of he eising REV LB models [8-5] for incompressible porous flows emplo he Bhanagar-Gross-Krook (BGK) model (also called he singe-relaaion-ime model) [6] o represen he collision process. Alhough he BGK model has become he mos popular one because of is simplici, here are several well-known criicisms on his model, such as he numerical insabili a low values of viscosi [7-9] and he inaccurac in reaing boundar condiions [3]. On he oher hand, i has been demonsraed ha he deficiencies of he BGK model can be addressed b emploing he muliple-relaaion-ime (MRT) model [3]. Hence, he aim of his paper is o develop a new MRT-LB model for incompressible porous flows a he REV scale based on he generalized model. In he model, a pressure-based MRT-LB equaion wih he eigh-b-eigh collision mari [3] is consruced in he framework of he sandard MRT-LB mehod. The remainder of his paper is organized as follows. In Secion, he MRT-LB model for incompressible porous flows a he REV scale is presened. In Secion 3, numerical ess of he
4 MRT-LB model are performed for he porous Poiseuille flow, porous Couee flow, lid-driven flow in a square porous cavi, and naural convecion flow in a square porous cavi. Finall, a brief conclusion is made in Secion 4.. MRT-LB model for incompressible flows in porous media. Macroscopic governing equaions The fluid flow is assumed o be wo-dimensional, laminar and incompressible. For isohermal incompressible porous flows a he REV scale, he generalized model proposed b Nihiarasu e al. [4, 7] is emploed in he presen sud. The dimensional governing equaions of he generalized model can be wrien as u, () u u F, () u u p e where is he average fluid densi, u and p are he volume-averaged fluid veloci and pressure, respecivel, is he porosi, and e is he effecive kineic viscosi. F = F, F denoes he oal bod force induced b he porous mari and oher eernal forces, which can be epressed as [6, ] F F u u u a, (3) K K where is he kineic viscosi of he fluid, K is he permeabili, F is he geomeric funcion, a is he bod force due o an eernal force, and u u u, in which u and u are he - and -componens of he macroscopic veloci u, respecivel. Based on Ergun s relaion [33], he geomeric funcion F and he permeabili K of he porous media can be epressed as [34] F , d p K 3 5, (4)
5 where d p is he solid paricle diameer. The flow governed b Eqs. () and () are characerized b he porosi and several dimensionless parameers: he Darc number Da, he viscosi raio J, and he Renolds number Re, which are defined as K e Da, J L, Re LU, (5) where L is he characerisic lengh, and U is he characerisic veloci.. MRT-LB model In his subsecion, a wo-dimensional MRT-LB model wih eigh velociies (DQ8 model) is presened o sud incompressible porous flows. According o Refs. [35, 36], he MRT-LB equaion wih an eplici reamen of he forcing erm can be wrien as eq,,,, Λ f e f M Λ m m M I S, (6) where M is he 8 8 orhogonal ransformaion mari, Λ is he non-negaive 8 8 diagonal relaaion mari, and I is he 8 8 uni mari. The boldface smbols, f, m, m eq, and S represen 8-dimensional (column) vecors:, f,, f,,, f8, f, f,,, f8 8, f e e e, m,, m,,, m8, m where T is he ranspose operaor, f, eq eq eq eq 8 T, T, m, m,, m,,, m, i S = S, S,, S 8 T, is he discree disribuion funcion corresponding o he T, T, discree veloci e and ime,, i m and eq, m are he veloci momens of he discree disribuion funcions f and he corresponding equilibrium momens, respecivel, and S i,,, 8 are componens of he forcing erm S. i
6 In he DQ8 model, he eigh discree velociies e i,,, 8 (see Fig. ) are defined as [3, i 37] where c, in which and cos i,sin i c, i ~ 4 e i, (7) cos i 9 4,sin i 9 4 c, i 5 ~ 8 are he ime sep and laice spacing, respecivel. The sound speed of he DQ8 model is c c 3. In he presen work, c is se o be, which leads o s. The ransformaion mari linearl ransforms he discree disribuion funcions f 8 (veloci space) o heir veloci momens 8 m (momen space): m Mf, f = M m. (8) For he DQ8 model, he eigh veloci momens m i,,, 8 corresponding o he i discree velociies are: where m m p is he pressure, m momenum j, j m, m, m3, m4, m5, m6, m7, m8 p, e, j, q, j, q, p, p T T, (9) e is relaed o energ, m3,5 j, are componens of he J, m4,6 q, are relaed o energ flu, and m7,8 p, are relaed o he diagonal and off-diagonal componens of he sress ensor. Wih he ordering of he veloci momens i m given above, he ransformaion mari M can be easil consruced ( c ) [3]: Μ. () The equilibrium momens eq m for he veloci momens m are given as follows:
7 m T eq 5 u eq eq eq eq eq p, e, j, q, j, q, p, p 3 3, () where eq 3 u u e p, 4 eq q c u, eq q cu, eq 3 u u p, eq 3 uu p 3. () To ge he correc generalized Navier-Sokes equaions () and (), he parameers are chosen as follows:, c c,, and 3 3. In he above equilibrium momens, we have emploed he incompressibili approimaion, i.e., he fluid densi and j, j J u. The relaaion mari Λ is given b Λ =diag s, s, s, s, s, s, s, s =diag, se,, sq,, sq, s, s. (3) The evoluion process of he MRT-LB equaion (6) consiss of wo seps: he collision sep and sreaming sep [7]. Usuall, he collision sep is implemened in he momen space: eq m m m m I S, (4) where m represen he veloci momens afer collision, while he sreaming sep is sill carried ou in he veloci space:,, f e f, (5) i i i where f M m.
8 In order o derive he correc equaions of hdrodnamics, he forcing erm S in he momen space should be chosen appropriael. For he DQ8 MRT-LB model, he forcing erm S can be chosen as 4u F 5 F F S F, (6) F uf uf uf u F where F, F F is given b Eq. (3). The fluid veloci u is defined as u e F, (7) 8 i fi i Noe ha he oal bod force F also conains he veloci u, so Eq. (7) is a nonlinear equaion for u. According o Eqs (3) and (7), he veloci u can be calculaed eplicil b v u l l l v, (8) where v, l and l are given b 8 v ei fi a, l i K, F l K. (9) According o Ref. [39], he fluid pressure p can be deermined b where 4 9 c p f s 8 s i i u, () and s u u 3. Here, u 3e u 4.5 e u.5 u wih 4 9, ~4 9 and 5~8 36. si i i i The effecive kineic viscosi e in he model is defined as
9 e s c () wih s7 s8 s. Through he Chapman-Enskog analsis [3, 35, 38] of he MRT-LB equaion (6) in he momen space, he generalized Navier-Sokes equaions () and () can be recovered eacl in he incompressible limi (see Appendi A for deails). If he eigh relaaion raes si i,,,8 are se o be a single value, i.e., Λ I, hen he presen MRT-LB model reduces o a BGK-LB model wih he following equilibrium disribuion funcion: f p e u e u u eq i i i i 4 cs cs cs cs. () I should be noed ha, as and Da, he presen MRT-LB model reduces o he incompressible MRT-LB model [3] for free-fluid flows wihou porous media. When F, he Brinkman-eended Darc equaion can be obained from he presen MRT-LB model. 3. Numerical ess In his secion, numerical simulaions of several pical wo-dimensional porous flows are carried ou o validae he proposed MRT-LB model. The esing problems include he Poiseuille flow in a channel filled wih porous media, he Couee flow beween wo parallel plaes filled wih porous media, he lid-driven flow in a square porous cavi, and he naural convecion flow in a square porous cavi. The presen numerical resuls are compared wih he analical and numerical soluions in previous sudies. In simulaions, we se, c =, and. The relaaion raes si i 8 are chosen as follows: s s3 s 5, s., s4 s 6., and s7 s8 s. Unless oherwise saed, he non-equilibrium erapolaion scheme [4] is emploed o rea he veloci and emperaure boundar condiions in simulaions. 3. Poiseuille flow in a channel filled wih porous media
10 The firs es problem is he Poiseuille flow in a D channel filled wih porous media. The heigh of he channel is L, and he flow is driven b an eernal force =, a along he channel a direcion. When he flow is full developed in he -direcion (along he channel direcion), he governing equaion of he flow can be epressed as wih u u L u e K F u u a K,,, and he -direcion veloci componen u is zero everwhere. The (3) Renolds number Re is defined as Re Lu, where u is he maimum veloci of he porous Poiseuille flow (wihou he nonlinear drag force, i.e., F ) along he cenerline of he channel []: where K. e ak L u cosh, (4) In simulaions, he porosi is se o be., he viscosi raio J is se o be, he Darc number Da changes from 6 o 3, and he Renolds number Re changes from. o. The relaaion rae s is se o be 5 3 (.6 ) wih a N N 8 8 square mesh. The nonequilibrium erapolaion scheme is emploed o he op and boom plaes for no-slip veloci boundar condiion, and he periodic boundar condiions are emploed o he inle and oule of he channel. A, he veloci momens are se o be heir equilibrium values, and he veloci field is iniialized o be wih a consan pressure p. The veloci profiles prediced b he presen MRT-LB model for differen Re and Da are shown in Fig.. The numerical soluions given b Guo and Zhao [] using he FD mehod are also included in he figure for comparison. From Fig. i is clearl seen ha he presen numerical resuls agree well wih he FD soluions repored in he lieraure. 3. Couee flow beween wo parallel plaes filled wih porous media
11 We now appl he MRT-LB model o he Couee flow in a D channel filled wih fluid-sauraed porous media. The flow is driven b he op plae moving along he -direcion wih a uniform veloci u, while he boom plae is fied. The Renolds number Re of his flow is defined b Re Lu. A he sead sae, he flow sill follows Eq. (4) ( a ) wih u, and u, L u. In simulaions, he porosi is se o be., he viscosi raio J is se o be, and he boundar and iniial condiions are he same as hose in he porous Poiseuille flow. For Re =.,, 5, and, he kineic viscosi is se o be.4,.8,.3, and.8 based on a N N 8 8 square mesh, respecivel. In Fig. 3, he veloci profiles prediced b he presen MRT-LB model are presened for various Re and Da. The FD soluions in Ref. [] are also included in Fig. 3 for comparison. Ecellen agreemen can be found beween he presen MRT-LB resuls and he FD soluions. The MRT-LB model is also applied o simulae he modified Couee flow in a D channel pariall filled wih a fluid-sauraed porous medium, which has been invesigaed in Refs. [,, 4]. A porous laer is posiioned in he lower par of he channel such ha here is a gap beween he op plae and he medium. The porosi is se o be. for L (porous region) and for L (free-fluid region). For small Re and Da, he nonlinear drag force can be negleced, and he Brinkman model is applied in he porous region. A he sead sae, he analical soluion of he flow can be epressed as [, 4]: u d d L L, (5) d ep L L where K, d Ku K, and d u K e. In simulaions, F is se o be zero, and he fluid kineic viscosi is se o be.6 wih a N N 8 8 laice. The relaaion rae s is deermined b s.5 J cs (for free-fluid region, J ). The
12 veloci profiles prediced b he presen MRT-LB model for J and 4 a Da. and Re. are ploed in Fig. 4. As shown, he presen numerical resuls agree well wih he analical ones. The disconinui of he veloci-gradien a he porous medium/free-fluid inerface is well capured b he presen MRT-LB model wihou an special reamen for he boundar condiion a he inerface in simulaions. 3.3 Lid-driven flow in a square porous cavi In his subsecion, we appl he MRT-LB model o simulae he incompressible flow in a square cavi filled wih a fluid-sauraed porous medium. The op wall of he cavi moves from lef o righ wih a uniform veloci u, while he boom, righ, and lef walls of he cavi are fied. The Renolds number Re of his flow is defined as Re Lu, where L is he heigh of he cavi. In simulaions, we se., Re, J =, and he kineic viscosi is se o be.56 based on a N N 8 8 square mesh. The nonequilibrium erapolaion scheme is adoped o rea i he veloci boundar condiions of f on he four walls. The veloci momens are se o be heir equilibrium values, and he veloci field is iniialized o be wih a consan pressure p a. The sreamlines prediced b he presen MRT-LB model for differen Darc numbers wih. and Re are illusraed in Fig. 5. From he figure we can see ha, as Da decreases, he vore inside he cavi becomes weaker, and he boundar laer near he op wall becomes hinner. This can be aribued o he lower permeabili ( K = DaL ) of he medium which resuls in lower fluid veloci. To be more informaive, he horizonal veloci componen u in he verical midplane L and he verical veloci componen u in he horizonal midplane L are ploed in Fig. 6. The FD soluions in Ref. [] are also ploed in Fig. 6 for comparison. Obviousl, he presen numerical resuls are in good agreemen wih he FD soluions.
13 As and Da ends o infini, he presen MRT-LB model reduces o he incompressible MRT-LB model for free-fluid flows wihou porous media. We now appl he MRT-LB model o he lid-driven cavi flow wihou porous media. In simulaions, we se.999, 6 Da, and u.. In Fig. 7, he veloci profiles hrough he cener of he cavi are ploed for differen Renolds numbers. The benchmark soluions of Ghia e al. [4] are also presened for comparison. Ecellen agreemen can be found beween he presen MRT-LB resuls and he benchmark soluions. 3.4 Naural convecion flow in a square cavi filled wih porous media Naural convecion flow in a square cavi filled wih a porous medium (see Fig. 8) has been sudied eensivel b man researchers [4, 7, 43, 44] based on he generalized model. For his problem, he veloci field is solved b he presen MRT-LB model and he emperaure field is solved b he DQ5 MRT-LB model [44] (see Appendi B for deails). The op and boom walls of he cavi are hermall insulaed, while he lef and righ verical walls are kep a consan bu differen emperaures T h and T c ( T h T ), respecivel. The Prandl number Pr and Raleigh number Ra c of his flow are defined as Pr and Ra T T T is 3 g TL Pr, respecivel, where h c he emperaure difference, is he hermal diffusivi of he fluid, g is he graviaional acceleraion, is he hermal epansion coefficien, and L is he disance beween he walls. The bod force a is defined as a j, where T T T = g T T is he reference emperaure, h c and j is he uni vecor in he -direcion. In simulaions, we se Pr, J, = (hermal capaci raio), e ( e is he effecive hermal diffusivi), ζ, ζ ζ T, ζ3 ζ4., Th, and T c =. According o Refs. [3, 45], he dimensionless relaaion imes of he veloci and emperaure fields can be deermined as
14 MaJL 3Pr, c Ra T cs cst.5, (6) J Pr respecivel, where Ma U c 3U is he Mach number ( U g TL is he characerisic s veloci, and usuall Ma.3 ), and c c 5 ( c st is he sound speed of he DQ5 model). The st grid sizes of,, and 5 5 are emploed for Da, 4, and 6, respecivel. Fig. 9 illusraes he sreamlines and isoherms for differen Raleigh numbers and Darc numbers (Darc-Raleigh number * Ra RaDa = ) wih Pr and.6. From he figure we can observe ha, as Da decreases, he veloci and hermal boundar laers become hinner near he ho and cold verical walls. As Da increases o, he sreamlines and isoherms are less crowded near he verical walls and more convecive miing occurs inside he cavi. To quanif he resuls, he average Nussel numbers of he lef verical wall are calculaed and lised in Table. The numerical resuls given b Nihiarasu e al. [4, 7] using he finie elemen mehod are also included in Table for comparison. To sum up, he numerical resuls of he presen MRT-LB model agree well wih hose resuls repored in previous sudies. 4. Conclusion In his paper, an MRT-LB model wih he eigh-b-eigh collision mari has been presened for simulaing incompressible flows in porous media a he REV scale. The ke poin of he MRT-LB model is o include he porosi ino he equilibrium momens, and add a forcing erm o he MRT-LB equaion in he momen space o accoun for he Darc (linear) and Forchheimer (nonlinear) drag forces of he solid mari based on he generalized model. Through he Chapman-Enskog analsis in he momen space, he generalized Navier-Sokes equaions can be recovered eacl wihou arificial compressible errors. Numerical simulaions of he porous Poiseuille flow, porous Couee flow,
15 lid-driven flow in square porous cavi, and naural convecion flow in a square porous cavi have been carried ou o demonsrae he presen MRT-LB model. The numerical resuls of he presen MRT-LB model agree well wih he analical soluions and/or oher numerical soluions repored in previous sudies. Acknowledgemens This work was suppored b he Naional Ke Basic Research Program of China (973 Program) (3CB834). References [] K. Vafai, Handbook of Porous Media, Marcel Dekker, New York,. [] I. Pop, D.B. Ingham, Convecive Hea Transfer: Mahemaical and Compuaional Modeling of Viscous Fluids and Porous Media, Pergamon, Oford,. [3] D.A. Nield, A. Bejan, Convecion in Porous Media, hird ed., Springer, New York, 6. [4] P. Nihiarasu, K.N. Seeharamu, T. Sundararajan, Naural convecive hea ransfer in a fluid sauraed variable porosi medium, In. J. Hea Mass Transfer 4 (997) [5] K. Vafai, C.L. Tien, Boundar and ineria effecs on flow and hea ransfer in porous media, In. J. Hea Mass Transfer 4 (98) [6] C.T. Hsu, P. Cheng, Thermal dispersion in a porous medium, In. J. Hea Mass Transfer 33 (99) [7] P. Nihiarasu, K. Ravindran, A new semi-implici ime sepping procedure for buoanc driven flow in a fluid sauraed porous medium, Compu. Meh. Appl. Mech. Eng. 65 (998)
16 [8] K. Vafai, S.J. Kim, On he limiaions of Brinkman-Forchheimer-Eended Darc equaion, In. J. Hea Fluid Flow 6 (995) -5. [9] U. Frisch, B. Hasslacher, Y. Pomeau, Laice-gas auomaa for he Navier-Sokes equaion, Phs. Rev. Le. 56 (986) 55. [] R. Benzi, S. Succi, M. Vergassola, The laice Bolzmann equaion: heor and applicaions, Phs. Rep. (99) [] S. Chen, G. D. Doolen, Laice Bolzmann mehod for fluid flows, Annu. Rev. Fluid Mech. 3 (998) [] Y. Wang, Y. L. He, T. S. Zhao, G. H. Tang, W. Q. Tao, Implici-eplici finie-difference laice Bolzmann mehod for compressible flows, In. J. Mod. Phs. C 8 (7) [3] Y.L. He, Y. Wang, Q. Li, Laice Bolzmann Mehod: Theor and Applicaions, Science Press, Beijing, 9. [4] Q. Li, K.H. Luo, X.J. Li, Laice Bolzmann mehod for relaivisic hdrodnamics: Issues on conservaion law of paricle number and disconinuiies, Phs. Rev. D 86 () [5] Y.L. He, Q. Liu, Q. Li, Three-dimensional finie-difference laice Bolzmann model and is applicaion o inviscid compressible flows wih shock waves, Phsica A 39 (3) [6] Y.X. Xia, Y.H. Qian, Laice Bolzmann simulaion for forced wo-dimensional urbulence, Phs. Rev. E 9() (4) 34. [7] S. Succi, Laice Bolzmann across scales: from urbulence o DNA ranslocaion, Eur. Phs. J. B 64 (8) [8] M.A.A. Spaid, F.R. Phelan, Jr., Laice Bolzmann mehods for modeling microscale flow in fibrous porous media, Phs. Fluids 9 (997)
17 [9] O. Dardis, J. McCloske, Laice Bolzmann scheme wih real numbered solid densi for he simulaion of flow in porous media, Phs. Rev. E 57 (998) [] D.M. Freed, Laice-Bolzmann mehod for macroscopic porous media modeling, In. J. Mod. Phs. C 9 (998) [] N.S. Mars, Improved approimaion of he Brinkman equaion using a laice Bolzmann mehod, Phs. Fluids 6 () [] Z. Guo, T.S. Zhao, Laice Bolzmann model for incompressible flows hrough porous media, Phs. Rev. E 66 () [3] Q. Kang, D. Zhang, S. Chen, Unified laice Bolzmann mehod for flow in muliscale porous media, Phs. Rev. E 66 () [4] S. Anwar, M.C. Sukop, Regional scale ransien groundwaer flow modeling using laice Bolzmann mehods, Compu. Mah. Appl. 58 (9) 5-3. [5] S.D.C. Walsh, H. Burwinkle, M.O. Saar, A new parial-bounceback laice-bolzmann mehod for fluid flow hrough heerogeneous media, Compu. Geosci. 35 (9) [6] P.L. Bhanagar, E.P. Gross, M. Krook, A model for collision processes in gases. I. Small ampliude processes in charged and neural one-componen ssems, Phs. Rev. 94 (954) [7] P. Lallemand, L.-S. Luo, Theor of he laice Bolzmann mehod: Dispersion, dissipaion, isorop, Galilean invariance, and sabili, Phs. Rev. E 6 () [8] D. d Humières, I. Ginzburg, M. Krafczk, P. Lallemand, L.-S. Luo, Muliple relaaion ime laice Bolzmann models in hree dimensions, Philos. Trans. R. Soc. Lond. A 36 () [9] Q. Li, Y.L. He, G.H. Tang, W.Q. Tao, Improved aismmeric laice Bolzmann scheme, Phs.
18 Rev. E 8 () [3] I. Ginzburg, D. d Humières, Mulireflecion boundar condiions for laice Bolzmann models, Phs. Rev. E 68 (3) [3] D. d Humières, Generalized laice-bolzmann equaions, in: B.D. Shizgal, D.P. Weaver (Eds.), Rarefied Gas Dnamics: Theor and Simulaions, in: Prog. Asronau. Aeronau., Vol. 59, AIAA, Washingon, DC, 99, pp [3] R. Du, B. Shi, Incompressible MRT laice Bolzmann model wih eigh velociies in D space, In. J. Mod. Phs. C (9) [33] S. Ergun, Fluid flow hrough packed columns, Chem. Eng. Prog. 48 (95) [34] K. Vafai, Convecive flow and hea ransfer in variable-porosi media, J. Fluid Mech. 47 (984) [35] M.E. MaCracken, J. Abraham, Muliple-relaaion-ime laice-bolzmann model for muliphase flow, Phs. Rev. E 7 (5) 367. [36] Q. Li, K.H. Luo, X.J. Li, Laice Bolzmann modeling of muliphase flows a large densi raio wih an improved pseudopoenial model, Phs. Rev. E, 87 (3) 533. [37] Y.H. Qian, D. d Humeier, P. Lallemand, Laice BGK models for Navier-Sokes equaion, Europhs. Le. 7 (99) [38] S. Chapman and T. G. Cowling, The Mahemaical Theor of Non-Uniform Gases, Cambridge Universi Press, London, 97. [39] Z. Guo, B. Shi, N. Wang, Laice BGK model for incompressible Navier-Sokes equaion, J. Compu. Phs. 65 () [4] Z.L. Guo, C.G. Zheng, B.C. Shi, Non-equilibrium erapolaion mehod for veloci and pressure
19 boundar condiions in he laice Bolzmann mehod, Chin. Phs. () 366. [4] N. Mars, D.P. Benz, E.J. Garboczi, Compuer simulaion sud of he effecive viscosi in Brinkman s equaion, Phs. Fluids 6 (994) [4] U. Ghia, K.N. Ghia, C.T. Shin, High-Re soluions for incompressible flow using he Navier-Sokes equaions and a muligrid mehod, J. Compu. Phs. 48 (98) [43] Z. Guo, T.S. Zhao, A laice Bolzmann model for convecion hea ransfer in porous media, Numer. Hea Transfer B 47 (5) [44] Q. Liu, Y.L. He, Q. Li, W.Q. Tao, A muliple-relaaion-ime laice Bolzmann model for convecion hea ransfer in porous media, In. J. Hea Mass Transfer 73 (4) [45] Q. Li, Y.L. He, Y. Wang, G.H. Tang, An improved hermal laice Bolzmann model for flows wihou viscous hea dissipaion and compression work, In. J. Mod. Phs. C 9 (8) 5-5. Appendi A: Chapman-Enskog analsis of he DQ8 MRT-LB model The Chapman-Enskog epansion mehod [3, 35, 38] is adoped o derive he generalized Navier-Sokes equaions () and () from he DQ8 MRT-LB model. To his end, he following epansions in ime and space are inroduced: n f f n e, e, i i i i i n! () () i i i i, (A.a) f f f f, (A.b), j j, S S, F F, (A.c) where is a small epansion parameer, S,, S S T 8, F F, F. Wih he above epansions, we can derive he following equaions from Eq. (6) as consecuive orders of he parameer in he momen space as
20 : eq m m, (A.a) where : j j : D MD M I C (, Λ D m Λ m I S, (A.b) Λ Λ m D I m D I S Λ m, (A.c) j ), D I jdiag e j, e j,, e8 j, j eij C M I M, Λ Λ, and m T, e, F, q, F, q, p, p. (A.3) C and C can be given eplicil b C, C 3. (A.4) From Eq. (A.b), he following equaions a he ime scale can be obained: 5 u p u u, 3 3 (A.5a) s p u u s e 3 3 S, (A.5b) u uu s3 u p s3 F 3 S, (A.5c) u u u u s4 u p s4 q 4 S, (A.5d) u u u s5 u p s5 F 5 S, (A.5e) u u u u s6 u p s6 q 6 S, (A.5f) u u 7 s u u s7 p S, (A.5g)
21 uu 8 s u u s8 p S. (A.5h) From Eq. (A.c), he following equaions a he ime scale can be obained: u p, (A.6a) s s s u F e p s p 4 s3 s s7 s8 S 3 S S7 S8 4, (A.6b) s s s u F p e s p 4 s5 s8 s s7 S 5 S8 S S7 4, (A.6c) Following he approach in Ref. [3], we can obain: u p, (A.7) equaion Combining Eqs. (A.5a) and (A.6a) wih Eq. (A.7) leads o he following incompressible coninui u u, (A.8) Noe ha e, p and p in Eqs. (A.6b) and (A.6c) are unknowns o be deermined. According o Eqs. (A.5b), (A.5g) and (A.5h), we can ge: s s e p u u 3 3 S, (A.9a) u u s7 s 7 p u u S, (A.9b) uu s8 s 8 p u u S. (A.9c) Neglecing he erms of order 3 O( u ) and higher-order erms of he form u ( u u ) j k k j, using Eqs. (A.5c) and (A.5e), we can obain:
22 u uf, (A.a) u uf, (A.b) u u uf uf. (A.c) Wih he above equaions, we can obain: 4 s uf u F s e u u, (A.a) 5 5 u F u F s 7 p u u s 7, (A.b) 3 s u 8 F u F s 8 p u u. (A.c) 3 Subsiuing Eq. (A.) ino Eq. (A.6), he following equaions a he ime scale can be derived: u u u u u 5 s 3 s7 u u, (A.a) 3 s8 u u u u u 3 s8 5 s u u. (A.b) 3 s7 Combining Eq. (A.) ( ime scale) wih Eqs. (A.5c) and (A.5e) ( ime scale), he following equaions can be derived ( ): u u u u p 5 s u u
23 u u u u F, (A.3a) 3 3 s7 s8 u u u u p u u 3 s8 u u u u F. (A.3b) 5 s 3 s7 Wih he aid of he incompressible coninui equaion u u, he incompressible generalized Navier-Sokes equaions () and () can be obained from Eq. (A.3). The effecive kineic viscosi is defined as c.5 wih s7 s8 s. e s Appendi B: DQ5 MRT-LB model The macroscopic emperaure equaion of naural convecion flow in porous media can be wrien as T u T e T, (B.) where is he hermal capaci raio, and e is he effecive hermal diffusivi. For he emperaure field, he DQ5 MRT-LB equaion is defined as eq g e, g, N,, n n, (B.) where N is a 5 5 orhogonal ransformaion mari, and is a diagonal relaaion mari. The boldface smbols, g, n, and eq n are 5-dimensional (column) vecors:, g,, g,,, g4, g, g,,, g4 4, g e e e, n,, n,,, n4, n eq eq eq eq 4 T, T, n, n,, n,,, n, T, T,
24 where g, is he emperaure disribuion funcion,, i n and eq, n are he veloci momens of he discree disribuion funcions g and he corresponding equilibrium momens, respecivel. In he DQ5 model, he five discree velociies e i,,, 4 are given b i,, i e i. (B.3) cos i,sin i c, i ~ 4 For he DQ5 model, he ransformaion mari N is given b N. (B.4) 4 The ransformaion mari N linearl ransforms he discree disribuion funcions g 5 o heir veloci momens 5 n : In he ssem of g equilibrium momens as [44] i i n Ng,, onl he emperaure g = N n. (B.5) 4 T n g is conserved quani. The i i n eq i,,, 4 for he veloci momens n i,,, 4 are defined n eq T, ut n, eq ut n, eq n eq 3 T, eq 4 n, (B.6) where is a consan ( 4 ). In he presen work, is se o be. The diagonal relaaion mari is given b: 3 4 =diag, ζ, ζ, ζ, ζ. (B.7) The effecive hermal diffusivi is defined as c.5 e e st T wih ζ ζ T and c 4 c. st
25 Figure Capions Fig.. Discree velociies of he DQ8 model. Fig.. Veloci profiles of he porous Poiseuille flow for differen Re and Da wih. and J. Solid lines represen he FD soluions [] and smbols represen he presen MRT-LB resuls. Fig. 3. Veloci profiles of he porous Couee flow for differen Re and Da wih. and J. Solid lines represen he FD soluions [] and smbols represen he presen MRT-LB resuls. Fig. 4. Veloci profiles of he porous Couee flow for differen viscosi raio J wih Re =. and Da =.. Solid lines represen he analical soluions from Eq. (5) and smbols represen he presen MRT-LB resuls. Fig. 5. Sreamlines of he lid-driven flow in a porous cavi for differen Da wih. and Re : (a) Da ; (b) Da 3 ; (c) Da 4. Fig. 6. Veloci profiles hrough he cener of he cavi: (a) horizonal veloci componen u in he verical midplane L ; (b) verical veloci componen u in he horizonal midplane L. Solid lines represen he FD soluions [] and smbols represen he presen MRT-LB resuls. Fig. 7. Veloci profiles hrough he cener of he cavi: (a) horizonal veloci componen u in he verical midplane L ; (b) verical veloci componen u in he horizonal midplane L. Solid lines represen he benchmark soluions [4] and smbols represen he presen MRT-LB resuls. Fig. 8. Naural convecion in a square cavi filled wih a porous medium. Fig. 9. Sreamlines and isoherms for differen Ra and Da wih Pr and.6 : (a) Da, 5 Ra ; (b) Da 4, 7 Ra ; (c) Da 6, 9 Ra.
26 Fig.. Discree velociies of he DQ8 model.
27 ..8 u/u.6.4. Da= 5, =. Re=. Re= Re=5 Re= Lines: FD soluions /L (a)..8 u/u.6.4. Re=., =. Da= 6 Da= 4 Da= 3 Lines: FD soluions /L (b) Fig.. Veloci profiles of he porous Poiseuille flow for differen Re and Da wih. and J. Solid lines represen he FD soluions [] and smbols represen he presen MRT-LB resuls.
28 ..8 u/u.6.4. Da=., =. Re=. Re=5 Re= Lines: FD soluions /L (a)..8.6 /L.4. Re=, =. Da= Da= 3 Da= 5 Lines: FD soluions....4 u/u.6.8. (b) Fig. 3. Veloci profiles of he porous Couee flow for differen Re and Da wih. and J. Solid lines represen he FD soluions [] and smbols represen he presen MRT-LB resuls.
29 ..8 /L.6.4. Re=., Da=. J= J=4 Lines: Analical soluions u/u Fig. 4. Veloci profiles of he porous Couee flow for differen viscosi raio J wih Re =. and Da =.. Solid lines represen he analical soluions from Eq. (5) and smbols represen he presen MRT-LB resuls.
30 (a) Da (b) Da (c) Da 4 Fig. 5. Sreamlines of he lid-driven flow in a porous cavi for differen Da wih. and Re : (a) Da ; (b) Da 3 ; (c) Da 4.
31 ..8 /L.6.4. Re=, =. Da= Da= 3 Da= 4 Lines: FD soluions u(l/, )/u (a).5..5 v(,l/)/u Re=, =. Da= Da= 3 Da= 4 Lines: FD soluions /L (b) Fig. 6. Veloci profiles hrough he cener of he cavi: (a) horizonal veloci componen u in he verical midplane L ; (b) verical veloci componen u in he horizonal midplane L. Solid lines represen he FD soluions [] and smbols represen he presen MRT-LB resuls.
32 ..8.6 /L.4. Da= 6, =.999 Re= Re=4 Re= Lines: Ghia e al. [4] u(l/, )/u (a).4. v(,l/)/u Da= 6, =.999 Re= Re=4 Re= Lines: Ghia e al. [4] /L (b) Fig. 7. Veloci profiles hrough he cener of he cavi: (a) horizonal veloci componen u in he verical midplane L ; (b) verical veloci componen u in he horizonal midplane L. Solid lines represen he benchmark soluions [4] and smbols represen he presen MRT-LB resuls.
33 Adiabaic Porous Medium Th g Tc Adiabaic Fig. 8. Naural convecion in a square cavi filled wih a porous medium.
34 (a) Da, 5 Ra (b) Da 4, 7 Ra (c) Da 6, 9 Ra Fig. 9. Sreamlines and isoherms for differen Ra and Da wih Pr and.6 : (a) Da, 5 Ra ; (b) Da 4, 7 Ra ; (c) Da 6, 9 Ra.
35 Table Capion Table. Comparisons of he average Nussel numbers for differen Ra, Da and wih Pr.. Da Table. Comparisons of he average Nussel numbers for differen Ra, Da and wih Pr.. Ra Ref.[4] Ref.[7] Presen Ref.[4] Ref.[7] Presen Ref.[4] Ref.[7] Presen
Flow-Induced Vibration Analysis of Supported Pipes with a Crack
Flow-Induced Vibraion Analsis of Suppored Pipes wih a Crack Jin-Huk Lee, Samer Masoud Al-Said Deparmen of Mechanical Engineering American Universi of Sharjah, UAE Ouline Inroducion and Moivaion Aeroacousicall
More informationTHE EFFECT OF SUCTION AND INJECTION ON UNSTEADY COUETTE FLOW WITH VARIABLE PROPERTIES
Kragujevac J. Sci. 3 () 7-4. UDC 53.5:536. 4 THE EFFECT OF SUCTION AND INJECTION ON UNSTEADY COUETTE FLOW WITH VARIABLE PROPERTIES Hazem A. Aia Dep. of Mahemaics, College of Science,King Saud Universiy
More informationPhysics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle
Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,
More informationImproved Approximate Solutions for Nonlinear Evolutions Equations in Mathematical Physics Using the Reduced Differential Transform Method
Journal of Applied Mahemaics & Bioinformaics, vol., no., 01, 1-14 ISSN: 179-660 (prin), 179-699 (online) Scienpress Ld, 01 Improved Approimae Soluions for Nonlinear Evoluions Equaions in Mahemaical Physics
More informationDiffusion & Viscosity: Navier-Stokes Equation
4/5/018 Diffusion & Viscosiy: Navier-Sokes Equaion 1 4/5/018 Diffusion Equaion Imagine a quaniy C(x,) represening a local propery in a fluid, eg. - hermal energy densiy - concenraion of a polluan - densiy
More informationREVERSE COMPUTATION OF FORCED CONVECTION HEAT TRANSFER USING ADJOINT FORMULATION
REVERSE COMPUTATION OF FORCED CONVECTION HEAT TRANSFER USING ADJOINT FORMULATION Kazunari Momose and Hideo Kimoo Graduae School of Engineering Science, Osaka Universiy Osaka 56-853, Japan E-mail: momose@me.es.osaka-u.ac.jp
More informationThe fundamental mass balance equation is ( 1 ) where: I = inputs P = production O = outputs L = losses A = accumulation
Hea (iffusion) Equaion erivaion of iffusion Equaion The fundamenal mass balance equaion is I P O L A ( 1 ) where: I inpus P producion O oupus L losses A accumulaion Assume ha no chemical is produced or
More informationUnsteady MHD Second Grade Fluids Flow in a Porous Medium with Ramped Wall Temperature
Unsead MHD Second Grade Fluids Flow in a Porous Medium wih Ramped Wall Temperaure ZULKHIBRI ISMAIL Universii Malasia Pahang Facul of Indusrial Sciences & Technolog Lebuh Raa Tun Razak 6300 Kuanan Pahang
More informationNumerical investigation of Ranque-Hilsch energy separation effect A.S. Noskov 1,a, V.N. Alekhin 1,b, A.V. Khait 1,a
Applied Mechanics and Maerials Online: 2013-01-11 ISSN: 1662-7482, Vol. 281, pp 355-358 doi:10.4028/www.scienific.ne/amm.281.355 2013 Trans Tech Publicaions, Swizerland Numerical invesigaion of Ranque-Hilsch
More informationTheory of! Partial Differential Equations!
hp://www.nd.edu/~gryggva/cfd-course/! Ouline! Theory o! Parial Dierenial Equaions! Gréar Tryggvason! Spring 011! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and
More informationAnd the solution to the PDE problem must be of the form Π 1
5. Self-Similar Soluions b Dimensional Analsis Consider he diffusion problem from las secion, wih poinwise release (Ref: Bluman & Cole, 2.3): c = D 2 c x + Q 0δ(x)δ() 2 c(x,0) = 0, c(±,) = 0 Iniial release
More informationThermal Modeling of a Honeycomb Reformer including Radiative Heat Transfer
Thermal Modeling of a Honeycomb Reformer including Radiaive Hea Transfer J Schöne *1, A Körnig 1, W Becer 1 and A Michaelis 1 1 Fraunhofer IKTS, Winerbergsraße 8, 0177 Dresden *Corresponding auhor: jaobschoene@isfraunhoferde
More informationThe motions of the celt on a horizontal plane with viscous friction
The h Join Inernaional Conference on Mulibody Sysem Dynamics June 8, 18, Lisboa, Porugal The moions of he cel on a horizonal plane wih viscous fricion Maria A. Munisyna 1 1 Moscow Insiue of Physics and
More informationTheory of! Partial Differential Equations-I!
hp://users.wpi.edu/~grear/me61.hml! Ouline! Theory o! Parial Dierenial Equaions-I! Gréar Tryggvason! Spring 010! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and
More informationCH.7. PLANE LINEAR ELASTICITY. Continuum Mechanics Course (MMC) - ETSECCPB - UPC
CH.7. PLANE LINEAR ELASTICITY Coninuum Mechanics Course (MMC) - ETSECCPB - UPC Overview Plane Linear Elasici Theor Plane Sress Simplifing Hpohesis Srain Field Consiuive Equaion Displacemen Field The Linear
More informationIB Physics Kinematics Worksheet
IB Physics Kinemaics Workshee Wrie full soluions and noes for muliple choice answers. Do no use a calculaor for muliple choice answers. 1. Which of he following is a correc definiion of average acceleraion?
More informationIntroduction to Physical Oceanography Homework 5 - Solutions
Laure Zanna //5 Inroducion o Phsical Oceanograph Homework 5 - Soluions. Inerial oscillaions wih boom fricion non-selecive scale: The governing equaions for his problem are This ssem can be wrien as where
More informationNavneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi
Creep in Viscoelasic Subsances Numerical mehods o calculae he coefficiens of he Prony equaion using creep es daa and Herediary Inegrals Mehod Navnee Saini, Mayank Goyal, Vishal Bansal (23); Term Projec
More informationME 391 Mechanical Engineering Analysis
Fall 04 ME 39 Mechanical Engineering Analsis Eam # Soluions Direcions: Open noes (including course web posings). No books, compuers, or phones. An calculaor is fair game. Problem Deermine he posiion of
More informationJournal of Chemical and Pharmaceutical Research, 2013, 5(9): Research Article. Thermal analysis of diesel engine piston
Available online www.jocpr.com Journal of Chemical and Pharmaceuical Research, 13, 5(9):388-393 Research Aricle ISSN : 975-7384 CODEN(USA) : JCPRC5 hermal analsis of diesel engine pison Honguan Zhang,
More informationNumerical Dispersion
eview of Linear Numerical Sabiliy Numerical Dispersion n he previous lecure, we considered he linear numerical sabiliy of boh advecion and diffusion erms when approimaed wih several spaial and emporal
More informationStructural Dynamics and Earthquake Engineering
Srucural Dynamics and Earhquae Engineering Course 1 Inroducion. Single degree of freedom sysems: Equaions of moion, problem saemen, soluion mehods. Course noes are available for download a hp://www.c.up.ro/users/aurelsraan/
More informationStability and Bifurcation in a Neural Network Model with Two Delays
Inernaional Mahemaical Forum, Vol. 6, 11, no. 35, 175-1731 Sabiliy and Bifurcaion in a Neural Nework Model wih Two Delays GuangPing Hu and XiaoLing Li School of Mahemaics and Physics, Nanjing Universiy
More informationExperiment 123 Determination of the sound wave velocity with the method of Lissajous figures
perimen 3 Deerminaion of he sound wave veloci wih he mehod of Lissajous figures The aim of he eercise To sud acousic wave propagaion in he air To deermine of he sound wave veloci in he air Mehodolog of
More informationESTIMATION OF DYNAMIC PANEL DATA MODELS WHEN REGRESSION COEFFICIENTS AND INDIVIDUAL EFFECTS ARE TIME-VARYING
Inernaional Journal of Social Science and Economic Research Volume:02 Issue:0 ESTIMATION OF DYNAMIC PANEL DATA MODELS WHEN REGRESSION COEFFICIENTS AND INDIVIDUAL EFFECTS ARE TIME-VARYING Chung-ki Min Professor
More information1. VELOCITY AND ACCELERATION
1. VELOCITY AND ACCELERATION 1.1 Kinemaics Equaions s = u + 1 a and s = v 1 a s = 1 (u + v) v = u + as 1. Displacemen-Time Graph Gradien = speed 1.3 Velociy-Time Graph Gradien = acceleraion Area under
More informationNUMERICAL SIMULATION OF A LIQUID SODIUM TURBULENT FLOW OVER A BACKWARD FACING STEP WITH A FOUR PARAMETER LOGARITHMIC TURBULENCE MODEL
h European Conference on Compuaional Mechanics (ECCM ) 7h European Conference on Compuaional Fluid Dynamics (ECFD 7) 11 June, Glasgow, UK NUMERICAL SIMULATION OF A LIQUID SODIUM TURBULENT FLOW OVER A BACKWARD
More informationAlmost Unbiased Estimator for Estimating Population Mean Using Known Value of Some Population Parameter(s)
Almos Unbiased Esimaor for Esimaing Populaion Mean Using Known Value of Some Populaion Parameers Rajesh Singh Deparmen of Saisics, Banaras Hindu Universi U.P., India rsinghsa@ahoo.com Mukesh Kumar Deparmen
More informationVehicle Arrival Models : Headway
Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where
More informationFormulation of the Stress Distribution Due to a Concentrated Force Acting on the Boundary of Viscoelastic Half-Space
Formulaion of he Sress Disribuion Due o a Concenraed Force Acing on he Boundar of Viscoelasic Half-Space Yun eng and Debao Zhou Deparmen of Mechanical and Indusrial Engineering Universi of Minnesoa, Duluh
More informationElementary Differential Equations and Boundary Value Problems
Elemenar Differenial Equaions and Boundar Value Problems Boce. & DiPrima 9 h Ediion Chaper 1: Inroducion 1006003 คณ ตศาสตร ว ศวกรรม 3 สาขาว ชาว ศวกรรมคอมพ วเตอร ป การศ กษา 1/2555 ผศ.ดร.อร ญญา ผศ.ดร.สมศ
More informationSub Module 2.6. Measurement of transient temperature
Mechanical Measuremens Prof. S.P.Venkaeshan Sub Module 2.6 Measuremen of ransien emperaure Many processes of engineering relevance involve variaions wih respec o ime. The sysem properies like emperaure,
More informationASTR415: Problem Set #5
ASTR45: Problem Se #5 Curran D. Muhlberger Universi of Marland (Daed: April 25, 27) Three ssems of coupled differenial equaions were sudied using inegraors based on Euler s mehod, a fourh-order Runge-Kua
More informationThe role of the error function in three-dimensional singularly perturbed convection-diffusion problems with discontinuous data
The role of he error funcion in hree-dimensional singularl perurbed convecion-diffusion problems wih disconinuous daa José Luis López García, Eser Pérez Sinusía Depo. de Maemáica e Informáica, U. Pública
More informationModule 2 F c i k c s la l w a s o s f dif di fusi s o i n
Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms
More informationIntegration of the equation of motion with respect to time rather than displacement leads to the equations of impulse and momentum.
Inegraion of he equaion of moion wih respec o ime raher han displacemen leads o he equaions of impulse and momenum. These equaions greal faciliae he soluion of man problems in which he applied forces ac
More informationHeat Transfer. Revision Examples
Hea Transfer Revision Examples Hea ransfer: energy ranspor because of a emperaure difference. Thermal energy is ransferred from one region o anoher. Hea ranspor is he same phenomena lie mass ransfer, momenum
More informationN. Sandeep 1 and V. Sugunamma 2
Journal of Applied Fluid Mechanics, Vol. 7, No., pp. 75-86, 4. Available online a www.jafmonline.ne, ISSN 735-357, EISSN 735-3645. Radiaion and Inclined Magneic Field Effecs on Unseady Hydromagneic Free
More informationThe Application of Optimal Homotopy Asymptotic Method for One-Dimensional Heat and Advection- Diffusion Equations
Inf. Sci. Le., No., 57-61 13) 57 Informaion Sciences Leers An Inernaional Journal hp://d.doi.org/1.1785/isl/ The Applicaion of Opimal Homoopy Asympoic Mehod for One-Dimensional Hea and Advecion- Diffusion
More informationForcing scheme in pseudopotential lattice Boltzmann model for multiphase flows
Forcing scheme in pseudopoenial laice Bolzmann model for muliphase flows Q. Li, K. H. Luo, *, and X. J. Li 3 School of Engineering Sciences, Universiy of Souhampon, Souhampon SO7 BJ, Unied Kingdom Cener
More information10. State Space Methods
. Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he
More informationVanishing Viscosity Method. There are another instructive and perhaps more natural discontinuous solutions of the conservation law
Vanishing Viscosiy Mehod. There are anoher insrucive and perhaps more naural disconinuous soluions of he conservaion law (1 u +(q(u x 0, he so called vanishing viscosiy mehod. This mehod consiss in viewing
More informationVerification of a CFD benchmark solution of transient low Mach number flows with Richardson extrapolation procedure 1
Verificaion of a CFD benchmark soluion of ransien low Mach number flows wih Richardson exrapolaion procedure S. Beneboula, S. Gounand, A. Beccanini and E. Suder DEN/DANS/DMS/STMF Commissaria à l Energie
More informationHaar Wavelet Operational Matrix Method for Solving Fractional Partial Differential Equations
Copyrigh 22 Tech Science Press CMES, vol.88, no.3, pp.229-243, 22 Haar Wavele Operaional Mari Mehod for Solving Fracional Parial Differenial Equaions Mingu Yi and Yiming Chen Absrac: In his paper, Haar
More informationOutline of Topics. Analysis of ODE models with MATLAB. What will we learn from this lecture. Aim of analysis: Why such analysis matters?
of Topics wih MATLAB Shan He School for Compuaional Science Universi of Birmingham Module 6-3836: Compuaional Modelling wih MATLAB Wha will we learn from his lecure Aim of analsis: Aim of analsis. Some
More informationAdvanced FDTD Algorithms
EE 5303 Elecromagneic Analsis Using Finie Difference Time Domain Lecure #5 Advanced FDTD Algorihms Lecure 5 These noes ma conain coprighed maerial obained under fair use rules. Disribuion of hese maerials
More informationMOMENTUM CONSERVATION LAW
1 AAST/AEDT AP PHYSICS B: Impulse and Momenum Le us run an experimen: The ball is moving wih a velociy of V o and a force of F is applied on i for he ime inerval of. As he resul he ball s velociy changes
More informationTurbulent Flows. Computational Modelling of Turbulent Flows. Overview. Turbulent Eddies and Scales
School of Mechanical Aerospace and Civil Engineering Turbulen Flows As noed above, using he mehods described in earlier lecures, he Navier-Sokes equaions can be discreized and solved numerically on complex
More informationFractional Method of Characteristics for Fractional Partial Differential Equations
Fracional Mehod of Characerisics for Fracional Parial Differenial Equaions Guo-cheng Wu* Modern Teile Insiue, Donghua Universiy, 188 Yan-an ilu Road, Shanghai 51, PR China Absrac The mehod of characerisics
More informationComputation of the Effect of Space Harmonics on Starting Process of Induction Motors Using TSFEM
Journal of elecrical sysems Special Issue N 01 : November 2009 pp: 48-52 Compuaion of he Effec of Space Harmonics on Saring Process of Inducion Moors Using TSFEM Youcef Ouazir USTHB Laboraoire des sysèmes
More informationdt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.
Mah Rahman Exam Review Soluions () Consider he IVP: ( 4)y 3y + 4y = ; y(3) = 0, y (3) =. (a) Please deermine he longes inerval for which he IVP is guaraneed o have a unique soluion. Soluion: The disconinuiies
More informationMulti-scale 2D acoustic full waveform inversion with high frequency impulsive source
Muli-scale D acousic full waveform inversion wih high frequency impulsive source Vladimir N Zubov*, Universiy of Calgary, Calgary AB vzubov@ucalgaryca and Michael P Lamoureux, Universiy of Calgary, Calgary
More information2002 November 14 Exam III Physics 191
November 4 Exam III Physics 9 Physical onsans: Earh s free-fall acceleraion = g = 9.8 m/s ircle he leer of he single bes answer. quesion is worh poin Each 3. Four differen objecs wih masses: m = kg, m
More informationSTUDY ON THE AIR MOVEMENT CHARACTER IN SOLAR WALL SYSTEM. Y Li, X Duanmu, Y Sun, J Li and H Jia
Proceedings: Building Simulaion 007 SUDY ON HE AIR MOVEMEN CHARACER IN SOLAR WALL SYSEM Y Li, X Duanmu, Y Sun, J Li and H Jia College of Archiecure and Civil Engineering, Beijing Universiy of echnology,
More informationarxiv: v1 [math.fa] 3 Jan 2019
DAMPED AND DIVERGENCE EXACT SOLUTIONS FOR THE DUFFING EQUATION USING LEAF FUNCTIONS AND HYPERBOLIC LEAF FUNCTIONS A PREPRINT arxiv:9.66v [mah.fa] Jan 9 Kazunori Shinohara Deparmen of Mechanical Sysems
More informationWEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x
WEEK-3 Reciaion PHYS 131 Ch. 3: FOC 1, 3, 4, 6, 14. Problems 9, 37, 41 & 71 and Ch. 4: FOC 1, 3, 5, 8. Problems 3, 5 & 16. Feb 8, 018 Ch. 3: FOC 1, 3, 4, 6, 14. 1. (a) The horizonal componen of he projecile
More informationProblem Set #1. i z. the complex propagation constant. For the characteristic impedance:
Problem Se # Problem : a) Using phasor noaion, calculae he volage and curren waves on a ransmission line by solving he wave equaion Assume ha R, L,, G are all non-zero and independen of frequency From
More informationChapter 15: Phenomena. Chapter 15 Chemical Kinetics. Reaction Rates. Reaction Rates R P. Reaction Rates. Rate Laws
Chaper 5: Phenomena Phenomena: The reacion (aq) + B(aq) C(aq) was sudied a wo differen emperaures (98 K and 35 K). For each emperaure he reacion was sared by puing differen concenraions of he 3 species
More informationKINEMATICS IN ONE DIMENSION
KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec
More informationVariational Iteration Method for Solving System of Fractional Order Ordinary Differential Equations
IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 1, Issue 6 Ver. II (Nov - Dec. 214), PP 48-54 Variaional Ieraion Mehod for Solving Sysem of Fracional Order Ordinary Differenial
More informationOn Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature
On Measuring Pro-Poor Growh 1. On Various Ways of Measuring Pro-Poor Growh: A Shor eview of he Lieraure During he pas en years or so here have been various suggesions concerning he way one should check
More informationHomotopy Perturbation Method for Solving Some Initial Boundary Value Problems with Non Local Conditions
Proceedings of he World Congress on Engineering and Compuer Science 23 Vol I WCECS 23, 23-25 Ocober, 23, San Francisco, USA Homoopy Perurbaion Mehod for Solving Some Iniial Boundary Value Problems wih
More informationSome Basic Information about M-S-D Systems
Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,
More informationLecture 4 Kinetics of a particle Part 3: Impulse and Momentum
MEE Engineering Mechanics II Lecure 4 Lecure 4 Kineics of a paricle Par 3: Impulse and Momenum Linear impulse and momenum Saring from he equaion of moion for a paricle of mass m which is subjeced o an
More informationKEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow
1 KEY Mah 4 Miderm I Fall 8 secions 1 and Insrucor: Sco Glasgow Please do NOT wrie on his eam. No credi will be given for such work. Raher wrie in a blue book, or on our own paper, preferabl engineering
More informationMotion along a Straight Line
chaper 2 Moion along a Sraigh Line verage speed and average velociy (Secion 2.2) 1. Velociy versus speed Cone in he ebook: fer Eample 2. Insananeous velociy and insananeous acceleraion (Secions 2.3, 2.4)
More informationSTATE-SPACE MODELLING. A mass balance across the tank gives:
B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing
More informationThe Contradiction within Equations of Motion with Constant Acceleration
The Conradicion wihin Equaions of Moion wih Consan Acceleraion Louai Hassan Elzein Basheir (Daed: July 7, 0 This paper is prepared o demonsrae he violaion of rules of mahemaics in he algebraic derivaion
More informationLecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.
Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in
More informationd 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3
and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or
More informationIMPROVED HYPERBOLIC FUNCTION METHOD AND EXACT SOLUTIONS FOR VARIABLE COEFFICIENT BENJAMIN-BONA-MAHONY-BURGERS EQUATION
THERMAL SCIENCE, Year 015, Vol. 19, No. 4, pp. 1183-1187 1183 IMPROVED HYPERBOLIC FUNCTION METHOD AND EXACT SOLUTIONS FOR VARIABLE COEFFICIENT BENJAMIN-BONA-MAHONY-BURGERS EQUATION by Hong-Cai MA a,b*,
More informationChapter 5 Kinematics
Chaper 5 Kinemaics In he firs place, wha do we mean b ime and space? I urns ou ha hese deep philosophical quesions have o be analzed ver carefull in phsics, and his is no eas o do. The heor of relaivi
More informationLinear Dynamic Models
Linear Dnamic Models and Forecasing Reference aricle: Ineracions beween he muliplier analsis and he principle of acceleraion Ouline. The sae space ssem as an approach o working wih ssems of difference
More informationHyperchaos Synchronization Between two Different Hyperchaotic Systems
ISSN 76-769, England, UK Journal of Informaion and Compuing Science Vo3, No., 8, pp. 73-8 Hperchaos Snchroniaion Beween wo Differen Hperchaoic Ssems Qiang Jia + Facul of Science, Jiangsu Universi, Zhenjiang,
More informationEXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE
Version April 30, 2004.Submied o CTU Repors. EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Per Krysl Universiy of California, San Diego La Jolla, California 92093-0085,
More informationUnsteady Mixed Convection Heat and Mass Transfer Past an Infinite Porous Plate with Thermophoresis Effect
Unseady Mixed Convecion Hea and Mass Transfer Pas an Infinie Porous Plae wih Thermophoresis Effec TARIQ AL-AZAB Mechanical Engineering Deparmen, Al-Al-Sal Communiy College Al-Balqa Applied Universiy P.O.Box
More information3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon
3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of
More informationwhere the coordinate X (t) describes the system motion. X has its origin at the system static equilibrium position (SEP).
Appendix A: Conservaion of Mechanical Energy = Conservaion of Linear Momenum Consider he moion of a nd order mechanical sysem comprised of he fundamenal mechanical elemens: ineria or mass (M), siffness
More informationModal identification of structures from roving input data by means of maximum likelihood estimation of the state space model
Modal idenificaion of srucures from roving inpu daa by means of maximum likelihood esimaion of he sae space model J. Cara, J. Juan, E. Alarcón Absrac The usual way o perform a forced vibraion es is o fix
More informationExact solution of the(2+1)-dimensional hyperbolic nonlinear Schrödinger equation by Adomian decomposition method
Malaa J Ma ((014 160 164 Exac soluion of he(+1-dimensional hperbolic nonlinear Schrödinger equaion b Adomian decomposiion mehod Ifikhar Ahmed, a, Chunlai Mu b and Pan Zheng c a,b,c College of Mahemaics
More informationModule 4: Time Response of discrete time systems Lecture Note 2
Module 4: Time Response of discree ime sysems Lecure Noe 2 1 Prooype second order sysem The sudy of a second order sysem is imporan because many higher order sysem can be approimaed by a second order model
More informationTHE FOURIER-YANG INTEGRAL TRANSFORM FOR SOLVING THE 1-D HEAT DIFFUSION EQUATION. Jian-Guo ZHANG a,b *
Zhang, J.-G., e al.: The Fourier-Yang Inegral Transform for Solving he -D... THERMAL SCIENCE: Year 07, Vol., Suppl., pp. S63-S69 S63 THE FOURIER-YANG INTEGRAL TRANSFORM FOR SOLVING THE -D HEAT DIFFUSION
More informationQ.1 Define work and its unit?
CHP # 6 ORK AND ENERGY Q.1 Define work and is uni? A. ORK I can be define as when we applied a force on a body and he body covers a disance in he direcion of force, hen we say ha work is done. I is a scalar
More informationAPPLICATION OF CHEBYSHEV APPROXIMATION IN THE PROCESS OF VARIATIONAL ITERATION METHOD FOR SOLVING DIFFERENTIAL- ALGEBRAIC EQUATIONS
Mahemaical and Compuaional Applicaions, Vol., No. 4, pp. 99-978,. Associaion for Scienific Research APPLICATION OF CHEBYSHEV APPROXIMATION IN THE PROCESS OF VARIATIONAL ITERATION METHOD FOR SOLVING DIFFERENTIAL-
More informationProgram: RFEM 5, RSTAB 8, RF-DYNAM Pro, DYNAM Pro. Category: Isotropic Linear Elasticity, Dynamics, Member
Verificaion Example Program: RFEM 5, RSTAB 8, RF-DYNAM Pro, DYNAM Pro Caegory: Isoropic Linear Elasiciy, Dynamics, Member Verificaion Example: 0104 Canilever Beam wih Periodic Exciaion 0104 Canilever Beam
More informationSection 3.8, Mechanical and Electrical Vibrations
Secion 3.8, Mechanical and Elecrical Vibraions Mechanical Unis in he U.S. Cusomary and Meric Sysems Disance Mass Time Force g (Earh) Uni U.S. Cusomary MKS Sysem CGS Sysem fee f slugs seconds sec pounds
More informationNovel FVM - LBM method for compressible flows
Dual Degree Projec Sage III Novel FVM - LBM mehod for compressible flows Under he guidance of Prof. Ami Agrawal and Prof. Bhalchandra Puranik Submied By: Arpi Agarwal (02D0036) 5 h Year Dual Degree Suden
More informationGENERALIZATION OF THE FORMULA OF FAA DI BRUNO FOR A COMPOSITE FUNCTION WITH A VECTOR ARGUMENT
Inerna J Mah & Mah Sci Vol 4, No 7 000) 48 49 S0670000970 Hindawi Publishing Corp GENERALIZATION OF THE FORMULA OF FAA DI BRUNO FOR A COMPOSITE FUNCTION WITH A VECTOR ARGUMENT RUMEN L MISHKOV Received
More informationRC, RL and RLC circuits
Name Dae Time o Complee h m Parner Course/ Secion / Grade RC, RL and RLC circuis Inroducion In his experimen we will invesigae he behavior of circuis conaining combinaions of resisors, capaciors, and inducors.
More informationDynamic Analysis of Damped Driven Pendulum using Laplace Transform Method
, ISSN 0974-570X (Online), ISSN 0974-578 (Prin), Vol. 6; Issue No. 3; Year 05, Copyrigh 05 by CESER PUBLICATIONS Dynamic Analysis of Damped Driven Pendulum using Laplace Transform Mehod M.C. Agarana and
More informationt 2 B F x,t n dsdt t u x,t dxdt
Evoluion Equaions For 0, fixed, le U U0, where U denoes a bounded open se in R n.suppose ha U is filled wih a maerial in which a conaminan is being ranspored by various means including diffusion and convecion.
More informationWeek 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem)
Week 1 Lecure Problems, 5 Wha if somehing oscillaes wih no obvious spring? Wha is ω? (problem se problem) Sar wih Try and ge o SHM form E. Full beer can in lake, oscillaing F = m & = ge rearrange: F =
More informationNMR Spectroscopy: Principles and Applications. Nagarajan Murali 1D - Methods Lecture 5
NMR pecroscop: Principles and Applicaions Nagarajan Murali D - Mehods Lecure 5 D-NMR To full appreciae he workings of D NMR eperimens we need o a leas consider wo coupled spins. omeimes we need o go up
More informationA NEW TECHNOLOGY FOR SOLVING DIFFUSION AND HEAT EQUATIONS
THERMAL SCIENCE: Year 7, Vol., No. A, pp. 33-4 33 A NEW TECHNOLOGY FOR SOLVING DIFFUSION AND HEAT EQUATIONS by Xiao-Jun YANG a and Feng GAO a,b * a School of Mechanics and Civil Engineering, China Universiy
More informationNotes on Kalman Filtering
Noes on Kalman Filering Brian Borchers and Rick Aser November 7, Inroducion Daa Assimilaion is he problem of merging model predicions wih acual measuremens of a sysem o produce an opimal esimae of he curren
More informationLecture 10: Wave equation, solution by spherical means
Lecure : Wave equaion, soluion by spherical means Physical modeling eample: Elasodynamics u (; ) displacemen vecor in elasic body occupying a domain U R n, U, The posiion of he maerial poin siing a U in
More information2001 November 15 Exam III Physics 191
1 November 15 Eam III Physics 191 Physical Consans: Earh s free-fall acceleraion = g = 9.8 m/s 2 Circle he leer of he single bes answer. quesion is worh 1 poin Each 3. Four differen objecs wih masses:
More informationBoundedness and Stability of Solutions of Some Nonlinear Differential Equations of the Third-Order.
Boundedness Sabili of Soluions of Some Nonlinear Differenial Equaions of he Third-Order. A.T. Ademola, M.Sc. * P.O. Arawomo, Ph.D. Deparmen of Mahemaics Saisics, Bowen Universi, Iwo, Nigeria. Deparmen
More informationSimple possibilities of thermal diffusivity estimation for small-sized samples, with a laser pulse heating and infrared cameras
0 h Inernaional Conference on Quaniaive InfraRed Thermograph Jul 7-30, 00, Québec (Canada) Simple possibiliies of hermal diffusivi esimaion for small-sied samples, wih a laser pulse heaing and infrared
More informationModelling traffic flow with constant speed using the Galerkin finite element method
Modelling raffic flow wih consan speed using he Galerin finie elemen mehod Wesley Ceulemans, Magd A. Wahab, Kur De Prof and Geer Wes Absrac A macroscopic level, raffic can be described as a coninuum flow.
More information