Stability of sub-Alfvénic plasma expansions