Introduction to the lattice Boltzmann method
|
|
- Mariah Snow
- 5 years ago
- Views:
Transcription
1 Introducton to LB Introducton to the lattce Boltzmann method Burkhard Dünweg Max Planck Insttute for Polymer Research Ackermannweg 10, D Manz, Germany Introducton Naver-Stokes Buldng an LB Equlbrum populatons Lnear collson Statstcal mechancs B. D. & A. J. C. Ladd, arxv: v2, Advances n Polymer Scence 221, 89 (2009
2 Lattce Boltzmann c small set of veloctes c h connects two stes lnearzed Boltzmann equaton (knetc theory of gases fully dscretzed stes r, lattce spacng a tme t, tme step h n ( r,t: real number, mass densty on ste r correspondng to velocty c Introducton to LB Introducton Naver-Stokes Buldng an LB Equlbrum populatons Lnear collson Statstcal mechancs n ( r + c h,t + h = n ( r,t = n ( r,t + {n ( r,t}
3 Lattce Boltzmann algorthm Introducton to LB ρ = n j = n c, u = j/ρ Π= n c c n eq neq (ρ, u = w ρ = ρ ( 1 + u c c 2 s + ( u c 2 2c 4 s neq c = j neq c c = ρcs 2 + ρ u u lnear relaxaton: = j L j(n j n eq j streamng u2, such that 2cs 2 Introducton Naver-Stokes Buldng an LB Equlbrum populatons Lnear collson Statstcal mechancs remarks: Π relaxaton rates determne vscostes thermal nose: Langevn stress term n Π
4 Goal of ths lecture Introducton to LB Introducton What s the connecton between lattce Boltzmann and Naver Stokes? Why 19 veloctes? Where does the equlbrum dstrbuton come from? What are the magc numbers w? Relaxaton rates vscostes? How to put n nose? Naver-Stokes Buldng an LB Equlbrum populatons Lnear collson Statstcal mechancs
5 Some phlosophcal remarks Introducton to LB n ( r + c h,t + h = n ( r,t = n ( r,t + {n ( r,t} ρ = n j = ρ u = n c = c = 0 Introducton Naver-Stokes Buldng an LB Equlbrum populatons Lnear collson Statstcal mechancs mass conservaton momentum conservaton localty rotatonal symmetry (lattce! Galle nvarance (fnte number of veloctes
6 only u c only u c s Ma = u/c s 1 low Mach number compressblty does not matter equaton of state does not matter choose deal gas! m p partcle mass: p = ρ m p k B T Introducton to LB Introducton Naver-Stokes Buldng an LB Equlbrum populatons Lnear collson Statstcal mechancs c 2 s = p ρ = 1 m p k B T p = ρc 2 s k B T = m p c 2 s
7 Necessty of thermal fluctuatons Ideal gas, temp. T, partcle mass m p, sound speed c s : k B T = m p c 2 s c s a/h (a lattce spacng, h tme step c s as small as possble Example (water: mass densty ρ = 10 3 kg/m 3 sound speed realstc: m/s sound speed artfcal: c s = 10 2 m/s temperature T 300K, k B T = J partcle mass: m p = kg macroscopc scale molecular scale lattce spacng a = 1mm a = 1nm tme step h = 10 5 s h = s mass of a ste m a = 10 6 kg m a = kg Boltzmann Bo = (m p /m a 1/2 Bo = (m p /m a 1/2 number = = 0.6 Introducton to LB Introducton Naver-Stokes Buldng an LB Equlbrum populatons Lnear collson Statstcal mechancs
8 Where we want to get Introducton to LB Introducton t ρ + α j α = 0 Naver-Stokes η αβγδ = t j α + β ( ρc 2 s δ αβ + ρu α u β = β σ αβ σ αβ = η αβγδ γ u δ (ζ 23 η δ αβ δ γδ + η (δ αγ δ βδ + δ αδ δ βγ Buldng an LB Equlbrum populatons Lnear collson Statstcal mechancs η shear vscosty ζ bulk vscosty
9 expanson Mult tme scale analyss: ε 1 (e. g. ε = 10 3 : r 1 = ε r nterpretaton: coarse graned ruler : 1µm nstead of 978nm t 1 = εt nterpretaton: coarse graned clock : 1ns nstead of 837ps t 2 = ε 2 t nterpretaton: yet more coarse graned clock : 1µs nstead of 976ns572ps t 1 to capture wave lke phenomena t 2 to capture dffusve phenomena Introducton to LB Introducton Naver-Stokes Buldng an LB Equlbrum populatons Lnear collson Statstcal mechancs locaton n space and tme: read off r 1, t 1, t 2
10 Introducton to LB r 1 fxed r vares wth ε LBE: r = ε 1 r 1 n = n (0 + εn (1 + O(ε 2 = (0 + ε (1 + ε 2 (2 + O(ε 3 (k = (k c = 0 Introducton Naver-Stokes Buldng an LB Equlbrum populatons Lnear collson Statstcal mechancs n ( r 1 + ε c h,t 1 + εh,t 2 + ε 2 h n ( r 1,t 1,t 2 = ( r 1,t 1,t 2 Taylor expanson wrt ε ε 0 corresponds to macroscopc lmt
11 ε orders Introducton to LB ε 0 : ε 1 : ε 2 : (0 = 0 ( t1 + c r1 n (0 = h 1 (1 t2 n (0 + h 2 ( t 1 + c r1 2 n (0 + ( t1 + c r1 n (1 = h 1 (2 Introducton Naver-Stokes Buldng an LB Equlbrum populatons Lnear collson Statstcal mechancs or t2 n (0 + 1 ( 2 ( t 1 + c r1 n (1 + n (1 = h 1 (2
12 Zeroth order Introducton to LB 0 = (0 = ({n (0 } {n (0 } collsonal nvarant, {n (0 } = n eq no spurous conservaton laws n (0 = n (0 (ρ, j ρ (0 = ρ, j (0 = j n eq = ρ n eq c = j hgher orders: take moments..., c... Introducton Naver-Stokes Buldng an LB Equlbrum populatons Lnear collson Statstcal mechancs
13 Zeroth moment: Mass conservaton Introducton to LB Now, t1 ρ + 1α j α = 0 t2 ρ = 0 α = ε 1α t = ε t1 + ε 2 t2 Introducton Naver-Stokes Buldng an LB Equlbrum populatons Lnear collson Statstcal mechancs Hence, contnuty equaton OK!!!
14 Frst moment: Momentum conservaton Introducton to LB wth t2 j α β t1 j α + 1β Π (0 αβ = 0 Π αβ = ( Π (1 αβ + Π(1 αβ = 0 n c α c β Introducton Naver-Stokes Buldng an LB Equlbrum populatons together: t j α + β Π eq αβ + 1 ( 2 β Π neq αβ + Π neq αβ = 0 Lnear collson Statstcal mechancs wth n neq = εn (1 comparson wth Naver Stokes: Π eq αβ = ρc2 s δ 1 αβ + ρu α u β 2 ( Π neq αβ + Π neq αβ = σ αβ
15 Second moment: A useful relaton wth or ( t1 Π (0 αβ + 1γΦ (0 αβγ = h 1 Π (1 αβ Π(1 αβ Π neq αβ later (from explct form of n eq : Φ αβγ = n c α c β c γ ( Π neq αβ = h t Π eq αβ + γφ eq αβγ Φ eq αβγ = ρc2 s (u αδ βγ + u β δ αγ + u γ δ αβ use contnuty and Euler ( t1!!! for t Π eq αβ = t(ρc 2 s δ αβ + ρu α u β result (neglectng terms O(u 3 : Π neq αβ Π neq αβ = hρc2 s ( α u β + β u α Introducton to LB Introducton Naver-Stokes Buldng an LB Equlbrum populatons Lnear collson Statstcal mechancs
16 Lnear stress relaxaton Introducton to LB later: together wth Π neq αβ Π αβ = Π αβ δ αβπ γγ Π neq αβ Π neq αα = γ s Πneq αβ = γ bπ neq αα Π neq αβ = hρc2 s ( αu β + β u α two equatons wth two unknowns Π neq, Π neq result: 1 ( Π neq 2 αβ + Π neq αβ = σ αβ wth η = hρc2 s γ s ζ = hρc2 s 1 γ s γ b 1 γ b Introducton Naver-Stokes Buldng an LB Equlbrum populatons Lnear collson Statstcal mechancs
17 So, what do we need? 1. n eq such that: n eq (and hopefully = ρ n eq c = j n eq c α c β = ρc 2 s δ αβ + ρu α u β n eq c α c β c γ = ρc 2 s (u α δ βγ + u β δ αγ + u γ δ αβ 2. such that n eq nvarant, and = 0 Π neq αβ = γ neq s Π αβ c = 0 Π neq αα = γ bπ neq αα Introducton to LB Introducton Naver-Stokes Buldng an LB Equlbrum populatons Lnear collson Statstcal mechancs
18 D3Q19 I: Equlbrum populatons ansatz: n eq (ρ, u = w ρ ( 1 + A u c + B( u c 2 + Cu 2 cubc symmetry: w = 1 w c α = 0 w c α c β = σ 2 δ αβ w c α c β c γ = 0 w c α c β c γ c δ = κ 4 δ αβγδ +σ 4 (δ αβ δ γδ + δ αγ δ βδ + δ αδ δ βγ Introducton to LB Introducton Naver-Stokes Buldng an LB Equlbrum populatons Lnear collson Statstcal mechancs
19 ρ, j, Π: n eq soluton: (ρ, u = w ρ ( 1 + u c c 2 s c 2 s = σ 2 + ( u c 2 2c 4 s u2 2cs 2 w = 1 σ 4 = σ2 2 κ 4 = 0 w = 1/3 (zero velocty w = 1/18 (6 nearest neghbors w = 1/36 (12 next nearest neghbors Introducton to LB Introducton Naver-Stokes Buldng an LB Equlbrum populatons Lnear collson Statstcal mechancs c 2 s = (1/3(a/h2 3rd moment OK
20 D3Q19 II: Lnear collson ansatz (n eq nvarant: n neq = j Γ j n neq j modes dagonal representaton: m k = e k n m 0 ρ (m 1,m 2,m 3 j m 4 span(ρ,π αα span(m 5,...,m 9 = span( Π αβ knetc modes m 10,...,m 19 m neq k conserved modes: egenvalue rrelevant: γ 0 =... = γ 3 = 1 bulk stress: γ 4 = γ b shear stresses: γ 5 =...γ 9 = γ s knetc modes: γ k = 0 (easest choce = γ k m neq k Introducton to LB Introducton Naver-Stokes Buldng an LB Equlbrum populatons Lnear collson Statstcal mechancs
21 Statstcal mechancs: Theory of fluctuatons B. D., U. Schller, A. J. C. Ladd, PRE 76, (2007 µ = m p /a 3, ν = n /µ, ν = w ρ/µ P ({ν } ( ν ν e ν δ µ ( ν ρ δ µ ν! Strlng: ν c j Introducton to LB Introducton Naver-Stokes Buldng an LB Equlbrum populatons S ({ν } = (ν lnν ν ν ln ν + ν ( ( Lnear collson Statstcal mechancs P exp[s ({ν }] δ µ ν ρ δ µ ν c j P = max S = max wth constrants soluton (neglectng terms O(u 3 : n eq as before
22 Fluctuaton statstcs Gaussan approxmaton around the mean: n eq P exp ( ρw ( P exp ( n neq 2 2µn eq ( n neq 2µρw 2 recall: unnormalzed modes: δ δ ( ( n neq n neq δ δ ( ( c n neq c n neq Introducton to LB Introducton Naver-Stokes Buldng an LB Equlbrum populatons Lnear collson Statstcal mechancs m k = e k n normalzed modes ˆm k m k : ˆm k = ê k n w µρ wth orthogonal matrx ê k
23 Fluctuaton statstcs n mode space, Monte Carlo Introducton to LB P r k Gaussan wth neq} ({ˆm k exp 1 2 k 4 ˆm neq k = γ k ˆm neq k + ϕ k r k ˆm neq 2 k Introducton Naver-Stokes Buldng an LB Equlbrum populatons Lnear collson r k = 0 r 2 k = 1 Statstcal mechancs detaled balance holds for ϕ k = ( 1 γk 2 1/2
24 Fluctuatng hydrodynamcs Introducton to LB appled to the stochastc term equatons of fluctuatng hydrodynamcs (Landau / Lfshtz 1957 t ρ + α j α = 0 t j α + β ( ρc 2 s δ αβ + ρu α u β = β η αβγδ γ u δ + β Q αβ Introducton Naver-Stokes Buldng an LB Equlbrum populatons Lnear collson Statstcal mechancs Q αβ = 0 Qαβ ( r,t Q γδ ( r,t = 2k B Tη αβγδ δ ( r r δ ( t t
Progress in the Understanding of the Fluctuating Lattice Boltzmann Equation
Progress n the Understandng of the Fluctuatng Lattce Boltzmann Equaton Burkhard Dünweg a, Ulf D. Schller a, Anthony J. C. Ladd b a Max Planck Insttute for Polymer Research, Ackermannweg 10, D-55128 Manz,
More informationLATTICE BOLTZMANN SIMULATIONS OF SOFT
LATTICE BOLTZMANN SIMULATIONS OF SOFT MATTER SYSTEMS Burkhard Dünweg 1 and Anthony J. C. Ladd 2 1 Max Planck Insttute for Polymer Research Ackermannweg 10, 55128 Manz, Germany duenweg@mpp-manz.mpg.de 2
More informationCoupling Molecular Dynamics and Lattice Boltzmann to Simulate Brownian Motion with Hydrodynamic Interactions
Couplng Molecular Dynamcs and Lattce Boltzmann to Smulate Brownan Moton wth Hydrodynamc Interactons Burkhard Dünweg Max Planck Insttute for Polymer Research Ackermannweg 10, 55128 Manz, Germany and Department
More informationNumerical Simulations of Particulate Suspensions via a Discretized Boltzmann Equation Part I. Theoretical Foundation
Under consderaton for publcaton n J. Flud Mech. 1 Numercal Smulatons of Partculate Suspensons va a Dscretzed Boltzmann Equaton Part I. Theoretcal Foundaton By ANTHONY J. C. LADD Lawrence Lvermore Natonal
More informationand Statistical Mechanics Material Properties
Statstcal Mechancs and Materal Propertes By Kuno TAKAHASHI Tokyo Insttute of Technology, Tokyo 15-855, JAPA Phone/Fax +81-3-5734-3915 takahak@de.ttech.ac.jp http://www.de.ttech.ac.jp/~kt-lab/ Only for
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 informationLesson 5: Kinematics and Dynamics of Particles
Lesson 5: Knematcs and Dynamcs of Partcles hs set of notes descrbes the basc methodology for formulatng the knematc and knetc equatons for multbody dynamcs. In order to concentrate on the methodology and
More informationTHEORY OF THE LATTICE BOLTZMANN METHOD FOR MULTI-PHASE AND MULTICOMPONENT FLUIDS
THEORY OF THE LATTICE BOLTZMANN METHOD FOR MULTI-PHASE AND MULTICOMPONENT FLUIDS A Thess Submtted to the Graduate Faculty of the North Dakota State Unversty of Agrculture and Appled Scence By Qun L In
More informationBasic concept of reactive flows. Basic concept of reactive flows Combustion Mixing and reaction in high viscous fluid Application of Chaos
Introducton to Toshhsa Ueda School of Scence for Open and Envronmental Systems Keo Unversty, Japan Combuston Mxng and reacton n hgh vscous flud Applcaton of Chaos Keo Unversty 1 Keo Unversty 2 What s reactve
More 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 informationClassical Mechanics ( Particles and Biparticles )
Classcal Mechancs ( Partcles and Bpartcles ) Alejandro A. Torassa Creatve Commons Attrbuton 3.0 Lcense (0) Buenos Ares, Argentna atorassa@gmal.com Abstract Ths paper consders the exstence of bpartcles
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 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 informationPhysics 141. Lecture 14. Frank L. H. Wolfs Department of Physics and Astronomy, University of Rochester, Lecture 14, Page 1
Physcs 141. Lecture 14. Frank L. H. Wolfs Department of Physcs and Astronomy, Unversty of Rochester, Lecture 14, Page 1 Physcs 141. Lecture 14. Course Informaton: Lab report # 3. Exam # 2. Mult-Partcle
More informationNumerical Simulations of Particulate Suspensions. via a Discretized Boltzmann Equation. Part I. Theoretical Foundation. (June 25, 1993) Abstract
UCRL-JC-113349 Numercal Smulatons of Partculate Suspensons va a Dscretzed Boltzmann Equaton Part I. Theoretcal Foundaton Anthony J. C. Ladd Lawrence Lvermore Natonal Laboratory, Lvermore, Calforna 94550.
More informationLecture 5: Ideal monatomic gas
Lecture 5: Ideal monatomc gas Statstcal mechancs of the perfect gas Ams: Key new concepts and methods: Countng states Waves n a box. Demonstraton that β / kt Heat, work and Entropy n statstcal mechancs
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 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 informationComputers and Mathematics with Applications
Computers and Mathematcs wth Applcatons 58 (009) 86 866 Contents lsts avalable at ScenceDrect Computers and Mathematcs wth Applcatons journal homepage: www.elsever.com/locate/camwa An Onsager-lke relaton
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 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 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 informationLecture 6/7 (February 10/12, 2014) DIRAC EQUATION. The non-relativistic Schrödinger equation was obtained by noting that the Hamiltonian 2
P470 Lecture 6/7 (February 10/1, 014) DIRAC EQUATION The non-relatvstc Schrödnger equaton was obtaned by notng that the Hamltonan H = P (1) m can be transformed nto an operator form wth the substtutons
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 informationSimulation 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 information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 6
REVIEW of Lecture 5 2.29 Numercal Flud Mechancs Fall 2011 Lecture 6 Contnuum Hypothess and conservaton laws Macroscopc Propertes Materal covered n class: Dfferental forms of conservaton laws Materal Dervatve
More informationNumerical 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 informationComputational Astrophysics
Computatonal Astrophyscs Solvng for Gravty Alexander Knebe, Unversdad Autonoma de Madrd Computatonal Astrophyscs Solvng for Gravty the equatons full set of equatons collsonless matter (e.g. dark matter
More informationNotes on Analytical Dynamics
Notes on Analytcal Dynamcs Jan Peters & Mchael Mstry October 7, 004 Newtonan Mechancs Basc Asssumptons and Newtons Laws Lonely pontmasses wth postve mass Newtons st: Constant velocty v n an nertal frame
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 informationπ e ax2 dx = x 2 e ax2 dx or x 3 e ax2 dx = 1 x 4 e ax2 dx = 3 π 8a 5/2 (a) We are considering the Maxwell velocity distribution function: 2πτ/m
Homework Solutons Problem In solvng ths problem, we wll need to calculate some moments of the Gaussan dstrbuton. The brute-force method s to ntegrate by parts but there s a nce trck. The followng ntegrals
More informationarxiv: v1 [physics.comp-ph] 20 Nov 2018
arxv:8.08530v [physcs.comp-ph] 20 Nov 208 Internatonal Journal of Modern Physcs C c World Scentfc Publshng Company Valdty of the Molecular-Dynamcs-Lattce-Gas Global Equlbrum Dstrbuton Functon M. Reza Parsa
More informationStatistical mechanics canonical ensemble
canoncal ensemble system n thermal equlbrum wth bath of system probablty of mcro state T = k/ bath system p = 1 Z e E Z = e E average energy of system he = p E = what about? P E E e P e E = @ ln Z @ replcate
More informationPhysics 212: Statistical mechanics II Lecture I
Physcs 212: Statstcal mechancs II Lecture I A theory s the more mpressve the greater the smplcty of ts premses, the more dfferent knds of thngs t relates, and the more extended ts area of applcablty. Therefore
More informationAxial Turbine Analysis
Axal Turbne Analyss From Euler turbomachnery (conservaton) equatons need to Nole understand change n tangental velocty to relate to forces on blades and power m W m rc e rc uc uc e Analye flow n a plane
More informationIntroduction to Statistical Methods
Introducton to Statstcal Methods Physcs 4362, Lecture #3 hermodynamcs Classcal Statstcal Knetc heory Classcal hermodynamcs Macroscopc approach General propertes of the system Macroscopc varables 1 hermodynamc
More informationPHY688, Statistical Mechanics
Department of Physcs & Astronomy 449 ESS Bldg. Stony Brook Unversty January 31, 2017 Nuclear Astrophyscs James.Lattmer@Stonybrook.edu Thermodynamcs Internal Energy Densty and Frst Law: ε = E V = Ts P +
More informationFluid-dynamic equations for reacting gas mixtures
Flud-dynamc equatons for reactng gas mxtures Marza Bs, Mara Gropp, Gampero Spga Dpartmento d Matematca, Unverstà d Parma Va D Azeglo 85, 43100 Parma, Italy e-mal: gampero.spga@unpr.t Abstract Startng from
More informationPHYS 215C: Quantum Mechanics (Spring 2017) Problem Set 3 Solutions
PHYS 5C: Quantum Mechancs Sprng 07 Problem Set 3 Solutons Prof. Matthew Fsher Solutons prepared by: Chatanya Murthy and James Sully June 4, 07 Please let me know f you encounter any typos n the solutons.
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 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 informationU.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017
U.C. Berkeley CS94: Beyond Worst-Case Analyss Handout 4s Luca Trevsan September 5, 07 Summary of Lecture 4 In whch we ntroduce semdefnte programmng and apply t to Max Cut. Semdefnte Programmng Recall that
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 informationJohn E. Thomas. Quark-gluon plasma T = K BIG BANG Computer simulation of RHIC collision. Ultracold atomic gas T = 10-7 K
Quantum hydrodynamcs n a strongly nteractng Ferm gas John E. Thomas Quark-gluon plasma T 10 1 K BIG BANG Computer smulaton of RHIC collson Ultracold atomc gas T 10-7 K JETLa Group Students: Yngy Zhang
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 informationLecture Notes on Linear Regression
Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume
More informationMechanics Physics 151
Mechancs Physcs 5 Lecture 3 Contnuous Systems an Fels (Chapter 3) Where Are We Now? We ve fnshe all the essentals Fnal wll cover Lectures through Last two lectures: Classcal Fel Theory Start wth wave equatons
More informationSome Aspects of the Boltzmann Equation for a Granular Gas
Some Aspects of the Boltzmann Equaton for a Granular Gas James W. Dufty Department of Physcs, Unversty of Florda, Ganesvlle, FL 326, USA Abstract. The Boltzmann equaton for a gas of smooth, nelastc hard
More informationThe Solution of the Two-Dimensional Gross-Pitaevskii Equation Using Lattice-Boltzmann and He s Semi-Inverse Method
Internatonal Journal of Mathematcal Analyss Vol., 7, no., 69-77 HIKARI Ltd, www.m-hkar.com https://do.org/.988/jma.7.634 The Soluton of the Two-Dmensonal Gross-Ptaevsk Equaton Usng Lattce-Boltzmann and
More informationIonization fronts in HII regions
Ionzaton fronts n HII regons Intal expanson of HII onzaton front s supersonc, creatng a shock front. Statonary frame: front advances nto neutral materal In frame where shock front s statonary, neutral
More informationThermodynamics General
Thermodynamcs General Lecture 1 Lecture 1 s devoted to establshng buldng blocks for dscussng thermodynamcs. In addton, the equaton of state wll be establshed. I. Buldng blocks for thermodynamcs A. Dmensons,
More informationThe Tangential Force Distribution on Inner Cylinder of Power Law Fluid Flowing in Eccentric Annuli with the Inner Cylinder Reciprocating Axially
Open Journal of Flud Dynamcs, 2015, 5, 183-187 Publshed Onlne June 2015 n ScRes. http://www.scrp.org/journal/ojfd http://dx.do.org/10.4236/ojfd.2015.52020 The Tangental Force Dstrbuton on Inner Cylnder
More information763622S ADVANCED QUANTUM MECHANICS Solution Set 1 Spring c n a n. c n 2 = 1.
7636S ADVANCED QUANTUM MECHANICS Soluton Set 1 Sprng 013 1 Warm-up Show that the egenvalues of a Hermtan operator  are real and that the egenkets correspondng to dfferent egenvalues are orthogonal (b)
More informationDUE: WEDS FEB 21ST 2018
HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant
More informationIrreversibility of Processes in Closed System
Unversty of Segen Insttute of Flud- & hermodynamcs 5 2/1 Irreversblty of Processes n Closed System m G 2 m c 2 2, p, V m g h h 1 mc 1 1 p, p, V G J.P. Joule Strrng experment v J.B. Fourer Heat transfer
More informationPoisson brackets and canonical transformations
rof O B Wrght Mechancs Notes osson brackets and canoncal transformatons osson Brackets Consder an arbtrary functon f f ( qp t) df f f f q p q p t But q p p where ( qp ) pq q df f f f p q q p t In order
More informationChapter 13: Multiple Regression
Chapter 13: Multple Regresson 13.1 Developng the multple-regresson Model The general model can be descrbed as: It smplfes for two ndependent varables: The sample ft parameter b 0, b 1, and b are used to
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 informationEffects of Polymer Concentration and Molecular Weight on the Dynamics of Visco-Elasto- Capillary Breakup
Effects of Polymer Concentraton and Molecular Weght on the Dynamcs of Vsco-Elasto- Capllary Breakup Mattheu Veran Advsor: Prof. Gareth McKnley Mechancal Engneerng Department January 3, Capllary Breakup
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 informationStudy Guide For Exam Two
Study Gude For Exam Two Physcs 2210 Albretsen Updated: 08/02/2018 All Other Prevous Study Gudes Modules 01-06 Module 07 Work Work done by a constant force F over a dstance s : Work done by varyng force
More informationTHERMAL DISTRIBUTION IN THE HCL SPECTRUM OBJECTIVE
ame: THERMAL DISTRIBUTIO I THE HCL SPECTRUM OBJECTIVE To nvestgate a system s thermal dstrbuton n dscrete states; specfcally, determne HCl gas temperature from the relatve occupatons of ts rotatonal states.
More informationA micro-macro decomposition based asymptotic-preserving scheme for the multispecies Boltzmann equation
A mcro-macro decomposton based asymptotc-preservng scheme for the multspeces Boltzmann equaton Sh Jn Yngzhe Sh Abstract In ths paper we extend the mcro-macro decomposton based asymptotc-preservng scheme
More informationQuantum Statistical Mechanics
Chapter 8 Quantum Statstcal Mechancs 8.1 Mcrostates For a collecton of N partcles the classcal mcrostate of the system s unquely specfed by a vector (q, p) (q 1...q N, p 1...p 3N )(q 1...q 3N,p 1...p 3N
More informationRelativistic Lattice Boltzmann Model with Improved Dissipation
Relatvstc Lattce Boltzmann Model wth Improved Dsspaton M. Mendoza, 1, I. Karln, 2, S. Succ, 3, 1, 4, and H. J. Herrmann 1 TH Zürch, Computatonal Physcs for ngneerng Materals, Insttute for Buldng Materals,
More informationAN INTRODUCTION TO ENTROPIC LATTICE BOLTZMANN SCHEME
AN INTRODUCTION TO ENTROPIC LATTICE BOLTZMANN SCHEME Francesca Tos 1 and Sauro Succ 2 1 Dp. d Matematca U.Dn Unverstá degl Stud d Frenze, Vale Morgagn, 67-a, 50134 Frenze, Italy tos@math.unf.t 2 Isttuto
More informationPolymer Chains. Ju Li. GEM4 Summer School 2006 Cell and Molecular Mechanics in BioMedicine August 7 18, 2006, MIT, Cambridge, MA, USA
Polymer Chans Ju L GEM4 Summer School 006 Cell and Molecular Mechancs n BoMedcne August 7 18, 006, MIT, Cambrdge, MA, USA Outlne Freely Jonted Chan Worm-Lke Chan Persstence Length Freely Jonted Chan (FJC)
More informationII.D Many Random Variables
II.D Many Random Varables Wth more than one random varable, the set of outcomes s an -dmensonal space, S x = { < x 1, x 2,, x < }. For example, descrbng the locaton and velocty of a gas partcle requres
More informationTurbulence. Lecture 21. Non-linear Dynamics. 30 s & 40 s Taylor s work on homogeneous turbulence Kolmogorov.
Turbulence Lecture 1 Non-lnear Dynamcs Strong non-lnearty s a key feature of turbulence. 1. Unstable, chaotc behavor.. Strongly vortcal (vortex stretchng) 3 s & 4 s Taylor s work on homogeneous turbulence
More informationQuantum Field Theory III
Quantum Feld Theory III Prof. Erck Wenberg February 16, 011 1 Lecture 9 Last tme we showed that f we just look at weak nteractons and currents, strong nteracton has very good SU() SU() chral symmetry,
More informationPhysics 240: Worksheet 30 Name:
(1) One mole of an deal monatomc gas doubles ts temperature and doubles ts volume. What s the change n entropy of the gas? () 1 kg of ce at 0 0 C melts to become water at 0 0 C. What s the change n entropy
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 informationParametric fractional imputation for missing data analysis. Jae Kwang Kim Survey Working Group Seminar March 29, 2010
Parametrc fractonal mputaton for mssng data analyss Jae Kwang Km Survey Workng Group Semnar March 29, 2010 1 Outlne Introducton Proposed method Fractonal mputaton Approxmaton Varance estmaton Multple mputaton
More informationThermodynamics and statistical mechanics in materials modelling II
Course MP3 Lecture 8/11/006 (JAE) Course MP3 Lecture 8/11/006 Thermodynamcs and statstcal mechancs n materals modellng II A bref résumé of the physcal concepts used n materals modellng Dr James Ellott.1
More informationHidden Markov Models & The Multivariate Gaussian (10/26/04)
CS281A/Stat241A: Statstcal Learnng Theory Hdden Markov Models & The Multvarate Gaussan (10/26/04) Lecturer: Mchael I. Jordan Scrbes: Jonathan W. Hu 1 Hdden Markov Models As a bref revew, hdden Markov models
More informationEfficient, General Point Cloud Registration with Kernel Feature Maps
Effcent, General Pont Cloud Regstraton wth Kernel Feature Maps Hanchen Xong, Sandor Szedmak, Justus Pater Insttute of Computer Scence Unversty of Innsbruck 30 May 2013 Hanchen Xong (Un.Innsbruck) 3D Regstraton
More informationLeast squares cubic splines without B-splines S.K. Lucas
Least squares cubc splnes wthout B-splnes S.K. Lucas School of Mathematcs and Statstcs, Unversty of South Australa, Mawson Lakes SA 595 e-mal: stephen.lucas@unsa.edu.au Submtted to the Gazette of the Australan
More informationQuantum Mechanics I - Session 4
Quantum Mechancs I - Sesson 4 Aprl 3, 05 Contents Operators Change of Bass 4 3 Egenvectors and Egenvalues 5 3. Denton....................................... 5 3. Rotaton n D....................................
More informationChapter 11 Angular Momentum
Chapter 11 Angular Momentum Analyss Model: Nonsolated System (Angular Momentum) Angular Momentum of a Rotatng Rgd Object Analyss Model: Isolated System (Angular Momentum) Angular Momentum of a Partcle
More informationIntroduction to Turbulence Modeling
Introducton to Turbulence Modelng Professor Ismal B. Celk West Vrgna nversty Ismal.Celk@mal.wvu.edu CFD Lab. - West Vrgna nversty I-1 Introducton to Turbulence CFD Lab. - West Vrgna nversty I-2 Introducton
More informationLagrangian Field Theory
Lagrangan Feld Theory Adam Lott PHY 391 Aprl 6, 017 1 Introducton Ths paper s a summary of Chapter of Mandl and Shaw s Quantum Feld Theory [1]. The frst thng to do s to fx the notaton. For the most part,
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 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 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 informationConvexity preserving interpolation by splines of arbitrary degree
Computer Scence Journal of Moldova, vol.18, no.1(52), 2010 Convexty preservng nterpolaton by splnes of arbtrary degree Igor Verlan Abstract In the present paper an algorthm of C 2 nterpolaton of dscrete
More informationGouy-Chapman model (1910) The double layer is not as compact as in Helmholtz rigid layer.
CHE465/865, 6-3, Lecture 1, 7 nd Sep., 6 Gouy-Chapman model (191) The double layer s not as compact as n Helmholtz rgd layer. Consder thermal motons of ons: Tendency to ncrease the entropy and make the
More informationOverview of the lattice Boltzmann method for nano- and microscale fluid dynamics in materials science and engineering
INSTITUTE OF PHYSICS PUBLISHING MODELLING AND SIMULATION IN MATERIALS SCIENCE AND ENGINEERING Modellng Smul. Mater. Sc. Eng. 12 (2004) R13 R46 PII: S0965-0393(04)77661-5 TOPICAL REVIEW Overvew of the lattce
More informationRigid body simulation
Rgd bod smulaton Rgd bod smulaton Once we consder an object wth spacal etent, partcle sstem smulaton s no longer suffcent Problems Problems Unconstraned sstem rotatonal moton torques and angular momentum
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 information8. Superfluid to Mott-insulator transition
8. Superflud to Mott-nsulator transton Overvew Optcal lattce potentals Soluton of the Schrödnger equaton for perodc potentals Band structure Bloch oscllaton of bosonc and fermonc atoms n optcal lattces
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 informationThe Governing Equations
The Governng Equatons L. Goodman General Physcal Oceanography MAR 555 School for Marne Scences and Technology Umass-Dartmouth Dynamcs of Oceanography The Governng Equatons- (IPO-7) Mass Conservaton and
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 informationThe Feynman path integral
The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space
More information10/9/2003 PHY Lecture 11 1
Announcements 1. Physc Colloquum today --The Physcs and Analyss of Non-nvasve Optcal Imagng. Today s lecture Bref revew of momentum & collsons Example HW problems Introducton to rotatons Defnton of angular
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 4 Feb 1997
Statc and Dynamc Propertes of Dsspatve Partcle Dynamcs arxv:cond-mat/970036v1 [cond-mat.stat-mech] 4 Feb 1997 C. Marsh Theoretcal Physcs, Oxford Unversty 1 Keble Road, Oxford OX1 3NP UK and G. Backx and
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 informationA conservative and entropy scheme for a simplified model of granular media
A conservatve and entropy scheme for a smplfed model of granular meda C. Buet 1, S. Corder 2, V. Dos Santos 3 1: CEA,91680 Bruyères le Chatel, France. emal: chrstophe.buet@cea.fr 2: Laboratore MAPMO,UMR
More informationLecture 7: Boltzmann distribution & Thermodynamics of mixing
Prof. Tbbtt Lecture 7 etworks & Gels Lecture 7: Boltzmann dstrbuton & Thermodynamcs of mxng 1 Suggested readng Prof. Mark W. Tbbtt ETH Zürch 13 März 018 Molecular Drvng Forces Dll and Bromberg: Chapters
More informationThe lattice Boltzmann method: Fundamentals and acoustics
Erlend Magnus Vggen The lattce Boltzmann method: Fundamentals and acoustcs Thess for the degree of Phlosophae Doctor Trondhem, February 2014 Norwegan Unversty of Scence and Technology Faculty of Informaton
More information