  The stationary hydrodynamic equations for the transonic accretion disks andflows around rotating black holes are presented by using the Kerr-Schildcoordinate where there is no coordinate singularity at the event horizon. Weuse two types of the causal viscosity prescription, and the boundary conditionsfor the transonic accretion flows are given at the sonic point. For one type ofthe causal viscosity prescription we also add the boundary conditions at theviscous point where the accreting radial velocity is nearly equal to theviscous diffusion velocity. Based on the formalism for the transonic accretiondisks, after we present the calculation method of the transonic solutions, thehorizon-penetrating transonic solutions which smoothly pass the event horizonare calculated for several types of the accretion flow models: the idealisothermal flows, the ideal and the viscous polytropic flows, the advectiondominated accretion flows (ADAFs) with the relativistic equation of state, theadiabatic accretion disks, the standard accretion disks, the supercriticalaccretion disks. These solutions are obtained for both non-rotating androtating black holes. The calculated accretion flows plunge into black holewith finite three velocity smaller than the speed of light even at the eventhorizon or inside the horizon, and the angular velocities of the accretion flowat the horizon are generally different from the angular velocity of theframe-dragging due to the black hole's rotation. These features contrast to theresults obtained by using the Boyer-Lindquist coordinate with the coordinatesingularity at the horizon.
We use two types of the causal viscosity prescription, and the boundary conditions for the transonic accretion flows are given at the sonic point. For one type of the causal viscosity prescription we also add the boundary conditions at the viscous point where the accreting radial velocity is nearly equal to the viscous diffusion velocity.