Surface figuring using chemically reactive plasma jet machining (PJM) is a promising non-conventional technology for deterministic ultra-precision machining of optical components. Based on chemical reactions between plasma generated radicals and the surface atoms this technology is capable to fabricate complex shaped free form surfaces. Since the material removal rate during PJM depends strongly on the surface temperature which itself is influenced by the jet heat flux to the surface, the arising nonlinear effects on the etch result have to be regarded. Conventionally applied dwell time calculation algorithms do not consider those effects leading to significant machining errors in some cases. In order to improve the machining procedure with respect to deterministic material removal yielding predictable results a process simulation model has been developed. This model considers spatio-temporal variations of surface temperature and temperature dependent material removal and is able to predict the final workpiece topography after machining.