Diagnostics of a Domestic Hot Water Storage Tank under Operating Conditions

: This paper presents a new method for the diagnostics of a hot water storage tank under operating conditions. Depending on the operating point of the tank, the method enables determination of thermal conductivity coefﬁcients of the coil heat exchanger, which allows us to determine the intensity of heat transfer between the transfer medium and water in the tank as well as of tank walls, which consequently enables determination of heat losses to the environment. Furthermore, the dynamic properties of the tank may also be determined by applying this method. The advantage of this method is possibility of analyzing changes in the material constants of the coil heat exchanger, tank walls, and dynamic properties of the tank as a function of mass ﬂow of the medium supplying the coil heat exchanger. The possibility of determining coefﬁcients of thermal conductivity as well as the inertia of tank and exchanger, based on temperature measurements acquired in operating conditions is a novelty in this paper. Knowing the variability of material constants and of dynamic properties of the tank as a function of medium ﬂow rate allows multicriteria optimization to be performed which, with a conventional design of the tank, results in a reduction of up to 10% in the time taken to prepare domestic hot water. by the relationship (39) with the determination coefﬁcient R 2 = 0.925. The distribution of the deviations proves the repeatability of the results obtained.


Introduction
One of the most important design parameters of a domestic hot water storage tank (DHWST) is the heat exchanger, which is installed inside the tank. The design parameters of the heat exchanger, such as heat transfer surface area in relation to the tank's volume and the material from which it is made, determine its power and heat transfer rate, which affects the time it takes to prepare hot water. The time for preparing hot water is also determined by the operating parameters, among others, temperature and volumetric flow of the heat transfer medium [1][2][3]. In practice, the hot water preparation time in the tank is expected to be as short as possible, which, irrespective of the number of consumer points, ensures that the system is comfortable to operate. One of the methods to reduce the time for preparation of hot water is to introduce stratification in the storage tank, which determines its size and shape and the shape and structural materials of the heat exchanger inside the tank [4,5]. Furthermore, improved stratification can be achieved by insulating the inner wall surface with an insulating material or other low thermal conductivity material compatible with the stored liquid [6][7][8]. At present, a frequent design solution for stratification of water in the storage tank is a tank with a built-in thermosiphon heat exchanger [9,10]. Due to the special design in the form of a thermosiphon inside the heat exchanger, hot water preparation time is reduced by up to 65% compared to a storage tank fitted with a conventional heat exchanger. However, this solution is significantly more expensive than a storage tank with a conventional heat exchanger, therefore optimization of the design of the coil heat exchanger installed in a hot water storage tank still poses a current and very important research problem, as well as a challenge for engineers [11][12][13]. In practice, domestic hot water storage tanks are varied in terms of the heat exchanger geometry and material in order to improve the heat exchange intensity and consequently reduce the time taken to prepare domestic hot water. For this reason, domestic hot water storage tanks are designed to be operated with a specific heat source, which determines the design, the developed heat transfer surface area of the exchanger and the materials used for its fabrication. It is assumed that in the case of condensing gas or oil boilers, manufacturers recommend the use of storage tanks with a built-in heat exchanger with the developed surface of 1 m 2 per 1 kW of the heating capacity of the device. In the case of a solar heating system, the developed heat transfer surface of the heat exchanger should be larger. It is assumed that per 1 kW of solar system heating capacity, the developed heat transfer surface should be 0.6 m 2 . In the case of a compressor heat pump, the developed heat exchange surface area of the heat exchanger should be even greater at 0.25 m 2 per 1 kW of heating capacity of the device [14]. Today, computational fluid dynamics modeling (CFD) is more and more often being used to design and analyze the thermal states of domestic hot water storage tanks. CFD analysis makes it possible to assess the impact of storage tank design parameters such as heat exchanger configuration [12], heat exchanger and tank wall material constants, operating conditions, and tank operating point [15,16] on the storage tank thermal performance. In practice, due to the variability of the operating point of the storage tank under operating conditions, at the design stage, even using CFD analyses, it is very difficult to predict how and within what limits the values of the material constants will change. The authors' operating experience shows that in practice, storage tanks are often improperly selected for the heat source they are operated with, which translates into a significant increase in hot water preparation time. This is particularly evident in systems where compressor heat pumps are used as heat sources. Therefore, it becomes important to determine the material constants of the key structural components, including the heat exchanger, under the operating conditions and their dynamic properties, which will make it possible to analyze the impact of changing the operating point of the heat exchanger on its operating parameters and allow multicriteria optimization to be carried out in order to improve the performance of the tank. In this regard, the method for diagnostics of hot water storage tank under operating conditions developed by the authors is a novel approach, enabling improvement of thermal properties of the tank.

Methodology
The most important operating parameters of the hot water tank are inertia of the tank and heat exchange surface of the exchanger. The intensity of the heat exchange process results from heat exchange surface and the material from which the exchanger is made, which is characterized by the heat conductivity coefficient expressed in [W/m 2 K]. The value of the heat conductivity coefficient is strongly related to operating parameters of the tank, including water temperature within the tank, temperature of the medium supplying the exchanger, and the volumetric flow. These parameters characterize the so-called operating point of the tank, which is dependent on them. The values of these parameters vary under operating conditions, therefore it is very important to determine the heat conductivity coefficient of the exchanger at a given operating point as a function of the volumetric flow of the medium supplying this exchanger. For this reason, the parameters mentioned above were selected for analysis, as described below.
The proposed diagnostics method is based on the comparison of two models: an analytical model, in the form of a system of differential equations, developed on the basis of an energy balance linking the characteristic parameters of the storage tank (thermal conductivity coefficients, tank volume, etc.) with an operational model, in the form of differential equations, developed on the basis of operational data using parametric identification. A prerequisite for the reliability of this method is the same order of analytical and parametric model. The System Identification library of the Matlab package was used to develop an operational model based on the operational data, therefore, the paper will explain the part of the methodology concerning the analytical model. The operational model of the storage tank based on energy balance was developed assuming that the water temperature in the tank is distributed uniformly over the entire height of the coil heat exchanger. The model therefore does not take into account possible stratification that may occur in the storage tank. However, in modern heating systems, due to the requirements of the users, internal circulation is ensured to compensate for stratification in the storage tank.
Assuming that the modeled hot water storage tank is a load for a solar collector array and is isolated from the solar segment by a coil heat exchanger placed inside the storage tank, the amount of heat supplied by the collector array to the tank can be determined from Equation (1).
Assuming that the distribution of water temperature in the storage tank over the height of the heat exchanger is uniform, the amount of heat received by the water in the storage tank can be calculated from Equation (2).
In a steady state, the amount of heat supplied by the heat exchanger is equal to the amount of heat received by the water in the tank. Equation (3) is, therefore, true.
Therefore, by substituting Equations (1) and (2) into Equation (3), it is possible to determine the outlet temperature of the medium from the heat exchanger according to Equation (4).
Assuming that the recorded curves of the inlet and outlet temperatures of the medium are dependent on the temperature of the water in the storage tank, Equation (4) can be simplified to the form 5.
According to Equation (5), the variation of the outlet temperature of the medium from the heat exchanger depends on the value of the U LW coefficient, the inlet temperature, the flow rate, and the specific heat of the medium in the system.
The temperature rise in a domestic hot water storage tank depends not only on the thermal capacity of the storage tank but also on the dynamic properties of the heat exchanger inside the storage tank. The dynamic properties of the heat exchanger can be determined from the differential equation.
By transforming Equation (6) into the form (7) and completing the Laplace transform, Equation (8) is obtained, which enables determination of the transition state of the outlet temperature of the medium from the coil.
The model of the coil heat exchanger developed on the basis of the energy balance consists of two paths: the main path related to the medium inlet temperature, described in Equation (9) by the transfer function G 1 (s), and the disturbance path related to the storage tank temperature, described in Equation (10) by the transfer function G 2 (s). The operational model of the heat exchanger can be represented by a block diagram (Figure 1). The model of the coil heat exchanger developed on the basis of the energy bala consists of two paths: the main path related to the medium inlet temperature, descri in Equation (9) by the transfer function G1(s), and the disturbance path related to the s age tank temperature, described in Equation (10)  The individual constituent paths of the model are described with the transfer fu tion, appropriate for a first-order inertial object. The gain factors of the individual p can be determined from Equations (11) and (12), and the equivalent time constant, wh is identical for both paths of the model, can be determined from Equation (13).   Equations (11) and (12) show that the gain factors of both model paths depend on flow and specific heat of the medium and the ULW coefficient. The equivalent time cons is related to the effective heat capacity of the heat exchanger and the flow rate and spe heat of the medium. The individual constituent paths of the model are described with the transfer function, appropriate for a first-order inertial object. The gain factors of the individual paths can be determined from Equations (11) and (12), and the equivalent time constant, which is identical for both paths of the model, can be determined from Equation (13).
Equations (11) and (12) show that the gain factors of both model paths depend on the flow and specific heat of the medium and the U LW coefficient. The equivalent time constant is related to the effective heat capacity of the heat exchanger and the flow rate and specific heat of the medium.
Assuming that there is no stratification in the storage tank over the height of the heat exchanger and that the amount of heat transferred by the coil is equal to the amount of heat received by the water in the storage tank, the energy balance of the tank can be described using Equation (14).
The temperature of the water in the storage tank can be calculated by transforming Equation (14) into the form: According to Equation (15), the water temperature rise in the storage tank depends on the inlet temperature of the medium to the coil and the specific heat and flow rate of the medium.
The water temperature rise in the storage tank depends on the amount of heat transferred through the heat exchanger and the heat losses through the tank walls. Knowing the amount of heat transferred by the heat exchanger into the storage tank and the heat losses from the tank to the environment, the temperature rise in the tank can be estimated using the Equation: (mc) ezas dT zas By transforming Equation (16) and rearranging it, we can obtain the Equation Equation (18) shows that the storage tank model developed on the basis of the energy balance consists of two paths: the main path described with Equation (19) by the transfer function G 1 (s), related to the medium inlet temperature T in to the coil, and the disturbance path described with Equation (20) by the transfer function G 2 (s) and related to the ambient temperature T amb . The operational model of the storage tank can be represented by a block diagram ( Figure 2).   According to the diagram in Figure 3, the DHW storage tank consists of two interde pendent objects: the coil heat exchanger and the water tank.
The inlet temperature of the medium to the heat exchanger and the ambient temper ature were taken as input signals, and the outlet temperature of the medium from th The transfer functions of both paths are appropriate for first-order inertial objects. The transfer functions describing the dynamic properties of the individual paths make it possible to determine the gain factors and the equivalent time constants. The main path gain factor can be determined from Equation (21), while the disturbance path gain factor can be determined from Equation (22). The equivalent time constant is defined by Equation (23).
The holistic model of the DHW storage tank can be represented with the following block diagram ( Figure 3).   According to the diagram in Figure 3, the DHW storage tank consists of two interdependent objects: the coil heat exchanger and the water tank.
The inlet temperature of the medium to the heat exchanger and the ambient temperature were taken as input signals, and the outlet temperature of the medium from the According to the diagram in Figure 3, the DHW storage tank consists of two interdependent objects: the coil heat exchanger and the water tank.
The inlet temperature of the medium to the heat exchanger and the ambient temperature were taken as input signals, and the outlet temperature of the medium from the exchanger was taken as the output signal. The water temperature in the heat exchanger is the variable that links the two parts of the storage tank model. According to the diagram of the storage tank shown in Figure 3, the parametric identification process would need to develop two independent models of the individual components of the storage tank, linked together by the water temperature in the storage tank. Consequently, in terms of its structure, the model would become complex. Assuming that the recorded curves of the inlet and outlet temperatures of the medium from the exchanger and the water temperature in the storage tank, as a function of time, are interdependent on each other and on the ambient temperature, the disturbance path can be neglected in both parts of the storage tank model. The operational model of the DHW storage tank is then simplified to the form represented with the block diagram in Figure 4. The dynamic properties of the main path of the operational model are described by the transfer function (24) and of the disturbance path-by the transfer function (25), the gain factors and the equivalent time constants of which can be estimated using Equations (26)-(29).
linked together by the water temperature in the storage tank. Consequently, in te its structure, the model would become complex. Assuming that the recorded curve inlet and outlet temperatures of the medium from the exchanger and the water te ture in the storage tank, as a function of time, are interdependent on each other the ambient temperature, the disturbance path can be neglected in both parts of th age tank model. The operational model of the DHW storage tank is then simplified form represented with the block diagram in Figure 4. The dynamic properties of th path of the operational model are described by the transfer function (24) and of t turbance path -by the transfer function (25)  The ULW and ULZ coefficients in the individual transfer functions of the paths model parts correspond to the heat loss coefficients. The ULW coefficient correspo the heat losses of the coil heat exchanger and the ULZ coefficient corresponds to t losses of the storage tank. The energy balance shows that the power supplied by t The U LW and U LZ coefficients in the individual transfer functions of the paths of the model parts correspond to the heat loss coefficients. The U LW coefficient corresponds to the heat losses of the coil heat exchanger and the U LZ coefficient corresponds to the heat losses of the storage tank. The energy balance shows that the power supplied by the heat exchanger should be equal to the power transferred to the storage tank. Therefore, the gain factors of both paths of the operational model of the D HW storage tank should be 1. The heat losses from the coil heat exchanger-U LW should be calculated with Equation (26). The power losses of the tank-U LZ should be calculated with Equation (28), and the thermal capacity of the tank with Equation (29). The heat transfer intensity of the heat exchanger to the storage tank, that is the exchanger thermal conductivity coefficient λ W , should be calculated using Equation (30).

Storage Tank Model Based on Operational Data
The diagnostics of a storage tank is illustrated by the operating data recorded during the operation of the solar heating system ( Figure 5). The system consisted of two flat solar collectors with a unit absorber area of 3 m 2 . The collectors were thermally loaded with a 100 dm 3 storage tank fitted with a coil heat exchanger with a volume of 2.3 m 3 and a developed heat transfer surface area of 0.75 m 2 ( Figure 6). The exchanger's declared power, according to the manufacturer, is 14 kW at the flow rate of the heating medium of 2.5 m 3 /h and the temperature of the heating medium of 70 • C, the temperature of feeding water of 10 • C, and the temperature of domestic hot water of 45 • C. The medium flow rate was forced by a pump with stepless adjustment between 2 and 9 l/min, and it was measured using a flow meter generating a voltage signal corresponding to the flow rate, in the range of 0-10 V DC. The temperature was measured using four-wire PT1000 sensors, which made it possible to compensate for the resistance of the wires. The system control algorithm was implemented in a PLC which used a TCP/IP protocol for communication with a computer and the installed SCADA LBX software (ver. 2.60 LAB-EL, Reguły Poland, 2020). This solution made it possible to monitor the system and record the measurement data. The measurement data were sampled and recorded with a 1-s interval.
(26). The power losses of the tank -ULZ should be calculated with Equation (28), and the thermal capacity of the tank with Equation (29). The heat transfer intensity of the hea exchanger to the storage tank, that is the exchanger thermal conductivity coefficient λW should be calculated using Equation (30).

Storage Tank Model based on Operational Data
The diagnostics of a storage tank is illustrated by the operating data recorded during the operation of the solar heating system ( Figure 5). The system consisted of two flat sola collectors with a unit absorber area of 3 m 2 . The collectors were thermally loaded with a 100 dm 3 storage tank fitted with a coil heat exchanger with a volume of 2.3 m 3 and a de veloped heat transfer surface area of 0.75 m 2 ( Figure 6). The exchanger's declared power according to the manufacturer, is 14 kW at the flow rate of the heating medium of 2.5 m 3 /h and the temperature of the heating medium of 70 °C, the temperature of feeding water o 10 °C, and the temperature of domestic hot water of 45 °C. The medium flow rate wa forced by a pump with stepless adjustment between 2 and 9 l/min, and it was measured using a flow meter generating a voltage signal corresponding to the flow rate, in the range of 0-10 V DC. The temperature was measured using four-wire PT1000 sensors, which made it possible to compensate for the resistance of the wires. The system control algo rithm was implemented in a PLC which used a TCP/IP protocol for communication with a computer and the installed SCADA LBX software (ver. 2.60 LAB-EL, Reguły Poland 2020). This solution made it possible to monitor the system and record the measuremen data. The measurement data were sampled and recorded with a 1-second interval.  The operating data recorded at various operating points of the tank are to be used to diagnose the storage tank. In case of the solar heating system under consideration, the operating point will vary depending on prevailing sunlight conditions and flow of the working medium. The impact of working medium mass flow in solar installation on heat transfer coefficients as well as on thermal inertia of DHW storage tank and coil heat exchanger were analyzed on data recorded during 15 days of the installation operation. During each experiment, the mass flow of the working medium was constant, and the values of solar radiation intensity were high and stable over time. The instantaneous variability of the intensity of solar radiation during the experiment had a negligible impact on the dynamic properties of the reservoir (very high reservoir time constant) and the thermal conductivity coefficients. All experiments were performed for three mass flows of the working medium through the coil exchanger: the minimum one was 0.038 ± 0.02 [kg/s], the average one was 0.09 ± 002 [kg/s] and the maximum one was 0.13 ± 02 [kg/s]. Examples of operational data from one of these 15 days are presented in Figure 7. During the experiment that day, the medium flow through the heat exchanger was 0.132 [l/s]. On this basis, the proposed diagnostic method was discussed, but it should be emphasized that in the case of operational data recorded for the remaining 14 days, the procedure was identical.
In the identification process, the DHW storage tank was treated according to the analytical model as a single-input dual-output object (Figure 4). The inlet temperature of the medium to the coil, which is also the outlet temperature of the medium from the collector array, was taken as the input signal. The outlet temperature of the medium from the coil, which is also the inlet temperature of the medium to the collector array, and the temperature of the water in the storage tank were taken as the output signals. For temperature measurements, class A Pt1000 sensors were used for which the measurement uncertainty in the temperature range -30 °C -300 °C was (± 0.15 + 0.002|t|). The operational model developed on the basis of such data enables the analysis of the dynamics of the coil heat The operating data recorded at various operating points of the tank are to be used to diagnose the storage tank. In case of the solar heating system under consideration, the operating point will vary depending on prevailing sunlight conditions and flow of the working medium. The impact of working medium mass flow in solar installation on heat transfer coefficients as well as on thermal inertia of DHW storage tank and coil heat exchanger were analyzed on data recorded during 15 days of the installation operation. During each experiment, the mass flow of the working medium was constant, and the values of solar radiation intensity were high and stable over time. The instantaneous variability of the intensity of solar radiation during the experiment had a negligible impact on the dynamic properties of the reservoir (very high reservoir time constant) and the thermal conductivity coefficients. All experiments were performed for three mass flows of the working medium through the coil exchanger: the minimum one was 0.038 ± 0.02 [kg/s], the average one was 0.09 ± 0.02 [kg/s] and the maximum one was 0.13 ± 0.2 [kg/s]. Examples of operational data from one of these 15 days are presented in Figure 7. During the experiment that day, the medium flow through the heat exchanger was 0.132 [l/s]. On this basis, the proposed diagnostic method was discussed, but it should be emphasized that in the case of operational data recorded for the remaining 14 days, the procedure was identical.
In the identification process, the DHW storage tank was treated according to the analytical model as a single-input dual-output object (Figure 4). The inlet temperature of the medium to the coil, which is also the outlet temperature of the medium from the collector array, was taken as the input signal. The outlet temperature of the medium from the coil, which is also the inlet temperature of the medium to the collector array, and the temperature of the water in the storage tank were taken as the output signals. For temperature measurements, class A Pt1000 sensors were used for which the measurement uncertainty in the temperature range −30 • C-300 • C was (±0.15 + 0.002|t|). The operational model developed on the basis of such data enables the analysis of the dynamics of the coil heat exchanger-main path and the analysis of the dynamics of the water temperature in the storage tank-disturbance path.  In the parametric identification process, an operational model of the storage tank was developed based on the recorded operational data. The static and dynamic properties of the DHW storage tank were best represented by the ARX autoregressive model, described by a system of differential Equation (31), as estimated in the validation process ( Figure 8). The determination coefficient of the model was 0.99.
The determination coefficient of the model was 0.99.

Analysis of Material Constants and Static and Dynamic Properties of the Storage Tank
In order to determine the dynamic properties of the coil and the storage tank, the operational model in the form of a system of differential Equation (31) was transformed using the Laplace transform to the form of transfer functions (32) and (34), and in order to compare it with the analytical model, on the basis of the step response curves of the main path and the disturbance path, it was approximated to the form of first-order transfer functions, which correspond to Formulas (33) and (35). Figure 9 shows the step response of the main path of the operational model described by the transfer function (32).

Analysis of Material Constants and Static and Dynamic Properties of the Storage Tank
In order to determine the dynamic properties of the coil and the storage tank, the operational model in the form of a system of differential Equation (31) was transformed using the Laplace transform to the form of transfer functions (32) and (34), and in order to compare it with the analytical model, on the basis of the step response curves of the main path and the disturbance path, it was approximated to the form of first-order transfer functions, which correspond to Formulas (33) and (35). Figure 9 shows the step response of the main path of the operational model described by the transfer function (32).    Comparing the transfer functions of the two paths of the storage tank model, it should be noted that they have an identical characteristic equation. Discrepancies in the transient curves between the paths of the operational model of the storage tank will result from the position of the zeros (numerator radicals) relative to the poles (denominator radicals) of the transfer function in the plane of the complex variable. The position of the zeros relative to the poles depends on the mutual trajectories of the output signals relative to the input signal. Figure 10 shows the step response of the disturbance path of the operational model of the DHW storage tank. On the basis of the transient curves, it should be concluded that the model's disturbance path has a non-minimum-phase character. In the first phase, the output signal (storage tank temperature) trends in the opposite direction to the intended one, and then aperiodically reaches a steady state. It is also worth noting that the gain factor of the disturbance path was 1.04 which, despite the good insulation of the storage tank, indicates power loss. In the case under consideration, the power losses from the storage tank calculated from the relationship (21)  from the position of the zeros (numerator radicals) relative to the poles (denominator radicals) of the transfer function in the plane of the complex variable. The position of the zeros relative to the poles depends on the mutual trajectories of the output signals relative to the input signal. Figure 10 shows the step response of the disturbance path of the operational model of the DHW storage tank. On the basis of the transient curves, it should be concluded that the model's disturbance path has a non-minimum-phase character. In the first phase, the output signal (storage tank temperature) trends in the opposite direction to the intended one, and then aperiodically reaches a steady state. It is also worth noting that the gain factor of the disturbance path was 1.04 which, despite the good insulation of the storage tank, indicates power loss. In the case under consideration, the power losses from the storage tank calculated from the relationship (21)

Analysis of the Variation of Material Constants and Storage Tank Operating Parameters as a Function of the Medium Flow
The gain factors and the equivalent time constant, and thus the material constants of both paths of the operational model of the DHW storage tank, described by Equations (24)

Analysis of the Variation of Material Constants and Storage Tank Operating Parameters as a Function of the Medium Flow
The gain factors and the equivalent time constant, and thus the material constants of both paths of the operational model of the DHW storage tank, described by Equations (24) and (25), are strongly dependent on the performance of the circulator pump. Having operational data recorded with different circulator pump output, the impact of flow rate on the static and dynamic properties of the DHW storage tank can be estimated. Figure 11 shows the impact of the circulator pump output   Figure 11. Impact of the circulator pump output on the value of the gain factor of the main path. Figure 12 shows the impact of the circulator pump output on the heat transfer r the heat exchanger. As in the case of the gain factor of the main path of the storage model, an increase in the output of the circulator pump causes an increase in the transfer rate of the heat exchanger. The relationship between the pump output an heat transfer intensity is linear and, in the case under consideration, can be described Equation (37).
= 7410 * ̇− 52.83  transfer rate of the heat exchanger. The relationship between the pump output an heat transfer intensity is linear and, in the case under consideration, can be described Equation (37).
= 7410 * ̇− 52.83 Analyzing the distribution of deviations for individual values of the heat tra coefficient, it can be observed that the highest error of 2.17% occurred for a flow ra 0.137 [l/s]. This demonstrates the repeatability of the results. Figure 13 shows the im of the circulator pump output on the value of the equivalent time constant of the exchanger. The graph shows that the higher the output of the circulator pump, the l The distribution of the deviations shows that the lowest repeatability of the results was obtained for the maximum flow rate of 0.137 [l/s], for which the value of the equivalent time constant was the lowest. It is likely that in addition to the output of the circulator pump, the value of the equivalent time constant was affected by the dynamic properties of the collector array, which were determined, among other things, by the variation in solar irradiance.
An increase in the circulator pump output translates into an increase in the heat transfer rate of the exchanger, which affects the dynamic properties of the storage tank. An increase in the flow rate of the medium through the heat exchanger results in a decrease in the value of the equivalent time constant (Figure 14). In the solar system under consideration, increasing the output of the circulator pump from the minimum to the maximum reduced the time constant of the storage tank by an average of 2300 [s]. This means that the water temperature rise in the storage was be faster at maximum output of the circulator pump. The change in the equivalent time constant of the storage tank in the system under consideration as a function of the circulator pump output can be described by the relationship (39) with the determination coefficient R 2 = 0.925. The distribution of the deviations proves the repeatability of the results obtained. The distribution of the deviations shows that the lowest repeatability of the res was obtained for the maximum flow rate of 0.137 [l/s], for which the value of the equ lent time constant was the lowest. It is likely that in addition to the output of the circul pump, the value of the equivalent time constant was affected by the dynamic prope of the collector array, which were determined, among other things, by the variatio solar irradiance. was obtained for the maximum flow rate of 0.137 [l/s], for which the value of the equ lent time constant was the lowest. It is likely that in addition to the output of the circul pump, the value of the equivalent time constant was affected by the dynamic prope of the collector array, which were determined, among other things, by the variatio solar irradiance.

Analysis of the Impact of the Circulator Pump Output on the Dynamic Propertie the Storage Tank
According to the above analyses and from the point of view of operation, the namic properties of the DHW storage tank were most affected by changing the outp the circulator pump. It is therefore important to determine the impact of the change in

Analysis of the Impact of the Circulator Pump Output on the Dynamic Propertie the Storage Tank
According to the above analyses and from the point of view of operation, the namic properties of the DHW storage tank were most affected by changing the outp the circulator pump. It is therefore important to determine the impact of the change in circulator pump output on the dynamic properties of the storage tank. Estimates ca made using the operational model of the storage tank defined by Equations (24) and Figure 16 shows

Analysis of the Impact of the Circulator Pump Output on the Dynamic Properties of the Storage Tank
According to the above analyses and from the point of view of operation, the dynamic properties of the DHW storage tank were most affected by changing the output of the circulator pump. It is therefore important to determine the impact of the change in the circulator pump output on the dynamic properties of the storage tank. Estimates can be made using the operational model of the storage tank defined by Equations (24) and (25). Figure 16 shows the step response of the water temperature in the storage tank with the circulator pump output varying between 0.039 [l/s] and 0.139 [l/s]. The longest time constant of the tank corresponded to the minimum flow rate of 0.039 [l/s]. The water in the storage tank heated up fastest at the maximum circulator pump output of 0.139 [l/s], which means that the thermal conductivity coefficient λ W for this flow rate will reach its maximum value. Analyzing the transient curves of the water temperature in the storage tank, it seems that the increase in the flow rate of the medium above 0.105 [l/s] does not affect the shortening of the value of the equivalent time constant, and thus, according to the relationship (30), the increase in the thermal conductivity coefficient of the exchanger λ W . In the system under consideration, the flow rate of the medium through the heat exchanger should not be greater than 0.105 [l/s]. The output of the circulator pump above 0.105 [l/s] will result in higher electricity consumption by the circulator pump, increasing the operating costs. maximum value. Analyzing the transient curves of the water temperature in the storage tank, it seems that the increase in the flow rate of the medium above 0.105 [l/s] does not affect the shortening of the value of the equivalent time constant, and thus, according to the relationship (30), the increase in the thermal conductivity coefficient of the exchanger λW. In the system under consideration, the flow rate of the medium through the heat exchanger should not be greater than 0.105 [l/s]. The output of the circulator pump above 0.105 [l/s] will result in higher electricity consumption by the circulator pump, increasing the operating costs.

Conclusion
The developed diagnostics method for a domestic hot water storage tank enables the analysis of changes in the material constants of the coil heat exchanger and the dynamic properties of the storage tank under operating conditions. The analytical model was derived from the energy balance of the DHW storage tank, which is subjected to the laws of

Conclusion
The developed diagnostics method for a domestic hot water storage tank enables the analysis of changes in the material constants of the coil heat exchanger and the dynamic properties of the storage tank under operating conditions. The analytical model was derived from the energy balance of the DHW storage tank, which is subjected to the laws of heat exchange. The discrete model was created for operational data from the measurements of the exchanger operating temperatures. The method consists in comparing both models (analog and discrete), so it is correct provided that both models are of the same order as described in the methodology. Therefore, the parameters of the DHW storage tank determined by this method are characteristic for a given work point. Having the operational data recorded at various operating points, it is possible to draw conclusions on the changes in the dynamic properties of the storage tank and exchanger as well as of their coefficients of thermal conductivity under operating conditions.
In the case studied, it was shown that the key operating parameter affecting the intensity of the heat transfer process and the dynamic properties of the domestic hot water storage tank was the flow rate of the working medium.