dump_shock_nemd

In shock wave piston simulations, it’s often crucial to compute thermo information at different regions, both before and after the shock wave passage.

Piston simulations commonly involve millions of atoms. Dumping all the virial and velocity data for each atom can lead to excessively large output files, making data processing cumbersome. The dump_shock_nemd command addresses this by calculating spatial thermo information during the simulation.

This feature calculates the spatial distribution of partical velocity (km/h), stress (GPa), temperature and density (g/cm3) in x-direction.

Syntax

dump_shock_nemd interval <time_interval> bin_size <size of each bin>

  • The interval parameter sets the output interval (number of steps).

  • The bin_size parameter, optional with a default value of 10, defines the thickness of each histogram bin, measured in Angstroms.

Examples

To output spatial thermo information every 1000 steps for a run in the x-direction, with a bin size of 20 Angstroms, include the following before the run keyword:

dump_shock_nemd interval 1000 bin_size 20