The “smooth” quantum hydrodynamic (QHD) model is derived specifically to handle in a mathematically rigorous way the discontinuities in the classical potential energy which occur at heterojunction barriers in quantum semiconductor devices. Smooth QHD model simulations of the resonant tunneling diode are presented which exhibit enhanced negative differential resistance when compared with simulations using the original O(ħ2) QHD model. In addition, smooth QHD simulations of a classical electron shock wave are presented which agree with classical hydrodynamic model simulations and which do not exhibit the spurious dispersive oscillations of the O(ħ2) QHD model.