We have developed a global fully three-dimensional magnetohydrodynamic solar wind model for the region that extends from the coronal base to 100 AU. The simulation domain consists of tree spherical shell subdomains with computational boundaries between them placed at 20 solar radii and 0.3 AU. The location of the first boundary ensures that the flow at the boundary is both supersonic and super-Alfvenic. A steady-state solution in the innermost (coronal) region is obtained by the time-relaxation method. The solution uses a tilted dipole model or solar magnetograms as the boundary condition at the coronal base and includes a flux of Alfven waves in the WKB approximation which provide additional acceleration for the coronal outflow in the open field regions. The intermediate region solution is constructed by the integration of steady-state equations along radius using a marching scheme. The outer region solution (0.3-100 AU) is obtained again by the time relaxation and takes into account turbulence transport and heating as well as heating, flow deceleration, and other effects due to the interstellar pickup protons treated as a separate fluid. We use the model to simulate the global steady-state structure of the solar wind from the coronal base to the heliospheric boundary and compare the results with Ulysses and Voyager observations
