Source:
http://people.bu.edu/robinc/kzk

This is set of files contains the time domain computer code developed at U 
Texas Austin, for solving the KZK equation for the case of a focused 
source.  The code accounts for nonlinearity, diffraction, thermoviscous 
absoprtion, and multiple relaxation frequencies.  Refer to the following 
articles for the details of the code:

Yang Sub-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, December 1993.

Yang Sub-Lee and Mark F. Hamilton, "Time-domain modeling of finite-
amplitude beams",  J. Acoust. Soc. Am., vol 97, pp 906--917 (1995).

Robin O. Cleveland, Mark F. Hamilton, and David T. Blackstock,  "Time-domain 
modeling of finite-amplitude waves in relaxing fluids", 
J. Acoust. Soc. Am. 99, 3312--3318 (1996).

Results from a code almost identical to this code have also been presented:

Michalakis A. Averkiou and Mark F. Hamilton, "Nonlinear distortion of short 
pulses radiated by plane and focused circular pistons", J. Acoust.  Soc.  
Am., 102, pp 2539-2548 (1997).


  The files that should be in this bundle are:

kzktexas2.f  the source code

GetCols.f   program that pulls the output file "pulse.out" apart so that 
            individual waveforms can be examined.

yslFig4.9   input file to generate Fig 4.9 from the Lee's dissertation

yslFig5.25  input file to generate Fig 5.25 from the Lee's dissertation

yslFig5.26i input file to generate Fig 5.26from the Lee's dissertation 
            except it reads the inital waveform in from a file rather 
            than generating it itself.

G5dp100tp100   The file with the initial pressure waveform for the above 
               input file. The cryptic name stands for Gain=5, 
               data points = 100,  time points=100.