We propose a new method, emission-line eclipse mapping, to map the velocity fields of an accretion disk in position space. Quiescent dwarf novae usually exhibit double-peaked emission line profiles because of disk rotation. Since a part of the disk having a different line-of-sight velocity is successively obscured by a companion in eclipsing systems, they show time-varying line profiles. We calculate time changes of the emission line profiles, assuming Keplerian rotation fields (v_\phi\propto r^(-1/2) with r being the distance from the disk center) and the emissivity distribution of j \propto r^(-3/2). We, then, apply the usual eclipse mapping technique to the light curves at each of 12-24 wavelengths across the line center to map the region with the same line-of-sight velocity. The reconstructed images typically exhibit a `two-eye' pattern for high line-of-sight velocities and we can recover the relation, v_\phi \propto d^(-1/2), on the assumption of the axisymmetric disk, where d is the separation between the two `eyes'. We will be able to probe the Keplerian rotation law, the most fundamental assumption adopted in many disk models, by high-speed spectroscopic observations with 8-m class telescopes.

Intermediate resolution phase-resolved spectra of WZ Sge were obtained on five consecutive nights (July 23 -- 27) covering the initial stage of the 2001 superoutburst. Double-peaked emission lines of He\textsc{II} at 4686 \AA, which were absent on July 23, emerged on July 24 together with emission lines of C\textsc{III} / N\textsc{III} Bowen blend. Analyses of the He\textsc{II} emission lines using the Doppler tomography revealed an asymmetric spiral structure on the accretion disk. This finding demonstrates that spiral shocks with a very short orbital period can arise during the initial stage of an outburst and may be present in all SU UMa stars.

The hydrodynamic behavior of an accretion disk in a close binary system is numerically simulated. Calculation is made for a region including the compact star and the gas-supplying companion. The equation of state is that of an ideal gas characterized by the specific heat ratio $\gamma$. Two cases with $\gamma$ of 1.01 and 1.2 are studied. Our calculations show that the gas, flowing from the companion via a Lagrangian L1 point towards the accretion disk, forms a fine gas beam (L1 stream), which penetrates into the disk. No hot spot therefore forms in these calculations. Another fact discovered is that the gas rotating with the disk forms, on collision with the L1 stream, a bow shock wave, which may be called an L1 shock. The disk becomes hot because the L1 shock heats the disk gas in the outer parts of the disk, so that the spiral shocks wind loosely even with $\gamma=1.01$. The L1 shock enhances the non-axisymmetry of the density distribution in the disk, and therefore the angular momentum transfer by the tidal torque works more effectively. The maximum of the effective $\alpha$ becomes $\sim 0.3$. The 'hot spot' is not formed in our simulations, but our results suggest the formation of the 'hot line', which is the L1 shock elongated along the penetrating L1 stream.

The history of hydrodynamic numerical simulations for accretion disks in close binary systems is reviewed, in which emphasis is placed, in particular, on the facts that spiral shock waves were numerically found in 1986 by researchers including one of the present authors and that spiral structure was discovered in IP Pegasi in 1997 by Steeghs et al. The results of our two and three-dimensional numerical simulations in recent years are then summarized, with comparison being made with observations.

We perform 2D and 3D numerical simulations of an accretion disc in a close binary system using the Simplified Flux vector Splitting (SFS) finite volume method. In our calculations, gas is assumed to be the ideal one, and we calculate the cases with gamma=1.01, 1.05, 1.1 and 1.2. The mass ratio of the mass-losing star to the mass-accreting star is one. Our results show that spiral shocks are formed on the accretion disc in all cases. In 2D calculations we find that the smaller gamma is, the more tightly the spiral winds. We observe this trend in 3D calculations as well in somewhat weaker sense. Recently, Steeghs, Harlaftis & Horne (1997) found the first convincing evidence for spiral structure in the accretion disc of the eclipsing dwarf nova binary IP Pegasi using the technique known as Doppler tomography. We may claim that spiral structure which we have discovered in numerical simulations before are now found observationally.

The first convincing piece of evidence of spiral structure in the accretion disc in IP Pegasi was found by Steeghs et al. (1997). We performed two kinds of 2D hydrodynamic simulations, a SFS finite volume scheme and a SPH scheme, with a mass ratio of 0.5. Both results agreed well with each other. We constructed Doppler maps and line flux-binary phase relations based on density distributions, the results agreeing well with those obtained by observation.

We perform 2D and 3D numerical simulations of an accretion disc in a close binary system using the Simplified Flux vector Splitting (SFS) finite volume method. In our calculations, gas is assumed to be the ideal one, and we calculate the cases with gamma=1.01, 1.05, 1.1 and 1.2. The mass ratio of the mass losing star to the mass accreting star is unity. Our results show that spiral shocks are formed on the accretion disc in all cases. In 2D calculations we find that the smaller gamma is, the more tightly the spiral winds. We observe this trend in 3D calculations as well in somewhat weaker sense.