A steady-state model was developed which gave detailed account of the thermal-hydraulic phenomena in down-comer-riser-drum loop of a natural circulation drum-type boiler. The model was composed of mass, energy and momentum conservation equations, together with algebraic relationships to describe the working fluids properties and phase change. The empirical relative distribution of heat flux along the height of the furnace was used to drive the model, the averaged heat flux in furnace was solved as unknown, and this fairly accords with the physical process. The model is capable of predicting working fluids temperature, pressure, circulation flow-rates and steam distributions in risers under various operating conditions. Results of the simulation provide a clear insight into two-phase flow and heat transfer conditions in evaporating tubes, and predicated circulation ratios are consistent with results obtained in other literatures.