This paper presents time-dependent solutions to the coupled dusty hydrodynamics equations describing the spherically symmetric expansion of cometary neutral gas in the vicinity of a cometary nucleus. The sublimation process is repressented by gas outflow from a dust-covered reservoir containing stationary gas whose pressure and density values are determined by the sublimating (Ts) and surface (T0) temperatures. The model resolves earlier ambiguities in determining gas production rates and provides analytic relations between Ts, T0, and the gas parameters at the sonic point. The time evolution of a cometary outburst was modeled. It was found that, as a result of the strong gas-dust interaction in the inner coma region, a 'slow' disturbance in both the dust and gas parameters will be created in addition to the familiar gas blast-wave solution. This new 'slow' disturbance, which propagates with a velocity of about 0.2 km/s, might be responsible for some of the observed slowly expanding cometary halos, such as the one which was recently identified using 1910 Mount Wilson high-resolution comet Halley photographs.