Modelling
Modelling
evaporator were carried out by Horuz et al. [5]. Water-cooled system cycle and the thermal duty of each component. For the
absorber and condenser were introduced in addition to the solution to converge, it requires careful selection of input par-
original air-cooled ones so as to examine the system perform- ameters, such as the UA or effectiveness values of the heat
ance under the conditions of a wider range of high-side press- transfer components. This makes the model difficult to use by
ures. The experimental investigations were carried out under inexperienced users. To obtain a reasonable prediction of
six different operating conditions including variable chilled system performance by the model, the correlated value of UA
water inlet temperature and flow rate, changeable chilled water or effectiveness for each component from test results is still
level in the evaporator drum, varied generator heat input, and needed.
different cooling water flow rate and inlet temperature to the Apart from the above system models formed by the inte-
water-cooled condenser. It was found that the cooling capacity, gration of simplified component models, the detailed com-
component heat transfer rates, and cooling coefficient of per- ponent models are still significant to describe operationally
formance (COP) were all increased with higher chilled water and assist the optimal designs for the components. Chua et al.
inlet temperature. In the mean time, the above three perform- [11] developed a detailed absorption generator model in which
ance parameters were all increased with elevated energy input the multi-stream non-isothermal heat and mass transfers were
to generator. In addition, the cooling capacity and COP were recognized and considered. Goel et al. [12] developed an
both augmented with increased condenser cooling water flow analytical model to simulate the combined heat and mass
rate but the COP was decreased with higher condenser cooling transfer process in a counter-current ammonia– water-based
water temperature. The solution high and low pressures which absorber. The model requires detailed description of the geo-
were quite important to the system performance were however metric characteristics of the absorber. It uses empirical corre-
not presented. lations to predict the heat and mass transfer coefficients and
The models of different types of absorption refrigeration accounts for both liquid and vapour phase mass transfer resist-
systems were reported in series by Grossman et al. [6-10] and ances. Similar investigations on absorbers were carried out by
were finalized into a modular simulation tool—ABSIM for Kaynakli et al. [13] and Kwon et al. [14]. All these models
advanced absorption systems. The codes underneath the require detailed description of the components’ geometric
system models are based on unit subroutines containing the data, which inversely limits their applications to different types
governing equations for the system’s components. Each com- of system.
ponent is assumed as one control volume and the applied gov- Any simulation model requires a method for calculating
erning equations for the control volume are similar to those the thermophysical properties of the fluid pair in the
used by Lazzarin et al. [3,4]. Based on user-supplied cycle system. A simple, practical and fairly accurate method to do
diagrams, working fluid and given operating conditions, this for the ammonia– water fluid pair was proposed by
the model calculates the properties at each state point in the Conde [15].
This paper presents a model of a low-temperature air- ‘1’. The weak solution is heated and distilled in the generator
cooled ammonia – water absorption chiller system. The system and as a result, a strong solution is formed at the bottom of
model is an integration of the main component sub-models. the generator and high concentration ammonia refrigerant at
To avoid the need to have detailed data on the geometric the top. The strong solution then exits the generator at ‘2’ and
characteristics of the components and ensure speed of solution, flows through an expansion device before entering the absorber
each component model was based on the lumped parameter where it is sprayed over the SHX. The refrigerant flows from
methodology. To facilitate the modelling processes, the genera- the top of the generator to the rectifier where more water in
tor and rectifier are combined into one control volume while the ammonia refrigerant condenses and returns back to the
the absorber is divided into two control volumes. In the mean generator. The nearly pure refrigerant then exits the rectifier at
time, extensive experiments were carried out including the ‘8’ and flows to the condenser for the repetition of the cycle.
variation of high and low solution pressures with the tempera- The model of each component was developed based on the
tures of air inlet and brine outlet. The test results are used to energy, mass and species conservation equations. For simpli-
validate the model and correct the heat transfer parameters city, each component was assumed as a single or two control
needed in the component models. The following sections volumes. The models of the more complicated components in
describe the modelling methodology adopted, validation the system are described as below.
against experimental measurements and model applications.
This is a first stage in the development of an integrated model
of a microturbine-based tri-generation system for low-
temperature food refrigeration applications. 2.1 Rectifier
The rectifier was assumed to be a single control volume. The
assumptions made in developing the model are:
2 MODELLING SCHEME WITH THE (1) Pure saturated refrigerant vapour at rectifier exit ( point 8
LUMPED METHOD in Figure 1).
(2) Saturated vapour solution at top of generator ( point 100 ).
A schematic diagram of the refrigeration system considered in (3) Saturated liquid solution from rectifier to generator ( point
this investigation is shown in Figure 1. Referring to the 10 ).
diagram, the system consists of eight main components: gen- (4) No heat loss in rectifier.
erator, rectifier, air-cooled condenser, sub-cooling heat exchan-
ger, liquid-cooling evaporator, absorber (solution-cooled), Conservation of mass for ammonia:
air-cooled heat exchanger and solution pump. After flowing ð1 þ RÞj100 ¼ j8 þ Rj10 ð1Þ
through the generator and rectifier, the nearly pure refrigerant
at ‘8’ enters the air-cooled condenser where it is condensed to Conservation of energy:
liquid at ‘9’. The refrigerant then passes through the connec- ð1 þ RÞh100 ¼ h8 þ Rh10 þ qr ð2Þ
tion pipes and arrives at ‘10’. From there the refrigerant flows
where
through a tube-in-tube heat exchanger where it exchanges heat
with the refrigerant vapour leaving the evaporator and is R¼m
_ 10 =m
_8
further subcooled at ‘11’. The refrigerant then flows through a
thermodynamic expansion valve at ‘12’ before entering the
evaporator coil where it evaporates while absorbing heat from
a brine solution. 2.2 Generator and rectifier
The refrigerant vapour exiting the evaporator at ‘13’ then The generator and rectifier are combined into one control
flows through the tube-in-tube heat exchanger where it is volume to facilitate the modelling. Assumptions for the gen-
superheated to ‘14’ before entering the absorber. The ‘strong’ erator are:
solution from the generator enters the absorber at ‘3’ where it
is first cooled by a solution heat exchanger (SHX) before (1) Saturated liquid for weak solution outlet of generator at
absorbing the refrigerant flowing from the evaporator. It ‘2’.
should be noted that in this paper ‘strong’ solution refers to a (2) No heat loss in generator.
solution that is rich in absorbent, whereas ‘weak’ solution
refers to a solution that is weak in absorbent and rich in Conservation of mass for ammonia:
refrigerant. The ‘weak’ solution formed in the absorber then
flows through the air-cooled heat exchanger where the heat f j1 ¼ ð f 1Þj2 þ j8 ð3Þ
generated by the absorption process is rejected to the ambient. Conservation of energy:
The solution then exits the heat exchanger at ‘5’ and is
pumped through the rectifier and absorber to the generator at qg ¼ h8 h2 þ qr þ f ðh2 h1 Þ ð4Þ
_ 3a j3a þ m
m _ 14 ¼ m
_ 4 j4 ð9Þ _ 12 ¼ m
m _8
Heat transfer:
Q_ e ¼ UA LMTD ð14Þ
where
DTe;max DTe;min
LMTD ¼
ln DTe;max =DTe;min
DTa;max ¼ maxðTcf in T13 ; Tcf out T12 Þ
DTa;min ¼ minðTcf in T13 ; Tcf out T12 Þ
and simulation over a range of chilled fluid outlet temperatures experiments, the ambient temperature which is quite impor-
and ambient temperatures. The scatters that can be observed in tant to the system performance varied between 15.58C and
the experimental results are mainly because of the variation in 17.58C. To investigate the performance of the system over a
ambient temperature during the tests and the slow response of wider range of ambient temperatures likely to be encountered
the system to variations in both the ambient temperature and in practice, simulations were carried out with ambient temp-
chilled fluid temperature. Despite the scatters, however, there eratures between 108C and 308C and chilled fluid flow temp-
are clear trends that show both cooling capacity and coling eratures between 2308C and 58C.
COP of the system to increase with increasing chilled fluid Figure 11 shows the variation of ammonia concentration in
outlet tempertaure at specific ambient temperature. This is the weak and strong solution with chilled fluid flow tempera-
mainly because of the higher evaporating temperature, which ture and ambient temperature. It can be seen that both the
leads to a higher refrigerant density and hence higher refriger- strong and weak solution concentrations increase with chilled
ant flow rate in the system. fluid temperature.
As the chilled fluid tmperature increases, the presures in
both the low and high pressure sides of the system increase as
4 MODEL APPLICATION shown in Figure 12. In addition, it shows that an increase in
the ambient temperature leads to an increase in the pressure in
The validated model can be used to explore the performance the condenser and the generator but has little effect on the
of the absorption unit at various operating conditions. In the pressure in the evaporator and absorber.
Figure 6. Tested ammonia water solution temperatures at generator and Figure 8. Comparison between simulated and experimental pressures in the
absorber inlets and outlets. low pressure side of the system.
Figure 7. Comparison between experimental and simulation results for the Figure 9. Comparison between experimental and simulation results for
high pressure side of the system. cooling capacity.
Figure 10. Comparison between test and simulation results for coefficient of Figure 13. Influence of chilled fluid flow temperature and ambient
performance. temperature on refrigerant and strong solution flow rates.
R ratio of return flow [2] Herold KE, Radermacher R, Klein S. Absorption Chillers and Heat Pumps.
T, t temperature (8C, K) CRC Press, 1996.
Tamb ambient air temperature (8C) [3] Lazzarin RM, Gasparella A, Longo GA. Ammonia-water absorption
machines for refrigeration: theoretical and real performances.
DT temperature difference (K)
Int. J. Refrigeration 1996;19:239– 46.
UA overall conductance (W/K)
[4] Lazzarin RM, Gasparella A, Romagnoni P. Experimental report on the
j concentration (% by mass) reliability of ammonia-water absorption chillers. Int. J. Refrigeration
1996;19:247– 56.
SUBSCRIPTS [5] Horuz I, Callander TMS. Experimental investigation of a vapor absorption
refrigeration system. Int. J. Refrigeration 2004;27:10– 16.
4l saturated liquid at point ‘4’ in Figure 1 [6] Grossman G, Perez-Blanco H. Conceptual design and performance analy-
a absorber, air sis of absorption heat pumps for waste heat utilization. Int. J. Refrigeration
cf chilled fluid 1982;5:361– 70.
[7] Grossman G, Gommed K, Gadoth D. A computer model for simulation of
e evaporator
absorption systems in flexible and modular form of absorption systems in
g generator
flexible and modular form (ORNL/Sub/90 ¼ -89673). Oak Ridge National
in inlet Laboratory, Oak Ridge, TN, 1991.
out outlet [8] Grossman G. Modular and flexible simulation of advanced absorption
r rectifier systems. In: Proceedings of the International Absorption Heat Pump
s strong solution (rich in absorbent) Conference, Vol. 31, AES, New Orleans, LA, 1994, 345– 52.
p pump [9] Grossman G, Wilk M. Advanced modular simulation of absorption
w weak solution (weak in absorbent) systems. Int. J. Refrigeration 1994;17:231– 44.
[10] Grossman G, Zaltash A. ABSIM-modular simulation of advanced absorp-
tion systems. Int. J. Refrigeration 2001;24:531 – 43.
[11] Chua HT, Toh HK, Malek A, Ng KC, Srinivasan K. A general thermodyn-
ACKNOWLEDGEMENTS amic framework for understanding the behaviour of absorption chillers.
Int. J. Refrigeration 2000;23:491– 507.
The authors wish to acknowledge the Food Technology Unit of [12] Goel N, Goswami DY. Analysis of a counter-current vapour flow absorber.
DEFRA for their financial support for this project and the con- Int. J. Heat Mass Transfer 2005;48:1283– 92.
tribution of the industrial collaborators, Bond Retail Services [13] Kaynakli O, Horuz I. Comparison of parallel and counter flow coil absor-
Ltd, Apex Air Conditioning, Doug Marriott Associates and ber performance. Int. J. Heat Mass Transfer 2006;33:211 –23.
Bowman Power. [14] Kwon K, Jeong S. Effect of vapor flow on the falling-film heat and mass
transfer of the ammonia/water absorber. Int. J. Refrigeration 2004;27:
955 – 64.
REFERENCES [15] Conde M. Thermophysical properties of fNH3 þ H2Og solution for the
industrial design of absorption refrigeration equipment. M.Conde
Engineering, 2006.
[1] Faramarzi PT, Walker DH. Investigation of secondary loop supermarket
[16] TRNSYS. A Transient System Simulation program. TRNSYS 16, 2005.
refrigeration systems. Consultant Report 500-04-013, March 2004.