KZKTexas2 Focused Beams

This codes solves the KZK equation on a uniform grid and is configured for focused sources.

• Read me file
• kzktexas2.f Fortran source code
• GetCols.f Fotran code that pulls the ascii output file "pulse.out" apart so that individual waveforms can be viewed
• kzktexas2.cpp C++ source code
• read_kzk.m Matlab file that opens up the binary output file from the cpp code and plots the individual waveforms.
• 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).

KZK Texas Time Domain page

RETURN to Robin's home page.