Massive stars may form in or be captured into active galactic nuclei (AGN) disks. Recent 1D studies employing stellar-evolution codes have demonstrated the potential for rapid growth of such stars through accretion up to a few hundred solar masses. We perform 3D radiation hydrodynamic simulations of moderately massive stars’ envelopes in order to determine the rate and critical radius R crit of their accretion process in an isotropic gas-rich environment in the absence of luminosity-driven mass loss. We find that in the “fast-diffusion” regime where characteristic radiative diffusion speed c/τ is faster than the gas sound speed c s , the accretion rate is suppressed by feedback from gravitational and radiative advection energy flux, in addition to the stellar luminosity. Alternatively, in the “slow-diffusion” regime where c/τ < c s , due to adiabatic accretion, the stellar envelope expands quickly to become hydrostatic and further net accretion occurs on thermal timescales in the absence of self-gravity. When the radiation entropy of the medium is less than that of the star, however, this hydrostatic envelope can become more massive than the star itself. Within this subregime, the self-gravity of the envelope excites runaway growth. Applying our results to realistic environments, moderately massive stars (≲100M ⊙) embedded in AGN disks typically accrete in the fast-diffusion regime, leading to a reduction of steady-state accretion rate 1–2 orders of magnitudes lower than expected by previous 1D calculations and R crit smaller than the disk scale height, except in the opacity window at temperature T ∼ 2000 K. Accretion in slow diffusion regime occurs in regions with very high density ρ ≳ 10−9 g cm−3, and needs to be treated with caution in 1D long-term calculations.