The article discusses the implementation of two-dimensional axisymmetric configurations for the assessment of pressurized thermal shock (PTS). It highlights the evaluation of new implemented models through Neptune computational fluid dynamics (CFD) code. It presents the best performing models for both test cases.