As discussed in the Lecture, collapsing structures virialize. For a collapsing gas cloud,
the potential energy thereby is converted into internal energy during the collapse. Here
we will try to perform the so called Evrard collapse, which is a standard test case
and very detailed described in
this article in section 3.3, see also figure 6.
So think about:
To run the experiment, we first need to
Now we can perform the simulation and analyze it.
PERIODIC
and
NOGRAVITY
and switch on
EVALPOTENTIAL
in the Config.sh
file before compiling.
box.param
, do not
forget to also set PeriodicBoundariesOn 0
before running the simulation.
cp -r $HOME/Hydro/PPM .
TimeBetSnapshot 0.1046949
to
get the same output times.
Again, we are producing more complex initial conditions and run larger simulations, especially we will
enlarge_boxes.pro
in IDL
setup_evcoll.pro
in IDL
show_collapse
in IDL
ifort -g -traceback -check all -fpe0 -o cloudsetup cloudsetup.f90
glass_100x100x100.txt
[warning: large file: 37 Mbytes],
which is a text version of the large glass file created with enlarge_boxes.pro
above,
with coordinates in the range of 0 . . . 1)
./cloudsetup
ifort -g -traceback -check all -fpe0 -o histogram histogram.f90
./histogram >histogram.txt
gnuplot histogram.plt
gv histogram.ps
cloud.ic
as the initial conditions file
and output
as the snapshot file base in the parameter file
ifort -g -traceback -check all -fpe0 -o readsnap readsnap.f90
for file in output_???; do ./readsnap $file >$file.txt; done
gnuplot evolution.plt
gv evolution.ps