A simulation model of brazed plate heat exchangers (BPHE) operating in steady-state in a cascade vapor compression refrigeration system has been developed. For this system the heat exchangers were the condenser, intermediate heat exchanger or cascade-condenser and evaporator. The model adopts a local analysis method, where the heat exchanger is divided into a series of small control volumes, to which the heat transfer rate equations and the fundamental of conservation of mass and energy equation. Local heat exchanger coefficients and friction factor are calculated using correlations available in literature, considering regions of single-phase or twophase flow for each one of the heat exchangers. These correlations have heat flux values ranging from 2,5kW/m2 to 185kW/m2, saturation temperatures from 5 degrees Celsius to 35 degrees Celsius and were applied to geometries with corrugation angle ranging from 20 degrees to 60 degrees. In order to calculate the thermal performance of the heat exchangers and the output conditions of the two fluids a computational program was developed in FORTRAN. This simulation considers known inlet conditions of the fluids and the geometry of the heat exchanger. The thermophysical properties of the refrigerants fluids were calculated using the version 7 of REFPROP, a package by NIST (National Institute of Standards and Technology), that allow for the simulation of heat exchangers with a wide operating range of refrigerants. The simulation results were compared with experimental data (condenser and evaporator) for R22 and R290 refrigerants, obtaining a good agreement. A sensibility analysis for heat exchangers, using the new R1234yf and R1234ze has also been carried out.