Response Analysis on Electrical Pulses under Severe Nuclear Accident Temperature Conditions Using an Abnormal Signal Simulation Analysis Module

1 Severe Accident and PHWR Safety Research Division, Korea Atomic Energy Research Institute, 150 Dukjin-dong, Yusong-gu, Daejeon 305-353, Republic of Korea 2 Integrated Safety Assessment Division, KAERI, 150 Dukjin-dong, Yuseong-gu, Daejeon 305-353, Republic of Korea 3 Nuclear Safety Research Headquarter, KAERI, 150 Dukjin-dong, Yuseong-gu, Daejeon 305-353, Republic of Korea 4 Department of Optical Communications Engineering, Gwangju University, 52 Hyuodeok-ro, Nam-gu, Gwangju 503-703, Republic of Korea 5 Department of Electrical Engineering, Korea University, 5ka, Anam-dong, Sungbuk-ku, Seoul 136-701, Republic of Korea


Introduction
To diagnose a severe nuclear accident, it is essential to determine a plant's status and to continuously monitor a plant's responses, expected from the actions taken by operators to mitigate a given accident [1].In addition, a correct interpretation of a plant's conditions during such an accident is of significant importance for a successful accident management [2].According to SECY 89-012, the United State Nuclear Regulatory Commission (U.S. NRC) requires that in a reasonable way instrumentation is identified as one of the key elements in utility accident management plants [3].Moreover, a review of a plant's accident management capabilities is regarded as a crucial element in achieving a regulatory closure for severe accident issues.During accidents, information and data from a plant's instruments are essential in assessing the status and response of that plant.There are some inherent uncertainties of the reliability of instrumentations expected while subjected to harsh environments occurring in severe nuclear accident conditions.Even for a severe accident condition, however, there are a few ways to obtain the relevant information and data.According to [4,5], for example, circuit simulation analysis and diagnosis methods are available to assess these instruments in detail even when they provide apparently abnormal readings.Then, some abnormal signal diagnostics and analysis could be made through modeling of an important circuit that is composed of passive elements (R-L-C), an input multifunction generator, and an output pulse measurement system.These simulations would be useful for investigating what the signal and circuit response characteristics would be, similar to the relevant symptoms that can be obtained under severe environmental conditions as a quality analysis method [5].
In this paper, a newly enhanced simulation module has been designed through the realization of a one body order system [6] via one program that links an abnormal signal simulation engine program (PSpice and Multi-SIM code) to a signal analysis tool program (LabVIEW code), which includes signal analysis functions from an important circuit model applicable under the conditions [7].In that case, the aforementioned abnormal signal simulation engine is used as a one order command which means a nonstop processing command to import the simulation data from an output of the engine code to an input of the main analysis tool program.The design of the above simulation module can be broken up into three main steps [8].The first step is a section of a comparative assessment of the instruments between the instrument specifications or EQ (environment qualification) requirement parameters and the environment parameters that can be obtained through a severe accident simulation code such as MAAP4 whose mean circumstance data is obtained through a simulation of various accident scenarios [9].From the results, if its comparative result is not satisfied, the second step is a section of the design of a circuit simulation to apply an abnormal signal simulation processing technique with an equivalent circuit, the third step is a section of a signal analysis of functions that have an additional 5 substeps.Consequently, three main functions exist and are characterized as an abnormal signal processing, an accident management, and an additional guide that are displayed on the initial screen of a simulation system ASSA module.Using the ASSA module, we could obtain the response characteristics from the output voltage levels of a pulse parameter according to a change of the resistance for each of R 2 , and R 3 , and the response characteristic from the output voltage level and time constants of the pulse parameter according to a change of the capacitance for each of C 1 , and C 2 .
As a result, when we can analyze the cause from an equivalent circuit, which consists of an R-C distribution parameter, from a response characteristic for a 4∼20 mA instrument during high temperature conditions of the severe accident conditions (SACs), it is also possible to detect and diagnose the abnormal signal patterns of a malfunctioning instrument.In addition, the pulse parameter estimation is operable to estimate a pulse parameter for a signal based on an analysis of the detected pulses and information pertaining to there being abnormal pulses in the signal as a quantitative analysis method.
When this simulation module is set up in a laboratory, we can obtain a variety of abnormal signals simulated under degraded conditions like higher temperature conditions of the SACs.A large variety of built-in models and a simulator adds flexibility to a simulation and includes the effect of harsh environments by using a special function to identify the components exceeding a manufactures' limits from a circuit modeling including the R-C passive elements.
We present the design steps with the diagnosis methods and the circuit simulation using a specific signal processing in Section 2, followed by a detailed description of work carried out during the past five years of our diagnostics and analysis work in Sections 3-5.Some concluding remarks are presented in Section 6 and a reference of papers resulting from the project is included.

Design Steps of Diagnosis Methods
and the Circuit Simulation with Specific Signal Processing

Design Steps of Diagnosis Methods for the Malfunctioning
Instruments.When an instrument which is providing information and data for severe accident management is apparently malfunctioning, a series of steps can be taken which include an indirect diagnosis of the instrument and a direct method of determining the value of the parameter.These methods are identified by five steps which are called operational aids.They are combinations of diagnostic actions and other means of measurement.These combinations can be applied when instrument readouts are suspect or faulty.The five types of aids are a diverse indication, parameter inference, a portable instrument, circuit diagnosis, and a portable circuit readout [4] (Figure 1).In this case, a circuit diagnosis from among the diagnostic evaluations belongs to action items for a malfunctioning instrument where a circuit diagnosis using a simulation tool code is possible.

Design
Steps of the Circuit Simulator.Figure 2 shows an initial screen consisting of functions for the ASSA module.It has three main functions which have not only an abnormal signal processing function, but also a severe accident management function based on EPRI TR-103412 and EPRI TR-102371 technical reports including a comparative assessment procedure for instruments [4,5], and an additional guide on the initial screen of an ASSA module.
Figure 3 shows a flowchart of the overall design steps for the simulator as an ASSA module system.The comparative assessment procedure for the instruments has the main information and data including the information needs for the instrument specifications, position information according to five areas in the containment building, and the circumstance parameters from the accident scenarios using the MAAP4 code. Figure 4(a) shows the three EQ conditions which include a normal condition, an abnormal condition, and an out of range condition.Figure 4   a scenario taken by MAAP4 code. Figure 5 shows a screen for a comparative assessment procedure as the first step for the ASSA module system.As a result, if the proper conditions are not satisfied by using the comparative assessment method, then proceed to a next step as the section of LabVIEW Tool code.As the next step, we need decision making for signal conditions which means three kinds of signal conditions, of which the first is normal condition signals, the second is abnormal condition signals, and the third is out of control condition signals.
In the case of abnormal condition signals, we are interested in the information on the patterns of the abnormal signals and the database for the signal patterns that also could be obtained by the engine code, which then need to be processed by any other five substeps from the analysis tool code, and finally we are able to obtain corrective signals as a final step, but not yet here.There are two codes to be used for the simulation, the designed simulator which is composed of the LabVIEW code as a main analysis tool code, and the PSpice and MultiSIM code as a simulation engine code.For these abnormal signals simulations, a new simulator through an analysis of a typical equivalent circuit modeling has been designed, which is composed of an analysis tool code as a main body and a simulation engine code whose  systems are imported to link an output file to an input file.

Functions and Equivalent Circuit of the Engine Code.
The engine code for the circuit simulation program was used to simulate the effects of an instrument loop degradation on a measurement condition.This engine code is useful for obtaining a correlation between the simulation circuit conditions and the actual circuit conditions.Although the goal was to develop a real model, the intricacies of developing the appropriate waveforms and incorporating the effects of the induced noise into the circuit turned a degraded circuit loop into lines and lines of netlists.A 4∼20 mA loop was used to measure the reactor coolant pressure the same as for most instruments.The overall circuit loop consists of a section of a control room, a section of a twisted pair shield cable, and a section of a transmitter.A transmitter circuit was treated as a loop for an important circuit for the passive elements R-L-C which could be affected by the severe accident conditions.In this case, the imported passive elements of the transmitter circuit were R 2 , R 3 , C 1 , and C 2 , as distributed parameters  of the circuit.Figure 6 represents an overall circuit model modified to run in the MultiSIM code.

Procedure for a One Body Order Simulation and Digital
Signal Processing.For the first step, during the operating procedure, the establishment of a simulation circuit modeling by the PSpice/MultiSIM code is attempted, where the PSpice/MultiSIM simulation output file is used for all the files; * .cir,* .sim,and * .netfiles.The next step needs the information on the circuit elements from the * .netfile, and then the * .netfile can be controlled by an element in the circuit, for using the Labview * .netfile condition, and then we take the csd file out of the * .netfile in the Labview and execute the system exec.vifile of the PSpice/MultiSIM condition.As the last step, we extract the output results from the csd file, and then the file can be displayed as some of the output abnormal signals.Finally, these signals can be processed by an analyzing tool.Thus, this system will be used for circumstance diagnostic for the abnormal signals, and this code will also allow us to use the enhanced  of the PSpice/MultiSIM engine code and the LabVIEW tool code, which can work as a one body order system for this simulation.The overall simulation result data can be used by one operation system as a one order command.So that the programs could be operated quickly and conveniently as an onsite instrument, in this simulation, we changed the element values in the circuit which is possible by directly changing the R-C element value according to the temperature condition in the initial screen menu of the tool code, so we do not need to use the functions of the engine code.
As a result, it is also possible to change the output signal patterns according to the changing element values in the tool code menu functions under SACs.It is also easy and convenient for a signal analysis for some types of noise patterns.The output results can be extended and analyzed for their response characteristics by the tool code, because the LabVIEW code has powerful analyzing functions and good interface while being onsite.Figure 7 shows the block diagram for a one body order simulation system.The loop circuit was composed of three sections which include the pressure transmitter section, the twisted shielded pair connection cable section, and the receiver section.The pressure transmitter has been simplified to obtain an equivalent representation of the resisters, capacitors, and inductors as distributive parameters of the circuit.It also includes voltage or current sources that can be varied through the circuit simulation code to obtain realistic circuit response characteristics from environmental changes in the containment.A specific parameter includes: an input leakage (shunt) resistance that is affected by damage to an interconnecting cable and; input series resistance to represent the spliced connections that might be affected by corrosion and whiskers phenomenon.Input capacitance is increased easily by a water intrusion into the cable or the transmitter housing.Capacitance for the ground is easily increased by the presence of water.A diode and load resistors represent the power supply current.A current generator represents the actual transmitter output signal.A typical cable is represented to describe the wire inductance and resistance, and the dielectric capacitance.The shield is represented as a third conductor.It is simplified to eliminate all the resistances, and to only model the inductance and capacitance.The voltage source drives a current to a voltage sensing resistor in series with the transmitter.Negative sides of the voltage source and the shields are attempted for the ground source in this analysis.Figure 8 shows a detail block diagram for the signal processing and analysis from the ASSA module system.Figure 9 is a screen of the ASSA for one body order simulation system.

Simulation and Response Characteristic Analysis from the Obtained Data
There are two types regarding information of pulse parameters which means changing a value for the output pulse signal with certain information to diagnosis a circuit condition such as a voltage level, a time constant including information of a rising time, falling time, and bandwidth compare to information for pulse parameters of the reference input signal.As a result, when we can analyze the cause of an equivalent circuit, which consists of an R-L-C distribution parameter, for a 4∼20 mA instrument of a severe accident condition from response characteristics, it is also the same meaning that is also possible to detect and diagnose the abnormal signal of a malfunctioning instrument from a simulation and response characteristic analysis.Figure 10 is patterns of output pulse signal of ASSA for a one body order simulation system.Method for the pulse parameter estimation: for pulse parameter, estimation for ASSA module has designed as the basic concept for circuit analysis procedure.In addition to this method, there is pulse detection logic, an abnormal pulse filter, and a pulse parameter estimation for pulse parameter estimation was already said and it is redundant.The pulse detection logic is operable to detect pulses and to estimate a pulse parameter of each detected pulse [10].The abnormal pulse filter is operable to analyze a pulse parameter of the detected pulses and to determine whether there are any abnormal pulses based on the analysis.The pulse parameter estimation is operable to estimate a pulse parameter for the signal based on an analysis of the detected pulses and information pertaining to there being any abnormal pulses in the signal.

Defining the Reference Signal from the Equivalent Circuit.
Figure 11 is reference input current (A) and output voltage (V) data for the pulse parameter as an initial value.This reference pulse signal was defined in the engine code as a linear approximation, with a 0-16 mA pulse with a 0 ms delay, 30 ms rising time, 30 ms falling time, and a 1 PPS repetition rate.The reference pulse also describes an output current pulse as measured across the current to voltage converter at the receiver.This is a reference input current I (A) and output voltage V (V) data of the pulse parameter as an initial value.In this figure, it appears as about 0.1 to 1.0 (V) pulse at a temperature of 27 ( • C).  5 (Ω)∼900 (Ω) and 1 (kΩ)∼10 (kΩ), respectively.Figure 14 shows the variations of output voltage according to the R 2 .

Results with
In the case of the R 2 simulation using the tool code, it can be seen that the resistor value changes for the high-voltage level data which means a good linearity (sensitivity) response characteristic over 1 (kΩ) to 10 (kΩ).In this resistance R 2 , the impulse response of an equivalent circuit for the 4-20 mA is its output when presented with an input standard signal called an impulse.The impulse response refers to the reaction of the circuit system which was consisted with that resistance R 2 to be obtain the response characteristic for the maximum resistance changing value facing each external temperature as the changing parameter, the response characteristic has an increasing characteristic in the circuit according to basic circuit theory, because the resistor (R 2 ) was consisted of the parallel component.Figures 15 and 16 show the output voltage level of the pulse parameter according to a change of the resistance R 3 which has a range of 5 (Ω)∼4 (kΩ) and 5 (kΩ)∼ 50 (kΩ), respectively.In the case of the R 3 simulation using the tool code, it can be seen that the resistor value changes for the high voltage level data which means a good linearity (sensitivity) response characteristic over 5 (Ω) to 50 (kΩ), in the case of an overall response characteristic with a decreasing response characteristic for an increasing resistance because the resistor (R 3 ) was consisted of the serial component in Figure 17 according to the basic circuit theory.Figure 17 shows a response characteristic of the output voltage for the pulse parameter for a change of the resistance of R 3 (5 (Ω)∼50 (kΩ)).

Results with Change of Capacitances.
Figure 18 shows the output voltage level data and the time constant for the pulse parameter according to a change of the capacitance of C 1 which has a range of 33 (nF)∼33 (mF).Figure 19 shows the response characteristics of the output voltage curve for the pulse parameter according to a change of the capacitance of C 1 .Figure 20 shows the response characteristics of the rising time (Tr) for the pulse parameter according to a change of the capacitance of C 1 .In the case of the C 1 simulation by using the tool code, it can be seen that the capacitance   value changes for the high voltage level data which means a good linearity (sensitivity) response characteristic over the same range, in the case of an overall response characteristic with a decreasing response characteristic for an increasing capacitance.The response characteristic of the rising time (Tr) for the pulse parameter which also means a good linearity (sensitivity) for the response characteristics over the same range, in the case of an overall response characteristic, reveals an increasing response characteristic for an increasing capacitance of C 1 .In this capacitance C 1 , the impulse response of an equivalent circuit for the 4-20 mA is its output when presented with an input standard signal called an impulse.The impulse response refers to the reaction of the circuit system which was consisted with that capacitance C 1 to obtain the response characteristic for the maximum capacitance changing value facing each external temperature as the changing parameter, the response characteristic of the output voltage curve is a decreasing characteristic, the response characteristic of the rising time (Tr) is an increasing characteristic in the circuit according to basic circuit theory, because the capacitance (C 1 ) was consisted of the parallel component.
Figure 21 shows the output voltage level data for the pulse parameter according to a change of the capacitance of C 2 which has a range of 10 (nF)∼10 (mF).Figure 22 shows the response characteristics of the output voltage level data for the pulse parameter according to a change of the capacitance of C 2 .Figure 23 shows the response characteristics of the rising time (Tr) for the pulse parameter according to a change  of the capacitance of C 2 .In the case of the C 2 simulation using the tool code, it can be seen that the capacitance value changes for the high voltage level data which means a good linearity (sensitivity) response characteristic over the same range, in the case of an overall response characteristic with a decreasing response characteristic for an increasing capacitance at 1.31 (mF), until about 1.31 (mF) it was a steady state.But the response characteristic of the rising time (Tr) for the pulse parameter which means a good linearity (sensitivity) for the response characteristics over the same range, in the case of an overall response characteristic reveals an increasing response characteristic for an increasing capacitance of C 2 in the circuit according to basic circuit theory.

Correlative Equation between R 2 and Temperature.
There is a correlative equation between the temperature and the amplitude of an output voltage in the equivalent circuit, which means that R 2 has been changing in the resistance range as a distributed parameter when the T temperature has been changing in the temperature range due to the accident conditions.Equation ( 1) is a correlative equation between the temperature and the amplitude of an output voltage.Figure 24 shows a relation curve between the T temperature and the amplitude of the V out according to (1).From this equation, for an inverse operation, it could also be said that T temperature has been changing in the temperature range when R 2 has been changing in the resistance as a distributed parameter in the equivalent circuit due to the accident conditions.In this case of the changing R 2 resistance value according to the changing temperature, experimental data is obtained by using in an electric oven as a real test.Equation ( 2) is also a correlative equation between the temperature and the amplitude of an output voltage, in the case of the changing R 2 resistance according to the T temperature.
Figure 25 shows the relation curve between the R 2 and the T temperature according to (2).Table 1 shows the changing R 2 resistance range according to the temperature changing range in a real test using in an electric oven: where T: changing temperature, T 0 : normal/initial temperature, R 0 : normal/initial resistance, R 2 (T): changing resistance at T temperature, c 1 : 1st-order temperature coefficient, i m : maximum input pulse current, v 0 : the normal/initial voltage, and R 3 , R 7 , R 8 : the circuit elements.

Relation between the C
1 and the Time Constant.The changing range for the time constant of the output voltage in the equivalent circuit when the C 1 has a maximum changing range from 33 nF to 33 mF, to be placed in extreme simulation conditions should be investigated.From the simulation result, it is also possible to explain an approximate equation for the relation between the C 1 and the time constant.In the case of above 33 uF, the time constant of the output voltage has a response characteristic as a steep  slope but it has response characteristic of almost no change when the C 1 has a changing range from 33 nF to 33 uF.Therefore, the time constant of the output voltage has a linear response characteristic when the C 1 has a changing range from 33 uF to 33 mF. Figure 26 shows a linear response characteristic between the rising time parameter and when C 1 has a changing range from 33 uF to 33 mF.Equation (3) defines a relation between the C 1 and the time constant as a rising time parameter from an output voltage shape.When the t is a time constant as a rising time parameter, c is the 1st-order temperature coefficient for the C 1 changing range.To obtain a rising time for R-C elements, Figure 27 shows changing output data for the pulse parameter for a changing range of the temperature of C 2 (27 ( • C) ∼ 87 ( • C)) according to correlation equation temperature coefficient to the time constant as a rising time parameter for R-C elements.Table 2 is output parameters of the pulse parameter according to a change of the capacitance of C 1 : where t: time constant as a rising time parameter, c: 1st-order temperature coefficient.

Data Analysis.
We could obtain frequency response characteristics of the multidimensions from the pulse parameter for a changing range of the capacitance C 1 .From this Fourier fast transform (FFT) analysis of the 3 dimensions plot, it will find the activity of distorting shape and it become loosing for a direct current (DC) property for a changing range.Therefore, it also will be able to provide correct  of the pulse parameter for a change of the capacitance C 1 , respectively.

Conclusions
In this paper, we provided a new simulator (called ASSA) developed through a circuit simulation with instrument diagnosis methods.Its underlying design requirements are fundamentally based on both EPRI TR-103412 and EPRI TR-102371 covering a variety of information for the performance assessment of plant instrumentation under severe accident conditions.As an alternative device to real systems, the simulator makes it possible for operators to take a correct action through a reasonable assessment of the performance of the plant instrumentation, even for severe accidents.The ASSA module system has been designed through a realization of a one body order system of one program linked to signal analysis tool programs and to an abnormal signal simulation engine as a one order command.As a result, the program will be able to analyze data quickly and conveniently, onsite.
Using the ASSA module system, we could obtain the response characteristics from the output voltage levels of a pulse parameter according to a change of the resistance for each of R 2 , and R 3 , and the response characteristic from the output voltage level data and time constant of the pulse parameter according to a change of the capacitance for each of C 1 , and C 2 .For the case of the pulse parameter following a change of the resistance, we obtained good response characteristics for the output voltage level data.For the case of the pulse parameter following a change of the capacitance, we could obtain much better response characteristics for the output voltage level data including the time constant patterns.With the ASSA module system, we additionally obtained a type of the varying signal patterns by a simple change of another element's R-C values.The forgoing function of the ASSA module makes it possible to change the abnormal signal patterns by changing the element R-C value and by using specific signal processing.In this paper, we can also obtain a correlative equation between the temperature and the amplitude of an output voltage in the equivalent circuit, which means that an R 2 has been changing in the resistance range as a distributed parameter when the T temperature has been changing in the temperature range by accident conditions.While the system module makes a signal analysis of noise patterns easy and convenient, it has not been built as a full simulation system as of yet.To make the ASSA module fully available for harsh environments, further research must be focused on developing an R-L-C (including active elements: a diode, N-P-N, P-N-P transistors) correlation equation for reflecting the effects of the circumstance parameters such as a high temperature, high pressure, high humidity, and high radiation.Since a variety of simulations to check on malfunctions or abnormal indications from instruments are possible through the ASSA module, it is possible to develop an alternative method to check on a degradation of the instruments with diagnosis logics for severe conditions and a database from the ASSA module analysis results.

Figure 1 :
Figure 1: Flowchart for the circuit diagnosis methods among the five types of aids.

Figure 2 :
Figure 2: Initial screen consisting of functions for ASSA module.

Figure 3 :
Figure 3: Flowchart of the overall design steps for the simulator as ASSA module.

Figure 4 :
Figure 4: Three EQ conditions and five phases of a scenario taken by MAAP4 code.

Figure 5 :
Figure 5: Screen of comparative assessment procedure as the first step for the ASSA module.

Figure 6 :
Figure 6: Overall circuit model modified to run in the MultiSIM code.

Figure 7 :
Figure 7: Block diagram for one body order simulation system.

Figure 8 :
Figure 8: Block diagram for the signal processing and analysis from ASSA module system.

Figure 9 :
Figure 9: A screen of the ASSA for a one body order simulation system.

Figure 10 :Figure 11 :
Figure 10: Patterns of output pulse signal of ASSA for a one body order simulation system.

Figure 17 :
Figure 17: Response characteristic of the output voltage for the pulse parameter for a change of the resistance of R 3 (5 (Ω)∼50 (kΩ)).

Figure 18 :
Figure 18: Output data for the pulse parameter for a change of the capacitance of C 1 (33 (nF)∼33 (mF)).

Figure 19 : 1 Figure 20 :
Figure 19: Response characteristic of the output voltage curve for the pulse parameter for a change of the capacitance C 1 (33 (nF)∼33 (mF)).

Figure 21 :Figure 22 :
Figure 21: Output data for the pulse parameter for a change of the capacitance of C 2 (10 (nF)∼10 (mF)).

2 Figure 23 :Figure 24 :
Figure 23: Response characteristic of the rising time for the pulse parameter for a change of the capacitance of C 2 (10 (nF)∼10 (mF)).

Figure 25 :
Figure 25: Relation curve between R 2 and the T temperature according to (2).

Figure 26 :
Figure 26: A linear response characteristic between the rising time and C 1 .

Figure 27 :
Figure 27: Output data for the pulse parameter for a change of the temperature of C 2 (27 ( • C) ∼ 87 ( • C)).

Figure 28 :Figure 29 :
Figure 28: Frequency response characteristic of two dimensions of the pulse parameter for a change of the capacitance C 1 .
(b) shows five phases of

Table 1 :
R 2 Resistance changing range according to the temperature changing range in a real test.

Table 2 :
Output data of the pulse parameter according to a change of the capacitance of C 1 .