high performance computing on graphics processing units: hgpu.org

A form of the conservation equations for fluid dynamics is presented, deduced using slightly less restrictive hypothesis than those necessary to obtain the well known Westervelt equation. This formulation accounts for full wave diffraction, nonlinearity, and thermoviscous dissipative effects. A CLAWPACK based, 2D finite volume method using the Roe linearization was implemented to obtain numerically the solution of the proposed equations. In order to validate the code, two different tests have been performed: one against a special Taylor shock-like analytic solution, the other against published results on a HIFU system, both with satisfactory results. The code is based on CLAWPACK and is written for parallel execution on a GPU, thus improving performance by a factor of over 60 when compared to the standard CLAWPACK code.