This paper presents the simulation models of the plate-fin, air-to-water (or water vapour) heat exchangers used as air-heating or air-cooling and dehumidifying coils in the HVAC (Heating, Ventilation and Air-Conditioning) systems. The thermal models are used to calculate the heat exchange between distributing air and coil pipes and outlet temperatures of air and heat or chilled fluid. The aerodynamic models are used to account for the pressure drop of the air crossing the coil tubes. They can also be used to optimize the structures of such coils. The models are based on principal laws of heat and mass conservation and fluid mechanics. They are transparent and easy to use.In our work, a coil is considered as an assembly of numbers of basic elements in which all the state variables are unique. Therefore we can conveniently simulate the coils with different structures and different geometric parameters. Two modular programs TRNSYS (Transient System Simulation) and ESACAP are utilized as supporting softwares which make the programming and simulation greatly simplified. The coil elements and a real coil were simulated. The results were compared with the data offered by the manufacturer (company SOFICA) and also with those obtained using critical methods such as NTU method, etc. and good agreement is attained. This paper presents the simulation models of the plate-fin, air-to-water (or water vapor) heat exchangers used as air-heating or air-cooling and dehumidifying coils in the HVAC (Heating, Ventilation and Air-Conditioning) systems. The thermal models are used to calculate the heat exchange between distributing air and coil pipes and outlet temperatures of air and heat or chilled fluid. The aerodynamic models are used to account for the pressure drop of the air crossing the coil tubes. They can also be used to optimize the structures of such coils. The models are based on principal laws of heat and mass conservation and fluid mechanics. They are transparent and easy to use. Our work, a coil is considered as an assembly of numbers of basic elements in which all the state variables are unique. Thus we can conveniently simulate the coils with different structures and different geometric parameters. Two modular programs TRNSYS (Transient System Simulation) and ESACAP are utilized as s The results were compared with the data offered by the manufacturer (company SOFICA) and also with those obtained using critical methods such as NTU method, etc . and good agreement is attained.