Standard Hydrodynamics Suite

This page lists a set of standardised hydrodynamics tests with analytic solutions.

The Sedov-Taylor blast wave test

Description

The Sedov blast (Sedov 1959) is a spherical explosion caused by the point injection of energy and in astronomy is representative of a supernova. The properties of the blast can be described by an analytic solution making it a simple test of hydrodynamic codes.

The Test

A box of side length 6 kpc containing a smooth distribution of gas (gamma=5/3) of density 1.24e7 M_sun/kpc³ (0.5 m_H/cm³) has 6.78e46 J of energy added to the centre. This corresponds to the energy of 1 supernova per 120 M_sun of gas for one particle at a resolution of 32³.

Examples of this initial condition in GADGET format are given for 32³, 64³ and 128³ particles with the energy added to the central 1, 8 and 64 particles respectively in a top-hat distribution. Units correspond to [L]=kpc, [M]=M_sun and [v]=km/s. This is only an example of a possible ic file, please feel free to create your own. Note, these ICs have zero temperatures; either set these to some minimum value, or use MinGasTemp in the parameter file.

After 10, 20 and 30 Myr the shock front should be at 0.8995, 1.1864 and 1.3958 kpc from the blast centre without cooling. The analytic solution for 10 Myr, 20 Myr and 30 Myr is in these files with the columns corresponding to radius (kpc), temperature (K), density (M_sun/kpc³) and radial velocity (km/s)

Results

Please return results for 0, 10, 20 and 30 Myr in either GADGET format or in ascii with the following 10 columns in this order: id,x,y,z,v_x,v_y,v_z,temperature,density,mass
id - particle or cell id
x,y,z - x,y,z position (or cell centre) in kpc
v_x,v_y,v_z - velocity of gas in x,y,z in km/s
temperature - temperature in K
density - density in M_sun/kpc³
mass - particle or cell mass in M_sun

Code: GADGET-2
User: Stuart Muldrew
IC: GADGET file given above
Data: 10 Myr, 20 Myr and 30 Myr in GADGET format for 128³ particles.
Plots: Denstiy evolution and density, velocity and temperature at 30 Myr.
Comments: 32 neighbours and ArtBulkViscCons of 2.0 used to obtain cleaner solution.

Code: Hydra 4.2
User: Stuart Muldrew
IC: GADGET file given above
Data: 10 Myr, 20 Myr and 30 Myr in ascii format for 128³ particles.
Plots: Denstiy evolution and density, velocity and temperature at 30 Myr.
Comments:

Code: GADGET-3
User: Chris Short
IC: GADGET file given above
Data: 10 Myr, 20 Myr and 30 Myr in GADGET format for 128³ particles.
Plots: Denstiy evolution and velocity and temperature at 30 Myr.
Comments: 64 neighbours and ArtBulkViscCons of 2.0.

Discussion forum for this suite

Recent changes RSS feed Creative Commons License Donate Driven by DokuWiki