AbstractSimulation of a hydrocarbons mixture separation by extractive distillation was based on binary vapor-liquid as well as liquid-liquid equilibrium data. Sulfolane was considered as extractive solvent for the selective toluene separation from a model mixture with heptane. For simulation of the chosen hydrocarbons mixture separation in the presence of an extractive solvent, the NRTL model was considered. A set of temperature-dependent binary NRTL parameters was evaluated independently by fitting the experimental vapor-liquid and liquid-liquid equilibrium data of the respective binary subsystems. In order to improve the description of the ternary vapor-liquid equilibrium, original NRTL model extended by the ternary contribution term was used. Parameters of the ternary contribution were obtained by direct fitting of available ternary liquid-liquid equilibrium data while employing the original binary NRTL parameters. Quality of the ternary vapor-liquid-liquid equilibrium description using the original and the extended excess Gibbs energy models was assessed by comparing the calculated compositions of conjugate liquid phases with experimental data. Using the extended NRTL model, mean deviation of the computed mole fractions decreased by approximately four times (9.73 × 10−3) compared to the value obtained using the original NRTL model. Both original and extended NRTL models were employed for the simulation of a model mixture separation by extractive distillation. At chosen experimental conditions, high purity distillate (x 1 > 0.999) was obtained. Results of the aromatics extractive distillation in the presence of sulfolane were compared to those obtained with N-methylpyrrolidone as the extractive solvent.