N-Body - attempt at spiral galaxy simulation

N_darkmatter = 25000
N_bulge = 15000
N_Disk = 60000

I say "attempt" because I am not sure how well I have constructed the initial conditions

Bulge - plummer sphere
Disk - attempted Kuzmin disk (q=1.8)
Dark Matter Halo = attempted isothermal sphere

Only the disk masses are plotted in this, all particles were advanced using Barnes Hut Treecode. Timesteps of 0.01 were used with a softening factor of 0.05. Unitless mass proportions for the bulge,disk, and DM halo are 0.2,1.0, and 15.0. The diameter of the disk is 10.0. The simulation advances from t = 0 to t = 170 in this video. The calculations are still running I paused it to output the plots as .png files with a gnuplot script and then combined them with mencoder.

Category:

Science & Technology

Tags:

• nbody
• n-body
• galactic dynamics
• astronomy
• astrophysics
• computer simulation
• computer
• simulation
• barnes
• hut
• treecode
• barnes hut treecode
• disk
• bulge
• halo
• dark matter
• gnuplot
• mencoder
• galaxy
• spiral galaxy
• fortran
• computer programming
• numerical simulation
• disk galaxy