As part of the Shallow Water Acoustics in a Random Medium (SWARM) experiment, a 16 element WHOI vertical line array (WVLA) was moored in 70 m of water off the New Jersey coast. A 400-Hz acoustic tomography source was moored some 32-km shoreward of this array, such that an acoustic path was created that was anti-parallel to the primary propagation direction for shelf-generated internal wave solitons. The presence of these soliton internal waves in the acoustic waveguide causes significant coupling of energy between propagating acoustic modes, creating fluctuations in modal intensities and modal peak arrival times, as well as time spreading of the pulses. Two methods by which acoustic propagation and scattering in soliton-filled waveguides can be modeled are presented here in order to understand and explain the scattering observed in the SWARM field data. The first method utilizes the Preisig and Duda [IEEE J. Ocean. Eng. 22, 256-269 (1997)] Sudden Interface Approximation (SIA) to represent the solitons. The second method, which is computationally slower, uses a finely meshed, "propagated" thermistor record to simulate the solitons in the SWARM experiment waveguide. Both numerical methods are found to generate scattering characteristics that are similar to the SWARM field data.