We run three-dimensional numerical simulations for the accretion flow around the white dwarf (WD) in the progenitor system of Tycho's supernova (SN). The mass of the WD, mass of the companion star, and the orbital period are set to be 1$M_{\odot}$, 1.6$M_{\odot}$, and 0.794 day, respectively, based on theoretical and observational researches of Tycho's SN remnant (SNR). We find that when the magnetic field in the accreted material is negligible, outflowing wind is concentrated near the equatorial plane. When the magnetic field has energy equipartition with internal energy, polar wind is comparable with the equatorial wind. A carefully chosen magnetic field between the above two cases ($B=5.44\times10^3 \rm{G}$) can roughly reproduce the latitude-dependent wind required to form the peculiar periphery of Tycho's SNR.

We apply our two-dimensional (2D), radially self-similar steady-state accretion flow model to the analysis of hydrodynamic simulation results of supercritical accretion flows. Self-similarity is checked and the input parameters for the model calculation, such as advective factor and heat capacity ratio, are obtained from time-averaged simulation data. Solutions of the model are then calculated and compared with the simulation results. We find that in the converged region of the simulation, excluding the part too close to the black hole, the radial distribution of azimuthal velocity $v_\phi$, density $\rho$ and pressure $p$ basically follows the self-similar assumptions, i.e. they are roughly proportional to $r^{-0.5}$, $r^{-n}$, and $r^{-(n+1)}$, respectively, where $n\sim0.85$ for the mass injection rate of $1000L_\mathrm{E}/c^2$, and $n\sim0.74$ for $3000L_\mathrm{E}/c^2$. The distribution of $v_r$ and $v_\theta$ agrees less with self-similarity, possibly due to convective motions in the $r\theta$ plane. The distribution of velocity, density and pressure in $\theta$ direction obtained by the steady model agrees well with the simulation results {within the calculation boundary of the steady model}. Outward mass flux in the simulations is overall directed toward polar angle of 0.8382 rad ($\sim 48.0^\circ$) for $1000L_\mathrm{E}/c^2$, and 0.7852 rad ($\sim 43.4^\circ$) for $3000L_\mathrm{E}/c^2$, and $\sim$94\% of the mass inflow are driven away as outflow, while outward momentum and energy fluxes are focused around the polar axis. Part of these fluxes lie in the region that are not calculated by the steady model, and special attention should be paid when the model is applied.

We investigate the optical variability of 7658 quasars from SDSS Stripe 82. Taking advantage of a larger sample and relatively more data points for each quasar, we estimate variability amplitudes and divide the sample into small bins of redshift, rest-frame wavelength, black hole mass, Eddington ratio and bolometric luminosity respectively, to investigate the relationships between variability and these parameters. An anti-correlation between variability and rest-frame wavelength is found. The variability amplitude of radio-quiet quasars shows almost no cosmological evolution, but that of radio-loud ones may weakly anti-correlate with redshift. In addition, variability increases as either luminosity or Eddington ratio decreases. However, the relationship between variability and black hole mass is uncertain; it is negative when the influence of Eddington ratio is excluded, but positive when the influence of luminosity is excluded. The intrinsic distribution of variability amplitudes for radio-loud and radio-quiet quasars are different. Both radio-loud and radio-quiet quasars exhibit a bluer-when-brighter chromatism. Assuming that quasar variability is caused by variations of accretion rate, the Shakura-Sunyaev disk model can reproduce the tendencies of observed correlations between variability and rest-frame wavelength, luminosity as well as Eddington ratio, supporting that changes of accretion rate plays an important role in producing the observed optical variability. However, the predicted positive correlation between variability and black hole mass seems to be inconsistent with the observed negative correlation between them in small bins of Eddington ratio, which suggests that other physical mechanisms may still need to be considered in modifying the simple accretion disk model.

In order to study the outflows from accretion disks, we solve the set of hydrodynamic equations for accretion disks in the spherical coordinates ($r\theta\phi$) to obtain the explicit structure along the $\theta$ direction. Using self-similar assumptions in the radial direction, we change the equations to a set of ordinary differential equations (ODEs) about the $\theta$-coordinate, which are then solved with symmetrical boundary conditions in the equatorial plane, and the velocity field is obtained. The $\alpha$ viscosity prescription is applied and an advective factor $f$ is used to simplify the energy equation.The results display thinner, quasi-Keplerian disks for Shakura-Sunyaev Disks (SSDs) and thicker, sub-Keplerian disks for Advection Dominated Accretion Flows (ADAFs) and slim disks, which are consistent with previous popular analytical models. However, an inflow region and an outflow region always exist, except when the viscosity parameter $\alpha$ is too large, which supports the results of some recent numerical simulation works. Our results indicate that the outflows should be common in various accretion disks and may be stronger in slim disks, where both advection and radiation pressure are dominant. We also present the structure dependence on the input parameters and discuss their physical meanings. The caveats of this work and possible improvements in the future are discussed.

We show that there exists a maximal possible accretion rate, beyond which global slim disk solutions cannot be constructed because in the vertical direction the gravitational force would be unable to balance the pressure force to gather the accreted matter. The principle for this restriction is the same as that for the Eddington luminosity and the corresponding critical accretion rate, which were derived for spherical accretion by considering the same force balance in the radial direction. If the assumption of hydrostatic equilibrium is waived and vertical motion is included, this restriction may become even more serious as the value of the maximal possible accretion rate becomes smaller. Previous understanding in the literature that global slim disk solutions could stand for any large accretion rates is due to the overestimation of the vertical gravitational force by using an approximate potential. For accretion flows with large accretion rates at large radii, outflows seem unavoidable in order for the accretion flow to reduce the accretion rate and follow a global solution till the central black hole.