  We present direct numerical simulations of inhomogeneous reduced
magnetohydrodynamic (RMHD) turbulence between the Sun and the Alfv\'en critical
point. These are the first such simulations that take into account the
solar-wind outflow velocity and the radial inhomogeneity of the background
solar wind without approximating the nonlinear terms in the governing
equations. RMHD turbulence is driven by outward-propagating Alfv\'en waves
($z^+$ fluctuations) launched from the Sun, which undergo partial non-WKB
reflection to produce sunward-propagating Alfv\'en waves ($z^-$ fluctuations).
We present ten simulations with different values of the correlation time
$\tau_{\rm c\,\sun}^+$ and perpendicular correlation length $L_{\perp \sun}$ of
outward-propagating Alfv\'en waves (AWs) at the coronal base. We find that
between 15% and 33% of the $z^+$ energy launched into the corona dissipates
between the coronal base and Alfv\'en critical point. Between 33% and 40% of
this input energy goes into work on the solar-wind outflow, and between 22% and
36% escapes as $z^+$ fluctuations through the simulation boundary at $r=r_{\rm
A}$. The $z^\pm$ power spectra scale like $k_\perp^{-\alpha^\pm}$, where
$k_\perp$ is the wavenumber in the plane perpendicular to $\vec{B}_0$. In our
simulation with the smallest value of $\tau_{\rm c\,\sun}^+$ ($\sim 2
\mbox{min}$) and largest value of $L_{\perp \sun}$ ($2\times 10^4 \mbox{km}$),
we find that $\alpha^+$ decreases approximately linearly with increasing
$\ln(r)$, reaching a value of 1.3 at $r=11.1 R_{\sun}$. Our simulations with
larger values of $\tau_{\rm c\,\sun}^+$ exhibit alignment between the contours
of constant $\phi^+$, $\phi^-$, $\Omega_0^+$, and $\Omega_0^-$, where
$\phi^\pm$ are the Els\"asser potentials and $\Omega_0^\pm$ are the outer-scale
parallel Els\"asser vorticities.
