Numerical values of the NRTL equation parameters for calculation of the vapour - liquid - liquid equilibria (VLLE) at atmospheric pressures have been presented for 5 ternary mixtures. These values were fitted to the experimental VLLE and vapour - liquid equilibrium (VLE) data to describe simultaneously, as accurately as possible, the VLE and the liquid - liquid equilibria (LLE). The coefficients of this model called further NRTL-VLL were used for simulations of n-propanol dehydration via heterogeneous azeotropic distillation. The calculations performed by a ChemCAD simulator were done for 4 mixtures using hydrocarbons, ether and ester as an entrainer. In majority simulations the top streams of the azeotropic column had composition and temperature similar to the corresponding experimental values of ternary azeotropes. The agreement between the concentrations of both liquid phases formed in a decanter and the experimental values of the LLE was good for all four simulations. The energy requirements were the most advantageous for the simulation with di-npropyl ether (DNPE) and isooctane. Simulations were performed also for one mixture using the NRTL equation coefficients taken from the ChemCAD database. In that case the compositions of the liquid organic phases leaving the decanter differed significantly from the experimental LLE data.