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.