Source: http://people.bu.edu/robinc/kzk
Dated: 1 September 1999

BurgersTX  version 1.1

  This set of files constitutes the time-domain computer code developed at 
The University of Texas at Austin to model the propagation of 
one-dimensional nonlinear sound waves in fluids.  The code is based on an 
augmented Burgers equation that accounts for nonlinearity, thermoviscous 
absorption, and absorption and dispersion due to an arbitrary number of 
relaxation phenomena, geometrical spreading, and a stratified 
(inhomogeneous) medium.  The principal references for this code are:

R. O. Cleveland, M. F. Hamilton, and D. T. Blackstock,  "Time-domain
modeling of finite-amplitude sound in relaxing fluids," J. Acoust. Soc.
Am., Vol. 99, 3312-3318 (1996).

R. O. Cleveland, "Propagation of sonic booms through a real stratified 
atmosphere," Ph.D. dissertation, The University of Texas at Austin (1995).

  A variant of this code (THOR) was used to model sonic boom propagation 
in the atmosphere:

Robin O. Cleveland, James P. Chambers, Henry E. Bass, Richard
Raspet, Mark F. Hmailton, David T. Blackstock, "Comparison
of computer codes for the propagation of sonic booms through
the atmosphere", J. Acoust. Soc. Am., vol 97, pp 906--917 (1995).

  The algorithm for the Burgers code was based on code to solve the KZK
equation (the KZK code is also available from this web site).   The 
references for the KZK code are:
 
Y.-S. Lee, "Numerical solution of the KZK equation for pulsed finite
amplitude sound beams in thermoviscous fluids," Ph.D. dissertation, The
University of Texas at Austin (1993).

Y.-S. Lee and M. F. Hamilton, "Time-domain modeling of pulsed
finite-amplitude sound beams," J. Acoust. Soc. Am., Vol. 97, pp. 906-917
(1995).


  The files that should be in this bundle are:

BurgersTX.f  The source code

RelaxFig1    Input file to generate Fig 1 from the Cleveland et al JASA paper

S2F10_4096   The data file with the initial pressure waveform for the above 
             input file.  The file describes a step-function with a total
			 number of samples of 4096 and with 10 points in the step.

RelaxFig2    Input file to generate Fig 2 from the Cleveland et al JASA paper.  

Nwave        Input file which shows the evolution of an N-wave.

N10_4096     The data file with the initial N-wave for the above input file.
             The file has 4096 points and there are 10 points each in the
			 head and tail shocks.