CN110617238B - Optimization design method of centrifugal pump impeller - Google Patents

Optimization design method of centrifugal pump impeller Download PDF

Info

Publication number
CN110617238B
CN110617238B CN201910914770.XA CN201910914770A CN110617238B CN 110617238 B CN110617238 B CN 110617238B CN 201910914770 A CN201910914770 A CN 201910914770A CN 110617238 B CN110617238 B CN 110617238B
Authority
CN
China
Prior art keywords
impeller
centrifugal pump
flow
design
parameters
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910914770.XA
Other languages
Chinese (zh)
Other versions
CN110617238A (en
Inventor
李晓俊
刘遥遥
郇钰莹
蔡云辉
陈波
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University of Technology ZJUT
Original Assignee
Zhejiang University of Technology ZJUT
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhejiang University of Technology ZJUT filed Critical Zhejiang University of Technology ZJUT
Priority to CN201910914770.XA priority Critical patent/CN110617238B/en
Publication of CN110617238A publication Critical patent/CN110617238A/en
Application granted granted Critical
Publication of CN110617238B publication Critical patent/CN110617238B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F04POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
    • F04DNON-POSITIVE-DISPLACEMENT PUMPS
    • F04D29/00Details, component parts, or accessories
    • F04D29/18Rotors
    • F04D29/22Rotors specially for centrifugal pumps
    • F04D29/2238Special flow patterns
    • F04D29/2244Free vortex

Landscapes

  • Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Structures Of Non-Positive Displacement Pumps (AREA)

Abstract

The invention relates to the field of hydraulic optimization of centrifugal pumps. The method can more accurately and quickly determine the optimal working condition point of the centrifugal pump impeller under the common constraint of the vortex structure and the entropy product to obtain the corresponding impeller design parameters, so that the impeller has better hydraulic performance and better stability. The technical scheme is as follows: an optimal design method of a centrifugal pump impeller is characterized by comprising the following steps: the method comprises the following steps: 1) firstly, specifically determining a group of parameters to be optimized according to design requirements and design experience; 2) modeling and grid division; 3) CFX simulation verification; 4) extracting a vortex core structure; 5) analyzing the entropy yield; 6) efficient point matching; 7) and optimizing the design parameters of the impeller by adopting a simulated annealing optimization algorithm.

Description

Optimization design method of centrifugal pump impeller
Technical Field
The invention relates to the field of hydraulic optimization of centrifugal pumps, in particular to an optimization design method for an impeller of a centrifugal pump.
Background
At present, centrifugal pumps are the most widely used pump type in the pump industry, and the usage amount of the centrifugal pumps accounts for about 70% of the whole pump industry. The viscous action of the fluid in the centrifugal pump and the special geometric structure of the centrifugal pump determine that the actual flow inside the centrifugal pump is three-dimensional, viscous and abnormal, and in addition, the action of various factors such as the fluid flowing around the flow vanes and rotating at high speed inevitably forms various unstable flows, which mainly comprise: the unsteady factors of the centrifugal pump include the common occurrence and obvious action of the unsteady factors such as the flow shedding, the backflow, the secondary flow, the wake flow and various vortexes, so that the flow field inside the centrifugal pump presents strong and complex unsteady flow characteristics. At the design operating point, although the centrifugal pump keeps high operation efficiency, the existence of unsteady flow characteristics can also lead to poor operation stability and high vibration noise of the centrifugal pump. The impeller is the most important overflowing part and acting part in the centrifugal pump, the design quality of the impeller directly influences the operation efficiency and the operation stability of the pump, and further influences the design targets of vibration, noise and the like of the pump, so that the impeller of the pump is of great significance in optimization design.
For the hydraulic design of centrifugal pump impellers, three major problems mainly need to be solved: (1) the efficiency is improved; (2) performance is improved; (3) and the operation stability of the pump is enhanced. The invention patent with the application number of 201110202524.5, 7/20/2011 discloses an optimal design method for an impeller of an anti-cavitation centrifugal pump, which is used for designing an optimal design method for the impeller for improving cavitation erosion performance by constructing optimization parameters and adopting an NSGA-II genetic algorithm as an optimization tool, and the method does not consider factors influencing the running stability of the impeller and the like. The key point is to fully know the internal fluid motion state of the centrifugal pump to ensure safe, stable and efficient operation of the centrifugal pump, wherein separation flow and vortex motion are the most common flow phenomena, and the hydraulic fluctuation can generate pressure pulsation on the centrifugal pump to influence the operation stability of the centrifugal pump. Therefore, on the premise of fully considering the efficiency of the impeller and the pressure pulsation factor, an impeller optimization design method is provided to improve the hydraulic performance and stability of the centrifugal pump.
Disclosure of Invention
The invention aims to overcome the defects of the background technology and provides an optimal design method of a centrifugal pump impeller, which can more accurately and quickly determine the optimal working condition point of the centrifugal pump impeller under the common constraint of a vortex structure and an entropy product to obtain corresponding impeller design parameters, so that the impeller has better hydraulic performance and better stability.
The technical scheme provided by the invention is as follows:
an optimal design method of a centrifugal pump impeller is characterized by comprising the following steps: the method comprises the following steps:
1) firstly, according to design requirements, a group of parameters to be optimized is specifically determined according to design experience: width of impeller inlet b1Impeller exit width b2Blade inlet angle beta1Blade exit angle beta2The number of blades z, this set of parameters to be optimized gives the constraints below. Other impeller parameters are used as design constants;
the constraint range of the parameters to be optimized is as follows:
Figure GDA0002778916230000021
in the above formula:
Km1-the impeller inlet axial surface coefficient,
ns-the specific speed of rotation of the pump,
Fr-a flow rate of the fluid,
h-the lift is set as the height of the pump,
n-the rotating speed of the impeller,
g-acceleration of gravity;
when the design parameters of the centrifugal pump impeller are optimized in the step, the distribution condition of the energy loss inside the centrifugal pump in the step 4) can be analyzed, so that a main area of the energy loss of the impeller is obtained, and a reference is provided for the adjustment of the design parameters of the impeller corresponding to the low-efficiency area of the impeller, so that the pertinence of the optimization of the design parameters is improved, and the verification times and the workload of the design parameters are reduced.
2) Modeling and meshing: modeling the calculation domain of the centrifugal pump impeller under a determined coordinate system (shown in figure 2) by using three-dimensional modeling software according to the parameters of the centrifugal pump impeller; utilizing ANSYS ICEM software to perform meshing on the three-dimensional model;
3) and (3) CFX simulation verification: the grid is led into CFX software, transient numerical simulation calculation of the centrifugal pump impeller under the full-flow working condition is carried out under given boundary conditions (including inlet pressure, outlet flow, wall surface and dynamic and static interface setting and impeller rotating speed), and the result can be automatically generated by the software; carrying out post-processing on the numerical simulation result to obtain a flow-lift and flow-efficiency performance curve of the centrifugal pump; comparing whether the lift and the efficiency under the drawn grid meet the design requirements or not according to the performance curve, and if the design requirements are not met, repeating the step 2) to divide the grid of the three-dimensional model again until the grid meeting the design requirements is obtained;
4) extracting a vortex core structure: based on the simulation result obtained in the step 3), extracting the vortex core structure in the calculation domain based on a Q method to obtain the data of the vortex core strength changing along with the flow, wherein the Q value represents the strength relation between the fluid infinitesimal rotation and the deformation;
5) entropy production analysis: based on the grid obtained in the step 3), calculating the distribution condition of the internal energy loss of the centrifugal pump according to the speed distribution of the CFD calculation result according to the formulas (1.1), (1.2) and (1.3); then carrying out volume integration on the formulas (1.2) and (1.3) in an impeller region by using a formula (1.4) to obtain the total entropy yield of the impeller;
Figure GDA0002778916230000031
Figure GDA0002778916230000041
Figure GDA0002778916230000042
in the above formula
Spro,t-total entropy production rate per unit volume of impeller,
Figure GDA0002778916230000043
-mean entropy production rate per unit volume of impeller,
Spro,d′-a rate of impeller pulsation entropy production per unit volume,
t-the temperature of the fluid, and,
p-the density of the fluid,
-a turbulent dissipation ratio of the fluid,
u, v, w-velocity components of the fluid in the direction X, Y, Z respectively,
Figure GDA0002778916230000044
representing the time-average velocity value, u ', v ', w ' representing the pulsation velocity value;
6) efficient point matching: carrying out non-dimensionalization on the data obtained in the step 4) and the step 5), and importing the data into data processing software for drawing to obtain a curve of the total entropy of the impeller and the vortex core strength changing along with the flow; wherein the non-dimensional flow corresponding to the minimum vortex core strength is Fr1The dimensionless flow corresponding to the minimum total entropy yield of the impeller is Fr2
7) Optimizing design parameters of the impeller by adopting a simulated annealing optimization algorithm to ensure that the dimensionless flow corresponding to the minimum vortex core strength and the dimensionless flow corresponding to the minimum impeller total entropy yield are both F'rAnd F'r=Fr1=Fr2At this point, the optimum design parameters for the impeller are found to be 1.
In the step 4), the step of extracting the vortex structure by using the Q method comprises the following steps:
Figure GDA0002778916230000045
Figure GDA0002778916230000046
Figure GDA0002778916230000047
the above equation is simplified to the following Q-values in cartesian coordinates:
Figure GDA0002778916230000051
wherein, the Q method defines the area where Q > 0 as the area where the vortex tube is; the deformation rate tensor and the vorticity in the above formula are calculated by editing a CEL expression in POST-processing software CFD-POST according to a numerical simulation result. Setting a threshold value, and extracting Fr=FrdThe infinitesimal iso-surface of (A) determines the vortex morphology, wherein FrdTo design the flow rate. After the Q value is determined, in CFD-POST, a vorticity isosurface under the Q value is made, the vorticity coefficient is adopted for coloring, then the strength of the vortex core is quantitatively calculated by software, and data are recorded.
In the step 6), the total entropy yield of the impeller is subjected to dimensionless treatment, wherein the numerical value of the entropy yield of each impeller is divided by the minimum value of the entropy yield of the impeller; carrying out non-dimensionalization treatment on the strength of the vortex core, namely dividing the numerical value of the strength of each vortex core by the minimum value of the strength of the vortex core; the flow is subjected to non-dimensionalization treatment by dividing the numerical value of each flow by the design flow Frd
The invention has the beneficial effects that:
1) the method reflects the change condition of the efficiency of the impeller by analyzing the entropy yield value inside the impeller; describing the change of the pressure pulsation by numerically processing the vortex core intensity; the impeller is optimally designed through the combination of entropy yield and vortex core strength, namely, the impeller is optimally designed by simultaneously considering the efficiency of the impeller and pressure pulsation, so that the impeller can be ensured to have better operating efficiency and operating stability.
2) The method adopts the Q method to extract the surface area of the vortex structure of the main vortex core area in the impeller of the centrifugal pump, and digitizes the vortex structure, thereby being beneficial to vividly understanding the change of pressure pulsation and providing a basis for the optimization of the impeller.
3) By adopting an energy entropy method and according to the distribution condition of energy loss in the centrifugal pump, the main area of energy loss of the impeller can be analyzed, and the reference is provided for adjusting the design parameters of the impeller corresponding to the low-efficiency area of the impeller.
4) Based on the Q method and the energy entropy theory combined constraint, the efficient working condition point of the centrifugal pump impeller is optimized and determined, so that the optimized centrifugal pump impeller parameters can be obtained more efficiently, and the optimization is more accurate and efficient than the traditional impeller.
5) The design parameters of the centrifugal pump impeller are optimally designed by adopting a simulated annealing algorithm, so that the optimal design parameters of the impeller can be more accurately and quickly obtained.
Drawings
Fig. 1 is a flow chart of an optimal design method of an impeller of a centrifugal pump according to the present invention.
Fig. 2 is a schematic view of the inlet and outlet angles of the impeller of the initial centrifugal pump.
Fig. 3 is a view from a-a of fig. 2 (which is a schematic view of the width of the inlet and outlet of the impeller of the initial centrifugal pump).
FIG. 4 is a grid diagram of the optimized centrifugal pump impeller.
FIG. 5 is a graph of total entropy yield and vortex core strength as a function of flow rate for an impeller of a centrifugal pump according to an embodiment of the present invention under initial parameters.
FIG. 6 is a graph of total entropy yield and vortex core strength as a function of flow for an impeller of a centrifugal pump according to an embodiment of the present invention under optimum design parameters.
Detailed Description
The following further description is made with reference to the embodiments shown in the drawings.
The method for optimally designing the centrifugal pump impeller shown in FIG. 1 is carried out according to the following steps:
1) firstly, according to design requirements, a group of parameters to be optimized is specifically determined according to design experience: width of impeller inlet b1Impeller exit width b2Blade inlet angle beta1Blade exit angle beta2The number of blades z (angles and port widths are indicated in fig. 2 and 3)) This set of parameters that need to be optimized and the constraint ranges given below. Other impeller parameters are used as design constants;
the design parameters of the centrifugal pump impeller to be optimized are constrained as follows:
Figure GDA0002778916230000061
in the above formula:
Km1-the impeller inlet axial surface coefficient,
ns-the specific speed of rotation of the pump,
Fr-a flow rate of the fluid,
h-the lift is set as the height of the pump,
n-the rotating speed of the impeller,
g-acceleration of gravity;
when the design parameters of the centrifugal pump impeller are optimized in the step, the distribution condition of the energy loss inside the centrifugal pump in the step 4) can be analyzed, so that a main area of the energy loss of the impeller is obtained, and a reference is provided for the adjustment of the design parameters of the impeller corresponding to the low-efficiency area of the impeller, so that the pertinence of the optimization of the design parameters is improved, and the verification times and the workload of the design parameters are reduced.
2) Modeling and meshing: modeling a calculation domain of the centrifugal pump impeller under a determined coordinate system (shown in figure 2; the Y axis of the coordinate system is vertical to the paper surface and faces outwards) by using three-dimensional modeling software according to basic parameters of the centrifugal pump impeller; the three-dimensional model is gridded using ANSYS ICEM software.
3) And (3) CFX simulation verification: and (3) introducing the grids into CFX software, and performing transient numerical simulation calculation on the centrifugal pump impeller under the full-flow working condition under the given boundary conditions including inlet pressure, outlet flow, wall surface and dynamic and static interface settings and impeller rotating speed, wherein the results are automatically generated by the software. Then carrying out post-processing on the numerical simulation result to obtain performance curves of the centrifugal pump, such as flow-lift, flow-efficiency and the like; comparing whether the lift and the efficiency under the drawn grid meet the design requirements or not according to the performance curve, and if the design requirements are not met, repeating the step 2) to perform grid division on the three-dimensional model again until the grid division meeting the design requirements is obtained;
4) vortex structure extraction: based on the simulation result obtained in the step 3), vortex structure extraction is carried out in a calculation domain based on a Q method to obtain data of vortex core strength changing along with flow, so as to realize numerical processing of the vortex core strength and visually describe the change of pressure pulsation; wherein, the Q value represents the strength relation between the rotation and deformation of the fluid infinitesimal;
the method for extracting the vortex structure comprises the following specific steps: determining the Q value according to the formulas (1.5), (1.6) and (1.7):
Figure GDA0002778916230000081
Figure GDA0002778916230000082
Figure GDA0002778916230000083
the above equation is simplified to the following Q-values in cartesian coordinates:
Figure GDA0002778916230000084
wherein, the Q method defines the area where Q > 0 as the area where the vortex tube is; the deformation rate tensor and the vorticity in the above formula are calculated by editing a CEL expression in POST-processing software CFD-POST according to a numerical simulation result. Setting a threshold value, and extracting Fr=FrdThe infinitesimal iso-surface of (A) determines the vortex morphology, wherein FrdTo design the flow rate. After the Q value is determined, in CFD-POST, a vorticity isosurface under the Q value is made, the vorticity coefficient is adopted for coloring, then the strength of the vortex core is quantitatively calculated by software, and data are recorded.
5) Entropy production analysis: based on the grid obtained in the step 3), calculating the distribution condition of the internal energy loss of the centrifugal pump according to the speed distribution of the CFD calculation result according to the formulas (1.1), (1.2) and (1.3) based on the energy entropy theory; the total entropy yield of the impeller is then obtained from equation (1.4) to reflect the variation in impeller efficiency.
Figure GDA0002778916230000085
Figure GDA0002778916230000091
In the above formula
Spro,t-total entropy production rate per unit volume of impeller,
Figure GDA0002778916230000092
-mean entropy production rate per unit volume of impeller,
Spro,d′-a rate of impeller pulsation entropy production per unit volume,
t-the temperature of the fluid, and,
p-the density of the fluid,
-a turbulent dissipation ratio of the fluid,
u, v, w-velocity components of the fluid in the direction X, Y, Z respectively,
Figure GDA0002778916230000093
representing the time-average velocity value, u ', v ', w ' representing the pulsation velocity value;
6) efficient point matching: carrying out non-dimensionalization on the data obtained in the step 4) and the step 5), and introducing the data into data processing software for drawing to obtain curves of the total entropy yield and the vortex core strength of the impeller along with the change of flow; wherein the non-dimensional flow corresponding to the minimum vortex core strength is Fr1The dimensionless flow corresponding to the minimum total entropy yield of the impeller is Fr2
In the step, the total entropy productivity of the impeller is subjected to dimensionless treatment, and the numerical value of the entropy productivity of each impeller is divided by the minimum value of the entropy productivity of the impellerA value; carrying out non-dimensionalization treatment on the strength of the vortex core, namely dividing the numerical value of the strength of each vortex core by the minimum value of the strength of the vortex core; the flow is subjected to non-dimensionalization treatment by dividing the numerical value of each flow by the design flow Frd
7) Optimizing design parameters of the impeller by adopting a simulated annealing optimization algorithm to ensure that the dimensionless flow corresponding to the minimum vortex core strength and the dimensionless flow corresponding to the minimum entropy yield are both F'r(F′rTo design the flow FrdDimensionless flow obtained after dimensionless treatment, hence F'r1), i.e. F'r=Fr1=Fr2At this point, the optimum design parameters for the impeller are found to be 1.
Example (b): the basic design parameters of a group of centrifugal pump impellers are empirically given as initial design parameters, and other parameters are constant values. The initial design parameters are shown in table 1.
TABLE 1
Figure GDA0002778916230000101
According to the initial design parameters, the specific implementation operation is carried out according to the following steps.
1. Establishing a specific three-dimensional model of the impeller by using three-dimensional modeling software SolidWorks and parameters in table 1, as shown in FIG. 2;
2. the model is subjected to grid division, the total number of grid nodes is 2556414, and the grid distribution is shown in fig. 4 (the case of optimal design parameters). And the boundary conditions shown in Table 2 are given for numerical simulation calculation;
TABLE 2
Figure GDA0002778916230000102
In the step, grid verification is required, namely, whether the lift and the efficiency obtained by comparison and calculation meet the design requirements or not is compared, if the lift and the efficiency do not meet the requirements, the grid is modified until the lift and the efficiency meet the design requirements;
3. after the grid meets the requirements, the grid quality is proved to be appropriate, the accuracy of the calculation result is reliable, and the extraction of the vortex core structure is carried out at the moment. Namely, the Q method is used to obtain the data of the vortex core intensity varying with the flow rate, the calculated vortex core intensity is 0.28756 at the flow rate, and then the flow rate is changed to continue the calculation. Finally obtaining a curve 2 shown in figure 5 after dimensionless processing;
4. according to the numerical simulation result obtained in the step 2, obtaining the relationship of the total entropy yield of the impeller along with the change of the flow according to the formulas (1.1), (1.2), (1.3) and (1.4), such as a curve 1 shown in fig. 5;
5. the flow rate corresponding to the minimum vortex core intensity on curve 1 in FIG. 5 is set as Fr1The flow value corresponding to the minimum value of the total entropy yield of the impeller on the curve 2 is set as Fr2. At this time, the dimensionless flow rate corresponding to the minimum vortex core intensity is Fr1When the total entropy production of the impeller is minimum, the corresponding dimensionless flow is Fr20.76. At this time Fr2>Fr1Therefore, parameter optimization is required.
6. Optimizing design parameters of the impeller by adopting a simulated annealing optimization algorithm and a simulated annealing optimization algorithm to ensure that dimensionless flow corresponding to the minimum vortex core strength and dimensionless flow corresponding to the minimum entropy yield are both F'r,F′r=Fr1=Fr2At this point, the optimum design parameters for the impeller are found to be 1.
After the centrifugal pump impeller is subjected to multiple times of optimization of design parameters (steps 1-6 above), the optimal design parameters (shown in table 3) are finally obtained, and a vortex core strength-flow curve and an impeller total entropy yield-flow curve of the centrifugal pump impeller are drawn again (see fig. 5), wherein at this time, the dimensionless flow corresponding to the minimum vortex core strength and the dimensionless flow corresponding to the minimum entropy yield both reach F'r,F′r=Fr1=Fr2The design parameter can be determined as the optimal design parameter of the impeller, namely 1.
TABLE 3
Figure GDA0002778916230000111

Claims (4)

1. An optimal design method of a centrifugal pump impeller is characterized by comprising the following steps: the method comprises the following steps:
1) firstly, according to design requirements, a group of parameters to be optimized is specifically determined according to design experience: width of impeller inlet b1Impeller exit width b2Blade inlet angle beta1Blade exit angle beta2The number z of the blades, the group of parameters needing to be optimized and the constraint range are given below, and other impeller parameters are used as design constants;
the constraint range of the parameters to be optimized is as follows:
Figure FDA0002778916220000011
30°<β1<40°
15°<β2<40°
Figure FDA0002778916220000012
in the above formula
Km1-the impeller inlet axial surface coefficient,
ns-the specific speed of rotation of the pump,
Fr-a flow rate of the fluid,
h-the lift is set as the height of the pump,
n-the rotating speed of the impeller,
g-acceleration of gravity;
2) modeling and meshing: modeling a calculation domain of the centrifugal pump impeller under the determined coordinate system by using three-dimensional modeling software according to parameters of the centrifugal pump impeller; utilizing ANSYS ICEM software to perform meshing on the three-dimensional model;
3) and (3) CFX simulation verification: importing the grids into CFX software, carrying out transient numerical simulation calculation on a centrifugal pump impeller under a full-flow working condition under a given boundary condition, and automatically generating a result by the software; carrying out post-processing on the numerical simulation result to obtain a flow-lift and flow-efficiency performance curve of the centrifugal pump; comparing whether the lift and the efficiency under the drawn grid meet the design requirements or not according to the performance curve, and if the design requirements are not met, repeating the step 2) to divide the grid of the three-dimensional model again until the grid meeting the design requirements is obtained;
4) extracting a vortex core structure: based on the simulation result obtained in the step 3), extracting the vortex core structure in the calculation domain based on a Q method to obtain the data of the vortex core strength changing along with the flow, wherein the Q value represents the strength relation between the fluid infinitesimal rotation and the deformation;
5) entropy production analysis: based on the grid obtained in the step 3), calculating the distribution condition of the internal energy loss of the centrifugal pump according to the speed distribution of the CFD calculation result according to the formulas (1.1), (1.2) and (1.3); then carrying out volume integration on the formulas (1.2) and (1.3) in an impeller region by using a formula (1.4) to obtain the total entropy yield of the impeller;
Figure FDA0002778916220000021
Figure FDA0002778916220000022
Figure FDA0002778916220000023
Figure FDA0002778916220000024
in the above formula
Spro,t-total entropy production rate per unit volume of impeller,
Figure FDA0002778916220000025
-mean entropy production rate per unit volume of impeller,
Spro,d′-a rate of impeller pulsation entropy production per unit volume,
t-the temperature of the fluid, and,
p-the density of the fluid,
-a turbulent dissipation ratio of the fluid,
u, v, w-velocity components of the fluid in the direction X, Y, Z respectively,
Figure FDA0002778916220000026
representing the time-average velocity value, u ', v ', w ' representing the pulsation velocity value;
6) efficient point matching: carrying out non-dimensionalization on the data obtained in the step 4) and the step 5), and introducing the data into data processing software for drawing to obtain curves of the total entropy yield and the vortex core strength of the impeller along with the change of flow; wherein the non-dimensional flow corresponding to the minimum vortex core strength is Fr1The dimensionless flow corresponding to the minimum total entropy yield of the impeller is Fr2
7) Optimizing design parameters of the impeller by adopting a simulated annealing optimization algorithm to ensure that the dimensionless flow corresponding to the minimum vortex core strength and the dimensionless flow corresponding to the minimum impeller total entropy yield are both F'rAnd F'r=Fr1=Fr2At this point, the optimum design parameters for the impeller are found to be 1.
2. The optimum design method of a centrifugal pump impeller according to claim 1, characterized in that: in the step 4), the step of extracting the vortex structure by using the Q method comprises the following steps: determining the Q value according to the formulas (1.5), (1.6) and (1.7):
tensor of deformation ratio-
Figure FDA0002778916220000031
Vorticity-
Figure FDA0002778916220000032
Q value-
Figure FDA0002778916220000033
The above equation is simplified to the following Q-values in cartesian coordinates:
Figure FDA0002778916220000034
the Q method defines the area with Q > 0 as the area where the vortex tube is located, and the deformation rate tensor and the vorticity in the above formula are calculated by editing a CEL expression in POST-processing software CFD-POST according to a numerical simulation result; extraction of Fr=FrdThe infinitesimal iso-surface of (A) determines the vortex morphology, wherein FrdTo design the flow rate; after the Q value is determined, in CFD-POST, a vorticity isosurface under the Q value is made, a vorticity intensity coefficient is adopted for coloring, then the strength of a vortex core is quantitatively calculated by software, and data are recorded.
3. The optimum design method of a centrifugal pump impeller according to claim 1, characterized in that: in the step 5), the total entropy yield of the impeller is subjected to dimensionless treatment, wherein the numerical value of the entropy yield of each impeller is divided by the minimum value of the entropy yield of the impeller; carrying out non-dimensionalization treatment on the strength of the vortex core, namely dividing the numerical value of the strength of each vortex core by the minimum value of the strength of the vortex core; the flow is subjected to non-dimensionalization treatment by dividing the numerical value of each flow by the design flow Frd
4. The method for optimally designing an impeller of a centrifugal pump according to claim 1, wherein the method comprises the steps of: in the step 3), the boundary conditions comprise inlet pressure, outlet flow, wall surface and dynamic and static interface settings and impeller rotating speed.
CN201910914770.XA 2019-09-26 2019-09-26 Optimization design method of centrifugal pump impeller Active CN110617238B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910914770.XA CN110617238B (en) 2019-09-26 2019-09-26 Optimization design method of centrifugal pump impeller

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910914770.XA CN110617238B (en) 2019-09-26 2019-09-26 Optimization design method of centrifugal pump impeller

Publications (2)

Publication Number Publication Date
CN110617238A CN110617238A (en) 2019-12-27
CN110617238B true CN110617238B (en) 2021-01-05

Family

ID=68924555

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910914770.XA Active CN110617238B (en) 2019-09-26 2019-09-26 Optimization design method of centrifugal pump impeller

Country Status (1)

Country Link
CN (1) CN110617238B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111523186B (en) * 2020-05-19 2024-01-19 重庆水泵厂有限责任公司 Optimization method of shape of water suction chamber for double water suction pump
CN112555451B (en) * 2020-10-29 2022-03-18 成都成高阀门有限公司 Pressing process parameter determination method for all-welded ball valve seat sealing groove and semi-finished valve seat
CN114692338B (en) * 2022-04-12 2024-04-26 浙江理工大学 Comprehensive optimization design method for cavitation and efficiency of low-temperature centrifugal pump
CN116186947B (en) * 2023-05-05 2023-08-04 陕西空天信息技术有限公司 Impeller simulation design method, device, equipment, medium and product

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003323468A (en) * 2002-05-01 2003-11-14 Tsurumi Mfg Co Ltd Method for designing balance of rigid rotor
JP2004185102A (en) * 2002-11-29 2004-07-02 Nippon Steel Corp Electromagnetic analysis device and method
CN105465037A (en) * 2015-12-01 2016-04-06 中国农业大学 Hydraulic power optimizing method and device for double-suction centrifugal pump impeller
CN105550404A (en) * 2015-12-07 2016-05-04 江苏大学 Method for analyzing hydraulic loss of centrifugal pump based on entropy theory
CN105604953A (en) * 2015-12-16 2016-05-25 江苏大学 Multistage-pump impeller stagger arranging method based on vibration optimization
CN106202795A (en) * 2016-07-21 2016-12-07 浙江理工大学 Centrifugal pump impeller and the method for designing thereof of constraint is combined based on entropy product and blade loading
CN206162530U (en) * 2016-07-21 2017-05-10 浙江理工大学 Centrifugal pump impeller based on entropy is produced and blade load allies oneself with contract bundle design method manufacturing
CN109960836A (en) * 2017-12-26 2019-07-02 中国科学院深圳先进技术研究院 A kind of method, apparatus and terminal device of turbulent wake emulation
CN109977345A (en) * 2019-01-29 2019-07-05 河海大学 A kind of method for numerical simulation of axial-flow pump clearance leakage of blade tip vortex cavitation

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003323468A (en) * 2002-05-01 2003-11-14 Tsurumi Mfg Co Ltd Method for designing balance of rigid rotor
JP2004185102A (en) * 2002-11-29 2004-07-02 Nippon Steel Corp Electromagnetic analysis device and method
CN105465037A (en) * 2015-12-01 2016-04-06 中国农业大学 Hydraulic power optimizing method and device for double-suction centrifugal pump impeller
CN105550404A (en) * 2015-12-07 2016-05-04 江苏大学 Method for analyzing hydraulic loss of centrifugal pump based on entropy theory
CN105604953A (en) * 2015-12-16 2016-05-25 江苏大学 Multistage-pump impeller stagger arranging method based on vibration optimization
CN106202795A (en) * 2016-07-21 2016-12-07 浙江理工大学 Centrifugal pump impeller and the method for designing thereof of constraint is combined based on entropy product and blade loading
CN206162530U (en) * 2016-07-21 2017-05-10 浙江理工大学 Centrifugal pump impeller based on entropy is produced and blade load allies oneself with contract bundle design method manufacturing
CN109960836A (en) * 2017-12-26 2019-07-02 中国科学院深圳先进技术研究院 A kind of method, apparatus and terminal device of turbulent wake emulation
CN109977345A (en) * 2019-01-29 2019-07-05 河海大学 A kind of method for numerical simulation of axial-flow pump clearance leakage of blade tip vortex cavitation

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
自吸泵内能量损失及非定常流动特性研究;舒欣,任芸,吴登昊,等;《水利学报》;20190831;1010-1019 *

Also Published As

Publication number Publication date
CN110617238A (en) 2019-12-27

Similar Documents

Publication Publication Date Title
CN110617238B (en) Optimization design method of centrifugal pump impeller
CN110909422B (en) Method for predicting and optimizing range of efficient working condition area of centrifugal pump impeller
Yang et al. Energy characteristics and optimal design of diffuser meridian in an electrical submersible pump
Pei et al. Multi-point optimization on meridional shape of a centrifugal pump impeller for performance improvement
Li et al. Numerical investigation of impeller trimming effect on performance of an axial flow fan
Zhou et al. CFD investigation and PIV validation of flow field in a compact return diffuser under strong part-load conditions
CN109578085B (en) Method for weakening unsteady acting force of turbine movable blade through guide blade inclination
CN112784370B (en) Design method of space guide vane of multistage centrifugal pump
Pei et al. Statistical analysis of pressure fluctuations during unsteady flow for low-specific-speed centrifugal pumps
Bai et al. Effect of blade wrap angle in hydraulic turbine with forward-curved blades
JP5770993B2 (en) A method for predicting and evaluating the stability of cavitation behavior of inducers or impellers
CN111832132B (en) Design method of hydraulic model of low-specific-speed high-speed centrifugal pump
Kim et al. Design optimization of mixed-flow pump in a fixed meridional shape
Luo et al. Investigation on the effect of the shaft transition form on the inflow pattern and hydrodynamic characteristics of the pre-shaft tubular pump device
Peng et al. Optimization design of multistage pump impeller based on response surface methodology
Lin et al. Research on the rotor-stator interaction of centrifugal pump based on sinusoidal tubercle volute tongue
CN114692338B (en) Comprehensive optimization design method for cavitation and efficiency of low-temperature centrifugal pump
Yi-bin et al. The numerical simulation of unsteady flow in a mixed flow pump guide vane
Pan et al. A comprehensive assessment on detached eddy simulation and delayed detached eddy simulation approaches in large-scale rotating separated flow in a centrifugal impeller under small flow rate
CN112761969A (en) Optimization method for structural parameters of centrifugal pump
Mahdi et al. Unsteady fluid flow analysis of tongue geometry in a centrifugal pump at design and off-design conditions
Wang et al. The numerical simulation of full flow field of Roto-Jet pump and analysis of energy losses
Weidong et al. Optimization design of stainless steel stamping multistage pump based on orthogonal test
Ali et al. Numerical analysis of flow phenomena in a centrifugal pump impeller of low specific speed
Wang et al. Comparative study on hydraulic performance of single-suction centrifugal pump and double-suction centrifugal pump

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant