ISENTROPE: calculate isentrope through a point

Overview

ISENTROPE integrates mechanical work along an isentropic compression or expansion path, using the scalar equation of state to calculate pressure. The equation of state is taken from the ARIADNE library.

A graphical interface to the program can be used on some systems.

Operation

From graphical interface:

   electra isentrope &

Direct execution: (we recommend using the expandec utility to expand input from multiple files with comments)

   expandec < isen.in | isentrope.x > isen.out

Graphical interface

Intuitive (?)

The title is a pull-down menu giving various options for running the program.
Looks for command-line options and a `.isentrope.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

All read from standard input.

equation_of_state (read to object of class (eos *), so includes type)
rho0 e0:	initial state (density, specific internal energy)
comp2:	final compression
nd:	number of densities to calculate states at on isentrope
acc:	integration accuracy, e.g. 1.0e-6
eps:	numerical bandwidth, e.g. 1.0e-8
verbose:	flag, 1 for diagnostic output
nsugsteps:	suggested integration steps between each output density, e.g. 10
maxsteps:	maximum integration steps to allow
mode:	output mode, graph or dump
abscissa:	state parameter for abscissa of graph (see below for list)
ordinate:	state parameter for ordinate of graph (see below for list)
upoffs:		offset to add to calculated particle speeds
upmult:		factor to scale calculated particle speeds (usually 1 or -1)

Recognised state parameters for abscissa and ordinate:

rho	density
v/v0	compression (1 at starting density)
v	specific volume
e	specific internal energy
p	pressure
up	particle speed
t	temperature

Output

Written in ASCII to standard output. If the "dump" option is selected, the output comprises:

# comment line identifying columns
For each shocked state (also calculated first for the unshocked state):
   mass_density
   specific_volume
   particle_speed
   pressure
   specific_internal_energy
   compression (v/v₀)
   temperature
   square_of_sound_speed

If the "graph" option is selected, only the requested columns are printed.

Files

isentrope.cc - source code
sample input file