HUGONIOT: calculate shock Hugoniot

Files

Graphical interface:
- hugoniot - shell script to set up environment variables and run graphical interface. Filenames and options can be supplied on the command line; these are passed to hugoniot.tcl. To work correctly, the electra script should be in the path (and should be set up to point to the right files).
- hugoniot.tcl - Tcl/Tk script
- alstein.eos - sample eqn of state data
- Tcl/Tk utility library
C++ program:
- hugoniot.cc - source code
- sample input file
- Ariadne material model library V3.0: equation of state parts
- WXC++ maths library: vector and ode integration parts

Graphical interface

Intuitive (?)

The title is a pull-down menu giving various options for running the program.
Looks for command-line options and a `.hugoniot.tcl_rc' file of directory defaults (Tcl/Tk utility library standards). If neither is found, settings revert to internal defaults. These may be regeneratewd by pressing `Defaults'.
`Save defaults' creates a new .rc file.
Buttons change the display to allow different groups of parameters to be edited.
On-line help.

Input to bare program

ASCII, read from standard input.

verbose		flag, 1 for diagnostic output
material_name
material_model (ARIADNE mat_type* object, so includes type)
initial_state
rotaxis rotang	axis (3-vector) and angle (radians) for rotating material to 1D loading plane
gradu		grad(u), 3x3 matrix
comp1 comp2	range of compressions to calculate over
nd		number of compression values
upoffs		offset to add to calculated shock speeds
upmult		factor to scale calculated shock speeds (usually 1 or -1)
usoffs		offset to add to calculated particle speeds
usmult		factor to scale calculated particle speeds (usually 1 or -1)
tol		solution accuracy, e.g. 1.0e-6
geomfac		geometrical energy increase factor in solution search, e.g. 1.1
denergy0	initial energy increment in solution search, e.g. 1.0e-6
maxits		maximum iterations in finding each solution

Output

ASCII, written to standard output. Some miscellaneous data are written at the start, to confirm some of the input settings. Each of these lines starts with # so it is filtered automatically e.g. by Gnuplot. The rest of the output consists of Hugoniot states, one per line.

For initial state and each state along the Hugoniot:
   mass_density
   specific_volume
   compression
   particle_speed
   shock_speed
   specific_internal_energy
   stress (Voigt notation, i.e. 6-vector)
   temperature
   bulk_sound_speed_squared
   specific_heat_capacity
   mean_pressure
   state