Method and device for calculating pump-induced pulsating pressure load of reactor and related product

Abstract

The invention relates to the technical field of nuclear engineering, and discloses a reactor pump-induced pulsating pressure load calculation method and device and related products. The method comprises the following steps: establishing a three-dimensional model of a reactor component to be calculated; simulation is carried out under the condition that main pump inlet and outlet pulsating pressure serves as a simulation excitation source, and frequency response data of all monitoring points are obtained; generating target frequency response data for any target reactor component based on the frequency response data of the associated monitoring point; determining target pressure loads respectively corresponding to a plurality of characteristic frequencies in the actual excitation source based on the target frequency response data; for any target frequency, calculating a corresponding target pressure transfer coefficient according to the amplitude and the target pressure load corresponding to the target frequency; and calculating the pulsating pressure load of the target reactor component based on the amplitude corresponding to each characteristic frequency and each target pressure transfer coefficient. According to the method, a refined reactor assembly three-dimensional model can be constructed, and on this basis, the accuracy of the pump-induced pulsating pressure load obtained through calculation is remarkably improved.

Description

Technical Field

[0001] The present disclosure relates to the field of nuclear engineering technology, and in particular to a method and device for calculating reactor pump-induced pulsation pressure load, and related products. Background Art

[0002] The nuclear-grade main pump, also known as the reactor coolant pump, is an important device in a nuclear power plant that drives the coolant flow and maintains the integrity of the pressure boundary. It plays a vital role in the safe operation of the nuclear power plant. In recent years, the industry has conducted a lot of research on the causes, characteristics, and spectrum information of the main pump pressure pulsation, but has paid less attention to the propagation mode of the main pump pressure pulsation in the reactor. This neglect may bring some serious consequences. For example, when the frequency of the pulsating pressure load generated by the main pump matches the modal frequency of the reactor components, that is, when resonance occurs, the resonance will cause the vibration amplitude of the reactor components to increase sharply, which will accelerate the fatigue damage of the components in the long run; in addition, within a specific reactor coolant volume, the main pump pulsating pressure load may induce standing waves and produce acoustic resonance, thereby aggravating the fatigue damage of the reactor components; in addition, if the pulsating pressure amplitude on the inside and outside of the reactor components is greatly different, the component will be subjected to a large pressure pulsation load, which will gradually cause fatigue damage, leading to crack initiation and rapid propagation. According to RG1.20 (a regulatory guide issued by the U.S. Nuclear Regulatory Commission, which provides guidance on the dynamic analysis of nuclear power plant safety systems), the reactor design should carry out an evaluation and analysis of pump-induced pulsating pressure loads. Therefore, the study of the calculation and evaluation of reactor pump-induced pulsating pressure loads is of great significance.

[0003] In the related technology, the foreign dedicated software ACSTIC2 (a computer code specifically used to predict and analyze pump-induced acoustic pressure in fluid handling systems) can be used to carry out in-pile acoustic analysis, that is, to analyze the transmission mode of the main pump pressure pulsation in the reactor and calculate the pulsating pressure load acting on the reactor components. However, the software ACSTIC2 can only perform a simple two-dimensional modeling of the reactor coolant system, and based on this, calculate and analyze the pulsating pressure load on the components. The final calculation result is not accurate enough and differs greatly from the actual pulsating pressure load. Therefore, it is urgent to solve this technical problem. Summary of the invention

[0004] In view of the above situation, the embodiments of the present disclosure provide a method, device and related products for calculating the reactor pump-induced pulsating pressure load, aiming to solve the above problem or at least partially solve the above problem.

[0005] In a first aspect, an embodiment of the present disclosure provides a method for calculating a reactor pump-induced pulsation pressure load, the method comprising:

[0006] In response to a model building request initiated by a user, a three-dimensional model of a reactor component to be calculated is established based on a target component in an actual reactor coolant system;

[0007] Under the condition that the pulsating pressure at the inlet and outlet of the main pump is used as the simulation excitation source, the three-dimensional model of the reactor assembly to be calculated is subjected to acoustic simulation calculation to obtain the frequency response data corresponding to each monitoring point in the model, wherein the response data in the frequency response data is the pulsating pressure load borne by the monitoring point;

[0008] For any target reactor component in the reactor assembly model to be calculated, generating target frequency response data based on frequency response data corresponding to an associated monitoring point of the target reactor component;

[0009] Based on the target frequency response data, determining target pressure loads corresponding to a plurality of characteristic frequencies in an actual excitation source;

[0010] For any target frequency among the multiple characteristic frequencies, a target pressure transfer coefficient of the target reactor component at the target frequency is calculated according to the amplitude corresponding to the target frequency and the target pressure load;

[0011] Based on the amplitude corresponding to each of the characteristic frequencies and each of the target pressure transfer coefficients, the pulsating pressure load borne by the target reactor component is calculated.

[0012] In a second aspect, an embodiment of the present disclosure further provides a device for calculating a reactor pump-induced pulsation pressure load, the device comprising:

[0013] A model building module, for responding to a model building request initiated by a user, and building a three-dimensional model of a reactor component to be calculated based on a target component in an actual reactor coolant system;

[0014] A simulation module is used to perform acoustic simulation calculation on the three-dimensional model of the reactor assembly to be calculated under the condition that the pulsating pressure at the inlet and outlet of the main pump is used as a simulation excitation source, and obtain frequency response data corresponding to each monitoring point in the model, wherein the response data in the frequency response data is the pulsating pressure load borne by the monitoring point;

[0015] A frequency response data generating module, for generating target frequency response data for any target reactor component in the reactor assembly model to be calculated based on the frequency response data corresponding to the associated monitoring point of the target reactor component;

[0016] a transfer coefficient calculation module, for determining target pressure loads corresponding to a plurality of characteristic frequencies in an actual excitation source respectively based on the target frequency response data; and for any target frequency among the plurality of characteristic frequencies, calculating a target pressure transfer coefficient of the target reactor component at the target frequency according to the amplitude corresponding to the target frequency and the target pressure load;

[0017] The load calculation module is used to calculate the pulsating pressure load borne by the target reactor component based on the amplitude corresponding to each characteristic frequency and each target pressure transfer coefficient.

[0018] In a third aspect, an embodiment of the present disclosure further provides an electronic device, comprising: a processor; and a memory arranged to store computer executable instructions, which, when executed, cause the processor to execute the steps of the above-mentioned method for calculating the reactor pump-induced pulsating pressure load.

[0019] In a fourth aspect, an embodiment of the present disclosure further provides a computer-readable storage medium, which stores one or more programs. When the one or more programs are executed by an electronic device including multiple application programs, the electronic device executes the steps of the above-mentioned reactor pump-induced pulsating pressure load calculation method.

[0020] By means of the above technical scheme, the reactor pump-induced pulsating pressure load calculation method, device and related products provided in the embodiment of the present disclosure can construct a refined three-dimensional model of the reactor component based on the target component in the actual reactor coolant system, compared with the two-dimensional load calculation method in the prior art; under the condition that the pulsating pressure at the inlet and outlet of the main pump is used as the simulation excitation source, the three-dimensional model is subjected to acoustic simulation calculation, and the response data corresponding to each monitoring point in the model is obtained as the frequency response data of the pulsating pressure load borne by the monitoring point, and on this basis, the target pressure transfer coefficient corresponding to the target reactor component under multiple characteristic frequencies of the actual excitation source is calculated; finally, according to the amplitude corresponding to each characteristic frequency and the corresponding target pressure transfer coefficient, the pulsating pressure load borne by the target reactor component is obtained by linear superposition. The scheme provided in this embodiment can significantly improve the accuracy of the pump-induced pulsating pressure load. In addition, the calculated pulsating pressure load can be used as the input of the subsequent stress fatigue analysis of the reactor component.

[0021] The above description is only an overview of the technical solution of the present disclosure. In order to more clearly understand the technical means of the present disclosure, it can be implemented according to the contents of the specification. In order to make the above and other purposes, features and advantages of the present disclosure more obvious and easy to understand, the specific implementation methods of the present disclosure are listed below. BRIEF DESCRIPTION OF THE DRAWINGS

[0022] The drawings described herein are used to provide a further understanding of the present disclosure and constitute a part of the present disclosure. The illustrative embodiments of the present disclosure and their descriptions are used to explain the present disclosure and do not constitute an improper limitation on the present disclosure. In the drawings:

[0023] Figure 1 A schematic flow chart of a method for calculating a reactor pump-induced pulsation pressure load provided in an embodiment of the present disclosure is shown;

[0024] Figure 2A An acoustic pressure response diagram of a monitoring point in the lower pile internals region provided by an embodiment of the present disclosure is shown;

[0025] Figure 2B It shows a sound pressure response diagram of a monitoring point in the descending section area provided by an embodiment of the present disclosure;

[0026] Figure 2C An acoustic pressure response diagram of a monitoring point in a fuel assembly area provided by an embodiment of the present disclosure is shown;

[0027] Figure 3 A schematic diagram of a perforated plate impedance model provided by an embodiment of the present disclosure is shown;

[0028] Figure 4 A schematic diagram showing the pulsating pressure load on the reactor components provided by the embodiment of the present disclosure is shown;

[0029] Figure 5 A schematic diagram of the structure of a reactor pump-induced pulsation pressure load calculation device provided in an embodiment of the present disclosure is shown;

[0030] Figure 6 A schematic structural diagram of an electronic device provided by an embodiment of the present disclosure is shown. DETAILED DESCRIPTION

[0031] In order to make the purpose, technical solutions and advantages of the present disclosure clearer, the technical solutions of the present disclosure will be clearly and completely described below in combination with the specific embodiments of the present disclosure and the corresponding drawings. Obviously, the described embodiments are only part of the embodiments of the present disclosure, not all of the embodiments. Based on the embodiments in the present disclosure, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present disclosure.

[0032] It should be noted that similar reference numerals and letters denote similar items in the following drawings, and therefore, once an item is defined in one drawing, further definition and explanation thereof is not required in subsequent drawings.

[0033] It should be noted that the terms "first", "second", etc. in the specification and claims of the present disclosure and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that such usage is interchangeable where appropriate, so that the embodiments of the present disclosure described herein can be implemented in an order other than those illustrated or described herein. In addition, the term "including" and its variants are to be interpreted as open-ended terms meaning "including but not limited to".

[0034] First, the professional terms involved in the embodiments of the present disclosure are explained:

[0035] 1) Reactor pump pulsating pressure load refers to the dynamic load generated by the periodic pressure fluctuations caused by the operation of the reactor coolant main pump acting on the internal components of the reactor coolant system.

[0036] As introduced above, the foreign dedicated software ACSTIC2 (a computer code specifically used to predict and analyze pump-induced acoustic pressure in fluid handling systems) can be used to carry out in-pile acoustic analysis, that is, to analyze the transmission mode of the main pump pressure pulsation in the reactor and calculate the pulsating pressure load acting on the reactor components. However, the software ACSTIC2 has defects such as being able to only perform simple two-dimensional modeling of the reactor coolant system and being unable to define complex acoustic boundaries. Therefore, the software is used to calculate and analyze the pulsating pressure load on the components, and the final calculation result is not accurate enough and differs greatly from the actual pulsating pressure load. Based on this, the present invention proposes a method, device and related products for calculating the pulsating pressure load caused by a reactor pump. The present disclosure is described in detail below through specific embodiments.

[0037] To facilitate understanding of this embodiment, a method for calculating a reactor pump-induced pulsating pressure load disclosed in an embodiment of the present disclosure is first introduced in detail. The execution subject of the method for calculating a reactor pump-induced pulsating pressure load provided in an embodiment of the present disclosure is generally a computer device with certain computing capabilities, and the computer device includes, for example: a terminal device or a server or other processing device, and the terminal device can be a user equipment (UE), a mobile device, a user terminal, or a terminal. In some possible implementations, the method for calculating a reactor pump-induced pulsating pressure load can be implemented by a processor calling a computer-readable instruction stored in a memory.

[0038] Figure 1 The schematic diagram of the process flow of the method for calculating the reactor pump pulsation pressure load provided by the embodiment of the present disclosure is shown. Figure 1 It can be seen that the embodiment of the present disclosure at least includes steps S101-S106:

[0039] S101: In response to a model building request initiated by a user, a three-dimensional model of a reactor component to be calculated is established based on a target component in an actual reactor coolant system.

[0040] This embodiment is implemented based on finite element software. Here, the finite element software may be, for example, OpenMC, ANSYS, etc., which is not limited in this embodiment. The target component may be selected according to actual needs. In a possible implementation, the target component includes at least one of the following: a reactor and a steam generator. During implementation, for example, if the target component includes a reactor, a three-dimensional model of the reactor component to be calculated may be established based on the reactor pressure vessel, the core and the internal components of the reactor, and the main pipeline. For another example, if the target component includes a steam generator, a three-dimensional model of the reactor component to be calculated may be established based on the steam generator and the main pipeline.

[0041] S102: Under the condition that the pulsating pressure at the inlet and outlet of the main pump is used as a simulation excitation source, an acoustic simulation calculation is performed on the three-dimensional model of the reactor assembly to be calculated to obtain frequency response data corresponding to each monitoring point in the model, wherein the response data corresponding to the frequency response data is the pulsating pressure load borne by the monitoring point.

[0042] This embodiment uses a simulated excitation source to replace the pulsating pressure of the main pump inlet and outlet. For example, if the target component includes a reactor, and there are 3 groups of main pumps in the actual reactor coolant system, 6 simulated excitation sources can be set as the pulsating pressure of each group of main pumps in the system. During implementation, the parameters of each simulated excitation source can be set according to the spectrum data provided by the main pump manufacturer. Specifically, in order to obtain the frequency response data of each monitoring point, for any simulated excitation source, the type, amplitude, phase, starting frequency, end frequency, frequency step and other parameters of the excitation source can be set. Exemplarily, the amplitude of the excitation source can be set to 1000, the phase to 0, the starting frequency to 0Hz, the end frequency to 500Hz, and the frequency step to 0.5Hz. In addition, the amplitude and phase of each simulated excitation source can be set according to actual needs. For the phase, generally speaking, when all the main pumps are in the same phase, the pump-induced pulsating load is the largest. Therefore, during implementation, this working condition can be used for conservative calculation. Finite element software can be used to set multiple monitoring points in the coolant of the flow channel.

[0043] During implementation, the three-dimensional model of the reactor component to be calculated can be meshed in response to a solver configuration request initiated by a user, and a solver can be configured to perform acoustic simulation calculations on the model based on the solver under the condition that the pulsating pressure at the inlet and outlet of the main pump is used as a simulation excitation source to obtain frequency response data corresponding to each monitoring point in the model.

[0044] S103: For any target reactor component in the reactor assembly model to be calculated, generate target frequency response data based on frequency response data corresponding to an associated monitoring point of the target reactor component.

[0045] It is understandable that the associated monitoring point of the target reactor component here is a monitoring point near the target reactor component. For example, the spatial distance between each monitoring point and the center of gravity of the target reactor component can be calculated, and the monitoring point with the smallest spatial distance can be selected as the associated monitoring point of the target reactor component.

[0046] Exemplarily, the associated monitoring point of the target reactor component is the aforementioned monitoring point A, and the frequency response data corresponding to the monitoring point A can be used as the target frequency response data corresponding to the target reactor component.

[0047] S104: Determine target pressure loads corresponding to a plurality of characteristic frequencies in an actual excitation source based on the target frequency response data.

[0048] Here, the characteristic frequencies in the actual excitation source may include shaft frequency, secondary shaft frequency, primary blade frequency, and secondary blade frequency. In implementation, a frequency response curve may be obtained based on discrete target frequency response data. For example, assuming that the number of monitoring points is 240, Figure 2A The acoustic pressure response diagram of a monitoring point in the lower pile internal component area provided by an embodiment of the present disclosure is shown. Figure 2B A sound pressure response diagram of a monitoring point in the descending section area provided by an embodiment of the present disclosure is shown. Figure 2C The acoustic pressure response diagram of a monitoring point in the fuel assembly area provided by the embodiment of the present disclosure is shown. Then, for any characteristic frequency, a frequency sweep can be performed within a certain range before and after the characteristic frequency, and the maximum amplitude is taken as the corresponding target pressure load.

[0049] S105: For any target frequency among the multiple characteristic frequencies, calculate a target pressure transfer coefficient of the target reactor component at the target frequency according to the amplitude and target pressure load corresponding to the target frequency.

[0050] For example, if the target frequency is the shaft frequency, the ratio of the target pressure load and the amplitude corresponding to the shaft frequency can be used as the target pressure transfer coefficient corresponding to the target reactor component at the shaft frequency. Similarly, the target pressure transfer coefficient corresponding to the target reactor component at other characteristic frequencies can be calculated.

[0051] S106: Based on the amplitude corresponding to each characteristic frequency and each target pressure transfer coefficient, the pulsating pressure load borne by the target reactor component is calculated.

[0052] In specific implementation, for example, for any characteristic frequency M1 in the actual excitation source, the amplitude corresponding to the characteristic frequency M1 can be multiplied by the corresponding target pressure transfer coefficient to obtain the product N1. Similarly, the products N2, N3 and N4 corresponding to other characteristic frequencies can be obtained. The products N1, N2, N3 and N4 are added together to obtain the pulsating pressure load borne by the target reactor components.

[0053] from Figure 1 It can be seen from the method shown that, compared with the two-dimensional load calculation method in the prior art, the embodiment of the present disclosure can construct a refined three-dimensional model of the reactor component based on the target component in the actual reactor coolant system; under the condition that the pulsating pressure at the inlet and outlet of the main pump is used as the simulation excitation source, the three-dimensional model is subjected to acoustic simulation calculation, and the response data corresponding to each monitoring point in the model is obtained as the frequency response data of the pulsating pressure load borne by the monitoring point, and on this basis, the target pressure transfer coefficient corresponding to the target reactor component under multiple characteristic frequencies of the actual excitation source is calculated; finally, according to the amplitude corresponding to each characteristic frequency and the corresponding target pressure transfer coefficient, the pulsating pressure load borne by the target reactor component is obtained by linear superposition. The solution provided in this embodiment can significantly improve the accuracy of the pump-induced pulsating pressure load. In addition, the calculated pulsating pressure load can be used as the input for subsequent stress fatigue analysis of reactor components.

[0054] In a possible implementation, the method further includes: in response to a geometric model establishment request initiated by a user, establishing a three-dimensional geometric model of the reactor component to be calculated based on the target component; in response to a parameter setting request initiated by a user, using the acquired acoustic parameter values ​​to set the acoustic parameters of the coolant in the three-dimensional geometric model of the reactor component to be calculated, the acoustic parameters including density and sound velocity; in response to a damping setting request initiated by a user, using the acquired attenuation coefficient to set the damping characteristics of the coolant in each flow channel in the three-dimensional geometric model of the reactor component to be calculated; in response to a boundary condition setting request initiated by a user, using the acquired boundary characteristic data to set the boundary conditions of the three-dimensional geometric model of the reactor component to be calculated, and obtaining the three-dimensional model of the reactor component to be calculated.

[0055] In this embodiment, the acoustic parameters of the coolant include but are not limited to density, sound velocity, etc. During implementation, the acoustic parameter values ​​of the coolant can be determined according to the actual hot operating conditions of the coolant. Specifically, they can be set according to the cold section and hot section of the coolant, respectively, and can also be set according to the average temperature and pressure of the reactor during operation.

[0056] In response to a damping setting request initiated by a user, the damping characteristic of the coolant in each flow channel in the three-dimensional geometric model of the reactor assembly to be calculated can be set using the acquired attenuation coefficient. When implementing the attenuation coefficient, for example, it can be set by a technician based on engineering experience.

[0057] In a possible implementation manner, the attenuation coefficient is generated according to the following method: using a preset model selection rule, according to the flow channel property data of each flow channel in the three-dimensional geometric model of the reactor assembly to be calculated, determining a target damping model, wherein the target damping model is one of a viscous damping model, a mass damping model, a stiffness damping model, a Rayleigh damping model, and a hysteresis damping model; for any target flow channel in the three-dimensional geometric model of the reactor assembly to be calculated, based on the target damping model, calculating the damping characteristic data corresponding to the target flow channel, and calculating the corresponding attenuation coefficient according to the damping characteristic data;

[0058] The step of calculating damping characteristic data corresponding to the target flow channel based on the target damping model, and calculating a corresponding attenuation coefficient according to the damping characteristic data, includes:

[0059] If the target damping model is the viscous damping model, it is determined whether the cross section of the target flow channel is a single flow channel. If the cross section of the target flow channel is a single flow channel, the damping characteristic data is calculated according to the following formula (1); if the cross section of the target flow channel is a multi-flow channel, the damping characteristic data is calculated according to the following formula (2); and the attenuation coefficient is calculated according to the following formula (3):

[0060]

[0061] If the target damping model is the mass damping model, the damping characteristic data is calculated according to the following formula (4), and the attenuation coefficient is calculated according to the following formula (5):

[0062]

[0063] If the target damping model is the stiffness damping model, the damping characteristic data is calculated according to the following formula (6), and the attenuation coefficient is calculated according to the above formula (5):

[0064]

[0065] If the target damping model is the Rayleigh damping model, the damping characteristic data is calculated according to the following formula (7), and the attenuation coefficient is calculated according to the above formula (5):

[0066]

[0067] If the target damping model is the hysteresis damping model, the damping characteristic data is determined according to the following formula (8), and the attenuation coefficient is calculated according to the above formula (5):

[0068] ζ=const(8)

[0069] Where R is the damping parameter, Δp is the pressure drop of the coolant flow, and q v is the coolant mass flow rate, ρ is the coolant density, l is the flow channel length, d is the flow channel equivalent diameter, k i is the resistance coefficient of the ith flow channel, A i is the area of ​​the i-th flow channel, N is the number of flow channels, N P is the preset attenuation coefficient, c represents the speed of sound, ζ represents the damping ratio, α and β represent the preset coefficients, ω represents the frequency, and const represents a constant.

[0070] After research, it was found that the momentum equation corresponding to the acoustic equation can be written as follows:

[0071]

[0072] Among them, ρ is the fluid density, u is the fluid velocity, and R is the transmission loss in the sound wave transmission, which determines the damping characteristics of the reactor sound field. According to the above definition, the acoustic characteristic equation (frequency response function) in the frequency domain can be obtained as follows:

[0073] ω 2 -iω·Rc 2 γ 2 =0 (10)

[0074] Wherein, γ represents the wave number of the sound wave, ω represents the frequency, and i represents the imaginary unit. It can be seen that the equation shown in formula (10) is similar to the frequency response function structure of the single-degree-of-freedom oscillator system in structural dynamics. Therefore, this embodiment refers to the method of structural mechanics and proposes an acoustic damping system, including a viscous damping model, a mass damping model, a stiffness damping model, a Rayleigh damping model, and a hysteresis damping model. In the above five damping models, the relationship between the damping parameter R, the damping ratio ζ, the speed of sound c, and the frequency ω is as shown in the above formulas (1), (2), (4), (6), and (7). In addition, for the viscous damping model, ζ=R / 2ωc. Here, for the preset coefficients α and β in the mass damping model, the stiffness damping model, and the Rayleigh damping model, the characteristic frequency and the corresponding damping ratio can be selected and brought into the formula for calculation.

[0075] During implementation, the preset model selection rules can be used to first determine the target damping model from the viscous damping model, mass damping model, stiffness damping model, Rayleigh damping model, and hysteresis damping model according to the flow channel property data of each flow channel in the three-dimensional geometric model of the reactor assembly to be calculated. Here, the flow channel property data includes but is not limited to: flow channel shape, flow channel size, whether there is a dead water area, flow channel area and other parameters, depending on the preset model selection rules. The preset model selection rules can be set according to actual needs, and the embodiments of the present disclosure are not limited to this.

[0076] In a possible implementation, the flow channel property data includes the flow channel shape and whether there is a dead water area. The preset model selection rule is used to determine the target damping model according to the flow channel property data of each flow channel in the three-dimensional geometric model of the reactor assembly to be calculated, including: based on each flow channel property data, judging whether there are a first flow channel and a second flow channel in the three-dimensional geometric model of the reactor assembly to be calculated; wherein the first flow channel has a dead water area, and the second flow channel includes at least one of the following: an elliptical flow channel, a spiral flow channel and an irregular flow channel; if the first flow channel and the second flow channel do not exist, selecting the viscous damping model as the target damping model; if the first flow channel or the second flow channel exists, selecting one from the mass damping model, the stiffness damping model, the Rayleigh damping model, and the hysteresis damping model as the target damping model.

[0077] This embodiment does not limit the type of the second flow channel, and it can be set according to actual conditions. For example, the second flow channel can also include a porous medium pore flow channel.

[0078] When the damping characteristic data of the first flow channel with a dead water area is calculated using the above formula (1) or (2), the pressure drop data of the coolant flow is missing, resulting in the calculated damping being 0, which is inaccurate; when the damping characteristic data of the second flow channel is calculated, the cross section of the second flow channel is complex and it is also inconvenient to calculate. Therefore, in this embodiment, when the first flow channel and the second flow channel do not exist in the model, the viscous damping model is selected as the target damping model; if the first flow channel or the second flow channel exists, one of the mass damping model, the stiffness damping model, the Rayleigh damping model, and the hysteresis damping model is selected as the target damping model. During implementation, the three-dimensional geometric model of the reactor assembly to be calculated, the mass damping model, the stiffness damping model, the Rayleigh damping model, and the hysteresis damping model can be sent to the terminal held by the technician, and the technician selects the damping model based on engineering experience, and sends the selection result to the execution subject of this embodiment, and the selection result is used as the target damping model.

[0079] For the above formula (10), considering γ as an unknown, the expression of the wave number when considering the damping effect can be obtained under the condition of small damping:

[0080]

[0081] Considering ω as an unknown, solving equation (11) yields:

[0082]

[0083] When the square root of equation (12) is 0, the acoustic damping is critical damping and the damping ratio ζ is 1. At this time:

[0084] R0=2cγ (13)

[0085] Under low damping conditions, γ can be written as follows:

[0086] ω=cγ (14)

[0087] According to the acoustic characteristic equation considering damping, the relationship between the damping ratio, R and frequency can be obtained:

[0088]

[0089] Substituting formula (15) into formula (11), we can get:

[0090]

[0091] In finite element software, the attenuation coefficient is used to represent the attenuation of pressure waves, where N P represents the attenuation coefficient defined in the finite element software, and its equation is:

[0092]

[0093] It can be deduced that:

[0094]

[0095] Therefore, for different damping types, the relationship between the attenuation coefficient and the damping ratio or R is as shown in the above equations (3) and (5), that is, for the viscous damping model:

[0096]

[0097] Mass Damping Model:

[0098]

[0099] Stiffness damping model:

[0100]

[0101] Rayleigh damping model:

[0102]

[0103] Hysteretic Damping Model:

[0104]

[0105] Table 1 below shows the R values ​​of the main regions provided in this embodiment:

[0106] Table 1 Damping characteristics of main areas in the reactor (viscous damping)

[0107] Reactor area R / s-1 Main pump outlet-pressure vessel 44.5 Pressure vessel-main pump inlet 10.7 Fuel assembly flow channel 11.8 Descending segment 13.0 Guide cylinder area 11.3 Upper internals 8.09

[0108] This embodiment is analogous to structural dynamics and creatively proposes a damping model system, which includes 5 models (viscous damping model, mass damping model, stiffness damping model, Rayleigh damping model, hysteresis damping model). During implementation, a more practical damping model can be selected for the model based on the attribute data of each flow channel in the model, thereby significantly improving the accuracy of the coolant damping characteristic setting of each flow channel in the model. At the same time, compared with the damping calculation method in the prior art that can only control one-dimensional sound waves, the damping calculation formula provided in this embodiment can control the three-dimensional sound field, and the calculation result is more accurate.

[0109] In response to a boundary condition setting request initiated by a user, the boundary condition of the three-dimensional geometric model of the reactor assembly to be calculated can be set using the acquired boundary characteristic data to obtain the three-dimensional model of the reactor assembly to be calculated. Here, the boundary conditions include but are not limited to the boundary conditions such as the perforated plate boundary and the boundary of different media. The boundary characteristic data is used to define the boundary conditions. During implementation, for the perforated plate boundary, the boundary characteristic data may include the impedance rate; for the boundary of different media, the boundary characteristic data may include the normal impedance of the interface.

[0110] Specifically, in a possible implementation manner, the boundary condition includes a perforated plate boundary condition, and the boundary characteristic data includes an impedance ratio; the impedance ratio is determined according to the following formula:

[0111]

[0112] Where η = a / 2, Z trans represents the transmission impedance, ρ0 represents the fluid density, c0 represents the fluid sound velocity, μ represents the dynamic viscosity, k represents the shear (viscosity) wave number, and t p represents the plate thickness, a represents the hole radius, η represents the hole end correction factor, σ represents the porosity, R f Indicates flow resistance.

[0113] This embodiment uses the impedance characteristic boundary to simulate the boundary conditions of the perforated plate structure. During implementation, the above formula (19) can be used to set internal impedance boundary conditions for the perforated plates in the model, such as the flow distribution cover, the lower core plate, the upper core plate, the grid plate, and the primary side of the compression cylinder. Figure 3 The schematic diagram of the impedance model of the perforated plate provided by the embodiment of the present disclosure is shown, wherein Lx and Ly represent the distance between the centers of adjacent holes in the horizontal direction and the distance between the centers of adjacent holes in the vertical direction, respectively.

[0114] During implementation, the flow channel with complex structure can also be simplified. For example, if the target component includes a reactor, since the flow channel of the fuel part is narrow and long, only its plane wave propagation mode can be considered, and other high-order sound waves are not considered.

[0115] This embodiment uses finite element software to construct a three-dimensional geometric model of the reactor assembly to be calculated, and sets the acoustic parameter values, coolant damping characteristics, and complex acoustic boundary conditions of the three-dimensional geometric model, and finally obtains the three-dimensional model of the reactor assembly to be calculated, thereby making the constructed three-dimensional model of the reactor assembly to be calculated more accurate and realistic, thereby ensuring the reliability and accuracy of the pressure load obtained by subsequent model-based calculations.

[0116] In a possible implementation manner, the associated monitoring points include inner associated monitoring points and outer associated monitoring points; and for any target reactor component in the reactor assembly model to be calculated, generating target frequency response data based on the frequency response data corresponding to the associated monitoring points of the target reactor component includes: acquiring a plurality of first pulsating pressure loads in the frequency response data of the inner associated monitoring points and a plurality of second pulsating pressure loads in the frequency response data of the outer associated monitoring points; calculating each target load difference according to each of the first pulsating pressure loads and the corresponding second pulsating pressure load; and generating the target frequency response data according to each target load difference.

[0117] Figure 4 A schematic diagram of the pulsating pressure load on the reactor components provided by the embodiment of the present disclosure is shown. Figure 4 This embodiment is described as an example.

[0118] For the target reactor component of the hanging basket assembly (part), in order to obtain its corresponding target frequency response data, we can first obtain multiple first pulsating pressure loads C11 (pressure pulsation inside the hanging basket) in the frequency response data of the inner associated monitoring point, and obtain multiple second pulsating pressure loads C21 (pressure pulsation outside the hanging basket) in the frequency response data of the outer associated monitoring point, and subtract the first pulsating pressure load C11 from the corresponding second pulsating pressure load C21 (the frequencies corresponding to C11 and C21 are equal) to obtain the target load differences c1, c2, c3, ..., c mWhere m is the frequency number of the excitation source. That is, the target load difference of the hanging basket assembly can be calculated according to the following formula:

[0119] P dif =P inner -P outer =(P inner-real -P outer-real )+(P inner-img -P outer-img )i (20)

[0121] Among them, P inner Represents the first pulsating pressure load, P outer Represents the second pulsating pressure load, P inner-real represents the real part of the first pulsating pressure load, P inner-img Represents the imaginary part of the first pulsating pressure load, P outer-real , represents the real part of the second pulsating pressure load, P outer-img Represents the imaginary part of the second pulsating pressure load.

[0122] Finally, according to the target load difference c1, c2, c3, ..., c m and their corresponding frequencies to generate target frequency response data corresponding to the target reactor components.

[0123] This embodiment uses the inner and outer associated monitoring points of the target reactor component to calculate the pressure load on the component in the target frequency response data. This method can make the pressure load more accurate, thereby making the target frequency response data more accurate, and further significantly improve the accuracy and reliability of the pulsating pressure transfer coefficient subsequently calculated based on the target frequency response data, as well as the pulsating pressure load on the target reactor component.

[0124] The embodiment of the present disclosure also provides a method for calculating the reactor pump-induced pulsation pressure load, comprising the following steps:

[0125] Step S1: In response to a geometric model establishment request initiated by a user, a three-dimensional geometric model of a reactor component to be calculated is established based on a reactor in an actual reactor coolant system.

[0126] Step S2: In response to a parameter setting request initiated by the user, the acoustic parameters of the coolant in the three-dimensional geometric model of the reactor assembly to be calculated are set using the acquired acoustic parameter values, where the acoustic parameters include density and sound velocity.

[0127] Step S3: In response to the damping setting request initiated by the user, the damping characteristics of the coolant in each flow channel in the three-dimensional geometric model of the reactor assembly to be calculated are set using the acquired attenuation coefficient.

[0128] Preferably, the attenuation coefficient is generated according to the following method: based on the attribute data of each flow channel, determine whether there are the first flow channel and the second flow channel in the three-dimensional geometric model of the reactor assembly to be calculated; wherein, there is a dead water area in the first flow channel, and the second flow channel includes at least one of the following: an elliptical flow channel, a spiral flow channel and an irregular flow channel; if the first flow channel and the second flow channel do not exist, select a viscous damping model as the target damping model; if the first flow channel or the second flow channel exists, select one from the mass damping model, the stiffness damping model, the Rayleigh damping model, and the hysteresis damping model as the target damping model; for any target flow channel in the three-dimensional geometric model of the reactor assembly to be calculated, based on the target damping model, calculate the damping characteristic data corresponding to the target flow channel, and calculate the corresponding attenuation coefficient according to the damping characteristic data; based on the target damping model, calculate the damping characteristic data corresponding to the target flow channel, and calculate the corresponding attenuation coefficient according to the damping characteristic data, including: if the target damping If the model is a viscous damping model, determine whether the cross section of the target flow channel is a single flow channel. If the cross section of the target flow channel is a single flow channel, calculate the damping characteristic data according to the above formula (1); if the cross section of the target flow channel is a multi-flow channel, calculate the damping characteristic data according to the above formula (2); and calculate the attenuation coefficient according to the above formula (3): if the target damping model is a mass damping model, calculate the damping characteristic data according to the above formula (4), and calculate the attenuation coefficient according to the above formula (5): if the target damping model is a stiffness damping model, calculate the damping characteristic data according to the above formula (6), and calculate the attenuation coefficient according to the above formula (5): if the target damping model is a Rayleigh damping model, calculate the damping characteristic data according to the above formula (7), and calculate the attenuation coefficient according to the above formula (5): if the target damping model is a hysteresis damping model, determine the damping characteristic data according to the above formula (8), and calculate the attenuation coefficient according to the above formula (5).

[0129] Step S4: In response to the perforated plate boundary condition setting request initiated by the user, the boundary conditions of the three-dimensional geometric model of the reactor assembly to be calculated are set using the acquired impedance ratio to obtain the three-dimensional model of the reactor assembly to be calculated. The impedance ratio is determined according to the above formula (19).

[0130] Step S5: Under the condition that the pulsating pressure at the inlet and outlet of the main pump is used as the simulation excitation source, acoustic simulation calculations are performed on the three-dimensional model of the reactor component to be calculated, and frequency response data corresponding to each monitoring point in the model are obtained. The response data in the frequency response data is the pulsating pressure load borne by the monitoring point.

[0131] Step S6: acquiring a plurality of first pulsating pressure loads in the frequency response data of the inner associated monitoring point and a plurality of second pulsating pressure loads in the frequency response data of the outer associated monitoring point.

[0132] Step S7: Calculate each target load difference according to each first pulsating pressure load and the corresponding second pulsating pressure load.

[0133] Step S8: Generate target frequency response data according to each target load difference.

[0134] Step S9: Based on the target frequency response data, determine the target pressure loads corresponding to the multiple characteristic frequencies in the actual excitation source.

[0135] Step S10: For any target frequency among the multiple characteristic frequencies, a target pressure transfer coefficient of a target reactor component at the target frequency is calculated according to the amplitude corresponding to the target frequency and the target pressure load.

[0136] Step S11: based on the amplitude corresponding to each characteristic frequency and each target pressure transfer coefficient, the pulsating pressure load borne by the target reactor component is calculated.

[0137] Those skilled in the art will appreciate that, in the above method of the specific implementation manner, the writing order of the steps does not mean a strict execution order, but does not constitute any limitation on the implementation process. The specific execution order of the steps should be determined by their functions and possible internal logic.

[0138] It should be noted that, in practical applications, all possible implementations described above can be arbitrarily combined to form possible embodiments of the present disclosure, which will not be described one by one here.

[0139] Based on the same concept, the embodiment of the present disclosure also provides a device for calculating the reactor pump-induced pulsation pressure load. Figure 5 FIG. 1 shows a schematic diagram of the structure of a reactor pump pulsation pressure load calculation device provided by an embodiment of the present disclosure, see FIG. Figure 5 As shown, the reactor pump induced pulsation pressure load calculation device 500 provided in the embodiment of the present disclosure includes:

[0140] The model building module 501 is used to respond to a model building request initiated by a user and to build a three-dimensional model of a reactor component to be calculated based on a target component in an actual reactor coolant system;

[0141] The simulation module 502 is used to perform acoustic simulation calculation on the three-dimensional model of the reactor assembly to be calculated under the condition that the pulsating pressure at the inlet and outlet of the main pump is used as the simulation excitation source, and obtain the frequency response data corresponding to each monitoring point in the model, wherein the response data in the frequency response data is the pulsating pressure load borne by the monitoring point;

[0142] The frequency response data generating module 503 is used to generate target frequency response data for any target reactor component in the reactor assembly model to be calculated based on the frequency response data corresponding to the associated monitoring point of the target reactor component;

[0143] The transfer coefficient calculation module 504 is used to determine the target pressure loads corresponding to the multiple characteristic frequencies in the actual excitation source based on the target frequency response data; for any target frequency among the multiple characteristic frequencies, calculate the target pressure transfer coefficient of the target reactor component at the target frequency according to the amplitude corresponding to the target frequency and the target pressure load;

[0144] The load calculation module 505 is used to calculate the pulsating pressure load borne by the target reactor component based on the amplitude corresponding to each of the characteristic frequencies and each of the target pressure transfer coefficients.

[0145] In a possible implementation manner, in the above device, the target component includes at least one of the following: a reactor, a steam generator.

[0146] In a possible implementation, in the above-mentioned device, the model building module 501 is also used for: responding to a geometric model establishment request initiated by a user, establishing a three-dimensional geometric model of the reactor component to be calculated based on the target component; responding to a parameter setting request initiated by a user, using the acquired acoustic parameter values ​​to set the acoustic parameters of the coolant in the three-dimensional geometric model of the reactor component to be calculated, the acoustic parameters including density and sound velocity; responding to a damping setting request initiated by a user, using the acquired attenuation coefficient to set the damping characteristics of the coolant in each flow channel in the three-dimensional geometric model of the reactor component to be calculated; responding to a boundary condition setting request initiated by a user, using the acquired boundary characteristic data to set the boundary conditions of the three-dimensional geometric model of the reactor component to be calculated, and obtaining the three-dimensional model of the reactor component to be calculated.

[0147] In a possible implementation manner, the device further includes a generation module, which is used to: determine a target damping model according to the flow channel property data of each flow channel in the three-dimensional geometric model of the reactor assembly to be calculated using a preset model selection rule, wherein the target damping model is one of a viscous damping model, a mass damping model, a stiffness damping model, a Rayleigh damping model, and a hysteresis damping model; for any target flow channel in the three-dimensional geometric model of the reactor assembly to be calculated, based on the target damping model, calculate the damping characteristic data corresponding to the target flow channel, and calculate the corresponding attenuation coefficient according to the damping characteristic data; the calculation of the damping characteristic data corresponding to the target flow channel based on the target damping model, and calculating the corresponding attenuation coefficient according to the damping characteristic data, includes:

[0148] If the target damping model is the viscous damping model, it is determined whether the cross section of the target flow channel is a single flow channel. If the cross section of the target flow channel is a single flow channel, the damping characteristic data is calculated according to the following formula (1); if the cross section of the target flow channel is a multi-flow channel, the damping characteristic data is calculated according to the following formula (2); and the attenuation coefficient is calculated according to the following formula (3):

[0149]

[0150] If the target damping model is the mass damping model, the damping characteristic data is calculated according to the following formula (4), and the attenuation coefficient is calculated according to the following formula (5):

[0151]

[0152] If the target damping model is the stiffness damping model, the damping characteristic data is calculated according to the following formula (6), and the attenuation coefficient is calculated according to the above formula (5):

[0153]

[0154] If the target damping model is the Rayleigh damping model, the damping characteristic data is calculated according to the following formula (7), and the attenuation coefficient is calculated according to the above formula (5):

[0155]

[0156] If the target damping model is the hysteresis damping model, the damping characteristic data is determined according to the following formula (8), and the attenuation coefficient is calculated according to the above formula (5):

[0157] ζ=const(8)

[0158] Where R is the damping parameter, Δp is the pressure drop of the coolant flow, and q v is the coolant mass flow rate, ρ is the coolant density, l is the flow channel length, d is the flow channel equivalent diameter, k i is the resistance coefficient of the ith flow channel,

[0159] A i is the area of ​​the i-th flow channel, N is the number of flow channels, N P is the preset attenuation coefficient, c represents the speed of sound, ζ represents the damping ratio, α and β represent the preset coefficients, ω represents the frequency, and const represents a constant.

[0160] In a possible implementation, in the above-mentioned device, the flow channel attribute data includes the flow channel shape and whether there is a dead water area; the attenuation coefficient generation module is also used to: based on each flow channel attribute data, determine whether there are a first flow channel and a second flow channel in the three-dimensional geometric model of the reactor assembly to be calculated; wherein, there is a dead water area in the first flow channel, and the second flow channel includes at least one of the following: an elliptical flow channel, a spiral flow channel and an irregular flow channel; if the first flow channel and the second flow channel do not exist, the viscous damping model is selected as the target damping model; if the first flow channel or the second flow channel exists, one is selected from the mass damping model, the stiffness damping model, the Rayleigh damping model, and the hysteresis damping model as the target damping model.

[0161] In a possible implementation manner, in the above device, the boundary condition includes a perforated plate boundary condition, and the boundary characteristic data includes an impedance ratio; the impedance ratio is determined according to the following formula:

[0162]

[0163] Where η = a / 2, Z trans represents the transmission impedance, ρ0 represents the fluid density, c0 represents the fluid sound velocity, μ represents the dynamic viscosity, k represents the shear (viscosity) wave number, and t p represents the plate thickness, a represents the hole radius, η represents the hole end correction factor, σ represents the porosity, R f Indicates flow resistance.

[0164] In a possible implementation manner, in the above-mentioned device, the associated monitoring points include inner associated monitoring points and outer associated monitoring points; the frequency response data generating module 503 is used to: obtain multiple first pulsating pressure loads in the frequency response data of the inner associated monitoring points and multiple second pulsating pressure loads in the frequency response data of the outer associated monitoring points; calculate a target load difference according to each of the first pulsating pressure loads and the corresponding second pulsating pressure loads; and generate the target frequency response data according to each target load difference.

[0165] It should be noted that any of the above-mentioned reactor pump-induced pulsating pressure load calculation devices can implement the above-mentioned reactor pump-induced pulsating pressure load calculation method in a one-to-one correspondence, which will not be repeated here.

[0166] Figure 6 FIG. 1 shows a schematic diagram of the structure of an electronic device provided by an embodiment of the present disclosure. Figure 6As shown, at the hardware level, the electronic device includes a processor, and optionally also includes an internal bus, a network interface, and a memory. Among them, the memory may include a memory, such as a high-speed random access memory (RAM), and may also include a non-volatile memory (non-volatile memory), such as at least one disk storage, etc. Of course, the electronic device may also include hardware required for other services.

[0167] The processor, network interface and memory can be interconnected through an internal bus, which can be an ISA (Industry Standard Architecture) bus, a PCI (Peripheral Component Interconnect) bus or an EISA (Extended Industry Standard Architecture) bus. The bus can be divided into an address bus, a data bus, a control bus, etc. For ease of representation, Figure 6 Only one bidirectional arrow is used in the diagram, but this does not mean that there is only one bus or only one type of bus.

[0168] The memory is used to store programs. Specifically, the program may include program code, and the program code includes computer operation instructions. The memory may include internal memory and non-volatile memory, and provides instructions and data to the processor.

[0169] The processor reads the corresponding computer program from the non-volatile memory into the memory and then runs it, forming a reactor pump-induced pulsation pressure load calculation device at the logical level. The processor executes the program stored in the memory and is specifically used to execute the above method.

[0170] The processor may be an integrated circuit chip with signal processing capabilities. In the implementation process, each step of the above method can be completed by an integrated logic circuit of hardware in the processor or an instruction in the form of software. The above processor may be a general-purpose processor, including a central processing unit (CPU), a network processor (NP), etc.; it may also be a digital signal processor (DSP), an application-specific integrated circuit (ASIC), a field-programmable gate array (FPGA) or other programmable logic devices, discrete gates or transistor logic devices, discrete hardware components. The disclosed methods, steps and logic block diagrams in the embodiments of the present disclosure may be implemented or executed. The general-purpose processor may be a microprocessor or the processor may also be any conventional processor, etc. The steps of the method disclosed in conjunction with the embodiments of the present disclosure may be directly embodied as being executed by a hardware decoding processor, or may be executed by a combination of hardware and software modules in a decoding processor. The software module may be located in a mature storage medium in the field such as a random access memory, a flash memory, a read-only memory, a programmable read-only memory or an electrically erasable programmable memory, a register, etc. The storage medium is located in the memory, and the processor reads the information in the memory and completes the steps of the above method in combination with its hardware.

[0171] The electronic device can execute the reactor pump pulsation pressure load calculation method provided by multiple embodiments of the present disclosure, and realize the reactor pump pulsation pressure load calculation device in Figure 5 The functions of the illustrated embodiments will not be described in detail in the embodiments of the present disclosure.

[0172] The embodiments of the present disclosure also provide a computer-readable storage medium, which stores one or more programs, and the one or more programs include instructions. When the instructions are executed by an electronic device including multiple application programs, the electronic device can execute the reactor pump-induced pulsating pressure load calculation method provided by multiple embodiments of the present disclosure.

[0173] Those skilled in the art will appreciate that the embodiments of the present disclosure may be provided as methods, systems, or computer program products. Therefore, the present disclosure may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present disclosure may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0174] The present disclosure is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present disclosure. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0175] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0176] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0177] In a typical configuration, a computing device includes one or more processors (CPU), input / output interfaces, network interfaces, and memory.

[0178] The memory may include non-permanent storage in a computer-readable medium, random access memory (RAM) and / or non-volatile memory in the form of read-only memory (ROM) or flash RAM. The memory is an example of a computer-readable medium.

[0179] Computer readable media include permanent and non-permanent, removable and non-removable media that can be implemented by any method or technology to store information. Information can be computer readable instructions, data structures, program modules or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, compact disk read-only memory (CD-ROM), digital versatile disk (DVD) or other optical storage, magnetic cassettes, magnetic tape, disk storage or other magnetic storage devices or any other non-transmission media that can be used to store information that can be accessed by a computing device. As defined herein, computer readable media does not include temporary computer readable media (transitory media), such as modulated data signals and carrier waves.

[0180] The disclosed embodiments also provide a computer program product, which carries a program code. The program code includes instructions that can be used to execute the steps of the method for calculating the reactor pump-induced pulsating pressure load described in the above method embodiment. For details, please refer to the above method embodiment, which will not be repeated here.

[0181] The computer program product may be implemented in hardware, software or a combination thereof. In one optional embodiment, the computer program product is implemented as a computer storage medium. In another optional embodiment, the computer program product is implemented as a software product, such as a software development kit (SDK).

[0182] It should also be noted that the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, commodity or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, commodity or device. In the absence of further restrictions, an element defined by the sentence "comprises a ..." does not exclude the presence of other identical elements in the process, method, commodity or device including the element.

[0183] It will be appreciated by those skilled in the art that the embodiments of the present disclosure may be provided as methods, systems or computer program products. Therefore, the present disclosure may take the form of a complete hardware embodiment, a complete software embodiment or an embodiment combining software and hardware. Moreover, the present disclosure may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0184] The above are only embodiments of the present disclosure and are not intended to limit the present disclosure. For those skilled in the art, the present disclosure may have various modifications and variations. Any modification, equivalent substitution, improvement, etc. made within the spirit and principle of the present disclosure shall be included in the scope of the claims of the present disclosure.

Claims

1. A method for calculating reactor pump-induced pulsation pressure load, characterized in that: The method comprises: In response to a model building request initiated by a user, a three-dimensional model of a reactor component to be calculated is established based on a target component in an actual reactor coolant system; Under the condition that the pulsating pressure at the inlet and outlet of the main pump is used as the simulation excitation source, the three-dimensional model of the reactor assembly to be calculated is subjected to acoustic simulation calculation to obtain the frequency response data corresponding to each monitoring point in the model, wherein the response data in the frequency response data is the pulsating pressure load borne by the monitoring point; For any target reactor component in the reactor assembly model to be calculated, generating target frequency response data based on frequency response data corresponding to an associated monitoring point of the target reactor component; Based on the target frequency response data, determining target pressure loads corresponding to a plurality of characteristic frequencies in an actual excitation source; For any target frequency among the multiple characteristic frequencies, a target pressure transfer coefficient of the target reactor component at the target frequency is calculated according to the amplitude corresponding to the target frequency and the target pressure load; Based on the amplitude corresponding to each of the characteristic frequencies and each of the target pressure transfer coefficients, the pulsating pressure load borne by the target reactor component is calculated.

2. The method according to claim 1, characterized in that The target component includes at least one of the following: a reactor, a steam generator.

3. The method according to claim 1, characterized in that: The method further comprises: In response to a geometric model establishment request initiated by a user, a three-dimensional geometric model of the reactor component to be calculated is established based on the target component; In response to a parameter setting request initiated by a user, the acoustic parameters of the coolant in the three-dimensional geometric model of the reactor assembly to be calculated are set using the acquired acoustic parameter values, wherein the acoustic parameters include density and sound velocity; In response to a damping setting request initiated by a user, using the acquired attenuation coefficient, setting damping characteristics of the coolant in each flow channel in the three-dimensional geometric model of the reactor assembly to be calculated; In response to a boundary condition setting request initiated by a user, the boundary conditions of the three-dimensional geometric model of the reactor assembly to be calculated are set using the acquired boundary characteristic data to obtain the three-dimensional model of the reactor assembly to be calculated.

4. The method according to claim 3, characterized in that The attenuation coefficient is generated according to the following method: Using a preset model selection rule, according to the flow channel property data of each flow channel in the three-dimensional geometric model of the reactor assembly to be calculated, a target damping model is determined, wherein the target damping model is one of a viscous damping model, a mass damping model, a stiffness damping model, a Rayleigh damping model, and a hysteresis damping model; For any target flow channel in the three-dimensional geometric model of the reactor assembly to be calculated, based on the target damping model, calculating the damping characteristic data corresponding to the target flow channel, and calculating the corresponding attenuation coefficient according to the damping characteristic data; The step of calculating damping characteristic data corresponding to the target flow channel based on the target damping model, and calculating a corresponding attenuation coefficient according to the damping characteristic data, includes: If the target damping model is the viscous damping model, it is determined whether the cross section of the target flow channel is a single flow channel. If the cross section of the target flow channel is a single flow channel, the damping characteristic data is calculated according to the following formula (1); if the cross section of the target flow channel is a multi-flow channel, the damping characteristic data is calculated according to the following formula (2); and the attenuation coefficient is calculated according to the following formula (3): If the target damping model is the mass damping model, the damping characteristic data is calculated according to the following formula (4), and the attenuation coefficient is calculated according to the following formula (5): If the target damping model is the stiffness damping model, the damping characteristic data is calculated according to the following formula (6), and the attenuation coefficient is calculated according to the above formula (5): If the target damping model is the Rayleigh damping model, the damping characteristic data is calculated according to the following formula (7), and the attenuation coefficient is calculated according to the above formula (5): If the target damping model is the hysteresis damping model, the damping characteristic data is determined according to the following formula (8), and the attenuation coefficient is calculated according to the above formula (5): ζ=const(8) Where R is the damping parameter, Δp is the pressure drop of the coolant flow, and q v is the coolant mass flow rate, ρ is the coolant density, l is the flow channel length, d is the flow channel equivalent diameter, k i is the resistance coefficient of the ith flow channel, A i is the area of ​​the i-th flow channel, N is the number of flow channels, N P is the preset attenuation coefficient, c represents the speed of sound, ζ represents the damping ratio, α and β represent the preset coefficients, ω represents the frequency, and const represents a constant.

5. The method according to claim 4, characterized in that The flow channel attribute data includes the flow channel shape and whether there is a dead water area; the use of a preset model selection rule to determine the target damping model according to the flow channel attribute data of each flow channel in the three-dimensional geometric model of the reactor assembly to be calculated includes: Based on the attribute data of each flow channel, determining whether there is a first flow channel and a second flow channel in the three-dimensional geometric model of the reactor assembly to be calculated; wherein the first flow channel has a dead water area, and the second flow channel includes at least one of the following: an elliptical flow channel, a spiral flow channel, and an irregular flow channel; If the first flow channel and the second flow channel do not exist, selecting the viscous damping model as the target damping model; If the first flow channel or the second flow channel exists, one is selected from the mass damping model, the stiffness damping model, the Rayleigh damping model, and the hysteresis damping model as the target damping model.

6. The method according to claim 3, characterized in that The boundary conditions include perforated plate boundary conditions, and the boundary characteristic data include impedance ratio; The resistivity is determined according to the following formula: Where η = a / 2, Z trans represents the transmission impedance, ρ0 represents the fluid density, c0 represents the fluid sound velocity, μ represents the dynamic viscosity, k represents the shear (viscosity) wave number, and t p represents the plate thickness, a represents the hole radius, η represents the hole end correction factor, σ represents the porosity, R f Indicates flow resistance.

7. The method according to any one of claims 1 to 6, characterized in that: The associated monitoring points include inner associated monitoring points and outer associated monitoring points; the generating target frequency response data for any target reactor component in the reactor assembly model to be calculated based on the frequency response data corresponding to the associated monitoring points of the target reactor component includes: Acquire a plurality of first pulsating pressure loads in the frequency response data of the inner associated monitoring point and a plurality of second pulsating pressure loads in the frequency response data of the outer associated monitoring point; Calculating each target load difference according to each of the first pulsating pressure loads and the corresponding second pulsating pressure loads; The target frequency response data is generated according to each of the target load differences.

8. A reactor pump induced pulsation pressure load calculation device, characterized in that: The device comprises: A model building module, for responding to a model building request initiated by a user, and building a three-dimensional model of a reactor component to be calculated based on a target component in an actual reactor coolant system; A simulation module is used to perform acoustic simulation calculation on the three-dimensional model of the reactor assembly to be calculated under the condition that the pulsating pressure at the inlet and outlet of the main pump is used as a simulation excitation source, and obtain frequency response data corresponding to each monitoring point in the model, wherein the response data in the frequency response data is the pulsating pressure load borne by the monitoring point; A frequency response data generating module, for generating target frequency response data for any target reactor component in the reactor assembly model to be calculated based on the frequency response data corresponding to the associated monitoring point of the target reactor component; a transfer coefficient calculation module, for determining target pressure loads corresponding to a plurality of characteristic frequencies in an actual excitation source respectively based on the target frequency response data; and for any target frequency among the plurality of characteristic frequencies, calculating a target pressure transfer coefficient of the target reactor component at the target frequency according to the amplitude corresponding to the target frequency and the target pressure load; The load calculation module is used to calculate the pulsating pressure load borne by the target reactor component based on the amplitude corresponding to each characteristic frequency and each target pressure transfer coefficient.

9. An electronic device, comprising: processor; as well as A memory arranged to store computer executable instructions, wherein the executable instructions, when executed, cause the processor to perform the steps of the method according to any one of claims 1 to 7.

10. A computer-readable storage medium storing one or more programs, characterized in that: When the one or more programs are executed by an electronic device including a plurality of application programs, the electronic device executes the steps of the method as claimed in any one of claims 1 to 7.

Images