Phase changes in multistream heat exchangers can be detected with a nonsmooth model.

Flash calculations with single-phase outlets can be handled with a nonsmooth model.

Complex simulations of a natural gas liquefaction process can be performed.

Simulations do not require the solution of a difficult optimization problem.

Automatic differentiation outperforms finite differencing in the solution procedure.