Reactor pump-induced pulsation pressure load calculation method, device and related products

Abstract

The present disclosure relates to the field of nuclear engineering technology and discloses a method, device, and related products for calculating reactor pump-induced pulsating pressure loads. The method comprises: establishing a three-dimensional model of the reactor assembly to be calculated; performing simulation using the pulsating pressure at the inlet and outlet of the main pump as the simulation excitation source to obtain frequency response data for each monitoring point; generating target frequency response data for any target reactor component based on the frequency response data of its associated monitoring points; determining target pressure loads corresponding to multiple characteristic frequencies in the actual excitation source based on the target frequency response data; calculating the corresponding target pressure transfer coefficient for any target frequency based on the amplitude and target pressure load corresponding to the target frequency; and calculating the pulsating pressure load of the target reactor component based on the amplitude and target pressure transfer coefficient corresponding to each characteristic frequency. The present disclosure can construct a refined three-dimensional model of the reactor assembly, based on which the accuracy of the calculated pump-induced pulsating pressure load can be significantly 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 pulsating pressure loads, and related products. Background Art

[0002] Nuclear-grade main pumps, also known as reactor coolant pumps, are crucial equipment in nuclear power plants, driving coolant flow and maintaining pressure boundary integrity. They are crucial for the safe operation of nuclear power plants. In recent years, the industry has conducted extensive research on the causes, characteristics, and spectral characteristics of main pump pressure pulsations. However, less attention has been paid to the propagation patterns of main pump pressure pulsations within the reactor. This neglect can have serious consequences. For example, when the frequency of the pulsating pressure load generated by the main pump matches the modal frequency of reactor components, creating resonance, the vibration amplitude of the reactor components increases dramatically, accelerating fatigue damage over the long term. Furthermore, within a specific reactor coolant volume, the main pump pulsating pressure load can induce standing waves, generating acoustic resonance, which exacerbates fatigue damage in reactor components. Furthermore, if the pulsating pressure amplitude differs significantly between the inside and outside of a reactor component, the component will experience significant pressure pulsation loads, leading to gradual fatigue damage and the initiation and rapid propagation of cracks. Furthermore, according to RG1.20 (a regulatory guide published by the U.S. Nuclear Regulatory Commission that provides guidance on the dynamic analysis of nuclear power plant safety systems), reactor designs should conduct an assessment and analysis of pump-induced pulsating pressure loads. Therefore, research on the computational assessment of reactor pump-induced pulsating pressure loads is of great significance.

[0003] In related technologies, ACSTIC2, a specialized software developed abroad (a computer code specifically for predicting and analyzing pump-induced acoustic pressure in fluid handling systems), can be used to conduct in-reactor acoustic analysis. This involves analyzing the transmission pattern of the main pump pressure pulsation within the reactor and calculating the pulsating pressure loads acting on reactor components. However, ACSTIC2 can only perform a simple two-dimensional modeling of the reactor coolant system, based on which it calculates and analyzes the pulsating pressure loads on the components. The final calculated results are inaccurate and differ significantly from the actual pulsating pressure loads. Therefore, there is an urgent need to address this technical issue. 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 pulsating pressure load, the method comprising:

[0006] In response to a model building request initiated by a user, a three-dimensional model of the 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, an acoustic simulation calculation is performed on the three-dimensional model of the reactor assembly to be calculated, and frequency response data corresponding to each monitoring point in the model is obtained. 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] Determining target pressure loads corresponding to a plurality of characteristic frequencies in an actual excitation source based on the target frequency response data;

[0010] 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 an amplitude corresponding to the target frequency and a target pressure load;

[0011] 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.

[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 is used to respond to a model building request initiated by a user and build a three-dimensional model of the reactor component to be calculated based on the target component in the actual reactor coolant system;

[0014] a simulation module for performing acoustic simulation calculations on the three-dimensional model of the reactor assembly to be calculated, using the pulsating pressure at the inlet and outlet of the main pump as a simulation excitation source, to 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, configured 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 points of the target reactor component;

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

[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 leveraging the aforementioned technical solutions, the reactor pump-induced pulsating pressure load calculation method, device, and related products provided in the disclosed embodiments can construct a refined three-dimensional model of the target reactor component in an actual reactor coolant system, compared to the two-dimensional load calculation methods used in the prior art. Using the main pump inlet and outlet pulsating pressure as the simulation excitation source, acoustic simulation calculations are performed on the three-dimensional model to obtain the response data corresponding to each monitoring point in the model, which is the frequency response data of the pulsating pressure load experienced by the monitoring point. Based on this, the target pressure transfer coefficient corresponding to the target reactor component at multiple characteristic frequencies of the actual excitation source is calculated. Finally, based on the amplitude and target pressure transfer coefficient corresponding to each characteristic frequency, the pulsating pressure load experienced by the target reactor component is obtained through linear superposition. The solution provided by this embodiment can significantly improve the accuracy of pump-induced pulsating pressure loads. In addition, the calculated pulsating pressure load can be used as input for 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 in accordance with 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 exemplary embodiments of the present disclosure and their descriptions are used to explain the present disclosure and do not constitute an improper limitation of the present disclosure. In the drawings:

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

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

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

[0026] Figure 2C shows an acoustic pressure response diagram of a monitoring point in a fuel assembly area according to an embodiment of the present disclosure;

[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 loads on reactor components provided by an embodiment of the present disclosure is shown;

[0029] Figure 5 A schematic structural diagram of a reactor pump-induced pulsation pressure load calculation device provided by 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] To make the objectives, technical solutions, and advantages of the present disclosure more clear, the technical solutions of the present disclosure will be clearly and completely described below in conjunction 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 of the present disclosure, all other embodiments obtained by ordinary technicians in this field without making creative efforts 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, it does not need to be further defined or explained in subsequent drawings.

[0033] It should be noted that the terms "first," "second," and the like in the specification and claims of the present disclosure and the accompanying drawings are used to distinguish similar objects and are not necessarily used to describe a specific order or precedence. 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 sequences other than those illustrated or described herein. In addition, the term "including" and its variations 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-induced 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 previously mentioned, 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 the defects of only being able to perform simple two-dimensional modeling of the reactor coolant system and being unable to define complex acoustic boundaries. Therefore, when using this software to 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. Based on this, the present invention proposes a method, device and related products for calculating the pulsating pressure load caused by the reactor pump. The present disclosure is described in detail below through specific embodiments.

[0037] To facilitate understanding of this embodiment, a detailed description of a method for calculating reactor pump-induced pulsating pressure loads disclosed in an embodiment of the present disclosure is first provided. The method for calculating reactor pump-induced pulsating pressure loads provided in this embodiment of the present disclosure is generally executed by a computer device with certain computing capabilities, such as a terminal device, a server, or other processing device. The terminal device may be a user equipment (UE), a mobile device, a user terminal, or a terminal. In some possible implementations, the method for calculating reactor pump-induced pulsating pressure loads may be implemented by a processor invoking computer-readable instructions stored in a memory.

[0038] Figure 1 The flow chart of the method for calculating the reactor pump-induced pulsating 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 includes at least 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 can be, for example, OpenMC, ANSYS, etc., which is not limited to this embodiment. The target component can be selected according to actual needs. In one 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 can be established based on the reactor pressure vessel, the core and the in-core components, 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 can 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. 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 at the inlet and outlet of the main pump. For example, if the target component includes a reactor and there are three groups of main pumps in the actual reactor coolant system, six simulated excitation sources can be set as the pulsating pressure at the inlet and outlet 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, ending frequency, frequency step and other parameters of the excitation source can be set. For example, the amplitude of the excitation source can be set to 1000, the phase to 0, the starting frequency to 0Hz, the ending 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. Regarding the phase, generally speaking, when all main pumps have the same phase, the pump-induced pulsating load is the largest. Therefore, during implementation, this working condition can be used for conservative calculations. Finite element software can be used to set multiple monitoring points in the coolant of the flow channel.

[0043] During implementation, in response to a solver configuration request initiated by a user, the three-dimensional model of the reactor component to be calculated can be meshed, 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 the 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 understood that the associated monitoring points of the target reactor component here are monitoring points 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. During 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 An acoustic pressure response diagram of a monitoring point in the lower stack internals region 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 value 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 to 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 coefficients corresponding to the target reactor components 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 component.

[0053] from Figure 1 As can be seen from the method shown, 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. Based on this, the target pressure transfer coefficient corresponding to the target reactor component at 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 by 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 input for subsequent stress fatigue analysis of the reactor component.

[0054] In a possible embodiment, 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 speed; 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 coolant's acoustic parameters include, but are not limited to, density and sound velocity. During implementation, the coolant's acoustic parameter values can be determined based on the actual coolant's hot operating conditions. Specifically, they can be set based on the coolant's cold and hot segments, or based on the average temperature and pressure during reactor operation.

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

[0057] In one possible embodiment, the attenuation coefficient is generated according to the following method: using a preset model selection rule, based on 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 hysteretic 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, damping characteristic data corresponding to the target flow channel is calculated, and based on the damping characteristic data, a corresponding attenuation coefficient is calculated;

[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 i-th 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 is the speed of sound, ζ is the damping ratio, α and β are the preset coefficients, ω is the frequency, and const is a constant.

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

[0071]

[0072] Where ρ is the fluid density, u is the fluid velocity, and R is the transmission loss during sound wave transmission, which determines the damping characteristics of the reactor acoustic field. Based on the above definitions, 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, a preset model selection rule can be used to determine a target damping model from among viscous damping models, mass damping models, stiffness damping models, Rayleigh damping models, and hysteretic damping models, based on 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, parameters such as flow channel shape, flow channel dimensions, the presence of dead water areas, and flow channel area, which are determined by the preset model selection rule. The preset model selection rule can be set according to actual needs and is not limited to this in the present embodiment.

[0076] In one possible embodiment, the flow channel attribute 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 attribute data of each flow channel in the three-dimensional geometric model of the reactor assembly to be calculated, including: based on each flow channel attribute data, judging 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, 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, 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, which can be set according to actual conditions. For example, the second flow channel can also include a porous medium pore flow channel.

[0078] When using the above formula (1) or (2), when calculating the damping characteristic data of the first flow channel in the dead water area, the lack of coolant flow pressure drop data results in the calculated damping being 0, which is inaccurate. When calculating the damping characteristic data of the second flow channel, the cross section of the second flow channel is more complex and the calculation is also inconvenient. 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. The technician selects the damping model based on engineering experience and sends the selection result to the execution body of this embodiment, which uses the selection result as the target damping model.

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

[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 obtain:

[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 areas 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 pile internals 8.09

[0108] This embodiment is analogous to structural dynamics and creatively proposes a damping model system, which includes five models (viscous damping model, mass damping model, stiffness damping model, Rayleigh damping model, and hysteresis damping model). During implementation, a damping model that is more in line with reality 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 to 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 three-dimensional sound fields, and the calculation results are more accurate.

[0109] In response to a user-initiated boundary condition setting request, the acquired boundary characteristic data can be used to set boundary conditions for the three-dimensional geometric model of the reactor assembly to be calculated, thereby obtaining the three-dimensional model of the reactor assembly to be calculated. Boundary conditions include, but are not limited to, perforated plate boundaries and boundaries between different media. Boundary characteristic data is used to define boundary conditions. In implementation, for perforated plate boundaries, boundary characteristic data may include resistivity; for boundaries between different media, boundary characteristic data may include the normal impedance of the interface.

[0110] Specifically, in one possible implementation, 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, t p represents the plate thickness, a represents the hole radius, η represents the hole end correction coefficient, σ represents the porosity, R f Indicates flow resistivity.

[0113] This embodiment uses impedance characteristic boundaries 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 FIG2 shows a schematic diagram of a perforated plate impedance model provided by an embodiment of the present disclosure, 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, complex flow channels can also be simplified. For example, if the target component includes a reactor, the long and narrow flow channel of the fuel section can be used to consider only its plane wave propagation mode, while other high-order acoustic waves can be ignored.

[0115] This embodiment uses finite element software to construct a three-dimensional geometric model of the reactor assembly to be calculated, and sets acoustic parameter values, coolant damping characteristics, and complex acoustic boundary conditions for the three-dimensional geometric model, ultimately obtaining a three-dimensional model of the reactor assembly to be calculated. This makes 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 subsequently calculated based on the model.

[0116] In one possible implementation, the associated monitoring points include inner associated monitoring points and outer associated monitoring points; and 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: obtaining a plurality of first pulsating pressure loads from the frequency response data of the inner associated monitoring points and a plurality of second pulsating pressure loads from the frequency response data of the outer associated monitoring points; calculating each target load difference based on each of the first pulsating pressure loads and the corresponding second pulsating pressure loads; and generating the target frequency response data based on each target load difference.

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

[0118] For the target reactor component, the 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 basket) from the frequency response data of the inner associated monitoring point, and obtain multiple second pulsating pressure loads C21 (pressure pulsation outside the basket) from the frequency response data of the outer associated monitoring point. Then, we 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 of the excitation source. 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 utilizes associated monitoring points on the inside and outside of a 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 precise, thereby making the target frequency response data more accurate, thereby significantly improving the accuracy and reliability of the pulsating pressure transfer coefficient subsequently calculated based on the target frequency response data, and ultimately the pulsating pressure load on the target reactor component.

[0124] The present disclosure also provides a method for calculating a reactor pump-induced pulsating 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 the reactor component to be calculated is established based on a reactor in an actual reactor coolant system.

[0126] Step S2: In response to the 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. 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 the first flow channel and the second flow channel exist 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 the 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 based on 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 based on 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 user's request to set the perforated plate boundary conditions, the obtained impedance ratio is used to set the boundary conditions of the three-dimensional geometric model of the reactor assembly to be calculated, thereby obtaining 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: 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.

[0132] Step S7: Calculate target load differences based on the first pulsating pressure loads and the corresponding second pulsating pressure loads.

[0133] Step S8: Generate target frequency response data based on the target load differences.

[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 based on the amplitude and target pressure load corresponding to the target frequency.

[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 understand that in the above method of the specific embodiment, the writing order of each step does not mean a strict execution order, but does not constitute any limitation on the implementation process. The specific execution order of each step should be determined by its function and possible internal logic.

[0138] It should be noted that, in practical applications, all possible implementation methods described above can be arbitrarily combined to form possible embodiments of the present disclosure, and 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 pulsating pressure load caused by a reactor pump. Figure 5 The schematic diagram of the structure of the reactor pump pulsation pressure load calculation device provided by the embodiment of the present disclosure is shown. Figure 5 As shown, the reactor pump-induced pulsation pressure load calculation device 500 provided by the embodiment of the present disclosure includes:

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

[0141] Simulation module 502 is configured to perform acoustic simulation calculations on the three-dimensional model of the reactor assembly to be calculated, using the pulsating pressure at the inlet and outlet of the main pump as a simulation excitation source, to obtain frequency response data corresponding to each monitoring point in the model, wherein the response data in the frequency response data represents the pulsating pressure load borne by the monitoring point;

[0142] A frequency response data generating module 503 is configured 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 points of the target reactor component;

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

[0144] The load calculation module 505 is configured 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, in the above device, the target component includes at least one of the following: a reactor, a steam generator.

[0146] In a possible embodiment, in the above-mentioned device, the model construction module 501 is further used to: in response to a geometric model establishment request initiated by a user, establish 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, use 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, and the acoustic parameters include density and sound speed; in response to a damping setting request initiated by a user, use 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, use the acquired boundary characteristic data to set the boundary conditions of the three-dimensional geometric model of the reactor component to be calculated, and obtain the three-dimensional model of the reactor component to be calculated.

[0147] In one possible embodiment, the device further includes a generation module, configured to: determine a target damping model based on 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 hysteretic damping model; for any target flow channel in the three-dimensional geometric model of the reactor assembly to be calculated, calculate the damping characteristic data corresponding to the target flow channel based on the target damping model, and calculate the corresponding attenuation coefficient based on the damping characteristic data; the calculation of the damping characteristic data corresponding to the target flow channel based on the target damping model and the calculation of the corresponding attenuation coefficient based on 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 i-th 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 is the speed of sound, ζ is the damping ratio, α and β are the preset coefficients, ω is the frequency, and const is a constant.

[0160] In one possible embodiment, 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 one possible implementation, 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, t p represents the plate thickness, a represents the hole radius, η represents the hole end correction coefficient, σ represents the porosity, R f Indicates flow resistivity.

[0164] In one possible implementation, in the above-mentioned device, the associated monitoring points include inner associated monitoring points and outer associated monitoring points; the frequency response data generation module 503 is configured to: obtain multiple first pulsating pressure loads from the frequency response data of the inner associated monitoring points and multiple second pulsating pressure loads from the frequency response data of the outer associated monitoring points; calculate a target load difference based on each of the first pulsating pressure loads and the corresponding second pulsating pressure loads; and generate the target frequency response data based on 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 described in detail 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, an internal bus, a network interface, and memory. The memory may include internal memory, such as high-speed random-access memory (RAM), and may also include non-volatile memory, such as at least one disk storage device. Of course, the electronic device may also include other hardware required for its 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. Buses can be divided into address buses, data buses, control buses, 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 one type of bus.

[0168] The memory is used to store programs. Specifically, the program may include program code, which includes computer operating 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 internal memory and then runs it, forming a device for calculating the reactor pump-induced pulsating pressure load at a logical level. The processor executes the program stored in the memory and is specifically used to perform the aforementioned method.

[0170] The processor may be an integrated circuit chip with signal processing capabilities. During implementation, each step of the above method can be completed by hardware integrated logic circuits in the processor or by software instructions. 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 gate or transistor logic devices, or discrete hardware components. The various methods, steps, and logic block diagrams disclosed in the embodiments of the present disclosure can be implemented or executed. The general-purpose processor may be a microprocessor or any conventional processor. The steps of the method disclosed in conjunction with the embodiments of the present disclosure can be directly implemented and executed by a hardware decoding processor, or by a combination of hardware and software modules in the decoding processor. The software module can be located in a storage medium mature in the art, such as random access memory, flash memory, read-only memory, programmable read-only memory, electrically erasable programmable memory, registers, 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 embodiment will not be described in detail in the embodiments of the present disclosure.

[0172] An embodiment of the present disclosure also proposes 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. Furthermore, 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 magnetic 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 produce 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 that can direct a computer or other programmable data processing device to work 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 The 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 operational steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing the instructions executed on the computer or other programmable device for implementing the process. Figure 1 a process or multiple processes and / or boxes Figure 1 A step that specifies a function in one or more boxes.

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

[0178] 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. Memory is an example of a computer-readable medium.

[0179] Computer-readable media includes permanent and non-permanent, removable and non-removable media that can be implemented by any method or technology to store information. The 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 disc read-only memory (CD-ROM), digital versatile disc (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 transitory computer-readable media (transitory media), such as modulated data signals and carrier waves.

[0180] An embodiment of the present disclosure also provides a computer program product, which carries program code. The instructions included in the program code can be used to execute the steps of the reactor pump-induced pulsating pressure load calculation method described in the above method embodiment. For details, please refer to the above method embodiment and will not be repeated here.

[0181] The computer program product may be implemented in hardware, software, or a combination thereof. In one embodiment, the computer program product is implemented as a computer storage medium. In another 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 "comprises," "includes," or any other variations thereof are intended to encompass non-exclusive inclusion, such that a process, method, commodity, or apparatus that includes a series of elements includes not only those elements but also other elements not explicitly listed, or includes elements inherent to such process, method, commodity, or apparatus. In the absence of further limitations, an element defined by the phrase "comprises a ..." does not exclude the presence of other identical elements in the process, method, commodity, or apparatus that includes the element.

[0183] Those skilled in the art will appreciate that embodiments of the present disclosure may be provided as methods, systems, or computer program products. Thus, the present disclosure may take the form of a fully hardware embodiment, a fully software embodiment, or an embodiment combining software and hardware. Furthermore, 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 magnetic disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0184] The above are merely examples of the present disclosure and are not intended to limit the present disclosure. Various modifications and variations are possible for those skilled in the art. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present disclosure are intended to be included within the scope of the claims of the present disclosure.

Claims

1. A method for calculating reactor pump-induced pulsating pressure load, characterized in that: The method comprises: In response to a model building request initiated by a user, a three-dimensional model of the 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, an acoustic simulation calculation is performed on the three-dimensional model of the reactor assembly to be calculated, and frequency response data corresponding to each monitoring point in the model is obtained. 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; Determining target pressure loads corresponding to a plurality of characteristic frequencies in an actual excitation source based on the target frequency response data; 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 an amplitude corresponding to the target frequency and a target pressure load; Calculating 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; Wherein, the target component includes at least one of the following: a reactor, a steam generator; If the target component includes the reactor, then establishing a three-dimensional model of the reactor component to be calculated based on the target component in the actual reactor coolant system includes: establishing a three-dimensional model of the reactor component to be calculated based on the reactor pressure vessel, the reactor core and the reactor internal components, and the main pipeline; If the target component includes the steam generator, then establishing the three-dimensional model of the reactor component to be calculated based on the target component in the actual reactor coolant system includes: establishing the three-dimensional model of the reactor component to be calculated based on the steam generator and the main pipeline; 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, using the acquired acoustic parameter values, setting acoustic parameters of the coolant in the three-dimensional geometric model of the reactor assembly 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, setting a damping characteristic of a 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.

2. The method according to claim 1, characterized in that The attenuation coefficient is generated according to the following method: Determining a target damping model based on 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 hysteretic 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): （1） （2） （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): （4） （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): （6） 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): （7） 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): （8） Where R is the damping parameter, is the pressure drop of the coolant flow, 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 i-th flow channel, A i is the area of the i-th flow channel, N is the number of flow channels, is the preset attenuation coefficient, c represents the speed of sound, is the damping ratio, and is the preset coefficient, is the frequency, Represents a constant.

3. The method according to claim 2, characterized in that The flow channel attribute data includes the flow channel shape and whether there is a dead water area; and determining the target damping model based on the flow channel attribute data of each flow channel in the three-dimensional geometric model of the reactor assembly to be calculated by using a preset model selection rule includes: Based on the attribute data of each flow channel, determining whether a first flow channel and a second flow channel exist 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.

4. The method according to claim 1, wherein The boundary conditions include perforated plate boundary conditions, and the boundary characteristic data include impedance ratio; The impedance ratio is determined according to the following formula: in, , represents the transmission impedance, represents the fluid density, represents the fluid sound velocity, represents the dynamic viscosity, represents the shear (viscous) wave number, Indicates the plate thickness, represents the hole radius, represents the hole end correction factor, represents the porosity, Indicates flow resistivity.

5. The method according to any one of claims 1 to 4, characterized in that: 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 frequency response data corresponding to the associated monitoring points of the target reactor component, including: 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 target load differences based on 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.

6. A reactor pump-induced pulsation pressure load calculation device, characterized in that: The device comprises: A model building module is used to build a three-dimensional model of the reactor component to be calculated based on the target component in the actual reactor coolant system in response to a model building request initiated by the user; a simulation module for performing acoustic simulation calculations on the three-dimensional model of the reactor assembly to be calculated, using the pulsating pressure at the inlet and outlet of the main pump as a simulation excitation source, to 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, configured 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 points of the target reactor component; a transfer coefficient calculation module, configured to determine, based on the target frequency response data, target pressure loads corresponding to a plurality of characteristic frequencies in the actual excitation source; and for any target frequency among the plurality of characteristic frequencies, calculate, based on the amplitude corresponding to the target frequency and the target pressure load, a target pressure transfer coefficient of the target reactor component at the target frequency; a load calculation module, configured 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; Wherein, the target component includes at least one of the following: a reactor, a steam generator; If the target component includes the reactor, the model building module is used to: establish a three-dimensional model of the reactor component to be calculated based on the reactor pressure vessel, reactor core, reactor internals, and main pipelines; If the target component includes the steam generator, the model building module is used to: establish a three-dimensional model of the reactor component to be calculated based on the steam generator and the main pipeline; The model building module is specifically used to: respond to a geometric model establishment request initiated by a user, and establish a three-dimensional geometric model of the reactor component to be calculated based on the target component; respond to a parameter setting request initiated by a user, and use 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, and the acoustic parameters include density and sound speed; respond to a damping setting request initiated by a user, and use 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; respond to a boundary condition setting request initiated by a user, and use the acquired boundary characteristic data to set the boundary conditions of the three-dimensional geometric model of the reactor component to be calculated, so as to obtain the three-dimensional model of the reactor component to be calculated.

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

8. 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 according to any one of claims 1 to 5.

Images