Numerical Simulation of Lid-Driven Cavity Flow Using the Lattice Boltzmann Method
|
|
- Kerrie Osborne
- 6 years ago
- Views:
Transcription
1 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 3, N. MUHAMAD AND K. SOPIAN 2 Department of Mechancal and Materal Engneerng 2 Solar Energy Research Insttute (SERI) Natonal Unversty of Malaysa 436 UKM Bang, Selangor MALAYSIA 3 Faculty of Mechancal Engneerng Technology Unversty of Malaysa 83 Skuda, Johor MALAYSIA * shahrr@eng.ukm.my Abstract: - Ths paper presents the smulaton of ld-drven cavty for deep and shallow flow usng the lattce Boltzmann method. The effect the Reynolds number on the flow pattern at aspect rato of 5,.5,.5, 4 was studed extensvely. The lattce Boltzmann method s one of the most recent smulaton technques based on molecular theory whch s found to be a very effcent numercal tool due to ts capablty to go deeper nto the partcle s doman, smulatng ther nteracton among groups of partcle and relatng the parameters back to macro condton. The comparson of the results was made wth CFD software, verson 6.. Excellent agreement was obtaned even relatvely wth coarse grds appled on lattce the Boltzmann smulaton scheme. Keywords: - Lattce Boltzmann method, dstrbuton functon, mcroscopc velocty, tme relaxaton, ld-drven cavty flow Introducton The lattce Boltzmann method (), a numercal method evolved from mathematcal statstcal approach has been well accepted as an alternatve numercal scheme n computatonal flud dynamcs feld. In comparson wth other numercal schemes, s a bottom up approach, derves the Naver- Stokes uaton from statstcal behavor of partcles dynamcs. The magnary propagaton and collson processes of flud partcles are reconstructed n the formulaton of scheme. These processes are represented by the evoluton of partcle dstrbuton functon, f ( x, t) whch descrbes the statstcal populaton of partcles at locaton x and tme t. The advantages of nclude smple calculaton procedure, sutablty for parallel computaton, ease and robust handlng of multphase flow, complex geometres, nterfacal dynamcs and others []. A few standard benchmark problems have been smulated by and the results are found to agree well wth the correspondng Naver- Stokes solutons [2]. The ld-drven cavty flow s one of the most mportant benchmarks for new numercal method to be developed. It represents the flow of a rectangular or square geometry where the flow s drven by a tangental moton wth constant velocty of a sngle ld, representng the Drchlet boundary condtons. Moreover, the drven-cavty flow exhbts a number of nterestng physcal features [3]. However, on the other hand, the smulaton of ld-drven flow nsde cavty s scarcely performed. Only a few studes have dealt wth the nfluence of the heght of cavty on the flow nsde. In comparson wth the flow nsde a square cavty, one new parameter has to be taken nto account, the aspect rato of the cavty, K = W/H, W beng the cavty wdth and H the depth. As shown by Hou et al. [3] and Shen et al. [4], the physcs of drven flow n a cavty ndcate the exstence of crtcal aspect rato at whch the corner eddes merge and form a prmary eddy. The dependence of the vortex structure on the aspect rato at dfferent Reynolds numbers was nvestgated by Pan and Acrvos [5]. A flow vsualzaton experment was conducted to nvestgate the effect of nerta force n the flow structure n the Reynolds number of range 2 Re 4. Due to expermental lmtaton, the shallow cavty flow should consst essentally of a ISSN: ISBN:
2 Proceedngs of the 3th WSEAS Internatonal Conference on APPLIED MATHEMATICS (MATH'8) sngle nvscd core of unform vortcty, whle the vortex structure below the prmary vortex n the deep cavty could not be captured. Therefore, the objectves of the present artcle are to extend the smulatons for cavty flows and to nterpret the results wth partcular focus on the structure of prmary and secondary eddes. 2 Numercal Method Hstorcally, s the logcal development of lattce gas automata (LGA) method [6]. Lke n LGA, the physcal space s dscretzed nto unform lattce nodes. Every node n the network s then connected wth ts neghbours through a number of lattce veloctes to be determned through the model chosen. The lattce Boltzmann uaton s gven by: f t + c f,α =Ω( f ) () x α where f s the dstrbuton functon for partcles wth velocty c at poston x and tme t. Equaton () conssts of two parts; propagaton (left-hand sde) whch refers to the propagaton of dstrbuton functon to the next node n the drecton of ts Ω (rght-hand sde) whch represent the collson of the partcle dstrbuton functon. In, magntude of c s set up so that n each tme step t, every dstrbuton functon propagates n a dstance of lattce nodes spacng x. Ths wll ensure that dstrbuton functon arrves exactly at the lattce nodes after t and colldes smultaneously. There are a few versons of collson operator probable velocty, and collson ( f) ( f) Ω publshed n lterature. However, the most well accepted verson due to ts smplcty and effcency was the Bhatnagar-Gross-Crook (BGK) collson model [7]. The uaton that represents ths model s gven by: ( x, c, t) f ( x, c, t) f Ω ( f) = (2) τ where f s the ulbrum dstrbuton functon and τ s the tme to reach ulbrum condton durng collson process and s called relaxaton tme. Equaton (2) of the BGK collson model descrbes that τ of non-ulbrum dstrbuton relaxes to ulbrum state wthn tme τ on every collson process. By replacng the BGK collson model nto the Boltzmann uaton, the BGK Boltzmann Equaton s obtaned, f f f f + cα = (3) t xα τ The general form of the lattce velocty model s expressed as DnQm where D represents spatal dmenson and Q s the number of connecton (lattce velocty) at every node. In ths paper, the nne-mcroscopc velocty or nne-bt model (D2Q) s used. The lattce geometry s shown n Fg. c 7 c 4 c 3 c 6 c c 2 c c 8 c 5 Fg. (D2Q) lattce geometry The ulbrum dstrbuton functon of the nne-bt model s [8]: f = ρω + 3( c u) + ( c u) u (4) 2 2 The weghts are: ω 4,, = ω 2~5 = ω 6~ = (5) 36 and the mcroscopc velocty components are: c = (6) The macroscopc quanttes can be calculated from: ρ (7) = = f u cf f (8) Through the Chapman-Enskog procedure, the ncompressble Naver-Stokes uatons derved from the ncompressble []: u = () u 2τ + u u= P+ 2 u () t 6 By comparng the momentum uaton obtaned from the second Newton s law wth uaton (), the relaton between tme relaxaton τ and shear vscosty ν s obtaned as follow: τ = 3 ν + 2 () Equaton () descrbes that the value of τ must be kept hgher than.5 n order to avod negatve value of knematc vscosty. Ths lmts our smulaton to low value of Reynolds number. However ths lmtaton can be partly solved by usng hgh number of nodes but wll lead to long computatonal tme. Ths s to be addressed n our future works. ISSN: ISBN:
3 Proceedngs of the 3th WSEAS Internatonal Conference on APPLIED MATHEMATICS (MATH'8) 3 Results and Dscusson Ths secton appled the lattce Boltzmann scheme dscussed n the preceedng secton to smulate the ld-drven cavty flow at dfferent Reynolds number. 3. Code Valdaton For the purpose of code valdaton, we carred out the smulaton of ld-drven flow n a square cavty, K = and compared wth benchmark result []. Ld-drven flow n a square cavty s well-known as a standard test case for the numercal schemes of flud flows. Supposed that the ld s located at the top boundary y = H, and moves wth a constant speed, U from left to rght. The Reynolds number for the system s gven by: UL Re = (2) ν Smulatons were done at Re= usng grd sze 2 2. dstrbuton are n good agreement wth the prevous work of Gha et. al []. 3.2 Deep Cavty Flow was carred out to study flud flow over a deep cavty wth aspect rato K = 5 and.5. D2Q model wth the BGK collson model used. Shown n Fg. 3 are the streamlnes for cavty of aspect rato.5. The three fgures correspond to Re=,4 and respectvely. These smulaton employed 2 nodes. The smulatons were consdered to have reached steady state when the r.m.s. change n horzontal and vertcal velocty decreased to magntude 6 or less. (a) (b) (c) Fg.3 Streamlne at steady state condton (a) Re = (b) Re = 4 and (c) Re = It may be seen from Fg. 3 that, for low Reynolds number smulaton Re=, a counter-rotatng vortex s formed below the movng ld. As the Reynolds number ncreases ( Re= 4), the center of the prmary eddy begns to move downwards wth respect to the top ld. For the case of hgh Reynolds number ( Re=), the prmary eddy s formed at the center of the geometry. The secondary vortex can also be clearly seen at low Raylegh number and ntally formed at the lower rght of the geometry. As we ncrease the Reynolds number, ths vortex shft to the center left of the geometry test case. Ths s perhaps due to the effect of vscous effect produced by the prmary vortex. Fg.2 (top) Streamlne at steady state condton. (bottom) Velocty profle at md-heght (y-velocty) and md-wdth (x-velocty) of cavty Fg. 2 shows plot of stream functon and velocty profle for the Reynolds numbers consdered. It s apparent that the flow structure and velocty 3.3 Shallow Cavty Flow In ths secton numercal soluton of lattce Boltzmann method for rectangular shallow cavtes s present. Fg.4 shows that for aspect rato D =.5 (a),(b) and (c), the prmary vortex center descend to the ccnter of cavty as the Reynolds number ncreased. Two secondary vortces also can be ISSN: ISBN:
4 Proceedngs of the 3th WSEAS Internatonal Conference on APPLIED MATHEMATICS (MATH'8) clearly seen at the lower corners of the cavty wth left vortex s bgger than the rght vortex. (a) (b) (c) Fg.4 Steamlne at steady state condton for K =.5 (a) Re = (b) Re = 4 and (c) Re = For aspect rato D = 4 (a),(b) and (c), as Reynolds number s ncreased, Fg.5 shows that the prmary vortex began to splt nto two vortces. No any secondary vortces n the corners has been notced. Re = Re = 4 Re = Fg.6 Comparson between and for K = 5 (a) (b Re = Re = 4 (c) Fg.5 Steamlne at steady state condton for K = 4 (a) Re = (b) Re = 4 and (c) Re = Fgs. 6, 7, 8 and show the comparson of centerlne velocty profle obtaned from and. In ths comparson, K = 5,.5,.5 and 4 for Re =, 4 and has been selected. It may be observed that the present results are n excellent agreement wth the results for every tested case. Re = Fg.7 Comparson between and for K = Re = Re = Re = Fg.8 Comparson between and for K =.5 ISSN: ISBN:
5 Proceedngs of the 3th WSEAS Internatonal Conference on APPLIED MATHEMATICS (MATH'8) Re = Re = Re = Fg. Comparson between and for K = 4 4 Concluson The flow structure of a two-dmensonal vscous flow n ld-drven deep cavtes has been studed numercally usng lattce Boltzmann scheme. The dynamc of the prmary and secondary vortces are well captured by these smulatons. Whenever comparson s possble, the present results are found to be n good agreement wth the results obtaned by CFD software, verson 6.. The effects of boundary condtons n the lateral drectons (threedmensonal) effect wll be consdered n our future study. Acknowledgments The authors would lke to acknowledge the Mnstry of Scence, Technology and Innovaton, Malaysa for sponsorng ths work. References: [] S. Chen, and G. Doolen,, Lattce Boltzmann Method for Flud Flows, Annual Revew of Flud Mechancs, Vol. 3, 8, pp [2] G. McNamaea, and B. Alder,, Analyss of Lattce Boltzmann Treatment of Hydrodynamcs., Physca A, Vol. 4, 3, pp [3] S. Hou, Q. Zou, S. Chen, G. Doolen and A. C. Cogley, Smulaton of cavty flow by the lattce Boltzmann Method., Journal of Computatonal Physcs, Vol. 8, 5, pp [4] C. Shen, and J. M. Floryan, Low Reynolds Number Flows over Cavtes, Physcs of Fluds, 85, Vol. 28, pp.-22. [5] F. Pan and A. Acrvos, Steady Flows n Rectangular Cavtes., Journal of Flud Mechancs, 67, Vol. 28, pp [6] U. Frsh, B. Hasslacher and Y. Pomeau, Lattce gas automata for the Naver-Stokes uaton, Physcs Rev. Letters, vol. 56, 86, pp [7] P. L. Bhatnagar, E. P. Gross and M. Krook, A Model for Collson Processes n Gasses, Phys. Rev. Vol. 4, 54, pp [8] C.S. Nor Azwad and T. Tanahash, Threedmensonal thermal lattce-boltzmann smulaton of natural convecton n a cubc cavty, Intl. J. Modern Physcs B, Vol. 2, 27, pp [] L. S. Luo and X. He, Lattce Boltzmann model for the ncompressble Naver-Stokes uaton, J. Stat. Phys. Vol. 88, 7, pp [] U. Gha, K. N. Gha, and C. T. Shn, Hgh Re Soluton for Incompressble Flow usng the Naver-Stokes Equatons and a Multgrd Method, Journal of Computatonal Physcs, Vol. 48, 82, pp [] C. S. Azwad and T. Tanahash, Smplfed Thermal Lattce Boltzmann n Incompressble Lmt., Internatonal Journal of Modern Physcs B, Vol. 2, 26, pp [2] DP. Zegler, Boundary Condtons For Lattce Boltzmann Smulatons. J. Stat Phys. Vol. 7, 3, pp [3] D.V. Patl, K.N. Lakshmsha and B. Rogg. Lattce Boltzmann smulaton of ld-drven flow n deep cavtes, Journal of Computers & Fluds, Vol. 35, 26, No., pp [4] M. Cheng, K.C. Hung, Vortex structure of steady flow n a rectangular cavty, Journal of Computers & Fluds, Vol. 35, 26, No., pp ISSN: ISBN:
Simulation of Incompressible Flows in Two-Sided Lid-Driven Square Cavities. Part II - LBM
Perumal & Dass CFD Letters Vol. 2(1) 2010 www.cfdl.ssres.net 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
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 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 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 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 informationMesoscale Numerical Method for Prediction of Thermal Fluid Flow through Porous Media
Nor Azwad Che Sdk, Mohd Irwan Mohd Azm Mesoscale Numercal Method for Predcton of Thermal Flud Flow through Porous Meda NOR AZWADI CHE SIDIK and MOHD IRWAN MOHD AZMI Department of Thermofluds Unverst Teknolog
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 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 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 information(IREMOS) On New Solutions for Non-Newtonian Visco-Elastic Fluid in Pipe by N. Moallemi, I. Shafieenejad, H. Davari, A. Fata
ISSN 974-98 Vol. 4 N. 4 August Internatonal Revew on Modellng and Smulatons Contents: (contnued from Part B) (IREMOS) PART C Emprcal Models for the Correlaton of Global Solar Radaton Under Malaysa Envronment
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationSimulation of Flow Pattern in Open Channels with Sudden Expansions
Research Journal of Appled Scences, Engneerng and Technology 4(19): 3852-3857, 2012 ISSN: 2040-7467 Maxwell Scentfc Organzaton, 2012 Submtted: May 11, 2012 Accepted: June 01, 2012 Publshed: October 01,
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 informationFlow equations To simulate the flow, the Navier-Stokes system that includes continuity and momentum equations is solved
Smulaton of nose generaton and propagaton caused by the turbulent flow around bluff bodes Zamotn Krll e-mal: krart@gmal.com, cq: 958886 Summary Accurate predctons of nose generaton and spread n turbulent
More 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 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 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 informationCONTROLLED FLOW SIMULATION USING SPH METHOD
HERI COADA AIR FORCE ACADEMY ROMAIA ITERATIOAL COFERECE of SCIETIFIC PAPER AFASES 01 Brasov, 4-6 May 01 GEERAL M.R. STEFAIK ARMED FORCES ACADEMY SLOVAK REPUBLIC COTROLLED FLOW SIMULATIO USIG SPH METHOD
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 informationLattice Boltzmann simulation of nucleate boiling in micro-pillar structured surface
Proceedngs of the Asan Conference on Thermal Scences 017, 1st ACTS March 6-30, 017, Jeju Island, Korea ACTS-P00545 Lattce Boltzmann smulaton of nucleate bolng n mcro-pllar structured surface Png Zhou,
More informationA NUMERICAL COMPARISON OF LANGRANGE AND KANE S METHODS OF AN ARM SEGMENT
Internatonal Conference Mathematcal and Computatonal ology 0 Internatonal Journal of Modern Physcs: Conference Seres Vol. 9 0 68 75 World Scentfc Publshng Company DOI: 0.4/S009450059 A NUMERICAL COMPARISON
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 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 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 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 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 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 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 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 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 informationTurbulence classification of load data by the frequency and severity of wind gusts. Oscar Moñux, DEWI GmbH Kevin Bleibler, DEWI GmbH
Turbulence classfcaton of load data by the frequency and severty of wnd gusts Introducton Oscar Moñux, DEWI GmbH Kevn Blebler, DEWI GmbH Durng the wnd turbne developng process, one of the most mportant
More 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 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 NEW FILTERED DYNAMIC SUBGRID-SCALE MODEL FOR LARGE EDDY SIMULATION OF INDOOR AIRFLOW
A NEW FILTERED DYNAMIC SUBGRID-SCALE MODEL FOR LARGE EDDY SIMULATION OF INDOOR AIRFLOW We Zhang and Qngyan Chen Buldng Technology Program Massachusetts Insttute of Technology 77 Mass. Ave., Cambrdge, MA
More 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 informationTurbulence and its Modelling
School of Mechancal Aerospace and Cvl Engneerng 3rd Year Flud Mechancs Introducton In earler lectures we have consdered how flow nstabltes develop, and noted that above some crtcal Reynolds number flows
More 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 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 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 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 informationSignificance of Dirichlet Series Solution for a Boundary Value Problem
IOSR Journal of Engneerng (IOSRJEN) ISSN (e): 5-3 ISSN (p): 78-879 Vol. 6 Issue 6(June. 6) V PP 8-6 www.osrjen.org Sgnfcance of Drchlet Seres Soluton for a Boundary Value Problem Achala L. Nargund* and
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 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 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 informationSimulation Study on Characteristics of the Vortex Structure in Human Mouth-throat Model in Cyclic Respiratory Pattern
Avalable onlne at www.scencedrect.com Energy Proceda 17 (2012 ) 1786 1792 2012 Internatonal Conference on Future Electrcal Power and Energy Systems Smulaton Study on Characterstcs of the Vortex Structure
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 informationTHE NEAR-WALL INFLUENCE ON THE FLOW AROUND A SINGLE SQUARE CYLINDER.
THE NEAR-WALL INFLUENCE ON THE FLOW AROUND A SINGLE SQUARE CYLINDER. Campregher, Rubens Faculty of Mechancal Engneerng, FEMEC Federal Unversty of Uberlânda, UFU 38400-902 Uberlânda - Brazl campregher@mecanca.ufu.br
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 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 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 informationPublication 2006/01. Transport Equations in Incompressible. Lars Davidson
Publcaton 2006/01 Transport Equatons n Incompressble URANS and LES Lars Davdson Dvson of Flud Dynamcs Department of Appled Mechancs Chalmers Unversty of Technology Göteborg, Sweden, May 2006 Transport
More 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 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 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 informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
More informationFUZZY FINITE ELEMENT METHOD
FUZZY FINITE ELEMENT METHOD RELIABILITY TRUCTURE ANALYI UING PROBABILITY 3.. Maxmum Normal tress Internal force s the shear force, V has a magntude equal to the load P and bendng moment, M. Bendng moments
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 informationComparative Studies of Law of Conservation of Energy. and Law Clusters of Conservation of Generalized Energy
Comparatve Studes of Law of Conservaton of Energy and Law Clusters of Conservaton of Generalzed Energy No.3 of Comparatve Physcs Seres Papers Fu Yuhua (CNOOC Research Insttute, E-mal:fuyh1945@sna.com)
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 informationDETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM
Ganj, Z. Z., et al.: Determnaton of Temperature Dstrbuton for S111 DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM by Davood Domr GANJI
More 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 information1. Why turbulence occur? Hydrodynamic Instability. Hydrodynamic Instability. Centrifugal Instability: Rayleigh-Benard Instability:
. Why turbulence occur? Hydrodynamc Instablty Hydrodynamc Instablty T Centrfugal Instablty: Ω Raylegh-Benard Instablty: Drvng force: centrfugal force Drvng force: buoyancy flud Dampng force: vscous dsspaton
More informationAssignment 5. Simulation for Logistics. Monti, N.E. Yunita, T.
Assgnment 5 Smulaton for Logstcs Mont, N.E. Yunta, T. November 26, 2007 1. Smulaton Desgn The frst objectve of ths assgnment s to derve a 90% two-sded Confdence Interval (CI) for the average watng tme
More informationNormally, in one phase reservoir simulation we would deal with one of the following fluid systems:
TPG4160 Reservor Smulaton 2017 page 1 of 9 ONE-DIMENSIONAL, ONE-PHASE RESERVOIR SIMULATION Flud systems The term sngle phase apples to any system wth only one phase present n the reservor In some cases
More informationNUMERICAL SIMULATION OF FLOW OVER STEPPED SPILLWAYS
ISSN: 345-3109 RCEE Research n Cvl and Envronmental Engneerng www.rcee.com Research n Cvl and Envronmental Engneerng 014 (04) 190-198 NUMERICAL SIMULATION OF FLOW OVER STEPPED SPILLWAYS Rasoul Daneshfaraz
More informationMagnetic Field Around The New 400kV OH Power Transmission Lines In Libya
ECENT ADVANCES n ENEGY & ENVIONMENT Magnetc Feld Around The New kv OH Power Transmsson Lnes In Lbya JAMAL M. EHTAIBA * SAYEH M. ELHABASHI ** Organzaton for Development of Admnstratve Centers, ODAC MISUATA
More information829. An adaptive method for inertia force identification in cantilever under moving mass
89. An adaptve method for nerta force dentfcaton n cantlever under movng mass Qang Chen 1, Mnzhuo Wang, Hao Yan 3, Haonan Ye 4, Guola Yang 5 1,, 3, 4 Department of Control and System Engneerng, Nanng Unversty,
More 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 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 informationNumerical simulation of flow past circular duct
Water Scence and Engneerng, 010, 3(): 08-16 do:10.388/.ssn.1674-370.010.0.009 http://.aterournal.cn e-mal: se008@vp.163.com Numercal smulaton of flo past crcular duct Ze-gao YIN* 1, Xan-e Cao 1, Hong-da
More information11. Dynamics in Rotating Frames of Reference
Unversty of Rhode Island DgtalCommons@URI Classcal Dynamcs Physcs Course Materals 2015 11. Dynamcs n Rotatng Frames of Reference Gerhard Müller Unversty of Rhode Island, gmuller@ur.edu Creatve Commons
More informationTemperature. Chapter Heat Engine
Chapter 3 Temperature In prevous chapters of these notes we ntroduced the Prncple of Maxmum ntropy as a technque for estmatng probablty dstrbutons consstent wth constrants. In Chapter 9 we dscussed the
More informationχ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body
Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown
More informationElectron-Impact Double Ionization of the H 2
I R A P 6(), Dec. 5, pp. 9- Electron-Impact Double Ionzaton of the H olecule Internatonal Scence Press ISSN: 9-59 Electron-Impact Double Ionzaton of the H olecule. S. PINDZOLA AND J. COLGAN Department
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 informationEnergy configuration optimization of submerged propeller in oxidation ditch based on CFD
IOP Conference Seres: Earth and Envronmental Scence Energy confguraton optmzaton of submerged propeller n oxdaton dtch based on CFD To cte ths artcle: S Y Wu et al 01 IOP Conf. Ser.: Earth Envron. Sc.
More 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 informationInteraction of a pair of horizontally aligned bubbles in gravity field
Interacton of a par of horzontally algned bubbles n gravty feld Han JIAO 1 ; Dongyan SHI 2 ; Zhka Wang 3 ; Hongqun LI 4 1 Harbn Engneerng Unversty, Chna 2 Harbn Engneerng Unversty, Chna 3 Harbn Engneerng
More informationChapter 3. r r. Position, Velocity, and Acceleration Revisited
Chapter 3 Poston, Velocty, and Acceleraton Revsted The poston vector of a partcle s a vector drawn from the orgn to the locaton of the partcle. In two dmensons: r = x ˆ+ yj ˆ (1) The dsplacement vector
More informationREAL TIME AIRFLOW SIMULATION IN BUILDINGS
REAL TIME AIRFLOW SIMULATION IN BUILDINGS Wangda Zuo, and Qngyan (Yan) Chen School of Mechancal Engneerng, Purdue Unversty, West Lafayette, USA ABSTRACT Real tme flow smulaton s crucal n emergency management
More informationLecture 5.8 Flux Vector Splitting
Lecture 5.8 Flux Vector Splttng 1 Flux Vector Splttng The vector E n (5.7.) can be rewrtten as E = AU (5.8.1) (wth A as gven n (5.7.4) or (5.7.6) ) whenever, the equaton of state s of the separable form
More informationComputational investigation of the external excitation frequency effect on liquid sloshing phenomenon
Computatonal nvestgaton of the external exctaton frequency effect on lqud sloshng phenomenon Abdallah Bouabd *, Zed Drss, Laboratory of Electro-Mechanc Systems (LASEM) Natonal School of Engneers of Sfax
More informationDynaflowC.F.D. conference
DynaflowC.F.D. conference Lattce Boltzmann Smulatons and Applcaton to Multphase Flow Sacha Jelc, Senor C.F.D. Thermal Engneer Exa GmbH January 13th 2011 Agenda Company Overvew Current Applcatons Overvew
More informationx = , so that calculated
Stat 4, secton Sngle Factor ANOVA notes by Tm Plachowsk n chapter 8 we conducted hypothess tests n whch we compared a sngle sample s mean or proporton to some hypotheszed value Chapter 9 expanded ths to
More informationAnalysis of Unsteady Aerodynamics of a Car Model with Radiator in Dynamic Pitching Motion using LS-DYNA
Analyss of Unsteady Aerodynamcs of a Car Model wth Radator n Dynamc Ptchng Moton usng LS-DYNA Yusuke Nakae 1, Jro Takamtsu 1, Hrosh Tanaka 1, Tsuyosh Yasuk 1 1 Toyota Motor Corporaton 1 Introducton Recently,
More information