Global solutions of optically thick advective accretion disks around black holes are constructed. The solutions are obtained by solving numerically a set of ordinary differential equations corresponding to a steady axisymmetric geometrically thin disk. We pay special attention to consistently satisfy the regularity conditions at singular points of the equations. For this reason we analytically expand a solution at the singular point, and use coefficients of the expansion in our iterative numerical procedure. We obtain consistent transonic solutions in a wide range of values of the viscosity parameter alpha and mass acretion rate. We compare two different form of viscosity: one takes the shear stress to be proportional to the pressure, while the other uses the angular velocity gradient-dependent stress. We find that there are two singular points in solutions corresponding to the pressure-proportional shear stress. The inner singular point locates close to the last stable orbit around black hole. This point changes its type from a saddle to node depending on values of alpha and accretion rate. The outer singular point locates at larger radius and is always of a saddle-type. We argue that, contrary to the previous investigations, a nodal-type inner singular point does not introduce multiple solutions. Only one integral curve, which corresponds to the unique global solution, passes simultaneously the inner and outer singular points independently of the type of inner singular point. Solutions with the angular velocity gradient-dependent shear stress have one singular point which is always of a saddle-type and corresponds to the unique global solution. The structure of accretion disks corresponding to both viscosities are similar.

We consider the effects of advection and radial gradients of pressure and radial drift velocity on the structure of optically thick accretion disks. We concentrate our efforts on highly viscous disks, $\alpha=1.0$, with large accretion rates. Contrary to disk models neglecting advection, we find that continuous solutions extending from the outer disk regions to the inner edge exist for all accretion rates we have considered. We show that the sonic point moves outward with increasing accretion rate, and that in the innermost disk region advection acts as a heating process that may even dominate over dissipative heating. Despite the importance of advection on it's structure, the disk remains geometrically thin.

We present solutions to the accretion disk structure equations in which advective cooling is accounted for in a self-consistent way. It is shown that for high rates of accretion, above a critical value, when without advection there are no solutions extending continuously from large to small radii, there exist global solutions with self-consistent advective cooling. These solutions are quite different fromthose solutions obtained earlier by other authors, where advective cooling was not treated self-consistently.

We present here the solution for accretion disk structure, which describes continuously the transition between optically thick and optically thin disk regions. We show that the disk structure equations without advection, used here, give two branches of solutions which do not intersect for $L<L_{b}\lesssim 0.6 L_{\rm Edd}$ for $\alpha=1.0$ and $M_{BH}=10^8 M_\odot$. For larger luminosities there are no global solutions of the equations without advection. We suggest, as one of the possibilities, that advection becomes important even before, at luminosities $\sim0.1L_{\rm Edd}$, when the inner disk regions become unstable due to radiation pressure dominance and induce the transition from optically thick to optically thin disk in the global solution.