Optimization of a user-defined fractal combustion model and its application to the assessment of the behaviour of a heavy-duty ng engine equipped with VVA under steady-state and transient conditions