CC BY
7
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
-□ □-
Розроблено ттегровану систему тдтримки функщонування теплонасосного енергопостачання на основi прогнозування змти температури м^цевог води. Змта витрати пари холодагента, числа обертiв електродвигуна компресора вiдбуваeться при вимi-рювант температури холодагента на виходi iз конденсатора, тиску випаровування, тиску конденсащг та частоти напруги.
Виконано комплексне математичне моделювання теплонасосног системи, що базуеться на ттегроватй системi тдтримки розряду грунту нарiвнi 10-8°С. Визначено витрату холодагента, потужтсть електродвигуна компресора, напругу, частоту напруги, число обертiв електродвигуна компресора, коефщ1ент продуктивностi теплонасосног системи для встановлених рiвнiв функщонування. Встановлено параметри конвективного тепло-обмту в конденсаторi, постшт часу та коефщгенти математич-них моделей динамжи змши температури мкцевог води, витрати пари холодагента, числа обертiв електродвигуна компресора.
Здобуто функщональну ощнку змти температури мкцевог води в дiапазонi 35-55 °С впродовж опалювального сезону, витрати пари холодагента, числа обертiв електродвигуна компресора. Визначення тдсумковог функциональной тформацп надае можливкть прийма-ти наступт випереджуючiргшення: на тдтримку змти тиску випаровування щодо змти витрати пари холодагента для цифрового управлтня; на тдтримку змти тиску випаровування щодо змти витрати пари холодагента та на змту частоти напруги щодо змти числа обертiв електродвигуна компресора для частотного управлтня. Тому, запропоновано прогнозування змти температури мкце-вог води на основi вимiрювання температури холодагента на виходi iз конденсатора. Саме ця ощнка у стввiдношеннi з вимiрюваним тиском випаровування, входить до складу аналхтичних визначень витрати холодагента та числа обертiв електродвигуна компресора. Здобуття таког ощнки та вимiрювання частоти напруги надае можливкть упереджено впливати на узгодження функщонування зовтшнього та внутршнього контурiв теплонасосног системи як при цифровому, так i частотному управлтт
Ключовi слова: теплонасосна система, частотне управлтня,
цифрове управлтня, тиск випаровування, тиск конденсащг -□ □-
UDC 621.31
| DOI: 10.15587/1729-4061.2018.139473 |
DEVELOPMENT OF ENERGY-SAVING TECHNOLOGY FOR
MAINTAINING THE
FUNCTIONING OF HEAT PUMP POWER SUPPLY
E. Chai kovskaya
PhD, Associate Professor, Senior Researcher Department of Theoretical, general and alternative energy Odessa National Polytechnic University Shevchenka ave., 1, Odessa, Ukraine, 65044 E-mail: eechaikovskaya@gmail.com
1. 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-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].
The outside air as a low-potential power source is available in heat pump power supply, but the change of air temperature in a wide range makes it difficult to maintain functioning of heat pump systems [1].
A heating system of the soil-water type requires the construction of special soil heat exchangers for heating the brine - 30 % ethylene glycol solution that is fed to the evaporator of a heat pump [2]. A vertical heat exchanger occupies a smaller area than a horizontal one but requires additional capital investments for drilling wells. It is possible to ensure heat extraction within 30-100 W per one meter of the length of a vertical heat exchanger depending on a soil type at the depth of 40-150 m, where the soil temperature
is 8-10 °C. The design of a vertical heat exchanger of the tube-in-tube type is most common. Brine, after transferring heat to a refrigerant in the evaporator of a heat pump, is fed by the circulation pump to the soil heat exchanger on one line, the other line provides a lift of the brine to the evaporator of a heat pump. Improvement of the tools for using the soil heat discharge is based on the development of new structures of vertical heat exchangers for the intensification of the heat exchange process. Thus, for example, in paper [4], the optimal structural parameters, such as length and diameter, were established based on mathematical modeling of a specially developed spiral heat exchanger. In paper [5], it was proposed to use the foundation for location of spiral heat exchangers with the purpose of optimization of a heat exchange surface. Heat pumps of the soil-water type are most common, but their use should not disturb the natural soil relaxation in the period when the heating does not function. Thus, an integrated system of maintaining soil heat discharge within 8-10 °C during the heating season as part of a dynamic system: exchanger-evaporator of the heat pump was proposed for this purpose. Based on prediction of a change in soil temperature, changes in brine flow rate were ensured based on the change of rotation rate of the motor of the circulation pump at measuring the brine temperature at the outlet of the evaporator. The exact initial and final term of soil heat discharge was established [6].
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.
2. 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.
where TSFHP(t) is the technological system of functioning of heat pump power supply; IDS(t) is the integrated dynamic subsystem (the soil heat exchanger - evaporator; evaporator - compressor; compressor - condenser); PIDS(t), PB(t) are the properties of the elements of the integrated dynamic subsystem, and of units of the technological subsystem, respectively; t is the time, s; z is the length coordinate, m; x(t) is the impacts; /(t) is the diagnosed parameters; K(t) is the coefficients of mathematical description; y(T, z) is the output parameters; d(T) is the dynamic parameters; Z(t), R(t) are the logical relations in IDS(t), TSFHP(t), respectively.
3. 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.
4. Materials and methods of the study regarding maintaining the functioning of heat pump power supply
4. 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:
/x0(t), ^(t), x2(t),1 IDS(t)(PIDS(t)l f (t),K(t),y(t,z), TSFHP (t) =•{ \d (t) I I. (1)
Z (t), PIDS (t)), R(t),(PBi (t) (x1(t), f (t), Kt (t), y (t)»,
— Dynamic subsystem:
the soil heat exchanger - evaporator; the evaporator - compressor; the compressor condenser
snse^)
Ô
ô ©
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
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.
4. 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]:
(IDS (t)(PIDS (t), MMIDS (t, z), CIDS(t), LCIDS(t) Xo(t), Xj(t), x2(t), / (t), K (t),\ y(t, z), d(t), FIIDS (t) /'
LMDIDS(t), MDIDS (t), NCF (t), SIDS(t), LSIDS(t) / (t), K (t), y(t, z), d (t), FIIDS (t)) PIDS (t))), R(t),
■^(txf (t);
K (t), y, (t)/
SFHP(t) =
(PB,(t)( ^V'\)CNCF(t)),
(2)
where SFHP(t) is the support of the functioning of heat pump power supply;
IDS (t) is the integrated dynamic subsystem that includes the dynamic system: the soil heat exchanger - evaporator, the evaporator - compressor, the compressor - condenser;
PIDS(t), PB(t) are the properties of the elements of the integrated dynamic subsystem, of units of the technological system, respectively;
MMIDS(t, 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;
MIIDS(t) is the boundary permissible change in soil temperature, refrigerant vapor flow rate, the number of rotation of the compressor motor, the local water temperature;
CIDS(t), MDIDS(t), SIDS(t) are the control of working capacity, decision making, state identification in the dynamic subsystem, respectively;
LCIDS(t), LMDIDS(t), LSIDS(t) are the logical relations in CIDS(t), MDIDS(t), SIDS(t), respectively;
FIIDS(t) is the functional summarizing information regarding decision making in the dynamic subsystem;
NCF(t), CNCF(t) are the new conditions of functioning, confirmation of new conditions of functioning from the units of the technological system;
R(t) 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.
x(t) - impacts;
/(t) - diagnosed parameters;
K(t) - coefficients of mathematical description;
y(T,z) - output parameters;
d(t) - dynamic parameters;
z - length coordinate, m;
t - time, s. Indices: i is the number of elements of SFHP(t); 0, 1, 2 are the initial mode, external, internal character of impacts.
Mathematical descriptions of (1, 2) allow maintaining functioning of the heat pump power supply using the following actions:
- control of working capacity (CIDS(t)) of the dynamic subsystem: the soil heat exchanger - evaporator, evaporator - compressor, compressor - condenser based on mathematical (MMIDS(t,z)) and logical (LCIDS(t)) modeling regarding obtaining the reference (MIIDS(t)) 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(t)) of the dynamic subsystem: the soil heat exchanger - evaporator, evaporator - compressor, compressor - condenser based on mathematical (MMIDS(t,z)) and logical (LCIDS(t)) modeling regarding obtaining the functional (FIIDS (t)) 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 (t)) with the use of functional information (FIIDS (t)), obtained based on logical modeling of obtained (LMDIDS(t));
- 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 (t));
- identification (SIDS(t)) of new conditions of functioning heat pump power supply (NCF(t)) based on logical modeling (LSIDS (t)) and confirmation of new conditions of functioning based on (R(t)) from the units of the technological system.
4. 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:
Wt
tw -G,, '
Kr e(1 - L/) (TwS + 1)P-
1(1 -
(3)
where
K =-
G,n
a r 0hr 0
a.„ A,
L= L +1
1 t G Cv * T
; Lr =—e = e(1 -Lr );
a,n h.
(TWS + 1)P-1, _£_, T = gwCw
P ' LJ " awok
P = T S + e* +1; T =
" m ' m
gmCm
a,„ nh„
where a is the heat transfer coefficient, kW/(m2-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, m2/m; tw, 0, a are the temperatures of local water, distribution wall, refrigerant, K; z is the coordinate of the condenser length, m; Tw, Tm 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; i is the time, s; S is the Laplace transform parameter; S=wj; ro 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 Kr, which assesses a change in the temperature of the refrigerant at the outlet of the condenser of the heat pump that is measured.
4. 4. 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:
W
G _pi Y I1 '
,_y .m
(4)
where
(T S + 1)B_ 1 x =_ f — ; y = -——---;
Xp Js dp r L—p '
T =
w __ 1
a, h
L =
GC„
a, h,
P = T S + £* +1; T =
I m 'm
gmC m
a, h
e* = e(1 _Lj); e =
h.
L * _—; L _ —,
r L +1 r a 0K 0 Lw
where a is the heat transfer coefficient, kW/(m2-K); p is the refrigerant pressure, MPa; fs is the intersection for refrigerant passing, m2; C is the specific thermal capacity, kJ(kg-K); G is the substance flow rate, kg/s; p is the refrigerant density, kg/m3; g is the specific weight of substance, kg/m; h is the specific surface, m2/m; z is the spatial coordinate of the condenser length, m; Tw, Tm are the time constants that characterize thermal accumulating capacity of local water, metal, s; S is the parameter of Laplace transform; S=roj; ro 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.
4. 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:
W- /,= f1 (1 _ «
-Y tU
(5)
K = 60/(1 _ s) ; x
Kf = ; xp
T— = —
g—C— a, h
L =
= _f dp; y=(T—S+1)P-1;
f dp ' Y L— p ;
GC„
a, h
P = Tm S + e* +1; Tm =
gmCm
e* = e(1 _ Lr *); e =
an h.
L*=
Lr +1
; Lr =
GrCr an h.
_ (TWS+1)p-1, Y1 _ p ' L '
r vn
where a is the heat transfer coefficient, kW/(m2-K); p is the refrigerant pressure, MPa; fs is the intersection for refrigerant passing, m2; f is the voltage frequency, Hz; pn is the number of poles pairs; C is the specific thermal capacity, kJ/(kg-K); G is the substance flow rate, kg/s; p is the density of refrigerant, kg/m3; g is the specific weight of substance, kg/m; h is the specific surface, m2/m; z is the coordinate of the condenser length, m; Tw, Tm are the time constants that characterize thermal accumulating capacity of local water, metal, s; S is the parameter of Laplace transform; S=roj; ro 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:
^ + (MA)Kre(1 _ O . ^ (A2+B12)
(6)
Coefficient Kr includes the temperature of a separating wall 0:
e = («—(o1 +o2)/2) + A(t1 +t2)/2)/(a,+ A),
(7)
where 01, o2 are the temperatures of refrigerant at the inlet and at the outlet of condenser, K, respectively;
A = 1/(8 m / X m + 1/a ),
(8)
where S is the thickness of the wall of the condenser, m; a is the heat transfer coefficient, kW/(m2-K); X is the thermal conductivity of metal of the condenser wall, kW/(m-K); t1, t2 are the temperatures of local water at the inlet and outlet of the condenser, K, respectively. Indices: r is the external flow -refrigerant; — is the internal flow - local water.
To use the true part of 01(ra), the following coefficients were obtained:
A1 = e * _T—Tmw2; A2 = e* +1; B1 = T— ew + T— w + Tm w;
A1A2 + B1B2
A2 + B22
B2 = TmW*' C1 =
L1 = 1 _ e<C •cos(_^D1); M1 = _e-çC'sin(_ÇD1).
D1 =
a2 b1 _ ab2
A22 + B22
(9)
(10)
(11)
The real part of transfer function (4) concerning evaluation of a change in refrigerant vapor flow rate was separated:
a
1
a
n
02 (w) = X pLw (C1L1) - (DM).
To use the real part of O2(w), the following coefficients were obtained:
(12) 5. Results of research into technology of maintaining the functioning of heat pump power supply
A =-Tmw2; A2 = £-TwTmw2; B1 = (£ + 1)w;
B2 = TwEW + TwW + TmW + £,
C = A1A2 + B1B2 , _ = A2B1 - A1B2 . C1 = 42 , D2 ' D1 =
A22 + B22
A22 + B22
L1 = 1 - e~zc • cos(-|D1); M1 = -e^ sin(-^D1).
(13)
(14)
(15)
The real part of transfer function (5) regarding evaluation of a change in the number of rotations of the compressor motor was separated:
O3 (w) = Kf x pLw (C1L1) - (D1M1).
To use the real part of O3(w), the following coefficients were obtained:
A =-Tmw2; A2 = £-TwTmw2; B1 = (£ + 1)w; B2 = Tw £W + Tw w + TmW + £, C = A1A2+B1B2 ; D = A2B1 - AB2;
C1 = /12 , D2 ' D1 = a2 , r,2 '
A22 + B2
A22 + B22
L1 = 1 - e~?C 1 cos(-ÇD1); M1 = -e~zc sin(-^D1).
(18)
(19)
n(x) = Gr (t, z)Kf (t) = 1 ~
= — JO3 (w)sin(tw / w)d w,
(22)
where tw is the temperature of local water, K; Gr is the refrigerant flow rate, kg/s; n is the number of rotations of the electric motor of the compressor, rpm.
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.
(16) 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
(17) pressure and voltage frequency are measured continuously.
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=roj), 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.
1 M
tw (t, z) = —J O1 (co)sin(tra / w)d w; (20)
2n 0
Gr (t, z) = — J O2 (co) sin( tw / w )dw; (21)
Boundary change of /,-at the outlet of the evaporator: 0°C...(-1,5) °C
Boundary change of parameters: to the condenser: 45.. .60 "C; K out of the condenser: 35.. .55 °C
First level: ievap=(-2,5)"C; /,"45-40,5 "C; /„=35-40 "C. Second level: /cVaP =(-3)°C; /,.=50^15,5 °C; /„,=35^15 °C. Third level: /cvap=(-3,5)°C; /,.=55-50,5 "C; /„=35-50 "C. Fourth level: /c»aP =H)"C; /,.=60-55,5 "C; /„.=35-55 "C.
Determining G'„ kg/s, Nc, kW, U, V,/ Hz, n, rpm, COP
Development and implementation: Wa , W:l_fr Wt _G |
Determining of parameters of heat exchange, time constants and coefficients of transfer functions
Reference dynamic characteristics, based on Wc _p , Wn_ft, Wt G ^
Fig. 2. Block diagram of the integrated mathematic modeling of heat pump power supply: U — voltage, V; f — voltage frequency, Hz; /— current, A; de, d) — external, internal diameter of the soil thermal exchanger, respectively, m; L — length of heat exchanger, m; tr, tw, tevap — temperature of refrigerant, local water, refrigerant evaporation, respectively, K; tr!n, twout — 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
Levels of functioning of the heat pump system -1
Reference dynamic characteristics based on WG , Wn ^, W\
Measurement of parameters: /?cvap (x), tIctnt(j),pcani{i),f{ï) (Cr1(t),(x0(t),x1(x)(/7 (x)(-)(+), jpcond(t)(+)(-),Gr(t)(0),«(t)(0),/(x)(-)(+))))JZci(t),
(fT2(T)(/d(T)((9(T)(+)(-)(/revap(T)(+)(-))))))),Zc2(T), (C/; (t)(K(T)(Xp (t)(-X+), Kf (-)(+), KT (x)(+)(-)))), Zc3 (T),
(CT4 ( T)(,v(x)(ör (x )(+)(-), n(i )(+)(-), C (x)(+)(-)))), Zi4(x),
(C7i(T)№XAGr(r) / AGrcesir(f)(+)(-),
A«(r)/A«cestf(r)(+)(-),
A'w(r) / A/wccstup (r)(+)(-)))),Zc5(x),
Obtaining summarizing information
Fig. 3. Block diagram of working capacity control of the heat pump system: tr out. tw, 8 are the refrigerant temperatures at the outlet of the condenser, of local water, of separating wall, respectively, K; pevap., pcond. are the evaporation pressure, condensation pressure, respectively, MPa; Gr 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; CTis the event control; Z is the logic relations; d is the dynamic parameters; x is the impacts; fd 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; t is the time, Indices: c is the working capacity control; e 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
Table 1
Mode parameters of the heat pump system
Levels of functioning Gr, kg/s Ne, kW U, V f, Hz n, rpm. COP
Level 1 0.0368 1. 23 266. 96 35.13 1053.9 5.13
Level 2 0.0380 1.39 302.17 39.8 1194 4.53
Level 3 0.0391 1.55 337.17 44.36 1330.8 4.06
Level 4 0.0398 1.73 380 50 1500 3.63
Note: Gr is the refrigerant flow rate, kg/s; Ne is the power of electric motor of compressor, kW; U is the voltage, V; f is the voltage frequency, HZ; n is the number of rotations of electric motor of compressor, rpm; COP is the coefficient of efficiency of heat pump system
Table 2
Parameters of heat exchange in the composition of mathematical models of dynamics
Levels of functioning Parameter
ar, W/(m2-K) a,, W/(m2-K) k, W/(m2-K)
Level 1 1193.16 873.18 496.38
Level 2 1245.17 972.35 536.83
Level 3 1319.06 1129.83 597.21
Level 4 1467.69 1423.96 706.81
Note: ar is the coefficient of convective heat exchange from the refrigerant to the wall of condenser, W/(m2K); aw is the coefficient of convective heat exchange from the condenser wall to local water, W/(m2K); k is the heat transfer coefficient, W/(m2K)
Table 3
Values of time constants and coefficients of mathematical models of dynamics
Levels of functioning Tw, s Tm, s L,,, m £ Lr., m L ' £* Z
Level 1 5.36 1.78 23.37 1.5883 2.60 0.28 1.1460 0.72
Level 2 5.02 1.60 20.98 1.5523 2.59 0.28 1.1177 0.67
Level 3 4.14 1.38 18.06 1.3570 2.57 0.28 0.9770 0.62
Level 4 3.28 1.90 14.31 1.1980 2.44 0.28 0.8506 0.60
Local water temperature at the established period of time is determined as follows:
W+i (t) = ^ + (Ki+1 (t) / Aiw,.^. (t) -
-Aiw (t) / Atw ^.(t))^-twi)), (23)
where tw is the local water temperature, °C; tw1, tw2 are the initial and final values of local water temperature, °C, respectively; i is the number of levels of functioning; t 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-105 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 °C=39.12 °C+(0.2585-0.2059)(55 °C-35 °C).
Levels of functioning of the heat pump system
>
Obtaining summarizing information
Yes
(CTc(t)(AG(t) / AGcestf (r) > AGcestl (r) / AGC e5t f (r))) < 1 > 0).
/'evapei+l
Digital control
(Cr.(r)(A/ï(r) / A«„slf(r) > A»„sll (r)/ A»„sl.r.(r))) < 1 > 0).
/(T)=/c,+ l(T)
Frequency control
(CTc( r)( AZW ( r) / A/w ccs up.W>A'„.c.ci i W/A>wxx s,.»p.M))<i>0).
(G7'(r)(AG(r) / AG. stf(r)<AGces L (rV AGc es t(r)))<l>0).
/^cvap(^) Pcvap e i-
Digital control
(C7;(r)(A«(r) / A«,csU-.(r) < A/îocsll (r) / AnccslJ. (r))) < 1 > 0).
yW=/e,-.(T)
Frequency control
(CTc(r)(A/w(r)/ A/.^ir) < A/wccill (r)/ Afw .(*"))) < 1 >0).
COP(r)=COP„y (t) —
Maintaining efficiency of heat pump power supply
n+i (x) = n, + ((^n+i (x)/ ^est.f.. (x)--An, (x)/ kncest.f (x))(n2- ni)), (25)
Fig. 4. Block diagram of maintaining efficiency of heat pump power supply: tr out, tw are the refrigerant temperatures at the outlet of the condenser, of local water, respectively, K; pevap., pcond. are the evaporation pressure, condensation
pressure, respectively, MPa; Gr 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
Refrigerant flow rate 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; nj, n2 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; t 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,540s 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.0369 kg/s+(0.9610-
-0.8291)(0.0398 kg/s-0.0368 kg/s),
1184 rpm=1136.7 rpm+(0.9395-
-0.8334)(1500 rpm-1053.9 rpm).
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-105 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.
(x) = G
-((AG,+i (x) / AG_t f. -AG, (x) / AG_t f.(x))(G2-Gi )),
(24)
where G is the refrigerant flow rate, kg/s; Gj, G2 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; t 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:
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 4
Integrated system of maintaining a change in the local water temperature
Table 5
Integrated system of a change in refrigerant flow rate and the number of rotations of the compressor
Time, t, 105 s A change in parameters At^(x)/ Atw. c. est. up. (t) tw(T), °C
0 tr in 45 C; tr out.=39.55 C; /=35.13 Hz 0 35
10.5 tr in 45 C, trout=39.55 C; /=34.55 Hz 0.0086 35.17
21 tr in=45 °C; tr out=39.6 °C; /=33.97 Hz 0.0218 35.43
31.5 tr in=45 °C; tr out=39.7 °C; /=33.38 Hz 0.0489 35.97
42 tr in=45 °C, tr out=39.8 °C; /=32.8 Hz 0.0744 36.48
52.5 trin=45 °C; tr out=39.9 °C; /=32.21 Hz 0.1007 37.01
63 trin=45°C; trout=40°C; /=31.63 Hz 0.1270 37.54
73.5 trin=45 °C; tr out=40.1 °C; /=31.05 Hz 0.1533 38.07
84 tr in=45 °C; tr out=40.3 °C; /=30.46 Hz 0.2059 39.12
94.5 Making decisi°n °n a change in v°ltage frequency: /=39.8 Hz; tr in=45 °C; tr out=40.5 °C 0.2585 40.17
105 trin=50 °C; tr out=44.2 °C; /=37.52 Hz 0.3282 41.56
115.5 Making decisicn °n a change in v°ltage frequency: /=44.36 Hz; tr in=50 C; tr out.=45 C 0.5184 45.36
126 tr in=55 °C; tr out=50.1 °C; /=30.46 Hz 0.7865 50.73
136.5 Making decisicn °n a change in v°ltage frequency: /=50 Hz; tr in=60 °C; tr out.=55.1 °C 0.9283 53.57
147 tr in=60 C; tr out=55.4 °C; /=50 Hz 0.9821 54.65
157.5 tr in=60 C; tr out=55.5 °C; /=50 Hz 1 55
161.6 tr in=60 °C; tr out=55.5 °C; /=50 Hz 1 55
Note: tr in, tr out, tw are the refrigerant temperature at the inlet of the condenser, at the outlet of the condenser, local water temperature, respectively, °C; f is the voltage frequency, Hz; t is the time, s. Index:c. est. up. is the constant estimated value of the parameter of the upper level of functioning
Time, t, 105 s A change in parameters AG(x)/ AG c.est.f.(T) An(x)/ An c .est. f. (T) G(t), kg/s n (t), rpm
0 pevap=0.266 MPa; Pc°nd=1.008 MPa; / =35.13 Hz 1 1 0.0368 1053.9
10.5 pevap=0.265 MPa; Pc°nd=1.009 MPa; /=34.55 Hz 0.9979 0.9815 0.0368 1062.1
21 pevap=0.265 MPa; Pc°nd=1.012 MPa; / =33.97 Hz 0.9932 0.9605 0.03681 1071.5
31.5 pevap=0.2645 MPa; Pc°nd=1.014 MPa; /=33.38 Hz 0.9899 0.9505 0.03682 1076
42 Pevap=0.264 MPa; Pc°nd=1.016 MPa; /=32.8 Hz 0.9866 0.9341 0.03683 1083.3
52.5 pevap=0.2635 MPa; Pc°nd=1.019 MPa; /=32.21 Hz 0.9820 0.9005 0.03684 1098.3
63 pevap=0.263 MPa; Pc°nd=1.022 MPa; /=31.63 Hz 0.9775 0.8988 0.03685 1099.1
73.5 pevap=0.2625 MPa; Pc°nd=1.028 MPa; /=31.05 Hz 0.9686 0.8567 0.03688 1117.9
84 pevap=0.262 MPa; Pc°nd=1.034 MPa; /=30.46 Hz 0.9610 0.8334 0.0369 1136.7
94.5 Making decisi°n °n a change in v°ltage fre-quency:/=39.8 Hz; Pevap=0.261 MPa; Pc°nd=1.1851 MPa 0.8291 0.9395 0.0373 1184
105 pevap=0.260 MPa; Pc°nd=1.325 MPa; /=37.52 Hz 0.7945 0.7685 0.0376 1260.3
115.5 Making decisi°n °n a change in v°ltage frequency: /=44.36 Hz; pevap=0.258 MPa; Pc°nd=1.3391 MPa 0.7841 0.9297 0.0378 1317.9
126 Pevap=0.254 MPa; Pc°nd=1.4298 MPa; /=30.46 Hz 0.7285 0.8005 0.038 1406.9
136.5 Making decisi°n °n a change in v°ltage frequency: /=50 Hz; pevap=0.252 MPa; Pc°nd=1.5206 MPa 0.7262 0 0.038 1500
147 pevap=0.252 MPa; pc°nd=1.5206 MPa; /=50 Hz 0 0 0.0398 1500
157.5 pevap=0.252 MPa; Pc°nd=1.5206 MPa; /=50 Hz 0 0 0.0398 1500
161.6 pevap=0.252 MPa; Pc°nd=1.5206 MPa; /=50 Hz 0 0 0.0398 1500
Note: Gr is the refrigerant flow rate. motor of the compressor, rpm; pevap. condensation pressure, MPa; f is the est. f. is the constant estimated value
kg/s; n is the number of rotations of the electric is the evaporation pressure, MPa;pcond. is the voltage frequency, Hz; t is the time, s. Index: c. of the parameter of the first level of functioning
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.
7. 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 °C...(-1,5) oC. 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, Kp, K/, Kx; 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,
œndensatfon pressure, v°ltage frequency that are measured with the reference values. Determining °f summarizing functfonal informatfon provides an °pp°rtunity t° make the foltawing advancing decisfons: t° maintain a change in evap°rati°n pressure as for a change in refrigerant vap°r
fl°w rate for digital œntrol: t° maintain a change in evap-°rati°n pressure as for a change °f refrigerant vap°r fl°w rate and a change °f v°ltage frequency as for a change °f the number °f rotatfons °f the electric m°t°r °f the c°mpress°r for frequency œntrol.
References
1. Jiang S. Air-Source Heat Pump Systems // Handbook of Energy Systems in Green Buildings. 2017. P. 1-44. doi: https:// doi.org/10.1007/978-3-662-49088-4_2-1
2. Rees S. J. An introduction to ground-source heat pump technology // Advances in Ground-Source Heat Pump Systems. 2016. P. 1-25. doi: https://doi.org/10.1016/b978-0-08-100311-4.00001-7
3. Chaikovskaya E. Devising an energy saving technology for a biogas plant as a part of the cogeneration system // Eastern-European Journal of Enterprise Technologies. 2015. Vol. 3, Issue 8 (75). P. 44-49. doi: https://doi.org/10.15587/1729-4061.2015.44252
4. Development of a Spiral Type Heat Exchanger for Ground Source Heat Pump System / Suzuki M., Yoneyama K., Amemiya S., Oe M. // Energy Procedia. 2016. Vol. 96. P. 503-510. doi: https://doi.org/10.1016Zj.egypro.2016.09.091
5. Liu W., Xu M. Research Progress of Pile Heat Exchangers in Ground Source Heat Pump System // Procedia Engineering. 2017. Vol. 205. P. 3775-3781. doi: https://doi.org/10.1016/jj.proeng.2017.10.141
6. Chaikovskaya E. The development of the method of maintaning the soil discharge in the heat pump energy supply // Technology audit and production reserves. 2016. Vol. 4, Issue 1 (30). P. 33-39. doi: https://doi.org/10.15587/2312-8372.2016.74705
7. Dincer I., Rosen M. A., Ahmadi P. Modeling and Optimization of Heat Pump Systems // Optimization of Energy Systems. 2017. P. 183-198. doi: https://doi.org/10.1002/9781118894484.ch6
8. Underwood C. P. Heat pump modelling // Advances in Ground-Source Heat Pump Systems. 2016. P. 387-421. doi: https:// doi.org/10.1016/b978-0-08-100311-4.00014-5
9. Li Y., Yu J. Theoretical analysis on optimal configurations of heat exchanger and compressor in a two-stage compression air source heat pump system // Applied Thermal Engineering. 2016. Vol. 96. P. 682-689. doi: https://doi.org/10.1016/j.applthermaleng.2015.11.132
10. Matuska T., Sourek B., Sedlar J. Heat Pump System Performance with PV System Adapted Control // Proceedings of EuroSun2016. 2016. doi: https://doi.org/10.18086/eurosun.2016.08.06
11. Yan G., Bai T., Yu J. Energy and exergy efficiency analysis of solar driven ejector-compressor heat pump cycle // Solar Energy. 2016. Vol. 125. P. 243-255. doi: https://doi.org/10.1016/jj.solener.2015.12.021
12. A Predictive Model for Smart Control of a Domestic Heat Pump and Thermal Storage / Van Leeuwen R. P., Gebhardt I., de Wit J. B., Smit G. J. M. // Proceedings of the 5th International Conference on Smart Cities and Green ICT Systems. 2016. doi: https://doi.org/10.5220/0005762201360145
13. Chaikovskaya E. E. Coordination energy production and consumption based on intellectual control heat and mass transfer processes // XV Minsk International Heat and Mass Transfer Forum: Heat and mass transfer processes in the energy and equipment. Energy savings. Minsk, 2016. P. 1-12.
14. Chaikovska Y. Y. Intellectual systems for support of energy systems functioning at level of decision-making // Power engineering: economics, technique, ecology. 2017. Issue 3. P. 114-118. doi: https://doi.org/10.20535/1813-5420.3.2017.117377
-□ □-
Дослiджено вплив установки надротор-них елементiв у виглядi статорног втулки з повздовжтми прямими пазами над передвклю-ченим колесом на характеристики модельного шнекововiдцентрового ступеня. Проведено фiзичний експеримент з використанням пла-нування для виршення проблеми оптимiзацiг геометричних параметрiв статорног втулки з повздовжтми прямими пазами в багатофак-торнш задачi покращення кавтацшно-ерозш-них характеристик шнекововiдцентрового ступеня за допомогою надроторних елементiв. Було визначено спектр збуджуючих частот коливань дослiджуваного шнекововiдцентрового ступеня вiд кавтацшних процеыв для можливостi використання параметра стiйкостi до кавта-цшног ерозп у якостi параметра оптимiзацiг. Експериментальним шляхом визначено опти-мальт розмiри надроторног втулки до^джу-ваного шнекововiдцентрового ступеня: Z=32, Ь=14,1=20,12=20. Це дозволило покращити кав^ тацшно-ерозшш якостi шнекововiдцентрового ступеня без змти габаритних розмiрiв та не погiршити натрш та енергетичш характеристики. Проведено додатковий фiзичний експе-римент за допомогою альтернативного методу визначення кавтацшно-ерозшних якостей для тдтвердження отриманих у до^дженн результатiв завдяки застосуванню надроторних елементiв у шнекововiдцентровому ступет. Використання надроторних елементiв у складi шнекововiдцентрових ступетв зазвичай обме-жувалося лише потребами тдвищення кавта-цшних якостей ступеня. В рамках до^дження, що описано в данш статтi, запропоновано використання цього елементу для боротьби з нега-тивними на^дками кавтацшног ерозп. Була тдтверджена можлив^ть такого використання та розроблен науково-методичн рекомен-дацп щодо проектування надроторних еле-ментiв у складi шнекововiдцентрового ступеня. Впровадження удосконалених перших шнеко-вовiдцентрових ступетв з надроторними еле-ментами в iснуючi конструкцп вiдцентрових насоЫв дозволить збшьшити наробтку до вiд-мови, що актуально для в^х галузей промисло-востi, де використовуються вiдцентровi насоси
Ключовi слова: вiдцентровий насос, шнеко-вовiдцентровий ступть, надроторн елементи,
кавiтацiйно-ерозiйнi характеристики -□ □-
UDC 621.65
I DOI: 10.15587/1729-4061.2018.139392 |
IMPROVEMENT OF CAVITATION EROSION CHARACTERISTICS OF THE CENTRIFUGAL INDUCER STAGE WITH THE INDUCER
BUSH
P. Tkach
Head of Hydraulic Part Research Department* Е-mail: pavlencio@gmail.com A. Yashchenko Head of Dynamics and Vibration Monitoring Laboratory* E-mail: yaschenko@vniiaen.sumy.ua O. G u s a k PhD, Associate Professor Department of Applied Fluid Aeromechanics** E-mail: info@teset.sumdu.edu.ua S. Khovanskyy PhD, Associate Professor Department of Applied Fluid Aeromechanics** E-mail: serg_83@ukr.net V. Panchenko Senior Lecturer
Department of Applied Fluid Aeromechanics ** E-mail: pan_va@ukr.net I. G rec h ka
PhD, Associate Professor Department of theory and computer aided design of mechanisms and machines National Technical University «Kharkiv Polytechnic Institute» Kyrpychova str., 2, Kharkiv, Ukraine, 61002 *JSC «VNIIAEN» "Research and Design Institute for Atomic and Power Pumpbuilding" 2-nd Zaliznychna str., 2, Sumy, Ukraine, 40003 **Sumy State University Rimskoho-Korsakova str., 2, Sumy, Ukraine, 40007
1. Introduction
Today one of the main problems is a lack of energy resources. Therefore, one of the main tasks of Ukraine and the world is to reduce inefficient energy consumption significantly. It is
well known that pumps are used almost in all industries, and according to various estimates consume from 20 % to 25 % of all electricity produced in the world, and in some industries, this value can range to 50 % [1]. One indicator of the pump energy efficiency is the cost of its life cycle, the analysis of
©
