functions in demo1.i -
demo1
|
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
|
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
|
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
|
sound Set up the initial conditions for evolve to launch a weak sound wave. | |
SEE ALSO: | shock, evolve |