A formation of the heat pump mathematical models A Y Marinchenko1 and A V Edelev1 1 Melentiev Energy Systems Institute of Siberian Branch of SB RAS, Lermontov St. 130, Irkutsk, Russia, 664033 E-mail: marinchenko@isem.irk.ru, flower@isem.irk.ru Abstract. The paper addresses the mathematical modeling heat pump main elements and heat pump based heat supply systems. It is used in the calculating expected heat supply costs on the basis of a heat supply system functioning simulation. Using these models, the current state of the heat pump NT-60-1 located in the Baikal Museum of the Scientific Center of the Siberian Branch of the Russian Academy of Sciences in Listvyanka was analyzed. Recommendations on improving the heat pump efficiency are given. 1. Introduction The existing Russian district heat supply systems use the thermal power plants that produce combined heat and electricity as heat sources. The boiler houses of various types are heat sources in other heat supply systems. One of the ways to increase the energy efficiency of district heating systems is to use a heat pump (HP) technology, which pick over heat from surface or underground water, industrial or domestic waste water, etc., and transfer it to the consumer with consuming a certain quantity of electricity. The fuel consumption to produce the same heat amount by the HP is lower than by the boiler house. Installing heat pump units in the existing boiler houses burning coal and fuel oil and located within areas with increased environmental requirements (Lake Baikal zone, etc.) will greatly reduce the fuel consumption. It results in the significant improvement of the environmental situation in a region. The most promising load for the HP is heating and hot water production. The heat is released in the form of hot water with the relatively low temperature (according to the temperature plan 95/70°C or lower). Under these conditions the HP provides significant fuel savings. However, the heat load depends on the outdoor temperature and varies significantly during day and according to the seasons of the year. Possible random deviations of various kinds lead to the necessity of using variable HP operation modes, redundancy and covering peak loads by other usually less capital-intensive heat sources. Significant opportunities to improve the efficiency of a HP based heat supply system are also available due to the optimization of equipment size and its operation. This work addresses the formation of mathematical models of the HP main elements and heat supply systems based on HP to use them in the calculating the expected heat supply system costs on the basis of the heat supply system functioning simulation. The rest of the paper is structured as follows. Section 2 and 3 give a brief overview of the use of the HP in heat supply systems and related work, respectively. Section 4 presents a framework for mathematical modeling HP main elements. Section 5 describes testing proposed mathematical models. Section 6 concludes the paper. Copyright © 2020 for this paper by its authors. Use permitted under Creative Commons License Attribution 4.0 International (CC BY 4.0). 2. The use of heat pumps in heat supply systems The HP is a device that uses the phase transitions of the working fluid (evaporation and condensation) to pick up energy from a low potential source (evaporation) and transfer it to a consumer with high temperature (condensation). The theoretical foundation of this process was presented in 1824 by S. Carnot. And in 1852 W. Thompson (Lord Kelvin) created the first working HP. A diagram of a heat supply system based on the HP is shown in Figure 1. 8 6 9 7 4 5 3 2 1 Figure 1. A diagram of a heat supply system based on the HP: 1 – pump; 2 – evaporator; 3 – compressor with electric drive; 4 – condenser; 5 – thermostatic valve; 6 – electric boiler; 7 – regulating water valves; 8 – heating devices; 9 – circulation pump of a heat supply system The principles of operation of the heat supply system in Figure 1 are as follows: 1. In heat exchanger 2 or evaporator thermal energy from an external heat source is transferred to the HP working fluid known as a refrigerant. 2. The refrigerant evaporates and moves to the compressor 3, which compresses it to a certain value. The refrigerant temperature, pressure, and enthalpy increase. 3. The compressed refrigerant moves through a heat exchanger 4 (condenser), in which the refrigerant condenses and transfers heat to water in a consumer’s side system. 4. The refrigerant enters a throttle valve 5, which reduces the refrigerant pressure, and then it moves to the evaporator 2. 5. After the condenser 4 water in a consumer’s side system flows to heating devices 8. If it is necessary, the electric boiler 6 heats water up to the required temperature. As noted above, it is possible to increase the efficiency of the HP based heat supply systems optimizing their design and operation. This optimization can be carried out on the basis of the methodology for calculating the expected heating costs based on the heat supply system functioning simulation. 3. Related work Significant economic and environmental advantages make the HP a promising technology in the field of heat supply for most regions and countries of the world. A number of various heat pumps types units are operating in many countries now [1-4]. Their number is increasing every year. The total number of installed HP units in the world exceeds 100 million. Many countries provide substantial support to households that use heat pumps through national energy conservation programs. It should be noted that the majority (about 68%) of new HP units have the capacity of up to 15 kW that corresponds to the heat requirements of a separate house with the heated living area of 140 m2. At the same time, there are quite big heat pumps units. Existing centralized district heat supply systems which use boiler houses as a heat source can be upgraded with combined installations, including the HP. A significant amount of work has been devoted to the methods of mathematical modeling and optimization of various HP types and heat-producing plants based on them [1-7]. Usually mathematical modeling is used to study of thermodynamic cycles [3, 8-11], determine the best working fluids for the HP [14-17], increase the thermodynamic [5, 7-10, 16-19] and economic efficiency of HP units [6, 7, 11, 19-22]. The specificity of our study is that the developed HP models will be used to investigate the application of environmental technologies in the Baikal natural territory infrastructure. 4. Mathematical models of the HP elements One of the main HP elements is a compressor, in which the process of compression of the working fluid occurs. In the design calculation of the compressor, the input flow rate, pressure and enthalpy of the working fluid, and its output pressure are known. The output enthalpy of the working fluid and maximum mechanical capacity of the compressor are calculated. The mathematical representation of the compression process of the working fluid is described by the following dependencies:  S 1  f spi P1 , H 1 ,  (1) H ad A  fips P2 , S1  , (2) H 2  H 1  H 1  H ad  ots , A (3)  N  G  H 1  H 2  m ,  (4) where S is entropy; N is capacity; ots is internal relative efficiency of the compressor;  m is mechanical efficiency of the compressor; H ad is enthalpy of the working fluid at the end of the ideal (adiabatic) process of compression to the outlet pressure, Р is the pressure; G is the flow rate; Н is enthalpy of the working fluid. Hereinafter, superscript 1 marks input parameters, and 2 marks output parameters. The HP condenser is to transfer heat from the working fluid to water in the consumer’s side system. In this case, the working fluid converts from a vaporous state to a liquid state, i.e., it condenses. In the condenser a shell and tube construction is most often used. In this case, the heated water is passed inside the pipes, and the working fluid is condensed on the outside of pipes in the inter-tube space. The system of equations describing the operation of the HP condenser includes the following dependencies:  heat balance equation  Gw  H w2  H w 1    ut  G1p  H 1p  H p2 ,   (5)  equation of material balance Gр2  G1p , (6)  heat transfer equation F  Gw  H w2  H w1   s , (7) Tlog  K where G is flow rate, H is enthalpy, ut is heat exchanger heat utilization coefficient, F is heating surface area, Tlog is logarithmic temperature head,  s is heat exchanger safety factor on the heating surface, K is heat transfer coefficient. Hereinafter, subscripts “p” and “w” denote the parameters related to the HP working fluid and water, respectively. Dependencies for determining the logarithmic temperature head are Tsat  f tpн Pp ,   (8) T  Tsat  Tk , 1 w (9)  w  T  f tpi P , H , 2 2 w 2 w (10) T  T  T   T  T  , sat 1 w sat w 2 (11)  T T  log 1 ln  sat  w  Tsat  T  2 w where Tsat is saturation temperature at an average pressure of the working fluid in the heat exchanger; Tk is end temperature pressure (at the outlet of water from the heater). The relationship between the enthalpy of water and temperature is  H w2  f ipt Pw2 ,Tw2 . (12) Dependence for determining the heat transfer coefficient is 1 K , (13) 1 D1 D1  D2   D1 1     2 D2 D1  D2   m 1 where D1 and D2 are the outer and inner diameters of the pipes, respectively, 1 is the heat transfer coefficient from the working fluid to the pipe,  2 is the heat transfer coefficient from the pipe to water, m is the thermal conductivity of the pipe metal. The heat transfer coefficient 1 is determined by the procedure described in [23] for the liquid condensation. Dependencies for determining the hydraulic resistance of a condenser along the path of the working fluid and water are P P  2     W  , 1 t tp t m cp 2 w (14) 2  g Vwcp w w P  P  0,33  2 1     W  , n tp n m 1 2 p (15) 2  g V p1 p p where g is the acceleration of gravity, W is the speed of water or steam, and V is the specific volume. The coefficients of hydraulic friction and local losses  tp ,  mn ,  tp ,  mt for water inside the pipes n t and for washing the pipe bundle are calculated in accordance with the standard calculation method: Gw1  Vwch Z tp    , (16) 4    Wwch  D2 2 F Ltp  , (17)   Z tn  D1 F  W p1  S1  D1  Z row  . (18)   D1  G1p  Vh1 where Z tp is number of pipes connected in parallel; Ltp is average length of one pipe; Z row is number of rows of pipes along the steam flow; S1 is transverse (along the working body) pipe pitch. The initial data for the modeling condenser are the flow rate, enthalpy and pressure of the heated water at the inlet, the pressure and enthalpy of the working fluid, the speeds of the working fluid and water, the inner and outer diameters of the pipes, and their transverse steps. The flow rate of the working fluid, the enthalpy and pressure of the water at the outlet, the enthalpy and pressure of the working fluid at the outlet, the surface area of the heating, the number of pipes, the number of rows of pipes along the working fluid, etc., are determined. The HP evaporator is to take heat from a low-temperature heat source at low temperature and transfer it to the working fluid; while the working fluid converts from the liquid phase to the gaseous one. The system of equations describing the heat transfer process in the evaporator is:  heat balance equation Gw  H w 1    H w2  ut  G1p  H p2  H 1p    , (19)  material balance equation Gр2  G1p , (20)  heat transfer equation F  Gw  H w1  H w2  з .  (21) Tlog  K Dependencies for determining the logarithmic temperature pressure are Tsat  f tpн Pp ,   (22) T  Tsat  Tk , 1 w (23)  w  T  f tpi P , H , 2 w 2 2 w (24) T  T  T   T  T  . 1 w sat 2 w sat (25)  T  Tsat  log 1 ln  w   T  Tsat  2 w Dependencies for determining the heat transfer coefficient, the hydraulic resistance of the evaporator along the path of the working fluid and water, its design and layout characteristics are similar to the dependences (14)-(18) used in the modeling the condenser. The heat transfer coefficient  1 is determined by the method described in [23] for the bubble boiling process. The initial data for the modeling evaporator are the enthalpy and pressure of the heated water at the inlet, the flow rate, the pressure and enthalpy of the working fluid, the speeds of the working fluid and water, the inner and outer diameters of the pipes, and their transverse steps. The water flow rate, enthalpy and pressure of the water at the outlet, the enthalpy of the working fluid at the outlet, the surface area of the heating, the number of pipes, the number of rows of pipes along the working fluid, etc. are determined. To calculate the thermodynamic and transport properties of the HP working fluid there is an approach that interpolates dependencies between given base (tabular) points at which the properties are known. The data for Freon R134a (http://www.coolprop.org/fluid_properties/fluids/index.html) were used as basic table values. 5. Testing mathematical models The developed mathematical model was used to analyse the current state of the HP NT-60-1 located in the Baikal Museum of the Scientific Center of the Siberian Branch of the Russian Academy of Sciences in Listvyanka and to develop recommendations for improving its efficiency. The source of low-grade heat is the running water of the pools with a discharge temperature depending on the season 4 - 9 ºС and a flow rate of about 20 m³/h. A general view of the HP NT-60-1 is shown in Figure 2. The scheme of the heat supply system of the Baikal Museum based on the NT-60-1 vapor compression HP is presented in Figure 3. Figure 2. The HP NT-60-1 for the heat supply of the Baikal Museum of the Scientific Center of the Siberian Branch of the Russian Academy of Sciences on the basis of renewable heat of Lake Baikal water The ozone-safe freon R-134a was used as the NT-60-1 working fluid, the HP scheme was developed taking into account the specifics of working with low-temperature renewable sources and provides higher conversion factors compared to traditional layout schemes. For this purpose an evaporator with an in-line boiling of the refrigerant was used, it is adapted for working conditions in cold Baikal water at low temperatures of the order of 2-4°C, and schemes for joint operation of the HP with other heat sources, including electric boilers, were developed. In addition, while working on cold heat sources in order to increase the specific HP heat output, it is planned to heat the low-grade heat carrier (Baikal water) in front of the evaporator due to the deep cooling of the refrigerant condensate in the subcooler. In this case, all or part of water passes through the supercooler 8. In the absence of the need for hot water supply, this operation scheme provides an increase in the overall efficiency of the vapor compression HP by minimizing irreversible losses in the throttle device. The estimated average annual HP coefficient of performance, which determines the ratio of heat generated to energy consumption of the compressor drive, was about 3.5. During the heating period the HP provides heating of water up to 60 ºС on the working fluid R-134a, which corresponds to the HP design specifications. The HP was also worked together with the electric boiler, previously used in the museum's heating center for heating. The implementation of the scheme of their joint application for heat supply where the HP works on the base load and of electric works only on peak load, reduces the total heating energy consumption by 2–3 times. Figure 3. The schematic diagram of the Baikal Museum’s heat supply system: 1 – pumping station; 2 – Museum’s aquariums; 3 – pumping installation; 4 – regulating water valves; 5 – evaporator with internal coolant boiling; 6 – reciprocating compressor with electric drive; 7 –condenser; 8 – supercooler; 9 – temperature-regulating valve; 10 – electric boiler; 11 – heating devices; 12 – circulation pump of the heating system Table 1 presents the nominal parameters of the HP NT-60-1 obtained from the operation manual. Table 2 presents the calculated HP parameters obtained using the mathematical model of the HP at the initial nominal parameters. As one can see from Table 2, the power consumption of the compressor of the HP is 24.9 kW at nominal temperatures at the outlet of the condenser and evaporator. Table 1. The passport parameters of the HP NT-60-1 Dimension Value Maximum power consumption kW 30.1 Nominal flow rate per condenser m3/h 5 3 Nominal evaporator water consumption m /h 27.8 Nominal water temperature at the condenser outlet 0 С 62±1 Nominal water temperature at the evaporator outlet 0 С 2.5±0.5 Coefficient of performance 3.5 Table 2. The HP NT-60-1 parameters calculated at nominal values Dimension Value Power consumption kW 24.9 Condenser water flow m3/h 5 Evaporator water consumption m3/h 27.8 Condenser inlet water temperature 0 С 30 Condenser outlet water temperature 0 С 60.9 Evaporator inlet water temperature 0 С 7 Evaporator outlet water temperature 0 С 2.48 Coefficient of performance 3.4 Table 3 shows the existing state of the NT-60-1. It should be noted that the water flow in the heat supply system was unknown, so in this case it was calculated according to the heat balance equation, taking into account losses in the HP. Also a fairly low value of the coefficient of performance of the HP indicates that its efficiency is rather low. Table 3. The real HP NT-60-1 parameters Dimension Value Power consumption kW 26 Condenser water flow m3/h 7 Evaporator water consumption m3/h 16 Condenser inlet water temperature 0 С 49 Condenser outlet water temperature 0 С 59 Evaporator inlet water temperature 0 С 6.4 Evaporator outlet water temperature 0 С 3.2 Coefficient of performance 1.54 Table 4 presents the NT-60-1 parameters calculated under the same conditions as in Table 3, except for the return water temperature, which was adopted at 42°C (7°C lower than before). As one can see from Table 4, it is allowed for almost the same electricity consumption value to double the efficiency of the NT-60-1 (coefficient of performance is 2.7) Table 4. The recommended HP NT-60-1 parameters Dimension Value Power consumption kW 25.6 Condenser water flow m3/h 7 Evaporator water consumption m3/h 16 Condenser inlet water temperature 0 С 42 Condenser outlet water temperature 0 С 59.9 Evaporator inlet water temperature 0 С 6.4 Evaporator outlet water temperature 0 С 1.8 Coefficient of performance 2.7 The main recommendation drown from the results of the experimental analysis is to increase the flow rate of low-grade water through the evaporator to nominal values. This will bring the parameters of NT-60-1 to ones given in the passport and significantly increase the efficiency of the HP. A decrease in the temperature at the inlet of the condenser to 30-35°C and at the output of the condenser to 500С will also have a positive effect, which will lead to a decrease in the compressor energy consumption without reducing the thermal efficiency of the HP. 6. Conclusions In the paper, we focused on the problem of the HP use study. It is a relevant problem. This is due to significant economic and environmental advantages the HP use. To this end, we have proposed the specialized mathematical models for modeling and studying the operation efficiency of HP under its practical use. Based on the experimental analysis, we have formed useful recommendations related to improving the functioning of this HP. Our future work related to applying the proposed models in calculating the expected heat supply system costs. Acknowledgments The study is supported by the Russian Foundation of Basic Research and Government of Irkutsk Region, project no. 20-47-380002. We sincerely thank the director, scientific leader, management, and employees of the Baikal Museum for advising and expert supporting in our study. References [1] Bonin J. 2015 Heat pump planning handbook (Routledge, New York) [2] Fischer D and Madani H 2017 On heat pumps in smart grids: A review Renewable and Sustainable Energy Reviews 70 342–357 [3] Arpagaus C, Bless F, Schiffmann J. and Bertsch S S 2016 Multi-temperature heat pumps: A literature review Int J Refrig 69 437–465. [4] Rees S ed. 2016 Advances in ground-source heat pump systems (Woodhead Publishing) [5] Wallerand A S, Kermani M, Kantor I and Maréchal F 2018 Optimal heat pump integration in industrial processes Applied Energy 219 68–92 [6] Averfalk H, Ingvarsson P, Persson U, Gong M and Werner S 2017 Large heat pumps in Swedish district heating systems Renewable and Sustainable Energy Reviews 79 1275–1284 [7] Simonetti R, Moretti L, Molinaroli L and Manzolini G 2020 Energetic and economic optimization of the yearly performance of three different solar assisted heat pump systems using a mixed integer linear programming algorithm Energy Convers. Manag. 206 112446 [8] Mota-Babiloni A, Mateu-Royo C, Navarro-Esbrí J, Molés F, Amat-Albuixech M and Barragán- Cervera Á 2018 Optimisation of high-temperature heat pump cascades with internal heat exchangers using refrigerants with low global warming potential Energy 165 1248–1258 [9] Culha O, Gunerhan H, Biyik E, Ekren O and Hepbasli A 2015 Heat exchanger applications in wastewater source heat pumps for buildings: A key review Energy Build. 104 215–232 [10] Zhang Q, Yang S, Shi P, Hou W, Zeng A, Ma Y and Yuan X 2020 Economically and thermodynamically efficient heat pump-assisted side-stream pressure-swing distillation arrangement for separating a maximum-boiling azeotrope Appl. Therm. Eng. 173 115228 [11] Blázquez C S, Verda V, Nieto I M, Martín A F and González-Aguilera D 2020 Analysis and optimization of the design parameters of a district groundwater heat pump system in Turin, Italy Renew. Energ. 149 374–383 [12] Bamigbetan O, Eikevik T M, Neksa P and Bantle M 2017 Review of vapour compression heat pumps for high temperature heating using natural working fluids Int. J. Refrig. 80 197–211 [13] Zühlsdorf B, Jensen J K and Elmegaard B 2019 Heat pump working fluid selection—economic and thermodynamic comparison of criteria and boundary conditions Int. J. Refrig. 98 500–513 [14] Uusitalo A, Turunen-Saaresti T, Honkatukia J, Tiainen J and Jaatinen-Värri A 2020 Numerical analysis of working fluids for large scale centrifugal compressor driven cascade heat pumps upgrading waste heat Applied Energy 269 115056 [15] Brodal E and Eiksund O 2020 Optimization study of heat pumps using refrigerant blends– Ejector versus expansion valve systems Int. J. Refrig. 111 136–146 [16] Willem H, Lin Y and Lekov A 2017 Review of energy efficiency and system performance of residential heat pump water heaters Energy Build. 143 191–201 [17] Kim B, Lee S H, Lee D and Kim Y 2020 Performance comparison of heat pumps using low global warming potential refrigerants with optimized heat exchanger designs Appl. Therm. Eng. 171 114990 [18] Peng X, Wang D, Wang G, Yang Y and Xiang S 2020 Numerical investigation on the heating performance of a transcritical CO2 vapor-injection heat pump system Appl. Therm. Eng. 166 114656 [19] Pinamonti M and Baggio P 2020 Energy and economic optimization of solar-assisted heat pump systems with storage technologies for heating and cooling in residential buildings Renew. Energ. 157 90–99 [20] Wang M, Deng C Wang Y and Feng X 2020 Exergoeconomic performance comparison, selection and integration of industrial heat pumps for low grade waste heat recovery Energy Convers. Manag. 207 112532 [21] Deng J, He S, Wei Q, Li J, Liu H, Zhang Z and Zhang H 2020 Field test and optimization of heat pumps and water distribution systems in medium-depth geothermal heat pump systems Energy Build. 209 109724 [22] Hou G, Taherian H, Li L, Fuse J and Moradi L 2020 System performance analysis of a hybrid ground source heat pump with optimal control strategies based on numerical simulations Geothermics 86 101849 [23] Isachenko V P, Osipova V, Semyonov S, Sukomel 2000 A Heat transfer (University Press of the Pacific)