KZKTexas2 Focused Beams
This codes solves the KZK equation on a uniform grid and is
configured for focused sources.
winkzk2.zip Zip file with all
the necessary C++ code for a Windows platform.
- Read me file
- kzktexas2.f Fortran source
- GetCols.f Fotran code that
pulls the ascii output file "pulse.out" apart so that
individual waveforms can be viewed
- kzktexas2.cpp C++ source
- read_kzk.m Matlab file that
opens up the binary output file from the cpp code and plots the
- yslFig4.9 Input file to generate
Fig 4.9 from the Lee's thesis. This is the focusing of a pulse
with nonlinear effects. The original figure is available in jpg
format Figure 4.9.
- yslFig5.25 Input file to generate
Fig 5.25 from Yang Sub Lee's thesis. This is focusing of a pulse
WITHOUT nonlinear distortion. The original figure is available
in jpg format Figure 5.25.
- yslFig5.26input Input
file to generate Fig 5.26 from Yang Sub Lee's thesis except in
this case the initial condition is read from a file rather than
generated by the code itself. The original figure is available
in jpg format Figure 5.26.
- G5dp100tp100 The file with
the initial pressure waveform for the above input file. This
is provide to demonstrate how to generate an arbitrary initial
waveform. The cryptic name stands for a gain of 5, 100 rho points
and 100 tau points.
Specific references to this code include:
- Michalakis A. Averkiou and Mark F. Hamilton, "Nonlinear
distortion of short pulses radiated by plane and focused circular
pistons", J. Acoust. Soc. Am., Vol. 102, pp 2539-2548 (1997).
- Michalakis A. Averkiou and Robin O. Cleveland, "Time
domain numerical modeling of the Dornier HM3 lithotripter,"
J. Acoust. Soc. Am. 106, pp 102-112 (1999).
Texas Time Domain page
RETURN to Robin's