Method suitable for starting MAAP-CANDU software under any accident condition

Abstract

The invention belongs to the technical field of nuclear power plant simulation, and particularly relates to a method suitable for starting MAAP-CANDU software under any accident condition. Comprising the following steps: step 1, establishing an MAAP model of a CANDU6 unit, and performing steady-state and typical serious accident working condition tests; 2, loading a model parameter file of the CANDU6 unit, and entering a dormant state after initialization; 3, parameter variables, generated by the full-range simulator, of system equipment under the accident condition and needed by initialization of MAAP-CANDU software are transmitted to the SIM program; 4, calculating initial values of state quantities of different working media; 5, if the software is initialized to crash, checking the error reporting module, checking the initial input parameter value of the module, and repeating the step 3; and 6, after the software is successfully initialized, smooth replacement calculation is carried out. The method has the advantages that the full-range simulator can schedule MAAP-CANDU software in real time, and the MAAP software and the CANDU software can calculate synchronously; the MAAP-CANDU software can be started from any accident working condition, and smooth replacement calculation can be realized.

Description

Technical Field

[0001] The invention belongs to the technical field of nuclear power plant simulation, and in particular relates to a method suitable for starting any accident condition of MAAP-CANDU software. Background Art

[0002] MAAP (The Modular Accident Analysis Program) is an integrated, modular program developed by the Electric Power Research Institute (EPRI) and FAI (Fauske & Associates) for severe accident analysis and calculation in nuclear power plants.

[0003] Based on the MAAP program, FAI and AECL have jointly developed a severe accident analysis software MAAP-CANDU suitable for the CANDU6 reactor type by integrating the core heating module program of OPG (Ontario Power Generation). The MAAP-CANDU software package contains a horizontally placed fuel channel model and CANDU reactor-specific equipment and systems, such as calandria, emergency core cooling system, moderator cooling system and end shield cooling system. MAAP-CANDU software can analyze and calculate severe accident phenomena such as evaporation of the main heat transfer system, core heating and melting, fuel channel failure, core disintegration, calandria failure, end shield failure and containment failure.

[0004] At present, the full-range simulator developed based on the real CANDU6 unit can simulate the normal operating conditions and design benchmark conditions of the CANDU6 unit, and plays a pivotal role in training operators of the CANDU6 unit. The full-range simulator of the CANDU6 unit cannot simulate severe accidents, so it is necessary to integrate MAAP-CANDU software to cope with the current dilemma.

[0005] Therefore, it is necessary to invent a method suitable for starting any accident condition of MAAP-CANDU software to perform follow-up calculations, expand the simulation range of the CANDU6 unit full-range simulator, and provide a platform for the training and drills of SAMG of CANDU6 units. Summary of the invention

[0006] The purpose of the present invention is to provide a method suitable for starting any accident condition of MAAP-CANDU software, by transferring important parameters of relevant modules such as the core, main heat transfer system, stabilizer, steam generator, degassing condenser, containment and calandria of a full-range simulator of a CANDU6 unit to the MAAP-CANDU software, and initializing the MAAP-CANDU software, so as to achieve the purpose of starting any accident condition of the MAAP-CANDU software, solve the problem that the simulation range of the full-range simulator of the CANDU6 unit is limited, expand the simulation model of the CANDU6 unit, and enable the MAAP-CANDU software to smoothly take over the calculation of serious accident phenomena of the CANDU6 unit.

[0007] The technical solution of the present invention is as follows: A method applicable to starting any accident condition of MAAP-CANDU software, comprising the following steps:

[0008] Step 1: Use MAAP-CANDU software to establish the MAAP model of the CANDU6 unit, and conduct steady-state and typical severe accident condition tests to verify the reliability of the simulation model;

[0009] Step 2: MAAP-CANDU software loads the model parameter file of the CANDU6 unit under the scheduling of the MAAPRinMem program and enters the dormant state after initialization;

[0010] Step 3: After the set time, the parameter variables of the system equipment under the accident condition required for the initialization of the MAAP-CANDU software generated by the full-range simulator are passed to the SIM program;

[0011] Step 4: After receiving the important parameters of the limited system equipment, the SIM program completes the calculation of the initial values ​​of the state quantities of different working fluids in the submodule;

[0012] Step 5: If the MAAP-CANDU software crashes during initialization, check the error module and verify the initial input parameter values ​​of the module, that is, the rationality of the parameter variable ranges of each module listed in step 3, and repeat step 3;

[0013] Step 6: After the MAAP-CANDU software is successfully initialized, a smooth succession calculation is performed. The process system model of the original full-range simulator provides the calculation boundary for the MAAP-CANDU software, and the MAAP-CANDU software responds to the original fault information.

[0014] The parameter variables in step 3 include: important time and control signals, core parameters, main heat transfer system parameters, pressurizer parameters, steam generator parameters, degassing condenser parameters, containment parameters, and calandria parameters.

[0015] The important time and control signals include the stack shutdown time and the pump shutdown time.

[0016] The core parameters include decay heat and core node temperature.

[0017] The main heat transfer system parameters include total pressure, steam partial pressure, hydrogen partial pressure, nitrogen partial pressure, oxygen partial pressure; air temperature, water temperature; cavitation fraction, node natural circulation state; node water quality.

[0018] The pressurizer parameters include total pressure, steam partial pressure, hydrogen partial pressure, nitrogen partial pressure, oxygen partial pressure, air temperature, water temperature, cylinder wall temperature, fluctuation pipe wall temperature, and pressurizer water quality.

[0019] The steam generator parameters include the pressure, steam partial pressure, hydrogen partial pressure, nitrogen partial pressure, oxygen partial pressure of the damaged loop steam generator, the air temperature and water temperature of the damaged loop steam generator, the water quality of the damaged loop steam generator, the cylinder wall temperature of the damaged loop steam generator, the pressure, steam partial pressure, hydrogen partial pressure, nitrogen partial pressure, oxygen partial pressure of the undamaged loop steam generator, the air temperature and water temperature of the undamaged loop steam generator, the water quality of the undamaged loop steam generator, and the cylinder wall temperature of the undamaged loop steam generator.

[0020] The degassing condenser parameters include total pressure, steam partial pressure, hydrogen partial pressure, nitrogen partial pressure, oxygen partial pressure, air temperature, water temperature, and water quality; the containment parameters include node total pressure, steam partial pressure, hydrogen partial pressure, nitrogen partial pressure, oxygen partial pressure, node air temperature, and node water quality; the calandria container includes total pressure, steam partial pressure, hydrogen partial pressure, nitrogen partial pressure, oxygen partial pressure, air temperature, water temperature, water quality, and calandria container wall temperature.

[0021] In step 4, the SIM program uses the line interpolation method to obtain the thermophysical property functions of light water and light water vapor respectively according to the thermophysical property parameter table, and then calculates the specific volume, specific enthalpy and water level of water at the startup time through the known temperature, pressure and water mass of the system equipment. The gas is calculated by the ideal gas equation to obtain the corresponding gas mass at the startup time. The heat sink temperature is calculated from the inner wall temperature through the uniform distribution function and the empirical coefficient. The SIM program completes the assignment of initial values ​​of the state quantities of different working fluids in each module, and then calls the original mass conservation ordinary differential equation and energy conservation ordinary differential equation of the MAAP-CANDU software to complete the solution of the state quantity at the next moment, that is, the initialization of the MAAP-CANDU software at the time of the accident condition is completed.

[0022] The step 4 specifically includes the following:

[0023] The line interpolation function is as follows:

[0024]

[0025] Among them, r is the enthalpy value; T is the system equipment parameter temperature; T1 and T2 are the corresponding pressures of the temperature in the thermophysical parameter table;

[0026] The water level calculation method except for the main heat transfer system is as follows:

[0027] L=Mv / (S)

[0028] Wherein, L is the water level; M is the water mass; v is the specific volume; S is the bottom area of ​​the equipment model;

[0029] The ideal gas equation is shown below:

[0030]

[0031] Where P is the partial pressure of the gas; V is the volume of the gas; T is the temperature; n is the molar amount of the gas substance; R is the molar gas constant; M is the mass of the gas; μ is the molar mass of the gas;

[0032] The solution formula for the heat sink temperature distribution obtained from the wall temperature is as follows:

[0033] T THS (i) = T THS -δ(T THS -T CONT )*i / X node

[0034] TTHS(i) is the temperature of the heat sink node; T THS is the known equipment wall temperature; Tcont is the containment node temperature of the equipment compartment; δ is the empirical coefficient factor; i is the node number of the radial division of the corresponding equipment wall; X node The number of nodes corresponding to the radial division of the equipment wall.

[0035] The beneficial effects of the present invention are: 1) the full-range simulator can realize real-time scheduling of the MAAP-CANDU software, and the two can be calculated synchronously; 2) the MAAP-CANDU software can be started from any accident condition, and smooth succession calculation can be realized; 3) after the MAAP-CANDU software is started from any accident condition, the subsequent accident process can be correctly simulated; 4) the calculation boundary of the MAAP-CANDU software can be coupled with the process system model of the simulator; 5) the MAAP-CANDU software can correctly respond to the fault information of the full-range simulator. BRIEF DESCRIPTION OF THE DRAWINGS

[0036] Figure 1 A flow chart of a method for starting any accident condition of MAAP-CANDU software provided by the present invention;

[0037] Figure 2Start the business process for any accident condition of MAAP-CANDU software;

[0038] Figure 3 This is the data processing flow chart of the SIM program. DETAILED DESCRIPTION

[0039] The present invention is further described in detail below with reference to the accompanying drawings and specific embodiments.

[0040] The present invention is based on CANDU6 unit data, uses MAAP-CANDU software to establish a MAAP model of the CANDU unit, and then realizes data communication between the MAAP-CANDU program and the simulator through the developed MAAPRinMem program. The MAAP-CANDU software accepts the scheduling of the simulator as a subroutine. At a certain artificially specified time, the simulator generates a set of key parameters of the system equipment under the accident condition, passes the parameters to the MAAP-CANDU software, calls the initialization program SIM developed specifically for the MAAP-CANDU software, initializes the MAAP-CANDU software, and enables the MAAP-CANDU software to realize the startup of the accident condition and to continue the calculation. After the MAAP-CANDU software is initialized, the model of the original full-range simulator replaced by the MAAP-CANDU software can be calculated for a short period of time and then officially exited, and other process system models that have not been replaced can provide calculation boundaries for the MAAP-CANDU software.

[0041] A method applicable to starting any accident condition of MAAP-CANDU software, comprising the following steps:

[0042] Step 1: Use MAAP-CANDU software to establish the MAAP model of the CANDU6 unit, and conduct steady-state and typical severe accident condition tests to verify the reliability of the simulation model;

[0043] Step 2: MAAP-CANDU software, under the scheduling of MAAPRinMem program (a data communication program specially developed for MAAP-CANDU software and simulator), loads the model parameter file of CANDU6 unit and enters the dormant state after initialization;

[0044] Step 3: After reaching the manually set time, the parameter variables of the system equipment under the accident condition required for the initialization of the MAAP-CANDU software generated by the full-range simulator are passed to the SIM program. The parameter variable range includes:

[0045] Important time and control signals: reactor shutdown time, pump shutdown time;

[0046] Core parameters: decay heat, core node temperature;

[0047] Main heat transfer system parameters: total pressure, steam partial pressure, hydrogen partial pressure, nitrogen partial pressure, oxygen partial pressure; air temperature, water temperature; cavitation fraction, node natural circulation state; node water quality;

[0048] Stabilizer parameters: total pressure, steam partial pressure, hydrogen partial pressure, nitrogen partial pressure, oxygen partial pressure, air temperature, water temperature, cylinder wall temperature, fluctuation pipe wall temperature, stabilizer water quality;

[0049] Steam generator parameters: pressure, steam partial pressure, hydrogen partial pressure, nitrogen partial pressure, oxygen partial pressure of the damaged loop steam generator, air temperature and water temperature of the damaged loop steam generator, water quality of the damaged loop steam generator, cylinder wall temperature of the damaged loop steam generator, pressure, steam partial pressure, hydrogen partial pressure, nitrogen partial pressure, oxygen partial pressure of the undamaged loop steam generator, air temperature and water temperature of the undamaged loop steam generator, water quality of the undamaged loop steam generator, cylinder wall temperature of the undamaged loop steam generator;

[0050] Degassing condenser parameters: total pressure, steam partial pressure, hydrogen partial pressure, nitrogen partial pressure, oxygen partial pressure, air temperature, water temperature, water quality;

[0051] Containment parameters: node total pressure, steam partial pressure, hydrogen partial pressure, nitrogen partial pressure, oxygen partial pressure, node air temperature, node water quality;

[0052] Calandria vessel parameters: total pressure, steam partial pressure, hydrogen partial pressure, nitrogen partial pressure, oxygen partial pressure, air temperature, water temperature, water quality, calandria vessel wall temperature;

[0053] Step 4: After receiving the important parameters of the limited system equipment, the SIM program will complete the calculation of the initial values ​​of the state quantities of different working fluids in the submodule.

[0054] The SIM program receives the important parameters of the limited system equipment and processes them as follows: Figure 3 As shown. The SIM program uses the line interpolation method to obtain the thermophysical property functions of light water and light water steam respectively according to the thermophysical property parameter table. Then the specific volume, specific enthalpy and water level of water at the startup moment are calculated by the known temperature, pressure and water mass of the system equipment. The gas is calculated by the ideal gas equation to obtain the corresponding gas mass at the startup moment. The heat sink temperature is calculated from the inner wall temperature through the uniform distribution function and the empirical coefficient. At this point, the SIM program completes the assignment of initial values ​​of the state quantities of different working fluids in each module. Subsequently, the original mass conservation ordinary differential equation and energy conservation ordinary differential equation of the MAAP-CANDU software are called to complete the solution of the state quantity at the next moment, that is, to complete the initialization of the MAAP-CANDU software at the moment of the accident condition.

[0055] The line interpolation function is as follows:

[0056]

[0057] Where r is the enthalpy (J / kg), specific volume (m 3 / kg) and other thermal properties; T is the system equipment parameter temperature, K; T1 and T2 are the temperatures inside and outside the thermal property parameter table (under corresponding pressure), K.

[0058] Schematic diagram of water level calculation method (except main heat transfer system):

[0059] L=Mv / (S)

[0060] L is the water level, m; M is the water mass, kg; v is the specific volume, m 3 / kg; S is the bottom area of ​​the equipment model (which changes with height).

[0061] The ideal gas equation is shown below:

[0062]

[0063] Where P is the gas partial pressure, Pa; V is the gas volume, m 3 ; T is temperature, K; n is the molar amount of gas substance, mol; R is the molar gas constant, J / (mol*K); M is the mass of gas, kg; μ is the molar mass of gas, kg / mol.

[0064] The solution formula for the heat sink temperature distribution obtained from the wall temperature is as follows:

[0065] T THS (i) = T THS -δ(T THS -T CONT )*i / X node

[0066] TTHS(i) is the temperature of the heat sink node, K; T THS is the known equipment wall temperature, K; Tcont is the containment node temperature of the equipment compartment, K; δ is the empirical coefficient factor; i is the node number of the radial division of the corresponding equipment wall; X node The number of nodes corresponding to the radial division of the equipment wall.

[0067] Step 5: If the MAAP-CANDU software crashes during initialization, check the error module and verify the initial input parameter values ​​of the module, that is, the rationality of the parameter variable ranges of each module listed in step 3, and repeat step 3;

[0068] Step 6: After the MAAP-CANDU software is successfully initialized, a smooth succession calculation is performed. The process system model of the original full-range simulator can provide a calculation boundary for the MAAP-CANDU software, and the MAAP-CANDU software responds to the original fault information.

[0069] A method applicable to starting up any accident condition of MAAP-CANDU software comprises the following steps: using MAAP-CANDU software to establish a MAAP model of a CANDU6 unit; under the scheduling of a MAAPRinMem program, the MAAP-CANDU software loads a model parameter file of the CANDU unit and enters a dormant state after the first step of initialization; after reaching a manually set time, important parameters of system equipment under accident conditions required for initialization of the MAAP-CANDU software generated by a full-range simulator are passed to a SIM program, and the initialization of the accident condition time of the MAAP-CANDU software is completed through the program, and the parameter variable range includes important time and control signals, important parameters of relevant modules such as a core, a main heat transfer system, a regulator, a steam generator, a degassing condenser, a containment, and a calandria; after the MAAP-CANDU software is successfully initialized, the calculation can be smoothly taken over, the process system model of the original full-range simulator can provide a calculation boundary for the MAAP-CANDU software, and the MAAP-CANDU software can respond to the original fault information.

[0070] The method of starting the MAAP-CANDU software under any accident condition based on this method is similar to the method of starting the MAAP-CANDU software under any accident condition in other software. If other people's products only contain some technical features, they are not infringing. Only when they cover all the above technical features will they be infringing.

Claims

1. A method applicable to the startup of any accident condition of MAAP-CANDU software, characterized in that: The steps include: Step 1: Use MAAP-CANDU software to establish the MAAP model of the CANDU6 unit, and conduct steady-state and typical severe accident condition tests to verify the reliability of the simulation model; Step 2: MAAP-CANDU software loads the model parameter file of the CANDU6 unit under the scheduling of the MAAPRinMem program and enters the dormant state after initialization; Step 3: After the set time, the parameter variables of the system equipment under the accident condition required for the initialization of the MAAP-CANDU software generated by the full-range simulator are passed to the SIM program; Step 4: After receiving the important parameters of the limited system equipment, the SIM program completes the calculation of the initial values ​​of the state quantities of different working fluids in the submodule; Step 5: If the MAAP-CANDU software crashes during initialization, check the error module and verify the initial input parameter values ​​of the module, that is, the rationality of the parameter variable ranges of each module listed in step 3, and repeat step 3; Step 6: After the MAAP-CANDU software is successfully initialized, a smooth succession calculation is performed. The process system model of the original full-range simulator provides the calculation boundary for the MAAP-CANDU software, and the MAAP-CANDU software responds to the original fault information.

2. A method for starting the MAAP-CANDU software in any accident condition as claimed in claim 1, characterized in that: The parameter variables in step 3 include: important time and control signals, core parameters, main heat transfer system parameters, pressurizer parameters, steam generator parameters, degassing condenser parameters, containment parameters, and calandria parameters.

3. A method for starting the MAAP-CANDU software in any accident condition as claimed in claim 2, characterized in that: The important time and control signals include the stack shutdown time and the pump shutdown time.

4. A method for starting the MAAP-CANDU software in any accident condition as claimed in claim 2, characterized in that: The core parameters include decay heat and core node temperature.

5. A method for starting the MAAP-CANDU software in any accident condition as claimed in claim 2, characterized in that: The main heat transfer system parameters include total pressure, steam partial pressure, hydrogen partial pressure, nitrogen partial pressure, oxygen partial pressure; air temperature, water temperature; cavitation fraction, node natural circulation state; node water quality.

6. A method for starting the MAAP-CANDU software in any accident condition as claimed in claim 2, characterized in that: The pressurizer parameters include total pressure, steam partial pressure, hydrogen partial pressure, nitrogen partial pressure, oxygen partial pressure, air temperature, water temperature, cylinder wall temperature, fluctuation pipe wall temperature, and pressurizer water quality.

7. A method for starting up any accident condition of MAAP-CANDU software as claimed in claim 2, characterized in that: The steam generator parameters include the pressure, steam partial pressure, hydrogen partial pressure, nitrogen partial pressure, oxygen partial pressure of the damaged loop steam generator, the air temperature and water temperature of the damaged loop steam generator, the water quality of the damaged loop steam generator, the cylinder wall temperature of the damaged loop steam generator, the pressure, steam partial pressure, hydrogen partial pressure, nitrogen partial pressure, oxygen partial pressure of the undamaged loop steam generator, the air temperature and water temperature of the undamaged loop steam generator, the water quality of the undamaged loop steam generator, and the cylinder wall temperature of the undamaged loop steam generator.

8. A method for starting the MAAP-CANDU software in any accident condition as claimed in claim 2, characterized in that: The degassing condenser parameters include total pressure, steam partial pressure, hydrogen partial pressure, nitrogen partial pressure, oxygen partial pressure, air temperature, water temperature, and water quality; the containment parameters include node total pressure, steam partial pressure, hydrogen partial pressure, nitrogen partial pressure, oxygen partial pressure, node air temperature, and node water quality; the calandria container includes total pressure, steam partial pressure, hydrogen partial pressure, nitrogen partial pressure, oxygen partial pressure, air temperature, water temperature, water quality, and calandria container wall temperature.

9. A method for starting any accident condition of MAAP-CANDU software as claimed in claim 1, characterized in that: In step 4, the SIM program uses the line interpolation method to obtain the thermophysical property functions of light water and light water vapor respectively according to the thermophysical property parameter table, and then calculates the specific volume, specific enthalpy and water level of water at the startup time through the known temperature, pressure and water mass of the system equipment. The gas is calculated by the ideal gas equation to obtain the corresponding gas mass at the startup time. The heat sink temperature is calculated from the inner wall temperature through the uniform distribution function and the empirical coefficient. The SIM program completes the assignment of initial values ​​of the state quantities of different working fluids in each module, and then calls the original mass conservation ordinary differential equation and energy conservation ordinary differential equation of MAAP-CANDU to complete the solution of the state quantity at the next moment, that is, completes the initialization of the MAAP-CANDU software at the time of the accident condition.

10. A method for starting up any accident condition of MAAP-CANDU software as claimed in claim 9, characterized in that: The step 4 specifically includes the following: The line interpolation function is as follows: Among them, r is the enthalpy value; T is the system equipment parameter temperature; T1 and T2 are the corresponding pressures of the temperature in the thermophysical parameter table; The water level calculation method except for the main heat transfer system is as follows: L=Mv / (S) Wherein, L is the water level; M is the water mass; v is the specific volume; S is the bottom area of ​​the equipment model; The ideal gas equation is shown below: Where P is the partial pressure of the gas; V is the volume of the gas; T is the temperature; n is the molar amount of the gas substance; R is the molar gas constant; M is the mass of the gas; μ is the molar mass of the gas; The solution formula for the heat sink temperature distribution obtained from the wall temperature is as follows: T THS (i) = T THS -δ(T THS -T CONT )*i / X node TTHS(i) is the temperature of the heat sink node; T THS is the known equipment wall temperature; Tcont is the containment node temperature of the equipment compartment; δ is the empirical coefficient factor; i is the node number of the radial division of the corresponding equipment wall; X node The number of nodes corresponding to the radial division of the equipment wall.

Images