A temperature feasibility constraint is an important part of multistream heat exchanger (MSHE) modeling. However, temperature feasibility of an MSHE model makes a numerical issue when a physical property package is used to obtain highly accurate temperature-enthalpy relationships in equation-oriented modeling environment. To resolve the issue, this study proposes a new MSHE model with enthalpy feasibility using the fact that enthalpy is a monotonically increasing function of temperature. A natural gas liquefaction process, called a single mixed refrigeration process, is optimized using the proposed MSHE model under an equation-oriented modeling environment with a physical property package as a case study. (c) 2018 Elsevier Ltd. All rights reserved.