A nonisothermal compositional model has been developed for examining nonaqueous phase liquid contamination and remediation scenarios. The governing mass balance equations and constraints have been presented, and various type of compositional formulations available have been examined. An efficient and robust formulation has been developed that addressed certain issues related to groundwater situations and overcomes numerical difficulties encountered with previous formulations. The proposed formulation collapses to that corresponding to the multiphase flow equation set when interface mass transfer is neglected. Numerical implementation of the formulation has been discussed, and example problems have been presented for benchmark and verification. The present formulation outperforms other currently used formulations for the isothermal case studied and performs at least as well as its counterparts for the nonisothermal case. Characteristics of the simulator that hinder convergence for difficult problems have been identified, and directions for even further improvements have been noted. phase mass transfer, can also be used to evaluate the effectiveness of gas venting and air sparging as feasible remedial measures for contaminated sites. Finally, models that also include nonisothermal effects can examine the effectiveness of processes such as hot water flooding and steam injection as viable remediation schemes. Permeability ko, m 2 Porosity •b Soil compressibility, Pa -1 Standard fluid molar densities Pw, mat m -3 ion, mat m-3 Pa, mat m -3 Component molecular weights Mw, kg mat -1 Mn, kg mat -1 Ma, kg mat-1 Component campressibitities