On the convergence, stability and computational speed of numerical schemes for 0-D IC engine cylinder modelling