Computer Program for Heating and Cooling Loads in Buildings

In the late 1960s, under the sponsorship of the Department of Housing and Urban Development, Tamami Kusuda began to develop a computer program to help architects and engineers predict the thermal performance of a building. Important enough at the time, this work was to balloon in significance in the context of the energy shortages and rapidly escalating energy prices of the 1970s. The computer program was named the National Bureau of Standards Load Determination (NBSLD) program [1,2]. It combined algorithms for transient conduction in the building structure, solar heat gains and radiant transfer, and convection between building surfaces and the room air to allow the prediction of temperatures and heating and cooling loads under dynamic conditions. Moreover, the algorithms employed in NBSLD were adopted by the American Society of Heating, Refrigerating and Air-Conditioning Engineers (ASHRAE) as a recommended procedure for computerized energy calculations. They were presented to the Society in a special publication [3] which was sold for many years.

NBSLD simulated a very simple building model because of the computer limitations of its time. The model included the walls, windows, ceiling, and floor of a single space which might represent a room or an entire building. The user described the geometry of this space, the thermal properties of the building surfaces, optical properties of windows, air infiltration, heat sources such as equipment and occupants, and weather data, includ-ing temperature, humidity, and solar radiation values for periods from one day to an entire year. The program converted the building data into the coefficients of a set of simultaneous equations which were solved for temperatures of all the building components as they responded to the internal and external heat sources over time. By fixing the space air temperature to a desired value, the solution of the simultaneous equations gave the amount of energy that had to be added or removed, i. e., heating or cooling loads, necessary to maintain that temperature.

NBSLD was the result of many years of leadership in this field by Kusuda. His earlier work included develop-ment of numerical methods to quickly compute transient conduction through building walls [4] and to determine the thermal properties of moist air [5]. He organized the first three symposia on the Use of Computers for Environmental Engineering Related to Buildings. The first was held at NBS in 1970 [6], the second in Paris in 1974, and the third in Banff, Canada, in 1978. NBSLD was developed with the cooperation of ASHRAE's Task Group on Energy Requirements, in particular the Subcommittee on Heating and Cooling Load Calculations which Kusuda chaired. The members of the subcommittee shared their specialized work and computer algorithms. Kusuda developed and imple-mented the NBSLD framework for putting the algorithms together, and, together with NBS staff, filled in the missing components of the complete model. The Bureau's experimental facilities were used for validation of the program [7,8].

Kusuda's work with the Task Group on Energy Requirements and subsequent technical committees led to his recognition as a fellow of ASHRAE, the society's highest honor. In 1969 the Automated Procedures Engineering Consultants (APEC) awarded Kusuda its first Honorary Membership to recognize his pioneering work on the use of computers for heating, cooling and ventilating engineers.

NBSLD and related research were applied to the energy crisis of the 1970s in the analysis of energy conservation options for specific buildings [9] and in the development of general guidelines [10,11]. Kusuda's work continued in related areas, such as using sunlight in place of electric lighting [12] and studying heat exchange to the ground from buildings and pipes [13]. Recognition for this work includes a Lifetime Achievement Award from the International Building Performance Simulation Association and Silver and Gold Medals from the Department of Commerce.

In addition to the research and energy conservation design done with the NBSLD program, the documenta-tion of its algorithms provided the basis for the develop-ment of other building energy simulation programs, both public and private. The algorithms were used almost immediately in the U. S. Army Construction Engineering Research Laboratory's (CERL) first build-ing energy analysis program [14]. In the following years CERL developed new program code to simulate multi-ple rooms and allow a greater degree of integration of the building loads with the air handling system and plant equipment such as boilers and air conditioners. CERL's BLAST program [15] was released in 1977. The gradual development of multizone airflow calculations has led to a tool [16] for studying indoor air quality and even designing smoke control systems.

At present, the Department of Energy is developing the EnergyPlus program [17], which merges the best features of BLAST and the DOE-2 program. In particu-lar, it is using the heat balance [18] concepts from BLAST which were derived from NBSLD. Because of the need to gather the latest such procedures in a single document, as was originally done by Kusuda, ASHRAE is currently working to provide improved algorithms based on a quarter century of research since the publication of the NBSLD algorithms.

Tamami Kusuda was born in Seattle, Washington. His family returned to Japan before World War II and lived there during the war (and the bombing). He received his B. S. degree from Tokyo University in 1947 and his M. S. and Ph. D. degrees from the University of Minnesota in 1952 and 1955.

Prepared by George Walton.

Bibliography

[1] Tamami Kusuda, NBSLD, the Computer Program for Heating and Cooling Loads in Buildings, Building Science Series 69, National Bureau of Standards, Washington, DC (1976).

[2] T. Kusuda, NBSLD, Computer Program for Heating and Cooling Loads in Buildings, NBSIR 74-574, National Bureau of Standards, Washington, DC (1974).

[3] ASHRAE Task Group on Energy Requirements, Procedure for Determining Heating and Cooling Loads for Computerizing Energy Calculations: Algorithms for Building Heat Transfer Subroutines, American Society of Heating, Refrigerating and Air-Conditioning Engineers, New York (1975).

[4] T. Kusuda, Thermal Response Factors for Multi-Layer Structures of Various Heat Conduction Systems, ASHRAE Trans. 75 (Part 1), 246-271 (1969).

[5] T. Kusuda, Algorithms for Psychrometric Calculations, (Skeleton Tables for the Thermodynamic Properties of Moist Air) , Building Science Series 21, National Bureau of Standards, Washington, DC (1970).

[6] T. Kusuda (ed.), Use of Computers for Environmental Engineer-ing Related to Buildings, Building Science Series 39, National Bureau of Standards, Washington, DC (1971).

[7] Bradley A. Peavy, Frank J. Powell, and Douglas M. Burch, Dynamic Thermal Performance of an Experimental Masonry Building, Building Science Series 45, National Bureau of Stan-dards, Washington, DC (1973).

[8] B. A. Peavy, D. M. Burch, F. J. Powell, and C. M. Hunt, Comparison of Measured and Computer-Predicted Thermal Performance of a Four Bedroom Wood-Frame Townhouse, Building Science Series 57, National Bureau of Standards, Washington, DC (1975).

[9] Tamami Kusuda, James E. Hill, Stanley T. Liu, James P. Barnett, and John W. Bean, Pre-Design Analysis of Energy Conservation Options for a Multi-Story Demonstration Office Building, Building Science Series 78, National Bureau of Standards, Washington, DC (1975).

[10] G. Kelly, T. Kusuda, and J. Hill, Potential for Energy Conserva-tion in Heating, Ventilating, and Air-conditioning Equipment for Buildings, in Energy Conservation Through Effective Energy Utilization, Jesse C. Denton, Stephen Webber, and John Moriarty, (eds.), NBS Special Publication 403, National Bureau of Standards, Washington, DC (1976) pp. 163-191.

[11] T. Kusuda, Use of Simulation Models of Buildings in Assessing Energy Conservation Strategies, in Energy Conservation Strate-gies in Buildings: Comfort, Acceptability, and Health, Jan A. J. Stolwijk (ed.), John B. Pierce Foundation, New Haven, CT (1978).

[12] T. Kusuda and J. W. Bean, Comparison Between a Simplified Daylighting Calculation Procedure and a Comprehensive Inter-reflection Model Calculation Procedure, in Proceedings of the Third International Symposium on the Use of Computers for Environmental Engineering Related to Buildings, Alberta, Canada, CNRC Publ. 17376, National Research Council of Canada, Ottawa, Canada (1978).

[13] T. Kusuda, M. Mizuno, and J. W. Bean, Seasonal Heat Loss Calculation for Slab-on-Grade Floors, NBSIR 81-2420, National Bureau of Standards, Washington, DC (1982).

[14] D. Hittle and B. Sliwinski, CERL Thermal Loads Analysis and Systems Simulation Program Volume 2: Program Reference Manual, CERL Interim Report E-81, US Army Construction Engineering Research Laboratory, Champaign, Illinois (1975).

[15] D. C. Hittle, The Building Loads Analysis and System Thermo-dynamics (BLAST) Program, CERL Technical Report E-119, US Army Construction Engineering Research Laboratory, Champaign, Illinois (1977).

[16] George N. Walton, CONTAM96 User Manual, NISTIR 6056, National Institute of Standards and Technology, Gaithersburg, MD (1997).

[17] D. B. Crawley, L. K. Lawrie, C. O. Pedersen, R. J. Liesen, D. E. Fisher, R. K. Strand, R. D. Taylor, F. C. Winkelmann, W. F. Buhl, A. E. Erdem, and Y. J. Huang, EnergyPlus: A New-Generation Building Energy Simulation Program, in Building Simulation ' 99, Volume I, Kyoto, Japan, International Building Performance Simulation Association (1999) pp. 81-88.

[18] R. Strand, F. Winkelmann, F. Buhl, J. Huang, R. Liesen, C. Pedersen, D. Fisher, R. Taylor, D. Crawley, and L. Lawrie, EnergyPlus: Enhancing and Extending the Capabilities of the Building Heat Balance Simulation Technique for Use in Energy-Plus, in Building Simulation ' 99, Volume I, September 1999, Kyoto, Japan, International Building Performance Simulation Association (1999) pp. 653-660.

Fig. 1. Analogous electric circuit for heat exchange process in a room.

Fig. 2. Physical model of a typical room represented by the circuit in Fig. 1.