To understand why supercritical accretion is feasible onto a neutron star, we carefully examine the accretion flow dynamics by 2.5-dimensional general relativistic (GR) radiation magnetohydrodynamic (RMHD) simulations, comparing the cases of accretion onto a non-magnetized neutron star (NS) and that onto a black hole (BH). Supercritical BH accretion is relatively easy, since BH can swallow excess radiation energy, so that radiation flux can be inward in its vicinity. This mechanism can never work for NS which has a solid surface. In fact, we find that the radiation force is always outward. Instead, we found significant reduction in the mass accretion rate due to strong radiation-pressure driven outflow.The radiation flux F rad is self-regulated such that the radiation force balances with the sum of gravity and centrifugal forces. Even when the radiation energy density much exceeds that expected from the Eddington luminosity E rad ≃ F rad τ /c > 10 2 L Edd /(4πr 2 c), the radiation flux is always kept below the certain value which makes it possible not to blow all the gas away from the disk. These effects make supercritical accretion feasible. We also find that a settling region, where accretion is significantly decelerated by radiation cushion, is formed around the NS surface. In the settling region, the radiation temperature and mass density roughly follow T rad ∝ r −1 and ρ ∝ r −3 , respectively. No settling region appears around the BH so that matter can be directly swallowed by the BH with supersonic speed.