A multi‐Gaussian beam model is described where the velocity field on the face of an ultrasonic immersion transducer is represented in terms of a large number of localized Gaussians. The versatility of this approach is demonstrated by using this type of model to represent both focused and planar probes, rectangular and circular probe shapes, and piston as well as non‐piston behavior. © 2004 American Institute of Physics