  The edge plasma turbulence and transport dynamics, as well as the divertor
power loads during the thermal quench phase of tokamak disruptions are
numerically investigated with BOUT++'s flux-driven, six-field electromagnetic
turbulence model. Here a transient yet intense particle and energy sources are
applied at the pedestal top to mimic the plasma power drive at the edge induced
by a core thermal collapse, which flattens core temperature profile.
Interesting features such as surging of divertor heat load (up to 50 times),
and broadening of heat flux width (up to 4 times) on the outer divertor target
plate, are observed in the simulation, in qualitative agreement with
experimental observations. The dramatic changes of divertor heat load and width
are due to the enhanced plasma turbulence activities inside the separatrix. Two
cross-field transport mechanisms, namely the $E\times B$ turbulent convection
and the stochastic parallel advection/conduction, are identified to play
important roles in this process. Firstly, elevated edge pressure gradient
drives instabilities and subsequent turbulence in the entire pedestal region.
The enhanced turbulence not only transports particles and energy radially
across the separatrix via $E\times B$ convection which causes the initial
divertor heat load burst, but also induces an amplified magnetic fluctuation
$\tilde{B}$. Once the magnetic fluctuation is large enough to break the
magnetic flux surface, magnetic flutter effect provides an additional radial
transport channel. In the late stage of our simulation, $|\tilde{B}_r/B_0|$
reaches to $10^{-4}$ level that completely breaks magnetic flux surfaces such
that stochastic field-lines are directly connecting pedestal top plasma to the
divertor target plates or first wall, further contributing to the divertor heat
flux width broadening.
