Standard Library

demo1

DOCUMENT demo1
  run the 1-D hydrocode demonstration.  An X window should pop up
  in which a movie of a wave propagating down a shock tube appears.
  Use the 'sound' and 'shock' commands to set two different interesting
  initial conditions; the default is 'sound'.

SEE ALSO: sound, shock, evolve

evolve

DOCUMENT evolve, time1
      or evolve, time1, time0
  Step the hydro calculation forward to TIME1,
  starting with the initial conditions in the
  RT, z, and v arrays at time TIME0 (default 0.0
  if omitted).  The calculation also depends on
  the constants M (mass/area/zone) and gammaM1
  (gamma-1 for the gamma-law equation of state).
  The pressure array p is updated in addition to
  the primary state arrays RT, z, and v.

  Boundary conditions are specified by setting
  either a boundary pressure or a boundary
  velocity at each end of the fluid column.
  bc0_v   - Boundary velocity at z(0), or []
            if z(0) has pressure BC.
  bc0_p   - Boundary pressure beyond z(0).
  bc0_time  - If bc0_v or bc0_p is an array,
              bc0_time is an array of the same
              length specifying the corresponding
              times for time dependent BCs.
  bc0_Z   - Acoustic impedance at z(0) if bc0_v
            is nil (default is 0).
  bcN_v, bcN_p, bcN_time, and bcN_Z have the same
  meanings for the z(n_zones) boundary.

  The worker routines OutputResults and
  TakeStep must be supplied.

shock

DOCUMENT sound
  Set up the initial conditions for evolve to launch a strong wave, which
  steepens into a shock as it propagates.

SEE ALSO: sound, evolve

sound

DOCUMENT sound
  Set up the initial conditions for evolve to launch a weak sound wave.

SEE ALSO: shock, evolve