Hi I am also trying to build a refrigerant cycle similar to the turbine system, and i am also using the DAE equations directly without the MTK components. When i reduce my DAE to steady state model, the NL solver can solve the system. But when i tried to solve the same steady state model using IDA solver, the solver failed to initialize even when i fed in the solution generated from NL solver as the initial condition. I am thinking this might be related to the coolprop functions i called. Is there any suggestions for building the DAE system using coolprop and include phase change scenario? I have set up this question DAE‑safe bounding for CoolProp property calls in Julia that contains the same problem. Thank you for your reply!