Abstract: | A new FDTD algorithm is proposed for analyzing ultrasonic pulse propagation in the human body, a problem connected with ESWL (extracorporeal shock wave lithotripsy). In this method, we do not use plane wave approximation but directly employ the original equations, taking into account the Lagrangian to derive new FDTD algorithms. This method is applied to an experimental setup and its numerical model that resembles an actual treatment situation in order to compare the sound pressure distributions obtained numerically with those obtained experimentally. The present method is shown to give clearly better results than the earlier method, from the point of view of the numerical reappearance of a strongly nonlinear waveform. © 2009 Wiley Periodicals, Inc. Electr Eng Jpn, 169(4): 29–36, 2009; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/eej.20819 |