# Energy Conversion and Transmission Characteristics Analysis of Ice Storage Air Conditioning System Driven by Distributed Photovoltaic Energy System

Authors: Yongfeng Xu; Ming Li; Reda Hassanien Emam Hassanien
Journal: International Journal of Photoenergy
 (2016)

---

## Abstract

In order to reduce the investment and operation cost of distributed PV energy system, ice storage technology was introduced to substitute batteries for solar energy storage. Firstly, the ice storage air conditioning system (ISACS) driven by distributed photovoltaic energy system (DPES) was proposed and the feasibility studies have been investigated in this paper. And then, the theoretical model has been established and experimental work has been done to analyze the energy coupling and transferring characteristics in light-electricity-cold conversion process. In addition, the structure optimization analysis was investigated. Results revealed that energy losses were high in ice making process of ice slide maker with only 17.38% energy utilization efficiency and the energy efficiency and exergy efficiency of ISACS driven by DPES were 5.44% and 67.30%, respectively. So the immersed evaporator and cointegrated exchanger were adopted for higher energy utilization efficiency and better financial rewards in structure optimization. The COP and exergy efficiency of ice maker can be increased to 1.48 and 81.24%, respectively, after optimization and the energy utilization efficiency of ISACS driven by DPES could be improved 2.88 times. Moreover, ISACS has the out-of-the-box function of ordinary air conditioning system. In conclusion, ISACS driven by DPES will have good application prospects in tropical regions without power grid.


---

## Body

## 1. Introduction

With the dramatic climate changes, the cooling demand has been increased and led to a rapid growth of energy consumption, which causes traditional fossil fuel energy shortage and great damage to climate and environment with the emissions of CO2 and harmful particles by extensive use of traditional fossil energy. Furthermore, a large number of the uses of the electric air conditioning can increase the tense situation between power supply of grid and demands of people. Therefore, refrigeration driven by solar energy becomes one of the promising approaches to reduce or partially replace the conventional refrigeration systems. Solar thermal refrigeration and solar photovoltaic refrigeration are two main working modes [1, 2]. Several studies have been conducted about different solar refrigeration options in the last years for the improvement and development of solar thermal refrigeration system such as operation efficiency and operating stability [3–10]. Compared to the solar thermal refrigeration, solar photovoltaic (PV) refrigeration has more advantages on refrigerating effect, stable operation, and energy utilization rate. Along with the decrease of the investment and operation cost and the gradual increase of the conversion efficiency of PV module, the photovoltaic refrigeration will be developed rapidly.There are two refrigeration models: thermoelectric refrigeration and vapor compression refrigeration can be driven by PV. As early as 2003, Dai et al. [11] researched thermoelectric refrigerator driven by solar cells and the results showed that the system has 0.3 COP and the refrigerator temperature could maintain 5°C~10°C. Due to the limitation of the working principle, thermoelectric refrigerator is only suitable for cold storage and cannot be used for freezing. So vapor compression refrigeration driven by PV caught researchers’ attention.Firstly, Aktacir [12] designed a multifunctional PV refrigerator and found that when indoor and outdoor average temperatures were 26.3°C and 24.9°C, the minimum temperature of the refrigerator reached −10.6°C. The research results showed that the refrigeration efficiency of the system needed to be improved. In order to improve PV refrigeration performance, Mba et al. [13] used MATLAB software to simulate PV refrigeration system operating process and analysis system operating characteristics in different conditions. The structure optimization and the great operation scheme of the system were mentioned in the research results. Furthermore, Tina and Grasso [14] designed a software program for monitoring and managing stand-alone PV refrigerator system in a remote area, which had a real-time monitor running status and automatically recorded data feeding back to the terminal equipment. The research work can provide references for the system operation stability and optimization. It was also reported that the required photovoltaic panel area of a solar electric-vapor compression refrigeration system increases as the evaporating temperature decreases and the coefficient of the performance variation of the cooling system decreases with the decrease of evaporating temperature which can provide an important reference on the structure optimization of PV refrigerator [15]. So in the investigation on the structure optimization of PV refrigeration, Kaplanis and Papanastasiou [16] improved the performance of a traditional refrigerator driven by PV through the structure rebuild and the performance of the optimized refrigerator could be improved greatly. In addition, in order to improve the operation stability of the PV refrigerator and reduce the cost of inverter, the AC compressor was replaced by DC compressor in the research work led by Ekren et al. [17]. The result showed that PV module conversion efficiency has a greater impact on the system exergy efficiency. Finally, for the comprehensive performance of PV refrigerator such as the energy conversion, management, and operation performance, powered on three conditions such as photovoltaic components, battery, and outage system, the results showed that the system COP gradually decreases from morning till night [18].According to comprehensive analysis, PV refrigeration system research is currently mainly concentrated on ice maker driven by PV. Batteries are essential component to store energy and to solve the intermittency of solar energy in PV refrigeration system. However, the use of batteries can increase the investing and running costs. So the PV refrigeration without batteries or batteries replacement technology was investigated. Axaopoulos and Theodoridis [19] designed a PV Ice maker without battery and studied its performance when the compressor operating efficiency was 9.2%; they found that this prototype has a good ice making capability and reliable operation as well as a great improvement in the startup characteristics of the compressors, which remain working even during days with low solar irradiation of 150 W/m2. But the refrigeration efficiency needed to be improved. The American SOLUS Refrigeration Company has developed a photovoltaic DC refrigerator substituting battery with water-propylene/ethylene glycol phase change material to store cold and reduce the system investing and running cost. Results revealed that the temperature inside the refrigerator was remaining stable at around 1.4°C when environment temperature was 32°C [20]. But the use of the water-propylene/ethylene glycol phase change material in the system increased the cost. Thus, the ice storage technology, stable performance, low cost, and large phase change latent heat have attracted researchers’ attention. Ice storage technology has a great role in saving building energy, transferring peak power to off-peak, improving grid load rate, and other aspects.The performances of ice storage were analyzed. Pu et al. [21] analyzed the effect of incorporating the ice thermal storage air conditioning system in power supply based on the cumulative exergy analysis method. They found that the total cumulative exergy consumption increased as the ice thermal storage (ITS) system applied. However, the average cumulative exergy variation decreased slightly as the cooling load of the ITS system increased. And then, the ice thermal energy storage (ITES) air conditioning system incorporating a phase change material (PCM) was analyzed from energy, exergy, economic, and environmental aspects [22]. Results showed that the electricity consumption of hybrid system was 6.7% and 17.1% lower than that of the simple ITES (without PCM) system and the conventional system, respectively. The economic performance of ITES was analyzed by Sanaye and Shirazi [23]. The results revealed that the amounts of electricity consumption and CO2 emission of ITS system were lower 9% and 9.8%, respectively, in comparison with those of a conventional system. Furthermore, the payback period of ITS system for extra capital cost was 3.43 years. Because the ITS had superior economic performance, the ITS was widely used in the large building cooling system throughout the world. Han et al. [24] conducted an analysis study to compare between the potential impacts of ice storage systems on load-shifting under a new credit-based incentive scheme and the existing incentive arrangement in Jiangsu, China. Results indicated that adopting different schemes of ice storage can make distinct impacts on load-shifting. Wang and Dennis [25] conducted an investigation for influencing factors on the energy saving performance of battery storage and phase change cold storage in a PV cooling system with TRNSYS in three distinct climates (Madrid, Shanghai, and Brisbane). Results showed that the CO2 clathrate hydrate was the best performing cold store and this approach provided credible savings compared to a system without energy storage.According to the above analysis, nowadays, the research and utilization of PV refrigeration and ice storage are relatively independent. In order to integrate the advantages of the two technologies together, the ice storage air conditioning system (ISACS) driven by distributed photovoltaic energy system (DPES) was established based on our previous research results [26]. This is the first time that the ice storage technology was employed in the PV system to substitute or partly replace batteries for the solar energy storage. As you know, it is clear that the tropical regions, such as Xishuangbanna in China and Bangkok in Thailand, have almost 300 days for cooling demand; thereby ice storage has a good application prospect in those regions. So our research work has certain significance. In our system, a few batteries were also used in order to provide the stability for system operation. At first, a theoretical model has been established and experimental work has been done to analyze the energy coupling and transferring characteristics in light-electricity-cold conversion process of ISACS driven by DPES. And then, the system optimization analysis was also investigated in order to achieve higher energy utilization efficiency and better financial rewards. In order to solve the problem of cold supplying process must lagging behind ice making process in traditional ice storage air conditioning system and improve the efficiency of traditional ice storage air conditioning system, system structure optimization analyses were investigated. So the immersed evaporator and cointegrated exchanger were suggested. Therefore, the present study can achieve the purpose of ISACS supplying cold efficiently with out-of-box functionality as the ordinary vapor compression air conditioning system. The research results can also provide some references for commercial application of ISACS driven by DPES without batteries in tropical regions without power grid.


## 2. Materials and Methods

### 2.1. ISACS Driven by DPES

#### 2.1.1. Configuration of ISACS Driven by DPES

Ice storage air conditioning system (ISACS) driven by distributed photovoltaic energy system (DPES) was mainly configured by DPES, ice maker, storage system, and air conditioning system. The working diagram of ISACS driven by DPES is shown in Figure1.Figure 1
Work diagram of ISACS driven by DPES with batteries.PV modules convert solar energy into electric energy which can be regulated by controller with maximum power point tracking to drive ice maker, ice storage system, and air conditioning system. In daytime, DPES receives solar energy and turns it into direct-current (DC) electric power which can be converted to alternating current (AC) electric power by inverter to drive AC compressor, water pump, ethylene glycol pump, and fan coil. To maintain the stability of electric energy supply, batteries were adopted and connected with controller to maintain the energy conversion and supply in the most optimized way. Ice maker and storage system were made up with AC compressor, condenser, expansion valve, disc evaporator, and ice storage tank. Circulating water can be frozen in a disc evaporator and the ice can be drooped into the ice storage tank when the hot refrigerant flowed into the evaporator which was controlled by the solenoid valve. Thereby, the ice maker worked as vapor compression refrigeration. In AC compressor, cryogenic R134a was compressed to high temperature and high pressure gas to be filtered in economizer and to release heat in condenser. Refrigerant was condensed to mild temperature and high pressure gas. When the gas inflows into throttle valve, it can be throttled to low temperature and low pressure liquid and then feeds into plate evaporator. And then, refrigerant flows into the other gas-liquid separator to be sucked into the compressor. Thereby, the refrigeration cycle will be completed. Air conditioning system was mainly made up of coil heat exchanger which was fixed in ice storage tank, ethylene glycol pump, solenoid valve, proportional control valve, and fan coil. Ethylene glycol is an adopted cold exchanging medium.According to the working principle diagram, 0.2 kW ISACS driven by DPES was established as shown in Figure2. The main component parameters of ISACS driven by DPES are shown in Table 1.Table 1
Main component parameters of ISACS driven by DPES.
Component Model Parameters DPES PV module JN-245 P m: 245 W, Vm: 34.5 V, Im: 7.10 A, VOC: 43.5 V, ISC: 8.18 A, module: length∗wide: 1640 mm∗990 mm, cells in series: length∗wide∗numbers: 155 mm∗155 mm∗60 Controller PL60 12–48 V 60 A charge, 30 A load Inverter Solar 48 V P: 3 kW, DC input voltage: 48 V,output voltage: 220 V, output frequency: 50 Hz Batteries SP12-65 Battery capacity: 12 V 65 Ah, four batteries in series ISACS Refrigerant R134a Molecular formula: CH2FCF3, boiling point: −26.1°C, critical temperature: 101.1°C Ice maker IM50 Ice production: 2.12 kg/h,P: 380 W Ice storage tank / Capacity: 20 cm∗20 cm∗20 cm Cold exchanging medium Ethylene glycol Melting point: −12.6°C, viscosity: 25.66 mPa·s Pump RS15-6 Power: 46–93 W, life: 6 m, maximum flow rate: 3.4 m3/h Fan coil / Fan type: YS 56-2, power: 180 W, voltage: 380 V, current: 0.53 A, speed: 2800 r/min, number of fins: 95, size: 23 cm∗8 cm∗20 cm, coil numbers: 26, coil inner diameter: 6 mmFigure 2
Pictures of 0.2 kW ISACS driven by DPES.DPES were made up of two 245Wp polycrystalline silicon PV in series and four batteries in series were used to store electric power. Refrigerant (R134a) temperatures and pressure were measured by T-type thermocouples and pressure transducers, respectively. Voltages and currents of PV modules were measured by a digital multimeter. The wind speed was measured by the wind speed transducer. Solar irradiation was measured by pyranometer. Compressor input power was measured by a wattmeter. Electromagnetic flow meter was used to measure refrigerant flow and cold exchanging medium flow. The parameters of all instruments are shown in Table 2.Table 2
Measuring instrument parameters and uncertainty.
Instrument Model Range Accuracy Application scope Maximum relative error Maximum absolute error Uncertainty(B class) Pyranometer Kipp & Zonen CMP-6 0–2000 (W/m2) ±5% 0–1000 (W/m2) ±10% ±100 W/m2 57.7348 W/m2 Thermocouples T −200 to 350 (°C) ±0.4% 0–150 (°C) ±0.93% ±1.4°C 0.8083°C Wind speed transducer EC-9S 0–70 (m/s) ±0.4% 0–10 (m/s) ±2.8% ±0.28 m/s 0.1617 m/s Electromagnetic flow meter KROHNE OPTIFLUX 5000 DN 5; 0–12 (m/s) ±0.15% 0–5 (m/s) ±0.36% ±0.018 m/s 0.0104 m/s Pressure transducer YOKOGAWA EJA430E 0.14–16 (MPa) ±0.055% 0–2 (MPa) ±0.44% ±0.0088 MPa 0.0051 MPa Electronic balance AHW-3 0–3 kg ±0.05 g 0–3 kg ±0.05% ±0.0015 kg 0.0009 MPa Wattmeter DELIXI DDS607 0–10000 kW·h ±0.01 kW·h 0–100 kW·h ±1% ±1 kW·h 0.5774 kW·h Digital multimeter FLUKE F-179 Voltage: 0–1000 V ±0.9% 0–380 V ±2.37% ±9.006 V 5.1996 V Current: 0–10 A ±1% 0–10 A ±1% ±0.1 A 0.0577 A


### 2.2. Theoretical Models on Energy Conversion and Transfer Characteristics of ISACS Driven by DPES

Theoretical model was established to analyze the energy conversion and transfer characteristics of ISACS driven by DPES as follows.


#### 2.2.1. Energy and Exergy Models of PV Modules

Energy balance equation of PV modules is expressed as(1)mpCpmoduleTp=QPV,in-QPV,loss-QPV,elect.Solar energy absorbed by PV modules can be estimated as(2)QPV,in=ατGSp.Electric power is expressed by(3)QPV,elect=VoutIout.Exergy and exergy losses of PV modules are given in [26]: (4)ΔEPV=QPV,loss1-TaTP+IscVoc-ImVm,EPV,in=GSp1-43TaTsun+13TaTsun4,EPV,in=EPV,out+ΔEPV.


#### 2.2.2. Energy and Exergy Models of Controller

Energy balance equation of controller can be estimated as follows:(5)QPV,elect=QC,out+QC,loss.QC,loss is the consumed energy by controller per unit time (W); usually QC,loss is 4% of QPV,elect.Exergy and exergy losses of controller are expressed as in [27]:(6)ΔEC=QC,loss1-TaTC,EC,in=EC,out+ΔEC.


#### 2.2.3. Energy and Exergy Models of Batteries

Energy balance equations of batteries are given by [28](7)QB,in=ξQB,out+QB,loss-ξQB,storage.ξ is a symbolic coefficient. When batteries supply power to load, ξ is assigned 1. On the contrary, ξ is assigned −1 when batteries are charged:(8)QB,loss=IB2r1+r2SOC+1r3-r4SOC.Here,(9)SOC=SOC0+∫0tIBdtBC.VF is the full charge rest voltage (V), SOC0 is the initial SOC value (1), and BC is batteries capacity (Ah).Exergy balance and exergy loss equations of batteries are written as [27](10)ΔEB=QB,loss1-TaTB,EB,in=EB,out+ΔEB.


#### 2.2.4. Energy and Exergy Models of Inverter

Energy balance empirical equations of inverter are given in [28]:(11)QI,loss=CPCR+CR-1CRQC,out+QB,out,QC,out+QB,out=QI,out+QI,loss.Exergy balance and exergy loss equations of inverter are expressed in [28]: (12)ΔEI=QI,loss1-TaTI,EC,out+EB,out=EI,out+ΔEI.


#### 2.2.5. Energy and Exergy Models of ISACS

ISACS are driven by stable electric energy outputted by DPES. The thermodynamic cycle refrigeration process of ice maker, refrigerant, R134a thermodynamic properties, andP-h changing curve of refrigerant R134a in refrigeration cycle are shown in Figure 3. It can be clearly seen in Figure 3 that area 1 is the saturated liquid, area 2 is the gas-liquid coexistence, and area 3 is the saturated vapor. Consequently, the refrigeration cycle was 1-2-2s-3-4-5-1 across the saturated vapor region and gas-liquid coexistence region. The process from 1-2 is showing an actual compressor working process, when the dry steam in state 1 is compressed to saturated vapor in state 2, which will be isentropic compressed to 2s under ideal condition. However, the process 2-2s-3-4 is the refrigerant condensing and releasing heat process in condenser. The process 2-2s is the natural cooling process and 2s-3-4 is the isobaric cooling process. The cryogen will be high pressure dry steam in state 3 and high pressure saturated liquid in state 4. The process 4-5 is the decompression cooling by throttle valve. Therefore, by flowing through throttle valve, the refrigerant temperature and pressure will be low and then it can flow into evaporator to absorb heat and refrigerate, which is the process 5-1.Figure 3
R134a thermodynamic curve.Energy and exergy analyses are shown as follows.(a)
The compressor energy balance equation is given by(13)m˙rfh1+W=m˙rfh2+QCP,loss.Compressor exergy model is given as(14)ECP,in=ΔECP+ECP,out,ΔECP=m˙rfh2-h1-Tas2-s1.(b)
The condenser energy balance equation is shown as(15)m˙rfh3=m˙rfh4+QE+QCO,loss,QE=m˙air,icemkCp,airTair,out-Tair,in.Condenser exergy model is expressed as(16)ΔECO=m˙rfh3-h4-Tas3-s4-QE1-Tair,inTair,out,ECO,in=ECO,out+ΔECO.(c)
It is isenthalpic throttling process to refrigerant in throttle valve, so(17)QCO,out=QTH,out.Throttle valve exergy model is given by(18)ΔETH=m˙rfTas5-s4,ETH,in=ETH,out+ΔETH.(d)
Evaporator energy balance equation is shown as(19)QEV,in=QEV,out+QAB+QEV,loss.QAB is evaporator absorption heat from water (W).Evaporator exergy model is given by(20)ΔEEV=m˙rfh1-h5-Tas1-s5,EAB,in=QAB1-TEVTa,EEV,in+EAB,in=EEV,out+ΔEEV.Energy balance equations in cold exchanging and supplying process are written as(21)Qice,storage=Qair,fancoil+ΔQice,loss,Qair,fancoil=∫0tm˙air,fancoilCp,airTair,fancoil-in-Tair,fancoil-outdt,Qice,storage=micehice.And exergy calculation equation is shown:(22)Eice,in=Qair1-Tair,fancoil-outTair,fancoil-in+ΔEice.The energy conversion and transmission efficiencyη and exergy efficiency ψ of each component of the system can be calculated by the following formulas:(23)η=QoutQin,ψ=EoutEin.


## 2.1. ISACS Driven by DPES

### 2.1.1. Configuration of ISACS Driven by DPES

Ice storage air conditioning system (ISACS) driven by distributed photovoltaic energy system (DPES) was mainly configured by DPES, ice maker, storage system, and air conditioning system. The working diagram of ISACS driven by DPES is shown in Figure1.Figure 1
Work diagram of ISACS driven by DPES with batteries.PV modules convert solar energy into electric energy which can be regulated by controller with maximum power point tracking to drive ice maker, ice storage system, and air conditioning system. In daytime, DPES receives solar energy and turns it into direct-current (DC) electric power which can be converted to alternating current (AC) electric power by inverter to drive AC compressor, water pump, ethylene glycol pump, and fan coil. To maintain the stability of electric energy supply, batteries were adopted and connected with controller to maintain the energy conversion and supply in the most optimized way. Ice maker and storage system were made up with AC compressor, condenser, expansion valve, disc evaporator, and ice storage tank. Circulating water can be frozen in a disc evaporator and the ice can be drooped into the ice storage tank when the hot refrigerant flowed into the evaporator which was controlled by the solenoid valve. Thereby, the ice maker worked as vapor compression refrigeration. In AC compressor, cryogenic R134a was compressed to high temperature and high pressure gas to be filtered in economizer and to release heat in condenser. Refrigerant was condensed to mild temperature and high pressure gas. When the gas inflows into throttle valve, it can be throttled to low temperature and low pressure liquid and then feeds into plate evaporator. And then, refrigerant flows into the other gas-liquid separator to be sucked into the compressor. Thereby, the refrigeration cycle will be completed. Air conditioning system was mainly made up of coil heat exchanger which was fixed in ice storage tank, ethylene glycol pump, solenoid valve, proportional control valve, and fan coil. Ethylene glycol is an adopted cold exchanging medium.According to the working principle diagram, 0.2 kW ISACS driven by DPES was established as shown in Figure2. The main component parameters of ISACS driven by DPES are shown in Table 1.Table 1
Main component parameters of ISACS driven by DPES.
Component Model Parameters DPES PV module JN-245 P m: 245 W, Vm: 34.5 V, Im: 7.10 A, VOC: 43.5 V, ISC: 8.18 A, module: length∗wide: 1640 mm∗990 mm, cells in series: length∗wide∗numbers: 155 mm∗155 mm∗60 Controller PL60 12–48 V 60 A charge, 30 A load Inverter Solar 48 V P: 3 kW, DC input voltage: 48 V,output voltage: 220 V, output frequency: 50 Hz Batteries SP12-65 Battery capacity: 12 V 65 Ah, four batteries in series ISACS Refrigerant R134a Molecular formula: CH2FCF3, boiling point: −26.1°C, critical temperature: 101.1°C Ice maker IM50 Ice production: 2.12 kg/h,P: 380 W Ice storage tank / Capacity: 20 cm∗20 cm∗20 cm Cold exchanging medium Ethylene glycol Melting point: −12.6°C, viscosity: 25.66 mPa·s Pump RS15-6 Power: 46–93 W, life: 6 m, maximum flow rate: 3.4 m3/h Fan coil / Fan type: YS 56-2, power: 180 W, voltage: 380 V, current: 0.53 A, speed: 2800 r/min, number of fins: 95, size: 23 cm∗8 cm∗20 cm, coil numbers: 26, coil inner diameter: 6 mmFigure 2
Pictures of 0.2 kW ISACS driven by DPES.DPES were made up of two 245Wp polycrystalline silicon PV in series and four batteries in series were used to store electric power. Refrigerant (R134a) temperatures and pressure were measured by T-type thermocouples and pressure transducers, respectively. Voltages and currents of PV modules were measured by a digital multimeter. The wind speed was measured by the wind speed transducer. Solar irradiation was measured by pyranometer. Compressor input power was measured by a wattmeter. Electromagnetic flow meter was used to measure refrigerant flow and cold exchanging medium flow. The parameters of all instruments are shown in Table 2.Table 2
Measuring instrument parameters and uncertainty.
Instrument Model Range Accuracy Application scope Maximum relative error Maximum absolute error Uncertainty(B class) Pyranometer Kipp & Zonen CMP-6 0–2000 (W/m2) ±5% 0–1000 (W/m2) ±10% ±100 W/m2 57.7348 W/m2 Thermocouples T −200 to 350 (°C) ±0.4% 0–150 (°C) ±0.93% ±1.4°C 0.8083°C Wind speed transducer EC-9S 0–70 (m/s) ±0.4% 0–10 (m/s) ±2.8% ±0.28 m/s 0.1617 m/s Electromagnetic flow meter KROHNE OPTIFLUX 5000 DN 5; 0–12 (m/s) ±0.15% 0–5 (m/s) ±0.36% ±0.018 m/s 0.0104 m/s Pressure transducer YOKOGAWA EJA430E 0.14–16 (MPa) ±0.055% 0–2 (MPa) ±0.44% ±0.0088 MPa 0.0051 MPa Electronic balance AHW-3 0–3 kg ±0.05 g 0–3 kg ±0.05% ±0.0015 kg 0.0009 MPa Wattmeter DELIXI DDS607 0–10000 kW·h ±0.01 kW·h 0–100 kW·h ±1% ±1 kW·h 0.5774 kW·h Digital multimeter FLUKE F-179 Voltage: 0–1000 V ±0.9% 0–380 V ±2.37% ±9.006 V 5.1996 V Current: 0–10 A ±1% 0–10 A ±1% ±0.1 A 0.0577 A


## 2.1.1. Configuration of ISACS Driven by DPES

Ice storage air conditioning system (ISACS) driven by distributed photovoltaic energy system (DPES) was mainly configured by DPES, ice maker, storage system, and air conditioning system. The working diagram of ISACS driven by DPES is shown in Figure1.Figure 1
Work diagram of ISACS driven by DPES with batteries.PV modules convert solar energy into electric energy which can be regulated by controller with maximum power point tracking to drive ice maker, ice storage system, and air conditioning system. In daytime, DPES receives solar energy and turns it into direct-current (DC) electric power which can be converted to alternating current (AC) electric power by inverter to drive AC compressor, water pump, ethylene glycol pump, and fan coil. To maintain the stability of electric energy supply, batteries were adopted and connected with controller to maintain the energy conversion and supply in the most optimized way. Ice maker and storage system were made up with AC compressor, condenser, expansion valve, disc evaporator, and ice storage tank. Circulating water can be frozen in a disc evaporator and the ice can be drooped into the ice storage tank when the hot refrigerant flowed into the evaporator which was controlled by the solenoid valve. Thereby, the ice maker worked as vapor compression refrigeration. In AC compressor, cryogenic R134a was compressed to high temperature and high pressure gas to be filtered in economizer and to release heat in condenser. Refrigerant was condensed to mild temperature and high pressure gas. When the gas inflows into throttle valve, it can be throttled to low temperature and low pressure liquid and then feeds into plate evaporator. And then, refrigerant flows into the other gas-liquid separator to be sucked into the compressor. Thereby, the refrigeration cycle will be completed. Air conditioning system was mainly made up of coil heat exchanger which was fixed in ice storage tank, ethylene glycol pump, solenoid valve, proportional control valve, and fan coil. Ethylene glycol is an adopted cold exchanging medium.According to the working principle diagram, 0.2 kW ISACS driven by DPES was established as shown in Figure2. The main component parameters of ISACS driven by DPES are shown in Table 1.Table 1
Main component parameters of ISACS driven by DPES.
Component Model Parameters DPES PV module JN-245 P m: 245 W, Vm: 34.5 V, Im: 7.10 A, VOC: 43.5 V, ISC: 8.18 A, module: length∗wide: 1640 mm∗990 mm, cells in series: length∗wide∗numbers: 155 mm∗155 mm∗60 Controller PL60 12–48 V 60 A charge, 30 A load Inverter Solar 48 V P: 3 kW, DC input voltage: 48 V,output voltage: 220 V, output frequency: 50 Hz Batteries SP12-65 Battery capacity: 12 V 65 Ah, four batteries in series ISACS Refrigerant R134a Molecular formula: CH2FCF3, boiling point: −26.1°C, critical temperature: 101.1°C Ice maker IM50 Ice production: 2.12 kg/h,P: 380 W Ice storage tank / Capacity: 20 cm∗20 cm∗20 cm Cold exchanging medium Ethylene glycol Melting point: −12.6°C, viscosity: 25.66 mPa·s Pump RS15-6 Power: 46–93 W, life: 6 m, maximum flow rate: 3.4 m3/h Fan coil / Fan type: YS 56-2, power: 180 W, voltage: 380 V, current: 0.53 A, speed: 2800 r/min, number of fins: 95, size: 23 cm∗8 cm∗20 cm, coil numbers: 26, coil inner diameter: 6 mmFigure 2
Pictures of 0.2 kW ISACS driven by DPES.DPES were made up of two 245Wp polycrystalline silicon PV in series and four batteries in series were used to store electric power. Refrigerant (R134a) temperatures and pressure were measured by T-type thermocouples and pressure transducers, respectively. Voltages and currents of PV modules were measured by a digital multimeter. The wind speed was measured by the wind speed transducer. Solar irradiation was measured by pyranometer. Compressor input power was measured by a wattmeter. Electromagnetic flow meter was used to measure refrigerant flow and cold exchanging medium flow. The parameters of all instruments are shown in Table 2.Table 2
Measuring instrument parameters and uncertainty.
Instrument Model Range Accuracy Application scope Maximum relative error Maximum absolute error Uncertainty(B class) Pyranometer Kipp & Zonen CMP-6 0–2000 (W/m2) ±5% 0–1000 (W/m2) ±10% ±100 W/m2 57.7348 W/m2 Thermocouples T −200 to 350 (°C) ±0.4% 0–150 (°C) ±0.93% ±1.4°C 0.8083°C Wind speed transducer EC-9S 0–70 (m/s) ±0.4% 0–10 (m/s) ±2.8% ±0.28 m/s 0.1617 m/s Electromagnetic flow meter KROHNE OPTIFLUX 5000 DN 5; 0–12 (m/s) ±0.15% 0–5 (m/s) ±0.36% ±0.018 m/s 0.0104 m/s Pressure transducer YOKOGAWA EJA430E 0.14–16 (MPa) ±0.055% 0–2 (MPa) ±0.44% ±0.0088 MPa 0.0051 MPa Electronic balance AHW-3 0–3 kg ±0.05 g 0–3 kg ±0.05% ±0.0015 kg 0.0009 MPa Wattmeter DELIXI DDS607 0–10000 kW·h ±0.01 kW·h 0–100 kW·h ±1% ±1 kW·h 0.5774 kW·h Digital multimeter FLUKE F-179 Voltage: 0–1000 V ±0.9% 0–380 V ±2.37% ±9.006 V 5.1996 V Current: 0–10 A ±1% 0–10 A ±1% ±0.1 A 0.0577 A


## 2.2. Theoretical Models on Energy Conversion and Transfer Characteristics of ISACS Driven by DPES

Theoretical model was established to analyze the energy conversion and transfer characteristics of ISACS driven by DPES as follows.


### 2.2.1. Energy and Exergy Models of PV Modules

Energy balance equation of PV modules is expressed as(1)mpCpmoduleTp=QPV,in-QPV,loss-QPV,elect.Solar energy absorbed by PV modules can be estimated as(2)QPV,in=ατGSp.Electric power is expressed by(3)QPV,elect=VoutIout.Exergy and exergy losses of PV modules are given in [26]: (4)ΔEPV=QPV,loss1-TaTP+IscVoc-ImVm,EPV,in=GSp1-43TaTsun+13TaTsun4,EPV,in=EPV,out+ΔEPV.


### 2.2.2. Energy and Exergy Models of Controller

Energy balance equation of controller can be estimated as follows:(5)QPV,elect=QC,out+QC,loss.QC,loss is the consumed energy by controller per unit time (W); usually QC,loss is 4% of QPV,elect.Exergy and exergy losses of controller are expressed as in [27]:(6)ΔEC=QC,loss1-TaTC,EC,in=EC,out+ΔEC.


### 2.2.3. Energy and Exergy Models of Batteries

Energy balance equations of batteries are given by [28](7)QB,in=ξQB,out+QB,loss-ξQB,storage.ξ is a symbolic coefficient. When batteries supply power to load, ξ is assigned 1. On the contrary, ξ is assigned −1 when batteries are charged:(8)QB,loss=IB2r1+r2SOC+1r3-r4SOC.Here,(9)SOC=SOC0+∫0tIBdtBC.VF is the full charge rest voltage (V), SOC0 is the initial SOC value (1), and BC is batteries capacity (Ah).Exergy balance and exergy loss equations of batteries are written as [27](10)ΔEB=QB,loss1-TaTB,EB,in=EB,out+ΔEB.


### 2.2.4. Energy and Exergy Models of Inverter

Energy balance empirical equations of inverter are given in [28]:(11)QI,loss=CPCR+CR-1CRQC,out+QB,out,QC,out+QB,out=QI,out+QI,loss.Exergy balance and exergy loss equations of inverter are expressed in [28]: (12)ΔEI=QI,loss1-TaTI,EC,out+EB,out=EI,out+ΔEI.


### 2.2.5. Energy and Exergy Models of ISACS

ISACS are driven by stable electric energy outputted by DPES. The thermodynamic cycle refrigeration process of ice maker, refrigerant, R134a thermodynamic properties, andP-h changing curve of refrigerant R134a in refrigeration cycle are shown in Figure 3. It can be clearly seen in Figure 3 that area 1 is the saturated liquid, area 2 is the gas-liquid coexistence, and area 3 is the saturated vapor. Consequently, the refrigeration cycle was 1-2-2s-3-4-5-1 across the saturated vapor region and gas-liquid coexistence region. The process from 1-2 is showing an actual compressor working process, when the dry steam in state 1 is compressed to saturated vapor in state 2, which will be isentropic compressed to 2s under ideal condition. However, the process 2-2s-3-4 is the refrigerant condensing and releasing heat process in condenser. The process 2-2s is the natural cooling process and 2s-3-4 is the isobaric cooling process. The cryogen will be high pressure dry steam in state 3 and high pressure saturated liquid in state 4. The process 4-5 is the decompression cooling by throttle valve. Therefore, by flowing through throttle valve, the refrigerant temperature and pressure will be low and then it can flow into evaporator to absorb heat and refrigerate, which is the process 5-1.Figure 3
R134a thermodynamic curve.Energy and exergy analyses are shown as follows.(a)
The compressor energy balance equation is given by(13)m˙rfh1+W=m˙rfh2+QCP,loss.Compressor exergy model is given as(14)ECP,in=ΔECP+ECP,out,ΔECP=m˙rfh2-h1-Tas2-s1.(b)
The condenser energy balance equation is shown as(15)m˙rfh3=m˙rfh4+QE+QCO,loss,QE=m˙air,icemkCp,airTair,out-Tair,in.Condenser exergy model is expressed as(16)ΔECO=m˙rfh3-h4-Tas3-s4-QE1-Tair,inTair,out,ECO,in=ECO,out+ΔECO.(c)
It is isenthalpic throttling process to refrigerant in throttle valve, so(17)QCO,out=QTH,out.Throttle valve exergy model is given by(18)ΔETH=m˙rfTas5-s4,ETH,in=ETH,out+ΔETH.(d)
Evaporator energy balance equation is shown as(19)QEV,in=QEV,out+QAB+QEV,loss.QAB is evaporator absorption heat from water (W).Evaporator exergy model is given by(20)ΔEEV=m˙rfh1-h5-Tas1-s5,EAB,in=QAB1-TEVTa,EEV,in+EAB,in=EEV,out+ΔEEV.Energy balance equations in cold exchanging and supplying process are written as(21)Qice,storage=Qair,fancoil+ΔQice,loss,Qair,fancoil=∫0tm˙air,fancoilCp,airTair,fancoil-in-Tair,fancoil-outdt,Qice,storage=micehice.And exergy calculation equation is shown:(22)Eice,in=Qair1-Tair,fancoil-outTair,fancoil-in+ΔEice.The energy conversion and transmission efficiencyη and exergy efficiency ψ of each component of the system can be calculated by the following formulas:(23)η=QoutQin,ψ=EoutEin.


## 2.2.1. Energy and Exergy Models of PV Modules

Energy balance equation of PV modules is expressed as(1)mpCpmoduleTp=QPV,in-QPV,loss-QPV,elect.Solar energy absorbed by PV modules can be estimated as(2)QPV,in=ατGSp.Electric power is expressed by(3)QPV,elect=VoutIout.Exergy and exergy losses of PV modules are given in [26]: (4)ΔEPV=QPV,loss1-TaTP+IscVoc-ImVm,EPV,in=GSp1-43TaTsun+13TaTsun4,EPV,in=EPV,out+ΔEPV.


## 2.2.2. Energy and Exergy Models of Controller

Energy balance equation of controller can be estimated as follows:(5)QPV,elect=QC,out+QC,loss.QC,loss is the consumed energy by controller per unit time (W); usually QC,loss is 4% of QPV,elect.Exergy and exergy losses of controller are expressed as in [27]:(6)ΔEC=QC,loss1-TaTC,EC,in=EC,out+ΔEC.


## 2.2.3. Energy and Exergy Models of Batteries

Energy balance equations of batteries are given by [28](7)QB,in=ξQB,out+QB,loss-ξQB,storage.ξ is a symbolic coefficient. When batteries supply power to load, ξ is assigned 1. On the contrary, ξ is assigned −1 when batteries are charged:(8)QB,loss=IB2r1+r2SOC+1r3-r4SOC.Here,(9)SOC=SOC0+∫0tIBdtBC.VF is the full charge rest voltage (V), SOC0 is the initial SOC value (1), and BC is batteries capacity (Ah).Exergy balance and exergy loss equations of batteries are written as [27](10)ΔEB=QB,loss1-TaTB,EB,in=EB,out+ΔEB.


## 2.2.4. Energy and Exergy Models of Inverter

Energy balance empirical equations of inverter are given in [28]:(11)QI,loss=CPCR+CR-1CRQC,out+QB,out,QC,out+QB,out=QI,out+QI,loss.Exergy balance and exergy loss equations of inverter are expressed in [28]: (12)ΔEI=QI,loss1-TaTI,EC,out+EB,out=EI,out+ΔEI.


## 2.2.5. Energy and Exergy Models of ISACS

ISACS are driven by stable electric energy outputted by DPES. The thermodynamic cycle refrigeration process of ice maker, refrigerant, R134a thermodynamic properties, andP-h changing curve of refrigerant R134a in refrigeration cycle are shown in Figure 3. It can be clearly seen in Figure 3 that area 1 is the saturated liquid, area 2 is the gas-liquid coexistence, and area 3 is the saturated vapor. Consequently, the refrigeration cycle was 1-2-2s-3-4-5-1 across the saturated vapor region and gas-liquid coexistence region. The process from 1-2 is showing an actual compressor working process, when the dry steam in state 1 is compressed to saturated vapor in state 2, which will be isentropic compressed to 2s under ideal condition. However, the process 2-2s-3-4 is the refrigerant condensing and releasing heat process in condenser. The process 2-2s is the natural cooling process and 2s-3-4 is the isobaric cooling process. The cryogen will be high pressure dry steam in state 3 and high pressure saturated liquid in state 4. The process 4-5 is the decompression cooling by throttle valve. Therefore, by flowing through throttle valve, the refrigerant temperature and pressure will be low and then it can flow into evaporator to absorb heat and refrigerate, which is the process 5-1.Figure 3
R134a thermodynamic curve.Energy and exergy analyses are shown as follows.(a)
The compressor energy balance equation is given by(13)m˙rfh1+W=m˙rfh2+QCP,loss.Compressor exergy model is given as(14)ECP,in=ΔECP+ECP,out,ΔECP=m˙rfh2-h1-Tas2-s1.(b)
The condenser energy balance equation is shown as(15)m˙rfh3=m˙rfh4+QE+QCO,loss,QE=m˙air,icemkCp,airTair,out-Tair,in.Condenser exergy model is expressed as(16)ΔECO=m˙rfh3-h4-Tas3-s4-QE1-Tair,inTair,out,ECO,in=ECO,out+ΔECO.(c)
It is isenthalpic throttling process to refrigerant in throttle valve, so(17)QCO,out=QTH,out.Throttle valve exergy model is given by(18)ΔETH=m˙rfTas5-s4,ETH,in=ETH,out+ΔETH.(d)
Evaporator energy balance equation is shown as(19)QEV,in=QEV,out+QAB+QEV,loss.QAB is evaporator absorption heat from water (W).Evaporator exergy model is given by(20)ΔEEV=m˙rfh1-h5-Tas1-s5,EAB,in=QAB1-TEVTa,EEV,in+EAB,in=EEV,out+ΔEEV.Energy balance equations in cold exchanging and supplying process are written as(21)Qice,storage=Qair,fancoil+ΔQice,loss,Qair,fancoil=∫0tm˙air,fancoilCp,airTair,fancoil-in-Tair,fancoil-outdt,Qice,storage=micehice.And exergy calculation equation is shown:(22)Eice,in=Qair1-Tair,fancoil-outTair,fancoil-in+ΔEice.The energy conversion and transmission efficiencyη and exergy efficiency ψ of each component of the system can be calculated by the following formulas:(23)η=QoutQin,ψ=EoutEin.


## 3. Results and Discussion

### 3.1. Calculation Results and Analysis

The related parameters are shown in Table3.Table 3
Related parameters of ISACS driven by DPES.
Parameters Value Outside Ts/K 5778 Ta/K 286.36–298 v/(m·s−1) 0.87–1.55 q/MJ·m−2 (08:00–16:00) 22.17 PV modules Q/kW·h (08:00–16:00) 2.76 Sp/m2 2.88 Sc/m2 3.24 VOC/V 89.35 ISC/A 8.18 Vout/V 50.90–54.20 Iout/A 0.50–8.50 Vm/V 34.50 Im/A 7.10 η0/% 17.50 α 0.92 τ 0.90 m/(kg·m−2) 23 Cpmodule/(J·kg−1·K−1) 1179.06 TP/K 287.04–325.15 Controller TC/K 303.15 Storage batteries C 48 V 65 Ah TB/K 313.15 Charging B 0.80 VF 56.53 VB/V 47.60–52.53 IB/A 8.00 r1 0.060 r2 0.041 r3 95.234 r4 51.856 Discharging B 0.70 VF 54.88 VB/V 54.88–44.60 IB/A 5.50–7.40 r1 0.052 r2 −0.012 r3 4.113 r4 −100.653 Inverter TI/K 318.15 Cp 10.045 CR 1.1885 Compressor m˙/(kg·s−1) 0.0127 Wp/W 328.31–340.25 Tin/K 268.15 Tout/K 313.15 Pin/kPa 243.71 Pout/kPa 770.21 hin/(kJ·kg−1) 395.01 hout/(kJ·kg−1) 425.00 sin/(kJ·kg−1·K−1) 1.7276 sout/(kJ·kg−1·K−1) 1.7500 Condenser Tin/K 313.15 Tout/K 303.15 Pin/kPa 700.00 Pout/kPa 700.00 hin/(kJ·kg−1) 425.00 hout/(kJ·kg−1) 241.80 sin/(kJ·kg−1·K−1) 1.7500 sout/(kJ·kg−1·K−1) 1.1437 m˙air,icemk/(kg·s−1) 0.21 Cp,air/(J·kg−1·K−1) 1000 Tair,in/K 286.36–296 Tair,out/K 291.36–301 Throttle valve Tin/K 303.15 Tout/K 263.15 Pin/kPa 700.00 Pout/kPa 200.00 hin/(kJ·kg−1) 241.80 hout/(kJ·kg−1) 241.80 sin/(kJ·kg−1·K−1) 1.1437 sout/(kJ·kg−1·K−1) 1.1500 Evaporator Tin/K 263.15 Tout/K 268.15 Pin/kPa 200.00 Pout/kPa 243.71 hin/(kJ·kg−1) 241.80 hout/(kJ·kg−1) 395.01 sin/(kJ·kg−1·K−1) 1.1500 sout/(kJ·kg−1·K−1) 1.7276 Cold exchanging and supplying Tice-tank/K 272.20–278.15 Tglycol/K 276.62–285.15 Tair,fancoil-in/K 293.45–289.95 Tair,fancoil-out/K 293.45–287.65 mICE/kg 16.98 hICE/(kJ·kg−1) 335 m˙air,fancoil/(kg·s−1) 0.12 tcold/s 1440008:00–17:00 is the working time from Monday to Friday. During working time, the cold demand in the house is very low, so the refrigerator can be driven by PV modules and then the cold can be stored in the ice storage tank. After work, the cold demands of the house increased and the cold stored in the daytime can be charged to service family. So the ISACS driven by DPES has great application prospects in domestic cooling field.The ISACS driven by DPES operated from 08:00 to 23:00. In daytime, solar energy was converted into electrical energy by PV modules. Before 11:00 AM, the irradiation was too low and the electricity generated by PV modules cannot drive ice maker. At the same time, batteries discharged to make up for the lack of electricity produced by PV modules and the ice maker could operate stably and reliably driven by photovoltaic battery hybrid energy supply system. And then, the irradiation gradually increased and the electricity increased along with irradiance. In hybrid energy supply system, the system output power was equal to the power rating of ice maker and was always the same but PV modules output power increased and batteries output power decreased accordingly. About 11 AM, the output power of PV modules was enough to drive ice maker so the output power of batteries is zero. After 11 AM, with the continuous increase of irradiance, not only was the electricity of PV modules used to drive ice maker but also the surplus electricity was stored in batteries. So the batteries were in charge. After noon, the irradiation from sun decreased step by step and so the electricity generated by PV modules reduced gradually. About 14:00, the electricity reducing step by step was only enough to drive ice maker. The charge state of batteries was end. Then, ice maker was driven by PV modules and batteries once again until the ice was enough at 16:00. All of the ice was stored in ice storage tank and the cold can be exchanged by the coil with refrigerating medium ethylene glycol flowing in it at 19:00. At the other end of the coil the cold of refrigerating medium can be blown into the air in the house by fan. So the temperature of house dropped step by step until it tended to balance. About 23:00, the cold of ice in the tank was released completely. And then, the fan and glycol pump would be shut down and the cold supply process was over. Through all of the process, energy is transformed from light to electricity by PV modules and it is stored with ice through phase change latent heat of water. In conversion and transmission process, the energy and exergy of system changed with the change of external environmental conditions. In order to describe system performance clearly and intuitively, the overall performance of the system was evaluated. So the energy efficiency and exergy efficiency of each component of the system in energy conversion and transmission process were calculated and the results are shown in Table4.Table 4
Calculated results of total energy and total exergy of each component of the system in energy conversion and transmission process.
Accepted energy/W Accepted exergy/W Useful power/W Useful exergy/W η/% ψ/% PV module 100007.5 112560.5 15889.79 95910.39 15.89 85.21 Controller 15889.79 95910.39 15254.2 95894.55 96.00 99.98 Batteries 15254.2 95894.55 15162.49 95887.94 99.40 99.99 Inverter 19772.17 95887.94 16222.1 95621.2 82.05 99.72 Distributed photovoltaic energy system 12.44 84.95 Compressor 156936.5 95621.2 151397.8 87286.14 96.47 91.28 Condenser 151397.8 87286.14 Air 51450 86445.17 90.88 99.04 Refrigeration 86136.41 Throttle valve 86136.41 86445.17 86136.41 85784.72 100 99.24 Evaporator 86136.41 85784.72 Output 140714.4 79810.74 336 93.04 Absorbing 54578 Ice 54578 79810.74 9484.22 76271.06 17.38 95.57 Ice maker system 51.20 79.77 Air conditioning system 94842.22 76271.06 81009.64 75747.45 85.42 99.31 ISACS driven by DPES 5.44 67.30Figure4 showed the change curves of energy efficiency and exergy efficiency of ice making and cold-storing process from 08:00 to 16:00 with ice maker driven by DPES. Figure 5 showed the change curves of system energy efficiency and exergy efficiency of fan coil driven by DPES from 19:00 to 23:00.Figure 4
Energy efficiency and exergy efficiency of ice maker driven by DPES.Figure 5
Energy efficiency and exergy efficiency of fan coil driven by DPES.In daytime, ice maker system was driven by DPES batteries and the system energy efficiency and exergy efficiency were 5.44% and 67.30%, respectively. It was found that the total energy efficiency and total exergy efficiency of DPES were 12.44% and 84.95%, respectively, and the total energy efficiency and total exergy efficiency of ice maker system were 51.20% and 79.77%, respectively. After 16:00, the total ice production was 16.98 kg. All of the ice was stored in the tank and air conditioner was open at 19:00 and sustained for about 4 hours with 85.42% total energy efficiency and 99.31% exergy efficiency.It was observed that in ice maker thermodynamic cycle refrigeration efficiency was only 17.38% and the wastage in cooling process was 45093.78 W. In the ice making process, the water can be pumped by a recycle pump and flowed through evaporator to absorb heat. A part of water can be frozen on the evaporator and the rest of water flowed back to the water storage. Thereby, in the water circulating process, the circulating water consumed much of the energy, as shown in Figure6.Figure 6
Ice slide maker (internal picture).In ice maker, compressor was the maximum exergy loss component with 8.78% exergy loss efficiency and it is clear that the evaporator is the most important component of ice maker; the total exergy of evaporator was 79810.74 W and the exergy efficiency was 93.04%. At the same time, the exergy flow from water to the evaporator was 76271.06 W with 95.57% exergy efficiency in ice making process.


### 3.2. Experimental Results and Analysis

#### 3.2.1. DPES

The output performance of PV modules was greatly affected by the external environment such as solar irradiance, wind speed, ambient temperature, and PV modules temperature. Therefore, the external environment should be tested through the experiment to calculate and analyze PV output performance. ISACS driven by DPES tested on October 22 in Kunming. Results showed that the irradiation quantityq was 22.17 MJ/m2 from 8:00 to 16:00 and the total amount of ice was 16.98 kg when the fan coil operated for nearly 4 hours. Meanwhile, PV modules output performances were tested and analyzed. Firstly, global solar radiation and ambient temperature were tested, as shown in Figure 7. PV modules temperature is a decisive component of output performance, which can be affected by wind speed changing at any time. Thus, wind speed and PV modules temperature can only be obtained by experimental measurement, as shown in Figure 8.Figure 7
Variations of global solar radiation and ambient temperature on 22nd Nov., 2015, in Kunming.Figure 8
Variations of wind speed and PV modules temperature on 22nd Nov., 2015, in Kunming.The voltage and current variations of PV modules were also tested on October 22 in Kunming. The variation curves are shown in Figure9. Batteries were adopted in testing process to ensure power output stability and store excess power generated by PV modules. The voltage and current of batteries were also measured and variation curves are also shown in Figure 10.Figure 9
Voltage and current variations of PV modules.Figure 10
Voltage and current variations of batteries.It was observed that the maximum tested current was 8.6 A at 12:35. The output voltage changes were 50.8 V–54.2 V and the maximum voltage appears at 12:36. In batteries discharging, ice maker can be driven for 10 hours by batteries only which were in full power state and disconnected with PV modules until the batteries electricity decreased and cannot drive the ice maker. The voltage declined from 54.88 V to 44.6 V step by step and the output current of PV modules was approximately 7 A with small fluctuations between 5.4 A and 7.4 A. In charging process for batteries, batteries disconnect with ice maker and connected with PV modules. Batteries current changed with solar irradiance and voltage increased gradually until batteries were in floating state.Electricity conversed from light by distributed photovoltaic energy system was transferred among batteries, ice maker, and air conditioning system. It was very important to optimize and match the energy among PV modules, storage batteries, and power consumption machine, which can provide an important reference for the future work. Therefore, the experimental study on ISACS driven by DPES was carried out in order to analyze the light-electricity-cold conversion process and the energy transfer characteristics. The results are shown in Figure11.Figure 11
The conversion and flow characteristics of power generated by PV modules among batteries, ice maker, and air conditioning system.There were ice making and ice melting abscission in ice maker operation process. In ice making process, the system current was 1.8 A, the sum of compressor working current and condenser fan current. However, in ice melting abscission process, the system current was 1.4 A, the sum of condenser fan current and solenoid valve operating current. On the other hand, the voltage of AC compressor was 220 V and it was stable all the time. Voltage and current of pump and fan were 220 V and 0.4 A and 380 V and 0.5 A, respectively. In DPES, photovoltaic modules and batteries were untied and used to drive ice maker. Before noon the PV modules output current increases along with solar irradiance and batteries output current gradually decreases. When PV modules output current increased to 6.5 A, batteries output current reduced to 0 A at 10:50 and ice maker can be completely driven by PV modules. Since then, ice maker was driven by PV module and the rest electricity can be stored in batteries. Batteries charging current increased step by step and reached the maximum of 2.0 A until the solar irradiance reached the peak values at 12:20. After noon, batteries charging current gradually reduces along with solar irradiance decreases. It was found that at 15:37 the charging current and PV modules output current were 0 A and 6.5 A. Therefore, ice maker was only driven by PV modules once again. Subsequently, the generated electricity by PV modules was not enough to drive the ice maker and batteries discharged to supplement the shortage. Batteries discharging current increases step by step along with the solar irradiance decreases. Batteries voltage reduced in discharging process and increased in charging process. Batteries separately derived pump and fan coil to supply cold from 19:00 to 23:00. Output current remained constant with 5 A and output voltage slightly decreased gradually. At 23:00, ice melts completely and then all machines were shut down. Until 07:30 of the next day, a new ice making and cold supplying cycle starts. Through the 27-hour experiment, batteries capacity decreased from 65 Ah to 9.76 Ah with a 55.24 Ah decrease. In the first day, batteries discharged 13.7 Ah in the morning and then were charged 7.39 Ah by PV modules. Batteries discharged 12.17 Ah from 15:37 to 19:00. And then ice maker stopped and cold exchanger ran, driven by batteries. The electricity consumption was 20.22 A from 19:00 to 23:00. In the next day, ice maker starts to work at 07:00 and electricity consumption in the morning discharging process and afternoon discharging process was 12.55 Ah and 11.61 Ah, respectively. Batteries were charged 7.62 Ah by PV modules from 10:11 to 15:36.


#### 3.2.2. Ice Maker and Cold Storage System

The ice production is shown in Figure12.Figure 12
Ice production of ice maker driven by DPES from 8:00 to 16:00 on 22nd Nov., 2015, in Kunming.It was found that the total amount of ice was 16.98 kg. Ice maker operation cycle time was 10 min and the ice production was about 0.35~0.36 kg every cycle. Ice making efficiency was about 2.12 kg/h. Another experimental test has been conducted at various temperatures of ice maker with every component in thermodynamic cycle of ice making process, as shown in Figure13. The inlet temperature and outlet temperature of evaporator, compressor, and condenser declined stably until ice making process was performed.Figure 13
Temperature variations of ice maker components during two refrigeration thermodynamic cycle processes.In ice making process, ice was pasted well together with five solid walls of grid plate evaporator and cannot be separated off by ice gravity. Thus, compressor must be shut down and a special solenoid valve as shown in Figure14 must be opened which was used for adjusting refrigerant flow and introduced the high temperature and high pressure refrigerant steam expelled from compressor into evaporator.Figure 14
Electromagnetic valve photo adopted in ice maker system.When the electromagnetic was opened, the evaporator plays the role of condenser to release heat to ice through solid walls and ice interfaces begin melting. Ice was divorced from evaporator and fell into tank. Meanwhile, evaporator inlet temperature and outlet temperature increased sharply until the solenoid valve was closed and compressor was turned on and a new ice making process begins. The refrigerant flowed out of evaporator into throttle valve and became low temperature and low pressure liquid. And then, the refrigerant flowed into condenser to absorb heat from outside and condenser inlet and outlet temperature declined sharply when low temperature refrigerant flowed into it. Condenser outlet temperature was higher than condenser inlet because refrigerant absorbed heat from outside through condenser, as shown in Figure12. Refrigerant next flowed into compressor and compressor inlet temperature increased a lot compared to ice making process, which can be changed with inlet and outlet temperature variation of evaporator and condenser. Compressor stops running in ice melting abscission process and the refrigerant gas cannot be compressed in compressor. Therefore, the compressor outlet temperature drops sharply when the solenoid valve opens and it returns smoothly when ice making process opened. Ice melting abscission process was extremely unfavorable to ISACS for a number of reasons as follows:(1)
Ice maker operation period extended for 200 s; furthermore, the ice melting abscission time was one-third time of an ice making cycle.(2)
Service life of compressor will be shortened by frequent start and stop compressor, which has a great impact on ice making process and was extremely unfavorable for energy supply process of DPES. Electricity consumption of solenoid valve was 3.0 × 10−4 kW h in ice melting abscission process and the total electricity consumption was 0.024 kW h from 7:30 to 19:00.


#### 3.2.3. Air Conditioning System

The performance of air conditioning system was tested at 19:00. At this time, the temperature of ice stored in ice storage tank and ice storage tank temperature were −3°C and −1°C, respectively. The indoor temperature was 20.5°C. The changes of temperature for exchanging cold and supplying cold processes are shown in Figure15.Figure 15
Temperatures change of exchanging cold and supplying cold process.Ice storage temperature remained constant at −1°C in ice absorbing and melting phase transition processes. However, after 22:20, all ice melted and the tank filled with water at low temperature. Subsequently, the temperature of water gradually increased when the water was employed to cool. Subsequently, it increased from 0°C to nearly 5°C within 40 min. At the end of the experiment, ice storage tank outlet temperature increased to 11.5°C from initial temperature of 3.4°C. In ice phase transition process, ice storage tank inlet and outlet temperatures should be constant values in theory; however, the temperature difference between inlet and outlet increased 7.15°C mainly as a result of the working temperature of pipeline pump. It was found that the fan outlet temperature decreased from 20.5°C to 15.1°C with a 5.4°C decline and the indoor temperature decreased from 20.5°C to 18°C with a 2.5°C decline within 4 hours. Cold power blown out by fan coil was 0.298 W and the cold exchanging and supplying efficiency of the air conditioning system was 90.4%. Furthermore, the ice storage tank has a good thermal insulation performance at 0.10 m thick polyurethane foam. The ice cold loss in ice storage tank was 218.59 kJ as 3.89% of input cold from 19:00 to 23:00 and the cold loss of fan coil was 313.92 kJ as 5.82% of input cold.


### 3.3. Model Validation

The system energy efficiency was 5.44% calculated as shown in Table4. In the experimental test process, the ice produce was 16.98 kg during 08:00~16:00 with the 22.17 MJ/m2 total irradiation. Because the ice maker was driven by PV modules and batteries, the ice produce was not only produced by PV modules. Therefore, the ice production of PV modules must be calculated. The rated power of ice maker was 380 W as shown in Table 1 which was not only made up with the power of compressor as shown in Table 3 but also included the powers of fan, solenoid valve, and water pump. The electricity consumption of ice maker was 3.04 kW·h during 08:00~16:00. But the electricity Q produced by PV modules was only 2.76 kW·h as shown in Table 3. So there was 90.79% power consumed by ice maker coming from PV modules. So the ice produced by PV modules was 15.42 kg. The ice phase change latent heat was 335 kJ/kg. So the energy stored by ice produced by PV modules was 5165.70 kJ. The total area of PV modules was 3.24 m2. So the total solar energy accepted by PV modules was 71830.80 kJ. And the cold exchange efficiency of air conditioning system was 85.42% as shown in Table 4. So the energy efficiency of ISACS driven by DPES tested in the experiment was 6.14%. The relative error of the calculated result was 11.4%. The model can be accepted in engineering application fields.


### 3.4. System Optimization

Analog computation and experimental tests results showed that there were some deficiencies of current system as follows:(1)
Usually the cold demand during the daytime is high, but the system can supply cold through ice melting made by ice maker driven by electricity generated by PV modules in daytime. Therefore, there is no matching between cold supply and demand in the same time.(2)
Energy loss is huge and ice making efficiency is low in ice making process of ice slide maker.Therefore, in order to improve the performance of ISACS driven by DPES and promote the project commercial promotion, some optimization and improvement measurements were proposed as follows:(a)
Evaporator immersion static refrigeration mode was adopted to replace ice harvester refrigeration mode. The optimized coil evaporator was immersed into water to absorb heat and make ice and all the energy was utilized, as shown in Figure16.(b)
The coil cold exchanger was cointegrated with coil evaporator. In refrigeration process, coil cold exchanger has the priority to get cold transferred from the coil evaporator next to the cold exchanger to supply cold for user. Surplus cold is used to make ice to store cold. Consequently, ISACS not only has the out-of-the-box function of ordinary air conditioning, but also effectively improves the appearing phenomenon of overcooling and remedies the disadvantage of cold supply after ice making process in traditional submerged ice making system. The top view of cointegration evaporators and cold exchanger immersed in the ice storage tank is shown in Figure17.Figure 16
Profile of evaporators and cold exchanger immersed in the ice storage tank.Figure 17
Top view of cointegration evaporators and cold exchanger immersed in the ice storage tank.Another immersion static ice maker with evaporator cointegrated with coil cold exchanger was constructed according to the optimization design shown in Figure18.Figure 18
The photo of evaporator immersion static ice maker.The performance of optimized ISACS driven by DPES can be tested and calculated as mentioned above. The system operated on May 12, 2016. The tested and calculated results were shown in Figure19 and Table 5.Table 5
The performance of two kinds of ISACS driven by DPES.
DPES Ice maker system Air conditioning system System Compressor Condenser Throttle valve Evaporator Ice Total η / % ψ / % η / % ψ / % η / % ψ / % η / % ψ / % η / % ψ / % η / % ψ / % η / % ψ / % η / % ψ / % η / % ψ / % 1 12.44 84.95 96.47 91.28 90.88 99.04 100 99.24 336 93.04 17.38 93.04 51.20 79.77 85.42 99.31 5.44 67.30 2 12.44 84.95 96.47 91.28 90.88 99.04 100 99.24 336 93.04 50.12 97.32 147.64 81.24 85.42 99.31 15.69 68.54 Note: 1 stands for the ISACS driven by DPES mentioned above; 2 stands for the optimized ISACS driven by DPES.Figure 19
Temperature changes of ice maker components after optimization.The water tank was full of 13.4 kg water. The system operated at 08:35 and all of the water changed into ice at about 10:40. And then the ice started subcooling and the ice surface became hard at 12:04. At the moment, the surface of evaporator temperature (the temperature of ice core) was −12.11°C and the water temperature was −8.09°C with a 4.20°C temperature difference between ice surface and ice core. At 12:04, the temperatures of evaporator inlet and outlet were, respectively, −20.49°C and −16.51°C and the compressor inlet temperature was −14.28°C. The energy efficiency and exergy efficiency of system components of optimized ISACS driven by DPES were calculated through above formulas and the performances of the system mentioned above and the optimized system were compared in Table5.By structure optimization, the ice making efficiency of evaporator immersion static ice maker was improved to 6.00 kg/h which was 2.88 times of ice slide maker with the same input power. Through theoretical calculation, it was observed that the ice making efficiency and exergy efficiency increased from 17.38% to 50.12% and from 95.57% to 97.32%, respectively. So the ice maker system energy efficiency was improved 2.88 times and exergy efficiency was improved from 76.85% to 78.52%. After the structure optimization, the energy efficiency of ISACS driven by DPES can be improved from 5.44% to 15.69% and the exergy efficiency can be improved from 67.30% to 68.54%. When evaporator cointegrated with coil cold exchanger, ISACS achieved out-of-box functionality as the ordinary vapor compression air conditioning system, which can effectively solve the problem of cold supplying process must lagging behind ice making process in traditional ice storage air conditioning system.


## 3.1. Calculation Results and Analysis

The related parameters are shown in Table3.Table 3
Related parameters of ISACS driven by DPES.
Parameters Value Outside Ts/K 5778 Ta/K 286.36–298 v/(m·s−1) 0.87–1.55 q/MJ·m−2 (08:00–16:00) 22.17 PV modules Q/kW·h (08:00–16:00) 2.76 Sp/m2 2.88 Sc/m2 3.24 VOC/V 89.35 ISC/A 8.18 Vout/V 50.90–54.20 Iout/A 0.50–8.50 Vm/V 34.50 Im/A 7.10 η0/% 17.50 α 0.92 τ 0.90 m/(kg·m−2) 23 Cpmodule/(J·kg−1·K−1) 1179.06 TP/K 287.04–325.15 Controller TC/K 303.15 Storage batteries C 48 V 65 Ah TB/K 313.15 Charging B 0.80 VF 56.53 VB/V 47.60–52.53 IB/A 8.00 r1 0.060 r2 0.041 r3 95.234 r4 51.856 Discharging B 0.70 VF 54.88 VB/V 54.88–44.60 IB/A 5.50–7.40 r1 0.052 r2 −0.012 r3 4.113 r4 −100.653 Inverter TI/K 318.15 Cp 10.045 CR 1.1885 Compressor m˙/(kg·s−1) 0.0127 Wp/W 328.31–340.25 Tin/K 268.15 Tout/K 313.15 Pin/kPa 243.71 Pout/kPa 770.21 hin/(kJ·kg−1) 395.01 hout/(kJ·kg−1) 425.00 sin/(kJ·kg−1·K−1) 1.7276 sout/(kJ·kg−1·K−1) 1.7500 Condenser Tin/K 313.15 Tout/K 303.15 Pin/kPa 700.00 Pout/kPa 700.00 hin/(kJ·kg−1) 425.00 hout/(kJ·kg−1) 241.80 sin/(kJ·kg−1·K−1) 1.7500 sout/(kJ·kg−1·K−1) 1.1437 m˙air,icemk/(kg·s−1) 0.21 Cp,air/(J·kg−1·K−1) 1000 Tair,in/K 286.36–296 Tair,out/K 291.36–301 Throttle valve Tin/K 303.15 Tout/K 263.15 Pin/kPa 700.00 Pout/kPa 200.00 hin/(kJ·kg−1) 241.80 hout/(kJ·kg−1) 241.80 sin/(kJ·kg−1·K−1) 1.1437 sout/(kJ·kg−1·K−1) 1.1500 Evaporator Tin/K 263.15 Tout/K 268.15 Pin/kPa 200.00 Pout/kPa 243.71 hin/(kJ·kg−1) 241.80 hout/(kJ·kg−1) 395.01 sin/(kJ·kg−1·K−1) 1.1500 sout/(kJ·kg−1·K−1) 1.7276 Cold exchanging and supplying Tice-tank/K 272.20–278.15 Tglycol/K 276.62–285.15 Tair,fancoil-in/K 293.45–289.95 Tair,fancoil-out/K 293.45–287.65 mICE/kg 16.98 hICE/(kJ·kg−1) 335 m˙air,fancoil/(kg·s−1) 0.12 tcold/s 1440008:00–17:00 is the working time from Monday to Friday. During working time, the cold demand in the house is very low, so the refrigerator can be driven by PV modules and then the cold can be stored in the ice storage tank. After work, the cold demands of the house increased and the cold stored in the daytime can be charged to service family. So the ISACS driven by DPES has great application prospects in domestic cooling field.The ISACS driven by DPES operated from 08:00 to 23:00. In daytime, solar energy was converted into electrical energy by PV modules. Before 11:00 AM, the irradiation was too low and the electricity generated by PV modules cannot drive ice maker. At the same time, batteries discharged to make up for the lack of electricity produced by PV modules and the ice maker could operate stably and reliably driven by photovoltaic battery hybrid energy supply system. And then, the irradiation gradually increased and the electricity increased along with irradiance. In hybrid energy supply system, the system output power was equal to the power rating of ice maker and was always the same but PV modules output power increased and batteries output power decreased accordingly. About 11 AM, the output power of PV modules was enough to drive ice maker so the output power of batteries is zero. After 11 AM, with the continuous increase of irradiance, not only was the electricity of PV modules used to drive ice maker but also the surplus electricity was stored in batteries. So the batteries were in charge. After noon, the irradiation from sun decreased step by step and so the electricity generated by PV modules reduced gradually. About 14:00, the electricity reducing step by step was only enough to drive ice maker. The charge state of batteries was end. Then, ice maker was driven by PV modules and batteries once again until the ice was enough at 16:00. All of the ice was stored in ice storage tank and the cold can be exchanged by the coil with refrigerating medium ethylene glycol flowing in it at 19:00. At the other end of the coil the cold of refrigerating medium can be blown into the air in the house by fan. So the temperature of house dropped step by step until it tended to balance. About 23:00, the cold of ice in the tank was released completely. And then, the fan and glycol pump would be shut down and the cold supply process was over. Through all of the process, energy is transformed from light to electricity by PV modules and it is stored with ice through phase change latent heat of water. In conversion and transmission process, the energy and exergy of system changed with the change of external environmental conditions. In order to describe system performance clearly and intuitively, the overall performance of the system was evaluated. So the energy efficiency and exergy efficiency of each component of the system in energy conversion and transmission process were calculated and the results are shown in Table4.Table 4
Calculated results of total energy and total exergy of each component of the system in energy conversion and transmission process.
Accepted energy/W Accepted exergy/W Useful power/W Useful exergy/W η/% ψ/% PV module 100007.5 112560.5 15889.79 95910.39 15.89 85.21 Controller 15889.79 95910.39 15254.2 95894.55 96.00 99.98 Batteries 15254.2 95894.55 15162.49 95887.94 99.40 99.99 Inverter 19772.17 95887.94 16222.1 95621.2 82.05 99.72 Distributed photovoltaic energy system 12.44 84.95 Compressor 156936.5 95621.2 151397.8 87286.14 96.47 91.28 Condenser 151397.8 87286.14 Air 51450 86445.17 90.88 99.04 Refrigeration 86136.41 Throttle valve 86136.41 86445.17 86136.41 85784.72 100 99.24 Evaporator 86136.41 85784.72 Output 140714.4 79810.74 336 93.04 Absorbing 54578 Ice 54578 79810.74 9484.22 76271.06 17.38 95.57 Ice maker system 51.20 79.77 Air conditioning system 94842.22 76271.06 81009.64 75747.45 85.42 99.31 ISACS driven by DPES 5.44 67.30Figure4 showed the change curves of energy efficiency and exergy efficiency of ice making and cold-storing process from 08:00 to 16:00 with ice maker driven by DPES. Figure 5 showed the change curves of system energy efficiency and exergy efficiency of fan coil driven by DPES from 19:00 to 23:00.Figure 4
Energy efficiency and exergy efficiency of ice maker driven by DPES.Figure 5
Energy efficiency and exergy efficiency of fan coil driven by DPES.In daytime, ice maker system was driven by DPES batteries and the system energy efficiency and exergy efficiency were 5.44% and 67.30%, respectively. It was found that the total energy efficiency and total exergy efficiency of DPES were 12.44% and 84.95%, respectively, and the total energy efficiency and total exergy efficiency of ice maker system were 51.20% and 79.77%, respectively. After 16:00, the total ice production was 16.98 kg. All of the ice was stored in the tank and air conditioner was open at 19:00 and sustained for about 4 hours with 85.42% total energy efficiency and 99.31% exergy efficiency.It was observed that in ice maker thermodynamic cycle refrigeration efficiency was only 17.38% and the wastage in cooling process was 45093.78 W. In the ice making process, the water can be pumped by a recycle pump and flowed through evaporator to absorb heat. A part of water can be frozen on the evaporator and the rest of water flowed back to the water storage. Thereby, in the water circulating process, the circulating water consumed much of the energy, as shown in Figure6.Figure 6
Ice slide maker (internal picture).In ice maker, compressor was the maximum exergy loss component with 8.78% exergy loss efficiency and it is clear that the evaporator is the most important component of ice maker; the total exergy of evaporator was 79810.74 W and the exergy efficiency was 93.04%. At the same time, the exergy flow from water to the evaporator was 76271.06 W with 95.57% exergy efficiency in ice making process.


## 3.2. Experimental Results and Analysis

### 3.2.1. DPES

The output performance of PV modules was greatly affected by the external environment such as solar irradiance, wind speed, ambient temperature, and PV modules temperature. Therefore, the external environment should be tested through the experiment to calculate and analyze PV output performance. ISACS driven by DPES tested on October 22 in Kunming. Results showed that the irradiation quantityq was 22.17 MJ/m2 from 8:00 to 16:00 and the total amount of ice was 16.98 kg when the fan coil operated for nearly 4 hours. Meanwhile, PV modules output performances were tested and analyzed. Firstly, global solar radiation and ambient temperature were tested, as shown in Figure 7. PV modules temperature is a decisive component of output performance, which can be affected by wind speed changing at any time. Thus, wind speed and PV modules temperature can only be obtained by experimental measurement, as shown in Figure 8.Figure 7
Variations of global solar radiation and ambient temperature on 22nd Nov., 2015, in Kunming.Figure 8
Variations of wind speed and PV modules temperature on 22nd Nov., 2015, in Kunming.The voltage and current variations of PV modules were also tested on October 22 in Kunming. The variation curves are shown in Figure9. Batteries were adopted in testing process to ensure power output stability and store excess power generated by PV modules. The voltage and current of batteries were also measured and variation curves are also shown in Figure 10.Figure 9
Voltage and current variations of PV modules.Figure 10
Voltage and current variations of batteries.It was observed that the maximum tested current was 8.6 A at 12:35. The output voltage changes were 50.8 V–54.2 V and the maximum voltage appears at 12:36. In batteries discharging, ice maker can be driven for 10 hours by batteries only which were in full power state and disconnected with PV modules until the batteries electricity decreased and cannot drive the ice maker. The voltage declined from 54.88 V to 44.6 V step by step and the output current of PV modules was approximately 7 A with small fluctuations between 5.4 A and 7.4 A. In charging process for batteries, batteries disconnect with ice maker and connected with PV modules. Batteries current changed with solar irradiance and voltage increased gradually until batteries were in floating state.Electricity conversed from light by distributed photovoltaic energy system was transferred among batteries, ice maker, and air conditioning system. It was very important to optimize and match the energy among PV modules, storage batteries, and power consumption machine, which can provide an important reference for the future work. Therefore, the experimental study on ISACS driven by DPES was carried out in order to analyze the light-electricity-cold conversion process and the energy transfer characteristics. The results are shown in Figure11.Figure 11
The conversion and flow characteristics of power generated by PV modules among batteries, ice maker, and air conditioning system.There were ice making and ice melting abscission in ice maker operation process. In ice making process, the system current was 1.8 A, the sum of compressor working current and condenser fan current. However, in ice melting abscission process, the system current was 1.4 A, the sum of condenser fan current and solenoid valve operating current. On the other hand, the voltage of AC compressor was 220 V and it was stable all the time. Voltage and current of pump and fan were 220 V and 0.4 A and 380 V and 0.5 A, respectively. In DPES, photovoltaic modules and batteries were untied and used to drive ice maker. Before noon the PV modules output current increases along with solar irradiance and batteries output current gradually decreases. When PV modules output current increased to 6.5 A, batteries output current reduced to 0 A at 10:50 and ice maker can be completely driven by PV modules. Since then, ice maker was driven by PV module and the rest electricity can be stored in batteries. Batteries charging current increased step by step and reached the maximum of 2.0 A until the solar irradiance reached the peak values at 12:20. After noon, batteries charging current gradually reduces along with solar irradiance decreases. It was found that at 15:37 the charging current and PV modules output current were 0 A and 6.5 A. Therefore, ice maker was only driven by PV modules once again. Subsequently, the generated electricity by PV modules was not enough to drive the ice maker and batteries discharged to supplement the shortage. Batteries discharging current increases step by step along with the solar irradiance decreases. Batteries voltage reduced in discharging process and increased in charging process. Batteries separately derived pump and fan coil to supply cold from 19:00 to 23:00. Output current remained constant with 5 A and output voltage slightly decreased gradually. At 23:00, ice melts completely and then all machines were shut down. Until 07:30 of the next day, a new ice making and cold supplying cycle starts. Through the 27-hour experiment, batteries capacity decreased from 65 Ah to 9.76 Ah with a 55.24 Ah decrease. In the first day, batteries discharged 13.7 Ah in the morning and then were charged 7.39 Ah by PV modules. Batteries discharged 12.17 Ah from 15:37 to 19:00. And then ice maker stopped and cold exchanger ran, driven by batteries. The electricity consumption was 20.22 A from 19:00 to 23:00. In the next day, ice maker starts to work at 07:00 and electricity consumption in the morning discharging process and afternoon discharging process was 12.55 Ah and 11.61 Ah, respectively. Batteries were charged 7.62 Ah by PV modules from 10:11 to 15:36.


### 3.2.2. Ice Maker and Cold Storage System

The ice production is shown in Figure12.Figure 12
Ice production of ice maker driven by DPES from 8:00 to 16:00 on 22nd Nov., 2015, in Kunming.It was found that the total amount of ice was 16.98 kg. Ice maker operation cycle time was 10 min and the ice production was about 0.35~0.36 kg every cycle. Ice making efficiency was about 2.12 kg/h. Another experimental test has been conducted at various temperatures of ice maker with every component in thermodynamic cycle of ice making process, as shown in Figure13. The inlet temperature and outlet temperature of evaporator, compressor, and condenser declined stably until ice making process was performed.Figure 13
Temperature variations of ice maker components during two refrigeration thermodynamic cycle processes.In ice making process, ice was pasted well together with five solid walls of grid plate evaporator and cannot be separated off by ice gravity. Thus, compressor must be shut down and a special solenoid valve as shown in Figure14 must be opened which was used for adjusting refrigerant flow and introduced the high temperature and high pressure refrigerant steam expelled from compressor into evaporator.Figure 14
Electromagnetic valve photo adopted in ice maker system.When the electromagnetic was opened, the evaporator plays the role of condenser to release heat to ice through solid walls and ice interfaces begin melting. Ice was divorced from evaporator and fell into tank. Meanwhile, evaporator inlet temperature and outlet temperature increased sharply until the solenoid valve was closed and compressor was turned on and a new ice making process begins. The refrigerant flowed out of evaporator into throttle valve and became low temperature and low pressure liquid. And then, the refrigerant flowed into condenser to absorb heat from outside and condenser inlet and outlet temperature declined sharply when low temperature refrigerant flowed into it. Condenser outlet temperature was higher than condenser inlet because refrigerant absorbed heat from outside through condenser, as shown in Figure12. Refrigerant next flowed into compressor and compressor inlet temperature increased a lot compared to ice making process, which can be changed with inlet and outlet temperature variation of evaporator and condenser. Compressor stops running in ice melting abscission process and the refrigerant gas cannot be compressed in compressor. Therefore, the compressor outlet temperature drops sharply when the solenoid valve opens and it returns smoothly when ice making process opened. Ice melting abscission process was extremely unfavorable to ISACS for a number of reasons as follows:(1)
Ice maker operation period extended for 200 s; furthermore, the ice melting abscission time was one-third time of an ice making cycle.(2)
Service life of compressor will be shortened by frequent start and stop compressor, which has a great impact on ice making process and was extremely unfavorable for energy supply process of DPES. Electricity consumption of solenoid valve was 3.0 × 10−4 kW h in ice melting abscission process and the total electricity consumption was 0.024 kW h from 7:30 to 19:00.


### 3.2.3. Air Conditioning System

The performance of air conditioning system was tested at 19:00. At this time, the temperature of ice stored in ice storage tank and ice storage tank temperature were −3°C and −1°C, respectively. The indoor temperature was 20.5°C. The changes of temperature for exchanging cold and supplying cold processes are shown in Figure15.Figure 15
Temperatures change of exchanging cold and supplying cold process.Ice storage temperature remained constant at −1°C in ice absorbing and melting phase transition processes. However, after 22:20, all ice melted and the tank filled with water at low temperature. Subsequently, the temperature of water gradually increased when the water was employed to cool. Subsequently, it increased from 0°C to nearly 5°C within 40 min. At the end of the experiment, ice storage tank outlet temperature increased to 11.5°C from initial temperature of 3.4°C. In ice phase transition process, ice storage tank inlet and outlet temperatures should be constant values in theory; however, the temperature difference between inlet and outlet increased 7.15°C mainly as a result of the working temperature of pipeline pump. It was found that the fan outlet temperature decreased from 20.5°C to 15.1°C with a 5.4°C decline and the indoor temperature decreased from 20.5°C to 18°C with a 2.5°C decline within 4 hours. Cold power blown out by fan coil was 0.298 W and the cold exchanging and supplying efficiency of the air conditioning system was 90.4%. Furthermore, the ice storage tank has a good thermal insulation performance at 0.10 m thick polyurethane foam. The ice cold loss in ice storage tank was 218.59 kJ as 3.89% of input cold from 19:00 to 23:00 and the cold loss of fan coil was 313.92 kJ as 5.82% of input cold.


## 3.2.1. DPES

The output performance of PV modules was greatly affected by the external environment such as solar irradiance, wind speed, ambient temperature, and PV modules temperature. Therefore, the external environment should be tested through the experiment to calculate and analyze PV output performance. ISACS driven by DPES tested on October 22 in Kunming. Results showed that the irradiation quantityq was 22.17 MJ/m2 from 8:00 to 16:00 and the total amount of ice was 16.98 kg when the fan coil operated for nearly 4 hours. Meanwhile, PV modules output performances were tested and analyzed. Firstly, global solar radiation and ambient temperature were tested, as shown in Figure 7. PV modules temperature is a decisive component of output performance, which can be affected by wind speed changing at any time. Thus, wind speed and PV modules temperature can only be obtained by experimental measurement, as shown in Figure 8.Figure 7
Variations of global solar radiation and ambient temperature on 22nd Nov., 2015, in Kunming.Figure 8
Variations of wind speed and PV modules temperature on 22nd Nov., 2015, in Kunming.The voltage and current variations of PV modules were also tested on October 22 in Kunming. The variation curves are shown in Figure9. Batteries were adopted in testing process to ensure power output stability and store excess power generated by PV modules. The voltage and current of batteries were also measured and variation curves are also shown in Figure 10.Figure 9
Voltage and current variations of PV modules.Figure 10
Voltage and current variations of batteries.It was observed that the maximum tested current was 8.6 A at 12:35. The output voltage changes were 50.8 V–54.2 V and the maximum voltage appears at 12:36. In batteries discharging, ice maker can be driven for 10 hours by batteries only which were in full power state and disconnected with PV modules until the batteries electricity decreased and cannot drive the ice maker. The voltage declined from 54.88 V to 44.6 V step by step and the output current of PV modules was approximately 7 A with small fluctuations between 5.4 A and 7.4 A. In charging process for batteries, batteries disconnect with ice maker and connected with PV modules. Batteries current changed with solar irradiance and voltage increased gradually until batteries were in floating state.Electricity conversed from light by distributed photovoltaic energy system was transferred among batteries, ice maker, and air conditioning system. It was very important to optimize and match the energy among PV modules, storage batteries, and power consumption machine, which can provide an important reference for the future work. Therefore, the experimental study on ISACS driven by DPES was carried out in order to analyze the light-electricity-cold conversion process and the energy transfer characteristics. The results are shown in Figure11.Figure 11
The conversion and flow characteristics of power generated by PV modules among batteries, ice maker, and air conditioning system.There were ice making and ice melting abscission in ice maker operation process. In ice making process, the system current was 1.8 A, the sum of compressor working current and condenser fan current. However, in ice melting abscission process, the system current was 1.4 A, the sum of condenser fan current and solenoid valve operating current. On the other hand, the voltage of AC compressor was 220 V and it was stable all the time. Voltage and current of pump and fan were 220 V and 0.4 A and 380 V and 0.5 A, respectively. In DPES, photovoltaic modules and batteries were untied and used to drive ice maker. Before noon the PV modules output current increases along with solar irradiance and batteries output current gradually decreases. When PV modules output current increased to 6.5 A, batteries output current reduced to 0 A at 10:50 and ice maker can be completely driven by PV modules. Since then, ice maker was driven by PV module and the rest electricity can be stored in batteries. Batteries charging current increased step by step and reached the maximum of 2.0 A until the solar irradiance reached the peak values at 12:20. After noon, batteries charging current gradually reduces along with solar irradiance decreases. It was found that at 15:37 the charging current and PV modules output current were 0 A and 6.5 A. Therefore, ice maker was only driven by PV modules once again. Subsequently, the generated electricity by PV modules was not enough to drive the ice maker and batteries discharged to supplement the shortage. Batteries discharging current increases step by step along with the solar irradiance decreases. Batteries voltage reduced in discharging process and increased in charging process. Batteries separately derived pump and fan coil to supply cold from 19:00 to 23:00. Output current remained constant with 5 A and output voltage slightly decreased gradually. At 23:00, ice melts completely and then all machines were shut down. Until 07:30 of the next day, a new ice making and cold supplying cycle starts. Through the 27-hour experiment, batteries capacity decreased from 65 Ah to 9.76 Ah with a 55.24 Ah decrease. In the first day, batteries discharged 13.7 Ah in the morning and then were charged 7.39 Ah by PV modules. Batteries discharged 12.17 Ah from 15:37 to 19:00. And then ice maker stopped and cold exchanger ran, driven by batteries. The electricity consumption was 20.22 A from 19:00 to 23:00. In the next day, ice maker starts to work at 07:00 and electricity consumption in the morning discharging process and afternoon discharging process was 12.55 Ah and 11.61 Ah, respectively. Batteries were charged 7.62 Ah by PV modules from 10:11 to 15:36.


## 3.2.2. Ice Maker and Cold Storage System

The ice production is shown in Figure12.Figure 12
Ice production of ice maker driven by DPES from 8:00 to 16:00 on 22nd Nov., 2015, in Kunming.It was found that the total amount of ice was 16.98 kg. Ice maker operation cycle time was 10 min and the ice production was about 0.35~0.36 kg every cycle. Ice making efficiency was about 2.12 kg/h. Another experimental test has been conducted at various temperatures of ice maker with every component in thermodynamic cycle of ice making process, as shown in Figure13. The inlet temperature and outlet temperature of evaporator, compressor, and condenser declined stably until ice making process was performed.Figure 13
Temperature variations of ice maker components during two refrigeration thermodynamic cycle processes.In ice making process, ice was pasted well together with five solid walls of grid plate evaporator and cannot be separated off by ice gravity. Thus, compressor must be shut down and a special solenoid valve as shown in Figure14 must be opened which was used for adjusting refrigerant flow and introduced the high temperature and high pressure refrigerant steam expelled from compressor into evaporator.Figure 14
Electromagnetic valve photo adopted in ice maker system.When the electromagnetic was opened, the evaporator plays the role of condenser to release heat to ice through solid walls and ice interfaces begin melting. Ice was divorced from evaporator and fell into tank. Meanwhile, evaporator inlet temperature and outlet temperature increased sharply until the solenoid valve was closed and compressor was turned on and a new ice making process begins. The refrigerant flowed out of evaporator into throttle valve and became low temperature and low pressure liquid. And then, the refrigerant flowed into condenser to absorb heat from outside and condenser inlet and outlet temperature declined sharply when low temperature refrigerant flowed into it. Condenser outlet temperature was higher than condenser inlet because refrigerant absorbed heat from outside through condenser, as shown in Figure12. Refrigerant next flowed into compressor and compressor inlet temperature increased a lot compared to ice making process, which can be changed with inlet and outlet temperature variation of evaporator and condenser. Compressor stops running in ice melting abscission process and the refrigerant gas cannot be compressed in compressor. Therefore, the compressor outlet temperature drops sharply when the solenoid valve opens and it returns smoothly when ice making process opened. Ice melting abscission process was extremely unfavorable to ISACS for a number of reasons as follows:(1)
Ice maker operation period extended for 200 s; furthermore, the ice melting abscission time was one-third time of an ice making cycle.(2)
Service life of compressor will be shortened by frequent start and stop compressor, which has a great impact on ice making process and was extremely unfavorable for energy supply process of DPES. Electricity consumption of solenoid valve was 3.0 × 10−4 kW h in ice melting abscission process and the total electricity consumption was 0.024 kW h from 7:30 to 19:00.


## 3.2.3. Air Conditioning System

The performance of air conditioning system was tested at 19:00. At this time, the temperature of ice stored in ice storage tank and ice storage tank temperature were −3°C and −1°C, respectively. The indoor temperature was 20.5°C. The changes of temperature for exchanging cold and supplying cold processes are shown in Figure15.Figure 15
Temperatures change of exchanging cold and supplying cold process.Ice storage temperature remained constant at −1°C in ice absorbing and melting phase transition processes. However, after 22:20, all ice melted and the tank filled with water at low temperature. Subsequently, the temperature of water gradually increased when the water was employed to cool. Subsequently, it increased from 0°C to nearly 5°C within 40 min. At the end of the experiment, ice storage tank outlet temperature increased to 11.5°C from initial temperature of 3.4°C. In ice phase transition process, ice storage tank inlet and outlet temperatures should be constant values in theory; however, the temperature difference between inlet and outlet increased 7.15°C mainly as a result of the working temperature of pipeline pump. It was found that the fan outlet temperature decreased from 20.5°C to 15.1°C with a 5.4°C decline and the indoor temperature decreased from 20.5°C to 18°C with a 2.5°C decline within 4 hours. Cold power blown out by fan coil was 0.298 W and the cold exchanging and supplying efficiency of the air conditioning system was 90.4%. Furthermore, the ice storage tank has a good thermal insulation performance at 0.10 m thick polyurethane foam. The ice cold loss in ice storage tank was 218.59 kJ as 3.89% of input cold from 19:00 to 23:00 and the cold loss of fan coil was 313.92 kJ as 5.82% of input cold.


## 3.3. Model Validation

The system energy efficiency was 5.44% calculated as shown in Table4. In the experimental test process, the ice produce was 16.98 kg during 08:00~16:00 with the 22.17 MJ/m2 total irradiation. Because the ice maker was driven by PV modules and batteries, the ice produce was not only produced by PV modules. Therefore, the ice production of PV modules must be calculated. The rated power of ice maker was 380 W as shown in Table 1 which was not only made up with the power of compressor as shown in Table 3 but also included the powers of fan, solenoid valve, and water pump. The electricity consumption of ice maker was 3.04 kW·h during 08:00~16:00. But the electricity Q produced by PV modules was only 2.76 kW·h as shown in Table 3. So there was 90.79% power consumed by ice maker coming from PV modules. So the ice produced by PV modules was 15.42 kg. The ice phase change latent heat was 335 kJ/kg. So the energy stored by ice produced by PV modules was 5165.70 kJ. The total area of PV modules was 3.24 m2. So the total solar energy accepted by PV modules was 71830.80 kJ. And the cold exchange efficiency of air conditioning system was 85.42% as shown in Table 4. So the energy efficiency of ISACS driven by DPES tested in the experiment was 6.14%. The relative error of the calculated result was 11.4%. The model can be accepted in engineering application fields.


## 3.4. System Optimization

Analog computation and experimental tests results showed that there were some deficiencies of current system as follows:(1)
Usually the cold demand during the daytime is high, but the system can supply cold through ice melting made by ice maker driven by electricity generated by PV modules in daytime. Therefore, there is no matching between cold supply and demand in the same time.(2)
Energy loss is huge and ice making efficiency is low in ice making process of ice slide maker.Therefore, in order to improve the performance of ISACS driven by DPES and promote the project commercial promotion, some optimization and improvement measurements were proposed as follows:(a)
Evaporator immersion static refrigeration mode was adopted to replace ice harvester refrigeration mode. The optimized coil evaporator was immersed into water to absorb heat and make ice and all the energy was utilized, as shown in Figure16.(b)
The coil cold exchanger was cointegrated with coil evaporator. In refrigeration process, coil cold exchanger has the priority to get cold transferred from the coil evaporator next to the cold exchanger to supply cold for user. Surplus cold is used to make ice to store cold. Consequently, ISACS not only has the out-of-the-box function of ordinary air conditioning, but also effectively improves the appearing phenomenon of overcooling and remedies the disadvantage of cold supply after ice making process in traditional submerged ice making system. The top view of cointegration evaporators and cold exchanger immersed in the ice storage tank is shown in Figure17.Figure 16
Profile of evaporators and cold exchanger immersed in the ice storage tank.Figure 17
Top view of cointegration evaporators and cold exchanger immersed in the ice storage tank.Another immersion static ice maker with evaporator cointegrated with coil cold exchanger was constructed according to the optimization design shown in Figure18.Figure 18
The photo of evaporator immersion static ice maker.The performance of optimized ISACS driven by DPES can be tested and calculated as mentioned above. The system operated on May 12, 2016. The tested and calculated results were shown in Figure19 and Table 5.Table 5
The performance of two kinds of ISACS driven by DPES.
DPES Ice maker system Air conditioning system System Compressor Condenser Throttle valve Evaporator Ice Total η / % ψ / % η / % ψ / % η / % ψ / % η / % ψ / % η / % ψ / % η / % ψ / % η / % ψ / % η / % ψ / % η / % ψ / % 1 12.44 84.95 96.47 91.28 90.88 99.04 100 99.24 336 93.04 17.38 93.04 51.20 79.77 85.42 99.31 5.44 67.30 2 12.44 84.95 96.47 91.28 90.88 99.04 100 99.24 336 93.04 50.12 97.32 147.64 81.24 85.42 99.31 15.69 68.54 Note: 1 stands for the ISACS driven by DPES mentioned above; 2 stands for the optimized ISACS driven by DPES.Figure 19
Temperature changes of ice maker components after optimization.The water tank was full of 13.4 kg water. The system operated at 08:35 and all of the water changed into ice at about 10:40. And then the ice started subcooling and the ice surface became hard at 12:04. At the moment, the surface of evaporator temperature (the temperature of ice core) was −12.11°C and the water temperature was −8.09°C with a 4.20°C temperature difference between ice surface and ice core. At 12:04, the temperatures of evaporator inlet and outlet were, respectively, −20.49°C and −16.51°C and the compressor inlet temperature was −14.28°C. The energy efficiency and exergy efficiency of system components of optimized ISACS driven by DPES were calculated through above formulas and the performances of the system mentioned above and the optimized system were compared in Table5.By structure optimization, the ice making efficiency of evaporator immersion static ice maker was improved to 6.00 kg/h which was 2.88 times of ice slide maker with the same input power. Through theoretical calculation, it was observed that the ice making efficiency and exergy efficiency increased from 17.38% to 50.12% and from 95.57% to 97.32%, respectively. So the ice maker system energy efficiency was improved 2.88 times and exergy efficiency was improved from 76.85% to 78.52%. After the structure optimization, the energy efficiency of ISACS driven by DPES can be improved from 5.44% to 15.69% and the exergy efficiency can be improved from 67.30% to 68.54%. When evaporator cointegrated with coil cold exchanger, ISACS achieved out-of-box functionality as the ordinary vapor compression air conditioning system, which can effectively solve the problem of cold supplying process must lagging behind ice making process in traditional ice storage air conditioning system.


## 4. Conclusions

(1)
The energy utilization ratio and exergy efficiency of ISACS driven by DPES were 5.44% and 67.30%, respectively. DPES conversion efficiency was 12.44% with 84.95% exergy efficiency and ice maker system energy efficiency was 51.20% with 79.77% exergy efficiency. (2)
In order to improve the refrigeration performance of ISACS, evaporator immersion static refrigeration mode was adopted to replace ice harvester refrigeration mode to achieve high efficiency refrigeration. (3)
The coil cold exchanger was cointegrated with coil evaporator. Consequently, ISACS not only has the out-of-the-box function of ordinary air conditioning system, but also effectively improves the appearance of overcooling phenomenon. (4)
By simulation, the system energy utilization efficiency could be improved from 5.44% to 15.69% with exergy efficiency increasing from 67.30% to 68.54%.



---

*Source: 4749278-2016-12-01.xml*