A growing sample of white dwarfs (WDs) with metal-enriched atmospheres are accompanied by excess infrared emission, indicating that they are encircled by a compact dusty disk of solid debris. Such `WD debris disks' are thought to originate from the tidal disruption of asteroids or other minor bodies, but the precise mechanism(s) responsible for transporting matter to the WD surface remains unclear, especially in those systems with the highest inferred metal accretion rates dM_Z/dt ~ 1e8-1e10 g/s. Here we present global time-dependent calculations of the coupled evolution of the gaseous and solid components of WD debris disks. Solids transported inwards (initially due to PR drag) sublimate at tens of WD radii, producing a source of gas that accretes onto the WD surface and viscously spreads outwards in radius, where it overlaps with the solid disk. If the aerodynamic coupling between the solids and gaseous disks is sufficiently strong (and/or the gas viscosity sufficiently weak), then gas builds up near the sublimation radius faster than it can viscously spread away. Since the rate of drag-induced solid accretion increases with gas density, this results in a runaway accretion process, during which the WD accretion rate reaches values orders of magnitude higher than can be achieved by PR drag alone. We explore the evolution of WD debris disks across a wide range of physical conditions and calculate the predicted distribution of observed accretion rates dM_Z/dt, finding reasonable agreement with the current sample. Although the conditions necessary for runaway accretion are at best marginally satisfied given the minimal level of aerodynamic drag between circular gaseous and solid disks, the presence of other stronger forms of solid-gas coupling---such as would result if the gaseous disk is only mildly eccentric---substantially increase the likelihood of runaway accretion.

Recent infrared observations have revealed presence of compact (radii < R_Sun) debris disks around more than a dozen of metal-rich white dwarfs (WD), likely produced by tidal disruption of asteroids. Accretion of high-Z material from these disks may account for the metal contamination of these WDs. It was previously shown using local calculations that the Poynting-Robertson (PR) drag acting on the dense, optically thick disk naturally drives metal accretion onto the WD at the typical rate \dot M_{PR} \approx 10^8 g/s. Here we extend this local analysis by exploring global evolution of the debris disk under the action of the PR drag for a variety of assumptions about the disk properties. We find that massive disks (mass > 10^{20} g), which are optically thick to incident stellar radiation inevitably give rise to metal accretion at rates \dot M > 0.2\dot M_{PR}. The magnitude of \dot M and its time evolution are determined predominantly by the initial pattern of the radial distribution of the debris (i.e. ring-like vs. disk-like) but not by the total mass of the disk. The latter determines only the disk lifetime, which can be several Myr or longer. Evolution of an optically thick disk generically results in the development of a sharp outer edge of the disk. We also find that the low mass (< 10^{20} g), optically thin disks exhibit \dot M << \dot M_{PR} and evolve on characteristic timescale \sim 10^5-10^6 yr, independent of their total mass.