Heavy ion (O+, O++, and N+) density enhancements in the outer plasmasphere have been observed using the retarding ion mass spectrometer instrument on the DE 1 satellite. These are seen at L shells from 2 to 5, with most occurrences in the L = 3 to 4 region; the maximum L shell at which these enhancements occur varies inversely with Dst. It is also known that enhancements of O+ and O++ overlie ionospheric electron temperature peaks. It is thought that these enhancements are related to heating of plasmaspheric particles through interactions with ring current ions. This was investigated using a time‐dependent one‐stream hydrodynamic model for plasmaspheric flows, in which the model flux tube is connected to the ionosphere. The model simultaneously solves the coupled continuity, momentum, and energy equations of a two‐ion (H+ and O+) quasi‐neutral, currentless plasma. This model is fully interhemispheric and diffusive equilibrium is not assumed; it includes a corotating tilted dipole magnetic field and neutral winds. First, diurnally reproducible results were found assuming only photoelectron heating of thermal electrons. For this case the modeled equatorial O+ density was below 1 cm−3 throughout the day. The O+ results also show significant diurnal variability, with standing shocks developing when production stops and O+ flows downward under the influence of gravity. Numerical tests were done with different levels of electron heating in the plasmasphere; these show that the equatorial O+ density is highly dependent on the assumed electron heating rates. Over the range of integrated plasmaspheric electron heating (along the flux tube) from 8.7 to 280×109 eV/s, the equatorial O+ density goes like the heating raised to the power 2.3.