There are actually a few aspects of that expression I don't understand, in fact I don't think they are right. It sounds like the general simplification is the "Stromgren sphere" type of simplification, where we imagine a very sharp boundary between the ionized gas and the non-ionized gas. That boundary is also a surface of force imbalance, because the ionized gas is both hotter and higher density (it just got ionized, so H turns into p+e). The high pressure inside the HII region presses on the essentially zero-pressure outside region, piling up a shell of mass that is moving supersonically relative to the external gas. The rate of change of momentum of that shell is the left-hand side, where note both its mass and velocity can change.
The right-hand side analyzes the momentum deposition into the shell. As the ionization front eats into that shell, a jetlike flow at speed u moves back into the HII region, and the recoil from that jet pushes on the shell (the second term on RHS). Also, there is the pressure from the HII region (the first term on the RHS). It all makes sense, but there seem to be some discrepancies.
First look at the pressure term, the first term on the RHS. That term basically comes from particles in the HII region bouncing off the shell, where c is the sound speed in the HII region. But the recoil of those particles will not be c, it will be more like c-rdot, just as is there in the second term. Also, the rate they encounter the shell will also be like c-rdot, not c. Presumably c >> u is implied, so they choose to worry about u-rdot but not c-rdot. Maybe that's OK, though I'm a little surprised that c >> u, so I would have thought you'd need both effects or neither.
A potentially more serious problem, it seems to me, is that they are assuming the same density in the HII region as in the jetlike flow coming off from the shell. What's the justification for that? The shell should be at higher density since it is being compressed. What's more, the u seems hard to connect with anything you'd actually see in the HII region, as it would quickly interact presumably (astronuc's question about u being a variable). Why wouldn't they want to characterize the recoil by the ionization rate instead? i would think you'd have an ionization front eating into the shell, controlled by how many ionizing photons are leaking out of the expanding HII region. If you knew that, then you'd know the rate that ionization was eating into the shell, and you know the recoil due to each ionization because the electron pretty much ends up moving backward into the HII region at something like c. That term would be quite independent of the density of the HII region.
Put differently, imagine there was no HII region at all, just a hollow void surrounded by a shell. You'd still have the ionizing photons leaking into the shell, and it would still cause ionization and recoil, but there rho=0. So I don't understand the second term on the RHS, it doesn't look right to me. I like the paper anyway though.