CN111651851A - Containment solving method and containment solver - Google Patents
Containment solving method and containment solver Download PDFInfo
- Publication number
- CN111651851A CN111651851A CN201910160710.3A CN201910160710A CN111651851A CN 111651851 A CN111651851 A CN 111651851A CN 201910160710 A CN201910160710 A CN 201910160710A CN 111651851 A CN111651851 A CN 111651851A
- Authority
- CN
- China
- Prior art keywords
- phase
- liquid phase
- gas
- term
- containment
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 239000012071 phase Substances 0.000 claims abstract description 160
- 239000007791 liquid phase Substances 0.000 claims abstract description 149
- 238000004134 energy conservation Methods 0.000 claims abstract description 40
- 238000004364 calculation method Methods 0.000 claims abstract description 33
- 230000001052 transient effect Effects 0.000 claims description 18
- 230000005484 gravity Effects 0.000 claims description 12
- 239000011159 matrix material Substances 0.000 claims description 12
- 230000001133 acceleration Effects 0.000 claims description 6
- 239000011800 void material Substances 0.000 claims description 6
- 239000007788 liquid Substances 0.000 claims 1
- 230000008569 process Effects 0.000 abstract description 12
- 239000012530 fluid Substances 0.000 abstract description 7
- 238000010586 diagram Methods 0.000 description 7
- 230000004888 barrier function Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 101100087530 Caenorhabditis elegans rom-1 gene Proteins 0.000 description 1
- 101100305983 Mus musculus Rom1 gene Proteins 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
Images
Abstract
The invention discloses a containment solving method and a containment solver. The containment solving method comprises the following steps: reading relevant information of the model; solving the momentum conservation equation of the gas phase and the liquid phase according to the related information to obtain the flow rate of the gas phase and the liquid phase; respectively judging whether the flow velocity directions of the gas phase and the liquid phase obtained by two adjacent calculations are consistent; if the two variables are consistent, the mass conservation equations of the steam, the non-condensable gas and the liquid phase and the energy conservation equations of the gas phase and the liquid phase are simultaneously established to obtain the values of the main variable and other variables; and correcting variables such as flow rate. The containment solution method provided by the invention adopts the multi-phase multi-flow field model, takes heat and mass transfer processes and phenomena among phases, wall surfaces and inside fluid into consideration, and discretely solves the mass conservation equation, the momentum equation and the energy conservation equation, so that the containment analysis program can calculate the thermal hydraulic process of the containment, and the containment solution method has the advantages of accuracy and reliability.
Description
Technical Field
The invention relates to the technical field of nuclear power, in particular to a containment solving method and a containment solver.
Background
As a third barrier for nuclear power plants, containment is also the last barrier of the defense-in-depth system, and its integrity and reliability are crucial to the safety of the nuclear power plant. The containment analysis program is a program for analyzing and calculating the thermal hydraulic response and integrity of the containment and is an important tool for containment thermal hydraulic design and safety analysis, so that the accuracy and reliability of the containment analysis program are very important.
Disclosure of Invention
The present invention is directed to solving, at least to some extent, one of the technical problems in the related art described above.
To this end, a first object of the present invention is to propose a containment solution method. The method adopts a multi-phase multi-flow field model, takes heat and mass transfer processes and phenomena among phases, wall surfaces and inside fluid into consideration, and carries out discrete solution on a mass conservation equation, a momentum equation and an energy conservation equation, so that the containment analysis program can calculate the thermal hydraulic process of the containment, and the method has the advantages of accuracy and reliability.
The second purpose of the invention is to provide a containment solver.
In order to achieve the above object, a first aspect of the present invention discloses a containment solution method, including the steps of: reading relevant information of the model; solving the momentum conservation equation of the gas phase and the liquid phase according to the related information to obtain the flow rate of the gas phase and the liquid phase; respectively judging whether the flow velocity directions of the gas phase and the liquid phase obtained by two adjacent calculations are consistent; if the two variables are consistent, the mass conservation equations of the steam, the non-condensable gas and the liquid phase and the energy conservation equations of the gas phase and the liquid phase are simultaneously established to obtain the values of the main variable and other variables; the values of variables such as flow rate are corrected.
According to the containment vessel solving method, a multi-phase multi-flow field model is adopted, heat and mass transfer processes and phenomena among phases, on wall surfaces and in fluid are considered, a mass conservation equation, a momentum equation and an energy conservation equation are subjected to discrete solution, calculation of a containment vessel analysis program on a thermal hydraulic process of a containment vessel is achieved, and the containment vessel solving method has the advantages of accuracy and reliability.
In some examples, the relevant information includes a model category, a model parameter value, an initial condition, and a boundary condition.
In some examples, solving by the conservation of momentum equation for the gas and liquid phases based on the correlation information to obtain the gas and liquid phase flow rates comprises: calculating terms of a gas phase and liquid phase momentum equation, wherein the gas phase and liquid phase momentum equation is of the form:
wherein each item of the gas phase and liquid phase momentum equation comprises a transient term, a convection term, a pressure term, a gravity term and a resistance term, wherein k is the liquid phase or the gas phase, αkIs the void fraction of the k phase, pkIs the density of the k phase, ukIs the velocity of the k phase, t is the time, p is the pressure, g is the acceleration of gravity, F'kResistance to k-phase;
and solving the momentum equations of the gas phase and the liquid phase to obtain the flow rates of the gas phase and the liquid phase.
In some examples, the simultaneous vapor, non-condensable gas and liquid phase conservation of mass equation and gas and liquid phase conservation of energy equation yields values of principal and other variables, including:
and calculating the mass conservation equations of the steam, the non-condensable gas and the liquid phase and the energy conservation equations of the gas phase and the liquid phase, wherein the mass conservation equations of the steam, the non-condensable gas and the liquid phase are in the form of:
wherein, each item of mass conservation equation of steam, non-condensable gas and liquid phase comprises transient term, convection term, mass source term and the like, wherein k is liquid phase, steam or non-condensable gas'kIs a mass source term of k phases;
the gas phase and liquid phase energy conservation equation is in the form of:
wherein each item of the gas phase and liquid phase energy conservation equation comprises a transient term, a convection term, a pressure work-doing term, a heat transfer term and the like, wherein k is a liquid phase or a gas phase, hkIs the specific enthalpy of the k phase, Φ'kIs the heat source term of the k phase, Q'kIs the heat transfer term of the k phase;
and simultaneously establishing a mass conservation equation of steam, noncondensable gas and liquid phase and an energy conservation equation of gas phase and liquid phase, solving based on a Newton-Raphson algorithm to obtain a full-field solution matrix, and calculating according to the full-field solution matrix to obtain main variables such as pressure and enthalpy and values of other variables.
In some examples, when the directions of the gas-phase flow velocity and the liquid-phase flow velocity obtained by two adjacent calculations are judged to be inconsistent, the momentum conservation equation is solved again until the directions of the flow velocities obtained by two adjacent calculations are consistent.
An embodiment of a second aspect of the present invention discloses a containment solver, including: the initialization module is used for reading in relevant information of the model; the solving module is used for solving the momentum conservation equation of the gas phase and the liquid phase according to the related information to obtain the flow velocity of the gas phase and the liquid phase, judging whether the flow velocity directions of the gas phase and the liquid phase obtained by two adjacent calculations are consistent, if so, combining the mass conservation equation of the steam, the non-condensable gas and the liquid phase and the energy conservation equation of the gas phase and the liquid phase to obtain the values of a principal variable and other variables, and correcting the values of the variables such as the flow velocity; and the output module is used for outputting the solving result.
According to the containment solver, a multiphase multi-flow field model is adopted, heat and mass transfer processes and phenomena among phases, wall surfaces and inside fluid are considered, a mass conservation equation, a momentum equation and an energy conservation equation are subjected to discrete solution, calculation of a containment analysis program on a thermal hydraulic process of a containment is achieved, and the containment solver has the advantages of accuracy and reliability.
In some examples, the relevant information includes a model category, a model parameter value, an initial condition, and a boundary condition.
In some examples, the solution module is to:
calculating terms of a gas phase and liquid phase momentum equation, wherein the gas phase and liquid phase momentum equation is of the form:
wherein each item of the gas phase and liquid phase momentum equation comprises a transient term, a convection term, a pressure term, a gravity term and a resistance term, wherein k is the liquid phase or the gas phase, αkIs the void fraction of the k phase, pkIs the density of the k phase, ukIs the velocity of the k phase, t is the time, p is the pressure, g is the acceleration of gravity, F'kResistance to k-phase;
and solving the momentum equations of the gas phase and the liquid phase to obtain the flow rates of the gas phase and the liquid phase.
In some examples, the solution module is to:
and calculating the mass conservation equations of the steam, the non-condensable gas and the liquid phase and the energy conservation equations of the gas phase and the liquid phase, wherein the mass conservation equations of the steam, the non-condensable gas and the liquid phase are in the form of:
wherein, each item of mass conservation equation of steam, non-condensable gas and liquid phase comprises transient term, convection term, mass source term and the like, wherein k is liquid phase, steam or non-condensable gas'kIs a mass source term of k phases;
the gas phase and liquid phase energy conservation equation is in the form of:
wherein each item of the gas phase and liquid phase energy conservation equation comprises a transient term, a convection term, a pressure work-doing term, a heat transfer term and the like, wherein k is a liquid phase or a gas phase, hkIs the specific enthalpy of the k phase, Φ'kIs the heat source term of the k phase, Q'kIs the heat transfer term of the k phase;
and simultaneously establishing a mass conservation equation of steam, noncondensable gas and liquid phase and an energy conservation equation of gas phase and liquid phase, solving based on a Newton-Raphson algorithm to obtain a full-field solution matrix, and calculating according to the full-field solution matrix to obtain main variables such as pressure and enthalpy and values of other variables.
In some examples, the solving module is configured to, when it is determined that the directions of the gas-phase flow velocity and the liquid-phase flow velocity obtained by two adjacent calculations are not consistent, re-solve the momentum conservation equation until the directions of the flow velocities obtained by two adjacent calculations are consistent.
Additional aspects and advantages of the invention will be set forth in part in the description which follows and, in part, will be obvious from the description, or may be learned by practice of the invention.
Drawings
The foregoing and additional aspects and advantages of the present invention will become apparent and readily appreciated from the following description of the embodiments, taken in conjunction with the accompanying drawings of which:
FIG. 1 is a flow diagram of a containment solution method according to one embodiment of the invention;
FIG. 2 is a flow diagram of a containment solver, according to one embodiment of the invention;
FIG. 3 is a grid schematic of a containment solution method according to one embodiment of the invention;
FIG. 4 is a schematic diagram illustrating a comparison of containment solution methods with the results of compartment temperature calculations of the related art according to an embodiment of the present invention;
FIG. 5 is a schematic diagram illustrating a comparison of containment solution methods with the calculated results of compartment pressures of the related art according to an embodiment of the present invention;
FIG. 6 is a block diagram of a containment solver, according to one embodiment of the invention.
Detailed Description
Reference will now be made in detail to embodiments of the present invention, examples of which are illustrated in the accompanying drawings, wherein like or similar reference numerals refer to the same or similar elements or elements having the same or similar function throughout. The embodiments described below with reference to the accompanying drawings are illustrative only for the purpose of explaining the present invention, and are not to be construed as limiting the present invention.
The containment solution method and the containment solver according to the embodiments of the present invention are described below with reference to the accompanying drawings.
FIG. 1 is a flow diagram of a containment solution method according to one embodiment of the invention.
As shown in FIG. 1 in conjunction with FIG. 2, a containment solution method according to an embodiment of the invention includes the following steps:
s101: and reading relevant information of the model.
Wherein the related information includes but is not limited to: model type, model parameter values, initial conditions, boundary conditions, and the like.
Specifically, when an input card is read, if the input card has an error, the error is reported, and error report information is provided for a user, and if the input card has no error, program variables are assigned according to the content of the input card, such as model type selection, model parameters, initial conditions, boundary conditions, and the like.
S102: and solving the momentum conservation equation of the gas phase and the liquid phase according to the related information to obtain the flow rates of the gas phase and the liquid phase.
For example: calculating the terms of the gas and liquid phase conservation of momentum equations, wherein the gas and liquid phase conservation of momentum equations are of the form:
wherein each term of the gas phase and liquid phase momentum conservation equation comprises a transient term, a convection term, a pressure term, a gravity term and a resistance term, wherein k is the liquid phase or the gas phase, αkIs the void fraction of the k phase, pkIs the density of the k phase, ukIs the velocity of the k phase, t is the time, p is the pressure, g is the acceleration of gravity, F'kResistance to k-phase;
and solving the gas-phase and liquid-phase momentum conservation equations to obtain the gas-phase and liquid-phase flow rates.
Namely: and calculating each item of the gas-phase and liquid-phase momentum conservation equations, solving the gas-phase and liquid-phase momentum conservation equations, and calculating to obtain the gas-phase and liquid-phase flow rates.
S103: and respectively judging whether the flow velocity directions of the gas phase and the liquid phase obtained by two adjacent calculations are consistent.
When the direction of the gas phase flow velocity and the liquid phase flow velocity obtained by two adjacent times of calculation is judged to be inconsistent, the momentum conservation equation is solved again until the flow velocity direction obtained by two adjacent times of calculation is consistent. Namely: and judging the speed direction, if the speed direction is opposite to the speed direction obtained by the last calculation, judging the upstream again according to a new speed value, and solving the momentum equation again until the speed directions obtained by the previous calculation and the next calculation are the same.
S104: if the two variables are consistent, the mass conservation equations of the steam, the non-condensable gas and the liquid phase and the energy conservation equations of the gas phase and the liquid phase are combined to obtain the values of the main variable and other variables.
For example: and calculating the mass conservation equations of the steam, the non-condensable gas and the liquid phase and the energy conservation equations of the gas phase and the liquid phase, wherein the mass conservation equations of the steam, the non-condensable gas and the liquid phase are in the form of:
wherein, each item of mass conservation equation of steam, non-condensable gas and liquid phase comprises transient term, convection term, mass source term and the like, wherein k is liquid phase, steam or non-condensable gas'kIs a mass source term of k phases;
the gas phase and liquid phase energy conservation equation is in the form of:
wherein each item of the gas phase and liquid phase energy conservation equation comprises a transient term, a convection term, a pressure work-doing term, a heat transfer term and the like, wherein k is a liquid phase or a gas phase, hkIs the specific enthalpy of the k phase, Φ'kIs the heat source term of the k phase, Q'kIs the heat transfer term of the k phase;
and simultaneously solving a mass conservation equation of the steam, the non-condensable gas and the liquid phase and an energy conservation equation of the gas phase and the liquid phase to obtain main variables such as pressure, enthalpy and the like and values of other variables.
Namely: calculating each item of a mass conservation equation and an energy conservation equation, then combining the mass conservation equation of steam, non-condensable gas and liquid phase and the energy conservation equation of gas phase and liquid phase, solving based on a Newton-Raphson algorithm to obtain a solution matrix of a full field, and then calculating according to the solution matrix of the full field to obtain main variables such as pressure, enthalpy and the like and values of other variables.
S105: the values of variables such as speed are corrected.
Specifically, the updated pressure value is used to correct variables such as the flow rate. Then, judging whether the calculation is converged according to a convergence criterion, if so, outputting the calculation result of the time step, and starting the calculation of the next time step; and if not, returning to start the next loop iterative computation.
Taking a simple example of simulating a containment mass-energy release accident as an example, as shown in FIG. 3, the calculationExample includes the boundary BC, the compartment rom 1. From 1 having a cross-sectional area of 100m2And the height is 2 m. The initial conditions include: the Room1 was filled with air at a pressure of 0.1MPa and a temperature of 400K. The boundary conditions are as follows: BC is the flow boundary, the boundary fluid is steam, the temperature is equal to 400K, and the flow is 1 kg/s.
Firstly, writing relevant information of a calculation example, including geometric parameters, initial conditions, boundary conditions, model parameters and the like, into an input card, solving according to a momentum equation to obtain flow rate, then substituting the flow rate into a mass conservation equation and an energy conservation equation, simultaneously solving to obtain variables such as temperature and pressure of a compartment, and comparing a calculation result with a containment program GOTHIC calculation result which is commonly used internationally, as shown in fig. 4 and 5.
According to the containment vessel solving method, a multi-phase multi-flow field model is adopted, heat and mass transfer processes and phenomena among phases, on wall surfaces and in fluid are considered, a mass conservation equation, a momentum equation and an energy conservation equation are subjected to discrete solution, calculation of a containment vessel analysis program on a thermal hydraulic process of a containment vessel is achieved, and the containment vessel solving method has the advantages of accuracy and reliability.
FIG. 6 is a block diagram of a containment solver, according to one embodiment of the invention. As shown in FIG. 6, a containment solver 600 according to an embodiment of the invention includes: an initialization module 610, a solving module 620, and an output module 630.
The initialization module 610 is configured to read in relevant information of the model; the solving module 620 is used for solving the momentum conservation equations of the gas phase and the liquid phase according to the related information to obtain the flow rates of the gas phase and the liquid phase, judging whether the flow rates of the gas phase and the liquid phase obtained by two adjacent calculations are consistent in direction, and if so, simultaneously solving the mass conservation equations of the steam, the non-condensable gas and the liquid phase and the energy conservation equations of the gas phase and the liquid phase to obtain the values of the principal variables and other variables and correcting the values of the variables such as the flow rate; the output module 630 is used for outputting the solution result.
In one embodiment of the invention, the relevant information comprises model class, model parameter values, initial conditions and boundary conditions.
In one embodiment of the invention, the solving module is configured to:
calculating terms of a gas phase and liquid phase momentum equation, wherein the gas phase and liquid phase momentum equation is of the form:
wherein each item of the gas phase and liquid phase momentum equation comprises a transient term, a convection term, a pressure term, a gravity term and a resistance term, wherein k is the liquid phase or the gas phase, αkIs the void fraction of the k phase, pkIs the density of the k phase, ukIs the velocity of the k phase, t is the time, p is the pressure, g is the acceleration of gravity, F'kResistance to k-phase;
and solving the momentum equations of the gas phase and the liquid phase to obtain the flow rates of the gas phase and the liquid phase.
In one embodiment of the invention, the solving module is configured to:
and calculating the mass conservation equations of the steam, the non-condensable gas and the liquid phase and the energy conservation equations of the gas phase and the liquid phase, wherein the mass conservation equations of the steam, the non-condensable gas and the liquid phase are in the form of:
wherein, each item of mass conservation equation of steam, non-condensable gas and liquid phase comprises transient term, convection term, mass source term and the like, wherein k is liquid phase, steam or non-condensable gas'kIs a mass source term of k phases;
the gas phase and liquid phase energy conservation equation is in the form of:
wherein each item of the gas phase and liquid phase energy conservation equation comprises a transient term, a convection term, a pressure work-doing term, a heat transfer term and the like, wherein k is a liquid phase or a gas phase, hkIs the specific enthalpy of the k phase, Φ'kIs the heat source term of the k phase, Q'kIs the heat transfer term of the k phase;
and simultaneously establishing a mass conservation equation of steam, noncondensable gas and liquid phase and an energy conservation equation of gas phase and liquid phase, solving based on a Newton-Raphson algorithm to obtain a full-field solution matrix, and calculating according to the full-field solution matrix to obtain main variables such as pressure and enthalpy and values of other variables.
In an embodiment of the invention, the solving module is configured to, when it is determined that the directions of the gas-phase flow velocity and the liquid-phase flow velocity obtained through two adjacent calculations are not consistent, re-solve the momentum conservation equation until the directions of the flow velocities obtained through two adjacent calculations are consistent.
According to the containment solver, a multiphase multi-flow field model is adopted, heat and mass transfer processes and phenomena among phases, wall surfaces and inside fluid are considered, a mass conservation equation, a momentum equation and an energy conservation equation are subjected to discrete solution, calculation of a containment analysis program on a thermal hydraulic process of a containment is achieved, and the containment solver has the advantages of accuracy and reliability.
It should be noted that a specific implementation manner of the containment solver in the embodiment of the present invention is similar to a specific implementation manner of the containment solution method in the embodiment of the present invention, and reference is specifically made to the description of the method portion, which is not described herein again.
In the description herein, references to the description of the term "one embodiment," "some embodiments," "an example," "a specific example," or "some examples," etc., mean that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the invention. In this specification, the schematic representations of the terms used above are not necessarily intended to refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Furthermore, various embodiments or examples and features of different embodiments or examples described in this specification can be combined and combined by one skilled in the art without contradiction.
In the description of the present invention, it should be noted that, unless otherwise specified and limited, the terms "mounted," "connected," and "connected" are used in a broad sense, and may be, for example, mechanically or electrically connected, or may be connected by two elements, directly or indirectly through an intermediate medium, and the specific meaning of the terms may be understood by those skilled in the art according to specific situations
Although embodiments of the present invention have been shown and described above, it is understood that the above embodiments are exemplary and should not be construed as limiting the present invention, and that variations, modifications, substitutions and alterations can be made to the above embodiments by those of ordinary skill in the art within the scope of the present invention.
Claims (10)
1. A containment solution method is characterized by comprising the following steps:
reading relevant information of the model;
solving the momentum conservation equation of the gas phase and the liquid phase according to the related information to obtain the flow rate of the gas phase and the liquid phase;
respectively judging whether the flow velocity directions of the gas phase and the liquid phase obtained by two adjacent calculations are consistent;
if the two variables are consistent, the mass conservation equations of the steam, the non-condensable gas and the liquid phase and the energy conservation equations of the gas phase and the liquid phase are simultaneously established to obtain the values of the main variable and other variables;
and correcting variables such as flow rate.
2. The containment solution method according to claim 1, wherein the relevant information includes model type, model parameter values, initial conditions, and boundary conditions.
3. The containment solution method according to claim 1, wherein solving the conservation of momentum equations for the gas and liquid phases based on the correlation information to obtain the gas and liquid phase flow rates comprises:
calculating terms of a gas phase and liquid phase momentum equation, wherein the gas phase and liquid phase momentum equation is of the form:
wherein each item of the gas phase and liquid phase momentum equation comprises a transient term, a convection term, a pressure term, a gravity term and a resistance term, wherein k is the liquid phase or the gas phase, αkIs the void fraction of the k phase, pkIs the density of the k phase, ukIs the velocity of the k phase, t is the time, p is the pressure, g is the acceleration of gravity, F'kResistance to k-phase;
and solving the momentum equations of the gas phase and the liquid phase to obtain the flow rates of the gas phase and the liquid phase.
4. The containment solution method according to claim 1, wherein the simultaneous steam, noncondensable gas and liquid phase conservation of mass equation and gas and liquid phase conservation of energy equation yields values of principal variables and other variables, including:
and calculating the mass conservation equations of the steam, the non-condensable gas and the liquid phase and the energy conservation equations of the gas phase and the liquid phase, wherein the mass conservation equations of the steam, the non-condensable gas and the liquid phase are in the form of:
wherein, each item of mass conservation equation of steam, non-condensable gas and liquid phase comprises transient term, convection term, mass source term and the like, wherein k is liquid phase, steam or non-condensable gas'kIs a mass source term of k phases;
the gas phase and liquid phase energy conservation equation is in the form of:
wherein each item of the gas phase and liquid phase energy conservation equation comprises a transient term, a convection term, a pressure work-doing term, a heat transfer term and the like, whereinK is a liquid or gas phase, hkIs the specific enthalpy of the k phase, Φ'kIs the heat source term of the k phase, Q'kIs the heat transfer term of the k phase;
and simultaneously establishing a mass conservation equation of steam, noncondensable gas and liquid phase and an energy conservation equation of gas phase and liquid phase, solving based on a Newton-Raphson algorithm to obtain a full-field solution matrix, and calculating according to the full-field solution matrix to obtain main variables such as pressure and enthalpy and values of other variables.
5. The containment solution method according to claim 1, wherein when the flow velocity directions of the gas phase and the liquid phase obtained by two adjacent calculations are judged to be inconsistent, the momentum conservation equation is solved again until the flow velocity directions obtained by two adjacent calculations are consistent.
6. A containment solver, comprising:
the initialization module is used for reading in relevant information of the model;
the solving module is used for solving the momentum conservation equation of the gas phase and the liquid phase according to the related information to obtain the flow velocity of the gas phase and the liquid phase, judging whether the flow velocity directions of the gas phase and the liquid phase obtained by two adjacent calculations are consistent, if so, combining the mass conservation equation of the steam, the non-condensable gas and the liquid phase and the energy conservation equation of the gas phase and the liquid phase to obtain the values of a principal variable and other variables, and correcting the values of the variables such as the flow velocity;
and the output module is used for outputting the solving result.
7. The containment solver of claim 6 wherein said related information includes model type, model parameter values, initial conditions and boundary conditions.
8. The containment solver of claim 6 wherein said solving module is configured to:
calculating terms of a gas phase and liquid phase momentum equation, wherein the gas phase and liquid phase momentum equation is of the form:
wherein each item of the gas phase and liquid phase momentum equation comprises a transient term, a convection term, a pressure term, a gravity term and a resistance term, wherein k is the liquid phase or the gas phase, αkIs the void fraction of the k phase, pkIs the density of the k phase, ukIs the velocity of the k phase, t is time, p is pressure, g is acceleration of gravity, Fk' is the resistance experienced by the k-phase;
and solving the momentum equations of the gas phase and the liquid phase to obtain the flow rates of the gas phase and the liquid phase.
9. The containment solver of claim 6 wherein said solving module is configured to:
and calculating the mass conservation equations of the steam, the non-condensable gas and the liquid phase and the energy conservation equations of the gas phase and the liquid phase, wherein the mass conservation equations of the steam, the non-condensable gas and the liquid phase are in the form of:
wherein, each item of mass conservation equation of steam, non-condensable gas and liquid phase comprises transient term, convection term, mass source term and the like, wherein k is liquid phase, steam or non-condensable gas'kIs a mass source term of k phases;
the gas phase and liquid phase energy conservation equation is in the form of:
wherein each item of the gas phase and liquid phase energy conservation equation comprises a transient term, a convection term, a pressure work-doing term, a heat transfer term and the like, wherein k is a liquid phase or a gas phase, hkIs the specific enthalpy of the k phase, Φ'kIs the heat source term of the k phase, Q'kIs the heat transfer term of the k phase;
and simultaneously establishing a mass conservation equation of steam, noncondensable gas and liquid phase and an energy conservation equation of gas phase and liquid phase, solving based on a Newton-Raphson algorithm to obtain a full-field solution matrix, and calculating according to the full-field solution matrix to obtain main variables such as pressure and enthalpy and values of other variables.
10. The containment solver of claim 6 wherein the solving module is configured to re-solve the momentum conservation equation when the directions of the gas-phase flow velocity and the liquid-phase flow velocity obtained by two adjacent calculations are determined to be inconsistent until the directions of the flow velocities obtained by two adjacent calculations are consistent.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910160710.3A CN111651851B (en) | 2019-03-04 | 2019-03-04 | Containment solving method and containment solver |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910160710.3A CN111651851B (en) | 2019-03-04 | 2019-03-04 | Containment solving method and containment solver |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111651851A true CN111651851A (en) | 2020-09-11 |
| CN111651851B CN111651851B (en) | 2024-01-30 |
Family
ID=72346128
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910160710.3A Active CN111651851B (en) | 2019-03-04 | 2019-03-04 | Containment solving method and containment solver |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111651851B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112613240A (en) * | 2020-11-26 | 2021-04-06 | 中国核电工程有限公司 | Calculation method for flow analysis in containment vessel under serious accident |
| CN112613158A (en) * | 2020-11-26 | 2021-04-06 | 中国核电工程有限公司 | Comprehensive analysis method for thermodynamic response of control body in containment vessel under serious accident |
| CN115906596A (en) * | 2022-11-18 | 2023-04-04 | 上海索辰信息科技股份有限公司 | Wall surface oil film calculation method |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20040128111A1 (en) * | 1998-03-24 | 2004-07-01 | Lang Fred D. | Method for detecting heat exchanger tube failures and their location when using input/loss performance monitoring of a recovery boiler |
| JP2015215772A (en) * | 2014-05-12 | 2015-12-03 | 株式会社東芝 | Heat transfer simulation device and heat transfer simulation method |
| CN106547944A (en) * | 2016-09-27 | 2017-03-29 | 大连理工大学 | A kind of dividing plate optimization method for suppressing compressible two phase flow to rock |
| CN107832573A (en) * | 2017-07-28 | 2018-03-23 | 中冶华天工程技术有限公司 | Predict the numerical computation method of flue gas recirculation sintering matter thermal coupling process |
| CN108536952A (en) * | 2018-04-03 | 2018-09-14 | 东北大学 | The computational methods of biphase gas and liquid flow gas holdup in a kind of determining ladle |
| CN207991832U (en) * | 2017-11-10 | 2018-10-19 | 国核华清（北京）核电技术研发中心有限公司 | The passive core cooling system bulk testing rack of pressurized-water reactor nuclear power plant |
-
2019
- 2019-03-04 CN CN201910160710.3A patent/CN111651851B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20040128111A1 (en) * | 1998-03-24 | 2004-07-01 | Lang Fred D. | Method for detecting heat exchanger tube failures and their location when using input/loss performance monitoring of a recovery boiler |
| JP2015215772A (en) * | 2014-05-12 | 2015-12-03 | 株式会社東芝 | Heat transfer simulation device and heat transfer simulation method |
| CN106547944A (en) * | 2016-09-27 | 2017-03-29 | 大连理工大学 | A kind of dividing plate optimization method for suppressing compressible two phase flow to rock |
| CN107832573A (en) * | 2017-07-28 | 2018-03-23 | 中冶华天工程技术有限公司 | Predict the numerical computation method of flue gas recirculation sintering matter thermal coupling process |
| CN207991832U (en) * | 2017-11-10 | 2018-10-19 | 国核华清（北京）核电技术研发中心有限公司 | The passive core cooling system bulk testing rack of pressurized-water reactor nuclear power plant |
| CN108536952A (en) * | 2018-04-03 | 2018-09-14 | 东北大学 | The computational methods of biphase gas and liquid flow gas holdup in a kind of determining ladle |
Non-Patent Citations (4)
| Title |
|---|
| WENQI ZHONG 等: "CFD simulation of dense particulate reaction system：Approaches, recent advances and applications", 《CHEMICAL ENGINEERING SCIENCE》, pages 16 - 43 * |
| 吕彦力 等: "旁路流对换热器壳侧气液两相流动特性的影响", 河南科技大学学报(自然科学版), pages 22 - 25 * |
| 张欣雨 等: "起伏管道空管投油气液两相流动的数值模拟", 中国科学院大学学报 * |
| 李华琪;江新标;陈立新;杨宁;胡攀;马腾跃;张良;: "空间堆堆芯热管蒸气流动计算方法研究", 核动力工程, no. 06, pages 1 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112613240A (en) * | 2020-11-26 | 2021-04-06 | 中国核电工程有限公司 | Calculation method for flow analysis in containment vessel under serious accident |
| CN112613158A (en) * | 2020-11-26 | 2021-04-06 | 中国核电工程有限公司 | Comprehensive analysis method for thermodynamic response of control body in containment vessel under serious accident |
| CN112613158B (en) * | 2020-11-26 | 2024-02-23 | 中国核电工程有限公司 | Comprehensive analysis method for controlling thermodynamic response in containment under serious accident |
| CN115906596A (en) * | 2022-11-18 | 2023-04-04 | 上海索辰信息科技股份有限公司 | Wall surface oil film calculation method |
| CN115906596B (en) * | 2022-11-18 | 2024-03-22 | 上海索辰信息科技股份有限公司 | Wall oil film calculation method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111651851B (en) | 2024-01-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111651851A (en) | Containment solving method and containment solver | |
| Mazzi et al. | Off-design dynamic model of a real Organic Rankine Cycle system fuelled by exhaust gases from industrial processes | |
| CN111723450A (en) | Safety analysis method and system for nuclear power plant | |
| CN109918787A (en) | The analogy method of aqueous vapor two-phase homogeneous flow in aqueduct based on finite volume method | |
| CN104834773A (en) | Simulation method for heat exchange performance of straight tube type once-through steam generator | |
| CN105512433A (en) | Fluid-solid node two-phase flow modeling method | |
| CN108763670A (en) | A kind of solution supercritical carbon dioxide reactor Brayton cycle transient process method | |
| Jordaan et al. | Numerical development of a coupled one-dimensional/three-dimensional computational fluid dynamics method for thermal analysis with flow maldistribution | |
| Zaversky et al. | Transient behavior of an active indirect two-tank thermal energy storage system during changes in operating mode–An application of an experimentally validated numerical model | |
| CN111125972B (en) | Hydraulic load analysis method for water loss accident of break of nuclear power plant | |
| Zhang et al. | Discrete unified gas-kinetic scheme for the conservative Allen-Cahn equation | |
| Dechamps | Modelling the transient behaviour of heat recovery steam generators | |
| CN113283189B (en) | Cross-dimension coupling analysis method for containment and passive containment air cooling system | |
| CN114357380A (en) | Thermal deviation calculation method based on real-time measurement data | |
| CN114388162A (en) | Helium-xenon cooling reactor control method and device and electronic equipment | |
| CN111723451B (en) | Sub-channel solving method for Newton-Rafferson algorithm | |
| Galindo-García et al. | CFD simulations of heat recovery steam generators including tube banks | |
| Ashaaria et al. | Fuzzy state space model for a pressurizer in a nuclear power plant | |
| Fu et al. | Development of Conservative Form of RELAP5 Thermal Hydraulic Equations: Part II—Numerical Approach and Code Results | |
| Altés Buch | Dynamic modeling of a steam Rankine Cycle for concentrated solar power applications | |
| Reznicek | Design and simulation of supercritical carbon dioxide recompression Brayton cycles with regenerators for recuperation | |
| Choi et al. | Integral behavior of the ATLAS facility for a 3-inch small break loss of coolant accident | |
| González-Mora et al. | Alternative Approach For Thermo-Hydraulic Modeling Of Direct Steam Generation In Parabolic Trough Solar Collectors | |
| Eoh et al. | Numerical simulation of a potential CO2 ingress accident in a SFR employing an advanced energy conversion system | |
| Ludington | Tools for supercritical carbon dioxide cycle analysis and the cycle's applicability to sodium fast reactors |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |