RESULTS (updated 2nd June 2008)

using PHANTOM SPH code
results of 256^3 particle hydro, 512^3 particle hydro, 256^3 particle MHD and 512^3 particle MHD are contained in the following file:





























SOME OLD STUFF (ignore, updated version in the pdf file)

velocity slices mid-way through the box


vx slices taken at z=0.5 (ie. half way through the box) in my 256^3 particle run (plots are xy): times shown are from t=0.02 to t=0.20 using an adaptive colour table
vslicexy.png

density slices mid-way through the box

below are density slices taken at z=0.5 (ie. half way through the box) in my 256^3 particle run (plots are xy): times shown are from t=0.02 to t=0.20 using an adaptive colour table.
Note that the crappyness in the initial conditions is NOT OF MY DOING! This can be seen by the fact that it improves as the run proceeds.
slicexy.png

here is the full horror of the initial conditions for density in the SPH data: ie. THE INITIAL CONDITIONS NEED IMPROVEMENT!
slicexy_00001.png

column density integrated through z

here are column density snapshot images from the 256^3 particle run (plots are xy) using a fixed colour table:
t=0.0
pricexy_00001.png
t=0.02
pricexy_00002.png
t=0.04
pricexy_00003.png
t=0.06
pricexy_00004.png
t=0.08
pricexy_00005.png
t=0.1
pricexy_00006.png
t=0.12
pricexy_00007.png
t=0.14
pricexy_00008.png
t=0.16
pricexy_00009.png
t=0.18
pricexy_00010.png
t=0.2
pricexy_00011.png

power spectra

power spectra in the 256^3 run, after interpolating from the particles to a 256^3 grid
(for the full comparison, I hope to calculate a power spectrum direct from the particles
so as not to lose information at high k due to the interpolation). Importantly, here the
interpolation of velocity is normalised, whilst the density interpolation is "raw".
Variations on this combination have also been done, but this seems to be the best approach
(well, the second best approach compared to getting the power spectrum direct from particles).