Monday, 07/24: Writing a pm-code - AIP

More documents

Recommendations

Info

PM: main <strong>code</strong> blocks 21 • Triangular Shaped Cloud (TSC): S(x) = 1 ∆x { 1 − |x|/∆x, |x| < ∆x 0, otherwise The fraction of particle’s mass assigned to a cell ijk is the shape function averaged over this cell: W (x p − x ijk ) = x ijk +∆x/2 ∫ dx ′ S(x p − x ′ ); x ijk −∆x/2 W (r p − r ijk ) = W (x p − x ijk )W (y p − y ijk )W (z p − z ijk ); The density in a cell ijk is then ρ ijk = N p ∑ p=1 m p W (r p − r ijk )
PM: main <strong>code</strong> blocks 22 In practice, of course, we loop over particles and assign the density to neighboring cells as opposed to summing over all particles for each cell as the straightforward reading of the above equation would suggest.
Page 1 and 2: Writing a PM code Andrey Kravtsov a
Page 3 and 4: equations I 2 Cosmological simulati
Page 5 and 6: N-body approach 4 Cosmological simu
Page 7 and 8: particle initial conditions 6 Cosmo
Page 9 and 10: equations of gasdynamics 8 Cosmolog
Page 11 and 12: equations summary 10 • Baryonic g
Page 13 and 14: dissipationless simulations 12 Part
Page 15 and 16: Model equations 14 Model equations
Page 17 and 18: Model equations 16 The quantities w
Page 19 and 20: Model equations 18 In dimensionless
Page 21: PM: main code blocks 20 Density Ass
Page 25 and 26: PM: main code blocks 24 Implementin
Page 27 and 28: PM: main code blocks 26 Solving the
Page 29 and 30: PM: main code blocks 28 ⇒ Given a
Page 31 and 32: PM: main code blocks 30 After n tim
Page 33 and 34: Zeldovich approximation (ZA) 32 Zel
Page 35 and 36: A test problem: 1D collapse of a pl
Page 41 and 42: Setting up cosmological ICs 40 The
Page 43 and 44: Setting up cosmological ICs 42 Appl
Page 45 and 46: project summary 44 Project summary
Page 47 and 48: References and links 46 PM in cosmo
Page 49 and 50: References and links 48 • Michael

Monday, 07/24: Writing a pm-code - AIP

You also want an ePaper? Increase the reach of your titles

Delete template?

Save as template?