A method and system for calculating the core gas-liquid two-phase nuclear-thermal coupling of a lead-based reactor

Abstract

The present application relates to a kind of lead-based reactor core gas-liquid two-phase nuclear heat coupling calculation method and system, the method includes: constructing core geometry model, generating reaction cross section data;Core geometry model includes fuel rod region and coolant region, and reaction cross section data includes multiple groups, multiple temperature points and multiple void fraction;According to reaction cross section data, based on the core neutron calculation model, the core power distribution is calculated;According to the core power distribution, based on the core gas-liquid two-phase thermal hydraulic calculation model, the thermal parameter of current time step is obtained, and thermal parameter includes fuel rod temperature, coolant temperature and void fraction.The calculation method provided by the present application can accurately simulate the influence of gas phase on neutron characteristics, and then accurately predict the core power distribution and thermal parameter under gas-liquid two-phase flow state, to cope with the power transient safety problem caused by gas-liquid two-phase flow that may occur in core under severe accident, to provide technical support for safety analysis and accident evaluation of lead-based reactor.

Description

Technical Field

[0001] This invention relates to the field of nuclear reactor design technology, and more specifically, to a calculation method and system for gas-liquid two-phase nuclear thermal coupling in lead-based reactor cores. Background Technology

[0002] With the development of nuclear energy technology, lead-based reactors, as an important type of fourth-generation advanced nuclear energy system, have attracted widespread attention due to their inherent advantages in safety, economy, and sustainability. Nuclear-thermal coupling calculations of lead-based reactor cores are a key technology for reactor design, safety analysis, and operation optimization, and are of great significance for accurately predicting the neutronics and thermal-hydraulic characteristics of the reactor core.

[0003] In existing technologies, nuclear-thermal coupling calculations in the field of lead-based reactors are limited to nuclear-thermal coupling calculations in the single-phase flow state of liquid metal. Under severe accident conditions in lead-based reactors, the reactor core may exhibit a gas-liquid two-phase flow state. The presence of the gas phase will significantly affect neutronics characteristics, thereby affecting power distribution and thermal-hydraulic characteristics; specifically, two-phase flow can induce transient changes in core power, threatening the reactor's safety boundary.

[0004] Therefore, there is an urgent need to develop a nuclear thermal coupling calculation method for the gas-liquid two-phase flow of lead-based reactors, which can accurately simulate the influence of gas phase on neutronics properties and precisely predict the core power distribution and thermal parameters under the gas-liquid two-phase flow state, providing technical support for the safety analysis and accident assessment of lead-based reactors. Summary of the Invention

[0005] This invention provides a calculation method for gas-liquid two-phase nuclear thermal coupling in lead-based reactor cores, which solves the problem that existing calculation methods for gas-liquid two-phase nuclear thermal coupling in lead-based reactor cores are only applicable to single-phase liquid metal flow and cannot cope with the power transient safety problem caused by gas-liquid two-phase flow that may occur in the reactor core under severe accidents.

[0006] This invention provides a calculation method for the gas-liquid two-phase nuclear thermal coupling of lead-based reactor cores to solve technical problems.

[0007] This invention provides a method for calculating the gas-liquid two-phase nuclear-thermal coupling of a lead-based reactor core, the method comprising:

[0008] S100, construct the core geometry model and generate reaction section data; wherein, the core geometry model includes the fuel rod region and the coolant region, and the reaction section data includes multiple groups, multiple temperature points and multiple cavitation fractions;

[0009] S120, Based on the reaction cross-section data, calculate the core power distribution using a core neutronics calculation model;

[0010] S140, based on the core power distribution and the core gas-liquid two-phase thermal-hydraulic calculation model, the thermal parameters for the current time step are obtained. These thermal parameters include fuel rod temperature, coolant temperature, and void fraction. The core gas-liquid two-phase thermal-hydraulic calculation model includes a thermal-hydraulic model and a transport model. The thermal-hydraulic model is used to calculate the fuel rod temperature and the coolant temperature. The transport model includes gas-liquid two-phase mass conservation equations, gas-liquid two-phase momentum conservation equations, and gas-liquid two-phase energy conservation equations, used to calculate the void fraction.

[0011] The mass conservation equation for a gas-liquid two-phase system is:

[0012] ;

[0013] In the formula, The density of the gas-liquid two-phase mixture is denoted as . For time, The velocity of the gas and liquid phases, For the fraction of voids, The density is the gas phase density. The density of the liquid phase;

[0014] The momentum conservation equation for a gas-liquid two-phase system is:

[0015] ;

[0016] In the formula, The density of the gas-liquid two-phase mixture is denoted as . For time, The velocity of the gas and liquid phases, For flow field pressure, The dynamic viscosity of the gas-liquid two-phase mixture is denoted as . For the fraction of voids, The dynamic viscosity of the gas phase. The dynamic viscosity of the liquid phase. It is the acceleration due to gravity;

[0017] The energy conservation equation for a gas-liquid two-phase system is:

[0018] ;

[0019] In the formula, Let be the density of the gas-liquid two-phase mixture. Enthalpy of a gas-liquid two-phase mixture For time, The velocity of the gas and liquid phases, Enthalpy of a gas-liquid two-phase mixture For the fraction of voids, For vapor phase enthalpy, For liquid phase enthalpy, Thermal conductivity, The characteristic temperature of the coolant. This is the rate at which heat is released from fuel to coolant per unit volume.

[0020] Optionally, step S140, which involves obtaining the thermal parameters for the current time step based on the core power distribution and the core gas-liquid two-phase thermal-hydraulic calculation model, includes:

[0021] S142, Based on the core power distribution and the thermal-hydraulic model, the fuel rod temperature and the coolant temperature are calculated.

[0022] S144, solve the gas-liquid two-phase mass conservation equation, gas-liquid two-phase momentum conservation equation, and gas-liquid two-phase energy conservation equation in sequence to calculate the cavitation fraction.

[0023] Optionally, step S100, the step of constructing the core geometric model and generating reaction section data, includes:

[0024] S102, Construct the core geometry model;

[0025] S104, insert bubbles into the coolant region to obtain a core geometry model containing bubbles; this setting is to simulate the gas-liquid two-phase flow state of the lead-based reactor under different operating conditions;

[0026] S106, Based on the core geometry model containing bubbles, the reaction cross section is calculated;

[0027] S108, homogenize the reaction cross section to generate the reaction cross section data.

[0028] Optionally, step S120, which involves calculating the core power distribution based on the reaction cross-section data and a core neutronics calculation model, includes:

[0029] S122, Fit the reaction cross-section data to obtain a continuous function of the reaction cross-section with respect to the cavitation fraction;

[0030] S124, based on the continuous function and the core neutron flux density, calculate the core neutron flux density; and based on the core neutron flux density, calculate the core power distribution.

[0031] Optionally, step S122, fitting the reaction cross-section data to obtain a continuous function of the reaction cross-section with respect to the cavitation fraction, includes:

[0032] S123, perform polynomial fitting on the reaction cross-section data to obtain a continuous function of the reaction cross-section with respect to the cavitation fraction; the continuous function is:

[0033] ;

[0034] In the formula, For the reaction cross section, This refers to the fraction of voids; These are the fitting coefficients. .

[0035] Optionally, after step S140, which involves obtaining the thermal parameters for the current time step based on the core power distribution and the core gas-liquid two-phase thermal-hydraulic calculation model, the method further includes:

[0036] S160, determine whether the core power distribution, the fuel rod temperature, the coolant temperature and the void fraction have converged. If so, proceed to the next time step and repeat steps S120 and S140.

[0037] Optionally, between step S140 and step S160, the method further includes:

[0038] S150: Determine whether the current time step has reached the preset time. If yes, output the result; otherwise, execute step S160.

[0039] This invention provides a method for calculating the gas-liquid two-phase nuclear-thermal coupling of lead-based reactor cores, which has at least the following beneficial technical effects:

[0040] By generating reaction cross-section data with multiple groups, multiple temperature points, and multiple void fractions, and transferring this data to the core neutronics calculation model, the core power distribution is calculated. This core power distribution is then transferred to the core gas-liquid two-phase thermal-hydraulic calculation model to obtain thermal parameters. This calculation method can accurately simulate the influence of gas phase on neutronics characteristics, and thus accurately predict the core power distribution and thermal parameters under gas-liquid two-phase flow conditions. This addresses the power transient safety issues that may arise in the core under severe accidents caused by gas-liquid two-phase flow, providing technical support for the safety analysis and accident assessment of lead-based reactors.

[0041] A second aspect of the present invention provides a lead-based reactor core gas-liquid two-phase nuclear thermal coupling calculation system, which applies the above-described lead-based reactor core gas-liquid two-phase nuclear thermal coupling calculation method, the system comprising:

[0042] The geometry model building module is used to build geometry models;

[0043] The cross-section data generation module is used to generate reaction cross-section data based on the core geometry model.

[0044] The core power distribution calculation module is used to calculate the core power distribution based on the reaction cross-section data and a core neutronics calculation model; and,

[0045] The thermal parameter calculation module obtains the thermal parameters for the current time step based on the core power distribution and the core gas-liquid two-phase thermal-hydraulic calculation model. The thermal parameters include fuel rod temperature, coolant temperature, and void fraction.

[0046] This invention provides a lead-based reactor core gas-liquid two-phase nuclear thermal coupling calculation system, which applies the lead-based reactor core gas-liquid two-phase nuclear thermal coupling calculation method described above and has all the advantages of the above coupling calculation method, which will not be repeated here. Attached Figure Description

[0047] Figure 1 A flowchart illustrating a lead-based reactor core gas-liquid two-phase nuclear-thermal coupling calculation method provided in this embodiment of the invention. Figure 1 ;

[0048] Figure 2 A flowchart illustrating the steps of constructing a core geometric model and generating reaction section data in a lead-based reactor core gas-liquid two-phase nuclear thermal coupling calculation method provided in an embodiment of the present invention.

[0049] Figure 3 A flowchart illustrating the steps of calculating the core power distribution based on the reaction cross-section data and the core neutronics calculation model in a lead-based reactor core gas-liquid two-phase nuclear thermal coupling calculation method provided in an embodiment of the present invention.

[0050] Figure 4 In a lead-based reactor core gas-liquid two-phase nuclear thermal coupling calculation method provided in an embodiment of the present invention, a flowchart is shown in which the thermal parameters of the current time step are obtained based on the core power distribution and the core gas-liquid two-phase thermal-hydraulic calculation model.

[0051] Figure 5 A flowchart illustrating a lead-based reactor core gas-liquid two-phase nuclear-thermal coupling calculation method provided in this embodiment of the invention. Figure 2 ;

[0052] Figure 6 The diagram shows the change of power distribution in a lead-based reactor core over time in a gas-liquid two-phase nuclear thermal coupling calculation method provided in this embodiment of the invention.

[0053] Figure 7 This is a block diagram of a lead-based reactor core gas-liquid two-phase nuclear thermal coupling calculation system provided in an embodiment of the present invention.

[0054] Explanation of reference numerals in the attached figures:

[0055] 10. Geometric model construction module; 20. Cross-section data generation module; 30. Core power distribution calculation module; 40. Thermal parameter calculation module. Detailed Implementation

[0056] To make the above-mentioned objects, features and advantages of the present invention more apparent and understandable, the following description is provided in conjunction with the accompanying drawings. Figures 1-7 Specific embodiments of the present invention will be described in detail below.

[0057] This invention provides a method for calculating the gas-liquid two-phase nuclear-thermal coupling of a lead-based reactor core, see appendix. Figure 1 The method includes:

[0058] S100, Construct the core geometry model and generate reaction cross-section data. In this step, the core geometry model includes the fuel rod region and the coolant region, and the reaction cross-section data includes multiple groups, multiple temperature points, and multiple void fractions. This step generates the reaction cross-section data based on OpenMC, an open-source Monte Carlo particle transport simulation software primarily used in nuclear reactor physics that supports geometric modeling. For details, see the appendix. Figure 2 The specific implementation process of step S100 is as follows:

[0059] Step S102: Construct the core geometry model; specifically, based on the Python API of OpenMC, construct a core geometry model consistent with the actual core structure; the core geometry model includes the fuel rod region and the coolant region; it should be noted that the core geometry model is a two-dimensional core geometry model, and the cross-section of the core geometry model is a regular hexagon.

[0060] Step S104: Insert bubbles into the coolant region to obtain a core geometry model containing bubbles; this setting is to simulate the gas-liquid two-phase flow state of the lead-based reactor under different operating conditions.

[0061] Step S106: Calculate the reaction cross section based on the core geometry model containing bubbles;

[0062] Step S108 involves homogenizing the reaction cross section to generate reaction cross section data, including multi-group, multi-temperature point, and multi-cavitation fraction data. Specifically, the reaction cross section is subjected to SPH correction to homogenize it, resulting in multi-group, multi-temperature point, and multi-cavitation fraction reaction cross section data. This setup yields a core geometry model containing bubbles and generates multi-group, multi-temperature point, and multi-cavitation fraction reaction cross section data, providing input for subsequent neutronics calculations. This allows for accurate simulation of the influence of gas phases on neutronics properties, and precise prediction of core power distribution and thermal parameters under gas-liquid two-phase flow conditions, providing technical support for the safety analysis and accident assessment of lead-based reactors.

[0063] S120, Based on the reaction cross-section data and the core neutronics calculation model, calculate the core power distribution; in this step, the reaction cross-section data is transferred to the Morphy neutronics calculation model to calculate the core power distribution; for details, see Appendix. Figure 3 The specific implementation process of step S120 is as follows:

[0064] Step S122 involves fitting the reaction cross-section data to obtain a continuous function of the reaction cross-section with respect to the cavitation fraction. It should be noted that the fitting method can be nonlinear fitting, polynomial fitting, or other existing fitting methods. For example, step S122 could be: S123, performing polynomial fitting on the reaction cross-section data to obtain a continuous function of the reaction cross-section with respect to the cavitation fraction; this continuous function is:

[0065] ;

[0066] In the formula, For the reaction cross section, This refers to the fraction of voids; These are the fitting coefficients. .

[0067] Step S124: Based on the aforementioned continuous function and the core neutron flux calculation model, calculate the core neutron flux density; and based on the core neutron flux density, calculate the core power distribution. It should be noted that this core neutron flux calculation model uses a neutron flux calculation model from the prior art, which will not be elaborated upon here.

[0068] S140, based on the core power distribution and the core gas-liquid two-phase thermal-hydraulic calculation model, the thermal parameters of the current time step are obtained, including fuel rod temperature, coolant temperature and void fraction.

[0069] The lead-based reactor core gas-liquid two-phase nuclear-thermal coupling calculation method provided in this invention generates reaction cross-section data with multiple groups, multiple temperature points, and multiple vacuolar fractions. This reaction cross-section data is then transferred to the core neutronics calculation model to calculate the core power distribution. This core power distribution is then transferred to the core gas-liquid two-phase thermal-hydraulic calculation model to obtain thermal parameters. This calculation method can accurately simulate the influence of gas phase on neutronics characteristics, and thus accurately predict the core power distribution and thermal parameters under gas-liquid two-phase flow conditions. This addresses the power transient safety issues that may arise from gas-liquid two-phase flow in the core under severe accidents, providing technical support for the safety analysis and accident assessment of lead-based reactors.

[0070] In this embodiment of the invention, in step S140, the core gas-liquid two-phase thermal-hydraulic calculation model includes a thermal-hydraulic model and a transport model. The thermal-hydraulic model is used to calculate the fuel rod temperature and coolant temperature. This thermal-hydraulic model adopts the thermal-hydraulic model in the prior art, which will not be described in detail here. The transport model includes the gas-liquid two-phase mass conservation equation, the gas-liquid two-phase momentum conservation equation, and the gas-liquid two-phase energy conservation equation, which are used to calculate the void fraction. It should be noted that the transport model is constructed based on the simplified assumptions of "no slip ratio in the gas-liquid two-phase system" and "no interphase heat transfer".

[0071] Specifically, since the velocities of the gas and liquid phases are equal, and the velocities of the gas phase and liquid phase are equal, the mass conservation equation for the gas and liquid phases is:

[0072] ;

[0073] In the formula, Let be the density of the gas-liquid two-phase mixture. For time, The velocity of the gas and liquid phases, For the fraction of voids, The density is the gas phase density. The density of the liquid phase;

[0074] Since there is no slip ratio between the gas and liquid phases, there is no relative motion between them. Therefore, the effects of interphase drag and lift can be ignored. Thus, the momentum conservation equation for the gas-liquid two-phase system is:

[0075] ;

[0076] In the formula, Let be the density of the gas-liquid two-phase mixture. For time, The velocity of the gas and liquid phases, For flow field pressure, The dynamic viscosity of the gas-liquid two-phase mixture is denoted as . For the fraction of voids, The dynamic viscosity of the gas phase. The dynamic viscosity of the liquid phase. It is the acceleration due to gravity;

[0077] Since there is "no interphase heat transfer" and "no slip ratio" between the gas and liquid phases, and their spatial distributions are synchronized, the energy conservation equation for the gas-liquid two-phase system can be simplified to:

[0078] ;

[0079] In the formula, Let be the density of the gas-liquid two-phase mixture. Enthalpy of a gas-liquid two-phase mixture For time, The velocity of the gas and liquid phases, Enthalpy of a gas-liquid two-phase mixture For the fraction of voids, For vapor phase enthalpy, For liquid phase enthalpy, Thermal conductivity, The characteristic temperature of the coolant. This is the rate at which heat is released from fuel to coolant per unit volume.

[0080] The lead-based reactor core gas-liquid two-phase nuclear thermal coupling calculation method provided in this embodiment of the invention can achieve refined simulation of gas-liquid two-phase flow and heat transfer in lead-based reactor cores by developing a deep fusion of transport model and thermal-hydraulic model.

[0081] In this embodiment of the invention, see appendix. Figure 4 Step S140, the step of obtaining the thermal parameters for the current time step based on the core power distribution and the core gas-liquid two-phase thermal-hydraulic calculation model, includes:

[0082] S142, based on the core power distribution and the thermal-hydraulic model, the fuel rod temperature and coolant temperature are calculated.

[0083] S144, solve the gas-liquid two-phase mass conservation equation, the gas-liquid two-phase momentum conservation equation, and the gas-liquid two-phase energy conservation equation in sequence to calculate the cavitation fraction.

[0084] In this embodiment of the invention, after step S140, which involves obtaining the thermal parameters for the current time step based on the core power distribution and the core gas-liquid two-phase thermo-hydraulic calculation model, please refer to the appendix. Figure 5 The methods also include:

[0085] S160: Determine whether the core power distribution, fuel rod temperature, coolant temperature, and void fraction have converged. That is, whether the change in core power distribution in adjacent time steps meets the preset power distribution difference, the change in fuel rod temperature in adjacent time steps meets the preset fuel rod temperature difference, the change in coolant temperature in adjacent time steps meets the preset coolant temperature difference, and the change in void fraction in adjacent time steps meets the preset void fraction difference. If yes, proceed to the next time step and repeat steps S120 and S140. If no, modify the reaction cross section data and recalculate the core power distribution based on the modified reaction cross section data and the core neutronics calculation model.

[0086] This setup is repeated iteratively until the changes in core power distribution, fuel rod temperature, and coolant temperature meet the convergence conditions, thereby achieving a precise coupled simulation of the physical processes and thermal-hydraulic processes of the core under gas-liquid two-phase conditions.

[0087] In this embodiment of the invention, see appendix. Figure 4 Between step S140 and step S160, the method further includes:

[0088] S150: Determine whether the current time step has reached the preset time. If yes, output the result; otherwise, execute step S160.

[0089] To further illustrate the present invention, the calculation method for gas-liquid two-phase nuclear thermal coupling of lead-based reactor core provided by the present invention will be described in more detail below with reference to the accompanying drawings and embodiments, but these should not be construed as limiting the scope of protection of the present invention.

[0090] Appendix Figure 6 This diagram illustrates the change in power distribution over time in a lead-based reactor core using a gas-liquid two-phase nuclear-thermal coupling calculation method provided in an embodiment of the invention. When bubbles enter the core, the core power distribution rapidly increases due to the increased fission cross-section and decreased absorption cross-section, reaching a maximum power of 180% of the original power. Later, the negative reactivity coefficient, resulting from the increased temperatures of both fuel rods and coolant, begins to take effect, and the power gradually decreases. The power distribution calculation results based on the gas-liquid two-phase nuclear-thermal coupling calculation method for lead-based reactor cores provided in this embodiment of the invention are consistent with expectations, namely, the power surge phenomenon caused by bubbles entering the core after a lead-based reactor accident.

[0091] This invention also provides a lead-based reactor core gas-liquid two-phase nuclear-thermal coupling calculation system, which applies the above-mentioned lead-based reactor core gas-liquid two-phase nuclear-thermal coupling calculation method, see appendix. Figure 7 The system includes a geometric model construction module 10, a cross-section data generation module 20, a core power distribution calculation module 30, and a thermal parameter calculation module 40. The geometric model construction module 10 is used to construct the geometric model; the cross-section data generation module 20 is used to generate reaction cross-section data based on the core geometric model; the core power distribution calculation module 30 is used to calculate the core power distribution based on the reaction cross-section data and the core neutronics calculation model; and the thermal parameter calculation module 40 is used to obtain the thermal parameters for the current time step based on the core power distribution and the core gas-liquid two-phase thermal-hydraulic calculation model. The thermal parameters include fuel rod temperature, coolant temperature, and void fraction.

[0092] The lead-based reactor core gas-liquid two-phase nuclear thermal coupling calculation system provided in this embodiment of the invention applies the above-mentioned lead-based reactor core gas-liquid two-phase nuclear thermal coupling calculation method and has all the beneficial effects of the calculation method, which will not be elaborated here.

[0093] While the present invention has been disclosed above, it is not limited thereto. Any person skilled in the art can make various modifications and alterations without departing from the spirit and scope of the invention; therefore, the scope of protection of the present invention should be determined by the scope defined in the claims.

Claims

1. A method for calculating the nuclear thermal coupling of gas-liquid two-phase phases in a lead-based reactor core, characterized in that, The method includes: S100, construct the core geometry model and generate reaction section data; wherein, the core geometry model includes the fuel rod region and the coolant region, and the reaction section data includes multiple groups, multiple temperature points and multiple cavitation fractions; S120, Based on the reaction cross-section data, calculate the core power distribution using a core neutronics calculation model; S140, based on the core power distribution and the core gas-liquid two-phase thermal-hydraulic calculation model, the thermal parameters for the current time step are obtained. These thermal parameters include fuel rod temperature, coolant temperature, and void fraction. The core gas-liquid two-phase thermal-hydraulic calculation model includes a thermal-hydraulic model and a transport model. The thermal-hydraulic model is used to calculate the fuel rod temperature and the coolant temperature. The transport model includes gas-liquid two-phase mass conservation equations, gas-liquid two-phase momentum conservation equations, and gas-liquid two-phase energy conservation equations, used to calculate the void fraction. The mass conservation equation for a gas-liquid two-phase system is: ； In the formula, The density of the gas-liquid two-phase mixture is denoted as . For time, The velocity of the gas and liquid phases, For the fraction of voids, The density is the gas phase density. The density of the liquid phase; The momentum conservation equation for a gas-liquid two-phase system is: ； In the formula, The density of the gas-liquid two-phase mixture is denoted as . For time, The velocity of the gas and liquid phases, For flow field pressure, The dynamic viscosity of the gas-liquid two-phase mixture is denoted as . For the fraction of voids, The dynamic viscosity of the gas phase. The dynamic viscosity of the liquid phase. It is the acceleration due to gravity; The energy conservation equation for a gas-liquid two-phase system is: ； In the formula, The density of the gas-liquid two-phase mixture is denoted as . Enthalpy of a gas-liquid two-phase mixture For time, The velocity of the gas and liquid phases, Enthalpy of a gas-liquid two-phase mixture For the fraction of voids, For vapor phase enthalpy, For liquid phase enthalpy, Thermal conductivity, The characteristic temperature of the coolant. This is the rate at which heat is released from fuel to coolant per unit volume.

2. The calculation method for gas-liquid two-phase nuclear thermal coupling in lead-based reactor cores according to claim 1, characterized in that, Step S140, which involves obtaining the thermal parameters for the current time step based on the core power distribution and the core gas-liquid two-phase thermal-hydraulic calculation model, includes the following steps: S142, Based on the core power distribution and the thermal-hydraulic model, the fuel rod temperature and the coolant temperature are calculated. S144, solve the gas-liquid two-phase mass conservation equation, gas-liquid two-phase momentum conservation equation, and gas-liquid two-phase energy conservation equation in sequence to calculate the cavitation fraction.

3. The calculation method for gas-liquid two-phase nuclear thermal coupling in lead-based reactor cores according to claim 1, characterized in that, Step S100, the step of constructing the core geometric model and generating reaction section data, includes: S102, Construct the core geometry model; S104, insert bubbles into the coolant region to obtain a core geometry model containing bubbles; S106, Based on the core geometry model containing bubbles, the reaction cross section is calculated; S108, homogenize the reaction cross section to generate the reaction cross section data.

4. The calculation method for gas-liquid two-phase nuclear thermal coupling of lead-based reactor cores according to any one of claims 1-3, characterized in that, Step S120, which involves calculating the core power distribution based on the reaction cross-section data and a core neutronics calculation model, includes the following steps: S122, Fit the reaction cross-section data to obtain a continuous function of the reaction cross-section with respect to the cavitation fraction; S124, based on the continuous function and the core neutron flux density, calculate the core neutron flux density; and based on the core neutron flux density, calculate the core power distribution.

5. The calculation method for gas-liquid two-phase nuclear thermal coupling of lead-based reactor cores according to claim 4, characterized in that, Step S122, which involves fitting the reaction cross-section data to obtain a continuous function of the reaction cross-section with respect to the cavitation fraction, includes: S123, perform polynomial fitting on the reaction cross-section data to obtain a continuous function of the reaction cross-section with respect to the cavitation fraction; the continuous function is: ； In the formula, For the reaction cross section, This refers to the fraction of voids; These are the fitting coefficients. .

6. The calculation method for gas-liquid two-phase nuclear thermal coupling of lead-based reactor cores according to any one of claims 1-3, characterized in that, After step S140, which involves obtaining the thermal parameters for the current time step based on the core power distribution and the core gas-liquid two-phase thermal-hydraulic calculation model, the method further includes: S160, determine whether the core power distribution, the fuel rod temperature, the coolant temperature and the void fraction have converged. If so, proceed to the next time step and repeat steps S120 and S140.

7. The calculation method for gas-liquid two-phase nuclear thermal coupling in lead-based reactor cores according to claim 6, characterized in that, Between step S140 and step S160, the method further includes: S150: Determine whether the current time step has reached the preset time. If yes, output the result; otherwise, execute step S160.

8. A lead-based reactor core gas-liquid two-phase nuclear thermal coupling calculation system, characterized in that, The system employing the lead-based reactor core gas-liquid two-phase nuclear-thermal coupling calculation method according to any one of claims 1-7 comprises: The geometric model building module (10) is used to build geometric models; The cross-section data generation module (20) is used to generate reaction cross-section data based on the core geometry model; The core power distribution calculation module (30) is used to calculate the core power distribution based on the reaction cross-section data and a core neutronics calculation model; and, The thermal parameter calculation module (40) obtains the thermal parameters of the current time step based on the core power distribution and the core gas-liquid two-phase thermal-hydraulic calculation model. The thermal parameters include fuel rod temperature, coolant temperature and void fraction.

Images