next up previous
Next: 4.1 Measurement Selection Up: Measurement Selection and Diagnosability Previous: 3 Fault Isolation

4 The Liquid Sodium Cooling System

 

Our transient-based diagnosis scheme has been successfully applied to a number of different hydraulic systems (see http://air.vuse.vanderbilt.edu:8080/ for a working version of the system) [7]. To investigate scalability of the measurement selection algorithm we applied it to the model of the secondary cooling loop of a fast breeder reactor. The need for a qualitative approach to fault detection and isolation in this system is motivated by its high-order (six), nonlinearity, and the non availability of precise and real-time numerical simulation models. The precision of flow sensors is limited by signal noise, and achieving hardware redundancy by installing flow sensors is expensive.

Heat from the reactor core is transported to the turbine by a primary and secondary cooling system. Liquid sodium is pumped through an intermediate heat exchanger to transport heat from the primary cooling loop to the feed water loop by means of a superheater and evaporator vessel (Fig. 7). Pump losses are modeled by tex2html_wrap_inline1062 . The coil in the intermediate heat exchanger that accounts for flow momentum build-up is represented by a fluid inertia, tex2html_wrap_inline1064 . The two sodium vessels are capacitances, tex2html_wrap_inline1066 and tex2html_wrap_inline1068 . An overflow column, tex2html_wrap_inline1070 , maintains a desired sodium level in the main motor. All connecting pipes are modeled as resistances.

   figure289
Figure 7: Secondary sodium cooling loop.

The main motor driver (Fig. 8) is a synchronous, tex2html_wrap_inline1072 motor, and as an assumption, its electrical field is considered to be present as soon as it is turned on. Therefore, dynamic electrical effects are not modeled, and the electrical part of the motor system can be represented as a source of mechanical energy with a given torque/angular velocity characteristic. The inertia of the rotor and the mass of transmission gear is modeled by tex2html_wrap_inline1074 , and the transmission ratio between motor and pump by n. Pump losses in the fluid connection between the motor and pump are modeled by a dissipation element, tex2html_wrap_inline1062 , and the pump inertia is represented as tex2html_wrap_inline1080 . The model of a centrifugal pump can be derived using conservation of power and momentum [6]. The pump is represented by tex2html_wrap_inline1082 which describes a modulated gyrator with modulus tex2html_wrap_inline1084 . If the pump veins are not curved, b = 0.

   figure299
Figure 8: Pump driven by an ac motor.

The derivation of the causal relations of the sodium pump are based on a modulation factor g between input angular velocity, tex2html_wrap_inline1090 , and output flow rate, tex2html_wrap_inline1092 , tex2html_wrap_inline1094 . This factor is directly proportional to tex2html_wrap_inline1090 and inversely proportional to tex2html_wrap_inline1092 , tex2html_wrap_inline1100 . The dependency of g on tex2html_wrap_inline1090 and tex2html_wrap_inline1106 can be explicitly modeled by edges between these variables and the affected variables. In case of the dynamic behavior, the affected variables are input torque, tex2html_wrap_inline1108 , and output pressure, tex2html_wrap_inline1110 , and the corresponding edges are added to the causal graph (Fig. 9).

   figure305
Figure 9: Temporal causal graph.

The dependency on system variables of the modulation factor results in nonlinear, quadratic, behavior tex2html_wrap_inline1112 , and, therefore, the relation on the edge between tex2html_wrap_inline1092 and tex2html_wrap_inline1108 is unknown, in general. A sensitivity analysis of this relation, shown in Fig. 10, reveals that depending on the values of tex2html_wrap_inline1090 and tex2html_wrap_inline1092 , the sensitivity of tex2html_wrap_inline1108 to tex2html_wrap_inline1092 is positive or negative. Given the nominal values of the steady state operation of the system, which is parameter dependent, the weight of tex2html_wrap_inline1126 can be determined as a positive (1) or negative (-1) influence. However, once a deviation occurs, tex2html_wrap_inline1090 and tex2html_wrap_inline1092 may differ from their nominal values and a different operating point may be reached. Since these new values are caused by failure, and, therefore, unknown, the influence may reverse and is unknown as well. Because this can only occur if tex2html_wrap_inline1108 is predicted to be high based on the proportional influence (-1 or 1), only a predicted decrease in tex2html_wrap_inline1108 is unambiguous, and, therefore, propagated. A predicted increase in tex2html_wrap_inline1108 is propagated as unknown.

   figure312
Figure 10: Detailed sensitivity analysis of tex2html_wrap_inline1142 .




next up previous
Next: 4.1 Measurement Selection Up: Measurement Selection and Diagnosability Previous: 3 Fault Isolation

Pieter J. Mosterman
Mon Aug 18 15:29:41 CDT 1997