 Open Access
 Total Downloads : 47
 Authors : Onovoh Johnpaul Nkemjika , Ujam Amechi Joseph
 Paper ID : IJERTV8IS050001
 Volume & Issue : Volume 08, Issue 05 (May 2019)
 Published (First Online): 04052019
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Evaluation of the Thermodynamic Properties of Air Cycle Refrigeration System
Ujam Amechi Joseph, Onovoh Johnpaul Nkemjika,
Department of Mechanical Engineering, Nnamdi Azikiwe University,
Awka, Nigeria.
Abstract: The effects of possible parameters affecting system performance air cycle refrigeration systems are examined through analysis of the thermodynamic properties. The study was carried out on varying number of air cycle refrigeration system used in air conditioning the cabin of an aircraft. Analysis was carried out for each of the system types to determine some parameters as cabin temperature (Tcb), Refrigerating effect (RE), compressor work (Wc) and coefficient of performance (COP) from the evaluated parameters. Results show that Bootstrap Air Cycle Refrigeration System has the best COP of 0.78 followed by the system working on reversed joule cycle (BellColeman cycle) with COP of 0.68, while the Reduced Ambient Air Cycle Refrigeration System produced the best refrigerating air temperature of 245.56K (27.440C) at entropy of 1.029127kJ/kgK that necessitated the cooling of the cabin from a temperature of 323K (500C) to 294K (210C). This is required comfort temperature for the cabin occupants. These have a whole lot of similarity with functioning data designed originally for the ACRSs meant for trains. Response surface methodology was done using Design Expert V11 (Trial Version) for optimization. The thermodynamic properties so determined in this study will be useful in designing Nigerian made aircraft and also serves as an important addition to the information storage bank of research institutes in the country.
Keywords— Refrigerating effects; compressor work; coefficient of performance; temperature; response surface methodology.

INTRODUCTION
In explaining refrigeration, it can be said to be a process of removing heat from an enclosed space or from a substance for the purpose of lowering the temperature. The flow of heat is in the direction of decreasing temperature, which is from high temperature areas to low temperature areas [1]. Reversed Carnot cycle can be use to demonstrate an ideal refrigerator which heat is taken in from a constant low temperature source at TL and is given out to a constant high temperature sink TH. Refrigeration in essence satisfies only one of the four requirements of the air conditioning, i.e. provision of lower temperature than the outdoor surroundings [2]. Air cycle refrigeration system (ACRS) is a gas cycle designed for the comfort of passengers in aircraft cabin as regards to this study. When the working fluid is a gas that is compressed and expanded but doesnt change phase, the refrigeration cycle is seen as a gas cycle [3]. Air is most often this working fluid. As there is no condensation and evaporation intended in a gas cycle, components corresponding to the condenser and evaporator in a vapour compression cycle are the hot and cold gas to gas heat exchangers in gas cycles. The air cycle is less efficient than the vapour compression cycle because the air cycle works on the reverse Brayton cycle instead of the reversed Rankine cycle [4].
For the purpose of this study, we will consider Aircycle refrigeration system as it relates to the cabin of an aircraft. Air does not undergo any phase change during the cycle; consequently, all the internal heat transfer processes are sensible heat transfer processes.
The significance of this research is to provide a dependable refrigerating property value of the ACRSs used which would be needed for the design of Nigerian customized ACRS considering the environmental factors. One main part of this study is resource evaluation and quantification as regards to refrigerating temperatures, coefficient of performance and thermal efficiency of ACRS.
Aircycle refrigeration system is widely used in some cold storage room, some industrial or commercial cooling towers, ships cabins and aircrafts cabins. For the purpose of this work, the analysis is basically for aircraft with a full passengers load and full crew within some 15 windows and 11% electrical load that flies at the speed of 10km/min at the altitude of 10km. The thermodynamic properties of the air at the altitude include 0.2650bar for pressure, 49.850C for temperature, 0.4134kg/m3 for density. This study excludes aircraft meant for wars, jets, and rockets etc. that travel at a very high speed.
1 1 Ram Air P6
Ram Air Ram Air P4 HE 2 P5 6
Po To
P1 P2 P3
Outside Air Primary Heat Exchanger 1 Secondary Cooling Compressor Compressor Turbine
Figure 1: Bootstrap Air Cycle Refrigeration System [5]

MATERIAL AND METHODS
Akanu Ibiam International Airport (AIIA) Enugu (ENU for IATA) Flight test facility was used for this study. The airlines used for the purpose of this study operate Akanu Ibiam International Airport (AIIA) Enugu (ENU for IATA). Dana airline (DAN for ICAO and AJ for IATA) and Discovery airline (DCV for ICAO and DO for IATA) and Arik airline (ARA for ICAO and W3 for IATA) were used. Six (6) air cycle refrigeration systems were all captured in different aircrafts which might be from the same airline as the case may be. The thermodynamic properties of the ACRS used in the aircraft were measured and sourced from the airline operatives. Except for Basic Air Cycle Refrigeration System with Evaporative Cooling of 2005 model seen in Discovery airline, every other ACRS in Dana and Arik airlines are of newer model (20092013). Reduced Ambient Air Cycle Refrigeration System, Regenerative Air Cycle Refrigeration System and Bootstrap Air Cycle Refrigeration System were studied from aircrafts under Dana airlines. Bell Coleman Cycle (System Working on Reversed Joule Cycle) and Basic Air Cycle Refrigeration System without Evaporative Cooling were installed in aircrafts under Arik airlines. However, the practical development of these thermodynamic properties was with the aid of measuring instruments. The measurements were performed continuously from the boarding (starting) time to the flight time. The measured data consisted of temperature, relative humidity and absolute pressure. The measured positions were mostly outside and ram air condition. The instruments used for measurement were the PRHT101 data recording instruments produced by Madgetech, an American company and the TRHP recording instruments produced by Qingsheng, a Chinese company. The instruments were placed in such a way that the sensor reads the average outside air condition at every altitude. These airlines have an air operator certificate issued by the civil aviation authority of Nigeria, (CAAN). The efficiencies of the compressors and turbines recorded were the ones generally used by the manufacturers of these aircrafts.

Evaluation of Temperature and Pressure at different Process
The relationship between the temperature and pressure at different processes are presented below are as follows
Take ca = Efficiency of auxiliary compressor, t = Efficiency of turbine/expander, To as Temperature of outside air, Po as Pressure of the outside air, P1 as Pressure of the air after ramming, P2 as Pressure of the air leaving the compressor, Pc as Cabin pressure (bar), Tc as Temperature of the cabin, T2 as Temperature of the leaving the compressor, T3 as Temperature of the air passed through the heat exchanger. All temperatures are in Kelvin and all pressures in bar.
Ramming action
Process 12 (Compression process)
T1 = To
1
(P1) (1)
Po
Process 23 (Condenser)
T2 = T1
1
(P2) (2)
P1
Temperature of air entering he expander at constant pressure is 500C Process 34 (Expansion process)
1
T3 = (P3) (3)
T4 P4
Where T4 and P4 are temperature and pressure entering the aircrafts cabin Process 341 (Secondary compression process)
T41
1
P4
3
3
= T ( ) (4)
P2
N/B: Discharge pressure (P4) of the air from the auxiliary compressor for Bootstrap air cycle refrigeration system =12bar
Process 45 (Heat exchanger)
ca
T 1 T3
= 4
T4T3
(5)
Process 56 (Expansion process)
T5 = T4(1 2) + T12 (6)
T 1
T 1
1
6 = (Pc)
T5 P4
(7)
6 5
6 5
T6 T5 t = T 1 T
(8)
Where T6 = Temperature of air entering the cabin
Refrigerating effect (RE) = Cp(Tc T6) (9)
Work of compression (Wc) = CP(T2 T1) (10)
Coefficient of performance (COP)
COP = Refrigerating effect
Work of compression
= RE
Wc
(11)

The Central Composite Design (CCD) of the thermal properties of Air cycle Refrigeration system
Version 11 CCD of response surface methodology (RSM) of Design Expert software was implemented in this study to design the experiment/optimization of the thermal properties of Air cycle Refrigeration system. The thermal properties were predicted using Central Composite Design (CCD). This was done to determine the best conditions for refrigeration in a cabin system. Also, this helps to examine the interactive effects of the three factors considered. In this study, the factors considered were air temperature (K), work of compression (kJ/kg) and refrigerating effect (kJ/kg) as independent variables while the coefficient of performance and entropy (kJ/kgK) were the dependent variables or responses. CCD involves varying the independent variables at five different levels (a, 1, 0, +1, +a). In this work, a set of 20 experiments were performed which consist of 8 factorial points, 6 centre points and 6 axial points. Moreover, this is because the replicates of factorial points and the replicates of axial points were two to increase the accuracy of the experiment and were performed in a randomized order. The range and levels of all individual variables used are given in Table 1. It is noteworthy to know that the optimal values of the test variables were first obtained in coded units and then converted to the uncoded units. The empirical equation is represented as shown below,
2 2 2
= 0 + + 2 + (12)
=1
=1
=
The regression analysis was performed to estimate the response function as a quadratic or second order polynomial: The importance of each term in the equation is to estimate the goodness of fit in each case. Response surfaces were drawn to determine the individual and interactive effects of the test variable on the coefficient of performance and entropy.
Table 1: Factors levels of independent variables in actual and coded form for COP and Entropy.
Factor Name
Units
Min. (a)
Max.(+ a)
Low(1)
High(+1)
Mean (0 Level)
Std. Dev.
Air Temperature(A)
K
232.15
248.96
1 235.56
+1 245.56
240.56
4.48
Refrigerating Effect (B)
kJ/kg
74.42
91.23
1 77.83
+1 87.83
82.83
4.48
Work of Compression (C)
kJ/kg
176.72
178.41
1 177.07
+1 178.07
177.57
0.4491


RESULTS AND DISCUSSION

Thermal Properties of ACRS
For each Air cycle refrigeration system, results were presented from the ambient condition of the air (ramming action). Temperature of air (refrigerant) that goes into the cabin (Tcb) was deduced, so also the Refrigerating effect (RE), Work of compression (WC) and the Coefficient of performance (COP) of the system. Plot was made to show the relationship between the ultimate Temperatures versus Entropies. Also the results obtained from the air cycle refrigeration systems showing the temperatures and pressures at every process were reported. In this research, it was also noticed that Bootstrap air cycle refrigeration system has the best COP of 0.78, followed by the system working on reversed joule cycle (Bell Coleman cycle with COP of 0.68. Reduced Ambient air refrigeration system has the least COP of 0.45.
6
5
4
3
2
1
6
5
4
3
2
1
0
50
100
150
Tcb(K)
200
250
300
0
50
100
150
Tcb(K)
200
250
300
Where
ACRSs
ACRSs
Figure 2: Bar chart showing Air Temperature output of ACRSs

System working on reversed Joule cycle (Bell Coleman cycle), 2. Basic air cycle refrigeration system without evaporative cooling, 3. Basic air cycle refrigeration system with evaporative cooling, 4. Bootstrap air cycle refrigeration system, 5. Regenerative air cycle refrigeration system and 6. Reduced ambient air cycle refrigeration system.
Figure 2 shows the difference in the temperatures obtained in all the ACRSs. This temperature is what balances the high temperature produced in the cabin of the aircraft by loads, passengers, electrical gadgets etc. The lower the temperature the more efficient and proficient outputs it yields. In the bar chart above, bar 4 which is the Bootstrap air cycle refrigeration system with temperature of 185.29K has the lowest temperature followed by the first bar which is BellColeman cycle. This is because the two ACRSs is widely used for aircrafts with smaller cabin. The ACRS with highest temperature is at sixth bar (Reduced ambient air cycle refrigeration system with largest temperature of 245.56K).
TS diagram for Bell Coleman cycle system is shown in figure 3 below. This was done at the critical conditions of 140.650C and 37.7436band showed an elaborate view that included enthalpy in kJ/kg and specific volume in m3/kg. At its ACRS temperature of 192.550C, the entropy extended from 80J/kgK to 2650J/kgK. It also indicated enthalpy state from 110kJ/kg and specific volume from 0.005 to 0.1m3/kg. The enthalpy reflects similar range in its pressureenthalpy plot. Reference [4] reported similar trend in the uses of Bell Coleman refrigeration system for air conditioning where the expanded air temperature is 289K and constant entropy from process 3 to 4.
Figure 3: TS diagram for Bell Coleman Cycle
Reference [6] reported in reversebootstrap air cycle refrigeration system driven by powered turbine similar trend being that this time it is for heat pump having inlet air temperature at 180K and outlet air temperature at 320K.
TS diagram for Bootstrap air cycle refrigeration system is shown in figure 5 below. This was done at the critical conditions of – 140.650C and 37.7436bar and showed an elaborate view that included enthalpy in kJ/kg and specific volume in m3/kg. At its ACRS temperature of 206.80C, the entropy extended from 300 to 3000J/kgK. It also indicated enthalpy state from 90kJ/kg and specific volume from 0.005 to 1.0m3/kg. The enthalpy reflects similar range in its pressureenthalpy plot.


CCD Analysis and Optimization of Coefficient of Performance and Entropy in ACRS
The COP depends mostly on the results if there is significant variation for combination of work of compression and refrigerating effect while entropy depends mostly on air temperature. The empirical relationship between the responses COP and entropy and the three variables i coded values obtained by using the statistical package DesignExpert 11 version for determining the levels of factors which gives optimum COP and entropy was given by the (13) and (14). The quadratic regression equations that fitted the data for all the new modeled thermal properties are illustrated also in (15) and (16).
Figure 4: TS diagram for Bootstrap Air Cycle Refrigeration System
Z = +0.342995+0.000373A+0.012391B0.004918C3.98720E20AB+9.91693E21AC
0.000037BC 8.65167E07A28.56618E07B2+0.000015C2 (13) S = +0.299235+0.021143A+0.000391B0.000018C0.00720AB+0.021693AC
3.214037 E21BC 0.0005167A28.56618E07B21.023015 E07C2 (14)
The accuracy and adequacy of the above proposed model was tested using the Design Expert sequential model sum of squares and the model test statistics. From the sequential test, it can be seen that the model Fvalues is 9415.8. The Model Fvalue of 9415.80 implies the model is significant and there is only a 0.01% chance that an Fvalue this large could occur due to noise.
Pvalues less than 0.0500 indicate model terms are significant because values greater than 0.1000 indicate the model terms are not significant.. In this case B, C, BC, CÂ² are significant model terms. Table 4.10 shows the regression coefficients and adjusted regression coefficient which are in close agreement with their respective predicted values. The Predicted RÂ² of 0.9990 is in reasonable agreement with the Adjusted RÂ² of 0.9998; i.e. the difference is less than 0.2. The same trend is noticed for entropy
when considered as response too. Adeq Precision measures the signal to noise ratio. A ratio greater than 4 is desirable. The ratio of 343.503 indicates an adequate signal; therefore the model can be used to navigate the design space.
The coefficient of variance (CV) is the ratio of the standard error of the estimate to the mean value of the observed response and is considered reproducible once it is not greater than 10%. In this work, the CV obtained for the entire model is 0.369% for COP and 0.443% for entropy.
Table 2: Regression coefficient of COP using quadratic model with respect to the independent variables
Source
Sum of Squares
df
Mean Square
Fvalue
pvalue
Block
0.0000
1
0.0000
Model
0.5041
9
0.0560
9415.80
< 0.0001
AAir Temperature
0.0000
1
0.0000
0.0000
1.0000
BRefrigerating Effect
0.4660
1
0.4660
78337.65
< 0.0001
CWork of Compression
0.0362
1
0.0362
6078.60
< 0.0001
AB
0.0000
1
0.0000
0.0000
1.0000
AC
0.0000
1
0.0000
0.0000
1.0000
BC
0.0018
1
0.0018
302.59
< 0.0001
AÂ²
8.889E06
1
8.889E06
1.49
0.2526
BÂ²
8.889E06
1
8.889E06
1.49
0.2526
CÂ²
0.0001
1
0.0001
18.30
0.0021
Residual
0.0001
9
5.949E06
Lack of Fit
0.0001
5
0.000025
4.22
0.02555
Cor Total
0.5042
19
The analysis of variance (ANOVA) results for the model terms are given in Tables 1. Analysis of variance (ANOVA) was applied for estimating the significance of the model at 5% significance level and a model is considered significant if the pvalue (significance probability value) is less than 0.05. From the pvalues presented in Table 2, it can be stated that B, C, BC, CÂ² are significant model term for COP. Due to these, the equation of the models is reduced thus:
Z = +0.342995 + 0.012391B – 0.004918C – 0.000037BC + 0.000015C2 (15) S = +0.299235 + 0.021143A – 0.00720AB + 0.021693AC – 0.0005167A2 (16)
In this study, the calculated data were also analyzed to check the correlation between the actual and predicted response. It can be seen from the figures that the data points on the plot were reasonably distributed near to the straight line, indicating a good relationship between the actual and predicted values of the COP and entropy separately. The underlying assumptions of the above analysis were appropriate. The result also suggests that the selected quadratic model was adequate in predicting the dependent variable or response variables for the data.
DesignExpertÂ® Software
Trial Version COP
Color points by value of COP:
0.35 0.97
1
0.9
0.8
Predicted
Predicted
0.7
0.6
0.5
0.4
0.3
Predicted vs. Actual
0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Actual
Figure 5: Predicted values versus the actual values for COP with respect to the independent variables
DesignExpertÂ® Software
Trial Version Entropy
Color points by value of Entropy:
0.611 1.122
1.2
1.1
Predicted
Predicted
1
0.9
0.8
0.7
0.6
Predicted vs. Actual
0.6 0.7 0.8 0.9 1 1.1 1.2
Actual
Figure 6: Predicted values versus the actual values for Entropy with respect to the independent variables
COP
Design points below predicted value 0.35 0.97
X1 = A: Air Temperature X2 = B: Refrigerating Effect
Actual Factor
C: Work of Compression = 164.71
1
0.9
0.8
0.7
COP
COP
0.6
0.5
0.4
0.3
138.4
131.67
124.94
118.21
111.48
104.75
98.02
245.56
236.95
228.34
219.73
211.12
202.51
B: Refrigerating Effect (kJ/kg9) 1.29
84.56
77.83
193.9
185.29
A: Air Temperature (K)
Figure 7: RSM 3D plot indicating interaction effects of factors (refrigerating effect and air temperature) on COP.
DesignExpertÂ® Software
Trial Version
Factor Coding: Actual

COP
COP
Design Points
0.35 0.97
X1 = A: Air Temperature X2 = B: Refrigerating Effect
Actual Factor
C: Work of Compression = 164.71
131.67
B: Refrigerating Effect (kJ/kg)
B: Refrigerating Effect (kJ/kg)
124.94
118.21
111.48
104.75
0.8
0.7
6
98.02 0.6
91.29
84.56
77.83
0.5
185.29 193.9 202.51 211.12 219.73 228.34 236.95 245.56
A: Air Temperature (K)
Figure 8: Response surface contour plot indicating interaction effects of factors (refrigerating effect and air temperature) on COP.
Entropy (kJ/kgK)
Design points below predicted value 0.611 1.122
X1 = A: Air Temperature X2 = B: Refrigerating Effect
Actual Factor
C: Work of Compression = 164.71
1.2
1.1
1
Entropy (kJ/kgK)
Entropy (kJ/kgK)
0.9
0.8
0.7
0.6
138.4
131.67
124.94
118.21
111.48
104.75
98.02
245.56
236.95
228.34
219.73
211.12
202.51
B: Refrigerating Effect (kJ/kg9) 1.29
84.56
77.83
193.9
185.29
A: Air Temperature (K)
Figure 9: RSM 3D plot indicating interaction effects of factors (refrigerating effect and air temperature) on Entropy.
DesignExpertÂ® Software
Trial Version
Factor Coding: Actual
Entropy (kJ/kgK) Design Points 0.611 1.122
138.4 Entropy (kJ/kgK)
131.67
124.94
X1 = A: Air Temperature X2 = B: Refrigerating Effect
Actual Factor
C: Work of Compression = 164.71
118.21
B: Refrigerating Effect (kJ/kg)
B: Refrigerating Effect (kJ/kg)
111.48
104.75
0.85 0.9 0.95 1 6
1.05
98.02
91.29
84.56
77.83
185.29 193.9 202.51 211.12 219.73 228.34 236.95 245.56
A: Air Temperature (K)
Figure 10: Response surface contour plot indicating interaction effects of factors (refrigerating effect and air temperature) on Entropy.
It can be seen that the error between the experimental and predicted is less than 1%, therefore it can be concluded that the generated model has sufficient accuracy to predict the response.
Table 3: Optimal regression coefficient of COP and Entropy with respect to the independent variables
Intercep t
Block [1]
A
B
C
AB
AC
BC
AÂ²
BÂ²
CÂ²
Entrop
y
1.02447
– 0.0180123
0.13
– 2.19168E
17
7.54624E
17
– 9.92733E
17
– 9.08243E
17
– 1.16238E
16
0.0651048
0.0115415
– 0.0115415
p values
< 0.0001
1.0000
1.0000
1.0000
1.0000
1.0000
< 0.0001
0.1518
0.1518
COP
0.66044
– 0.0012261
7
– 2.52439E
17
0.184721
– 0.051455
7
3.76436E
17
1.97326E
18
0.015
– 0.00078567
4
– 0.00078567
4
0.0027498
6
p values
1.0000
< 0.0001
< 0.0001
1.0000
1.0000
< 0.0001
0.2526
0.2526
0.0021
235.56 245.56 77.83 87.83
A:Air Temperature = 245.56 B:Refrigerating Effect = 87.83
177.07 178.07 1 1.04
C:Work of Compression = 178.07 Entropy = 1.03079
Desirability = 1.000
Solution 1 out of 100
Figure 11: The desirability range of the independent variables


CONCLUSION
The results of the study showed the level of effectiveness of the different ACRSs analyzed in this work considering Nigerian environment. The condition and thermodynamic properties of air used as a case study was gotten from catalogs of the airline operatives at AIIA based on year round mean conditions at 45o latitude and varies with the time of the year and the weather patterns.
However, to show the relevance of this work, it corresponds with the aim and objectives of the International Refrigeration and Air conditioning conference held at Purdue July 12 to 15, 2010 about aircrafts cabin. They noted that the cabin temperature would rise up to 500C (323K) within some 15 windows, 11% electrical load of lightening of air craft that flies at about 10km and 12km. This has a whole lot to do with scope and limitation of this research. From the foregoing, Bootstrap air refrigeration cycle is the efficient ACRS for this aircrafts cabin. This is because Bootstrap air refrigeration system produced a higher COP than other systems. Moreover, it has been long recognized that other air cycle refrigeration systems yield a close COP and Refrigerating effect. COP is one of the responses and refrigerating effect is one of the factors used in optimization of the thermal properties of ACRS using response surface methodology in this work. The COP, RE, Wc and Tcb of the different Air cycle refrigeration system has been successfully identified and optimized through response surface methodology in this research which left Bootstrap air cycle refrigeration system the better in terms of COP and RE where as Reduced Ambient ACRS is better in Tcb.
REFERENCES

P. Liping, Q. Yue, L. Dong, and L. Meng, Thermal comfort assessment in civil aircraft cabins. Chinese Journal of Aeronautics. Beihang University. Volume 27, Issue 2, Pages 210216, April 2014.

A. Ujam, Refrigeration and air conditioning (B. Eng, M. Eng, Eng Mnse, MNSE, MASME) 2009.

L. Shengjun, Z. Zhenying, and T. Lili, Thermodynamic analysis of the actual air cycle refrigeration system. Hebei United University, Tangshan 063009, China. Earthquake engineering research center of Hebei Province, Tangshan 063009, China. Tangshan College, Tangshan 063000, China, 2011.

U. Shet, T. Sundararajan, and J. Mallikarjuna, Note on Refrigeration cycle. Indian institute of Technology Madras, 2012.

P. Desai, Refrigeration and Air conditioning. Watana press 01435 Khana publishers. 2B, Nath Market, Nai Sarak, New Delhi11007 (India), pp. 1 8, 2004.

Z. Deng, (2009) Research on Reversebootstrap Air Cycle Refrigeration System Driven by Powered Turbine. Master's thesis, Nanjing University of Aeronautics and Astronautics, 2009.