Green’s functions effectively simulate shear waves produced by an applied acoustic radiation force in elastic and viscoelastic soft tissue. In an effort to determine the optimal parameters for these simulations, the convergence of Green's function-based calculations is evaluated for realistic spatial distributions of the initial radiation force “push.” The input to these calculations is generated by FOCUS, the “Fast Object-oriented C + + Ultrasound Simulator,” which computes the approximate intensity fields generated by a Phillips L7-4 ultrasound transducer array for both focused and unfocused beams. The radiation force in the simulation model, which is proportional to the simulated intensity, is applied for 200 μs, and the resulting displacements are calculated with the Green's function model. Simulation results indicate that, for elastic media, convergence is achieved when the intensity field is sampled at roughly one-tenth of the wavelength of the compressional component that delivers the radiation force “push.” Aliasing and oscillation artifacts are observed in the model for an elastic medium at lower sampling rates. For viscoelastic media, spatial sampling rates as low as two samples per compressional wavelength are sufficient due to the low-pass filtering effects of the viscoelastic medium. [Supported in part by NIH Grants R01EB012079 and R01DK092255.]