Investigation of a New Monte Carlo Method for the Transitional Gas Flow
|
|
- Russell McLaughlin
- 6 years ago
- Views:
Transcription
1 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 Monte Carlo method (DSMC) s well developed foarefed gas flow n transton flow regme when 0.0<Kn<. However such a smulaton for a complex 3D vacuum system s stll a challengng task because of the huge demand on the memory and long computatonal tme. On the other hand f Kn>0 the gas flow s free molecular and can be smulated by the Test Partcle Monte Carlo method (TPMC) wthout any problem even for a complex 3D vacuum system. In ths paper we wll nvestgate the approach to extend the TPMC to transton flow regme by consderng the collson between gas molecules as an nteracton between a probe molecule and the gas background. Recently ths collson mechansm has been mplemented nto ProVac3D a new TPMC smulaton program developed by KIT. The prelmnary smulaton result shows a correct nonlnear ncreasng of the gas flow. However there s stll a quanttatve dscrepancy wth the expermental data whch means further mprovement s needed. Keywords: Test partcle Monte Carlo smulaton molecular collson transtonal gas flow PACS: n 05.0.Ln INTRODUCTION The gas flow n a vacuum system s characterzed by the Knudsen number Kn=λ/d wth λ beng the mean free path of the molecules and d the sze of the vacuum chamber. For an ultra hgh vacuum or a hgh vacuum system the mean free path s much greater than the sze of the vacuum chamber so Kn>0. In ths case the gas molecules seldom collde wth each other and the gas flow s called as free molecular flow whch can be smulated by the Test Partcle Monte Carlo method (TPMC). If Kn<0.0 the gas can be usually consdered as contnuous flud and the flow s called as vscous flow whch can be smulated by Computatonal Flud Dynamcs (CFD). The gas flow regme when 0.0<Kn< s called as transton flow regme whch can be smulated by knetc theory or the Drect Smulaton Monte Carlo method (DSMC). However the smulaton of a complex 3D vacuum system wth DSMC or knetc theory s stll a challenge task because of the huge demand on the memory and long computatonal tme. Unfortunately the gas flows n some ITER applcatons for example the gas flow nsde the neutralzer of the Neutral Beam Inecton system (NBI) the gas flow before the cryogenc pump n pellet necton scenaro s wthn the transton flow regme. So that t s reasonable to explore a new smulaton Monte Carlo method for the transtonal gas flow. DEVELOPMENT OF PROVAC3D INCLUDING MOLECULAR COLLISIONS ProVac3D standng for 3D densty PROfle n the VACuum system s a TPMC smulaton program developed by the Karlsruhe Insttute of Technology (KIT). In our prevous works [-3] we have well descrbed ts man deas and successfully cross checked and used t n dfferent applcatons n free molecular flow regme. By recordng the tme of flght of molecules we can calculate the gas densty dstrbuton. The natural dea to extend ProVac3D nto transton flow regme s to nclude the collson between gas molecules by consderng approprate nteracton between the probe molecule and the gas background whch s descrbed n Fgure. Suppose P(x) r s the probablty that the probe molecule meets a target molecule of the gas background before t has passed a dstance x r t reads as r P( x) = exp( σ n( x) dx) ()
2 FIGURE. The collson between the probe molecule and target molecule. where n(x r ) s the densty of the background gas molecules and σ the collson cross secton. The ntegral s along the path of the probe molecule. From ths relaton one can derve Equaton (2) to determne the collson tme t of the probe molecule wth one of the molecule n background. t log( rnd) ( )( ) (2) n x y z v probe ubulk dt = σ 0 where v probe s the speed of the probe molecule u bulk the bulk speed of the background gas rnd a random number. Suppose τ s the possble collson tme wth the wall f no collson wth the molecule of the background s taken nto account. When t < τ the collson wth the target molecule wll happen at the poston In such a case the velocty hard sphere collson model. v wth ( ) 2 relatve bulk velocty u r v v (3) r of the probe molecule after the collson can be calculated for the frst attempt by the v = and r v = 2 x = x0 + v r t. r [ v e + ( v + v )] relatve x) = sn( θ ) cos( ϕ ) y) = sn( θ ) sn( ϕ ) z) = cos( θ ) v r beng the velocty of the randomly chosen target molecule superposed by v ( x) = v sn( θ ) cos( ϕ v ( y) = v sn( θ v ( z) = v cos( θ In Equatons (4) and (5) the azmuth angles of the target molecule θ ϕ ) are randomly chosen n full space ( however the azmuth angles of the probe molecule after the collson ( θ ϕ ) are randomly chosen n the half-space determned by theelatve velocty. On the other hand when t > τ the collson wth the wall wll happen. In such a case Cercgnan-Lamps boundary condtons have been ncluded [45]. )sn( ϕ z. x y (4) (5) 2
3 NORMALIZATION PROCEDURE AND CONVERGENCE Above mentoned collson mechansm has been recently mplemented nto ProVac3D. Naturally t can be used to the applcaton n whch a dlute gas speces s movng nsde the gas background of another gas speces [6]. However the emphass of the present work s to smulate the transton flow of sngle gas speces through a crcular tube connectng wth a dosng dome wth gas densty n and pumpng dome wth neglgble gas densty n2<<n as llustrated n Fgure 2. FIGURE 2. ProVac3D smulaton model. The gas flow s smulated as the functon of the gas densty n n the dosng dome by an teraton process n = n0 + k n( k = 023 L). (6) In order to start from a free molecular background n 0 wth k=0 must be small enough. For our smulatons n s chosen at least one magntude smaller than n 0 so that the collsons between the probe molecules can be always neglgble. For each teraton step k a hgh number of probe partcles representng n has been smulated. As sad the smulaton doman s dvded nto cells n ProVac3D smulaton [-3]. If s the cell ndex we can calculate the average tme of flght of all test partcles () ; here subscrpt k means the k th teraton step. Because even there are collsons the densty s stll proportonal to fts k () and nversely proportonal to the cell volume V () we can normalze to obtan the densty dstrbuton untl the k th teraton step s fnshed fts k n k ( ) = C ( fts0 ( ) n0 + fts ( ) n + L + ftsk ( ) n) / V ( ) (7) wth C beng the global normalzaton constant. Please note that the densty dstrbuton () s then used n Equaton (2) to determne the new collson tme n n k the next (k+) th teraton step. We have checked ths normalzaton procedure wth dfferent cell shapes cell szes molecular collson cross sectons and ncrease step n and proven that t s convergent. SIMULATION RESULTS AND COMPARISON WITH EXPERIMENTAL DATA The smulaton s carred out for a crcular tube of L=0.570 m and D=0.06 m. For ths crcular tube of L/D=9.75 t s not only hard to calculate wth knetc theory because the length-to-dameteato s too small but also hard to calculate wth DSMC because the length-to-dameteato s too large. The gas smulated s ntrogen at 5 C. The collson cross secton s calculated from the expermental data of the vscosty [7]. The experment was carred out n the TRANSFLOW test rg [8]. The smulaton doman s only the tube tself f we neglect the collsons nsde the dosng dome. In such a case C = 4 va wth v beng the thermal molecular speed and A the area of the tube cross secton. Consequently the gas flow rate Q can be determned by the transmsson probabltes w k n each teraton step 3
4 Q = C L k ( w0 n0 + w n + + w n Other macroscopc parameters such as bulk velocty and temperature etc. can be obtaned n the smlar way. In the smulaton n 0 s chosen as 0 8 correspondng to an ntal P= Pa and an ntal Kn=00. respectvely. n s chosen as 0 7. After k=9990 teraton steps the fnal n s 0 2 whch s 000 tmes bgger than the ntal one correspondng to a fnal P=3.98 Pa and a fnal Kn=0. respectvely. So the gas flow evolves from the free molecular flow to the begnnng of the transtonal flow. For each teraton step 0 5 test partcles are used whch means that one test partcle represents about ntrogen molecules. There are 5 meshes along the radal drecton and 40 meshes along the axal drecton of the tube. In total the smulaton doman s dvded nto 200 cells. The smulaton s fnshed roughly n 2 hours by a desktop PC wth a CPU at 2.67 GHz. ). (8) FIGURE 3. Comparson of expermental result and ProVac3D smulaton for a crcular tube of L/D=9.75. The ProVac3D smulaton result n Fgure 3 shows a correct nonlnear ncreasng of the gas flow as P ncreases but there s stll a quanttatve dscrepancy wth the expermental data. Frst of all we have neglected the pressure n the pumpng dome. From the experment we know that the pressure n the pumpng dome though always beng at least one order of magntude lower than that n the dosng dome for ths partcular tube wll ntroduce a small back streamng. Secondly the entrance effect or the collsons n the dosng dome may play a role. We wll mprove our model to nclude these effects. Addtonally the smulated relatve axal densty and the axal bulk velocty nsde the tube are shown n Fgure 4 and Fgure 5 respectvely. We can see that the results for n=0 8 and n=0 2 are dfferent whch s qualtatvely reasonable. Unfortunately we have nether the expermental data nor other numercal results to compare wth. FIGURE 4. Relatve densty nsde the tube for a crcular tube of L/D=
5 FIGURE 5. Bulk velocty n the axal drecton for a crcular tube of L/D=9.75. CONCLUSIONS We have expanded our ProVac3D Monte Carlo smulaton program to explore a new approach for smulatng the transtonal gas flow. In ths new approach the collsons between molecules are represented by the collsons of the probe partcle wth the target partcle n the background. One convergent normalzaton scheme s suggested. Wth ths procedure the gas flow through a crcular tube s smulated as a functon of the ncreasng molecular densty n the dosng dome. The prelmnary smulaton result shows the correct nonlnear ncreasng of the gas flow rate. However there s stll a quanttatve dscrepancy wth the expermental data whch means further mprovement s needed. Obvously ths new smulaton approach has several advantages. Frst of all lke usual TPMC the test partcle s smulated one by one. So we wll not need the huge memory lke DSMC whch makes the smulaton of a 3D problem possble. Actually ProVac3D has mplemented the enttes to model a complcated 3D vacuum system. The present work shows the computaton tme s much faster than DSMC. Secondly as well known DSMC s very hard to be parallelzed. However the parallelzaton of TPMC s straghtforward. Ths work for ProVac3D s n progress. Thrdly the gas flow s smulated by an teraton process n the suggested normalzaton scheme. Ths means that we can even nterrupt the smulaton process and later on contnue the smulaton of the further evoluton of the gas flow by usng the obtaned result as the background. Nevertheless the drawback of ths normalzaton scheme s also clear snce we always need to start from the free molecular flow and there are many smulaton steps needed to reach the transton flow regme. REFERENCES. X. Luo M. Dremel and Ch. Day ProVac3D and Applcaton to the Neutral Beam Inecton System of ITER n 26 th Internatonal Symposum on Rarfed Gas Dynamcs edted by T. Abe AIP Conference Proceedngs 084 Amercan Insttute of Physcs Melvlle NY 2009 pp M. Dremel Chr. Day S. Hanke and X. Luo Cryopump desgn development for the ITER Neutral Beam Inectors Fuson Engneerng and Desgn (2009). 3. X. Luo Chr. Day 3D Monte Carlo vacuum modelng of the neutral beam necton system of ITER Fuson Engneerng and Desgn n Press. 2. C. Cercgnan and M. Lamps Knetc models for gas-surface nteractons Trans. Theory Stat. Phys. 0-04(97) 5. F. Sharpov Applcaton of the Cercgnan-Lamps Scatterng Kernel to Channel Gas Flows n 22 nd Internatonal Symposum on Rarfed Gas Dynamcs edted by T. J. Bartel and M. A. Galls AIP Conference Proceedngs 585 Amercan Insttute of Physcs Melvlle NY 200 pp S. Longo and P. Domede A Monte Carlo model for seeded atomc flows n the transton regme J. Comput. Phys (2009). 7. W. A. Cole and W. A. Wakeham The vscosty of ntrogen oxygen and ther bnary mxtures n the lmt of zero densty J. Phys. Chem. Ref. Data (985). 8. S. Varouts S. Nars V. Hauer Chr. Day and D. Valougeorgs Computatonal and expermental study of gas flows through long channels of varous cross sectons n the whole range of the Knudsen number J. Vac. Sc. Technol. A (2009). 5
Numerical Simulation of Lid-Driven Cavity Flow Using the Lattice Boltzmann Method
Proceedngs of the 3th WSEAS Internatonal Conference on APPLIED MATHEMATICS (MATH'8) Numercal Smulaton of Ld-Drven Cavty Flow Usng the Lattce Boltzmann Method M.A. MUSSA, S. ABDULLAH *, C.S. NOR AZWADI
More 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 informationGaussian Mixture Models
Lab Gaussan Mxture Models Lab Objectve: Understand the formulaton of Gaussan Mxture Models (GMMs) and how to estmate GMM parameters. You ve already seen GMMs as the observaton dstrbuton n certan contnuous
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 informationAppendix B: Resampling Algorithms
407 Appendx B: Resamplng Algorthms A common problem of all partcle flters s the degeneracy of weghts, whch conssts of the unbounded ncrease of the varance of the mportance weghts ω [ ] of the partcles
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 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 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 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 informationMarkov Chain Monte Carlo Lecture 6
where (x 1,..., x N ) X N, N s called the populaton sze, f(x) f (x) for at least one {1, 2,..., N}, and those dfferent from f(x) are called the tral dstrbutons n terms of mportance samplng. Dfferent ways
More informationSUPPLEMENTARY INFORMATION
do: 0.08/nature09 I. Resonant absorpton of XUV pulses n Kr + usng the reduced densty matrx approach The quantum beats nvestgated n ths paper are the result of nterference between two exctaton paths of
More informationNeutral-Current Neutrino-Nucleus Inelastic Reactions for Core Collapse Supernovae
Neutral-Current Neutrno-Nucleus Inelastc Reactons for Core Collapse Supernovae A. Juodagalvs Teornės Fzkos r Astronomjos Insttutas, Lthuana E-mal: andrusj@tpa.lt J. M. Sampao Centro de Físca Nuclear da
More informationEPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski
EPR Paradox and the Physcal Meanng of an Experment n Quantum Mechancs Vesseln C Nonnsk vesselnnonnsk@verzonnet Abstract It s shown that there s one purely determnstc outcome when measurement s made on
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 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 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 informationThe optimal delay of the second test is therefore approximately 210 hours earlier than =2.
THE IEC 61508 FORMULAS 223 The optmal delay of the second test s therefore approxmately 210 hours earler than =2. 8.4 The IEC 61508 Formulas IEC 61508-6 provdes approxmaton formulas for the PF for smple
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 informationDesign and Optimization of Fuzzy Controller for Inverse Pendulum System Using Genetic Algorithm
Desgn and Optmzaton of Fuzzy Controller for Inverse Pendulum System Usng Genetc Algorthm H. Mehraban A. Ashoor Unversty of Tehran Unversty of Tehran h.mehraban@ece.ut.ac.r a.ashoor@ece.ut.ac.r Abstract:
More informationProblem Set 9 Solutions
Desgn and Analyss of Algorthms May 4, 2015 Massachusetts Insttute of Technology 6.046J/18.410J Profs. Erk Demane, Srn Devadas, and Nancy Lynch Problem Set 9 Solutons Problem Set 9 Solutons Ths problem
More informationResource Allocation with a Budget Constraint for Computing Independent Tasks in the Cloud
Resource Allocaton wth a Budget Constrant for Computng Independent Tasks n the Cloud Wemng Sh and Bo Hong School of Electrcal and Computer Engneerng Georga Insttute of Technology, USA 2nd IEEE Internatonal
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 informationPhysics for Scientists and Engineers. Chapter 9 Impulse and Momentum
Physcs or Scentsts and Engneers Chapter 9 Impulse and Momentum Sprng, 008 Ho Jung Pak Lnear Momentum Lnear momentum o an object o mass m movng wth a velocty v s dened to be p mv Momentum and lnear momentum
More informationThe non-negativity of probabilities and the collapse of state
The non-negatvty of probabltes and the collapse of state Slobodan Prvanovć Insttute of Physcs, P.O. Box 57, 11080 Belgrade, Serba Abstract The dynamcal equaton, beng the combnaton of Schrödnger and Louvlle
More informationLecture Note 3. Eshelby s Inclusion II
ME340B Elastcty of Mcroscopc Structures Stanford Unversty Wnter 004 Lecture Note 3. Eshelby s Incluson II Chrs Wenberger and We Ca c All rghts reserved January 6, 004 Contents 1 Incluson energy n an nfnte
More informationNumerical simulation of a binary gas flow inside a rotating cylinder
Journal of Mechancal Scence and Technology 3 (9) 848~86 Journal of Mechancal Scence and Technology wwwsprngerlnkcom/content/738-494x DOI 7/s6-8-- Numercal smulaton of a bnary gas flow nsde a rotatng cylnder
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 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 informationSupporting Information
Supportng Informaton The neural network f n Eq. 1 s gven by: f x l = ReLU W atom x l + b atom, 2 where ReLU s the element-wse rectfed lnear unt, 21.e., ReLUx = max0, x, W atom R d d s the weght matrx to
More 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 informationPhysics 114 Exam 2 Fall 2014 Solutions. Name:
Physcs 114 Exam Fall 014 Name: For gradng purposes (do not wrte here): Queston 1. 1... 3. 3. Problem Answer each of the followng questons. Ponts for each queston are ndcated n red. Unless otherwse ndcated,
More informationKernel Methods and SVMs Extension
Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general
More 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 informationNon-equilibrium structure of the vapor-liquid interface of a binary fluid
Non-equlbrum structure of the vapor-lqud nterface of a bnary flud Aldo Frezzott Dpartmento d Matematca del Poltecnco d Mlano - Pazza Leonardo da Vnc 32 - I233 Mlano - Italy Abstract. The evaporaton of
More informationmodeling of equilibrium and dynamic multi-component adsorption in a two-layered fixed bed for purification of hydrogen from methane reforming products
modelng of equlbrum and dynamc mult-component adsorpton n a two-layered fxed bed for purfcaton of hydrogen from methane reformng products Mohammad A. Ebrahm, Mahmood R. G. Arsalan, Shohreh Fatem * Laboratory
More informationVQ widely used in coding speech, image, and video
at Scalar quantzers are specal cases of vector quantzers (VQ): they are constraned to look at one sample at a tme (memoryless) VQ does not have such constrant better RD perfomance expected Source codng
More informationThe Study of Teaching-learning-based Optimization Algorithm
Advanced Scence and Technology Letters Vol. (AST 06), pp.05- http://dx.do.org/0.57/astl.06. The Study of Teachng-learnng-based Optmzaton Algorthm u Sun, Yan fu, Lele Kong, Haolang Q,, Helongang Insttute
More informationWilbur and Ague 4 WILBUR AND AGUE; APPENDIX DR1. Two-dimensional chemical maps as well as chemical profiles were done at 15 kv using
DR2006139 Wlbur and Ague 4 WILBUR AND AGUE; APPENDIX DR1 MINERAL ANALYSES Two-dmensonal chemcal maps as well as chemcal profles were done at 15 kv usng the JEOL JXA-8600 electron mcroprobe at Yale Unversty
More informationMarkov Chain Monte Carlo (MCMC), Gibbs Sampling, Metropolis Algorithms, and Simulated Annealing Bioinformatics Course Supplement
Markov Chan Monte Carlo MCMC, Gbbs Samplng, Metropols Algorthms, and Smulated Annealng 2001 Bonformatcs Course Supplement SNU Bontellgence Lab http://bsnuackr/ Outlne! Markov Chan Monte Carlo MCMC! Metropols-Hastngs
More informationImplicit Integration Henyey Method
Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure
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 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 informationRobert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations
Quantum Physcs 量 理 Robert Esberg Second edton CH 09 Multelectron atoms ground states and x-ray exctatons 9-01 By gong through the procedure ndcated n the text, develop the tme-ndependent Schroednger equaton
More informationONE-DIMENSIONAL COLLISIONS
Purpose Theory ONE-DIMENSIONAL COLLISIONS a. To very the law o conservaton o lnear momentum n one-dmensonal collsons. b. To study conservaton o energy and lnear momentum n both elastc and nelastc onedmensonal
More informationPhysics 607 Exam 1. ( ) = 1, Γ( z +1) = zγ( z) x n e x2 dx = 1. e x2
Physcs 607 Exam 1 Please be well-organzed, and show all sgnfcant steps clearly n all problems. You are graded on your wor, so please do not just wrte down answers wth no explanaton! Do all your wor on
More informationPY2101 Classical Mechanics Dr. Síle Nic Chormaic, Room 215 D Kane Bldg
PY2101 Classcal Mechancs Dr. Síle Nc Chormac, Room 215 D Kane Bldg s.ncchormac@ucc.e Lectures stll some ssues to resolve. Slots shared between PY2101 and PY2104. Hope to have t fnalsed by tomorrow. Mondays
More information1 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations
Physcs 171/271 -Davd Klenfeld - Fall 2005 (revsed Wnter 2011) 1 Dervaton of Rate Equatons from Sngle-Cell Conductance (Hodgkn-Huxley-lke) Equatons We consder a network of many neurons, each of whch obeys
More informationThree-Phase Distillation in Packed Towers: Short-Cut Modelling and Parameter Tuning
European Symposum on Computer Arded Aded Process Engneerng 15 L. Pugjaner and A. Espuña (Edtors) 2005 Elsever Scence B.V. All rghts reserved. Three-Phase Dstllaton n Packed Towers: Short-Cut Modellng and
More informationBrownian-Dynamics Simulation of Colloidal Suspensions with Kob-Andersen Type Lennard-Jones Potentials 1
Brownan-Dynamcs Smulaton of Collodal Suspensons wth Kob-Andersen Type Lennard-Jones Potentals 1 Yuto KIMURA 2 and Mcho TOKUYAMA 3 Summary Extensve Brownan-dynamcs smulatons of bnary collodal suspenton
More informationChapter - 2. Distribution System Power Flow Analysis
Chapter - 2 Dstrbuton System Power Flow Analyss CHAPTER - 2 Radal Dstrbuton System Load Flow 2.1 Introducton Load flow s an mportant tool [66] for analyzng electrcal power system network performance. Load
More informationA Comparative Investigation into Aerodynamic Performances of Two Set Finned Bodies with Circular and Non Circular Cross Sections
A Comparatve Investgaton nto Aerodynamc Performances of Two Set Fnned Bodes wth Crcular and Non Crcular Cross Sectons MAHMOUD MANI and SHADI MAHJOOB Aerospace Engneerng Department Amrkabr Unversty of Technology
More 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 informationA Computational Viewpoint on Classical Density Functional Theory
A Computatonal Vewpont on Classcal Densty Functonal Theory Matthew Knepley and Drk Gllespe Computaton Insttute Unversty of Chcago Department of Molecular Bology and Physology Rush Unversty Medcal Center
More informationAn Improved Model for the Droplet Size Distribution in Sprays Developed From the Principle of Entropy Generation maximization
ILASS Amercas, 9 th Annual Conference on Lqud Atomzaton and Spray Systems, oronto, Canada, May 6 An Improved Model for the Droplet Sze Dstrbuton n Sprays Developed From the Prncple of Entropy Generaton
More informationConsistency & Convergence
/9/007 CHE 374 Computatonal Methods n Engneerng Ordnary Dfferental Equatons Consstency, Convergence, Stablty, Stffness and Adaptve and Implct Methods ODE s n MATLAB, etc Consstency & Convergence Consstency
More informationMoments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is.
Moments of Inerta Suppose a body s movng on a crcular path wth constant speed Let s consder two quanttes: the body s angular momentum L about the center of the crcle, and ts knetc energy T How are these
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 informationMODELING TRAFFIC LIGHTS IN INTERSECTION USING PETRI NETS
The 3 rd Internatonal Conference on Mathematcs and Statstcs (ICoMS-3) Insttut Pertanan Bogor, Indonesa, 5-6 August 28 MODELING TRAFFIC LIGHTS IN INTERSECTION USING PETRI NETS 1 Deky Adzkya and 2 Subono
More informationConstitutive Modelling of Superplastic AA-5083
TECHNISCHE MECHANIK, 3, -5, (01, 1-6 submtted: September 19, 011 Consttutve Modellng of Superplastc AA-5083 G. Gulano In ths study a fast procedure for determnng the constants of superplastc 5083 Al alloy
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 informationLINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity
LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 30 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 2 Remedes for multcollnearty Varous technques have
More informationPhysics 5153 Classical Mechanics. Principle of Virtual Work-1
P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal
More informationLinear Momentum. Center of Mass.
Lecture 6 Chapter 9 Physcs I 03.3.04 Lnear omentum. Center of ass. Course webste: http://faculty.uml.edu/ndry_danylov/teachng/physcsi Lecture Capture: http://echo360.uml.edu/danylov03/physcssprng.html
More informationEffect of loading frequency on the settlement of granular layer
Effect of loadng frequency on the settlement of granular layer Akko KONO Ralway Techncal Research Insttute, Japan Takash Matsushma Tsukuba Unversty, Japan ABSTRACT: Cyclc loadng tests were performed both
More informationComputational issues surrounding the management of an ecological food web
Computatonal ssues surroundng the management of an ecologcal food web Wllam J M Probert, Eve McDonald-Madden, Nathale Peyrard, Régs Sabbadn AIGM 12, ECAI2012 Montpeller, France Ratonale Ecology has many
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 informationA Particle Filter Algorithm based on Mixing of Prior probability density and UKF as Generate Importance Function
Advanced Scence and Technology Letters, pp.83-87 http://dx.do.org/10.14257/astl.2014.53.20 A Partcle Flter Algorthm based on Mxng of Pror probablty densty and UKF as Generate Importance Functon Lu Lu 1,1,
More information1 American Institute of Aeronautics and Astronautics Approved for public release; distribution unlimited
4nd AIAA Aerospace Scences Meetng and Exhbt 5-8 January 004, Reno, Nevada AIAA 004-35 4 nd AIAA Aerospace Scences Meetng and Exhbt Reno, Nevada 5-8 January 004 DEVELOPMENT OF A TWO-WAY COUPLED MODEL FOR
More informationJournal of Biomechanical Science and Engineering
Scence and Engneerng Smulaton Study on Effects of Hematocrt on Blood Flow Propertes Usng Partcle Method * Ken-ch TSUBOTA, ** Shgeo WADA *** and Takam YAMAGUCHI ** **Department of Boengneerng and Robotcs,
More informationSimulation and experiment of the effect of clearance of impeller wear-rings on the performance of centrifugal pump
IOP Conference Seres: Earth and Envronmental Scence Smulaton and experment of the effect of clearance of mpeller wear-rngs on the performance of centrfugal pump To cte ths artcle: S X Chen et al 01 IOP
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationGeneralized Linear Methods
Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set
More informationINTERROGATING THE FLOW BEHAVIOUR IN A NOVEL MAGNETIC DESICCANT VENTILATION SYSTEM USING COMPUTATIONAL FLUID DYNAMICS (CFD)
INTERROGATING THE FLOW BEHAVIOUR IN A NOVEL MAGNETIC DESICCANT VENTILATION SYSTEM USING COMPUTATIONAL FLUID DYNAMICS (CFD) Auwal Dodo*, Valente Hernandez-Perez, Je Zhu and Saffa Rffat Faculty of Engneerng,
More information9 Derivation of Rate Equations from Single-Cell Conductance (Hodgkin-Huxley-like) Equations
Physcs 171/271 - Chapter 9R -Davd Klenfeld - Fall 2005 9 Dervaton of Rate Equatons from Sngle-Cell Conductance (Hodgkn-Huxley-lke) Equatons We consder a network of many neurons, each of whch obeys a set
More informationIntroduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:
CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and
More informationConvergence of random processes
DS-GA 12 Lecture notes 6 Fall 216 Convergence of random processes 1 Introducton In these notes we study convergence of dscrete random processes. Ths allows to characterze phenomena such as the law of large
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 informationSecond Order Analysis
Second Order Analyss In the prevous classes we looked at a method that determnes the load correspondng to a state of bfurcaton equlbrum of a perfect frame by egenvalye analyss The system was assumed to
More information6.3.4 Modified Euler s method of integration
6.3.4 Modfed Euler s method of ntegraton Before dscussng the applcaton of Euler s method for solvng the swng equatons, let us frst revew the basc Euler s method of numercal ntegraton. Let the general from
More informationNote 10. Modeling and Simulation of Dynamic Systems
Lecture Notes of ME 475: Introducton to Mechatroncs Note 0 Modelng and Smulaton of Dynamc Systems Department of Mechancal Engneerng, Unversty Of Saskatchewan, 57 Campus Drve, Saskatoon, SK S7N 5A9, Canada
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 informationA PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS
HCMC Unversty of Pedagogy Thong Nguyen Huu et al. A PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS Thong Nguyen Huu and Hao Tran Van Department of mathematcs-nformaton,
More informationTitle: Radiative transitions and spectral broadening
Lecture 6 Ttle: Radatve transtons and spectral broadenng Objectves The spectral lnes emtted by atomc vapors at moderate temperature and pressure show the wavelength spread around the central frequency.
More information2 More examples with details
Physcs 129b Lecture 3 Caltech, 01/15/19 2 More examples wth detals 2.3 The permutaton group n = 4 S 4 contans 4! = 24 elements. One s the dentty e. Sx of them are exchange of two objects (, j) ( to j and
More informationPhysics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1
P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the
More informationLab 2e Thermal System Response and Effective Heat Transfer Coefficient
58:080 Expermental Engneerng 1 OBJECTIVE Lab 2e Thermal System Response and Effectve Heat Transfer Coeffcent Warnng: though the experment has educatonal objectves (to learn about bolng heat transfer, etc.),
More informationStructure and Drive Paul A. Jensen Copyright July 20, 2003
Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.
More informationDeterministic and Monte Carlo Codes for Multiple Scattering Photon Transport
Determnstc and Monte Carlo Codes for Multple Scatterng Photon Transport Jorge E. Fernández 1 1 Laboratory of Montecuccolno DIENCA Alma Mater Studorum Unversty of Bologna Italy Isttuto Nazonale d Fsca Nucleare
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 informationECEN 5005 Crystals, Nanocrystals and Device Applications Class 19 Group Theory For Crystals
ECEN 5005 Crystals, Nanocrystals and Devce Applcatons Class 9 Group Theory For Crystals Dee Dagram Radatve Transton Probablty Wgner-Ecart Theorem Selecton Rule Dee Dagram Expermentally determned energy
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 informationOn the correction of the h-index for career length
1 On the correcton of the h-ndex for career length by L. Egghe Unverstet Hasselt (UHasselt), Campus Depenbeek, Agoralaan, B-3590 Depenbeek, Belgum 1 and Unverstet Antwerpen (UA), IBW, Stadscampus, Venusstraat
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 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 informationWorkshop: Approximating energies and wave functions Quantum aspects of physical chemistry
Workshop: Approxmatng energes and wave functons Quantum aspects of physcal chemstry http://quantum.bu.edu/pltl/6/6.pdf Last updated Thursday, November 7, 25 7:9:5-5: Copyrght 25 Dan Dll (dan@bu.edu) Department
More informationUniversity of Washington Department of Chemistry Chemistry 453 Winter Quarter 2015
Lecture 2. 1/07/15-1/09/15 Unversty of Washngton Department of Chemstry Chemstry 453 Wnter Quarter 2015 We are not talkng about truth. We are talkng about somethng that seems lke truth. The truth we want
More informationAcceleration of Gas Bubble-Free Surface Interaction Computation Using Basis Preconditioners
Acceleraton of Gas Bubble-Free Surface Interacton Computaton Usng Bass Precondtoners K. L. Tan, B. C. Khoo and J. K. Whte ABSTRACT The computaton of gas bubble-free surface nteracton entals a tme-steppng
More informationSTATISTICAL MECHANICS
STATISTICAL MECHANICS Thermal Energy Recall that KE can always be separated nto 2 terms: KE system = 1 2 M 2 total v CM KE nternal Rgd-body rotaton and elastc / sound waves Use smplfyng assumptons KE of
More informationCHAPTER 7 ENERGY BALANCES SYSTEM SYSTEM. * What is energy? * Forms of Energy. - Kinetic energy (KE) - Potential energy (PE) PE = mgz
SYSTM CHAPTR 7 NRGY BALANCS 1 7.1-7. SYSTM nergy & 1st Law of Thermodynamcs * What s energy? * Forms of nergy - Knetc energy (K) K 1 mv - Potental energy (P) P mgz - Internal energy (U) * Total nergy,
More informationSupplementary Notes for Chapter 9 Mixture Thermodynamics
Supplementary Notes for Chapter 9 Mxture Thermodynamcs Key ponts Nne major topcs of Chapter 9 are revewed below: 1. Notaton and operatonal equatons for mxtures 2. PVTN EOSs for mxtures 3. General effects
More information