Abstract
New FDTD algorithm is proposed for analyzing ultrasonic pulse propagation in the human body, the problem being connected with ESWL (Extracorporeal Shock Wave Lithotripsy). In this method, we do not use plane wave approximation but employ directly the original equations taking account of Lagrangian to derive new FDTD algorithms. This method is applied to an experimental setup and its numerical model that resemble actual treatment situation to compare sound pressure distributions obtained numerically with those obtained experimentally. It is shown that the present method gives clearly better results than the earlier method, in the viewpoint of numerical reappearance of strongly nonlinear waveform.
