  We report simulation results for turbulent magnetic reconnection obtained
using a newly developed Reynolds-averaged magnetohydrodynamics model. We find
that the initial Harris current sheet develops in three ways, depending on the
strength of turbulence: laminar reconnection, turbulent reconnection, and
turbulent diffusion. The turbulent reconnection explosively converts the
magnetic field energy into both kinetic and thermal energy of plasmas, and
generates open fast reconnection jets. This fast turbulent reconnection is
achieved by the localization of turbulent diffusion. Additionally, localized
structure forms through the interaction of the mean field and turbulence.
