Numerical simulations are used to investigate soliton-like propagation in tapered silicon core optical fibers. The simulations are based on a realistic tapered structure with nanoscale core dimensions and a decreasing anomalous dispersion profile to compensate for the effects of linear and nonlinear loss. An intensity misfit parameter is used to establish the optimum taper dimensions that preserve the pulse shape while reducing temporal broadening. Soliton formation from Gaussian input pulses is also observed--further evidence of the potential for tapered silicon fibers to find use in a range of signal processing applications.