A physical mechanism that drives FU Orionis-type outbursts is reconsidered. We study the effect of inner part of a circumstellar disk covering a region from near the central star to the radius of approximately $5$ AU (hereafter, the inner disk). Using the fluctuated mass accretion rate onto the inner disk $\dot{M}_{\rm out}$, we consider the viscous evolution of the inner disk and the time variability of the mass accretion rate onto the central star $\dot{M}_{\rm in}$ by means of numerical calculation of an unsteady viscous accretion disk in a one-dimensional axisymmetric model. First, we calculate the evolution of the inner disk assuming an oscillating $\dot{M}_{\rm out}$. It is shown that the time variability of $\dot{M}_{\rm in}$ does not coincide with $\dot{M}_{\rm out}$ due to viscous diffusion. Second, we investigate the properties of spontaneous outbursts with temporally constant $\dot{M}_{\rm out}$. Outburst occur only in a limited range of mass accretion rates onto the inner disk $10^{-10}<\dot{M}_{\rm out}< 3\times 10^{-6}~{\rm M}_{\odot} {\rm yr}^{-1}$ due to gravo-magneto limit cycle (GML). Finally, we discuss the case with a combination of episodic $\dot{M}_{\rm out}$ and accretion outbursts cause by the GML in the inner disk. The GML can drive accretion outbursts onto the star even for the case of fluctuating $\dot{M}_{\rm out}$, although fluctuations of $\dot{M}$ decay during transmitting the inner disk inwards. We newly identified two modes of outburst which are spontaneous one and stimulated one. In a stimulated mode of outburst, $\dot{M}_{\rm out}$ does appear directly in $\dot{M}_{\rm in}$ (the latter defining the stellar accretion luminosity). In a spontaneous mode of outburst, $\dot{M}_{\rm out}$ appears as the interval between outbursts.

The growing process of both a young protostar and a circumstellar disk is investigated. Viscous evolution of a disk around a single star is considered with a model where a disk increases its mass by dynamically accreting envelope and simultaneously loses its mass via viscous accretion onto the central star. We focus on the circumstellar disk with high mass accretion rate onto the disk $\dot{M}=8.512c_{\rm s}^3/G$ as a result of dynamical collapse of rotating molecular cloud core. We study the origin of the surface density distribution and the origin of the disk-to-star mass ratio by means of numerical calculations of unsteady viscous accretion disk in one-dimensional axisymmetric model. It is shown that the radial profiles of the surface density $\Sigma$, azimuthal velocity $v_{\phi}$, and mass accretion rate $\dot{M}$ in the inner region approach to the quasi-steady state. Profile of the surface density distribution in the quasi-steady state is determined as a result of angular momentum transport rather than its original distribution of angular momentum in the cloud core. It is also shown that the disk mass becomes larger than the central star in the long time limit as long as temporary constant mass flux onto the disk is assumed. After the mass infall rate onto the disk declines owing to the depletion of the parent cloud core, the disk-to-star mass ratio $M_{\rm disk}/M_*$ decreases. The disk-to-star mass ratio becomes smaller than unity after $t> 10^5 \rm yr$ and $t>10^6 \rm yr$ from the beginning of the accretion phase in the case with $\alpha_0 =1 {\rm and} 0.1$, respectively, where $\alpha_0 $ is the constant part of viscous parameter. In the case with $\alpha_0 \leq 10^{-2}$, $M_{\rm disk}/M_*$ is still larger than unity at $2 \rm Myr$ from the beginning of the accretion phase.