Natural gas liquefaction is an energy intensive process where the feed is cooled from ambient temperature down to cryogenic temperatures. Different liquefaction cycles exist depending on the application, with dual mixed refrigerant processes normally considered for the large-scale production of Liquefied Natural Gas (LNG). Large temperature spans and small temperature differences in the heat exchangers make the liquefaction processes difficult to analyze. Exergetic losses from irreversible heat transfer increase exponentially with a decreasing temperature at subambient conditions. Consequently, an accurate and robust simulation tool is paramount to allow designers to make correct design decisions. However, conventional process simulators, such as Aspen Plus, suffer from significant drawbacks when modeling multistream heat exchangers. In particular, no rigorous checks exist to prevent temperature crossovers. Limited degrees of freedom and the inability to solve for stream variables other than outlet temperatures also makes such tools inflexible to use, often requiring the user to resort to a manual iterative procedure to obtain a feasible solution. In this article, a nonsmooth, multistream heat exchanger model is used to develop a simulation tool for two different dual mixed refrigerant processes. Case studies are presented for which Aspen Plus fails to obtain thermodynamically feasible solutions.