Inverse problem with several anomalies. Heating and DHW installation We consider again the installation of Example E.9.1, to which the two anomalies of Example E.9.2 are incorporated. As we have said, one of them is in the RS radiator system, with a 10% decrease in its efficiency and the other is in the HX heat exchanger, with a 35% decrease in its heat transfer coefficient. The objective of this example is to solve the discrimination problem of intrinsic anomalies by means of the method of characteristic equations in parallel with the method of malfunctions and dysfunctions.
Since 13 components have been considered, 13 characteristic curves must be defined. For achieving this, the different Types existing in TRNSYS v 17 were used, the characteristic curves of two of these components having been shown in the two previous examples (Example E.9.4 , Example E.9.3). The simulation of the installation was carried out hour by hour, with the demand profiles that have been represented in Fig. E.9.2. The cumulative values of the exergy of the flows in the reference conditions and for the state with the two anomalies introduced in the equipment described above are presented in Table E.9.6. In order not to include external induced anomalies (such as control inter-vention, changes in external conditions, etc.), the reference state coincides with the free condition.
When there are several anomalies in an installation, the anomaly in one of the i-th components leads to the unit consumption variation in the component itself \Delta k_{i,i n t} (intrinsic malfunction) and in the rest of the j-th equipment (j = 1,n, j ≠ i) it leads to induced malfunctions \Delta k_{j,i nd} as well as a variation in its products \Delta P_{j} (dysfunctions). The objective is to discriminate between these induced malfunctions and dysfunctions, in order to evaluate the increase in the resources consumption due to the intrinsic mal-function of the i-th component. By using the MD method, the extra consumption of resources associated with the increase in the products of the j -th equipment can be obtained through the dysfunctions D F_{ji} but additional information is needed to eval-uate the effect of the other induced malfunctions.
Complementing the MD diagnosis with that of the CC diagnosis, the component with the most important intrinsic malfunction can be identified, that which causes the greatest impact on fuel. Once this intrinsic malfunction is eliminated, the analysis is carried out again by applying the MD method a second time. The impact on fuel calculated in the first analysis \Delta F_{T}^{1} less the impact on fuel made in the second \Delta F_{T}^{2} (when the main anomaly has been eliminated) is the fuel saving obtained when that anomaly has been eliminated, that is to say
On the other hand, this fuel saving is the sum of the intrinsic malfunction M F_{i,i n t}^{1^{st}} plus the effects induced in that first analysis \left(\sum_{j}D F_{j i}^{1^{st}}+\sum_{j}M F_{j i}^{i n d,1^{st}}\right) plus the total product variation of the plant \Delta P_{s}^{1,2} and the dysfunction that is generated between both states. In this way, induced malfunctions can be calculated. The method can be repeated as many times as intrinsic malfunctions exist and the inverse diagnosis problem solved, see Fig E.9.12. For more details of the method see Picallo et al. 2016 [E.1].
The values obtained are presented in Table E.9.8. This table refers to the first anal-ysis, which contains the two malfunctions, and then the second is presented, which has the most influential malfunction removed. In column M F_{i n t}^{1^{st}} the intrinsic malfunctions are represented, while column M F_{i nd}^{1^{st}} represents the values of the induced malfunc-tions, as the characteristic curves of the components are not represented by horizontal lines. The sum of both columns corresponds to the total malfunctions of each compo-nent, while the last column represents the dysfunctions. In view of the table we can make the following comments:
Since we identify two components with intrinsic malfunctions, the next step is to eliminate the intrinsic malfunction of greatest weight (that of RS). Once that intrinsic malfunction has been eliminated, and this component is brought back to the reference conditions, we perform the simulation again, in order to evaluate the decrease in the impact on fuel between the first and the second simulations.
Table E.9.9 shows the results obtained by applying the method of characteristic equations to this second simulation. In view of the table we can make the following comments:
The general results corresponding to the anomalies which were deliberately caused are collected in Table E.9.10. The rows M F_{i n t},\sum_{}^{}M F_{i n d},\mathrm{and}\sum D F are the intrinsic and induced malfunctions and the dysfunctions of the components with anomalies; while row D F_{0}+\Delta P_{s} indicates the effects of the anomalies on the final production. Finally, D F_{anomaly} summarizes the impact on fuel of each anomaly.
Thus, the anomaly in RS generates an extra consumption of 6886 MJ of which 7599 MJ are due to the anomaly proper, and the remaining -714 MJ are due to the decrease in the final product. On the other hand, the anomaly in HX generates an over-consumption of -590 MJ (underconsumption in this case), of which 277 MJ are due to the anomaly in the heat exchanger and the remaining -867 MJ are explained by the reduction in the DHW production.
Table E.9.6 Accumulated values of the flows exergy during the heating period.
[GJ] | B1 | B2 | B3 | B4 | B5 | B6 | B7 | B8 | B9 | B10 | B11 | |
FREE | 122.9 | 100.1 | 372.3 | 351.8 | 192.0 | 180.3 | 153.3 | 38.9 | 29.8 | 182.9 | 169.2 | |
REAL | 122.9 | 99.2 | 369.8 | 348.4 | 190.7 | 179.1 | 151.9 | 38.8 | 29.8 | 181.7 | 166.7 | |
B12 | B13 | B14 | B15 | B16 | B17 | B18 | B19 | B20 | B21 | B22 | B23 | B24 |
57.6 | 122.7 | 111.6 | 37.2 | 28.3 | 0.2 | 6.5 | 2.3 | 149.1 | 0.04 | 0.2 | 5.8 | 0.03 |
57.7 | 121.4 | 109.2 | 36.2 | 27.5 | 0.2 | 6.4 | 2.3 | 155.3 | 0.05 | 0.2 | 5.6 | 0.03 |
Table E.9.8 MF, DF values deduced using the characteristic curves method.
Characteristic curves | ||||
M F_{i n t}^{1^{st}} | M F_{i nd}^{1^{st}} | D F^{1^{st}} | ||
① | CB | _ | -1214 | 6467 |
② | HC | _ | -450 | 500 |
③ | D1 | _ | _ | _ |
④ | V1 | _ | _ | _ |
⑤ | M1 | _ | _ | _ |
⑥ | HX | 323 | -117 | -22 |
⑦ | V2 | _ | _ | _ |
⑧ | M2 | _ | -136 | 42 |
⑨ | M3 | _ | _ | _ |
⑩ | RS | 1212 | -119 | _ |
⑪ | T | _ | -40 | -6 |
⑫ | V3 | _ | -1 | -15 |
⑬ | D2 | _ | _ | _ |
Table E.9.9 MF, DF and Ps deduced using the method of characteristic curves in the second simulation.
Characteristic curves | MF and DF diagnosis | |||||
M F_{i n t}^{2^{nd}} | M F_{i nd}^{2^{nd}} | D F_{0}^{2^{nd}} | D F_{}^{2^{nd}} | D P_{s}^{2^{nd}} | ||
① | CB | _ | -2048 | -754 | 2197 | _ |
② | HC | _ | -143 | 1 | 82 | _ |
③ | D1 | _ | _ | _ | _ | _ |
④ | V1 | _ | _ | _ | _ | _ |
⑤ | M1 | _ | _ | _ | _ | _ |
⑥ | HX | 317 | -118 | -6 | -9 | _ |
⑦ | V2 | _ | _ | _ | _ | _ |
⑧ | M2 | _ | -45 | -11 | 59 | _ |
⑨ | M3 | _ | _ | _ | _ | _ |
⑩ | RS | _ | 18 | _ | _ | _ |
⑪ | T | _ | -33 | -12 | _ | _ |
⑫ | V3 | _ | -1 | -10 | _ | -76 |
⑬ | D2 | _ | _ | _ | _ | _ |
Table E.9.10 General results of the diagnosis.
RS anomaly | HX anomaly | |
M F_{i n t} | 1212 | 317 |
\textstyle\sum M F_{i n d} | -695 | -1270 |
\sum D F | 7082 | 1230 |
D F_{0}+\Delta P_{s} | -714 | -867 |
\Delta F_{a n o m a l y} | 6886 | -590 |