CN106321465A - Kinetics characterization method for rotating stall of centrifugal pump impeller - Google Patents
Kinetics characterization method for rotating stall of centrifugal pump impeller Download PDFInfo
- Publication number
- CN106321465A CN106321465A CN201610705946.7A CN201610705946A CN106321465A CN 106321465 A CN106321465 A CN 106321465A CN 201610705946 A CN201610705946 A CN 201610705946A CN 106321465 A CN106321465 A CN 106321465A
- Authority
- CN
- China
- Prior art keywords
- rotating stall
- impeller
- centrifugal pump
- fractal
- stall
- 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.)
- Pending
Links
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04D—NON-POSITIVE-DISPLACEMENT PUMPS
- F04D15/00—Control, e.g. regulation, of pumps, pumping installations or systems
- F04D15/0088—Testing machines
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04D—NON-POSITIVE-DISPLACEMENT PUMPS
- F04D29/00—Details, component parts, or accessories
- F04D29/18—Rotors
- F04D29/22—Rotors specially for centrifugal pumps
- F04D29/2205—Conventional flow pattern
- F04D29/2216—Shape, geometry
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
- Control Of Non-Positive-Displacement Pumps (AREA)
Abstract
The invention discloses a kinetics characterization method for rotating stall of a centrifugal pump impeller. The kinetics characterization method mainly comprises the following steps: performing numerical simulation or experimental measurement on a centrifugal pump to obtain a pressure time series of a point at the outlet of the impeller, and re-building a low-dimensional dynamic system with kinetics properties, that is, performing phase-space reconstruction on the pressure time series to rebuild the low-dimensional dynamic system; and analyzing the fractal characteristics of rotating stall, and measuring the flow complexity after rotating stall by use of a given fractal dimension. According to the method, based on numerical simulation of flow before and after rotating stall of the centrifugal pump impeller, the phase-space reconstruction study method is used in combination with the C-C algorithm to study the pressure time series to obtain that some fractal structures exist in the rotating stall state, and the fractal dimensions of the fractal structures can more accurately capture the dynamics characteristics of rotating stall. The quantitative changes can be effectively used for characterizing the flow complexities at different positions near the outlet of the impeller after rotating stall.
Description
Technical field
The present invention relates to fractal dynamics field in rotating machinery phase space reconfiguration and nonlinear kinetics, especially one
Plant the Kinetic Characterization method of centrifugal pump impeller rotating stall.
Background technology
When centrifugal pump runs under low flow rate condition, than design conditions, great changes have taken place in internal flow, often occurs losing
The Complex Flows such as speed, Secondary Flow, import and export backflow.When particularly flow is reduced to certain limit, due to the flow inlet angle of attack relatively
Greatly, blade suction side flow separation increases, and forms the vortex of clogging impeller runner, and impeller enters stall conditions, blocks accordingly
The vortex of impeller passage is referred to as stall group.In rotating stall, stall group not only can destroy the uniformity of interior flow field, also can produce
Raw extra dynamic loading, can cause the fatigue at blade high stress point, fracture time serious.When stall frequencies and centrifugal pump group
When frequency is consistent, even can cause the vibration of whole unit.The existence of these flow instabilities phenomenons not only can cause whole machine
Group efficiency declines, it is also possible to make centrifugal pump operation stability and safety be affected by serious but, the efficiency of centrifugal pump and
Stablize and limited by rotating stall.That is, by control rotating stall be improve centrifugal pump internal flow power performance
A subject matter.Therefore, the dynamic characteristic of research rotating stall rotates for the flowing controlled within compressor impeller and loses
Speed is critically important.At present, the spectroscopic analysis methods for rotating stall pressure fluctuation is the main method of extraction system feature.
Specifically, the stall of centrifugal pump and flow separation can affect blade tip clearance by numerical simulation study, the method that liquid flows into
The flowing at place, thus control rotating stall.
Centrifugal pump internal flow unstability is studied by other research worker, but few people pass through dynamics research
Rotating stall.Generally, the rotating stall of centrifugal pump is to be obtained by the time series analyzing dynamical system, but hides
Kinetic character in systems can not be directly displayed by the time series studying these state variables, therefore loses for rotating
The description of speed is not enough.
Seasonal effect in time series phase space reconfiguration proposed in 1980, and then the method combines other technologies by other research worker
Being applied to more field, such as, chaos attractor is used to predict the rotating stall of compressor.Research is also had to compare stall
The correlation dimension of pressure signal front and back, relevant dimension is counted as the internal characteristics of fan stall simultaneously.
Chaotic Time Series Analysis is the important development of Nonlinear Time Series Analysis, and it has been widely used in
The directions such as biomedicine, meteorology, finance.Phase space reconfiguration is as the committed step of Chaotic Time Series Analysis, and it embeds
The computational methods of dimension m and delay time T are in constantly development.1986, Broomhead and King in Practical Calculation first
Propose selected (m-1) τ-value, while increasing m, reduce the value (keeping (m-1) τ is constant) of τ, choose optimal m and τ
The method of value.1996, D.Kugiumtzis proposed time delay τdChoose should be independent of Embedded dimensions m, and should
Studied by other research worker according to centrifugal pump internal flow unstability, but few people are rotated by dynamics research and lose
Speed.Generally, the rotating stall of centrifugal pump is to be obtained by the time series analyzing dynamical system, but is hidden in system
In kinetic character can not be directly displayed by the time series studying these state variables.Rely in time delay window τ w=
(m-1) τ, and to obtain τ w >=τ p, τ p here according to experiment be the mean orbit cycle, can be estimated by seasonal effect in time series fluctuation
Out.But above method or be readily obtained the conclusion of mistake, or amount of calculation is excessive, it is not easy to implement.1999, H.S.Kim,
R.Eykholt and J.D.Salas proposes C-C algorithm, its compared with before method have calculate accurately, feature that amount of calculation is little.
C-C algorithm association integration is able to rely on m；And time window τ w depends on m, and τ changes with m.
Scholar is not yet had to utilize the method for phase space reconfiguration to be centrifuged pump in nonlinear dynamics theory at present
Kinetic Characterization.
Summary of the invention
For solving above-mentioned technical problem, the present invention provides a kind of Kinetic Characterization method of centrifugal pump impeller rotating stall,
Utilizing phase space reconfiguration in nonlinear dynamics theory, dynamic (dynamical) feature can obtain from univariate time series,
Thus be used for characterizing rotating stall.
The technical scheme that this method uses is as follows:
Step one: the method utilizing numerical simulation or experiment to measure, changes the flow of centrifugal pump until rotating mistake
The phenomenon of speed, obtains the pressure time sequence of monitoring point, impeller a certain runner exit position one during rotating stall；
Step 2: phase space reconfiguration is applied to the pressure time sequence of this point, obtain reconstruct comprises kinetic property
Low cone motive system；
Step 3: utilize low cone motive system to analyze the fractal property of this rotating stall further, with fractal structure
Complex flow after Analysis of Fractal Dimension rotating stall.
Monitoring site in step one in a certain runner of impeller, with impeller outlet distance be impeller radius 2%～
5%, with back shroud distance be impeller outlet width 30%～70%, in the runner middle position of adjacent blades.
In step 2 during utilizing time delays method to carry out phase space reconfiguration, C-C algorithm is for the crucial ginseng obtained
Number, such as time delay and the Embedded dimensions of phase space reconstruction.Phase space reconfiguration is utilized, from monotropic in nonlinear dynamics theory
The time series of amount obtains dynamic (dynamical) feature, thus analyzes rotating stall.
Beneficial effects of the present invention: measure, with combining C-C based on to the numerical simulation of centrifugal impeller rotating stall or experiment
The technique study of the phase space reconfiguration research of algorithm obtains pressure time sequence.There are some fractal knots under rotating stall state
Structure, the fractal dimension of fractal structure can catch the dynamic characteristic of rotating stall more accurately.The fractal dimension of diverse location
Showing, near impeller outlet with Radius, owing to impeller outlet is with Similarity of Flow at Radius, fractal dimension is also
Similar, position is further away from impeller outlet, and fractal dimension is the fewest.These quantitative changes can be used for characterizing stall posterior lobe effectively
Wheel near exit various location complex flow.
Accompanying drawing explanation
Fig. 1 is the method step flow chart of the present invention.
Fig. 2 is the time series chart under stall conditions.
Fig. 3 is the seasonal effect in time series phase space reconfiguration under steady statue and the seasonal effect in time series phase space weight under stall conditions
The comparison diagram of structure.
Fig. 4 is the correlation dimension figure of phase-space reconstruction under stall conditions.
Detailed description of the invention
It is an object of the invention to provide a kind of Kinetic Characterization method of centrifugal pump impeller rotating stall, at Nonlinear Dynamic
Theory of mechanics utilizes phase space reconfiguration, dynamic (dynamical) feature can be obtained from univariate time series, thus analyze rotation
Turn stall.It appeared that fractal structure from the phase space of reconstruct, fractal dimension can be used to describe answering of parting structure simultaneously
Polygamy.In some sense, fractal structure is considered as the dynamic characteristic of a kind of centrifugal impeller rotating stall, and is conducive to
Characterize and analyze rotating stall.
Below as a example by a certain specific embodiment of centrifugal pump stall, in conjunction with the flow process of Fig. 1, selected specific design parameter is right
The present invention is further described:
Step one: utilize the means that numerical simulation or experiment are measured, when changing flow to 100kg/s, this centrifugal pump leaf
Wheel rotates the phenomenon of stall, obtains the pressure time sequence of monitoring point, impeller a certain runner exit position one during rotating stall
Row (see Fig. 2), this monitoring site in a certain runner of impeller, with impeller outlet distance be impeller radius 2%～5%, with rear
30%～70% that cover plate distance is impeller outlet width, in adjacent blades middle position.
Step 2: phase space reconfiguration is applied to these pressure time sequences, and reconstruct one and comprise kinetic property
Low cone motive system, during phase space reconfiguration, C-C algorithm for obtain key parameter.Study the science at non-linear dynamic
Utilize phase space reconfiguration in Lun, from univariate time series, obtain dynamic (dynamical) feature, thus analyze rotating stall.
As in figure 2 it is shown, during to 3000 step, already at steady statue, therefore (time sequence can be taken in conjunction with C-C algorithm
Row length N=3000, Embedded dimensions m=2, radius ri=i σ/2, i=1,2,3,4) optimal time obtained postpones τ=9, time
Between delay aperture τW=11, according to formula τW=(m-1) τ obtains Embedded dimensions m=2.22.Again pressure time sequence is carried out empty mutually
Between rebuild, obtain a series of fractal structure.
As it is shown on figure 3, the seasonal effect in time series phase space reconfiguration under steady statue and the seasonal effect in time series phase under stall conditions
Space Reconstruction is very different, and has the most similar figure in stall figure, in nonlinear kinetics, and the figure of these self similarities
Shape is fractal structure, and it is one of dynamic characteristic of rotating stall, and after display stall, the time series of impeller outlet is not
Unordered, can be used to characterize and analyze rotating stall.
Step 3: low cone motive system is made up of fractal structure, utilizes low cone motive system to analyze this point further and rotates
The fractal property of stall, with the correlation dimension d in the fractal dimension of fractal structure2Sign rotating stall characteristic such as Fig. 4:
T is the time delay of phase-space reconstruction.
Claims (4)
1. a Kinetic Characterization method for centrifugal pump impeller rotating stall, comprises the steps:
Step one: measured the means of rotating stall by numerical simulation or experiment, obtain during rotating stall in the arbitrary runner of impeller
The pressure time sequence of the monitoring point at exit position；
Step 2: phase space reconfiguration is applied to the pressure time sequence of this point, obtain reconstruct comprises the low of kinetic property
Cone motive system；
Step 3: utilize low cone motive system to analyze the fractal property of this rotating stall further, fractal with fractal structure
Complex flow after dimensional analysis rotating stall.
The Kinetic Characterization method of a kind of centrifugal pump impeller rotating stall the most according to claim 1, it is characterised in that step
Monitoring site described in rapid one in the arbitrary runner of impeller, with impeller outlet distance be impeller radius 2% ~ 5%, with impeller bonnet
30% ~ 70% that plate distance is impeller outlet width, in adjacent blades middle position.
The Kinetic Characterization method of a kind of centrifugal pump impeller rotating stall the most according to claim 1, step 2 is specially
The method using time delays carries out phase space reconfiguration, wherein utilizes C-C algorithm to obtain the embedding of time delay and phase space reconstruction
Enter dimension.
The Kinetic Characterization method of a kind of centrifugal pump impeller rotating stall the most according to claim 1, step 2 is specially
In nonlinear dynamics theory, utilize phase space reconfiguration, from univariate pressure time sequence, obtain fractal structure, profit
By the fractal dimension phenetic analysis rotating stall of fractal structure.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610705946.7A CN106321465A (en) | 2016-08-22 | 2016-08-22 | Kinetics characterization method for rotating stall of centrifugal pump impeller |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610705946.7A CN106321465A (en) | 2016-08-22 | 2016-08-22 | Kinetics characterization method for rotating stall of centrifugal pump impeller |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106321465A true CN106321465A (en) | 2017-01-11 |
Family
ID=57741997
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610705946.7A Pending CN106321465A (en) | 2016-08-22 | 2016-08-22 | Kinetics characterization method for rotating stall of centrifugal pump impeller |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106321465A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109030849A (en) * | 2018-05-22 | 2018-12-18 | 大连透平机械技术发展有限公司 | The monitoring method and device of compressor rotating stall |
| CN112417774A (en) * | 2020-10-14 | 2021-02-26 | 沈阳透平机械股份有限公司 | Method and device for identifying rotating stall of compressor and computer equipment |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH06101695A (en) * | 1992-05-28 | 1994-04-12 | Nikkiso Co Ltd | Axial diffuser and method for running centrifugal pump provided with axial diffuser |
| EP0571739B1 (en) * | 1992-05-28 | 1997-10-15 | Nikkiso Co., Ltd. | Axial diffusor for a centrifugal pump |
| CN1675470A (en) * | 2002-08-23 | 2005-09-28 | 约克国际公司 | System and method for detecting rotating stall in a centrifugal compressor |
| WO2009041851A1 (en) * | 2007-09-24 | 2009-04-02 | Central Institute Of Aviation Motors (Ciam) | Method for monitoring the operating modes of a compressor and a device for carrying out said method |
| CN202266441U (en) * | 2011-03-03 | 2012-06-06 | 刘滨军 | Rotary drum type vortex pump, slit taking device and slit dredger |
| CN104154043A (en) * | 2014-08-07 | 2014-11-19 | 江苏大学 | Axial flow pump without blade top leakage vortex |
| CN105805046A (en) * | 2016-01-25 | 2016-07-27 | 江苏大学 | Device for improving internal flow stability of diagonal flow pump |
-
2016
- 2016-08-22 CN CN201610705946.7A patent/CN106321465A/en active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH06101695A (en) * | 1992-05-28 | 1994-04-12 | Nikkiso Co Ltd | Axial diffuser and method for running centrifugal pump provided with axial diffuser |
| EP0571739B1 (en) * | 1992-05-28 | 1997-10-15 | Nikkiso Co., Ltd. | Axial diffusor for a centrifugal pump |
| CN1675470A (en) * | 2002-08-23 | 2005-09-28 | 约克国际公司 | System and method for detecting rotating stall in a centrifugal compressor |
| WO2009041851A1 (en) * | 2007-09-24 | 2009-04-02 | Central Institute Of Aviation Motors (Ciam) | Method for monitoring the operating modes of a compressor and a device for carrying out said method |
| CN202266441U (en) * | 2011-03-03 | 2012-06-06 | 刘滨军 | Rotary drum type vortex pump, slit taking device and slit dredger |
| CN104154043A (en) * | 2014-08-07 | 2014-11-19 | 江苏大学 | Axial flow pump without blade top leakage vortex |
| CN105805046A (en) * | 2016-01-25 | 2016-07-27 | 江苏大学 | Device for improving internal flow stability of diagonal flow pump |
Non-Patent Citations (1)
| Title |
|---|
| 王乐等: "《离心压缩机叶轮旋转失速的相空间重构及分形特征》", 《航空动力学报》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109030849A (en) * | 2018-05-22 | 2018-12-18 | 大连透平机械技术发展有限公司 | The monitoring method and device of compressor rotating stall |
| CN109030849B (en) * | 2018-05-22 | 2021-06-11 | 大连透平机械技术发展有限公司 | Method and device for monitoring rotating stall of compressor |
| CN112417774A (en) * | 2020-10-14 | 2021-02-26 | 沈阳透平机械股份有限公司 | Method and device for identifying rotating stall of compressor and computer equipment |
| CN112417774B (en) * | 2020-10-14 | 2024-02-23 | 沈阳透平机械股份有限公司 | Identification method and device for rotating stall of compressor and computer equipment |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Semlitsch et al. | Flow phenomena leading to surge in a centrifugal compressor | |
| Wang et al. | Application of the two-phase three-component computational model to predict cavitating flow in a centrifugal pump and its validation | |
| CN109185211B (en) | A kind of pumping plant unit pressure pulsation prediction method instantaneously calculated based on Fluid Mechanics Computation | |
| Liu et al. | Experimental investigation and numerical analysis of unsteady attached sheetcavitating flows in a centrifugal pump | |
| Haldeman et al. | Experimental and computational investigation of the time-averaged and time-resolved pressure loading on a vaneless counter-rotating turbine | |
| Deskos et al. | Turbulence-resolving simulations of wind turbine wakes | |
| Dong et al. | Study on unstable characteristics of centrifugal pump under different cavitation stages | |
| CN106321465A (en) | Kinetics characterization method for rotating stall of centrifugal pump impeller | |
| Wang et al. | Applicability of eddy viscosity turbulence models in low specific speed centrifugal pump | |
| CN105889136B (en) | A kind of centrifugal pump improved method for being used for reducing the pulsation of centrifugal pump internal pressure | |
| Mao et al. | Applications of different turbulence models in simulations of a large annular volute-type pump with the diffuser | |
| Andersen et al. | Reduced order model of the inherent turbulence of wind turbine wakes inside an infinitely long row of turbines | |
| Xu et al. | Pressure and velocity fluctuations characteristics of the tip clearance flow in an axial compressor stage at the near-stall condition | |
| Zhou et al. | Hydraulic design and performance analysis of low specific speed centrifugal pump | |
| CN110259702A (en) | A kind of centrifugal pump operational effect determination method | |
| Angelini et al. | Identification of losses in turbomachinery with machine learning | |
| Gölcü | Neural network analysis of head-flow curves in deep well pumps | |
| CN106934134A (en) | A kind of method for characterizing the centrifugation unsteady intensity of pump cavitation | |
| Li et al. | Cavitation simulation and NPSH prediction of a double suction centrifugal pump | |
| Tan et al. | PIV Measurements of the Flow in the Tip Region of a Compressor Rotor | |
| Chenaux et al. | Aeroelastic investigation of an annular transonic compressor cascade: Numerical sensitivity study for validation purposes | |
| Tan et al. | Visualization and time resolved PIV measurements of the flow in the tip region of a subsonic compressor rotor | |
| Zhao et al. | Surge prediction in a gas-liquid mixed-flow pump based on multifractal detrended fluctuation analysis | |
| Zhang et al. | Characteristic analysis on the pressure fluctuation in the impeller of a low specific speed mixed flow pump | |
| Taghavi-Zenouz et al. | Multiple spike stall cells in low speed axial compressor rotor blade row |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170111 |
|
| WD01 | Invention patent application deemed withdrawn after publication |