A numerical model of four-wave mixing between picosecond pulses in a tensile-strained bulk semiconductor optical amplifier is presented. The model utilizes a modified Schrödinger equation to model pulse propagation. The Schrödinger equation parameters such as the material gain first and second order dispersion, linewidth enhancement factors and optical loss coefficient were obtained using a previously developed steady-state model. The predicted four-wave mixing pulse characteristics show reasonably good agreement with experimental pulse characteristics obtained using Frequency Resolved Optical Gating. In particular simulations predict a large increase in the pulse pedestals of the FWM converted pulse, which is verified by experiment. © 2008 Springer Science+Business Media, LLC.