  The Westervelt equation is a model for the propagation of finite amplitude
ultrasound. The method of discrete exterior calculus can be used to solve this
equation numerically. A significant advantage of this method is that it can be
used to find numerical solutions in the discrete space manifold and the time,
and therefore is a generation of finite difference time domain method. This
algorithm has been implemented in C++.
