Nonequilibrium molecular dynamics simulations are applied to the investigation of step-flow kinetics at crystal-melt interfaces of silicon, modeled with the Stillinger-Weber potential ͓Phys. Rev. B 31, 5262 ͑1985͔͒.Step kinetic coefficients are calculated from crystallization rates of interfaces that are vicinals of the faceted ͑111͒ orientation. These vicinal interfaces contain periodic arrays of bilayer steps, and they are observed to crystallize in a step-flow growth mode at undercoolings lower than 40 K. Kinetic coefficients for both ͓110͔ and ͓121͔ oriented steps are determined for several values of the average step separation, in the range of 7.7-62.4 Å. The values of the step kinetic coefficients are shown to be highly isotropic, and are found to increase with increasing step separation until they saturate at step separations larger than ϳ50 Å. The largest step kinetic coefficients are found to be in the range of 0.7-0.8 m / ͑sK͒, values that are more than five times larger than the kinetic coefficient for the rough ͑100͒ crystal-melt interface in the same system. The dependence of step mobility on step separation and the relatively large value of the step kinetic coefficient are discussed in terms of available theoretical models for crystal growth kinetics from the melt.