Ondol heating, a kind of radiant floor heating, is a main method used in housing units in Korea. Building energy simulation including ondol and relevant facilities has not been performed due to its complexity. For evaluating energy consumption and indoor temperature variation, a new method should be proposed. At the present work, a dynamic simulation on ondol heating was tried by combining TRNSYS and EES. Characteristic functions for a pump, hot water coils and a gas boiler were simultaneously solved by EES, and calculated flow rates and supply temperature of hot water were provided as inputs of the active layer of TYPE 56 in TRNSYS. The results by the simulation on a typical housing unit in Korea shows a good trend in a viewpoint of actual behavior of ondol heating.