Aerosol removal result determination method and device coupled with multiple influence factors

Abstract

The invention relates to an aerosol removal result determination method and device coupled with multiple influence factors. The method comprises the following steps: acquiring thermal hydraulic parameters in a containment; determining a stress equation of the spray liquid drops in the containment, wherein the stress equation is used for indicating the stress equation of the spray liquid drops under the action of entrainment and settlement phenomena in the containment flow field; on the basis of the stress equation and the thermal hydraulic parameters, the spraying liquid drop state is obtained; determining the spatial distribution of the aerosol particles, wherein the spatial distribution of the aerosol particles is used for indicating the spatial distribution of the aerosol particles under the action of entrainment and settlement phenomena in the containment flow field; and based on the state of the spraying liquid drops and the spatial distribution of the aerosol particles, determining the removal efficiency of the spraying liquid drops to the aerosol particles so as to determine the removal result of the spraying liquid drops to the aerosol. The technical effects that the aerosol particle removal result is more accurate by considering the aerosol formation entrainment effect of the flow field in the containment and the existence of a natural sedimentation mechanism and coupling multiple influence factors are achieved.

Description

Technical Field

[0001] The present invention relates to the technical field of nuclear emergency treatment, and in particular to a method and a device for determining an aerosol removal result by coupling multiple influencing factors. Background Art

[0002] During a serious accident of a pressurized water reactor, a primary circuit breach or a benchmark accident such as a steam generator can cause core damage and release fission products into the containment, resulting in radioactive contamination. The spray system is one of the engineering safety facility systems. One of the functions of the spray system is to reduce the emission dose of radionuclides by spraying water, and to remove radioactive aerosol particles by different collection mechanisms through the action of radioactive aerosol particles. Therefore, the study of the mechanism of removing radioactive aerosols by spraying the containment can cope with the emergency response after the accident.

[0003] In the actual process of spraying the containment to remove radioactive aerosols, after the spraying is turned on, the spray droplets will have a cooling and decompression effect on the temperature and pressure inside the containment, thereby causing natural convection / forced convection inside the containment, resulting in local circulation of the flow field. The air flow movement in the free air space will entrain the radioactive aerosols to circulate, and at the same time, wall settlement will occur under gravity conditions.

[0004] In the related art, there is a problem that the aerosol particle removal results are inaccurate because the convection phenomenon of the flow field in the containment under the effect of spray cooling and pressure reduction is not taken into account, which will cause the aerosol to be entrained, and the radioactive aerosol will have a natural sedimentation mechanism.

[0005] The above problems need to be solved urgently. Summary of the invention

[0006] The invention discloses a method and a device for determining an aerosol removal result by coupling multiple influencing factors, aiming to solve the technical problems existing in the prior art.

[0007] The present invention adopts the following technical solutions:

[0008] On the one hand, the present invention provides a method for determining aerosol removal results by coupling multiple influencing factors, which includes: in an accident scenario, obtaining thermal-hydraulic parameters corresponding to computational domain nodes in a containment; determining a force equation of spray droplets in the containment, wherein the force equation is used to indicate the force equation of the spray droplets under the action of entrainment and sedimentation phenomena in the containment flow field; obtaining a spray droplet state based on the force equation and the thermal-hydraulic parameters; determining a spatial distribution of aerosol particles based on a variable transport model of aerosol particles in the containment, wherein the spatial distribution of aerosol particles is used to indicate the spatial distribution of aerosol particles under the action of entrainment and sedimentation phenomena in the containment flow field; determining a removal efficiency of the aerosol particles by the spray droplets based on the spray droplet state and the spatial distribution of the aerosol particles; and determining a removal result of the aerosol particles by the spray droplets based on the removal efficiency of the aerosol particles.

[0009] Optionally, in an accident scenario, obtaining thermal-hydraulic parameters corresponding to computational domain nodes in the containment shell includes: obtaining the total duration from the occurrence of the accident scenario to the end of the accident scenario, and the initial spatial position distribution of aerosol particles in the containment shell during the occurrence of the accident scenario; dividing the total duration into multiple time steps, dividing the computational domain of the containment shell into multiple computational domain nodes, and dividing the aerosol particles into multiple particle size intervals; selecting a first time step, and iteratively calculating the multiple computational domain nodes in the containment shell to obtain the thermal-hydraulic parameters corresponding to the first time step, wherein the first time step is the current time step among the multiple time steps.

[0010] Optionally, the first time step is selected, and the multiple computational domain nodes in the containment are iteratively calculated to obtain the thermal-hydraulic parameters corresponding to the first time step, including: determining the phase change heat of the spray droplets during evaporation and boiling on the multiple computational domain nodes in the containment to obtain a heat transfer model, wherein the heat transfer model is used to indicate the relationship between the mass change rate of the spray droplets and the vaporization heat phase change enthalpy; determining the boiling phase change of the spray droplets on the multiple computational domain nodes in the containment to obtain a mass transfer model, wherein the mass transfer model is used to indicate the phase change caused by boiling of the spray droplets when the temperature of the spray droplets reaches the boiling point or the water vapor pressure on the surface of the spray droplets is greater than the ambient pressure; determining the mass transfer heat transfer model based on the heat transfer model and the mass transfer model; solving the mass transfer heat transfer model to obtain the thermal-hydraulic parameters.

[0011] Optionally, determining the phase change heat of the spray droplets during evaporation and boiling to obtain a heat transfer model includes: the heat transfer model is calculated as follows:

[0012]

[0013] Among them, m p is the mass of the spray droplets, c p is the specific heat capacity of the spray droplets, T p is the temperature of the spray drop, t is the time, h fg is the phase change enthalpy of vaporization.

[0014] Optionally, determining the boiling phase change of the spray droplets to obtain a mass transfer model includes: the mass transfer model is calculated as follows:

[0015]

[0016] Among them, h fg is the phase change enthalpy of vaporization heat, h is the convection heat transfer coefficient, T ∞ is the temperature of high-temperature and high-pressure water vapor, T p is the temperature of the spray droplets.

[0017] Optionally, the variable transport model of aerosol particles in the containment is used to determine the spatial distribution of aerosol particles, including: determining the variable transport model of aerosol particles in the containment; when the aerosol particles are entrained in the free air space of the containment, adding a drag term to the component momentum transport equation in the variable transport model to obtain the entrainment effect equation of the aerosol particles; when the aerosol particles are subject to sedimentation in the free air space of the containment, adding a gravity term to the component momentum transport equation in the variable transport model to obtain the gravity effect equation of the aerosol particles; based on the variable transport model, the entrainment effect equation and the gravity effect equation, determining the spatial distribution of the aerosol particles under the effects of entrainment and sedimentation in the containment flow field.

[0018] Optionally, the adding of a drag term to the component momentum transport equation in the variable transport model to obtain the entrainment effect equation of the aerosol particles includes: the drag in the drag term obeys Stokes' law, and the drag is calculated as follows:

[0019]

[0020] Among them, F D is the drag coefficient, F D (uu P ) is the drag force, u is the velocity of the non-condensable gas containing water vapor, u P is the flow velocity of aerosol particles, μ is the fluid dynamic viscosity, C c is a function related to the mean free path of gas molecules and the radius of aerosol particles, ρ P is the density of aerosol particles, d pis the diameter of aerosol particles.

[0021] Optionally, the gravity term is added to the component momentum transport equation in the variable transport model to obtain the gravity effect equation of the aerosol particles, including: the gravity in the gravity term is calculated as follows:

[0022]

[0023] Among them: F g is the force in the direction of gravity after considering buoyancy, ρ P is the density of aerosol particles, ρ gas is the density of the mixture of aerosol particles and spray droplets.

[0024] Optionally, the removal efficiency of the aerosol particles by the spray droplets is determined based on the state of the spray droplets and the spatial distribution of the aerosol particles, including: determining a single droplet inertial capture efficiency when the aerosol particles are captured by the spray droplets by an inertial collision mechanism; determining a single droplet interception capture efficiency when the aerosol particles are captured by the spray droplets by an interception mechanism; determining a single droplet Brownian diffusion capture efficiency when the aerosol particles are captured by the spray droplets by a Brownian diffusion mechanism; determining the removal efficiency of the aerosol particles by the spray droplets based on the single droplet inertial capture efficiency, the single droplet interception capture efficiency and the single droplet Brownian diffusion capture efficiency.

[0025] Optionally, when the aerosol particles are captured by the spray droplets by the inertial collision mechanism, determining the single droplet inertial capture efficiency includes: calculating the single droplet inertial capture efficiency as follows:

[0026]

[0027] Among them, η imp is the inertial capture efficiency of a single droplet, St is the Stokes number, ρ p is the density of aerosol particles, d P is the diameter of the aerosol particles, d P is the diameter of the spray droplets, U G is the velocity of the free space gas in the containment, U D is the velocity of the spray droplets, μ G is the dynamic viscosity of the gas in the free space within the containment.

[0028] Optionally, when the aerosol particles are captured by the spray droplets by an interception mechanism, determining the single droplet interception capture efficiency includes: calculating the single droplet interception capture efficiency as follows:

[0029]

[0030] Among them, η R is the single droplet interception and capture efficiency, d P is the diameter of the aerosol particles, d D is the diameter of the spray droplets.

[0031] Optionally, when the aerosol particles are captured by the spray droplets using the Brownian diffusion mechanism, determining the single droplet Brownian diffusion capture efficiency includes: calculating the single droplet Brownian diffusion capture efficiency as follows:

[0032]

[0033] Among them, η dff is the single droplet interception and capture efficiency, ρ G is the density of the free space gas in the containment, k is the Boltzmann constant, μ G is the dynamic viscosity of the free space gas in the containment, C is the Cunningham correction coefficient, and T is the absolute temperature of the free space gas in the containment.

[0034] According to another aspect of an embodiment of the present invention, there is also provided an aerosol removal result determination device coupled with multiple influencing factors, including: a parameter acquisition module, used to obtain the thermal-hydraulic parameters of the containment corresponding to the computational domain nodes in the containment under an accident scenario; a force equation determination module, used to determine the force equation of the spray droplets in the containment, wherein the force equation is used to indicate the force equation of the spray droplets under the action of entrainment and sedimentation phenomena in the containment flow field; a spray droplet state module, used to obtain the spray droplet state based on the force equation and the thermal-hydraulic parameters. a particle spatial distribution module, for determining the spatial distribution of aerosol particles based on a variable transport model of aerosol particles in the containment, wherein the spatial distribution of aerosol particles is used to indicate the spatial distribution of aerosol particles under the action of entrainment and sedimentation phenomena in the containment flow field; a removal efficiency determination module, for determining the removal efficiency of the aerosol particles by the spray droplets based on the state of the spray droplets and the spatial distribution of the aerosol particles; a removal result module, for determining the removal result of the aerosol by the spray droplets based on the removal efficiency of the aerosol particles.

[0035] According to another aspect of an embodiment of the present invention, a non-volatile storage medium is provided, wherein the non-volatile storage medium stores a plurality of instructions, wherein the instructions are suitable for being loaded by a processor and executing any one of the methods for determining an aerosol removal result by coupling multiple influencing factors.

[0036] According to another aspect of an embodiment of the present invention, a computer program product is further provided, comprising a computer program, wherein when the computer program is executed by a processor, the steps of any one of the methods for determining an aerosol removal result by coupling multiple influencing factors are implemented.

[0037] The technical solution adopted by the present invention can achieve at least one of the following beneficial effects:

[0038] In an embodiment of the present invention, thermal-hydraulic parameters corresponding to computational domain nodes in a containment are obtained under an accident scenario; force equations of spray droplets in the containment are determined, wherein the force equations are used to indicate the force equations of the spray droplets under the effects of entrainment and sedimentation in the containment flow field; the state of spray droplets is obtained based on the force equations and the thermal-hydraulic parameters; the spatial distribution of aerosol particles is determined based on a variable transport model of aerosol particles in the containment, wherein the spatial distribution of aerosol particles is used to indicate the spatial distribution of aerosol particles under the effects of entrainment and sedimentation in the containment flow field; based on the state of spray droplets and the spatial distribution of aerosol particles, the removal efficiency of the aerosol particles by the spray droplets is determined; based on the removal efficiency of the aerosol particles, the removal result of the aerosol particles by the spray droplets is determined. The purpose of obtaining the state of spray droplets and the spatial distribution of aerosol particles under the action of entrainment and sedimentation phenomena in the containment flow field is achieved, thereby realizing the technical effect of taking into account that the convection phenomenon of the flow field in the containment will form an entrainment effect on the aerosol, and that there is a natural sedimentation mechanism for radioactive aerosols, and coupling multiple influencing factors to make the obtained aerosol particle removal results more accurate, and further solving the technical problem of inaccurate aerosol particle removal results due to failure to take into account the convection phenomenon of the flow field in the containment under the action of spray cooling and pressure reduction, which will form an entrainment effect on the aerosol, and that there is a natural sedimentation mechanism for radioactive aerosols. BRIEF DESCRIPTION OF THE DRAWINGS

[0039] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the following briefly introduces the drawings required for describing the embodiments, which constitute a part of the present invention. The exemplary embodiments of the present invention and their descriptions explain the present invention and do not constitute improper limitations on the present invention. In the drawings:

[0040] Figure 1 It is a flow chart of a method for determining an aerosol removal result by coupling multiple influencing factors in Example 1 of the present invention;

[0041] Figure 2 It is a distribution diagram of aerosol particles and spray droplets in a containment vessel of a method for determining radioactive aerosol diffusivitation in Example 1 of the present invention;

[0042] Figure 3 is a force diagram of aerosol particles in a method for determining radioactive aerosol diffusivitation motion in Example 1 of the present invention;

[0043] Figure 4 is an aerosol particle spray capture diagram of a radioactive aerosol diffusophore motion determination method in Example 1 of the present invention;

[0044] Figure 5 It is a flow chart of an optional method for determining aerosol removal results by coupling multiple influencing factors in Example 2 of the present invention;

[0045] Figure 6 It is a structural schematic diagram of an aerosol removal result determination device coupled with multiple influencing factors in Example 3 of the present invention. DETAILED DESCRIPTION

[0046] In order to make the purpose, technical solution and advantages of the present invention clearer, the technical solution of the present invention will be clearly and completely described below in conjunction with the specific embodiments of the present invention and the corresponding drawings. In the description of the present invention, it should be noted that the term "or" is usually used in the sense of including "and / or", unless the content clearly indicates otherwise.

[0047] In the description of the present invention, it should be noted that, unless otherwise clearly specified and limited, the terms "installed", "connected", and "connected" should be understood in a broad sense. For example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or a magnetic connection; it can be a direct connection, or it can be indirectly connected through an intermediate medium, or it can be the internal connection of two elements. For ordinary technicians in this field, the specific meanings of the above terms in the present invention can be understood according to specific circumstances. In addition, in the description of this application, the terms "first", "second", etc. are only used to distinguish the description, and cannot be understood as indicating or implying relative importance. In the description of the present invention, the meaning of "multiple" is at least two, such as two, three or more, etc., unless otherwise clearly and specifically limited.

[0048] Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0049] First, to facilitate understanding of the embodiments of the present invention, some terms or nouns involved in the present invention are explained below:

[0050] Coupling refers to the phenomenon that two or more systems or two forms of motion influence each other and even unite through interaction.

[0051] Aerosol is a colloidal dispersion system formed by small solid or liquid particles dispersed and suspended in a gas medium, also known as a gas dispersion system. Its dispersed phase is small solid or liquid particles, usually between 0.001 and 100 μm in size, but may also vary depending on the source and formation reasons.

[0052] In order to solve the problems existing in the prior art, the embodiments of the present application provide a method and device for determining an aerosol removal result by coupling multiple influencing factors.

[0053] Example 1

[0054] This embodiment provides a method for determining aerosol removal results by coupling multiple influencing factors, such as Figure 1 As shown, Figure 1 : is a flow chart of a method for determining an aerosol removal result by coupling multiple influencing factors in Example 1 of the present invention, the method comprising:

[0055] Step S102, in an accident scenario, obtaining thermal hydraulic parameters corresponding to the computational domain nodes in the containment;

[0056] Optionally, different types of components in the flow field of the containment include: non-condensable gas containing water vapor, radioactive aerosol particle components and spray droplets, and the physical modeling method of the Euler system is used for the non-condensable gas of water vapor, the physical modeling method of the Euler system is used for the radioactive aerosol particle components, and the physical modeling method of the Lagrangian system is used for the spray droplets. Specifically, the physical modeling method of the Euler system is a method based on node consideration of variable transport, and the particle tracking method of the Lagrangian system is a method based on calculation of spray droplet particle variables.

[0057] Optionally, the temperature and pressure reduction effect of the spray droplets on the containment environment is calculated by constructing a physical model of the interaction between different components. The calculation method is to use the heat and mass transfer model of the spray droplets and high-temperature water vapor for the spray droplets and the surrounding high-temperature and high-pressure water vapor under the Lagrangian system to calculate the temperature and pressure reduction process in the containment. The heat and mass transfer model represents the interaction between the non-condensable gas component containing water vapor and the spray droplet component. The specific calculation method is: the heat transfer in the containment is represented by calculating the phase change heat during the evaporation and boiling of the spray droplets, and the mass transfer in the containment is represented by calculating the boiling phase change of the spray droplets.

[0058] Optionally, based on the heat and mass transfer model, the thermal-hydraulic parameters corresponding to the computational domain nodes in the containment are calculated.

[0059] In some preferred embodiments, in an accident scenario, thermal-hydraulic parameters corresponding to the computational domain nodes in the containment are obtained, including: obtaining the total time from the occurrence of the accident scenario to the end of the accident scenario, and the initial spatial position distribution of aerosol particles in the containment during the accident scenario; dividing the total time into multiple time steps, dividing the computational domain of the containment into multiple computational domain nodes, and dividing the aerosol particles into multiple particle size intervals; selecting a first time step, and iteratively calculating the multiple computational domain nodes in the containment to obtain the thermal-hydraulic parameters corresponding to the first time step, wherein the first time step is the current time step among the multiple time steps.

[0060] Optionally, since the environment inside the containment is very complex under actual cooling conditions, the comprehensive effect of steam condensation on the removal of aerosol particles is relatively strong, and the diffusion and swimming movement of aerosol particles in actual conditions is directly related to the water vapor condensation condition on the containment wall, and is not a single spatial distribution, and the distribution of aerosol particles in the containment is not uniform, so it is necessary to divide the calculation domain inside the containment into multiple calculation domain nodes, and divide the aerosol particles into multiple particle size intervals, so as to effectively alleviate the influence of the very complex environment inside the containment under actual cooling conditions on the accuracy of the thermal-hydraulic parameter results finally obtained.

[0061] Optionally, the accident process is time-dependent. As time goes by, the accident scenario will change. The accident scenario is divided into multiple stages (time steps) according to the timeline, and the thermal-hydraulic parameters of each stage and the aerosol particle removal results corresponding to the thermal-hydraulic parameters are obtained respectively. The aerosol particle removal results of multiple stages are sorted out, thereby realizing the entire development process of aerosol particle removal, and effectively analyzing the aerosol particle removal results based on the entire development process.

[0062] Optionally, the first time step is taken as the current time step, and the corresponding thermal-hydraulic parameters in the first time step and the aerosol particle removal results determined based on the thermal-hydraulic parameters are obtained through an iterative loop; the next time step is updated as the current time step, and the iterative execution continues. If the calculation time physical time t (total duration) has been reached, the iterative calculation is stopped. At this time, the temperature, pressure and humidity field of the containment after spraying and the refined results of the removal of aerosol particles by spraying under the coupling effects of airflow entrainment and natural sedimentation have been obtained.

[0063] Optionally, the computational domain division of the containment should ensure that the multiple computational domain nodes are encrypted in the area where diffusion and phoresis mainly occur, and judged by the corresponding judgment criteria to ensure the applicability of the subsequent wall model. After the accident, a large amount of water vapor is generated to form a condensation process on the wall of the containment. The condensed water vapor will remove the aerosol particles. Therefore, the applicability of the wall model is the applicability of the computational domain of the containment. Figure 2 As shown, Figure 2 This is a distribution diagram of aerosol particles and spray droplets in a containment vessel of a method for determining radioactive aerosol diffusivitation in Example 1 of the present invention.

[0064] Optionally, node independence verification and time step independence verification are required for the computational domain node division process and the total time division process; among them, node independence verification mainly evaluates the sensitivity of the calculation results to the density of the computational domain nodes. In numerical calculations, the density of the computational domain nodes will directly affect the accuracy and efficiency of the calculation results. Time step independence verification evaluates the sensitivity of the calculation results to the size of the time step. In transient calculations, the selection of the time step affects the convergence and accuracy of the calculation.

[0065] Optionally, the aerosol is divided into different particle size intervals according to a source term particle size distribution density function, and the actual source term particle size distribution is fitted using a probability density function, which uses a log-normal distribution function or a Rosen-Rammler distribution function (Rosin-Rammler).

[0066] In some preferred embodiments, a first time step is selected, and multiple computational domain nodes within the containment are iteratively calculated to obtain thermal-hydraulic parameters corresponding to the first time step, including: determining the phase change heat of the spray droplets during evaporation and boiling at multiple computational domain nodes within the containment to obtain a heat transfer model, wherein the heat transfer model is used to indicate the relationship between the mass change rate of the spray droplets and the vaporization heat phase change enthalpy; determining the boiling phase change of the spray droplets at multiple computational domain nodes within the containment to obtain a mass transfer model, wherein the mass transfer model is used to indicate the phase change caused by boiling of the spray droplets when the temperature of the spray droplets reaches the boiling point or the water vapor pressure on the surface of the spray droplets is greater than the ambient pressure; determining the mass transfer heat transfer model based on the heat transfer model and the mass transfer model; solving the mass transfer heat transfer model to obtain thermal-hydraulic parameters.

[0067] Optionally, since the cooling and depressurizing effect of the spray droplets on the containment under the high-temperature and high-pressure steam environment is mainly the phenomenon of the latent heat of phase change taken away by the evaporation and boiling of the spray droplets, it is necessary to calculate the phase change heat during the evaporation and boiling of the spray droplets, that is, to calculate the phase change heat through the heat transfer model. The mass change rate of the spray droplets changes due to the mass transfer phenomenon caused by the condensation of high-temperature and high-pressure water vapor in the air on the spray droplets and the evaporation and boiling of the spray droplets. In order to obtain the mass change rate of the spray droplets, it is necessary to determine the mass transfer model of the spray droplets, that is, the mass transfer model based on the evaporation and condensation of the spray droplets and the mass transfer model of the boiling of the spray droplets. Among them, the mass transfer model of the boiling of the spray droplets means that when the temperature of the spray droplets reaches the boiling point, or the vapor pressure of the water on the surface of the spray droplets is greater than the ambient pressure, the spray droplets boil, and at this time, the heat obtained by the spray droplets from the humid air will be completely used for phase change.

[0068] In some preferred embodiments, the phase change heat of the spray droplets during evaporation and boiling is determined to obtain a heat transfer model, including: the heat transfer model is calculated as follows:

[0069]

[0070] Among them, m p is the mass of the spray droplets, c p is the specific heat capacity of the spray droplets, T p is the temperature of the spray drop, t is the time, h fg is the phase change enthalpy of vaporization.

[0071] In some preferred embodiments, the boiling phase change of the spray droplets is determined to obtain a mass transfer model, including: the mass transfer model is calculated as follows:

[0072]

[0073] Among them, h fg is the phase change enthalpy of vaporization heat, h is the convection heat transfer coefficient, T ∞ is the temperature of high-temperature and high-pressure water vapor, T p is the temperature of the spray droplets.

[0074] Optionally, the mass transfer model corresponding to the evaporation and condensation of the spray droplets is calculated as follows:

[0075]

[0076] Among them, k c is the mass transfer coefficient, A p is the surface area of ​​the spray droplets, ρ is the density of high-temperature and high-pressure water vapor air, B m is the spallation mass number.

[0077] Specifically, the calculation formula for the spallation mass number is as follows:

[0078]

[0079] Among them, Y i,s is the mass fraction of water vapor on the surface of the spray droplet, Y i,∞ is the mass fraction of water in the humid air corresponding to water vapor.

[0080] Optionally, the mass transfer model for the boiling of the spray droplets is when the temperature of the droplets reaches the boiling point T bp When the vapor pressure of water on the surface of the spray droplets is greater than the ambient pressure, the droplets boil, and all the heat obtained by the droplets from the humid air will be used for phase change.

[0081] The mass transfer model corresponding to the boiling of spray droplets is calculated as follows:

[0082]

[0083] Among them, h fg is the phase change enthalpy of vaporization heat, h is the convection heat transfer coefficient, T ∞ is the temperature of high-temperature and high-pressure water vapor, T p is the temperature of the spray droplets.

[0084] Specifically, the final required mass transfer model is determined based on the mass transfer model corresponding to the evaporation and condensation of the spray droplets and the mass transfer model corresponding to the boiling of the spray droplets, as well as the state of the spray droplets.

[0085] Step S104, determining the force equation of the spray droplets in the containment, wherein the force equation is used to indicate the force equation of the spray droplets under the action of entrainment and sedimentation phenomena in the containment flow field;

[0086] Step S106, obtaining the spray drop state based on the force equation and thermal hydraulic parameters;

[0087] Step S108, determining the spatial distribution of aerosol particles based on the variable transport model of aerosol particles in the containment, wherein the spatial distribution of aerosol particles is used to indicate the spatial distribution of aerosol particles under the effects of entrainment and sedimentation in the containment flow field;

[0088] Optionally, the variable transport model of radioactive aerosol particles needs to consider the drag force and the gravity direction force considering the buoyancy. The calculation process needs to combine thermal hydraulic parameters, as well as spray droplet state parameters such as spray droplet position and spray droplet size temperature.

[0089] Optionally, by constructing physical models between different components (using the physical modeling method of the Euler system for non-condensable gases such as water vapor, the physical modeling method of the Euler system for radioactive aerosol particle components, and the physical modeling method of the Lagrangian system for spray droplets), the entrainment and sedimentation of radioactive aerosol particles in the containment flow field are considered, which mainly considers the air entrainment effect and gravity sedimentation of radioactive aerosol particles in the free air space. The specific consideration method is to add the force under the corresponding transport mechanism to the variable transport model of the radioactive aerosol particle components under the Euler system to enable the radioactive aerosol particles to achieve the corresponding mechanism transport behavior.

[0090] In some preferred embodiments, the spatial distribution of aerosol particles is determined based on the variable transport model of aerosol particles in the containment, including: determining the variable transport model of aerosol particles in the containment; when the aerosol particles are entrained in the free air space of the containment, adding a drag term to the component momentum transport equation in the variable transport model to obtain the entrainment effect equation of the aerosol particles; when the aerosol particles are subjected to sedimentation in the free air space of the containment, adding a gravity term to the component momentum transport equation in the variable transport model to obtain the gravity effect equation of the aerosol particles; based on the variable transport model, the entrainment effect equation and the gravity effect equation, determining the spatial distribution of aerosol particles under the effect of entrainment and sedimentation in the containment flow field. Figure 3 As shown, Figure 3 It is a force diagram of aerosol particles in a method for determining radioactive aerosol diffusivitation motion in Example 1 of the present invention.

[0091] In some preferred embodiments, a drag term is added to the component momentum transport equation in the variable transport model to obtain the entrainment effect equation of the aerosol particles, including: the drag in the drag term obeys Stokes' law, and the drag is calculated as follows:

[0092]

[0093] Among them, F D is the drag coefficient, F D (uu P ) is the drag force, u is the velocity of the non-condensable gas containing water vapor, u P is the flow velocity of aerosol particles, μ is the fluid dynamic viscosity, C c is a function related to the mean free path of gas molecules and the radius of aerosol particles, ρ P is the density of aerosol particles, d p is the diameter of aerosol particles.

[0094] Optionally, the drag term should be added separately in three dimensions (the force equation of the non-condensable gas of water vapor, the force equation of the radioactive aerosol particle components, and the force equation of the spray droplets) to realize the calculation of the entrainment effect of unsteady airflow on radioactive aerosol particles and eliminate the interference of other gas components on the force of radioactive aerosol particles.

[0095] In some preferred embodiments, a gravity term is added to the component momentum transport equation in the variable transport model to obtain the gravity effect equation of the aerosol particles, including: the gravity in the gravity term is calculated as follows:

[0096]

[0097] Among them: F g is the force in the direction of gravity after considering buoyancy, ρ P is the density of aerosol particles, ρ gas is the density of the mixture of aerosol particles and spray droplets.

[0098] Step S110, determining the removal efficiency of the spray droplets on the aerosol particles based on the spray droplet state and the spatial distribution of the aerosol particles;

[0099] Optionally, by constructing physical models between different components (using the physical modeling method of the Euler system for non-condensable gases such as water vapor, using the physical modeling method of the Euler system for radioactive aerosol particle components, and using the physical modeling method of the Lagrangian system for spray droplets), the entrainment and sedimentation phenomena of radioactive aerosol particles in the containment flow field are considered, and the main considerations are: spray droplets capture aerosols by inertial collision, interception, and Brownian diffusion mechanisms for aerosol particles of different particle sizes. Specifically, the capture of radioactive aerosol particles by the inertial collision mechanism is achieved by constructing a single droplet inertial capture efficiency model for spray droplets and surrounding radioactive aerosol components under the Lagrangian system. The capture of radioactive aerosol particles by the interception mechanism is achieved by constructing a single droplet interception capture efficiency model for spray droplets and surrounding radioactive aerosol components under the Lagrangian system. The capture of radioactive aerosol particles by Brownian diffusion mechanism is achieved by constructing a single droplet Brownian diffusion capture efficiency model for the spray droplets and surrounding radioactive aerosol components under the Lagrangian system. Figure 4 As shown, Figure 4 This is an aerosol particle spray capture diagram of a method for determining radioactive aerosol diffusivitation movement in Example 1 of the present invention.

[0100] In some preferred embodiments, the removal efficiency of aerosol particles by spray droplets is determined based on the state of spray droplets and the spatial distribution of aerosol particles, including: when aerosol particles are captured by spray droplets by inertial collision mechanism, determining the single droplet inertial capture efficiency; when aerosol particles are captured by spray droplets by interception mechanism, determining the single droplet interception capture efficiency; when aerosol particles are captured by spray droplets by Brownian diffusion mechanism, determining the single droplet Brownian diffusion capture efficiency; based on the single droplet inertial capture efficiency, the single droplet interception capture efficiency and the single droplet Brownian diffusion capture efficiency, the removal efficiency of aerosol particles by spray droplets is determined.

[0101] In some preferred embodiments, when aerosol particles are captured by spray droplets by an inertial collision mechanism, determining the single droplet inertial capture efficiency includes: calculating the single droplet inertial capture efficiency as follows:

[0102]

[0103] Among them, η imp is the inertial capture efficiency of a single droplet, St is the Stokes number, ρ p is the density of aerosol particles, d P is the diameter of the aerosol particles, d P is the diameter of the spray droplets, U G is the velocity of the free space gas in the containment, U D is the velocity of the spray droplets, μ G is the dynamic viscosity of the gas in the free space within the containment.

[0104] In some preferred embodiments, when aerosol particles are captured by spray droplets by an interception mechanism, determining the single droplet interception capture efficiency includes: calculating the single droplet interception capture efficiency as follows:

[0105]

[0106] Among them, η R is the single droplet interception and capture efficiency, d P is the diameter of the aerosol particles, d D is the diameter of the spray droplets.

[0107] In some preferred embodiments, when aerosol particles are captured by spray droplets using a Brownian diffusion mechanism, determining the single droplet Brownian diffusion capture efficiency includes: calculating the single droplet Brownian diffusion capture efficiency as follows:

[0108]

[0109] Among them, η dff is the single droplet interception and capture efficiency, ρ Gis the density of the free space gas in the containment, k is the Boltzmann constant, μ G is the dynamic viscosity of the free space gas in the containment, C is the Cunningham correction coefficient, and T is the absolute temperature of the free space gas in the containment.

[0110] Step S112, based on the removal efficiency of aerosol particles, determining the removal result of aerosol particles by the spray droplets.

[0111] Optionally, the removal results of aerosol particles by spray droplets are mainly affected by factors such as the characteristics of the spray droplets (such as size and concentration) and operating conditions (such as wind speed and spray intensity). Based on the single droplet inertial capture efficiency, single droplet interception capture efficiency and single droplet Brownian diffusion capture efficiency obtained above, the removal efficiency of aerosol particles by spray droplets is determined. Within the total duration of the accident scenario, the removal efficiency of aerosol particles by spray droplets corresponding to each time step is determined, so that the removal results of aerosol particles in the current accident scenario can be determined. Based on the removal results, the aerosol particles can be effectively analyzed to determine the reaction process of the current accident scenario and reduce the radiation hazard of aerosol particles.

[0112] Through the above steps S102 to S112, the purpose of obtaining the state of spray droplets and the spatial distribution of aerosol particles under the action of entrainment and sedimentation phenomena in the containment flow field is achieved, thereby realizing the technical effect of taking into account that the convection phenomenon of the flow field in the containment will form an entrainment effect on the aerosol, and that there will be a natural sedimentation mechanism for radioactive aerosols, and coupling multiple influencing factors to make the obtained aerosol particle removal result more accurate, thereby solving the technical problem of inaccurate aerosol particle removal result due to failure to take into account the convection phenomenon of the flow field in the containment under the action of spray cooling and pressure reduction, which will form an entrainment effect on the aerosol, and that there will be a natural sedimentation mechanism for radioactive aerosols.

[0113] Example 2

[0114] Based on the above embodiments and optional embodiments, the present invention also proposes an optional implementation mode: Figure 5 is a flow chart of an optional method for determining aerosol removal results by coupling multiple influencing factors in Example 2 of the present invention, such as Figure 5 As shown, the method includes:

[0115] Step S1, dividing the containment calculation domain of the accident object into a number of calculation domain nodes, dividing the physical time t during the accident into different time steps (t0, t1...tn), and dividing the aerosol particles into different particle size intervals;

[0116] Step S2, iteratively solving the transport equation, component transport model and simplified mass transfer and heat transfer model for the containment calculation domain nodes at the current time step (first time step) to obtain the thermal hydraulic parameters of the containment under condensation conditions at the current time step;

[0117] Step S3, solving the force equation of the spray droplets based on the Lagrangian system at the current time step to obtain the spatial position distribution of the spray droplets, and updating the temperature and mass of the spray droplets based on the thermal hydraulic parameters calculated by the heat and mass transfer model, as well as updating the state of the spray droplets under the current containment environment;

[0118] Step S4, iteratively solving the variable transport model of the radioactive aerosol particle component for the containment calculation domain node at the current time step, and obtaining the spatial distribution of the radioactive aerosol in different particle size intervals of the radioactive aerosol particle component at the current time step; at the same time, in the calculation domain node where the radioactive aerosol particle component and the spray droplet exist at the same time, the spray capture model of the radioactive aerosol particle component is solved, and the solution result is added as a source term into the transport equation of the radioactive aerosol particle component, so as to realize the capture calculation of the radioactive aerosol under the spraying action, and update the mass of the radioactive aerosol in different particle size intervals and the spatial distribution of the aerosol particles after transport;

[0119] Step S5, update the next time step as the current time step, iterate the thermal hydraulic parameters, spray droplet state and spatial distribution of aerosol particles, and stop the iterative calculation when the calculation time physical time t (total duration) has been reached. At this time, the refined results of aerosol removal by spraying under the coupling effect of temperature, pressure and humidity field of spray droplets and airflow entrainment and natural sedimentation of the containment have been obtained.

[0120] Through the above steps S1 to S5, the non-condensable gas containing water vapor, the radioactive aerosol components and the spray droplets are modeled using the Eulerian system, the Eulerian system and the Lagrangian system, respectively, and three interaction models are established, among which the heat and mass transfer model of the spray droplets and the high-temperature water vapor can consider the cooling and pressure reduction effect of the spray droplets on the containment environment; the method of adding the force under the corresponding transport mechanism to the variable transport model of the radioactive aerosol components can consider the entrainment and sedimentation phenomena of the radioactive aerosol in the flow field; the mechanism model of the spray capture aerosol can consider the removal of radioactive aerosols of different particle sizes by the spray.

[0121] At the same time, the coupling can be solved iteratively over time on multiple computational domain nodes to consider the cooling and depressurization effects of spraying droplets on the containment environment, the entrainment and sedimentation of radioactive aerosols in the flow field, and the removal of radioactive aerosols by spraying, thereby achieving a detailed spatiotemporal analysis of the radioactive aerosol spray removal mechanism under the coupling effects of multiple influencing factors (entrainment, sedimentation, and spraying).

[0122] Example 3

[0123] According to an embodiment of the present invention, there is also provided an embodiment of a device for implementing the above-mentioned method for determining aerosol removal results by coupling multiple influencing factors. Figure 6 is a schematic diagram of the structure of an aerosol removal result determination device coupled with multiple influencing factors in Example 3 of the present invention, such as Figure 6 As shown, the above-mentioned aerosol removal result determination device coupled with multiple influencing factors includes: a parameter acquisition module 401, a force equation determination module 402, a spray droplet state module 403, a particle space distribution module 404, a removal efficiency determination module 405 and a removal result module 406, wherein:

[0124] The parameter acquisition module 401 is used to obtain the thermal hydraulic parameters of the containment corresponding to the computational domain nodes in the containment under the accident scenario;

[0125] A force equation determination module 402, connected to the parameter acquisition module 401, is used to determine the force equation of the spray droplets in the containment, wherein the force equation is used to indicate the force equation of the spray droplets under the effect of entrainment and sedimentation in the containment flow field;

[0126] The spray drop state module 403 is connected to the force equation determination module 402 and is used to obtain the spray drop state based on the force equation and thermal hydraulic parameters;

[0127] The particle spatial distribution module 404 is connected to the spray droplet state module 403 and is used to determine the spatial distribution of aerosol particles based on the variable transport model of aerosol particles in the containment, wherein the spatial distribution of aerosol particles is used to indicate the spatial distribution of aerosol particles under the effect of entrainment and sedimentation in the containment flow field;

[0128] The removal efficiency determination module 405 is connected to the particle spatial distribution module 404 and is used to determine the removal efficiency of the spray droplets on the aerosol particles based on the spray droplet state and the spatial distribution of the aerosol particles;

[0129] The removal result module 406 is connected to the removal efficiency determination module 405 and is used to determine the removal result of the aerosol by the spray droplets based on the removal efficiency of the aerosol particles.

[0130] It should be noted that the above modules can be implemented by software or hardware. For example, for the latter, it can be implemented in the following ways: the above modules can be located in the same processor; or the above modules can be located in different processors in any combination.

[0131] It should be noted that the above-mentioned measurement parameter acquisition module 401, force equation determination module 402, spray droplet state module 403, particle space distribution module 404, removal efficiency determination module 405 and removal result module 406 correspond to steps S102 to S112 in the embodiment, and the examples and application scenarios implemented by the above-mentioned modules and corresponding steps are the same, but are not limited to the contents disclosed in the above-mentioned embodiment. It should be noted that the above-mentioned modules can be run in a computer terminal as part of the device.

[0132] It should be noted that the optional or preferred implementation of this embodiment can refer to the relevant description in the embodiment, which will not be repeated here.

[0133] The above-mentioned aerosol removal result determination device coupled with multiple influencing factors can also include a processor and a memory. The above-mentioned measurement parameter acquisition module 401, force equation determination module 402, spray droplet state module 403, particle space distribution module 404, removal efficiency determination module 405 and removal result module 406 are all stored in the memory as program modules, and the processor executes the above-mentioned program modules stored in the memory to realize the corresponding functions.

[0134] The processor includes a kernel, which retrieves the corresponding program module from the memory. The kernel may be one or more. The memory may include non-permanent memory in a computer-readable medium, random access memory (RAM) and / or non-volatile memory, such as read-only memory (ROM) or flash RAM. The memory includes at least one memory chip.

[0135] According to an embodiment of the present application, an embodiment of a non-volatile storage medium is also provided. Optionally, in this embodiment, the non-volatile storage medium includes a stored program, wherein when the program is running, the device where the non-volatile storage medium is located is controlled to execute any of the above-mentioned methods for determining aerosol removal results by coupling multiple influencing factors.

[0136] Optionally, in this embodiment, the non-volatile storage medium may be located in any computer terminal in a computer terminal group in a computer network, or in any mobile terminal in a mobile terminal group, and the non-volatile storage medium includes a stored program.

[0137] Optionally, when the program is running, the device where the non-volatile storage medium is located is controlled to perform the following functions: in an accident scenario, obtain the thermal-hydraulic parameters corresponding to the computational domain nodes in the containment; determine the force equation of the spray droplets in the containment, wherein the force equation is used to indicate the force equation of the spray droplets under the action of entrainment and sedimentation in the containment flow field; obtain the state of the spray droplets based on the force equation and the thermal-hydraulic parameters; determine the spatial distribution of aerosol particles based on the variable transport model of aerosol particles in the containment, wherein the spatial distribution of aerosol particles is used to indicate the spatial distribution of aerosol particles under the action of entrainment and sedimentation in the containment flow field; determine the removal efficiency of aerosol particles by spray droplets based on the state of spray droplets and the spatial distribution of aerosol particles; determine the removal result of aerosol particles by spray droplets based on the removal efficiency of aerosol particles.

[0138] According to an embodiment of the present application, an embodiment of a processor is also provided. Optionally, in this embodiment, the processor is used to run a program, wherein when the program is run, any one of the above-mentioned methods for determining aerosol removal results by coupling multiple influencing factors is executed.

[0139] According to an embodiment of the present application, an embodiment of a computer program product is also provided. Optionally, in this embodiment, the computer program product includes a computer program, which, when executed by a processor, implements any of the steps of the aerosol removal result determination method coupled with multiple influencing factors.

[0140] Optionally, the above-mentioned computer program product, when executed on a data processing device, is suitable for executing an initialization program having the following method steps: in an accident scenario, obtaining thermal-hydraulic parameters corresponding to the computational domain nodes in the containment; determining the force equation of the spray droplets in the containment, wherein the force equation is used to indicate the force equation of the spray droplets under the action of entrainment and sedimentation phenomena in the containment flow field; obtaining the state of the spray droplets based on the force equation and the thermal-hydraulic parameters; determining the spatial distribution of aerosol particles based on a variable transport model of aerosol particles in the containment, wherein the spatial distribution of aerosol particles is used to indicate the spatial distribution of aerosol particles under the action of entrainment and sedimentation phenomena in the containment flow field; determining the removal efficiency of aerosol particles by spray droplets based on the state of spray droplets and the spatial distribution of aerosol particles; and determining the removal result of aerosol particles by spray droplets based on the removal efficiency of aerosol particles.

[0141] An embodiment of the present invention provides an electronic device, which includes a processor, a memory, and a program stored in the memory and executable on the processor, wherein the processor implements the following steps when executing the program: in an accident scenario, obtaining thermal-hydraulic parameters corresponding to computational domain nodes in a containment; determining a force equation of spray droplets in the containment, wherein the force equation is used to indicate the force equation of spray droplets under the effects of entrainment and sedimentation in a containment flow field; obtaining a state of spray droplets based on the force equation and the thermal-hydraulic parameters; determining a spatial distribution of aerosol particles based on a variable transport model of aerosol particles in the containment, wherein the spatial distribution of aerosol particles is used to indicate the spatial distribution of aerosol particles under the effects of entrainment and sedimentation in a containment flow field; determining a removal efficiency of aerosol particles by spray droplets based on the state of spray droplets and the spatial distribution of aerosol particles; and determining a removal result of aerosol particles by spray droplets based on the removal efficiency of aerosol particles.

[0142] The above sequence of the embodiments of the present invention is for description only and does not represent the superiority or inferiority of the embodiments.

[0143] In the above embodiments of the present invention, the description of each embodiment has its own emphasis. For parts that are not described in detail in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

[0144] In the several embodiments provided in this application, it should be understood that the disclosed technical content can be implemented in other ways. Among them, the device embodiments described above are only schematic. For example, the division of the above modules can be a logical function division. There may be other division methods in actual implementation, such as multiple modules or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some interfaces, modules or indirect coupling or communication connection of modules, which can be electrical or other forms.

[0145] The modules described above as separate components may or may not be physically separated, and the components shown as modules may or may not be physical modules, that is, they may be located in one place or distributed on multiple modules. Some or all of the modules may be selected according to actual needs to achieve the purpose of the present embodiment.

[0146] In addition, each functional module in each embodiment of the present invention may be integrated into one processing module, or each module may exist physically separately, or two or more modules may be integrated into one module. The above-mentioned integrated module may be implemented in the form of hardware or in the form of software functional modules.

[0147] If the above-mentioned integrated module is implemented in the form of a software function module and sold or used as an independent product, it can be stored in a computer-readable non-volatile storage medium. Based on this understanding, the technical solution of the present invention is essentially or the part that contributes to the prior art or all or part of the technical solution can be embodied in the form of a software product. The computer software product is stored in a non-volatile storage medium, including a number of instructions to enable a computer device (which can be a personal computer, server or network device, etc.) to perform all or part of the steps of the various embodiments of the present invention. The aforementioned non-volatile storage medium includes: U disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), mobile hard disk, magnetic disk or optical disk and other media that can store program codes.

[0148] The above are only preferred embodiments of the present invention. It should be pointed out that, for ordinary technicians in this technical field, several improvements and modifications can be made without departing from the principles of the present invention. These improvements and modifications should also be regarded as the scope of protection of the present invention.

Claims

1. A method for determining aerosol removal results by coupling multiple influencing factors, characterized in that: include: Under accident scenarios, obtain the thermal-hydraulic parameters corresponding to the computational domain nodes in the containment; Determining a force equation of the spray droplets in the containment, wherein the force equation is used to indicate the force equation of the spray droplets under the effects of entrainment and sedimentation in the containment flow field; Based on the force equation and the thermal hydraulic parameters, a spray drop state is obtained; Determining the spatial distribution of aerosol particles based on a variable transport model of aerosol particles in the containment, wherein the spatial distribution of aerosol particles is used to indicate the spatial distribution of aerosol particles under the effects of entrainment and sedimentation in the containment flow field; Determining the removal efficiency of the aerosol particles by the spray droplets based on the state of the spray droplets and the spatial distribution of the aerosol particles; Based on the removal efficiency of the aerosol particles, the removal result of the aerosol particles by the spray droplets is determined.

2. The method for determining aerosol removal results by coupling multiple influencing factors according to claim 1, characterized in that: In the accident scenario, obtaining the thermal hydraulic parameters corresponding to the computational domain nodes in the containment shell includes: Obtain the total duration from the start of the accident scenario to the end of the accident scenario, as well as the initial spatial position distribution of aerosol particles in the containment during the accident scenario; Dividing the total time into a plurality of time steps, dividing the computational domain of the containment into a plurality of computational domain nodes, and dividing the aerosol particles into a plurality of particle size intervals; A first time step is selected, and the multiple computational domain nodes in the containment are iteratively calculated to obtain thermal-hydraulic parameters corresponding to the first time step, wherein the first time step is a current time step among the multiple time steps.

3. The method for determining aerosol removal results by coupling multiple influencing factors according to claim 2, characterized in that: The selecting of the first time step, performing iterative calculation on the plurality of computational domain nodes in the containment, and obtaining thermal hydraulic parameters corresponding to the first time step, comprises: Determine the phase change heat of the spray droplets during evaporation and boiling on the plurality of computational domain nodes in the containment to obtain a heat transfer model, wherein the heat transfer model is used to indicate the relationship between the mass change rate of the spray droplets and the vaporization heat transfer phase change enthalpy; Determine the boiling phase change of the spray droplets on the multiple computational domain nodes in the containment to obtain a mass transfer model, wherein the mass transfer model is used to indicate the phase change caused by boiling of the spray droplets when the temperature of the spray droplets reaches the boiling point or the water vapor pressure on the surface of the spray droplets is greater than the ambient pressure; Determining a mass transfer and heat transfer model based on the heat transfer model and the mass transfer model; The mass transfer and heat transfer model is solved to obtain the thermal hydraulic parameters.

4. The method for determining aerosol removal results by coupling multiple influencing factors according to claim 3, characterized in that: The step of determining the phase change heat of the spray droplets during evaporation and boiling to obtain a heat transfer model comprises: The heat transfer model calculation is as follows: Among them, m p is the mass of the spray droplets, c p is the specific heat capacity of the spray droplets, T p is the temperature of the spray drop, t is the time, h fg is the phase change enthalpy of vaporization.

5. The method for determining aerosol removal results by coupling multiple influencing factors according to claim 3, characterized in that: The step of determining the boiling phase change of the spray droplets to obtain a mass transfer model comprises: The mass transfer model calculation is as follows: Among them, h fg is the phase change enthalpy of vaporization heat, h is the convection heat transfer coefficient, T ∞ is the temperature of high-temperature and high-pressure water vapor, T p is the temperature of the spray droplets.

6. The method for determining aerosol removal results by coupling multiple influencing factors according to claim 1, characterized in that: The variable transport model of aerosol particles in the containment is used to determine the spatial distribution of aerosol particles, including: Determine the variable transport model of aerosol particles within the containment; In the case where the aerosol particles are entrained in the free gas space of the containment, a drag term is added to the component momentum transport equation in the variable transport model to obtain the entrainment effect equation of the aerosol particles; When the aerosol particles are subjected to sedimentation in the free gas space of the containment, a gravity term is added to the component momentum transport equation in the variable transport model to obtain a gravity effect equation of the aerosol particles; Based on the variable transport model, the entrainment effect equation and the gravity effect equation, the spatial distribution of the aerosol particles under the effects of entrainment and sedimentation phenomena in the containment flow field is determined.

7. The method for determining aerosol removal results by coupling multiple influencing factors according to claim 6, characterized in that: The drag term is added to the component momentum transport equation in the variable transport model to obtain the entrainment effect equation of the aerosol particles, including: The drag force in the drag term obeys Stokes' law and is calculated as follows: Among them, F D is the drag coefficient, F D (uu P ) is the drag force, u is the velocity of the non-condensable gas containing water vapor, u P is the flow velocity of aerosol particles, μ is the fluid dynamic viscosity, C c is a function related to the mean free path of gas molecules and the radius of aerosol particles, ρ P is the density of aerosol particles, d p is the diameter of aerosol particles.

8. The method for determining aerosol removal results by coupling multiple influencing factors according to claim 6, characterized in that: The gravity term is added to the component momentum transport equation in the variable transport model to obtain the gravity effect equation of the aerosol particles, including: The gravity in the gravity term is calculated as follows: Among them: F g is the force in the direction of gravity after considering buoyancy, ρ P is the density of aerosol particles, ρ gas is the density of the mixture of aerosol particles and spray droplets.

9. The method for determining aerosol removal results by coupling multiple influencing factors according to claim 1, characterized in that: The step of determining the removal efficiency of the aerosol particles by the spray droplets based on the state of the spray droplets and the spatial distribution of the aerosol particles comprises: determining a single droplet inertial capture efficiency when the aerosol particles are captured by the spray droplets by an inertial collision mechanism; Determining the single droplet interception capture efficiency when the aerosol particles are captured by the spray droplets by an interception mechanism; determining a single droplet Brownian diffusion capture efficiency when the aerosol particles are captured by the spray droplets by a Brownian diffusion mechanism; The removal efficiency of the aerosol particles by the spray droplets is determined based on the single droplet inertial capture efficiency, the single droplet interception capture efficiency, and the single droplet Brownian diffusion capture efficiency.

10. The method for determining aerosol removal results by coupling multiple influencing factors according to claim 9, characterized in that: The method of determining the inertial capture efficiency of a single droplet when the aerosol particles are captured by the spray droplets by an inertial collision mechanism comprises: The calculation of the single droplet inertial capture efficiency is as follows: Among them, η imp is the inertial capture efficiency of a single droplet, St is the Stokes number, ρ p is the density of aerosol particles, d P is the diameter of the aerosol particles, d P is the diameter of the spray droplets, U G is the velocity of the free space gas in the containment, U D is the velocity of the spray droplets, μ G is the dynamic viscosity of the gas in the free space within the containment.

11. The method for determining aerosol removal results by coupling multiple influencing factors according to claim 9, characterized in that: The method of determining the single droplet interception and capture efficiency when the aerosol particles are captured by the spray droplets by the interception mechanism comprises: The single droplet interception capture efficiency is calculated as follows: Among them, η R is the single droplet interception and capture efficiency, d P is the diameter of the aerosol particles, d D is the diameter of the spray droplets.

12. The method for determining aerosol removal results by coupling multiple influencing factors according to claim 9, characterized in that: The method of determining the single droplet Brownian diffusion capture efficiency when the aerosol particles are captured by the spray droplets using the Brownian diffusion mechanism comprises: The single droplet Brownian diffusion capture efficiency is calculated as follows: Among them, η dff is the single droplet interception and capture efficiency, ρ G is the density of the free space gas in the containment, k is the Boltzmann constant, μ G is the dynamic viscosity of the free space gas in the containment, C is the Cunningham correction coefficient, and T is the absolute temperature of the free space gas in the containment.

13. An aerosol removal result determination device coupled with multiple influencing factors, characterized in that: include: The parameter acquisition module is used to obtain the thermal hydraulic parameters of the containment corresponding to the computational domain nodes in the containment under accident scenarios; A force equation determination module, used to determine the force equation of the spray droplets in the containment, wherein the force equation is used to indicate the force equation of the spray droplets under the effect of entrainment and sedimentation in the containment flow field; A spray drop state module, used for obtaining the spray drop state based on the force equation and the thermal hydraulic parameters; A particle spatial distribution module, used to determine the spatial distribution of aerosol particles based on a variable transport model of aerosol particles in the containment, wherein the spatial distribution of aerosol particles is used to indicate the spatial distribution of aerosol particles under the effects of entrainment and sedimentation in the containment flow field; a removal efficiency determination module, configured to determine the removal efficiency of the spray droplets on the aerosol particles based on the state of the spray droplets and the spatial distribution of the aerosol particles; The removal result module is used to determine the removal result of the aerosol by the spray droplets based on the removal efficiency of the aerosol particles.

14. A non-volatile storage medium, characterized in that: The non-volatile storage medium stores a plurality of instructions, and the instructions are suitable for being loaded by a processor and executed by the method for determining an aerosol removal result by coupling multiple influencing factors as described in any one of claims 1 to 12.

15. A computer program product comprising a computer program, characterized in that When the computer program is executed by a processor, the steps of the method for determining an aerosol removal result by coupling multiple influencing factors as described in any one of claims 1 to 12 are implemented.

Images