Outflows from young stellar objects: bringing numerical simulations closer to observations of Herbig-Haro objects