Fast simulation of second harmonic ultrasound fields