  We studied the final phases of galactic mergers, focusing on interactions
between supermassive black holes (SMBHs) and the interstellar medium in a
central sub-kpc region, using an N-body/hydrodynamics code. This numerical
experiment aims to understand the fate of the gas supplied by mergers of two or
more galaxies with SMBHs, whose masses are $10^7 M_\odot$. We observed that the
mass accretion rate to one SMBH exceeds the Eddington accretion rate when the
distance between two black holes (BHs) rapidly decreases. However, this rapid
accretion phase does not last for more than $10^7$ yrs, and it drops to $\sim$
10% of the Eddington rate in the quasi-steady accretion phase. The rapid
accretion is caused by the angular momentum transfer from the gas to the
stellar component, and the moderate accretion in the quasi-steady phase is
caused {by turbulent viscosity and gravitational torque in the disk. The second
merger event enhances the mass accretion to the BHs; however, this phase takes
place on a similar timescale to the first merger event. We also found that the
AGN feedback and the mass accretion to BHs can coexist in the central region of
merged galaxies, if the amount of feedback energy is given as $(2 \times
10^{-4} - 2 \times 10^{-3} )\dot{M} c^2$, where $\dot{M}$ is the accretion rate
to $r= 1$ pc. The accretion rate is suppressed by $\sim$ 1/50 in the
quasi-steady accretion phase for $0.02 \dot{M} c^2$. The fraction of the gas
that finally falls to each BH is approximately 5-7% of the supplied total gas
mass ($10^8 M_\odot$), and 15-20% of the gas forms a circumnuclear gas inside
100 pc. This remnant gas heavily obscures the luminous phase of the active
galactic nuclei (AGN) during merger events, and the moderate AGN feedback does
not alter this property.
