3.1. Description of a Process Line for Waste Heat Recovery
To maintain the high quality of the resulting film at such a high production rate, it is necessary to maintain a constant temperature of the first roll—the calender [
29]. The need to cool this roller and the large amount of waste heat generated here prompted researchers from the Department of Fundamentals of Engineering at the Warsaw University of Life Sciences to design a line to recover this heat and use it in the technological process of stretch film production. The technological line consists of two hydraulic systems separated by a plate heat exchanger (
Figure 7). Heat pump evaporators or fan coolers supply the primary side of the exchanger as the source of process cooling. The load for the secondary side of the heat exchanger (plate exchanger) is the hydraulic installation enabling the reception of the mentioned technological heat from the first calender cylinder.
The even surface temperature of the rotating cylinder is ensured by the appropriate water flow (inlet and outlet temperatures: 15 and 16 °C, respectively) of about 150 m
3/h. During film production, about 176 kW of waste heat (low temperature) is generated, which would be dispersed in the atmospheric air. The innovation of the process is using this heat for drying regranulate thanks to the increase in the temperature level of waste heat. The increase in the temperature level of low-temperature waste heat is achieved using a specially designed heat recovery system based on two compressor heat pumps powered by ecological refrigerants, i.e., agents with low ODP and GWP coefficients. The main technological problem is ensuring a stable and constant inlet water (coolant) temperature in the cylinder. Therefore, two parallel-connected buffer tanks (Z1 and Z2) with a volume of 2000 dm
3 each (simultaneously increasing the water charge in the installation) were used as cold storage. The water pump P1, with smooth regulation of its capacity, enables modification of the water inlet temperature to the calender (three-way mixing valve ZM1 with smooth regulation). The essential source of cooling (the first stage of technological heat recovery) is the compressor heat pump PC1 designed by employees of the Department of Fundamentals of Engineering and Energy at the Institute of Mechanical Engineering at the Warsaw University of Life Sciences (nominal power according to the PNEN 14511 standard (temperatures of the lower and upper heat source are 10 and 35 °C, respectively) is 220 kW), the lower heat source of which is the energy necessary to cover the cooling demand of the calender. PC1 uses semi-hermetic piston compressors, whereby the adjustment of its cooling power (to the current cooling demand of the calender) is achieved thanks to the inverter (regulation of the compressor supply voltage frequency: 30–50 Hz) regulating the compressor capacity. The PC1 evaporator is hydraulically coupled to buffer tanks Z3 and Z4 (2000 dm
3 each) and thanks to the plate heat exchanger (separation between glycol and water on the primary and secondary sides of the heat exchanger, respectively) with Z1 and Z2. Water pump P2 (with smooth regulation of its capacity) causes water flow between the secondary side of the exchanger and Z1 and Z2, while pump P3 (with smooth regulation of its capacity—ensuring a constant temperature gradient: 3K) causes glycol flow between buffer tanks Z3 and Z4 and the primary side of the plate heat exchanger. The pump P4 (with smooth regulation of its capacity: maintaining a constant difference between the inlet and outlet temperature of glycol to the evaporator 3 K, ensuring high efficiency of PC1) causes glycol flow between buffer tanks Z3 and Z4 and the evaporator of PC1. The condenser of PC1, depending on the operating mode (dependent on the current heat demand of the high-temperature technological process of preliminary preparation of regranulate), discharges heat to buffer tanks Z5 and Z6 (2000 dm
3) or fan cooler DR1. The PC1 operating mode is set by changing the position of three-way switching valves ZP1 and ZP2. The first mode of operation of the heat pump PC1 (valves PZ1 and ZP2 in position B) consists of supplying heat to buffer tanks Z5 and Z6 (heat load for the evaporator of the high-temperature heat pump PC2). However, if PC2 is not started or when the temperature in Z5 and Z6 has reached the set point, the following two modes of operation are possible. When the ambient temperature is higher than 3 °C, the medium-temperature heat is dispersed in the environment using the DR1 fan cooler; the buffer tank Z7 increases the glycol charge in the installation (it has a beneficial effect on the length of the heat pump operation cycle). However, when the glycol temperature in Z5 and Z6 reaches the set point, and the ambient temperature is lower than 3 °C, then the ZP1 and ZP2 valves are in position A, and the glycol flow between Z3 and Z4 and DR1 is forced using the P4 pump. This solution allows for significant savings in the electrical energy required to drive PC1 (in this mode, it is switched off). The high-temperature compressor heat pump PC2 enables the generation of high-temperature technological heat. The lower source of PC2 is Z5 and Z6, and the heat load of its condenser is the buffer tank Z8 (high-temperature heat source for the V-type heat exchanger). The heat exchanger is a source of heat used in the preliminary drying process of the regranulate. Details about the technological line for waste heat recovery can be found in our earlier publication [
30].
Based on the technical data from the Project Leader, it was estimated that the demand for technological cooling of the first calender roll would be 176 kW. This value was the basis for the design calculations of individual components of the technological heat recovery line. The demand for technological cooling of the first calender roll with a power of 176 kW and a temperature of 15 °C is the design value of the cooling power of a low-temperature, single-stage heat pump. The heating power and electrical energy required to drive the compressor will be variable and dependent on other operating parameters, including the necessary temperature of the upper source of the first stage of the heat recovery system and the volumetric efficiency of the compressor, which, by the design assumptions, can be regulated using frequency converters.
Analyzing the technical documentation of compressors available on the market in terms of selected refrigerants and the value of the nominal calculated cooling capacity of the first calender roll, it was decided that the most advantageous solution from the technological and economic point of view (taking into account investment and operating costs) would be to design a refrigeration device based on two semi-hermetic Bitzer piston compressors. The use of such a solution meets all the design assumptions. It allows for performing service activities that affect the length of the compressor’s life cycle, including replacing operating fluids (oil) and oil filters. It is also possible to regenerate the compressor if necessary.
3.3. Simulation Studies of the Refrigeration System in the Matlab Simulink Thermolib Package
The simulation studies were aimed at determining the following for the assumed device operating envelope: heating power, cooling power, electric energy consumed by the compressor, COP coefficient, heat transfer rate (in W/K) of the evaporator and condenser necessary for complete evaporation and condensation of the refrigerant and determining the required flow of the working medium being the heat carrier for the assumed design work gradients on the primary side of the 3 K evaporator and the secondary side of the 5 K condenser. The simulation studies began with the design of the refrigeration system according to the adopted single-stage system. At key points of the designed refrigeration system, the mass flow, temperature, pressure and, based on humidity, the state of matter of the refrigerant were examined.
From the design point of view, it was essential to determine the value of the total heat loss coefficient of the evaporator and condenser as a function of the mass flow of the medium on the primary side in the case of the evaporator and on the secondary side in the case of the condenser. In the case of the evaporator, the medium on the primary side was water; in the case of the condenser, the medium on the secondary side was a 40% ethylene glycol solution with a solidification temperature of −35 °C. From the operational point of view, it is essential to determine the mass flow on the primary side of the evaporator and the secondary side of the evaporator at which, for the assumed design temperature difference of the medium (evaporator 3 K, condenser 5 K), the medium will change its state of aggregation—in the case of the evaporator to gas and in the case of the condenser to liquid.
For this reason, PI controllers were designed to regulate the flow on the evaporator’s primary and the condenser’s secondary sides, the input signal of which was the control deviation between the set temperature gradient and the measured value. Simulation studies were started by finding the most unfavorable operating point of the compressor subcritical heat pump; i.e., finding such an operating point defined by the operating envelope at which the heating power and cooling power have the highest values. By determining for this operating point the value of the total heat loss coefficient of the exchanger for which the refrigerant will change its state of matter in the case of the evaporator to gas and in the case of the condenser to liquid, for the assumed temperature gradient, it is possible to calculate the required minimum flow number of the evaporator and condenser. For all other operating points located in the device’s operating envelope, the selected heat exchange surface of the evaporator and condenser will be greater than the minimum required, as a result of which the refrigerant will continuously be evaporated in the case of the evaporator and always condensed in the case of the condenser. In the case of the 1234yf agent, the most unfavorable operating point (operating point characterized by the highest cooling power and heating power) turned out to be the operating point for the highest saturation temperature of the refrigerant and the lowest condensation temperature. For the agent to change its state of matter to gas, the value of the total heat transfer coefficient of the evaporator was 132,000 W/K, while in the case of the condenser, the value of this coefficient was 84,000 W/K.
Figure 8 shows the cooling capacity of a low-temperature heat pump for a wide range of condensation (from 25 to 50 °C) and saturation (from −6 to 10 °C) temperatures. It can be seen that the cooling capacity increases with increasing saturation temperatures and decreasing concentration temperatures. The available cooling power corresponds to the technological cold for cooling the first calender roll. The nominal operating point is 4 °C and 35 °C (saturation and condensation temperature, respectively). The deviation of the parameters visible in
Figure 8 and
Figure 9 (2 °C and 30 °C) results from the isentropic efficiency of the compressor (the manufacturer’s catalogue data were used in the calculations).
The dependence of heating power on the refrigerant’s saturation temperature is presented similarly (
Figure 9). The device generates the highest heating power value of 405 kW for the refrigerant’s highest design saturation temperature (10 °C) and the lowest condensation temperature (25 °C). For this operating point, the power of both compressors is 39 kW. In the case of an increase in the condensation temperature from 25 °C to 50 °C at a constant saturation temperature of the refrigerant of 10 °C, the heating power received from the condenser will decrease from 405 to 354 kW, and the compressor power will increase to 87.2 kW. The lowest value of heating power is obtained for the saturation temperature of the refrigerant of −6 °C and its condensation temperature of 50 °C. The heating power value for this operating point is 183.8 kW with a demand for compressor power of 38.8 kW.
A specific trend can be observed when analyzing the power variation in compressors controlled at 100% as a function of the refrigerant saturation temperature. In the entire range of the refrigerant saturation temperature for the lowest value of the refrigerant saturation temperature, the total power of both compressors operating at maximum volumetric efficiency is practically constant at 40 kW. The higher the refrigerant condensation temperature, the higher the electricity consumption, and the higher the refrigerant saturation temperature (
Figure 10). At the highest tested condensation (50 °C) and saturation (10 °C) temperatures of the refrigerant, the demand for compressor power is almost 90 kW.
As a result (
Figure 11), the most significant COP variability over the entire range of refrigerant saturation temperature is achieved by the designed device at the lowest condensation temperature constant over time, which is 25 °C. For a refrigerant saturation temperature of 10 °C and a condensation temperature of 25 °C, the maximum COP value of 10.37 is achieved. The minor COP variability over the entire range of refrigerant saturation temperature occurs for the highest condensation temperature of 50 °C: the average 3.8. For the nominal operating point (4 °C and 35 °C, saturation temperature and condensation temperature, respectively), the COP is 5.5.
The COP value is determined by the temperature of the lower and upper heat sources, which determines the saturation and condensation temperature of the refrigerant. As was mentioned, in the case of a low-temperature heat pump, the COP is determined by the required temperature of the technological medium cooling the first calender roll. The nominal temperature of the medium is 16 °C, and the heat exchanger was designed for the difference between the temperature of the medium and the refrigerant saturation temperature of 10 K. Therefore, the nominal saturation temperature is 6 °C. The temperature of the upper heat source depends on the required temperature of the technological process: whether the heat is dissipated in the atmosphere or used for drying. The computer simulations shown in
Figure 8,
Figure 9,
Figure 10 and
Figure 11 concern a wide range of temperatures (condensation: from 25 to 50 °C and saturation from −6 to 10 °C) for which COP values were obtained (even COP > 10).
Analyzing the quality of regulation obtained using the designed PI controllers, it is visible that in the case of the evaporator, the maximum regulation deviation is 0.35 °C. The smallest is close to 0 (
Figure 12). In the condenser case, the maximum regulation deviation is 0.35 °C, and the minimum equals 0 (
Figure 13). The reason for the deviation of the parameters in
Figure 12 and
Figure 13 results from the mass flow of the working medium. The temperature delta was not achieved exactly as 3.0 and 5.0 (3.12, 3.15, 3.25, 3.37 and 4.83, 5.11, 5.37, respectively).
The simulation studies allowed the calculation of the required flow of the working medium on the evaporator’s primary and the condenser’s secondary sides. For both the evaporator and the condenser, the required flow of the medium increases linearly for a constant saturation temperature of the refrigerant and an increase in the lower source temperature. In the case of the evaporator, the most enormous volumetric flow is 2217.5 dm
3/min, and in the case of the condenser, the maximum flow value is 1582.1 dm
3/min. The calculated values of the mass flow through the evaporator and condenser allowed the selection of circulation pumps, forcing the working medium’s flow through the evaporator’s primary and the condenser’s secondary sides (
Figure 14 and
Figure 15).
The fact that the piston compressors can regulate volumetric efficiency in the range of 30–50 Hz allows the device’s efficiency to be adjusted to the current demand for technological cooling.
3.4. Selection of Structural Elements and Dimensioning of Pipes for the Refrigeration System of a Low-Temperature Compressor Heat Pump
Based on the simulation studies, a refrigeration system for a low-temperature subcritical compressor heat pump was developed (
Figure 16).
The discharge line is marked in red, the suction line is in liquid, and the liquid line is in yellow. The oil line is marked in brown. The refrigeration system is based on two semi-hermetic piston compressors connected in parallel (tandem). In the event of the need to perform service activities, a rota-lock valve is installed on each compressor’s suction and discharge nozzles. Since the system can be operated depending on the demand for technological cooling with one or two compressors switched on, a non-return valve is installed on the discharge nozzle of each compressor. A high and low-pressure transducer is installed in the discharge and suction lines, respectively, providing information on the current value of the refrigerant pressure using an analogue voltage signal in the range of 0–10 V DC. High and low-pressure analogue transducers also constitute the first level of protection of the refrigeration system against too-high and too-low refrigerant pressure. In addition, the second and third levels of security in mechanical pressure switches were used to protect against too-low and too-high pressures. The second level protects individual compressors, while the third level protects both. Plate heat exchangers were used as the evaporator and condenser. A filter drier was installed in the liquid line to protect the motor expansion valve, and a sight glass was used to assess the refrigerant’s state. Due to the possibility of operating one or two compressors and due to the possibility of regulating their capacity, a liquid refrigerant reservoir was installed in the fluid line. It is a refrigerant storage facility that adjusts the appropriate amount of refrigerant depending on the temperature of the lower and upper source and the volumetric capacity of the compressors. A safety valve was used to protect the system against high pressure. The safety valve was connected in parallel to the evaporator by connecting the inlet nozzle of the valve with a pipe to the high-pressure liquid line and the outlet nozzle with a pipe to the suction line. It should be noted that it was connected between two siphons that also function as a liquid separator. This solution protects the compressor from being flooded with a liquid agent in the event of low superheating of the agent and a change in the evaporator’s cooling capacity. The compressor oil pans were connected with an oil line, in which a sight glass was installed. The oil line is designed to equalize the oil level in the oil pans if only one compressor is switched on. This solution ensures correct oil distribution in the device’s refrigeration system. To protect the suction line from the condensation of water vapor contained in the air and, as a result, from frosting, a 19 mm thick rubber insulation was used. The diameters of individual sections of the refrigeration system were selected based on Equation (2).
The design problem related to the selection of diameters is that depending on the temperature of the lower and upper source of the specified work envelope, the volume flow of the refrigerant is variable due to the change in the density of the refrigerant. The diameter of the discharge line was estimated assuming that the refrigerant flow rate would be variable within the range of 8–20 m/s, the high-pressure liquid line within the range of 0.5–1.5 m/s, the low-pressure liquid line within the range of 0.4–0.5 m/s, and the refrigerant flow rate in the suction line would be greater than 8 m/s. The selected diameters of the individual sections of the refrigeration system were marked on the diagram (
Figure 16).
SolidWorks software was used to model the low-temperature compressor heat pump. The construction is based on a steel frame made of rectangular steel bars. The frame is placed on legs made of round steel bars (
Figure 17).
To reduce vibrations caused by the compressors, the cooling system is placed on springs where guides are mounted. The designed device is 2270 mm long, 1127 mm wide and 1350 mm high. The actual dimensions of the device may differ from those of the design because the width and length of the device primarily result from the technology used by the manufacturer to make bends, siphons, and counter-siphons. In the case of using ready-made bends or elbows, regardless of the connection technology, the actual dimension should be consistent with the design. In the case of using pipe bending technology, the dimension of the device will be larger than the design one. The front of the device is the control cabinet, which contains the electrical installation with protections, the controller and the operator panel.
3.5. High-Temperature Heat Pump
An analogous procedure was followed for the high-temperature subcritical compressor heat pump.
Knowing the volumetric capacity of the compressor of 156 m
3/h, and for the assumed temperatures, the value of the specific volume of steam at the suction nozzle (p-h diagram of the R1234yf refrigerant), the mass flow of the refrigerant was calculated using Equation (1) (
Table 2).
Based on the mass flow of the refrigerant for the assumed design values of the refrigerant saturation and condensation temperature, similar to the low-temperature heat pump, thermodynamic calculations of the refrigeration system were performed in the Matlab & Simulink Thermolib Basic package to determine, among others, heating power, cooling power, compressor power, and COP coefficient.
Figure 18,
Figure 19,
Figure 20,
Figure 21,
Figure 22,
Figure 23,
Figure 24 and
Figure 25 show the operating parameters of both compressors’ high-temperature heat pump at 100% control (supply voltage frequency 50 Hz).
Figure 18 shows a graph of the cooling power of the designed low-temperature heat pump for the saturation temperature of the refrigerant in the range of 1–30 °C, which corresponds to the lower source temperature in the range of 26–40 °C, for refrigerant condensation temperatures in the range of 75–85 °C. The graph shows that the highest value of cooling power necessary to evaporate the refrigerant (the most unfavorable operating point of the device from the design point) occurs for the highest saturation temperature of the refrigerant of 30 °C and the lowest condensation temperature of the refrigerant of 75 °C, and the value of cooling power necessary to evaporate the refrigerant is 218 kW. As can be seen from the graph, the cooling power of the designed heat pump, even with 100% control of both compressors, can be adjusted to the current demand for technological cooling by regulating the value of the condensation temperature of the refrigerant. For example, at a constant saturation temperature of 30 °C, increasing the value of the condensation temperature of the refrigerant from 75 °C to 85 °C, the cooling power decreases from 218 to 190.4 kW. The lowest value of cooling power is obtained for the refrigerant condensation temperature of 85 °C at its saturation temperature of 16 °C. The cooling power value for this operating point is 130.2 kW.
The dependence of heating power on the refrigerant’s saturation temperature is presented similarly (
Figure 19). The device generates the highest heating power value of 315.9 kW for the refrigerant’s highest design saturation temperature (30 °C) and the lowest condensation temperature (75 °C). For this operating point, the power of both compressors is 85 kW. In the case of an increase in the condensation temperature from 75 °C to 85 °C at a constant saturation temperature of the refrigerant of 30 °C, the heating power received from the condenser will decrease from 315.9 kW to 282.9 kW, and the power of the compressors will increase to 94.2 kW. The lowest value of heating power is obtained for the saturation temperature of the refrigerant of 16 °C and its condensation temperature of 85 °C. The heating power value for this operating point is 188.1 kW with a compressor power demand of 59.8 kW.
A particular trend can be observed by analyzing the variation in the compressor power controlled at 100% as a function of the refrigerant saturation temperature. In the entire range of the refrigerant saturation temperature, the compressor power increases regardless of the upper source temperature, while in the saturation temperature range of 16–24 °C, it increases more steeply than in the saturation temperature graph of 24–30 °C (
Figure 20).
As a result, the variation in the COP coefficient over the entire range of the refrigerant saturation temperature of the designed device for the design temperatures of the upper source is similar. The maximum value of the COP coefficient of 3.71 is achieved for the refrigerant saturation temperature of 30 °C and the condensation temperature of 75 °C. The lowest value of the COP coefficient (2.78) occurs for the highest condensation temperature of 85 °C and the lowest saturation temperature of 16 °C (
Figure 21).
Analyzing the quality of regulation obtained using the designed PI controllers, it is visible that in the case of the evaporator, the regulation deviation is maximally 0.19 °C, and the smallest is close to 0 (
Figure 22). In the case of the condenser; the regulation deviation is at the level of ±0.003 °C (
Figure 23). The reason for the deviation of the parameters in
Figure 22 and
Figure 23 results from the mass flow of the working medium. The temperature delta was not achieved exactly as 3.0 and 5.0 (3.05, 3.19 and 5.002, respectively).
The simulation studies allowed the calculation of the required flow of the working medium on the evaporator’s primary and the condenser’s secondary sides. For both the evaporator and the condenser, the required flow of the medium increases linearly for a constant saturation temperature of the refrigerant and an increase in the lower source temperature. In the case of the evaporator, the enormous volumetric flow is 1300 dm
3/min (
Figure 24). In the case of the condenser, the maximum flow value is 1130.1 dm
3/min (
Figure 25). The calculated values of the mass flow through the evaporator and condenser allowed the selection of circulation pumps forcing the flow of the working medium through the primary side of the evaporator and the secondary side of the condenser.
The piston compressors can regulate volumetric efficiency in 30–50 Hz. It allows the device’s efficiency to meet the current demand for technological cooling.