- Generate glass-like initial conditions for Smoothed Particle Hydrodynamics
- Relaxation of the particle distribution is done using al algorithm based on Weighted Voronoi Tesselations (Diehl et al. 2012)
- Additional particle reshuffling can be enabled to improve over- and undersampled maxima/minima
- A full suite of (analytical) testproblems is included
Tailor the Makefile (or CMakeLists.txt) to your machine, and run
make or cmake
Adjust the parameterfile ics.par and run the code. See the examples folder for example parameterfiles.
OMP_NUM_THREADS=4 ./WVTICs ics.par
| Problem | 1E4 | 1E5 | 1E6 |
|---|---|---|---|
| Constant Density | 1.7 | 0.8 | 0.35 |
| TopHat | 3.5 | 1.7 | 0.8 |
| SawTooth | 0.15 | 0.08 | 0.04 |
| Sinus | 1.3 | 0.6 | 0.3 |
| Mach 2 | 0.4 | 0.19 | x |
| Mach 3 | 0.4 | 0.19 | x |
| Mach 4 | 0.4 | 0.19 | x |
| Sod Shock | 5.0e-6 | 2.0e-6 | x |
| Sedov Blast | 1.7 | 0.8 | 0.35 |
| Kelvin-Helmholtz | 15.0 | 9.0 | 12.0 |
| Keplarian Ring | 14.0 | 6.0 | x |
| Cold Blob | 0.5 | 0.23 | 0.1 |
| Evrard Collapse | 1.7 | 0.8 | 0.35 |
| Zel’dovic Pancake | 1.7 | 0.8 | 0.35 |
| Box | 14.0 | 6.5 | 3.2 |
| Gresho-Vortex | 11.5 | 5.5 | x |
| Exponential Disk | 80.0 | x | x |
| Boss-Bodenheimer | 0.6 | 0.3 | 0.15 |
| Beta Profile | 1.7 | 0.8 | 0.35 |
| Fast Rotor | 12.0 | 6.0 | 3.5 |
| Strong Blast | 12.0 | 6.0 | 3.5 |
| Orszang-Tang Vort. | 12.0 | 6.0 | 3.5 |
| Lin. Alfvén Wave | 0.2 | x | x |
Contributions are welcome and will be credited
Please cite:
MIT
-
For Python
- See the Python files in the
scriptsfolder
- See the Python files in the
-
For IDL
- Do we want to add example scripts?
(c) Donnert, Arth, Steinwandel, Halbesma, Boess 2016-9