The present work is concerned with the development of a simulation model for shell and coil heat exchangers with first application in water-chillers. The model considers the heat exchanger as divided into two zones: two-phase region (boiling) and superheating region. Different refrigerant properties and heat transfer mechanisms are thus taken into account for each zone. For prescribed heat exchanger's geometry and water and efrigerant inlet conditions, a program in FORTRAN calculates the evaporator's performance which includes: refrigerant and water outlet conditions, evaporator's thermal capacity, and tube distribution for each zone. The REFPROP code programmed in FORTRAN was used for getting the local refrigerant properties, and the new kattan et al. s model for refrigerants boiling, including the new flow pattern map, was used for predicting the internal convective coefficient of heat transfer at the boiling region. A simple correlation was obtained for predicting the water convective coefficient depending on Reynolds number for the specific geometry at tested ranges of temperature. An evaporator was manufactured and an accurate test rig that included two controlled temperature baths and a data acquisition system was used for obtaining used data. The model is able to deal with a number of pure refrigerants and refrigerant blends.