  We conducted global hydrodynamic simulations of protoplanetary disk evolution
with an adaptive Shakura-Sunyaev {\alpha} prescription to represent the layered
disk structure, and starting with the collapse phase of the molecular cloud.
With the canonical values of model parameters, self-consistent dead zones
formed at the scale of a few au. The instabilities associated with the dead
zone and corresponding outbursts, similar to FUor eruptions, were also observed
in the simulations.
