CN114065493A - Method and system for optimizing stability of analysis program of reactor system - Google Patents
Method and system for optimizing stability of analysis program of reactor system Download PDFInfo
- Publication number
- CN114065493A CN114065493A CN202111310736.5A CN202111310736A CN114065493A CN 114065493 A CN114065493 A CN 114065493A CN 202111310736 A CN202111310736 A CN 202111310736A CN 114065493 A CN114065493 A CN 114065493A
- Authority
- CN
- China
- Prior art keywords
- program
- time
- abnormal
- reactor system
- calculation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Withdrawn
Links
- 238000000034 method Methods 0.000 title claims abstract description 22
- 230000002159 abnormal effect Effects 0.000 claims abstract description 39
- 238000012544 monitoring process Methods 0.000 claims abstract description 18
- 238000005457 optimization Methods 0.000 claims abstract description 5
- 238000004364 calculation method Methods 0.000 claims description 20
- 230000000704 physical effect Effects 0.000 claims description 3
- 239000011800 void material Substances 0.000 claims description 3
- 230000003247 decreasing effect Effects 0.000 claims description 2
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000004422 calculation algorithm Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 230000006399 behavior Effects 0.000 description 1
- 230000005494 condensation Effects 0.000 description 1
- 238000009833 condensation Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000008020 evaporation Effects 0.000 description 1
- 238000001704 evaporation Methods 0.000 description 1
- 238000002789 length control Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000005514 two-phase flow Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
- G21C—NUCLEAR REACTORS
- G21C17/00—Monitoring; Testing ; Maintaining
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/30—Nuclear fission reactors
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Plasma & Fusion (AREA)
- High Energy & Nuclear Physics (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
The invention relates to the technical field of reactor thermal hydraulic power and safety analysis, and particularly discloses a method and a system for optimizing stability of an analysis program of a reactor system. The method comprises the following steps: s1, distributing a plurality of sections of memory areas, and storing global parameters before the current computing time; s2, monitoring the running state of the program, and returning the program calculation state to the time point nearest to the abnormal time when the abnormal running of the program is monitored; and S3, reducing the time step of the program, restarting the program by using the global parameter at the latest moment, if the program runs abnormally, sequentially selecting a more previous time point, reducing the time step again sequentially, restarting by using the new global parameter, and gradually increasing the reduced time step sequentially until the time step returns to the initial set value after normal monitoring. The method can obviously increase the stability of the system program, realizes the stability optimization of the program, and is beneficial to improving the efficiency of safety analysis.
Description
Technical Field
The invention belongs to the technical field of reactor thermal hydraulic power and safety analysis, and particularly relates to a method and a system for optimizing stability of an analysis program of a reactor system.
Background
The nuclear reactor safety analysis is mainly used for modeling a nuclear power plant system by utilizing a safety analysis program and carrying out related analysis and calculation. The nuclear reactor accident phenomenon is complex, and can involve the severe change process of parameters such as phase change, evaporation and condensation, and the severe change of the parameters can cause the unstable calculation process of the analysis program of the reactor system, and even the problems of program error report or deviation of the calculation result from the real phenomenon occur. The stability of the program directly determines the accuracy and the analysis efficiency of the safety analysis of the reactor.
Due to the complexity of the physical phenomena themselves during a reactor accident, severe fluctuations of the parameters during the simulation cannot be avoided. Although a great deal of research is carried out on two-phase flow algorithms at home and abroad, an absolutely stable and reliable algorithm framework is sought to solve the problem that a system program is complicated and calculation errors occur, but the absolutely stable algorithm framework is not found so far. At present, time step length is usually manually corrected in safety analysis, and the problem is solved based on a breakpoint restart mode, and the mode is not beneficial to checking calculation errors deviating from a real phenomenon on one hand, and greatly influences accident analysis efficiency on the other hand.
Disclosure of Invention
The invention aims to provide a method and a system for optimizing the stability of an analysis program of a reactor system, which solve the problem of the computational stability of the analysis process of the reactor system.
The technical scheme of the invention is as follows: a method for stable optimization of an analysis program of a reactor system, the method comprising the steps of:
s1, distributing a plurality of sections of memory areas in the analysis program of the reactor system, and storing global parameters of a plurality of characteristic state points before the current calculation time;
s2, monitoring the running state of the analysis program of the reactor system, and returning the program calculation state to the time point nearest to the abnormal time when the abnormal running of the program is monitored;
s3, reducing the time step of the analysis program of the reactor system, restarting the program by using the global parameter at the latest moment, if the abnormal operation of the program is continuously detected, sequentially selecting a more previous characteristic state time point, reducing the time step again in sequence, restarting the program by using the global parameter at the characteristic state time point, and gradually increasing the reduced time step until the time step returns to the initial set value after the program is monitored to be calculated normally.
The step S3 specifically includes:
s3.1, reducing the time step length to 1/A of the original step length, and automatically restarting a program by using the global parameters stored in the memory of the time point, wherein A is a coefficient more than 1;
s3.2, monitoring whether the program is abnormal near the original breakpoint time, if so, returning the program calculation state to the previous characteristic state point time, reducing the time step length to 1/A of the original step length, restarting the program by using the global parameters of the corresponding characteristic state point time, and repeating the steps until the program is normally calculated;
and S3.3, monitoring the time T of the physical problem calculated by the analysis program of the reactor system, if the time T of the calculated physical problem is greater than the threshold time delta T, multiplying the reduced time step length by A in sequence until the time step length is restored to the initial set value, finishing the program breakpoint continued calculation, and continuing to perform subsequent calculation.
The step 1 specifically comprises:
in a system analysis program, a plurality of sections of memory areas are distributed by adding a multi-dimensional array variable, and the memory of each section of memory area is used for storing global parameters of a plurality of characteristic state points before the current moment.
The global parameters include pressure, temperature, flow rate, void fraction, non-condensable gas fraction, temperature of thermal build-up, and neutronics parameters of all control volumes.
The specific steps of determining the abnormal time in step S2 are as follows:
monitoring the running state of an analysis program of a reactor system, and judging the abnormal time when a breakpoint or an abnormal time deviated from the real condition occurs in the running of the program when a physical property variable in a certain control body is out of range, a non-numerical value condition occurs in a certain parameter or a dividend is 0.
The time step reduction scaling factor a is 1.9.
The threshold time Δ T is the sum of the abnormal time Ter and a set time dT.
A reactor system program stability optimization system comprises a temporary memory allocation module, a breakpoint and calculation abnormity monitoring module and a time step control module;
the temporary memory allocation module allocates a plurality of memory areas in a multi-dimensional array variable increasing mode; the memory area sequentially stores global parameters of a plurality of characteristic state point moments before the current calculation moment;
the breakpoint and calculation abnormity monitoring module can monitor the running state of an analysis program of the reactor system, record the time of the abnormal running state when the program generates an abnormal time, and return the program calculation state to the time point nearest to the abnormal time;
the time step control module can adjust and reduce the time step when receiving a program operation abnormal signal, and restart the program by using the time point nearest to the abnormal time; if the program operation abnormal signal is continuously received, continuously selecting a more previous characteristic state time point, adjusting and reducing the time step length, restarting the program by using the global parameter of the latest characteristic state time point, and repeating the steps until the program operation abnormal signal is not received, and gradually adjusting and increasing the time step length to the initial set value.
The time step control module adjusts the decreasing time step and the increasing time step to be reciprocal to each other.
And the time step length control module adjusts and increases the time step length after judging that the time for calculating the physical problem by the program is more than the threshold time.
The invention has the following remarkable effects: the stability optimization method and the system for the reactor system analysis program can obviously improve the stability of the system program, particularly the stability of the program in complex behaviors.
Drawings
Fig. 1 is a schematic diagram of a method for optimizing stability of an analysis program of a reactor system according to the present invention.
Detailed Description
The invention is described in further detail below with reference to the figures and the embodiments.
As shown in fig. 1, a method for optimizing stability of an analysis program of a reactor system specifically includes the following steps:
s1, distributing a plurality of sections of memory areas in the analysis program of the reactor system, and storing global parameters of a plurality of characteristic state points before the current calculation time;
in a system analysis program, distributing a plurality of sections of memory areas by adding a multi-dimensional array variable, wherein the memory of each section of memory area is used for storing global parameters of a plurality of characteristic state points before the current moment, and the global parameters comprise pressure, temperature, flow rate, void fraction, non-condensable gas fraction, temperature of heat construction and neutron parameters of all control bodies; for example, the multidimensional array variable may store global parameters at a time 5 feature state points before the current time;
s2, monitoring the running state of the analysis program of the reactor system, and returning the program calculation state to the time point nearest to the abnormal time when the abnormal running of the program is monitored;
monitoring the running state of an analysis program of a reactor system, judging the breakpoint or abnormal time deviating from the real situation when a physical property variable in a certain control body is out of range, a certain parameter has a non-numerical value condition or the dividend has a value of 0, recording the time Ter at the moment, and returning the program calculation state to the time point nearest to the abnormal time;
s3, reducing the time step of the analysis program of the reactor system, restarting the program by using the global parameter at the latest moment, if the abnormal operation of the program is continuously detected, sequentially selecting a more previous characteristic state time point, reducing the time step again in sequence, restarting the program by using the global parameter at the characteristic state time point, and gradually increasing the reduced time step until the time step returns to the initial set value after the program is monitored to be calculated normally;
s3.1, reducing the time step length to 1/1.9 of the original step length, and automatically restarting a program by using the global parameters stored in the memory of the time point;
s3.2, monitoring whether the program is abnormal near the original breakpoint time, if so, returning the program calculation state to the previous characteristic state point time, reducing the time step length to 1/1.9 of the original step length, then restarting the program by using the global parameters of the corresponding characteristic state point time, and repeating the steps until the program is normally calculated;
and S3.3, monitoring the time T of the physical problem calculated by the reactor system analysis program, if the calculated physical problem time T is greater than a threshold time delta T, sequentially multiplying the reduced time step by 1.9, and completing program breakpoint continued calculation and continuing to perform subsequent calculation until the time step is restored to an initial set value, wherein the threshold time delta T is the sum of the abnormal time Ter and a set time dT.
Claims (10)
1. A method for optimizing stability of an analysis program of a reactor system is characterized by comprising the following steps:
s1, distributing a plurality of sections of memory areas in the analysis program of the reactor system, and storing global parameters of a plurality of characteristic state points before the current calculation time;
s2, monitoring the running state of the analysis program of the reactor system, and returning the program calculation state to the time point nearest to the abnormal time when the abnormal running of the program is monitored;
s3, reducing the time step of the analysis program of the reactor system, restarting the program by using the global parameter at the latest moment, if the abnormal operation of the program is continuously detected, sequentially selecting a more previous characteristic state time point, reducing the time step again in sequence, restarting the program by using the global parameter at the characteristic state time point, and gradually increasing the reduced time step until the time step returns to the initial set value after the program is monitored to be calculated normally.
2. The method for optimizing stability of the analysis program of the reactor system according to claim 1, wherein the step S3 specifically includes:
s3.1, reducing the time step length to 1/A of the original step length, and automatically restarting a program by using the global parameters stored in the memory of the time point, wherein A is a coefficient more than 1;
s3.2, monitoring whether the program is abnormal near the original breakpoint time, if so, returning the program calculation state to the previous characteristic state point time, reducing the time step length to 1/A of the original step length, restarting the program by using the global parameters of the corresponding characteristic state point time, and repeating the steps until the program is normally calculated;
and S3.3, monitoring the time T of the physical problem calculated by the analysis program of the reactor system, if the time T of the calculated physical problem is greater than the threshold time delta T, multiplying the reduced time step length by A in sequence until the time step length is restored to the initial set value, finishing the program breakpoint continued calculation, and continuing to perform subsequent calculation.
3. The method for optimizing stability of the analysis program of the reactor system according to claim 1, wherein the step 1 specifically comprises:
in a system analysis program, a plurality of sections of memory areas are distributed by adding a multi-dimensional array variable, and the memory of each section of memory area is used for storing global parameters of a plurality of characteristic state points before the current moment.
4. The method of claim 1, wherein the global parameters include pressure, temperature, flow rate, void fraction, non-condensable gas fraction, temperature of thermal build, and neutronics parameters of all control volumes.
5. The method for optimizing stability of the analysis program of the reactor system according to claim 1, wherein the step S2 of determining the abnormal time specifically comprises:
monitoring the running state of an analysis program of a reactor system, and judging the abnormal time when a breakpoint or an abnormal time deviated from the real condition occurs in the running of the program when a physical property variable in a certain control body is out of range, a non-numerical value condition occurs in a certain parameter or a dividend is 0.
6. The method of claim 2, wherein the time step reduction scaling factor A is 1.9.
7. The method as claimed in claim 2, wherein the threshold time Δ T is the abnormal time Ter plus a set time dT.
8. A reactor system program stability optimization system is characterized by comprising a temporary memory allocation module, a breakpoint and calculation abnormity monitoring module and a time step control module;
the temporary memory allocation module allocates a plurality of memory areas in a multi-dimensional array variable increasing mode; the memory area sequentially stores global parameters of a plurality of characteristic state point moments before the current calculation moment;
the breakpoint and calculation abnormity monitoring module can monitor the running state of an analysis program of the reactor system, record the time of the abnormal running state when the program generates an abnormal time, and return the program calculation state to the time point nearest to the abnormal time;
the time step control module can adjust and reduce the time step when receiving a program operation abnormal signal, and restart the program by using the time point nearest to the abnormal time; if the program operation abnormal signal is continuously received, continuously selecting a more previous characteristic state time point, adjusting and reducing the time step length, restarting the program by using the global parameter of the latest characteristic state time point, and repeating the steps until the program operation abnormal signal is not received, and gradually adjusting and increasing the time step length to the initial set value.
9. The system of claim 8, wherein the step size control module adjusts the decreasing step size and the increasing step size to be reciprocal to each other.
10. The system of claim 8, wherein the step control module increases the step by a time determined by the time to calculate the physical problem for the program greater than a threshold time.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111310736.5A CN114065493A (en) | 2021-11-04 | 2021-11-04 | Method and system for optimizing stability of analysis program of reactor system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111310736.5A CN114065493A (en) | 2021-11-04 | 2021-11-04 | Method and system for optimizing stability of analysis program of reactor system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114065493A true CN114065493A (en) | 2022-02-18 |
Family
ID=80274470
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111310736.5A Withdrawn CN114065493A (en) | 2021-11-04 | 2021-11-04 | Method and system for optimizing stability of analysis program of reactor system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114065493A (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109002624A (en) * | 2018-07-26 | 2018-12-14 | 上海交通大学 | Supersonic speed rigidity combustion flows double adaptive decouples Optimized Simulated method and system |
| CN109189607A (en) * | 2018-08-13 | 2019-01-11 | 南京翼辉信息技术有限公司 | A kind of method that application program breakpoint restores |
| CN113408147A (en) * | 2021-07-15 | 2021-09-17 | 中国科学院近代物理研究所 | Reactor fuel performance analysis and calculation method, system, storage medium and equipment |
-
2021
- 2021-11-04 CN CN202111310736.5A patent/CN114065493A/en not_active Withdrawn
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109002624A (en) * | 2018-07-26 | 2018-12-14 | 上海交通大学 | Supersonic speed rigidity combustion flows double adaptive decouples Optimized Simulated method and system |
| CN109189607A (en) * | 2018-08-13 | 2019-01-11 | 南京翼辉信息技术有限公司 | A kind of method that application program breakpoint restores |
| CN113408147A (en) * | 2021-07-15 | 2021-09-17 | 中国科学院近代物理研究所 | Reactor fuel performance analysis and calculation method, system, storage medium and equipment |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112855408B (en) | Early warning method and early warning device for hydroelectric generating set | |
| CN104811344A (en) | Network dynamic service monitoring method and apparatus | |
| CN112083698B (en) | Production line real-time control method and device based on digital twin model | |
| CN108980628A (en) | A kind of water supply network booster alarm method | |
| US20240003122A1 (en) | Boom correction method and device for working machine | |
| CN109492974A (en) | The more Weibull assembly of elements spare parts demand amounts of larger cargo ships entirety alternate maintenance determine method | |
| US6477218B1 (en) | Control system of nuclear power plant, and control method thereof | |
| CN114065493A (en) | Method and system for optimizing stability of analysis program of reactor system | |
| WO2022267209A1 (en) | Parallel valve set control method and apparatus, controller, and storage medium | |
| CN116293410B (en) | Analysis method and device based on intelligent electronic scale | |
| US9099207B2 (en) | System, method, and program for monitoring reactor core | |
| WO2024087422A1 (en) | Power tilt prediction method and apparatus, device, storage medium, and program product | |
| CN116246802A (en) | Reactor control method, apparatus, computer device and storage medium | |
| CN112736935B (en) | Method for online checking model parameters of power system stabilizer by using PSS compensation angle | |
| CN115410731A (en) | Reactor entering feasibility analysis method, device and equipment for repairing assembly in reactor | |
| CN114896102A (en) | Data protection time point recovery method and system | |
| CN110162853B (en) | Mechanical experiment data preprocessing method based on mechanical model | |
| CN111794813B (en) | Method and device for monitoring operation performance of steam turbine and electronic equipment | |
| CN111445074B (en) | Low-temperature economizer energy saving optimization method, device, equipment and storage medium | |
| CN110763951B (en) | Fastest change detection model construction method for power transmission line power failure detection and identification | |
| CN114673603B (en) | Engine control system safety monitoring method, device, computer equipment and medium | |
| CN116931437B (en) | Control method and system of biomass gasification carbonization furnace | |
| CN117584790B (en) | Capacity-free charging pile control system | |
| CN118170124B (en) | Fault analysis method, device, equipment and storage medium of instrument control system | |
| CN117592355A (en) | Start-stop control method and device for metal stack, computer equipment and storage medium |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WW01 | Invention patent application withdrawn after publication | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20220218 |