DEVELOPMENT OF ENERGY- SAVING TECHNOLOGY FOR MAINTAINING THE FUNCTIONING OF HEAT PUMP POWER SUPPLY

10. Elnajjar E., Hamdan M. O., Selim M. Y. E. Experimental investigation of dual engine performance using variable LPG composition fuel // Renewable Energy. 2013. Vol. 56. P. 110–116. doi: https://doi.org/10.1016/j.renene.2012.09.048 11. Research into emissions of nitrogen oxides when converting the diesel engines to alternative fuels / Kryshtopa S., Panchuk M., Dolishnii B., Kryshtopa L., Hnyp M., Skalatska O. // Eastern-European Journal of Enterprise Technologies. 2018. Vol. 1, Issue 10 (91). P. 16–22. doi: https://doi.org/10.15587/1729-4061.2018.124045 12. Numerical study on the combustion process of a biogas spark-ignition engine / Carrera J., Riesco J., Martínez S., Sánchez F., Gallegos A. // Thermal Science. 2013. Vol. 17, Issue 1. P. 241–254. doi: https://doi.org/10.2298/tsci111115152c


Introduction
Under conditions for saving natural fuel and decreasing harmful emissions into the atmosphere, heat pump power supply with the use of renewable power sources is gaining further development [1][2][3].Thus, for example, the heat pump, using fermented must as a low-potential energy source, is recommended in order to maintain the operation of the biogas plant as a part of a cogeneration power system.Due to additional biogas production, the proposed technology provides an opportunity to increase marketability of a biogas plant and decrease the cost value of production of electricity and heat in the range of 20-30 % [3].
Under condition of the change of parameters of renewable power sources, a relevant task regarding further development of technologies of heat pump power supply is to maintain efficiency of heat pump systems based on coordination of production and power flow rate in terms of energy saving.
It is known that frequency control is based on changing the rotation rate of a compressor due to a change in supply voltage frequency but measures the change in pressure of a refrigerant in the evaporator.Digital regulation of the spiral compressor is based on modification of the period of compression of refrigerant vapor based on a specially developed electromagnetic valve, but also assesses a change in evaporation pressure.Maintaining the productivity of the heat pump system requires ensuring that changes in pressure evaporation in the evaporator should correspond to a change in pressure of refrigerant condensation in the condenser, which is difficult to execute based on measuring a change in refrigerant pressure in the evaporator.
For this purpose, it is necessary to predict a change in the temperature of local water when measuring the refrigerant temperature at the outlet of the condenser.Making advancing decisions on a change in refrigerant flow rate and the number of rotations of the motor of a compressor when measuring evaporation pressure, condensation pressure and voltage frequency will ensure a compliance of a change in evaporation pressure in the evaporator with a change in refrigerant condensation pressure in the condenser.

Literature review and problem statement
The known methods of improvement of heat pump systems -economic, exergy, thermal economy -allow determining the optimal operating conditions at the static level, which complicates the coordination of production and energy consumption under operation conditions.Thus, in paper [7], a comparative analysis of the energy and exergy methods in terms of optimization of performance coefficient of the heat pump system was made.The preference was given to the exergy method based on assessment of the influence of refrigerant flow rate on the power of the compressor.In paper [8], the authors performed mathematical modeling of components of the heat pump system regarding further modeling at the level of the heat pump system, but at the static level.Moreover, the dynamic approach to mathematical modeling of the heat pump system was considered in this paper, but only for the purpose of optimization of designing heat pump systems and stabilization of control systems.Paper [9] focuses on the dynamic approach to optimization of heat pump systems.The impact of intensification of heat exchange on the performance of the heat pump system, but with the use of mathematical models with concentrated parameters, was established.Assessment of a change in parameters only in time does not provide the possibility of predicting changes of parameters of the technological process.
Due to complexity of the analysis of heat pump systems at the static level regarding the performance optimization, the authors propose the use of combined systems of heat pump power supply with connection, for example, of photo cells [10].The influence of a change of solar radiation on maintaining the power of a compressor regarding a decrease in the term of using an external electric network, also at the static level, was described in paper [10].The technological scheme of heat pump power supply regarding the use of a heat pump heater from a solar engine was proposed in paper [11].The possibility to raise the suction pressure of a compressor and increase the efficiency of operation of the heat pump system with the use of energy and exergy methods, also at the static level, was determined.In addition to the static models, the models [12] that set the goal of control of heat pump systems with the use of intelligent networks are also known.But the means, proposed in this work, require measurement of process parameters.This approach cannot be used for making advancing decisions to maintain the operation of heat pump systems.
According to the technological scheme, it is possible to separate two circuits of heat pump power supply: external and internal.The external circuit provides extraction of heat from a low-potential power source regarding evaporation of the refrigerant in the evaporator of the heat pump.The internal circuit is made up from the evaporator, the compressor and the condenser.It is the evaporator and the condenser, located in the internal circuit, which are connecting elements of the external and internal circuits.A change in vapor flow rate that is fed to suction in the compressor of the heat pump depends both on a change in parameters of a low-potential power source and changes in flow rate.There are various means to maintain the performance of heat pump systems in terms of the impact on the compressor operation.Thus, the cyclic mode applies the on-off principle, which requires additional electricity consumption to start the compressor.Frequency control is based on changing compressor rotation rate due to a change in supply voltage frequency regarding changes in performance but measures only a change in refrigerant pressure in the evaporator.Digital control of the spiral compressor is based on the modification in the period of time of refrigerant vapor compression based of a specially developed electromagnetic valve in respect to the axial alignment of spirals, but it also assesses a change in evaporation pressure.Maintaining the performance of the heat pump system requires ensuring the compliance of a change in evaporation pressure in the evaporator with a change in refrigerant con-densation pressure in the condenser of the heat pump.With the view to energy saving, maintaining of functioning of heat pump power supply should be based on prediction of a change in local water temperature when measuring the temperature the refrigerant at the outlet of the condenser.It is necessary to make advancing decisions regarding a change in refrigerant flow rate and rotation rate of the compressor motor when measuring evaporation pressure, condensation pressure and voltage frequency when it comes to both frequency and digital control.This substantiates the need for research in this direction.

The aim and objectives of the study
The aim of present research is to develop the energy-saving technology of maintaining functioning of heat pump power supply as a component of the technological system.
The set goal will be accomplished if the following tasks are solved: -substantiation of the necessity of analytical estimates of a change in the local water temperature, refrigerant vapor flow rate.the number of rotations of the compressor motor regarding making advancing decisions on a change in performance of the heat pump system; -measuring the refrigerant temperature at the outlet of the condenser, evaporation pressure, condensation pressure, and voltage frequency; -development of the structural scheme and integrated mathematical and logical modeling as for obtaining the reference and functional estimation of a change in the local water temperature, refrigerant vapor flow rate, the number of rotations of the compressor motor; -to offer an integrated system of maintaining functioning of the heat pump power supply at the decision-making level.

1. Mathematical substantiation of the architecture of the technological system
Based on the methodological and mathematical substantiation of the architecture of technological systems [13,14], the architecture of the technological system of functioning of heat pump power supply was proposed.Its basis is the integrated dynamic subsystem, which includes the following dynamic systems: the soil heat exchanger -evaporator, the evaporator -compressor, the compressor -condenser.Other units that are the parts of the technological system are the units of charge, discharge and the unit of functional evaluation of efficiency coefficient, which are in the coordinated interaction with the dynamic subsystem (Fig. 1).
Using formula (1), mathematical substantiation of the architecture of the technological system of heat pump power supply functioning was described: where TSFHP(τ) is the technological system of functioning of heat pump power supply; IDS(τ) is the integrated dynamic subsystem (the soil heat exchanger -evaporator; evaporator -compressor; compressor -condenser); PIDS(τ), PB(τ) are the properties of the elements of the integrated dynamic subsystem, and of units of the technological subsystem, respectively; τ is the time, s; z is the length coordinate, m; x(τ) is the impacts; f(τ) is the diagnosed parameters; K(τ) is the coefficients of mathematical description; y(τ, z) is the output parameters; d(τ) is the dynamic parameters; Z(τ), R(τ) are the logical relations in IDS(τ), TSFHP(τ), respectively.
Indices: i is the number of elements of the technological system; 0, 1, 2 are the initial stationary mode, external and internal character of impacts.

2. Mathematical substantiation of maintaining the functioning of heat pump power supply
The mathematical description (2) of maintaining the functioning of heat pump power supply, based on mathematical substantiation of the architecture of technological systems, methodology of mathematical description of dynamics of power systems, the method of cause-effect relations graph was proposed [13,14]: , ( ), ( , ), ( ), ( ) where SFHP(τ) is the support of the functioning of heat pump power supply; IDS (τ) is the integrated dynamic subsystem that includes the dynamic system: the soil heat exchanger -evaporator, the evaporator -compressor, the compressor -condenser; PIDS(τ), PB(τ) are the properties of the elements of the integrated dynamic subsystem, of units of the technological system, respectively; MMIDS(τ, z) is the mathematical modeling of soil temperature dynamics, refrigerant vapor flow consumption, the number of rotations of the compressor motor, the local water temperature; CIDS(τ), MDIDS(τ), SIDS(τ) are the control of working capacity, decision making, state identification in the dynamic subsystem, respectively; LCIDS(τ), LMDIDS(τ), LSIDS(τ) are the logical relations in CIDS(τ), MDIDS(τ), SIDS(τ), respectively; FIIDS(τ) is the functional summarizing information regarding decision making in the dynamic subsystem; NCF(τ), CNCF(τ) are the new conditions of functioning, confirmation of new conditions of functioning from the units of the technological system; R(τ) is the logical relations between the dynamic subsystem and the units of discharge, charge, functional estimation of performance coefficient, which are included in the technological system of functioning of heat pump power supply.
Mathematical descriptions of (1, 2) allow maintaining functioning of the heat pump power supply using the following actions: -control of working capacity (CIDS(τ)) of the dynamic subsystem: the soil heat exchanger -evaporator, evaporator -compressor, compressor -condenser based on mathematical (MMIDS(τ,z)) and logical (LCIDS(τ)) modeling regarding obtaining the reference (MIIDS(τ)) estimation of a change in soil temperature, refrigerant flow rate, the number of rotations of the electric motor of the compressor, local water temperature; -control of working capacity (CIDS(τ)) of the dynamic subsystem: the soil heat exchanger -evaporator, evaporator -compressor, compressor -condenser based on mathematical (MMIDS(τ,z)) and logical (LCIDS(τ)) modeling regarding obtaining the functional (FIIDS (τ)) estimation of a change in soil temperature, refrigerant flow rate, the number of rotations of the electric motor of the compressor, local water temperature; -decision making (MDIDS (τ)) with the use of functional information (FIIDS (τ)), obtained based on logical modeling of obtained (LMDIDS(τ)); -decision making on a change of voltage frequency with the use of the functional estimation of a change in soil temperature, refrigerant flow rate, the number of rotations of the electric motor of the compressor, local water temperature (FIIDS(τ)); -identification (SIDS(τ)) of new conditions of functioning heat pump power supply (NCF(τ)) based on logical modeling (LSIDS (τ)) and confirmation of new conditions of functioning based on (R(τ)) from the units of the technological system.

3. Mathematical modeling of dynamics of change in the local water temperature
According to formulas (1), (2), it was proposed to predict changes in local water temperature when measuring the refrigerant temperature at the outlet of the condenser of the heat pump.Transfer function for the channel: "local water temperature -refrigerant flow rate" was obtained as a result of solving the system of nonlinear differential equations.A change in the local water temperature both in time and along the condenser length coordinate at a change in refrigerant flow rate is represented as follows: where where a is the heat transfer coefficient, kW/(m 2 •K); C is the specific thermal capacity, kJ/(kg•K); G is the substance flow rate, kg/s; g is the specific weight of substance, kg/m; h is the specific surface, m 2 /m; t w ,, σ are the temperatures of local water, distribution wall, refrigerant, K; z is the coordinate of the condenser length, m; T w , T m are the time constants, that characterize thermal accumulating capacity of local water, metal, s; m is the indicator of dependence of heat transfer coefficient on flow rate; is the time, s; S is the Laplace transform parameter; S=ωϳ; ω is the frequency, 1/s.Indices: 0 is the original stationary mode; 1 is the inlet to the condenser; w is the internal flow -local water; m is the metal wall; r is the external flow -refrigerant.
Transfer function for the channel: "local water temperature -refrigerant vapor flow rate" was obtained based on the solution of the system of nonlinear differential equations using the Laplace transform.A system of differential equations includes the equation of state as an assessment of a physical model of the dynamic system: compressor -condenser of the heat pump.The system of differential equations also includes the equation of energy of transferring and receiving media -the refrigerant and local water, respectively, and the equation of thermal balance for the wall of the condenser.
The equation of energy of the receiving medium was developed with representation of a change on the local water temperature both in time and along the spatial coordinate, which coincides with the direction of the flow of medium motion.The equation of energy of transferring medium includes coefficient K r , which assesses a change in the temperature of the refrigerant at the outlet of the condenser of the heat pump that is measured.

Mathematical modeling of dynamics of change in the refrigerant vapor flow rate
According to formulas (1), ( 2), the estimation of a change in refrigerant flow rate as for feed to compressor suction when measuring evaporation pressure at the outlet of the evaporator and condensing pressure at the outlet of the condenser of the heat pump was proposed.For this purpose, the system of differential equations, which includes the equation of refrigerant energy, local water, and the equation of thermal balance for the wall of the condenser was supplemented by equation of refrigerant continuity.The result of solving the system of differential equations using the Laplace transform is the transfer function for the channel: "flow rate of refrigerant -refrigerant pressure", which assesses a change in flow rate of refrigerant at a change of evaporation pressure and condensation pressure: where ; ; where a is the heat transfer coefficient, kW/(m 2 •K); p is the refrigerant pressure, МPа; f s is the intersection for refrigerant passing, m 2 ; C is the specific thermal capacity, kJ(kg•K); G is the substance flow rate, kg/s; ρ is the refrigerant density, kg/m 3 ; g is the specific weight of substance, kg/m; h is the specific surface, m 2 /m; z is the spatial coordinate of the condenser length, m; T w , T m are the time constants that characterize thermal accumulating capacity of local water, metal, s; S is the parameter of Laplace transform; S=ωϳ; ω is the frequency, 1/s.Indices: 0 is the original stationary mode; 1 is the inlet to the condenser; w is the internal flow -local water; m is the metal wall; r is the external flow -refrigerant.

5. Mathematical modeling of dynamics of change in the number of rotations of compressor motor
According to formulas (1), ( 2), estimation of a change in the number of rotations of the compressor motor of the heat pump was proposed.Parameters that are measured are the following: evaporation pressure at the outlet of the evaporator, condensation pressure at the outlet of the condenser and voltage frequency.The transfer function for the channel: "the number of rotations of the motor of the heat pump compressor -voltage frequency" takes the following form: where a is the heat transfer coefficient, kW/(m 2 •K); p is the refrigerant pressure, МPа; f s is the intersection for refrigerant passing, m 2 ; f is the voltage frequency, Hz; p n is the number of poles pairs; C is the specific thermal capacity, kJ/(kg•K); G is the substance flow rate, kg/s; ρ is the density of refrigerant, kg/m 3 ; g is the specific weight of substance, kg/m; h is the specific surface, m 2 /m; z is the coordinate of the condenser length, m; T w , T m are the time constants that characterize thermal accumulating capacity of local water, metal, s; S is the parameter of Laplace transform; S=ωϳ; ω is the frequency, 1/s.Indices: 0 is the original stationary mode; 1 is the inlet to the condenser; w is the internal flow -local water; m is the metal wall; r is the external flow -refrigerant.
The true part of transfer function (3) regarding estimation of a change in local water temperature was separated: Coefficient K r includes the temperature of a separating wall θ: where s 1 , s 2 are the temperatures of refrigerant at the inlet and at the outlet of condenser, K, respectively; where δ is the thickness of the wall of the condenser, m; a is the heat transfer coefficient, kW/(m 2 •K); λ is the thermal conductivity of metal of the condenser wall, kW/(m•K); t 1 , t 2 are the temperatures of local water at the inlet and outlet of the condenser, K, respectively.Indices: r is the external flowrefrigerant; w is the internal flow -local water.
To use the true part of O 1 (w), the following coefficients were obtained: The real part of transfer function (4) concerning evaluation of a change in refrigerant vapor flow rate was separated: To use the real part of O 2 (w), the following coefficients were obtained: The real part of transfer function ( 5) regarding evaluation of a change in the number of rotations of the compressor motor was separated: To use the real part of O 3 (w), the following coefficients were obtained: Implementation of transfer functions (3) to ( 5), obtained based on the use of the operational method for solving the systems of nonlinear differential equations, maintain the parameter of Laplace transform -S(S=ωϳ), where ω is the frequency, 1/s.To transfer from frequency region to time region, real parts (6), ( 12), ( 16), obtained as a result of mathematical treatment of transfer functions, were separated.These parts are included in integrals (20)-( 22), providing the opportunity to obtain dynamic characteristics of the local water temperature, refrigerant flow rate, the number of rotations of the compressor motor, respectively, using the reverse Fourier transform.
where t w is the temperature of local water, K; G r is the refrigerant flow rate, kg/s; n is the number of rotations of the electric motor of the compressor, rpm.

Results of research into technology of maintaining the functioning of heat pump power supply
According to formulas (1), ( 2), an integrated mathematical modeling of the heat pump power supply with the use of the developed block diagram was performed (Fig. 2).
According to the proposed block diagram (Fig. 2), Tables 1-3 shows the results of an integrated mathematical modeling of heat pump power supply.Thermal and physical properties of refrigerant R 134a were used.Estimation of the boundary change of the refrigerant temperature at the outlet of the evaporator of the heat pump was established based on maintaining soil heat discharge at the level 10-8 °C [6].
Based on the proposed mathematical substantiation of maintaining functioning of heat pump power supply (1), ( 2), the block diagrams (Fig. 3, 4) regarding working capacity and maintaining the efficiency of the heat pump system was developed.
The comprehensive integrated system of maintaining functioning of the heat pump power supply was developed (Tables 4, 5).The system is based on prediction of a change in local water temperature regarding a change in refrigerant vapor flow rate, the number of rotations of the electric motor of the compressor.Temperature of the refrigerant at the outlet of the condenser, evaporation pressure, condensation pressure and voltage frequency are measured continuously.i is the number of levels of functioning; τ is the time, s.Index: c. est.up. is the constant estimated value of the parameter of the upper level of functioning.Thus, for example, in the period of time of 94.5•10 5 s (2625 hours) after the beginning of a heating season, which was selected for the city of Kyiv (Ukraine) and is 4,448 hours, absolute value of local water temperature using formula (23) is: 40.17 °С=39.12°С+(0.2585-0.2059)(55°С-35 °С).Fig. 3. Block diagram of working capacity control of the heat pump system: t r out.t w , q are the refrigerant temperatures at the outlet of the condenser, of local water, of separating wall, respectively, K; p evap., p cond. are the evaporation pressure, condensation pressure, respectively, МPа; G r is the refrigerant flow rate, kg/s; f is the voltage frequency, Hz; n is the number of rotations of the electric motor of the compressor, rpm; CT is the event control; Z is the logic relations; d is the dynamic parameters; x is the impacts; f d is the diagnosed parameters; y is the original parameters; K is the coefficients of mathematical description; COP is the coefficient of efficiency of the heat pump; i is the time, Indices: с is the working capacity control; е is the reference value of a parameter, c. est.f. is the constant, estimated value of the parameter of the first level of functioning; c. est.up.is the constant, estimated value of the parameter of the upper level of functioning; 0, 1, 2 is the initial stationary mode, external, internal parameters; 3 is the coefficients of equations of dynamics; 4 is the essential diagnosed parameters; 5 is the dynamic parameters Refrigerant flow rate in the established period of time is determined as follows: where G is the refrigerant flow rate, kg/s; G 1 , G 2 is the initial and final values of refrigerant flow rate, kg/s, respectively (Table 5); i is the number of levels of functioning of heat pump power supply; τ is the time, s.Index: c. est.f. is the constant estimated value of the parameter of the first level of functioning The number of rotations of the electric motor of the compressor in the established period of time is determined as follows: ) ) where n is the number of rotations of the electric motor of the compressor, rpm; n 1 , n 2 are the initial and final values of the number of rotations of the electric motor of the compressor, rpm, respectively; i is the number of levels of functioning of the heat pump power supply; τ is the time, s.Index: c. est.f. is the constant estimated value of the parameter of the first level of functioning.Thus, for example, in the period of time of 94,5•10 5 s (2625 hours) after the beginning of a heating season, which was selected for the city of Kyiv (Ukraine) and is 4448 hours, absolute values of refrigerant flow rate, of number of rotations of the electric motor of the compressor with the use of formulas (24), (25) are: 0.0373 kg/s=0.0369kg/s+(0.9610--0.8291)(0.0398kg/s-0.0368kg/s), 1184 rpm=1136.7 rpm+(0.9395--0.8334)(1500rpm-1053.9rpm).
Graphic dependence of a change in local water temperature within the heating period regarding making decisions on the change in efficiency of the heat pump power supply is shown in Fig. 5.
Thus, for example, in the period of time 94.5•10 5 s (2625 hours) regarding heating local water from 40.17 ºC to 45.36 °C at digital control, it is necessary to make the advancing decision to increase refrigerant flow rate up to the level of 0.0373 kg/s.Frequency control needs making additional decision to increase voltage frequency by 13.2 %, to increase the number of rotations of the electric motor of the compressor up to the level of 1184 rpm. .Block diagram of maintaining efficiency of heat pump power supply: t r out, t w are the refrigerant temperatures at the outlet of the condenser, of local water, respectively, K; p evap., p cond. are the evaporation pressure, condensation pressure, respectively, МPа; G r is the refrigerant flow rate, kg/s; f is the voltage frequency, Hz; n is the number of rotations of the electric motor of the compressor, rpm; COP is the coefficient of efficiency of the heat pump system; i is the time, Indices: e is the reference value of the parameter; c. est.f. is the constant, estimated value of the parameter of the first level of functioning; c. est.up. is the constant, estimated value of the parameter of the upper level of functioning; l., n. l. are the level, new level of parameters, respectively Such actions will make it possible to increase vapor compression in the compressor by increasing the vapor temperature at the inlet to the condenser up to the level of 50 °C.Correctness of subsequent local water heating is proved by the use of an integrated system of maintaining the discharge of soil, the temperature of which in this period of time is 9.21 °C [6].
6. Discussion of results of studying the energy-saving technology for maintaining the functioning of heat pump power supply The integrated system of maintaining the efficiency of the heat pump power supply based on coordination of power production and consumption under power saving conditions.
During control of heat pump power supply, measurement of refrigerant pressure in the heat pump evaporator is usually used.But the use of such measurement is substantiated only under condition of evaluation of a change in refrigerant temperature and under condition when coordination of power production and consumption is not considered.When trying to go beyond these limits in order to enhance the efficiency of the heat pump power supply, the on-off functioning mode is used.When using frequency or digital control, an additional complication of electric circuits is used.There occur objective difficulties, associated with considerable costs of both electrical energy for starting processes and additional capital investments.The way of overcoming these difficulties is presented.It is based on the fact that prediction of a change in local water temperature with measurement of the refrigerant temperature at the outlet of the condenser of the heat pump was proposed.This estimation in the ratio of the measured evaporation pressure is a part of the proposed analytical determining of the refrigerant flow rate as for the digital control and the number of rotations of the electric motor of the compressor for frequency control.This estimation provides an opportunity to influence in advance the coordination of functioning of the external and internal circuits of the heat pump system, i.e. energy production and consumption.The outcome of the conducted research is the results of a comprehensive mathematical and logical modeling of the heat pump power supply using the boundary change in refrigerant temperature at the outlet of the heat pump evaporator, determined on the basis of the integrated system of maintaining soil heat discharge at the level of 10-8 °C.Thus, we established the relationship of a change in refrigerant evaporation temperature in the evaporator as a part of the internal circuit of the heat pump system and an external circuit as for maintaining of a change in soil temperature within the heating period.
Making advancing decisions on a change in vapor flow rate, influencing the action of an electromagnetic valve for digital control, allows ensuring a change of compression of refrigerant suction steam into the compressor and a change in compression of vapor injection into the condenser.Making advancing decisions to change voltage frequency as for a change of the number of rotations of the electric motor of the compressor also allows setting the compliance of compression of vapor suction with a change in injection compression at frequency control.Establishment new vapor parameters regarding the transition to a new level of operation maintains the new coefficient of efficiency of the heat pump system.The results of the study are the continuation of the research in the direction of coordination of production and energy flow rate [3,6,13,14].The presented results can be used in the development of intelligent systems of functioning of controllers of heat pump systems with both frequency and digital control.The development of this research is the planned testing of the research under the conditions of using heat pump systems of different types and different capacities as a part of the proposed technological system.This is due to the use of different sources of low-potential energy for determining various means of maintaining of a change of parameters during the technological process as, for example, in energy saving technology of biogas production based on the heat pump power supply, the energy source for which is fermented must [3].When using heat pump systems of various power and improving intelligent systems of controllers functioning, it is necessary to establish the means of control: frequency or digital, in order not to cause unnecessary additional capital investment.

Conclusions
1. Prediction of a change in local water temperature at measuring the refrigerant temperature at the outlet of the condenser was proposed.Making advancing decisions on changing the refrigerant flow rate and the number of rotations of the electric motor of the compressor at measuring evaporation pressure, condensation pressure and voltage frequency ensures compliance of a change in evaporation pressure in the evaporator with a change in refrigerant condensation in the condenser.
2. The comprehensive mathematical modeling of the heat pump system, based on an integrated system of maintaining the soil heat discharge at the level of 10-8 °C.The boundary change of the refrigerant temperature at the outlet of the evaporator was established: 0 o С…(-1,5) o С. Mode parameters of the heat pump system, parameters of heat exchange in the condenser, time constants and coefficients of mathematical models of dynamics for the established levels of functioning were determined.The standard dynamic estimations of a change in local water temperature, refrigerant vapor flow rate, the number of rotations of the electric motor of the compressor were obtained.Logical modeling of control of working capacity of the heat pump system, which follows the causation principle, was performed.
The logical unit has the components that evaluate: a change of refrigerant temperature at the outlet of the condenser, evaporation pressure, condensation pressure, voltage frequency that are measured; a change in the temperature of the condenser wall; a change of coefficients of mathematical models of dynamics, ϰ p , K f , K х ; a change in local water temperature, refrigerant flow rate, the number of rotations of the electric motor of the compressor; a change of dynamic parameters; the resulting unit of working capacity control for obtaining a functional estimation of a change of local water temperature, refrigerant vapor flow rate, the number of rotations of the electric motor of the compressor.
3. Maintaining the functioning of the heat pump power supply for prediction of a change in local water temperature at continuous measurement of the refrigerant temperature at the outlet of the condenser, evaporation pressure, condensation pressure and voltage frequency was proposed.Maintaining a change in performance of the heat and pump system is based on a comparison of evaporation pressure, condensation pressure, voltage frequency that are measured with the reference values.Determining of summarizing functional information provides an opportunity to make the following advancing decisions: to maintain a change in evaporation pressure as for a change in refrigerant vapor flow rate for digital control: to maintain a change in evaporation pressure as for a change of refrigerant vapor flow rate and a change of voltage frequency as for a change of the number of rotations of the electric motor of the compressor for frequency control.

Fig. 1 .
Fig. 1.The architecture of the technological system of functioning of heat pump power supply: 1 -charge unit; 2 -unit of functional estimation of performance coefficient; 3 -discharge unit

Fig. 2 .
Fig. 2. Block diagram of the integrated mathematic modeling of heat pump power supply: U -voltage, V; f -voltage frequency, Hz; I -current, А; d e , d i -external, internal diameter of the soil thermal exchanger, respectively, m; L -length of heat exchanger, m; t r , t w , t evap -temperature of refrigerant, local water, refrigerant evaporation, respectively, K; t r in , t w out -temperature of refrigerant at the inlet, local water at the outlet, respectively, K; Gr -refrigerant flow rate, kg/s; n -number of rotation of the electric motor of the compressor, rpm; COP -coefficient of performance of the heat pump system

Fig. 5 .
Fig. 5. Maintaining efficiency of heat pump power supply, where 1, 2, 3, 4 are the first, second, third and fourth levels of functioning, respectively

Table 1
Mode parameters of the heat pump system

Table 2
Parameters of heat exchange in the composition of mathematical models of dynamics Note: α r is the coefficient of convective heat exchange from the refrigerant to the wall of condenser, W/(m 2 •K); α w is the coefficient of convective heat exchange from the condenser wall to local water, W/(m 2 •K); k is the heat transfer coefficient, W/(m 2 •K)

Table 3
Values of time constants and coefficients of mathematical models of dynamics

Table 5
Integrated system of a change in refrigerant flow rate and the number of rotations of the compressor Gr is the refrigerant flow rate.kg/s; n is the number of rotations of the electric motor of the compressor, rpm; pevap. is the evaporation pressure, MPa; pcond. is the condensation pressure, МPа; f is the voltage frequency, Hz; τ is the time, s.Index: c. est.f. is the constant estimated value of the parameter of the first level of functioning