Dynamic compensation method for correcting deviation of inherent performance prediction of annular pressurized water centrifugal pump

Abstract

The present application relates to the chemical pump engineering technical field, provide a kind of annular pressurized water centrifugal pump inherent performance prediction deviation dynamic offset type correction method, specific steps are: obtaining annular pressurized water centrifugal pump characteristic parameter and carrying out spatial dispersion to its calculation domain;Through conventional steady-state numerical simulation, obtain preliminary performance prediction value;Introduce energy offset correction factor to reflect the dynamic distribution characteristics of inherent performance prediction deviation with flow;Energy offset correction factor is used to determine energy offset correction source item to depict the equivalent momentum required to compensate inherent performance prediction deviation;Steady-state numerical simulation is carried out using energy offset correction source item to obtain performance correction value.The present application can effectively adaptively offset correction to its inherent performance prediction deviation by introducing energy offset correction factor, solves the existing numerical simulation method in predicting annular pressurized water centrifugal pump performance due to inherent deviation and seriously affects the engineering problem of predication, provides technical support for accurately grasping product quality.

Description

Technical Field

[0001] This invention relates to the field of chemical pump engineering technology, and in particular to a dynamic compensation correction method for predicting inherent performance deviations of an annular centrifugal water pump. Background Technology

[0002] Annular centrifugal pumps (or annular flow channel centrifugal pumps) offer engineering advantages such as compact structure, easy manufacturing, and reliable strength, making them widely used in the chemical industry for transporting high-temperature and high-pressure media, as well as in other industrial applications involving the transport of impurities. To ensure the good hydraulic performance of annular centrifugal pumps, numerical simulation methods are indispensable in modern hydraulic design for quickly and accurately predicting their performance. Compared to prototype testing, this important auxiliary method provides rapid feedback on pump performance, thereby reducing economic costs and shortening the development cycle.

[0003] Specifically, using the design flow rate as a reference condition, steady-state CFD calculations with the same grid are conducted within the conventional hydraulic design flow range of -50% to +130% to predict the head, shaft power, and efficiency of annular centrifugal pumps. This is currently the most common method in hydraulic design, and can even be considered a near-standard operating mode. Its purpose is to quickly understand basic performance and ensure efficient product development. However, existing practical experience shows that annular centrifugal pumps have unique characteristics of high flow rate, high speed, and strong turbulence. The performance values ​​of annular centrifugal pumps predicted using this method have a large inherent deviation from experimental results, especially in head and shaft power, where there is a constant positive deviation, and this deviation increases with increasing flow rate. Currently, there is no effective solution to this problem in related technologies, making this inherent performance prediction deviation seriously affect designers' prediction of annular centrifugal pump performance, thereby reducing the quality and efficiency of hydraulic design.

[0004] To address this engineering problem, there is an urgent need to propose an effective method to correct inherent performance prediction biases, so as to better serve the hydraulic design and product development of annular centrifugal pumps. Summary of the Invention

[0005] This invention provides a dynamic compensation correction method for predicting the inherent performance deviation of an annular centrifugal pump, which solves the engineering problem of the inherent deviation affecting the prediction when predicting the performance of an annular centrifugal pump in existing methods, thereby providing technical support for quickly and accurately grasping product quality.

[0006] This invention provides a dynamic compensation correction method for predicting inherent performance deviations of an annular centrifugal water pump, comprising the following steps: Step 1: Determine the relevant characteristic parameters of the annular centrifugal pump, including at least the design flow rate. Q 0. Rated speedn and impeller outlet diameter D 2; Step 2: Spatial discretize the computational domain of the annular centrifugal pump; wherein, the computational domain of the annular centrifugal pump includes the impeller domain and the annular pressure water domain; Step 3: Use SST k-oh The turbulence model uses the mass flow rate at the inlet as the inlet boundary condition and the static pressure at the outlet as the outlet boundary condition, for a flow rate of 0.5. Q 0, 0.6 Q 0, 0.7 Q 0, 0.8 Q 0, 0.9 Q 0, 1.0 Q 0, 1.1 Q 0, 1.2 Q 0 and 1.3 Q Steady-state numerical simulations were performed under nine operating conditions to obtain predicted values ​​for the head, shaft power, and efficiency of the annular centrifugal pump under each operating condition. Step 4: Introduce energy compensation correction factors Oh This is used to reflect the dynamic distribution characteristics of the prediction deviation of the inherent performance of the annular centrifugal pump with flow rate; wherein, the energy compensation correction factor corresponds to the impeller domain in the annular centrifugal pump. Oh 叶 The determination requires two active parameters determined through power compensation iterative calculations and two directly specified inert parameters. The two active parameters are low-order values. a and high-order value b The two inert parameters are slope factors. c and Mileage Points d ; Corresponding to the energy compensation correction factor of the annular pressure water area in the annular pressure centrifugal pump. Oh 环 Determined through iterative calculation of head compensation; Step 5: Utilize the energy compensation correction factor from Step 4 Oh 叶 and Oh 环 To determine the energy compensation correction source terms corresponding to the impeller domain and the annular pressure water domain, respectively. M 叶 and M 环 It is used to characterize the equivalent momentum required to compensate for the inherent performance prediction deviation of the annular centrifugal pump, thereby realizing the adaptive compensation correction of the energy field inside the annular centrifugal pump. Step Six: Continue with SST from Step Three k-oh Based on the turbulence model, the energy compensation correction source term in step five is... M 叶 andM 环 Add these values ​​to the momentum equations for the impeller domain and the annular pressure domain, respectively. Use the mass flow rate inlet from step three as the outlet boundary condition and the static pressure outlet as the outlet boundary condition. For 0.5... Q 0~1.3 Q Steady-state numerical simulations are performed on any flow conditions of interest within the range of 0 to obtain the corrected values ​​of head, shaft power, and efficiency of the annular centrifugal pump under the corresponding operating conditions.

[0007] According to the dynamic compensation correction method for the inherent performance prediction deviation of the annular centrifugal pump provided by the present invention, step two specifically includes the following steps: Tetrahedral meshes are used to spatially discretize the impeller domain and annular pressure water area of ​​the annular centrifugal pump. The grid size of the impeller domain is taken as the impeller outlet diameter. D 2 / 40, the grid size of the annular pressure water area is taken as the impeller outlet diameter. D 1 / 20 of 2.

[0008] According to the dynamic compensation correction method for the inherent performance prediction deviation of the annular centrifugal pump provided by the present invention, the energy compensation correction factor corresponding to the impeller domain in step four is... Oh 叶 The calculation formula is as follows: In the formula, a The lower-order value, b For the aforementioned high-order value, c The slope factor is mentioned above. d For the aforementioned mileage points, For Strauha numbers.

[0009] According to the dynamic compensation correction method for the inherent performance prediction deviation of the annular pressurized water centrifugal pump provided by the present invention, the energy compensation correction source terms corresponding to the impeller domain and the annular pressurized water domain in step five are as follows: M 叶 and M 环 The calculation formulas are as follows: ; ; In the formula, r For water flow density, V 叶x , V 叶y , V 叶z For the three-dimensional components of the relative velocity in the impeller domain, V 环x ,V 环y , V 环z For the three-dimensional components of the absolute velocity of the annular pressure water body, i , j , k is the basis vector of the three-dimensional rectangular coordinate system.

[0010] According to the dynamic compensation correction method for the inherent performance prediction deviation of the annular centrifugal pump provided by the present invention, the low-position value is determined by power compensation iterative calculation in step four. a The specific process is as follows: Follow the method in step six for 0.5 Q Steady-state calculations were performed under operating condition 0 to determine the energy compensation correction factor for the impeller domain. Oh 叶 The initial value is 0, and it increases by 0.0005 every 50 iterations until the calculated shaft power is equal to 0.5 in step three. Q The predicted shaft power value for condition 0 decreases by 1%, and at this time... Oh 叶 The iteration result is assigned to the low-order value. a Used to calculate the impeller domain energy compensation correction factor. Oh 叶 .

[0011] According to the dynamic compensation correction method for the inherent performance prediction deviation of the annular centrifugal pump provided by the present invention, the high-order value is determined by power compensation iterative calculation in step four. b The specific process is as follows: Following the method in step six, a steady-state trial calculation is performed on the optimal operating condition predicted in step three, and an energy compensation correction factor is assigned to the impeller domain. Oh 叶 High value in b The initial value is set to 0, and then increased by 0.0005 every 50 iterations until the calculated shaft power decreases by 7% compared to the predicted shaft power under the optimal operating condition in step three. b The iterative result is used as the high-order value required in step five to calculate the impeller domain energy compensation correction factor. Oh 叶 .

[0012] According to the dynamic compensation correction method for the inherent performance prediction deviation of the annular pressurized centrifugal pump provided by the present invention, in step four, the energy compensation correction factor corresponding to the annular pressurized water area is determined by iterative calculation of head compensation. Oh 环 The specific process is as follows: Following the method in step six, a steady-state trial calculation is performed on the optimal operating condition predicted in step three, and an energy compensation correction factor is applied to the annular pressure water area. Oh 环 The initial value is 0, and it increases by 0.0005 every 50 iterations until the calculated head decreases by 7% compared to the head prediction of the optimal working condition in step three. Oh 环 The iterative results serve as the energy compensation correction factor for the annular pressure water area required in step five. Oh 环 .

[0013] According to the dynamic compensation correction method for the prediction deviation of the inherent performance of the annular centrifugal pump provided by the present invention, the slope factor in step four... c =147, Mileage Points

[0014] Through any of the above embodiments, the present invention has at least the following beneficial effects: The present invention provides a dynamic compensation correction method for the inherent performance prediction deviation of annular pressurized centrifugal pumps. By introducing an energy compensation correction factor, for a given annular pressurized centrifugal pump, it can adaptively compensate and correct the inherent performance prediction deviation within the flow range of interest in conventional hydraulic design. This achieves dynamic adjustment and accurate feedback of the hydraulic performance prediction value of the annular pressurized centrifugal pump. It solves the engineering problem in related technologies where the inherent deviation of numerical simulation methods in predicting the performance of annular pressurized centrifugal pumps seriously affects the performance prediction. Thus, it provides technical support for quickly and accurately grasping product quality and improving hydraulic design efficiency. Attached Figure Description

[0015] To more clearly illustrate the technical solutions in this invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of this invention. For those skilled in the art, other drawings can be obtained from these drawings without creative effort.

[0016] Figure 1 This is a flowchart illustrating the implementation of the dynamic compensation correction method for predicting the inherent performance deviation of an annular centrifugal pump provided by the present invention. Figure 2 This is a schematic diagram of the flow calculation domain of an embodiment of the annular centrifugal pump provided by the present invention; Figure 3 This is a comparison chart of the head prediction results of an embodiment of the annular centrifugal pump provided by the present invention; Figure 4 This is a comparison chart of shaft power prediction results for an embodiment of the annular centrifugal pump provided by the present invention. Detailed Implementation

[0017] To make the objectives, technical solutions, and advantages of this invention clearer, the technical solutions of this invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some, not all, of the embodiments of this invention. All other embodiments obtained by those skilled in the art based on the embodiments of this invention without creative effort are within the scope of protection of this invention.

[0018] It should be noted that SST in the description of this invention k-oh The turbulence model is the most widely used two-equation eddy viscosity model in the numerical simulation of centrifugal pumps. Those skilled in the art, after understanding the technical solution of this invention, can utilize SST... k-oh Turbulence model implementation, therefore for SST k-oh The specific formation of the turbulence model will not be elaborated upon.

[0019] The following is combined Figure 1 This invention describes a dynamic compensation correction method for predicting inherent performance deviations in annular centrifugal pumps. This compensation correction method is applied to the performance prediction of annular centrifugal pumps. Through this compensation correction method, accurate performance prediction can be achieved, and prediction deviations can be reduced. Specifically, it includes the following steps: Step 1: Determine the relevant characteristic parameters of the annular centrifugal pump, including at least the design flow rate. Q 0. Rated speed n and impeller outlet diameter D 2; Step 2: Spatial discretize the computational domain of the annular centrifugal pump; wherein, the computational domain of the annular centrifugal pump includes the impeller domain and the annular pressure water domain; Step 3: Use SST k-oh The turbulence model uses the mass flow rate at the inlet as the inlet boundary condition and the static pressure at the outlet as the outlet boundary condition, for a flow rate of 0.5. Q 0, 0.6 Q 0, 0.7 Q 0, 0.8 Q 0, 0.9 Q 0, 1.0 Q 0, 1.1 Q 0, 1.2 Q 0 and 1.3 Q Steady-state numerical simulations were performed under nine operating conditions to obtain predicted values ​​for the head, shaft power, and efficiency of the annular centrifugal pump under each operating condition. Step 4: Introduce energy compensation correction factors OhThis is used to reflect the dynamic distribution characteristics of the prediction deviation of the inherent performance of the annular centrifugal pump with flow rate; wherein, the energy compensation correction factor corresponds to the impeller domain in the annular centrifugal pump. Oh 叶 The determination requires two active parameters determined through power compensation iterative calculations and two directly specified inert parameters. The two active parameters are low-order values. a and high-order value b The two inert parameters are slope factors. c and Mileage Points d ; Corresponding to the energy compensation correction factor of the annular pressure water area in the annular pressure centrifugal pump. Oh 环 Determined through iterative calculation of head compensation; Step 5: Utilize the energy compensation correction factor from Step 4 Oh 叶 and Oh 环 To determine the energy compensation correction source terms corresponding to the impeller domain and the annular pressure water domain, respectively. M 叶 and M 环 It is used to characterize the equivalent momentum required to compensate for the inherent performance prediction deviation of the annular centrifugal pump, thereby realizing the adaptive compensation correction of the energy field inside the annular centrifugal pump. Step Six: Continue with SST from Step Three k-oh Based on the turbulence model, the energy compensation correction source term in step five is... M 叶 and M 环 These are added to the momentum equations of the impeller domain and the annular pressure domain, respectively, using the mass flow rate inlet and static pressure outlet from step three as the inlet and outlet boundary conditions, respectively, for 0.5 Q 0~1.3 Q Steady-state numerical simulations are performed on any flow conditions of interest within the range of 0 to obtain the corrected values ​​of head, shaft power, and efficiency of the annular centrifugal pump under the corresponding operating conditions.

[0020] In step one of the above embodiments, the parameters are determined and obtained by using a model of a real-world annular centrifugal pump. For example, the rotational speed is one of the characteristic parameters. n and impeller outlet diameter D 2. This can be determined during the design phase of the annular centrifugal pump.

[0021] In step two, the computational domain for the annular pressurized centrifugal pump specifically comprises the impeller domain and the annular pressurized water area, and the annular pressurized centrifugal pump can be spatially discretized using tetrahedral or hexahedral meshes. It should be understood that the impeller domain and the annular pressurized water area are inherent computational domains of the annular pressurized centrifugal pump.

[0022] In step three, SST can be preset in the system. k-oh The turbulence model invokes SST via a call when executing the method of this invention. k-oh The turbulence model was used to obtain preliminary predictions and subsequent corrections for the annular pressurized centrifugal pump.

[0023] The following specific examples illustrate this: Example 1 Step 101: Obtain the relevant characteristic parameters of the annular centrifugal water pump, including the design flow rate. Q 0. Rated speed n and impeller outlet diameter D 2. Among them, design flow rate Q 0 = 0.639m 3 / s, rated speed n =151.84s -1 Impeller outlet diameter D 2 = 0.463m.

[0024] Step 102: Spatial discretization of the computational domain of the annular centrifugal pump is performed using a tetrahedral mesh; wherein, the mesh size of the impeller domain is taken as the impeller outlet diameter. D 2 / 40, the grid size of the annular pressure water area is taken as the impeller outlet diameter. D 1 / 20 of 2; Step 103: Use SST, the most widely used method in numerical simulation of centrifugal pumps. k-oh The turbulence model uses the mass flow rate at the inlet as the inlet boundary condition and the static pressure at the outlet as the outlet boundary condition, for a flow rate of 0.5. Q 0, 0.6 Q 0, 0.7 Q 0, 0.8 Q 0, 0.9 Q 0, 1.0 Q 0, 1.1 Q 0, 1.2 Q 0 and 1.3 Q Nine operating conditions were simulated in steady state to obtain the initial predicted values ​​of head, shaft power and efficiency of the annular centrifugal pump under each operating condition. Step 104: Introduce energy compensation correction factor OhThis is used to reflect the dynamic distribution characteristics of the prediction deviation of the inherent performance of an annular centrifugal pump with flow rate; among which, the energy compensation correction factor corresponds to the impeller domain. Oh 叶 The determination requires two active parameters determined through power compensation iterative calculations and two directly specified inert parameters. The two active parameters are the low-order values. a and high-order value b The two inert parameters are the slope factor. c and Mileage Points d ; energy compensation correction factor corresponding to the annular pressure water area Oh 环 Determined through iterative calculation of head compensation; Step 105: Utilize the energy compensation correction factor from step 104 Oh 叶 and Oh 环 To determine the energy compensation correction source terms corresponding to the impeller domain and the annular pressure water domain, respectively. M 叶 and M 环 Among them, the M 叶 and M 环 It is used to characterize the equivalent momentum required to compensate for the inherent performance prediction deviation of the annular centrifugal pump, thereby realizing the adaptive compensation correction of the energy field inside the annular centrifugal pump. Step 106: Continue with the SST from step 103 k-oh Based on the turbulence model, the energy compensation correction source term in step 105 is... M 叶 and M 环 The momentum equations for the impeller domain and the annular pressure water domain were added respectively to correct the computational domain. The mass flow rate inlet and static pressure outlet from step 103 were used as the inlet and outlet boundary conditions for 0.5... Q 0~1.3 Q Corrected steady-state numerical simulations are performed on any flow condition of interest within the range of 0 to obtain the corrected values ​​of head, shaft power and efficiency of the annular centrifugal pump under the corresponding operating conditions.

[0025] In the above embodiment, the energy compensation correction factor corresponding to the impeller domain in step 104 Oh 叶 The calculation formula is as follows: In the formula, a These are the low-order values ​​of the activity parameters. b These are the high-order values ​​of the activity parameters.c The slope factor is an inert parameter. d Odometer points for inert parameters For Strauha numbers.

[0026] In the above embodiments, the energy compensation correction source terms corresponding to the impeller domain and the annular pressure water domain in step 5 are... M 叶 and M 环 The calculation formulas are as follows: ; ; In the formula, r For water flow density, V 叶x , V 叶y , V 叶z For the three-dimensional components of the relative velocity in the impeller domain, V 环x , V 环y , V 环z For the three-dimensional components of the absolute velocity of the annular pressure water body, i , j , k is the basis vector of the three-dimensional rectangular coordinate system.

[0027] In the above embodiment, the low-bit value is determined by power compensation iterative calculation in step 104. a The specific process is as follows: Follow the method in step 106 to process 0.5 Q Steady-state calculations were performed under operating condition 0 to determine the energy compensation correction factor for the impeller domain. Oh 叶 The initial value is 0, and it increases by 0.0005 every 50 iterations until the calculated shaft power is equal to 0.5 as described in step 103. Q The predicted shaft power value for condition 0 decreases by 1%, and at this time... Oh 叶 The iteration result is assigned to the low-order value. a Used to calculate the impeller domain energy compensation correction factor. Oh 叶 Specifically, based on the characteristic parameters determined above, for this embodiment, the final determined low-order value... a =0.007.

[0028] In the above embodiment, the high-order value is determined by power compensation iterative calculation in step 104. b The specific process is as follows: Following the method in step 106, a steady-state trial calculation is performed on the optimal operating condition (the predicted flow point with the highest efficiency) obtained in step 103, and an energy compensation correction factor is applied to the impeller domain. Oh 叶 High value in b The initial value is set to 0, and then increased by 0.0005 every 50 iterations until the calculated shaft power decreases by 7% compared to the predicted shaft power value for the optimal operating condition obtained in step 103. At this point... b The iteration result is used as the high-order value required in step 105. b Used to calculate the impeller domain energy compensation correction factor. Oh 叶 Specifically, based on the characteristic parameters determined above, for this embodiment, the final determined high-order value is... b =0.03.

[0029] In the above embodiment, step 104 determines the energy compensation correction factor corresponding to the annular pressure water area through iterative calculation of head compensation. Oh 环 The specific process is as follows: Following the method in step 106, a steady-state trial calculation is performed on the optimal operating condition (the predicted flow point with the highest efficiency) obtained in step 103, and an energy compensation correction factor is applied to the annular pressure water area. Oh 环 The initial value is 0, and it increases by 0.0005 every 50 iterations until the calculated head decreases by 7% compared to the head prediction value of the optimal working condition obtained in step 103. At this point... Oh 环 The iterative results serve as the energy compensation correction factor for the annular pressure water area required in step 105. Oh 环 Specifically, based on the characteristic parameters determined above, the energy compensation correction factor for the annular pressure water area is finally determined for this embodiment. Oh 环 =0.003.

[0030] Furthermore, in step 104, the preferred value for the specified inertial parameter slope factor is... c =147, the preferred value for mileage points is...

[0031] Understandably, the grid size of the impeller domain and the grid size of the annular pressure water domain can be appropriately adjusted around the recommended values ​​based on actual needs to better meet the prediction requirements of different scenarios.

[0032] To further illustrate the solution of the present invention, a comparison is made between relevant prediction methods and the prediction method described below: like Figure 3 As shown, this is a comparison chart of the head prediction results for an embodiment of the annular centrifugal pump provided by the present invention. Figure 4 The figure shown is a comparison of shaft power prediction results for an embodiment of the annular centrifugal pump provided by the present invention. Analysis shows that the head value of the annular centrifugal pump predicted using existing methods (…) H 0 represents the head value at the design flow rate measured in a high-precision experiment, and the shaft power value. P The shaft power values ​​(0 for the design flow rate measured in high-precision experiments) show a significant positive deviation from the experimental values, and this deviation increases with increasing flow rate, reaching a maximum of nearly 20%. This clearly affects designers' predictions of the hydraulic performance of the annular centrifugal pump. In contrast, the head and shaft power values ​​of the annular centrifugal pump calculated using the method of this invention are much closer to the high-precision experimental values, especially within the range of -20%. Q 0~+130% Q Within the flow rate range of 0. As can be seen from the above comparison, this invention effectively corrects the inherent performance prediction bias of existing methods and significantly improves the prediction accuracy of the performance of annular centrifugal pumps.

[0033] Through the above description of the embodiments, those skilled in the art can clearly understand that the present invention can adaptively compensate and correct the inherent performance prediction deviation within the flow range of interest in conventional hydraulic design for a given annular centrifugal pump. It realizes dynamic adjustment and accurate feedback of the hydraulic performance prediction value of the annular centrifugal pump, and solves the engineering problem that the existing numerical simulation method seriously affects the performance prediction due to inherent deviations when predicting the performance of annular centrifugal pumps. Thus, it provides technical support for quickly and accurately grasping product quality and improving hydraulic design efficiency.

[0034] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A dynamic compensation correction method for predicting inherent performance deviations of an annular centrifugal water pump, characterized in that, Includes the following steps: Step 1: Determine the relevant characteristic parameters of the annular centrifugal pump, including at least the design flow rate. Q 0. Rated speed n and impeller outlet diameter D 2; Step 2: Spatial discretize the computational domain of the annular centrifugal pump; Step 3: Use SST k-ω The turbulence model uses the mass flow rate at the inlet as the inlet boundary condition and the static pressure at the outlet as the outlet boundary condition, for a flow rate of 0.

5. Q 0, 0.6 Q 0, 0.7 Q 0, 0.8 Q 0, 0.9 Q 0, 1.0 Q 0, 1.1 Q 0, 1.2 Q 0 and 1.3 Q Steady-state numerical simulations were performed under nine operating conditions to obtain predicted values ​​for the head, shaft power, and efficiency of the annular centrifugal pump under each operating condition. Step 4: Introduce energy compensation correction factors Ω This is used to reflect the dynamic distribution characteristics of the prediction deviation of the inherent performance of the annular centrifugal pump with flow rate; wherein, the energy compensation correction factor corresponds to the impeller domain in the annular centrifugal pump. Ω 叶 The determination requires two active parameters determined through power compensation iterative calculations and two directly specified inert parameters. The two active parameters are low-order values. a and high-order value b The two inert parameters are slope factors. c and Mileage Points d ; Corresponding to the energy compensation correction factor of the annular pressure water area in the annular pressure centrifugal pump. Ω 环 Determined through iterative calculation of head compensation; Step 5: Utilize the energy compensation correction factor from Step 4 Ω 叶 and Ω 环 To determine the energy compensation correction source terms corresponding to the impeller domain and the annular pressure water domain, respectively. M 叶 and M 环 It is used to characterize the equivalent momentum required to compensate for the inherent performance prediction deviation of the annular centrifugal pump, thereby realizing the adaptive compensation correction of the energy field inside the annular centrifugal pump. Step Six: Continue with SST from Step Three k-ω Based on the turbulence model, the energy compensation correction source term in step five is... M 叶 and M 环 Add these values ​​to the momentum equations for the impeller domain and the annular pressure domain, respectively, using the mass flow rate inlet from step three as the inlet boundary condition and the static pressure outlet as the outlet boundary condition, for 0.5... Q 0~1.3 Q Steady-state numerical simulations were performed on any flow rate condition of interest within the range of 0, thereby obtaining the corrected values ​​of head, shaft power, and efficiency of the annular centrifugal pump under the corresponding operating conditions; The energy compensation correction source terms corresponding to the impeller domain and the annular pressure water region in step five. M 叶 and M 环 The calculation formulas are as follows: ； ； In the formula, ρ For water flow density, V 叶x , V 叶y , V 叶z For the three-dimensional components of the relative velocity in the impeller domain, V 环x , V 环y , V 环z The three-dimensional components of the absolute velocity of the annular pressure water body. i , j , k is the basis vector of the three-dimensional rectangular coordinate system.

2. The dynamic compensation correction method for the prediction deviation of the inherent performance of the annular centrifugal pump according to claim 1, characterized in that, Step two specifically includes the following steps: Tetrahedral meshes are used to spatially discretize the impeller domain and annular pressure water area of ​​the annular centrifugal pump. The grid size of the impeller domain is taken as the impeller outlet diameter. D 2 / 40, the grid size of the annular pressure water area is taken as the impeller outlet diameter. D 1 / 20 of 2.

3. The dynamic compensation correction method for the prediction deviation of the inherent performance of the annular centrifugal pump according to claim 1, characterized in that, The energy compensation correction factor corresponding to the impeller domain in step four. Ω 叶 The calculation formula is as follows: In the formula, a The lower-order value, b For the aforementioned high-order value, c The slope factor is mentioned above. d For the aforementioned mileage points, Let be the Strauha number, where Q The flow rate of the annular centrifugal pump under the operating conditions of interest.

4. The dynamic compensation correction method for the prediction deviation of the inherent performance of the annular centrifugal pump according to claim 3, characterized in that, In step four, the low-order value is determined through power compensation iterative calculation. a The specific process is as follows: Follow the method in step six for 0.5 Q Steady-state calculations were performed under operating condition 0 to determine the energy compensation correction factor for the impeller domain. Ω 叶 The initial value is 0, and it increases by 0.0005 every 50 iterations until the calculated shaft power is equal to 0.5 in step three. Q The predicted shaft power value for condition 0 decreases by 1%, and at this time... Ω 叶 The iteration result is assigned to the low-order value. a Used to calculate the impeller domain energy compensation correction factor. Ω 叶 .

5. The dynamic compensation correction method for the prediction deviation of the inherent performance of the annular centrifugal pump according to claim 3, characterized in that, In step four, the higher bit value is determined through power compensation iterative calculation. b The specific process is as follows: Following the method in step six, a steady-state trial calculation is performed on the optimal operating condition predicted in step three, and an energy compensation correction factor is assigned to the impeller domain. Ω 叶 High value in b The initial value is set to 0, and then increased by 0.0005 every 50 iterations until the calculated shaft power decreases by 7% compared to the predicted shaft power under the optimal operating condition in step three. b The iterative result is used as the high-order value required in step five to calculate the impeller domain energy compensation correction factor. Ω 叶 .

6. The dynamic compensation correction method for the prediction deviation of the inherent performance of the annular centrifugal pump according to claim 1, characterized in that, In step four, the energy compensation correction factor corresponding to the annular pressure water area is determined through iterative calculation of head compensation. Ω 环 The specific process is as follows: Following the method in step six, a steady-state trial calculation is performed on the optimal operating condition predicted in step three, and an energy compensation correction factor is applied to the annular pressure water area. Ω 环 The initial value is 0, and it increases by 0.0005 every 50 iterations until the calculated head decreases by 7% compared to the head prediction of the optimal working condition in step three. At this point... Ω 环 The iterative results serve as the energy compensation correction factor for the annular pressure water area required in step five. Ω 环 .

7. The dynamic compensation correction method for the prediction deviation of the inherent performance of the annular centrifugal pump according to claim 1, characterized in that, Slope factor in step four c =147, Mileage Points

Images