Optimisation and Cost Analysis of a Lithium Bromide Absorption Solar Cooling System

Similar documents
Life cycle analysis of thermosyphon solar water heaters

ANNUAL ENERGY PERFORMANCE OF SOLAR THERMAL SYSTEMS IN BRAŞOV, ROMANIA

Modeling and analyzing solar cooling systems in Polysun

Optimization of Thermosiphon Solar Water Heating systems working in a Mediterranean Environment of Cyprus

Renewables. Vacuum Tube Solar Systems. Solar Energy to the Power of

BUILDING INTEGRATION OF RENEWABLE ENERGY SYSTEMS TOWARDS ZERO OR NEARLY ZERO ENERGY BUILDINGS

Conversion Design of Solar Houses from Active to Passive

Solar Energy Technologies

Dimensioning a small-sized PTC solar field for heating and cooling of a hotel in Almería (Spain)

New developed solar thermal systems for heating and cooling Budapest, 16 th April 2009 Tsekouras Panagiotis Mech. Engineer NTUA Centre for Renewable

Performance Evaluation of Solar Parabolic Trough for Cloths Laundry application

Sustainable Building Façade and Advanced Fenestration Systems

Solar Air-Conditioning Systems in small - medium scale applications

Although electricity is efficient at heating water it is expensive and not available everywhere.

The impact of hot and cold storages on a solar absorption cooling system for an office building

Hot water persons. Map section

Simulation And Optimization of a Solar Absorption Cooling System Using Evacuated Tube Collectors

FLATE Hillsborough Community College - Brandon (813)

Experimental Study to Evaluate the Performance of Iraqi Passive House in Summer Season

Performance of a Solar Heating System with Photovoltaic Thermal Hybrid Collectors and Heat Pump

Thermal Characteristics of a Vernacular Building Envelope

Comfort Cooling Application Using Fixed Focus Solar Parabolic Dish Concentrator Integrated with Double Effect Vapor Absorption Machine

SOLAR COOLING WITH SMALL SIZE CHILLER: STATE OF THE ART

Performance evaluation of hybrid solar parabolic trough concentrator systems in Hong Kong

Performance Modeling of a Solar Driven Absorption Cooling System for. Carnegie Mellon University s Intelligent Workplace

LOW-ENERGY HOUSING ESTATE SUNDAYS

Solar Oven. Curriculum Levels 5-6 Science. Activity Description

Camel Solar brings you some of the most advanced solar thermal systems available on the market today through innovation, research and development.

UNIT FOUR SOLAR COLLECTORS

use Solar thermal for hot water production DESIGN APPROACHES

Future tree planting on the two acre lot will offset remaining lifetime carbon emissions.

Available online at ScienceDirect. Energy Procedia 91 (2016 )

Performance evaluation of solar absorption cooling system of Bahal (Haryana)

Operation of Swimming Pool Solar Heating at Blackmount School.

Renewable energy from the sun Daikin solar thermal systems for the installer

SOLAR ENGINEERING OF THERMAL PROCESSES

Performance Improvement in Solar Evacuated Glass Tube Collector

Experimental investigation of single-phase and twophase closed thermosyphon solar water heater systems

use Solar thermal for hot water production

EFFICIENT SOLAR COOLING

Microgeneration Installation Standard : MCS 024 SOLAR DOMESTIC HOT WATER ENERGY CALCULATION. to be used in conjunction with MIS 3001:

Héctor Rubio Plana Gas Natural, SDG S.A. 18th. June 2009

Solar Energy Today s Clean Energy Choice

Solar Cooling Technologies for Southern Climates - A System Comparison

SOLAR PHOTOVOLTAIC ENERGY

HEATING WITH A HEAT PUMP

Report On Task 57, Subtask B1: Development of Accelerated Ageing Tests for Evacuated Tube Collectors

- 2 - SME Q1. (a) Briefly explain how the following methods used in a gas-turbine power plant increase the thermal efficiency:

Building integration of solar renewable energy systems towards zero or nearly zero energy buildings

Verified net Zero Energy Building with air source heat pumps for SME

Climate responsive design

Available online at ScienceDirect. Energy Procedia 91 (2016 )

SorTech package solution description

Investigations of thermal self-sufficient Residential Buildings with Solar Energy Systems

COMPARISONS BETWEEN TRNSYS SOFTWARE SIMULATION AND F-CHART PROGRAM ON POOL HEATING SYSTEM

Research on integrated solar and geothermal energy engineering design in hot summer and cold winter area

Environmental Benefits of Domestic Solar Water Heating Systems

Indo-Swiss Building Energy Efficiency Project. Case Study: Aranya Bhawan, Jaipur

Figure 1: Uncovered and covered PVT collectors

Technical Guide to Solar Thermal. Better World Energy. How to harness the sun s energy and create hot water with a Solar Thermal System

Energy Technology & Conservation. Week_05. Instructor: Mr. Adnan Qamar. Mechanical Engineering Department

Chromasun Micro-Concentrator

Solar Absorption Aqua-Ammonia Absorption system simulation base on Climate of Malaysia

A Solar Wall System Utilized in Rural Houses of Northeast China

On the Night Heat Losses in Thermosyphon Solar Water Heaters

DOUBLE-EFFECT SOLAR ABSORPTION THERMAL ENERGY STORAGE

Microgeneration Installation Standard : MCS 024 SOLAR DOMESTIC HOT WATER ENERGY CALCULATION. to be used in conjunction with MIS 3001:

ENERGY EVALUATION BY SIMULATION FOR EFFECTIVE USE OF SEWAGE HEAT

Solar Thermal Systems- A reliable technology for Green Buildings

Solar Thermal. [Strategy]

Solar Cooling for Southern Climates, Double Effect Absorption Chillers with High Concentrating Collectors Versus Standard Single Effect Systems

Energy Resources. Assosciate Professor

The Double Glass Skin Façade: a sustainable option for the Curtain Wall in Mediterranean Climate Countries?

Gas Natural, SDG S.A. 18th. June 2009

Design and optimization of solar water heating system for a five star hotel in Varanasi

Simulation Studies on Performance of Solar Cooling System in UAE Conditions

Simulation and Optimization of Vacuum Tube Solar Collector Water Heating System in Iran

Assignment 3 Integrated Energy Design

Australian Solar Cooling Interest Group (ausscig) Conference Financial analysis of solar cooling systems in Australia

Project PolarSol koti 41\193\- Map section

EVALUATION OF THE SOLAR INCOME FOR BRAŞOV URBAN AREA

Innovative Solar thermal systems for heating and cooling

Project Report: SunMaxx Commercial Case Studies

27/10/15. Rod Hughes Architect Director 2030 architects ltd. Penrith

*Corresponding authors. Tel.: x415, address: (N. Zhu)

Continental J. Engineering Sciences 6 (1): 24-30, 2011 Wilolud Journals, 2011 ` Printed in Nigeria

Original Article. World Journal of Engineering Research and Technology WJERT.

Project Report: SunMaxx Commercial Case Studies

Project Report: SunMaxx Commercial Case Studies

COMPARATIVE STUDIES ON THERMAL EFFICIENCY OF SINGLE AND DOUBLE GLAZED FLAT PLATE SOLAR WATER HEATER

Experimental investigation of photovoltaic-thermal hybrid solar collector

SolarCity Partnership. Toronto Fire Station # kwt Solar Water Heating Installation. Technology. Monitoring. Best Practices

A Review of Solar Collectors in Solar Thermal Applications

Solar-Thermal Feasibility Study Farm #5 Boutique Winery

SOLAR ENERGY INTRODUCTION:

Tasks of the planner

BASIC EUROPEAN GLAZED COLLECTORS & SOLAR DOMESTIC

Solar-Assisted Space Heating of a Highly Insulated Energy Efficient House

Analysis of Solar Thermal Cooling System Using TRANSOL

Nearly Zero Energy Building in Lecco Modern technological building compared to an ancient Villa

Transcription:

Optimisation and Cost Analysis of a Lithium Bromide Absorption Solar Cooling System Georgios A. Florides and Soteris A. Kalogirou Mechanical Engineering Department, Higher Technical Institute, Nicosia, Cyprus. Corresponding email: SKalogirou@hti.ac.cy SUMMARY This paper presents the optimisation of the various components of a lithium bromide (LiBr) absorption solar cooling system such as the type, slope and area of solar collector and storage tank size. The collector types considered are the flat plate, compound parabolic and evacuated tube collectors. The optimisation is based on an energy benefit analysis, i.e., the amount of useful energy collected against the life cycle cost of the solar system. The above analysis considers the current prices for fuel and equipment costs, which have increased greatly during the last 5 years. A solar absorption system consisting of a lithium bromide-water unit, a solar collector and a storage tank is modelled with the TRNSYS computer program using a typical meteorological year (TMY) for a hot climate (Nicosia, Cyprus). For this analysis a typical house is modelled for a full year. The solar system can be used during summer to provide part of the heat required by the absorption unit and during winter to provide part of the heating load. The results indicate that due to the present high cost of fuel a large part of the building load can be covered with solar energy. INTRODUCTION Solar energy is in abundance in Cyprus. In summer, mean monthly temperatures for Nicosia at 14. hours in July, are 35.4 C with the temperature sometimes reaching 43 C. Therefore, there is a need to lower the indoor temperature considerably in order to be able to provide comfort. Solar cooling of buildings seems to be one of the most attractive solutions. This is an application in which the demand for cooling energy closely matches the availability of solar energy, not only to the seasonal but also to the daily variation. Many researchers have developed solar assisted absorption refrigeration systems. Most of them have been produced as experimental units and computer codes were written to simulate the systems. Relevant studies are shown in [1-4]. The objective of this paper is to model a complete system, consisting of a number of solar collectors, storage tank, a boiler and a LiBr-water absorption refrigerator, which will cover a typical house load during the whole year. Conclusions are drawn by comparing the economics of the system between the years 21 and 27 when there was a dramatic change in the fuel price. For this analysis the typical house indicated in Fig. 1 is considered. The house thermal load is minimised by considering an insulated roof, insulated walls, double glazed windows, internal shading and night ventilation (3 ach) in summer. The above factors were studied and found to be economically viable [5]. The double-walls are made of.1 m hollow brick and.2m

plaster on each side and a layer of.5 m insulation in between. The roof is constructed from fair-faced.15 m heavy concrete,.5 m polystyrene insulation,.7 m screed and.4 m asphalt, covered with aluminum paint of.55 solar absorptivity. The above construction requires an annual cooling load at 25 C of 17,6 kwh with a peak load of 1.3 kw and an annual heating load at 21 C of 3,53 kwh with a peak load of 5.5 kw. Fig. 1 Model house The TRNSYS program is used to model the complete system (house load estimation with solar powered heating and absorption cooling), together with the weather values of a typical meteorological year (TMY) file for Nicosia, Cyprus. The solar powered system consists of an array of solar collectors, boiler, storage tank, 11 kw absorption cooling unit, pumps and thermostats. Using this approach a system optimisation is performed in order to select the right equipment, i.e., the collector type, the storage tank volume and the collector slope angle and area. The collector area is decided by performing an economic analysis of the system. Also the longterm integrated system performance and the dynamic system behaviour is evaluated and an economic comparison is made between the years 21 and 27. CHARACTERISTICS OF AN 11 KW WATER-LIBR ABSORPTION CHILLER The characteristics needed in TRNSYS deck file are the generator load, the mass flow of the heating water to the generator heat exchanger and its input and output temperatures. Referring to previous work [6] the calculated 11 kw unit characteristics are as shown in Table 1. Also other parameters needed for the TRNSYS input file were calculated based on this previous work.

Table 1 Water-LiBr absorption refrigeration system calculations based on a generator temperature of 75 C and a solution heat exchanger exit temperature of 55 C Point H (kj/kg) %LiBr m (Kg/s) P (kpa) T ( C) (X) Remarks 1 83.5691.93 34.9 55 2 83.5691 4.82 34.9 55 3 124.7.5691 4.82 55 55 Sub-cooled liquid 4 183.2.5217 4.82 75 6 5 137.8.5217 4.82 51.5 6 6 137.8.5217.93 44.5 6 7 2612.474 4.82 7 Superheated Steam 8 131..474 4.82 31.5 Saturated liquid 9 131..474.93 6 1 2511.8.463.93 6 Saturated vapour 11 23.45.11.93 6 Saturated liquid Description Symbol kw Capacity (evaporator output power) Absorber heat, rejected to the environment Heat input to the generator Condenser heat, rejected to the environment Q 11. e Q 14.1 a Q 14.9 g Q 11.8 c Coefficient of performance COP.74 THE COMPLETE SYSTEM CHARACTERISTICS The complete system besides the absorption refrigerator consists of a number of solar collectors, a thermally insulated vertical storage tank, a conventional boiler and interconnecting piping. A schematic of the system showing also the simulation program information flow is shown in Fig. 2. The system was modelled with the TRNSYS simulation program. The program consists of many subroutines that model subsystem components. The type number of every TRNSYS subroutine used to model each component is also shown in Fig. 2. TYPE 11(2) Diverter TYPE 6 Absorption Refrigerator TYPE 2 Thermostat TYPE 3 Pump 2 TYPE 11(3) Mixer TYPE 38 Algebraic Tank TYPE 2 Thermostat TYPE 1 Solar Collector TYPE 19 House TYPE 3 Pump 1 Fig. 2 Circuit diagram and TRNSYS types used for modelling the system.

The construction and type of the solar collectors is important and relevant to the operation and efficiency of the whole system. Three types of solar collectors, modelled with TRNSYS Type 1, are evaluated in this study as follows: Flat Plate Collectors: These are predominantly used for domestic hot water production. The external casing of the collectors considered in this study, which are manufactured locally, is made of high corrosion resistant galvanised steel sheet, sprayed with aluzinc paint. The casing is covered with a 4 mm thick single glass and sealed with a rubber gasket. The absorber plate is made of copper. High radiation absorption is achieved by the use of black fine matt finish on the copper surface, which has a high absorption coefficient. The underside of the absorber plate and the side casing are well insulated to reduce conduction losses with 5 mm and 3 mm fibreglass insulation respectively. Each flat plate collector consist of 12 evenly spaced parallel copper pipes, 15 mm in diameter, embossed by semi-circular grooves formed in the flat plate absorber. Conventional flat plate collectors are developed for use in sunny and warm climates. During cold, cloudy and windy days their performance is greatly reduced. Furthermore, condensation and moisture can cause early deterioration of internal materials resulting in reduced performance and system failure. Compound Parabolic Concentrating Collectors (): These collectors use curved reflecting surfaces to concentrate sunlight onto a small absorber area. s are used for higher water temperature applications than the flat plate collectors. Such a focusing collector performs very well in direct sunlight but, depending on the concentration ratio, does not perform well under cloudy or hazy skies because only a few rays are captured and reflected onto the absorber. A compound parabolic collector system can work either as a stationary system, or it can track the sun. In stationary systems, much of the sunlight that hits the concentrator's reflector often misses the absorber. Tracking devices allow the collector to follow the sun's movement across the sky. This ensures that the concentrator always faces toward the sun. Such systems are usually employed for higher temperature applications. Since residential applications only require medium temperatures, stationary systems are usually employed. Stationary systems are also less expensive and are easier to install and maintain. Concentrating collectors work best in climates with a high amount of direct solar radiation as in Cyprus. Evacuated Tube Collectors: These collectors are highly efficient, made of an absorber pipe enclosed within a larger glass tube. The absorber pipe may also be attached to a black copper fin that fills the tube (absorber plate). The space between the glass and the absorber is evacuated. Only radiant heat is transmitted from the absorber pipe to the outside of the collector and in this way the efficiency of the collector is substantially increased. The tubes are usually mounted into a manifold and the heated liquid is either used directly or circulates through a heat exchanger and gives off its heat to water. The coated surface of the absorber pipe has high absorption (> 92 %) of solar radiation and low emittance (< 6 %) for the infrared heat radiation. The absorber coating also has high resistance to long-term vapour condensation and high operating temperature. The performance of the system employing these three types of collectors is investigated in order to select the most suitable for the present application. The final selection is made by considering the financial viability of the system.

Hot water is stored in a TRNSYS Type 38 tank. The vertical cylinder construction is made of copper and is thermally insulated with polyurethane. Also the tank is protected by a galvanised outer shell.6 mm thick. The backup boiler (TRNSYS Type 6) is assumed to have a maximum heating rate of 18 kw and a set upper temperature of 93 C. A number of thermostats (TRNSYS Type 2) are also used in order to control: a. The flow to the solar panels, allowing the fluid to circulate only when the temperature of the fluid returning from the collectors to the storage tank is higher than that of the fluid delivered to the load; and b. The operation of the boiler, allowing the boiler to operate only when the temperature of the fluid delivered to the load is below an optimum value. In this case the boiler will keep the water temperature delivered to the absorption cooler always above 85 C. SYSTEM OPTIMISATION AND ENERGY FLOWS A number of simulations were carried out in order to optimise the various factors affecting the performance of the system. The parameters considered are as follows: The collector slope angle: The solar heat gain from the system for various collector slope angles was examined. The optimum angle in the Cyprus environment is: i. 25-3 for the flat plate collector ii. 3 for the compound parabolic collector, and iii. 3 for the evacuated tube solar collector This is due to the solar altitude angle, which for the latitude of Cyprus (35 ) can reach to 78 during noon in June. Also, because of the load characteristics, with the total cooling loads being about 6 times bigger than the heating loads, the optimum angle should be such that the collectors are absorbing greater heat during summer. Storage tank size: This factor also plays a role in the optimisation of the system. Simulations show that for all three types of collectors, a smaller tank size results in slightly less energy consumption by the boiler and slightly less energy collected by the solar collectors. The storage tank size is therefore decided only from the length of time intervals between firing the boiler. Between these intervals the tank should be able to supply the system with the needed water mass flow at the correct temperature in the summer. As it results, when the storage tank size is small the temperature in the cylinder cannot be kept above 85 C for long periods of time. Assuming that the boiler is firing every 12 minutes, which may be a reasonable interval, the optimum storage tank size needed for the system would be: i..8 m 3 for the flat plate collector ii..6 m 3 for the compound parabolic collector and iii..6 m 3 for the evacuated tube solar collector

Finally the effect of the collector area is evaluated against the boiler heat required. As it is expected the greater the collector area the less the boiler heat needed as indicated in Fig. 3 and the more the collected heat as indicated in Fig. 4. To determine the optimum collector area the cost of the solar system must be compared against the fuel saved due to its use. The cost of the various types of collectors is shown in Table 2. Table 2 Cost of solar collectors Collector type 21 Cost ( / m 2 ) Flat plate collector 186 Compound parabolic collector 35 Evacuated tube collector 425 27 Cost ( / m 2 ) 22 37 46 The economic method used, is the life cycle analysis [7]. For solar systems the following equation can be used: System annual cost = Extra mortgage payment + Maintenance cost Extra fuel savings Extra electricity savings Extra tax savings (1) Boiler heat requied (kwh) 3 28 26 24 22 2 18 16 14 12 flat plate evacuated tube 1 2 3 Collector area (m 2 ) Fig. 3 Effect of the collector area on the boiler heat required by the system Savings ( ) 5-5 - -15-2 -25-3 -35-4 -45 Flat plate Evacuated tube 5 1 15 2 25 3 35 Collector Area m 2 Fig. 4 Effect of the collector area on the collector heat gain

The scenario considered is that 3% of the initial cost of the system is paid at the beginning and the rest is paid in equal instalments in the next 1 years. The evaluation is based on the price of the diesel in 21, which was.29 per lt and the current price which is.73 per lt. The results of the economic analysis are presented in Fig. 5. As it is observed the only economically viable solution in 21 was to use the with a collector area of 15 m 2. The flat plate collector and the evacuated tube collector were not suitable for this application then due to the low price of fuel. Now, in 27, the situation has changed and and evacuated tube collectors can be used economically as indicated in Fig. 6. Savings ( ) - -2-3 -4-5 Flat plate Evacuated tube 5 1 15 2 25 3 35 Collector Area m 2 Fig. 5 Collector area against life cycle savings in for a fuel price of.29 per lt and a 2 year period. Savings ( ) 9 8 7 6 5 4 3 2 - -2 Flat p late Evacuated tube 5 1 15 2 25 3 35 Collector Area m 2 Fig. 6 Collector area against life cycle savings in for a fuel price of.73 per lt and a 2- year period.

SYSTEM LONG-TERM PERFORMANCE AND ECONOMIC ANALYSIS The specifications of the final system obtained from the optimisation study are shown in Table 3. The energy flows of the system are shown in Figs. 7 and 8 for the years 21 and 27 respectively. The cooling load of the building reaches a maximum monthly value of 42 kwh (in July), whereas the maximum monthly heating load occurs during January and is equal to 125 kwh. The heat required from the conventional boiler is also shown in the above figures. The maximum monthly load supplied by the solar system is 15 kwh and as can be seen from the difference of the curves for the cooling load and boiler heat, nearly all collector heat can be utilised for cooling or heating purposes. Table 3 The final optimum system in 21 and 27 Item Year 21 Year 27 Collector type Collector area Collector slope Storage tank size 15 m 2 3 6 lt 35 m 2 3 6 lt Energy flow (kwh) 6 5 4 3 2 collector heat boiler heat cooling load heating load JAN FEB MAR APR MAY JUN JUL AUG SEP OCT NOV DEC Months Fig. 7 System energy flows for 21 Energy flow (kwh) 5 45 4 35 3 25 2 15 5 collector heat boiler heat cooling load heating load JAN FEB MAR APR MAY JUN JUL AUG SEP OCT NOV DEC Months Fig. 8 System energy flows for 27

The annual cooling load of 17,6 kwh is covered with a total supply of 15,22 kwh of boiler heat, supplemented by 85 kwh of solar heat, offered by the solar system. The annual heating load of 353 kwh is covered with a total supply of 288 kwh of boiler heat and 15 kwh of solar heat. In 21, the total life cycle cost of a complete system comprising the collector and the absorption unit, for a 2-year period was 22,7 with an estimated price for the absorption unit together with its accessories equal to 8,14. In order to have the same total life cycle cost the price of the absorption unit together with its accessories should not be higher than 3,4. This could be possible only when absorption units were mass-produced. Today (in 27) the final optimum system, as obtained from the complete system simulations, consists of 35 m 2 compound parabolic collectors. The only limitation is the available area on the house roof, which cannot be completely blocked. The price of diesel in 27, which is.73 per lt indicates that it is economical to replace about 2, kwh with solar energy collected with 35 m 2 of and the life cycle savings of such a system are 7,7. The total life cycle cost of a complete system comprising the collector and the absorption unit, for a 2- year period is 38,2. CONCLUSIONS The final optimum system for 27 as obtained from the complete system simulations, consisted of 35 m 2 compound parabolic collector tilted at 3 from horizontal and 6 liters hot water storage tank. The typical insulated house considered, requires an annual cooling load at 25 C of 17,6 kwh with a peak load of 1.3 kw and an annual heating load at 21 C of 3,53 kwh with a peak load of 5.5 kw. This annual load could be met by using about 28,1 kwh of boiler heat. The price of diesel in 21, which was.29 per lt indicated that it was economical to replace about 1, kwh with solar energy collected with 15 m 2 of and the life cycle savings of such a system were 32. The life cycle savings for the present system and 27 fuel price (.73 per lt) give life cycle savings of 7,7. The increase of the cost of fuel compels consumers to move to other energy forms in order to meet their needs. So gradually renewable energy forms can economically replace part of the fuel energy as the fuel price increases. REFERENCES 1. Ghaddar, N., Shihab, M. and Bdeir, F. 1997. Modelling and Simulation of Solar Absorption System Performance in Beirut, Renewable Energy, 1 (4), pp. 539-558. 2. Erhard, A. and Hahne, E. 1997. Test and Simulation of a Solar-Powered Absorption Cooling Machine, Solar Energy, 59 (4-6), pp. 155-162. 3. Hammad, M. and Zurigat, Y. 1998. Performance of a Second Generation Solar Cooling Unit, Solar Energy, 62 (2), pp. 79-84. 4. Ardehali, M.M., Shahrestani, M. and Adams, C. 27. Energy simulation of solar assisted absorption system and examination of clearness index effects on auxiliary heating, Energy Conversion and Management, 48(3), pp. 864 87. 5. Florides, G., Tassou, S., Kalogirou, S. and Wrobel, L. 22. Measures Used to Lower Building Energy Consumption and their Cost Effectiveness, Applied Energy, Vol. 73, No. 3-4, pp. 299-328. 6. Florides, G., Kalogirou, S., Tassou, S. and Wrobel, L. 22. Modelling, Simulation and Warming Impact Assessment of a Domestic-Size Absorption Solar Cooling System, Applied Thermal Engineering, Vol. 22, No. 12, pp. 1313-1326. 7. Kalogirou, S. 1996. Economic Analysis of Solar Energy Systems Using Spreadsheets, Proceedings of the Forth World Renewable Energy Congress, Denver, Colorado, USA, Vol. 2, pp. 133-137.