CN108763670A - A kind of solution supercritical carbon dioxide reactor Brayton cycle transient process method - Google Patents
A kind of solution supercritical carbon dioxide reactor Brayton cycle transient process method Download PDFInfo
- Publication number
- CN108763670A CN108763670A CN201810463358.6A CN201810463358A CN108763670A CN 108763670 A CN108763670 A CN 108763670A CN 201810463358 A CN201810463358 A CN 201810463358A CN 108763670 A CN108763670 A CN 108763670A
- Authority
- CN
- China
- Prior art keywords
- reactor
- carbon dioxide
- supercritical carbon
- fluid
- node
- 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
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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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/06—Power analysis or power 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/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
- Engine Equipment That Uses Special Cycles (AREA)
Abstract
The invention discloses a kind of solution supercritical carbon dioxide reactor Brayton cycle transient process methods, initialize supercritical carbon dioxide reactor parameter first；Then the quality, interior energy of control volume and the flow taken in current iteration step are calculated and obtains the transient state core power of reactor；Then it using the state of circulating coolant as the boundary condition of reactor fuel rod Calculation of Heat Transfer, calculates Temperature Distribution in fuel rod and updates pressure, flow, the enthalpy parameter taken over；It calculates the change of internal energy rate of each fluid node and is judged, finally judge whether that reaching final calculates the time, to form the transient process of supercritical carbon dioxide Brayton cycle.This invention simplifies the solution flows of supercritical carbon dioxide reactor Brayton cycle transient process, in turn ensure computational accuracy.
Description
Technical field
It is analyzed the invention belongs to nuclear reactor safety and calculates technical field, and in particular to a kind of solution supercritical carbon dioxide
Reactor Brayton cycle transient process method.
Background technology
Supercritical carbon dioxide reactor has many advantages, such as that the thermal efficiency is high, system simplification is compact, is made using Brayton cycle
For energy conversion system, whole system is made of the part such as reactor, regenerator, cooler, compressor, turbine, is to have very much
The forth generation nuclear reactor concept of foreground.Supercritical carbon dioxide reactor uses direct cooling cycle, eliminates steam generation
The equipment such as device, large-scale heat exchanger, simple system.Mechanical equipment in circulation loop is smaller, the startup of capital equipment and
The motor-driven response time substantially reduces, and is conducive to flexibly the operating status of adjustment reactor cycle in time.Based on above-mentioned advantage, surpass
Critical carbon dioxide reactor concept can be not only used for the high efficiency commercial reactor of exploitation different capacity level, can also become core
Submarine, core aircraft carrier etc. require that equipment volume is small, the powered equipment of changed power timely military equipment, have wide city
Field application prospect.
In the Brayton cycle of supercritical carbon dioxide reactor, the supercritical carbon dioxide at suction port of compressor, which is in, to be faced
Near boundary's point.Supercritical carbon dioxide is acutely and rapid in the variation of Near The Critical Point physical property, leads to the cooling in Brayton cycle
Agent physical property will undergo strong nonlinear change, and the characteristic of this coolant physical property acute variation had both improved overcritical titanium dioxide
The cycle efficieny of carbon reactor Brayton cycle, also the Transient Numerical Simulation entirely to recycle bring challenge.Physical property it is violent
Variation can cause to solve can encounter numerical computations and be difficult to restrain in supercritical carbon dioxide reactor Brayton cycle transient process
The problem of, increase the difficulty of transient analysis software development.
However through investigation, rarely have open achievement and patent introduction anti-for solving supercritical carbon dioxide both at home and abroad at present
Answer the numerical method of heap Brayton cycle transient process.
Invention content
In view of the above-mentioned deficiencies in the prior art, the technical problem to be solved by the present invention is that providing, a kind of solution is super to face
Boundary's carbon dioxide reaction heap Brayton cycle transient process method carries out quasi-static place to the turbine machinery in Brayton cycle
Reason, simplifies solution procedure, improves solving speed.
The present invention uses following technical scheme：
A kind of solution supercritical carbon dioxide reactor Brayton cycle transient process method, initializes overcritical two first
Carbon dioxide reaction heap parameter；Then the transient state core power of reactor is calculated；Then using the state of circulating coolant as
The boundary condition of reactor fuel rod Calculation of Heat Transfer calculates Temperature Distribution in fuel rod；Update the pressure taken over, flow, enthalpy
Parameter；It calculates the change of internal energy rate △ U/U of each fluid node and is compared with specified value, finally judge whether to reach final meter
Evaluation time repeats the reactor that above step obtains supercritical carbon dioxide reactor Brayton cycle in each time step
The flow regime of power rating, fuel rod state of temperature, turbine Mechanical Running Condition and cyclic internal cooling but agent is super to be formed
The transient process of critical carbon dioxide Brayton cycle
Specifically, steps are as follows：
S1, initialize the core power of supercritical carbon dioxide reactor, Brayton cycle coolant channel fluid state,
Calculate time step parameter；
The quality of supercritical carbon dioxide and energy conservation equation and the momentum taken over are kept in S2, discrete Brayton cycle
Permanent equation；Quality, interior energy and the flow taken over of control volume in current iteration step is calculated；
S3, by the enthalpy and specific volume of each node in recycling, according to instantaneous in CO 2 fluid physical property packet calculate node
Pressure；
S4, by calculate Point reactor nutron kinetics equation, the transient state core power of reactor is calculated；
Boundary condition of the state as reactor fuel rod Calculation of Heat Transfer of S5, circulating coolant, are changed by heat conduction, convection current
Heat calculates Temperature Distribution each in fuel rod；
The parameters such as pressure, flow, the enthalpy that S6, update are taken over；
S7, the change of internal energy rate △ U/U for calculating each fluid node, if the △ U/U maximum values of each node are more than specified value,
Then think that current iteration step calculates not restrain, time step halves, and returns to step S3 and recalculates；If the △ U/ of each fluid node
U maximum values are less than specified value, then it is assumed that the calculating convergence of current iteration step, into next step；
S8, judge whether to reach the final calculating time, if reached, terminate calculating；If do not reached, redirect
To step S2, the calculating of next time step is carried out.
Further, step S2 is specially：The quality and the conservation of energy of supercritical carbon dioxide in discrete Brayton cycle
Equation and the momentum conservation equation taken over, using Gauss-Saden that iteration to being recycled about the solution of Jacobin matrix J
Mass incremental, energy increment and the flow increment Δ Y respectively taken over of each node of fluid, to which next iteration step be calculated
Node quality, energy and the flow taken over, the calculating of mass incremental, energy increment and flow increment Δ Y is as follows：
Δ Y=J-1·B
Wherein, Δ t is time step, and nv is number of nodes, and nj is take over quantity, YnFor the argument value at current time, n
Indicate current time.
Further, the energy conservation equation of CO 2 fluid and mass-conservation equation are applied on control volume, meter
The variation of the quality, interior energy and pressure of control volume is calculated, momentum conservation equation is applied to take over, and calculating flows in and out control
The coolant flow of body, the conservation of mass of all control volumes, the difference equation of the conservation of energy and the conservation of momentum taken over is poor
Equations simultaneousness is divided to obtain following formula together：
Above formula is substituted into after the right first item is carried out first order Taylor series expansion and the time-derivative in left side is carried out discrete
Obtain the mass incremental and energy increment and take over flow increment Δ Y of node.
Further, the quality of CO 2 fluid, momentum and energy conservation equation are respectively：
Wherein, ρ is density；W is fluid flow；V is fluid flow rate；P is pressure；DhFor equivalent diameter；F is frictional resistance
Coefficient；gzFor acceleration of gravity；U is fluid interior energy；H is fluid enthalpy；Q is fluid endogenous pyrogen, and j is the section for flowing into present node
Point quantity, z is height of node.
Further, it when fluid flows through gas-turbine tool, uses and the node pressure that turbine machinery is brought is calculated with drag
Variation and energy variation：
Δ P=P1 T(Rp-1)
Wherein, Δ P is turbine machinery pressure change, P1 TFor turbine mechanical inputs stagnation pressure, RpFor turbine machinery pressure ratio, by property
Energy curve is inquired to obtain, Wc.v.For turbine mechanical convection body amount of work,For turbine mechanical flow, ηadFor the thermal insulation of turbine machinery
Efficiency is inquired to obtain by performance curve,Enthalpy is exported for turbine machinery,Enthalpy is exported for turbine machinery；
Turbine machinery is the node pressure variation and energy variation that fluid is brought, and will be respectively acting on fluid fundamental equation
Momentum conservation equation and energy conservation equation.
Further, in step S2, the pipeline in Brayton cycle is indicated with node, spatially fluid recirculation system
Discrete using staggered-mesh, momentum grid cell deviates quality and energy grid unit half space increment.
Further, in step S4, the reactor Point reactor nutron kinetics equation of six groups of delayed neutrons is as follows：
Wherein：N (t) is reactor core fission power；T is the time；ρ (t) is global reactivity；β is total delayed neutron fraction；Λ
It is prompt neutron for the time；λiFor the decay coefficient of i-th group of delayed neutron；Ci(t) it is the fission power of i-th group of delayed neutron.
Further, in step S5, the Temperature Distribution in fuel rod is as follows：
Wherein：ρ is density；cpFor specific heat；T is temperature；T is the time；K is thermal coefficient；Q is endogenous pyrogen power density.
Further, the efficiency and pressure ratio for solving turbine machinery come interpolation using the performance curve of turbine machinery, in transient state
When analyzing program solution, turbine machinery is one and takes over a node that the pressure liter being calculated by pressure ratio or pressure drop will act on
In the work(that the rotation axis on the equation of momentum of turbine machinery entrance sleeve, being calculated does turbine machinery internal flow, with source
The form of item increases the energy conservation equation in turbine mechanical fluid node.
Compared with prior art, the present invention at least has the advantages that：
A kind of solution supercritical carbon dioxide reactor Brayton cycle transient process method of the present invention, with fluid homogeneous flow
Based on model and reactor point pile neutron equation, in conjunction with the physical property of supercritical carbon dioxide, friction and heat transfer characteristic, face super
Boundary CO 2 fluid fundamental equation carries out time and spatially discrete, uses transformation period step-length or for Jacobian matrix
Increase the modes such as multiplication factor to promote convergence of this method near carbon dioxide critical point, efficiently solves Bretton and follow
The numerical divergence problem and turbine machinery transient process Solve problems that ring transient state can encounter when solving.
Further, state of the supercritical carbon dioxide in Brayton cycle is essentially single-phase, using homogeneous phase model
The fundamental equation of fluid is calculated, the number of fundamental equation can be effectively reduced, is reduced required alternate when fundamental equation calculates
Constitutive relationship not only simplifies calculating step, and ensures that this method can accurately obtain the state change of supercritical carbon dioxide
Change.
Further, it uses and shortens the time step calculated or increase the modes such as multiplication factor for Jacobian matrix, it can
Effectively to solve the problems, such as that the numerical solution that physical property mutation is brought when supercritical carbon dioxide intends Near The Critical Point is not restrained.
Further, fuel Doppler effect, cooling succinctly can effectively be solved using Point reactor nutron kinetics equation
The reactor capability feedback that agent temperature effect and control rod action are brought.
Further, one-dimension conduct equation can accurately calculate the distribution of temperature radially in fuel rod, be followed for Bretton
Ring transient analysis provides the fuel pellet and wall surface temperature result of variations that researcher is concerned about.
Further, the turbine machinery performance curve being input from the outside provides the pressure ratio and efficiency of turbine machinery, can be fast
The transient operating condition of fast accurate simulation gas-turbine tool, meanwhile, this method can unify the substantially square of fluid in turbine machinery
Journey form and the fluid fundamental equation form in common pipe are the key that realize that Closed Brayton Power Cycle transient process calculates.
Further, accurate supercritical carbon dioxide physical property, heat exchange and friction law formula, ensure that this method just
True property.
In conclusion the present invention is based on simplified fluid fundamental equation, Point reactor nutron kinetics equation, rational gas-turbines
Tool model, accurate supercritical carbon dioxide physical property and heat exchange friction law formula and effective method of value solving, not only simplify
The solution flow of supercritical carbon dioxide reactor Brayton cycle transient process in turn ensures computational accuracy, ensures to be based on being somebody's turn to do
The transient process of method solves the transient physical process that can accurately react supercritical carbon dioxide reactor Brayton cycle.
Below by drawings and examples, technical scheme of the present invention will be described in further detail.
Description of the drawings
Fig. 1 is flow chart of the present invention；
Fig. 2 is the discrete figure of Brayton cycle coolant channel of the present invention；
Fig. 3 is turbine machinery model schematic.
Specific implementation mode
The present invention provides a kind of solution supercritical carbon dioxide reactor Brayton cycle transient process methods, first just
The core power of beginningization supercritical carbon dioxide reactor, calculates time step at Brayton cycle coolant channel fluid state
Parameter；Then quality, interior energy and the flow taken over of control volume in current iteration step is calculated；Pass through the wink in node
When pressure；Obtain the transient state core power of reactor；Then using the state of circulating coolant as reactor fuel rod in terms of heat transfer
The boundary condition of calculation calculates Temperature Distribution in fuel rod；Update pressure, flow, the enthalpy parameter taken over；Calculate each fluid node
Change of internal energy rate △ U/U, if the change of internal energy rate maximum value of each node be more than specified value, then it is assumed that current iteration step calculate
It does not restrain, time step halves, and recalculates the instantaneous pressure in node；If the maximum value of the change of internal energy rate of each fluid node
Less than specified value, then it is assumed that the calculating convergence of current iteration step finally judges whether to reach the final calculating time.More than being based on
Method for solving can obtain the reactor capability shape of supercritical carbon dioxide reactor Brayton cycle in each time step
The flow regime of state, fuel rod state of temperature, turbine Mechanical Running Condition and cyclic internal cooling but agent, to form overcritical two
The transient process of carbonoxide Brayton cycle.
The present invention can accurately and effectively solve the transient process that supercritical carbon dioxide reactor Bretton follows, and can be used for
The crash analysis of supercritical carbon dioxide reactor Brayton cycle and control strategy research.
Referring to Fig. 1, a kind of solution supercritical carbon dioxide reactor Brayton cycle transient process method of the present invention, packet
Include following steps：
S1, initialize the core power of supercritical carbon dioxide reactor, Brayton cycle coolant channel fluid state,
Calculate the parameters such as time step；
The quality of supercritical carbon dioxide and energy conservation equation and the momentum taken over are kept in S2, discrete Brayton cycle
Permanent equation；Quality, interior energy and the flow taken over of control volume in current iteration step is calculated；
Quality, momentum and the energy conservation equation of CO 2 fluid be respectively：
Wherein：ρ is density, kg/m3；A is area, m2；W is fluid flow, kg/s；V is fluid flow rate, m/s；P is pressure
Power, MPa；DhFor equivalent diameter, m；F is coefficient of frictional resistance；gzFor acceleration of gravity, m/s2；U is fluid interior energy, kJ/kg；h
For fluid enthalpy, kJ/kg；Q is fluid endogenous pyrogen, kW/kg；J is the number of nodes for flowing into present node, and z is height of node；
When fluid flows through gas-turbine tool (gas-turbine and compressor), it will use and calculate what steamer machinery was brought with drag
Node pressure changes and energy variation：
Δ P=P1 T(Rp-1)
Wherein, Δ P is turbine machinery pressure change, P1 TFor turbine mechanical inputs stagnation pressure, RpFor turbine machinery pressure ratio, by property
Energy curve is inquired to obtain, Wc.v.For turbine mechanical convection body amount of work,For turbine mechanical flow, ηadFor the thermal insulation of turbine machinery
Efficiency is inquired to obtain by performance curve,Enthalpy is exported for turbine machinery,Enthalpy is exported for turbine machinery.
Turbine machinery is the node pressure variation and energy variation that fluid is brought, and will be respectively acting on fluid fundamental equation
Momentum conservation equation and energy conservation equation.
In numerical computations, the pipeline in Brayton cycle is all indicated with node.Spatially fluid recirculation system
Discrete using staggered-mesh, momentum grid cell deviates quality and energy grid unit half space increment.Fig. 2 gives two
Intersection mesh generation between adjacent node K and node L.It is take over j between node K and node L.Energy conservation equation and matter
It measures conservation equation to be applied on node, the variation of the quality, interior energy and pressure of calculate node, momentum conservation equation is applied to connect
Guan Shang calculates the coolant flow for flowing in and out control volume, Wj-1, Wj and Wj+1 as shown.
By the conservation of mass of all nodes, the difference equation of the conservation of energy and the conservation of momentum difference equation simultaneous taken over
Together, following form can be obtained：
First order Taylor series expansion is carried out to first item on the right of formula (4), can be obtained：
f(Yn+1)=f (Yn)+J(Yn)·(Yn+1-Yn) (5)
Formula (5) substitutes into formula (4) and can obtain：
It is discrete to the time-derivative progress on the left of formula (6), it can obtain：
Formula (7), which is substituted into formula (6), can finally obtain：
J Δs Y=B (8)
Wherein,
Δ Y=[Δ W1... Δ Wj... Δ Wnj,ΔM1... Δ MK... Δ Mnv,ΔU1... Δ UK... Δ Unv]T
Wherein, Δ t is time step, and nv is number of nodes, and nj is take over quantity, YnFor the argument value at current time, n
Indicate current time.
Using Gauss-Saden your iteration can solution formula (8), obtain the mass incremental of each node of circulation of fluid, energy increases
Amount and the flow increment respectively taken over, to be calculated next iteration step node quality, energy and the flow taken over.
S3, by the enthalpy and specific volume of each node in recycling, according to instantaneous in CO 2 fluid physical property packet calculate node
Pressure；
S4, by calculate Point reactor nutron kinetics equation, the transient state core power of reactor is calculated；
Consider that the reactor Point reactor nutron kinetics equation of six groups of delayed neutrons is：
In formula：N (t) is reactor core fission power, W；T is time, s；ρ (t) is global reactivity, $；β is total delayed neutron
Share；Λ is prompt neutron for time, s；λiFor the decay coefficient of i-th group of delayed neutron, s-1；Ci(t) it is i-th group of delayed neutron
Fission power, W.
Boundary condition of the state as reactor fuel rod Calculation of Heat Transfer of S5, circulating coolant, are changed by heat conduction, convection current
Heat calculates Temperature Distribution each in fuel rod；
In formula：ρ is density；cpFor specific heat；T is temperature；T is the time；K is thermal coefficient；Q is endogenous pyrogen power density.
The parameters such as pressure, flow, the enthalpy that S6, update are taken over；
S7, the change of internal energy rate △ U/U for calculating each fluid node, if the maximum value of the change of internal energy rate of each node is more than
Specified value, then it is assumed that current iteration step, which calculates, not to be restrained, and time step halves, and returns to step S3 and recalculates；If each fluid section
The maximum value of the change of internal energy rate of point is less than specified value, then it is assumed that the calculating convergence of current iteration step, into next step；
S8, judge whether to reach the final calculating time, if reached, terminate calculating；If do not reached, redirect
To step S2, the calculating of next time step is carried out.
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is
A part of the embodiment of the present invention, instead of all the embodiments.The present invention being described and shown in usually here in attached drawing is real
Applying the component of example can be arranged and be designed by a variety of different configurations.Therefore, the present invention to providing in the accompanying drawings below
The detailed description of embodiment be not intended to limit the range of claimed invention, but be merely representative of the selected of the present invention
Embodiment.Based on the embodiments of the present invention, those of ordinary skill in the art are obtained without creative efforts
The every other embodiment obtained, shall fall within the protection scope of the present invention.
Supercritical carbon dioxide heat exchange, which calculates, uses disclosed Gnielinski relational expressions, disclosed in the use of friction law formula
Zigrang-Sylvester relational expressions.
Referring to Fig. 3, for the solution of the turbine mechanical (compressor and gas-turbine) in Brayton cycle, it is believed that transient state mistake
The change of fluid of turbine machine internal is in quasi-static state in journey.Under the assumptions, using the performance curve of turbine machinery come
Interpolation solves the efficiency and pressure ratio of turbine mechanical (compressor and gas-turbine).In transient analysis program solution, by turbine machinery
One is regarded as to take over and a node.The pressure liter (or pressure drop) being calculated by pressure ratio will act on turbine machinery entrance sleeve
The equation of momentum on, and the work(that the rotation axis being calculated does turbine machinery internal flow, it will increased in the form of source item
In the energy conservation equation of turbine mechanical fluid node.
The above content is merely illustrative of the invention's technical idea, and protection scope of the present invention cannot be limited with this, every to press
According to technological thought proposed by the present invention, any change done on the basis of technical solution each falls within claims of the present invention
Protection domain within.
Claims (10)
1. a kind of solution supercritical carbon dioxide reactor Brayton cycle transient process method, which is characterized in that initial first
Change supercritical carbon dioxide reactor parameter；Then the transient state core power of reactor is calculated；Then with circulating coolant
Boundary condition of the state as reactor fuel rod Calculation of Heat Transfer, calculate Temperature Distribution in fuel rod；Update take over pressure,
Flow, enthalpy parameter；It calculates the change of internal energy rate △ U/U of each fluid node and is compared with specified value, finally judged whether
Reach final and calculate the time, repetition above step obtains supercritical carbon dioxide reactor Bretton in each time step and follows
The flow regime of reactor capability state, fuel rod state of temperature, turbine Mechanical Running Condition and the cyclic internal cooling of ring but agent,
To form the transient process of supercritical carbon dioxide Brayton cycle.
2. a kind of solution supercritical carbon dioxide reactor Brayton cycle transient process method according to claim 1,
It is characterized in that, being as follows：
S1, it the core power for initializing supercritical carbon dioxide reactor, Brayton cycle coolant channel fluid state, calculates
Time step parameter；
The quality of supercritical carbon dioxide and energy conservation equation and the conservation of momentum side taken in S2, discrete Brayton cycle
Journey；Quality, interior energy and the flow taken over of control volume in current iteration step is calculated；
S3, by the enthalpy and specific volume of each node in recycling, according to the instantaneous pressure in CO 2 fluid physical property packet calculate node；
S4, by calculate Point reactor nutron kinetics equation, the transient state core power of reactor is calculated；
Boundary condition of the state as reactor fuel rod Calculation of Heat Transfer of S5, circulating coolant, pass through heat conduction, heat convection meter
Calculate Temperature Distribution each in fuel rod；
The parameters such as pressure, flow, the enthalpy that S6, update are taken over；
S7, the change of internal energy rate △ U/U for calculating each fluid node recognize if the △ U/U maximum values of each node are more than specified value
It calculates and does not restrain for current iteration step, time step halves, and returns to step S3 and recalculates；If the △ U/U of each fluid node are most
Big value is less than specified value, then it is assumed that the calculating convergence of current iteration step, into next step；
S8, judge whether to reach the final calculating time, if reached, terminate calculating；If do not reached, step is jumped to
Rapid S2 carries out the calculating of next time step.
3. a kind of solution supercritical carbon dioxide reactor Brayton cycle transient process method according to claim 2,
It is characterized in that, step S2 is specially：In discrete Brayton cycle the quality of supercritical carbon dioxide and energy conservation equation with
And the momentum conservation equation taken over, using Gauss-Saden, your iteration is each to obtaining circulation of fluid about the solution of Jacobin matrix J
Mass incremental, energy increment and the flow increment Δ Y respectively taken over of node, to which the node of next iteration step be calculated
Quality, energy and the flow taken over, the calculating of mass incremental, energy increment and flow increment Δ Y is as follows：
Δ Y=J-1·B
Wherein, Δ t is time step, and nv is number of nodes, and nj is take over quantity, YnFor the argument value at current time, n is indicated
Current time.
4. a kind of solution supercritical carbon dioxide reactor Brayton cycle transient process method according to claim 3,
It is characterized in that, the energy conservation equation of CO 2 fluid and mass-conservation equation are applied on control volume, control is calculated
The variation of the quality, interior energy and pressure of body, momentum conservation equation are applied to take over, and calculating flows in and out the cold of control volume
But agent flux, by the conservation of mass of all control volumes, the difference equation of the conservation of energy and the conservation of momentum difference equation taken over
Simultaneous obtains following formula together：
Above formula is substituted into after the right first item is carried out first order Taylor series expansion and discrete obtain is carried out to the time-derivative in left side
The mass incremental and energy increment and take over flow increment Δ Y of node.
5. a kind of solution supercritical carbon dioxide reactor Brayton cycle transient process method according to claim 4,
It is characterized in that, the quality of CO 2 fluid, momentum and energy conservation equation are respectively：
Wherein, ρ is density；W is fluid flow；V is fluid flow rate；P is pressure；DhFor equivalent diameter；F is coefficient of frictional resistance；
gzFor acceleration of gravity；U is fluid interior energy；H is fluid enthalpy；Q is fluid endogenous pyrogen, and j is the number of nodes for flowing into present node
Amount, z is height of node.
6. a kind of solution supercritical carbon dioxide reactor Brayton cycle transient process method according to claim 5,
It is characterized in that, when fluid flows through gas-turbine tool, use with drag calculate node pressure variation that turbine machinery is brought and
Energy variation：
Δ P=P1 T(Rp-1)
Wherein, Δ P is turbine machinery pressure change, P1 TFor turbine mechanical inputs stagnation pressure, RpIt is bent by performance for turbine machinery pressure ratio
Line is inquired to obtain, Wc.v.For turbine mechanical convection body amount of work,For turbine mechanical flow, ηadFor the adiabatic effect of turbine machinery
Rate is inquired to obtain by performance curve,Enthalpy is exported for turbine machinery,Enthalpy is exported for turbine machinery；
Turbine machinery is the node pressure variation and energy variation that fluid is brought, and will be respectively acting on the momentum of fluid fundamental equation
Conservation equation and energy conservation equation.
7. a kind of solution supercritical carbon dioxide reactor Brayton cycle transient process side according to claim 2 or 5
Method, which is characterized in that in step S2, the pipeline in Brayton cycle is indicated with node, and spatially fluid recirculation system is adopted
Discrete with staggered-mesh, momentum grid cell deviates quality and energy grid unit half space increment.
8. a kind of solution supercritical carbon dioxide reactor Brayton cycle transient process method according to claim 2,
It is characterized in that, in step S4, the reactor Point reactor nutron kinetics equation of six groups of delayed neutrons is as follows：
Wherein：N (t) is reactor core fission power；T is the time；ρ (t) is global reactivity；β is total delayed neutron fraction；Λ is wink
Send out neutron generation time；λiFor the decay coefficient of i-th group of delayed neutron；Ci(t) it is the fission power of i-th group of delayed neutron.
9. a kind of solution supercritical carbon dioxide reactor Brayton cycle transient process method according to claim 2,
It is characterized in that, in step S5, the Temperature Distribution in fuel rod is as follows：
Wherein：ρ is density；cpFor specific heat；T is temperature；T is the time；K is thermal coefficient；Q is endogenous pyrogen power density.
10. a kind of solution supercritical carbon dioxide reactor Brayton cycle transient process side according to claim 1 or 2
Method, which is characterized in that the efficiency and pressure ratio for solving turbine machinery come interpolation using the performance curve of turbine machinery, in transient analysis
When program solution, turbine machinery is taken over for one and will act on gas with a node, the pressure liter being calculated by pressure ratio or pressure drop
On the equation of momentum of wheel mechanical entrance sleeve, the work(that the rotation axis that is calculated does turbine machinery internal flow, with source item
Form increases the energy conservation equation in turbine mechanical fluid node.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810463358.6A CN108763670B (en) | 2018-05-15 | 2018-05-15 | Method for solving Brayton cycle transient process of supercritical carbon dioxide reactor |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810463358.6A CN108763670B (en) | 2018-05-15 | 2018-05-15 | Method for solving Brayton cycle transient process of supercritical carbon dioxide reactor |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108763670A true CN108763670A (en) | 2018-11-06 |
| CN108763670B CN108763670B (en) | 2021-04-20 |
Family
ID=64007763
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810463358.6A Active CN108763670B (en) | 2018-05-15 | 2018-05-15 | Method for solving Brayton cycle transient process of supercritical carbon dioxide reactor |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108763670B (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109783904A (en) * | 2018-12-28 | 2019-05-21 | 西安交通大学 | A kind of width parameter area carbon dioxide physical property method for solving |
| CN111680405A (en) * | 2020-05-25 | 2020-09-18 | 西安理工大学 | Natural circulation capacity hydraulic characteristic calculation method |
| CN111723451A (en) * | 2019-03-04 | 2020-09-29 | 国家电投集团科学技术研究院有限公司 | Sub-channel solving method of Newton-Rafferson algorithm |
| CN112071457A (en) * | 2020-08-07 | 2020-12-11 | 西安交通大学 | Load tracking method for supercritical carbon dioxide direct cooling reactor system |
| CN112733469A (en) * | 2020-12-30 | 2021-04-30 | 西安交通大学 | Supercritical fluid applicability evaluation method of Brayton cycle system of special reactor |
| CN113255268A (en) * | 2021-05-21 | 2021-08-13 | 北京华大九天科技股份有限公司 | Method for detecting and repairing transient analysis non-convergence in circuit simulation |
| CN115374724A (en) * | 2022-08-29 | 2022-11-22 | 中国核动力研究设计院 | Fluid critical heat flux density analysis method in natural circulation system under motion condition |
| CN117521560A (en) * | 2024-01-03 | 2024-02-06 | 中国核动力研究设计院 | Modeling method and device for supercritical carbon dioxide turbine model and computing equipment |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1775732B1 (en) * | 2005-08-31 | 2013-07-31 | Westinghouse Electric Sweden AB | A method of estimating dryout properties in a nuclear light water reactor |
| CN103902784A (en) * | 2014-04-11 | 2014-07-02 | 华北电力大学 | Safety analysis calculating device for transient nuclear heat coupling of supercritical water reactor |
| CN105653869A (en) * | 2016-01-05 | 2016-06-08 | 中国核动力研究设计院 | Three-dimensional transient performance analysis method for supercritical water reactor core |
| CN106202866A (en) * | 2016-06-24 | 2016-12-07 | 西安交通大学 | One stablizes accurate reactor physics thermal technology's coupling calculation |
-
2018
- 2018-05-15 CN CN201810463358.6A patent/CN108763670B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1775732B1 (en) * | 2005-08-31 | 2013-07-31 | Westinghouse Electric Sweden AB | A method of estimating dryout properties in a nuclear light water reactor |
| CN103902784A (en) * | 2014-04-11 | 2014-07-02 | 华北电力大学 | Safety analysis calculating device for transient nuclear heat coupling of supercritical water reactor |
| CN105653869A (en) * | 2016-01-05 | 2016-06-08 | 中国核动力研究设计院 | Three-dimensional transient performance analysis method for supercritical water reactor core |
| CN106202866A (en) * | 2016-06-24 | 2016-12-07 | 西安交通大学 | One stablizes accurate reactor physics thermal technology's coupling calculation |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109783904A (en) * | 2018-12-28 | 2019-05-21 | 西安交通大学 | A kind of width parameter area carbon dioxide physical property method for solving |
| CN111723451A (en) * | 2019-03-04 | 2020-09-29 | 国家电投集团科学技术研究院有限公司 | Sub-channel solving method of Newton-Rafferson algorithm |
| CN111680405A (en) * | 2020-05-25 | 2020-09-18 | 西安理工大学 | Natural circulation capacity hydraulic characteristic calculation method |
| CN112071457B (en) * | 2020-08-07 | 2023-04-07 | 西安交通大学 | Load tracking method for supercritical carbon dioxide direct cooling reactor system |
| CN112071457A (en) * | 2020-08-07 | 2020-12-11 | 西安交通大学 | Load tracking method for supercritical carbon dioxide direct cooling reactor system |
| CN112733469A (en) * | 2020-12-30 | 2021-04-30 | 西安交通大学 | Supercritical fluid applicability evaluation method of Brayton cycle system of special reactor |
| CN112733469B (en) * | 2020-12-30 | 2022-09-30 | 西安交通大学 | Supercritical fluid applicability evaluation method for Brayton cycle system of special reactor |
| CN113255268A (en) * | 2021-05-21 | 2021-08-13 | 北京华大九天科技股份有限公司 | Method for detecting and repairing transient analysis non-convergence in circuit simulation |
| CN113255268B (en) * | 2021-05-21 | 2022-05-24 | 北京华大九天科技股份有限公司 | Method for detecting and repairing transient analysis non-convergence in circuit simulation |
| CN115374724A (en) * | 2022-08-29 | 2022-11-22 | 中国核动力研究设计院 | Fluid critical heat flux density analysis method in natural circulation system under motion condition |
| CN115374724B (en) * | 2022-08-29 | 2023-08-22 | 中国核动力研究设计院 | Fluid critical heat flow density analysis method in natural circulation system under motion condition |
| CN117521560A (en) * | 2024-01-03 | 2024-02-06 | 中国核动力研究设计院 | Modeling method and device for supercritical carbon dioxide turbine model and computing equipment |
| CN117521560B (en) * | 2024-01-03 | 2024-03-26 | 中国核动力研究设计院 | Modeling method and device for supercritical carbon dioxide turbine model and computing equipment |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108763670B (en) | 2021-04-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108763670A (en) | A kind of solution supercritical carbon dioxide reactor Brayton cycle transient process method | |
| CN111680458B (en) | Thermodynamic hydraulic transient calculation method suitable for sodium water direct-current steam generator | |
| Liu et al. | Core thermal-hydraulic evaluation of a heat pipe cooled nuclear reactor | |
| Li et al. | Dynamic simulation of two-tank indirect thermal energy storage system with molten salt | |
| Ming et al. | Dynamic modeling and validation of the 5 MW small modular supercritical CO2 Brayton-Cycle reactor system | |
| CN103902784A (en) | Safety analysis calculating device for transient nuclear heat coupling of supercritical water reactor | |
| CN109241662A (en) | A kind of electric heating interconnection integrated energy system optimal energy flow calculation methodologies | |
| CN106709133A (en) | Method for soft measurement of nuclear power station reactor core temperature fields on basis of neutral network surface fitting | |
| CN114444413B (en) | Plate-shaped fuel reactor core sub-channel level three-dimensional thermal hydraulic analysis method | |
| Ali et al. | Overall efficiency analysis of an innovative load-following nuclear power plant-thermal energy storage coupled cycle | |
| CN111680417B (en) | Thermal design method suitable for sodium water direct-current steam generator | |
| CN114330169B (en) | Physical and thermal coupling analysis method for small modularized super-safe gas cooled reactor | |
| Zhang et al. | Transfer function modeling and simulation of HPR1000 | |
| Lu et al. | Dynamic modeling and numerical investigation of novel pumped thermal electricity storage system during startup process | |
| Ferroni et al. | Exergy analysis of a PWR core heat transfer | |
| CN113793711A (en) | Method for analyzing coupling heat transfer characteristics of lithium-cooled nuclear reactor and Stirling generator | |
| Prajapati et al. | Power conversion from spherical tokamak test reactor with helium-cooled and water-cooled blanket | |
| Liao et al. | Dynamic Modeling and Characteristic Analysis of Helium Brayton Cycle System | |
| Khan et al. | Performance optimization of power generation system for helium gas/lead lithium dual-coolant blanket system | |
| Deng et al. | Numerical analysis and optimization of the charging process on a shell‐and‐tube latent heat thermal energy storage unit for a solar power plant with direct steam generation | |
| Zhang et al. | Characteristics and optimization of SCO2 Brayton cycle system for high power sodium-cooled fast reactor on Mars | |
| D’OLIVEIRA et al. | Numerical analysis on hot water storage tank integrated with phase change materials | |
| Tang et al. | Calculation and Analysis on Performance of Boiler Thermal Sensor Based on Topological Structure of Fluid Network | |
| Ali et al. | The development of a thermal-hydraulic and nonlinear dynamic system for molten salt reactors | |
| Nejedly et al. | Studies on components for a molten salt reactor |
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 |