Rare-earth-doped optical fibers are one of the most promising solid-state lasers. In these fiber lasers, a cladding-pumping scheme using double-clad fibers is utilized to increase the overall conversion efficiency of pumping light. To maintain acceptable beam quality, the low-numerical aperture large-mode-area fibers is effective for the double-clad fibers because the effects of stimulated Raman scattering can be reduced via the corresponding reduction in the power density in the large fiber core. For the large-mode-area double-clad fibers, fiber fuse propagation was investigated theoretically by the explicit finite-difference method using the thermochemical SiO x production model. In the calculation, we assumed the fiber to be in an atmosphere and that part (40 µm in length) of the core was heated to a temperature of 2,923 K. The threshold power for the double-clad fiber with the core radius of 10 µm was 1.6 W at 1.080 µm and it was close to the experimental value. The power dependence of the velocity of fiber fuse propagation was calculated for the double-clad fibers with the core radius of 10 and 15 µm. The calculated velocities were in fair agreement with the experimental values observed in the input power range from 1 kW to 3.5 kW at 1.080 µm.