CN107273570A - A kind of blade pump cavitation Inductive noise Numerical Prediction Method - Google Patents

A kind of blade pump cavitation Inductive noise Numerical Prediction Method Download PDF

Info

Publication number
CN107273570A
CN107273570A CN201710325748.2A CN201710325748A CN107273570A CN 107273570 A CN107273570 A CN 107273570A CN 201710325748 A CN201710325748 A CN 201710325748A CN 107273570 A CN107273570 A CN 107273570A
Authority
CN
China
Prior art keywords
cavitation
model
pump
file
sound
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
Application number
CN201710325748.2A
Other languages
Chinese (zh)
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.)
Jiangsu University
Original Assignee
Jiangsu University
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 Jiangsu University filed Critical Jiangsu University
Priority to CN201710325748.2A priority Critical patent/CN107273570A/en
Publication of CN107273570A publication Critical patent/CN107273570A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Structures Of Non-Positive Displacement Pumps (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The present invention provides a kind of blade pump cavitation Inductive noise Numerical Prediction Method, for calculating interior outfield radiated noise of the vane pump in low pressure or big flow occasion generation cavitation, including:Convection current field domain and pump case structure are modeled, and convection current field domain carries out mesh generation, check mesh quality, and according to pump carrying out practically parameter setting pressure inlets, speed export boundary condition sets solver parameter, and calculate the non-permanent file of cavitation;Carry out stationary iteration calculating and obtain cavitation numerical value destination file under each cavitation number, and model progress turbulence model and cavitation model are verified;The initial value that permanent destination file as the cavitation unsteady computation of cavitation model will be added carries out Non-stationary iterative method calculating, obtains under each cavitation number comprising pulsation and gridding information file;In acoustics software for calculation blade pump cavitation Inductive noise is calculated using the coupling algorithm of acoustic analogy theory and FInite Element and boundary element method.The present invention can relatively accurately carry out cavitation Inductive noise forecast.

Description

A kind of blade pump cavitation Inductive noise Numerical Prediction Method
Technical field
Patent of the present invention belongs to fluid machinery cavitation noise calculating field, is a kind of numerical value of blade pump cavitation Inductive noise Forecasting procedure.
Background technology
Vane pump, as a kind of transducing head applied to national economy every field, is to be related to national economy and country The grand strategy equipment of safety.The key issue core that cavitation develops as vane pump engineering field, restricts it towards highly reliable Property, the direction of high stability are developed.Due to the complexity of cavitation mechanism, numerical computation method, which gradually substitutes prototype experiment, turns into sky Change one of Main Means of research.
And the numerical simulation of cavitation is influenceed by multiphase flow model, turbulence model, cavitation model and every fluid parameter, because This not only needs to select suitable model for prediction blade pump cavitation Inductive noise, is made an uproar with greater need for the main frequency range concentrated to energy Sound Accurate Prediction.Research for cavitation noise at present focuses primarily upon single spherical bubble collapse and produces impulse radiation noise reason By prediction, this method can not well with real data contrast verification, for cavitation it is serious when hole movement cause pressure low Cavitation Inductive noise produced by frequency vibration is swung can not be calculated accurately.And as the high-frequency noise energy attenuation produced by bubble collapse very It hurry up, low-frequency noise energy transmission is stronger, it is bigger to vane pump reliability of operation and stability harm, therefore blade pump cavitation Inductive noise numerical prediction is a research puzzle.
Patent proposition of the present invention is a kind of can to calculate acoustic fluid domain and thin-walled by what boundary element method was coupled with FInite Element The method of acoustic power flow effect calculates blade pump cavitation Inductive noise between structure, solves the numerical prediction to cavitation Inductive noise and asks Topic, design and cavitation monitoring tool to altitude performance vane pump are of great significance.
The content of the invention
The purpose of the present invention:For problem and shortage present in current blade pump cavitation Forecasting Methodology, it is contemplated that A kind of Numerical Prediction Method this method of more accurate blade pump cavitation Inductive noise is provided, passes through CFD numerical computations software moulds Intend vane pump cavitation characterization and extract fixed dipole and rotating dipole sound source file, application boundary member method calculates cavitation induction Radiated sound field in noise, and application boundary member and Coupled with Finite Element algorithm calculate the outfield that acoustic power flow effect produces and radiate and make an uproar Sound, obtains internal field sound pressure level and radiated level curve, external radiation acoustical power frequency respective function and sound field under each cavitation number Point to distribution.
The technical scheme is that:A kind of blade pump cavitation Inductive noise Numerical Prediction Method includes:(1) convection current field domain And pump case structure is modeled, convection current field domain carries out mesh generation, mesh quality is checked, according to pump carrying out practically parameter setting Pressure inlets, speed export boundary condition sets solver parameter, and calculate the non-permanent file of cavitation.(2) stationary iteration is carried out Calculating obtains cavitation numerical value destination file under each cavitation number, and model progress turbulence model and cavitation model are verified.(3) The initial value that permanent destination file as the cavitation unsteady computation of cavitation model will be added carries out Non-stationary iterative method calculating, obtains Pulsation and gridding information file are included under each cavitation number.(4) acoustic analogy theory and FInite Element are applied in acoustics software for calculation And the coupling algorithm of boundary element method calculates blade pump cavitation Inductive noise.
A kind of blade pump cavitation Inductive noise Numerical Prediction Method, comprises the following steps:Step 1) build enclosed visualization Experimental bench, external characteristics and cavitation performance curve are obtained by testing, and obtain a cavitation stage internal field vacuole distribution, and threedimensional model is built Vertical and mesh generation:Convection current field domain and pump case structure are modeled, and convection current field domain carries out mesh generation, check mesh quality;
Step 2) non-cavitating numerical computations file is set using CFD software:By step 1) grid import CFD software, enter Row bound condition setting, turbulence model selection, solver parameter setting calculates non-cavitating numerical computations file;
Step 3) turbulence model is modified until reaching precision of prediction:By contrast test and calculated value external characteristics it Between difference, determine turbulence model if appropriate for this status condition, if difference is outside allowed band, return to step 2) reselect Turbulence model is modified to turbulence model until precision of prediction is reached, until difference is in allowed band, then into step 4);
Step 4) the permanent numerical computations file of cavitation is set:Gas-liquid two-phase physical parameter, setting initial component border are set Condition, sets inlet-pressure boundary condition, cavitation model selection;
Step 5) cavitation model is modified until reaching precision of prediction:By contrasting cavitation performance parameter and test value Difference, determines cavitation model if appropriate for operating mode instantly, the return to step 4 if not meeting) it is further by high speed photography test results Amendment gasification coefficient, condensation coefficient and density revision function optimization cavitation model are until reach precision of prediction, until meeting calculating Required precision enters step 6);
Step 6) unsteady cavitation numerical value calculation document is set:Fixed dipole and rotating dipole sound source, which are extracted, to be set, The permanent original document of cavitation is iterated to calculate into unsteady numerical value, cavitation form is obtained;
Step 7) internal field noise is calculated using LMS Virtual.Lab software BEM modules:Fixed and rotating dipole sound source Input, boundary element and site mess generation, carry out fluent material attribute definition, and boundary condition is set, and sound absorption attribute is set, and is started Sound field and site distribution are calculated;
Step 8) modal calculation is carried out using LMS Virtual.Lab softwares FEM model:BEM/FEM coupling calculating is obtained Outfield radiated noise, including FEM meshing, modal calculation, damping addition, coupled vibro-acoustic definition and calculating, obtain certain empty Change external radiation acoustical power frequency respective function and sound field under coefficient and point to distribution.
In such scheme, the step 1) three-dimension modeling and mesh generation specifically include following steps:
Step 11) application 3 d modeling software carries out three-dimensional modeling to pump fluid domain, and three-dimensional modeling is carried out to the pump housing, it is defeated Go out the model file that suffix is stp;
Step 12) mesh generation is carried out to fluid domain and housing respectively, and mesh quality is checked, to ensure computational accuracy, Adjust mesh quality and be more than 0.2, minimum grid angle is more than 18 degree, and exports the grid file that suffix is cfx.
In such scheme, the step 2) set non-cavitating numerical computations file to specifically include following step using CFD software Suddenly:
Step 21) by step 1) grid file import in CFD software and carry out computation model setting, select Isothermal heat transfer models and k-Epsilon turbulence models close Navier Stokes equation;Fluid media (medium) is chosen and root Factually actual value sets fluid media (medium) density and viscosity;Set and import and export boundary condition, environmental pressure is set according to actual value, Pump impeller rotating speed is set according to actual conditions, rotary speed unit is rev min^-1;Rotation computational fields wall is set to no sliding Wall, sets rotary shaft wall relative velocity consistent with impeller direction of rotation to be set to 0, uses rotating coordinate systems;Static meter Domain wall is calculated to be set to, without sliding wall, be set to global rest frame;Sound computational fields interface uses MRF Multiple references System, transformed coordinate system is set to Frozen Rotor, and Mesh connection mode is set to GGI;Computational fields interface connects mould silently Type is set to generic connectivity, and Mesh connection mode is set to GGI;Convective term is set to high resolution mode, turbulence number in solver parameter Value item is set to single order windward, solves step number and is set according to actual conditions, and it is 10 that convergence criterion residual error, which is set to average value RMS value, ^-4, is monitored using CEL language write-in lift efficiency expression formula, the non-cavitating numerical computations file that output suffix is def;
Step 22) non-cavitating file calculates and obtains efficiency lift performance curve, the destination file that output suffix is res.
Wherein pump Outer performance calculation expression formula is:
Efficiency Efficiency:
massFlow()@OUTLET*g*Head/(Torque*2950*2*pi/60)
Lift Head:
(massFlowAve(Total Pressure in Stn Frame)@OUTLET-massFlowAve(Total Pressure in Stn Frame)@INLET)/g/997
Torque Torque:
torque_z()@BLADE+torque_z()@SHROUD+torque_z()@HUB
In such scheme, the step 3) modified turbulent model, revised turbulent viscosity is:
Wherein, CμFor model constants;K is tubulence energy;ε is turbulence dissipation rate;F is filter function;ρ is that gas-liquid mixed is close Degree;ρvFor density of gas phase;ρmFor density of liquid phase.
In such scheme, the step 4) set the permanent numerical computations file of cavitation to specifically include following steps:
Step 41) gaseous medium is added into computational fields and medium physical parameter is set according to actual conditions, selection is equal Phase flow model, mass transfer is set to Cavitation, and inputs the saturated vapour pressure of the medium;
Step 42) inlet boundary volume fraction is set, liquid volume fraction is set to 1, gaseous state is set to 0;
Step 43) it regard the res files calculated in the case of non-cavitating as the cavitation under original document one atmospheric pressure of progress Calculate, and back pressure is reduced after the completion of calculating, until occurring cavitation phenomenon, export the permanent result text under required cavitation number Part.
In such scheme, the step 5) cavitation model amendment, empty verbiage model final form is:
pturb=0.39 ρ k
pv=(psat+pturb/2)
Wherein, UFor free stream velocity;L is characterized length;t=L/UIt is characterized time scale;Cdest=9 × 105; Cprod=3 × 104;αvFor gaseous phase volume fraction;PvFor vapor pressure;Psat、PturbFor pressure correcting coefficient.
In such scheme, the step 6) set unsteady cavitation numerical value calculation document to specifically include following steps:
Step 61) by computation model in the permanent calculation document of cavitation be changed to transient state calculate, set transient state calculate total time and Unit interval step-length;
Step 62) coordinate system is revised as instantaneous rotor stator in the surface model of data boundary;
Step 63) the unsteady solver parameter of modification, setting transient state to calculate, data are set and middle spacer step flow field is preserved;
Step 64) set fixed dipole and rotating dipole sound source to extract, pressure fluctuation monitoring point is set, and is applied CEL language writes acoustic power level equation, and exports unsteady computation file;
Step 65) carry out unsteady meter under different cavitation numbers using each res files of the permanent calculating of cavitation as original document Calculate, and it is cgns files to export the suffix comprising pulsation information,
Wherein acoustic power level equation expression formula is:
Mt=sqrt (ke)/c
Wa=0.1*Wden*ed*Mt^5
Wref=6.7*10e-19 [W m^-3]
Lw=10*log (Wa/Wref)
C=1500 [m s^-1].
In formula, Mt is the scale of turbulence;Ke is Turbulent Kinetic;Wa is acoustical power;Wref is reference;Wden is liquid phase Density;C is the velocity of sound in liquid phase;Lw is acoustic power level.
In such scheme, the step 7) it is middle using in the fixed dipole of LMS Virtual.Lab software BEM modules calculating Field noise comprises the following steps:
Step 71) the cgns files comprising fluid grid and Fluid pressure information are imported into acoustic boundary element module calculating pressure Power is pulsed and preserved;
Step 72) LMS Virtual.Lab Software Create boundary element mesh is used, the interior setting analysis model of BEM modules is Direct, interior, element, import boundary element mesh and site grid, and defining trellis-type and carrying out pre-treatment makes net Lattice unit normal direction points to sound field homonymy;
Step 73) fluent material and attribute in vane pump are defined, pressure result file is imported, pulsation is mapped into boundary element Grid;
Step 74) acceleartion boundary condition and sound absorption attribute are defined, start sound field and site distribution is calculated, obtain certain cavitation Several lower internal field site sound pressure levels and radiated curve.
In such scheme, the step 72) the middle specific step using LMS Virtual.Lab Software Create boundary element mesh It is rapid as follows:
Step 721) by file importing LMS Virtual.Lab software of the pump water model suffix for stp, and according to acoustics Size divides surface grids, and acoustic mesh size meets f with calculating frequencymax=c/6L, wherein fmaxFor highest frequency can be calculated, with Calculate time step and set relevant, c is sound spread speed in fluid media (medium), L is grid cell length;
Step 722) enter border element module, and choose the entitled bdf of surface grids output suffix grid file;
Step 723) open bdf grid files, by different characteristic angle will import and export face, import and export pipeline, impeller and Pump chamber is separated for follow-up data mapping.
In such scheme, the step 7) it is middle using in LMS Virtual.Lab software BEM modules calculating rotating dipole Field noise comprises the following steps as follows:
Step 75) the cgns files comprising fluid grid and Fluid pressure information are imported into acoustic boundary element module, by week To defining and calculate concentrfated load on blade;
Step 76) generation boundary element mesh, on border element module by analysis model be set to direct, interior, Element, imports boundary element mesh and site grid, defines trellis-type and carrying out pre-treatment refers to grid cell normal direction To sound field with survey;
Step 77) fluent material and attribute in pump are defined, import blade surface pressure result file;
Step 78) fan sound source and sound absorption attribute are defined, carry out sound-filed simulation and site is calculated, obtain internal field under certain cavitation number Rotating dipole pressure frequency response function.
In such scheme, the step 8) modal calculation, BEM/ are carried out using LMS Virtual.Lab softwares FEM model It is as follows that FEM coupling calculating obtains outfield radiated noise step:
Step 81) file that pump housing suffix is stp is imported, define case material attribute and preserve;
Step 82) divide finite element grid, the grid defined in finite-element module and structure attribute, addition containment surfaces definition Constraint, and carry out finite element structure modal calculation;
Step 83) boundary element mesh and pump body structure mode, addition modal damping and sensing are imported in the element module of border Property site, and define each trellis-type, carry out pre-treatment operation;
Step 84) fluent material attribute and structure attribute are defined, import volute pressure pulsation result and carry out data mapping;
Step 85) acceleartion boundary condition, sound absorption attribute are defined, carry out coupled sound fields distribution calculating and site is calculated, obtain External radiation acoustical power frequency respective function and sound field point to distribution under to certain cavitation coefficient, cavitation factor, Toma coefficient.
Compared with prior art, the beneficial effects of the invention are as follows:
1. the present invention is modified by the inside and outside attribute testing of vane pump to turbulence model and cavitation model, can be accurate Cavitation cycle form is simulated, is that blade pump cavitation Inductive noise numerical forecast provides safeguard.
2. the present invention proposes that a kind of method that can calculate acoustic power flow effect between acoustic fluid domain and thin-wall construction is used for Blade pump cavitation Inductive noise is calculated, the numerical prediction to cavitation Inductive noise is solved the problems, such as.
3. the present invention is according to cavitation form, pressure fluctuation distribution and changing rule, pair of the variation characteristic of radiated sound field Than the accuracy for illustrating this forecasting procedure, design and cavitation diagnosis to altitude performance pump provide foundation.
Brief description of the drawings
Fig. 1 is the flow chart of an embodiment of the present invention;
Fig. 2 is Leaf of embodiment of the present invention pump case and motor schematic three dimensional views;
Fig. 3 is Leaf pump fluid domain schematic three dimensional views of the embodiment of the present invention;
Fig. 4 is experimental simulation comparison diagram under each cavitation number after cavitation model optimization of the embodiment of the present invention;
Fig. 5 is fluid domain boundary element mesh figure in the embodiment of the present invention;
Fig. 6 is internal field noise pattern under different cavitation numbers in the embodiment of the present invention;
Fig. 7 is centre field of embodiment of the present invention radiated sound field distribution map;
Fig. 8 is the outfield directive property distribution map that cavitation number of the embodiment of the present invention is 0.31.
In figure, 1, motor;2nd, axle;3rd, blade pump casing;4th, suction flange;5th, outlet(discharge) flange;6th, outlet extending section water body; 7th, spiral case water body;8th, impeller water body;9th, import extension water body.
Embodiment
The present invention is described in further detail with reference to the accompanying drawings and detailed description, but protection scope of the present invention It is not limited to this.
Fig. 1 show the flow chart of blade pump cavitation Inductive noise Numerical Prediction Method of the present invention, and the present embodiment is adopted Exemplified by a marine centrifugal pump, its main design parameters:Flow Q=20m3/h, lift H=24.5m, rotating speed 2950r/min, Pumped (conveying) medium is seawater, impeller inlet diameter 52mm, impeller outlet diameter 142mm, the number of blade 6, impeller outlet width 8mm, snail Shell base circle diameter (BCD) 152mm, spiral case inlet width 21mm.Enclosed visualized experiment platform is built, external characteristics and sky are obtained by testing Change performance curve, obtain a centrifugation pump cavitation stage internal field vacuole distribution.
The blade pump cavitation Inductive noise Numerical Prediction Method, is mainly used in calculating vane pump in low pressure or big flow The interior outfield radiated noise for occurring cavitation is closed, is comprised the following steps:
Step 1) three-dimension modeling and mesh generation:
Three-dimensional modeling is carried out to pump fluid domain using UG 3D sculpting softwares, and three-dimensional modeling is carried out to pump body structure, and Export the model file that suffix is stp.Be illustrated in figure 2 visible motor 1 in blade pump case and motor schematic three dimensional views, figure, Axle 2, blade pump casing 3, suction flange 4 and outlet(discharge) flange 5;Being illustrated in figure 3 can in vane pump fluid domain schematic three dimensional views, figure See outlet extending section water body 6, spiral case water body 7, impeller water body 8 and import extension water body 9.Mesh generation is carried out to fluid domain, And mesh quality is checked, mesh quality is more than 0.2, and minimum grid angle is more than 18 degree, number of grid 2524115.Export suffix For cfx grid file.
Step 2) non-cavitating numerical computations file is set using CFD software:
Grid is imported in ANSYS CFX softwares and computation model setting is carried out, selection Isothermal heat transfer models and K-Epsilon turbulence models;Fluid media (medium) chooses 25 DEG C of Water at, Media density 998kg/m3;Set and import and export perimeter strip Part, pressure inlets are set to 1atm, and flowexit is set to 20m3/ h, wheel speed is set to 2950rev min^-1;Rotation is calculated Domain wall is set to, without sliding wall, set rotary shaft wall relative velocity consistent with impeller direction of rotation to be set to 0, use Rotating coordinate systems;Static computational fields wall is set to, without sliding wall, be set to global rest frame;Sound computational fields Interface uses MRF Moving references, and transformed coordinate system is set to Frozen Rotor, and Mesh connection mode is set to GGI;It is quiet Quiet computational fields interface link model is set to generic connectivity, and Mesh connection mode is set to GGI;Convective term is set in solver parameter For high resolution mode, turbulent flow numerical term is set to single order windward, solves step number and is set according to actual conditions, convergence criterion residual error is set It is 10^-4 for average value RMS value;It is monitored using CEL language write-in lift efficiency expression formula, output suffix is def's Non-cavitating numerical computations file.
Step 3) amendment turbulence model:
By the inside and outside characterisitic parameter of contrast and experiment value size, turbulence model is determined if appropriate for this status condition, if by mistake Difference is more than 5%, then turbulence model is modified until reaching precision of prediction.
Step 4) the permanent numerical computations file of cavitation is set:
Gaseous medium vapour is added into computational fields, molal weight 1.0kg kmol^-1 are set, homogeneous flow mould is selected Type, mass transfer is set to Cavitation, and inputs the saturated vapour pressure 3574pa of the medium.Inlet boundary volume integral is set Number, is set to 1, gaseous state is set to 0 by liquid volume fraction.The res files calculated in the case of non-cavitating are carried out as original document Cavitation under one atmospheric pressure is calculated, and the permanent result reduced after the completion of calculating under back pressure, cavitation number required for output is literary Part.
Step 5) amendment cavitation model:
By contrasting cavitation performance parameter and experiment value difference, cavitation model is determined if appropriate for operating mode instantly, by repairing Positive gasification coefficient, condensation coefficient and density revision function optimization cavitation model are until reach precision of prediction.
The distribution of cavitation vacuole is simulated after optimization and high-speed photography experiment shooting vacuole distribution is as shown in Figure 4.As can be seen that logical Cavitation distribution and experiment appropriateness substantially that the method is simulated are crossed, vacuole primarily occur ins blade suction surface close to impeller inlet position Put, as cavitation develops, gradually extend to blade outer rim, after cavitation develops to a certain extent, each blade vacuole distribution is connected, Runner is blocked, hydraulic performance is influenceed, steeply rises pump operation vibration noise.
Step 6) unsteady cavitation numerical value calculation document is set:
Computation model in the permanent calculation document of cavitation is changed into transient state to calculate, sets transient state to calculate total time T=0.4s, it is single Position time step is set to Δ T=0.00005787s.Coordinate system is revised as instantaneous rotor stator in the surface model of data boundary. Unsteady solver parameter is changed, sets transient state to calculate spacer step flow field in data and setting and preserves t=0.05s.Set fixed even Extremely son and rotating dipole sound source is extracted, and sets pressure fluctuation monitoring point, and application CEL language write-in acoustic power level equations, and Export unsteady computation file.It is non-fixed under different cavitation numbers that each res files of the permanent calculating of cavitation are carried out as original document Often calculate, and it is cgns files to export the suffix comprising pulsation information.
Step 7) sound field in fixed dipole cavitation Inductive noise is calculated using LMS Virtual.Lab software BEM modules:
Cgns files comprising fluid grid and Fluid pressure information are imported into acoustic boundary element module and calculate pressure fluctuation And preserve.Pump water model suffix is imported into LMS Virtual.Lab softwares for stp file, and divided according to acoustics size Surface grids, it is desirable to which acoustic mesh size is set to 10mm, absolute tolerance is set to 1mm.As shown in Figure 5.Into border element module, and select The entitled bdf of middle surface grids output suffix grid file, opens grid file, by different characteristic angle by import and export face, enter Outlet conduit, pump impeller, turbine and pump chamber are separated for follow-up data mapping.In BEM modules set analysis model be direct, Interior, element, import boundary element mesh, and site grid is arranged at 4 times away from import and export face calibers, define grid class Type and carry out pre-treatment make grid cell normal direction point to sound field homonymy.Using acceleration as boundary condition, border is imported and exported Plate face sound absorption attribute is defined, remaining surface is total reflection wall, and acoustic impedance is 1.5 × 106kg·m-2·s-1, the velocity of sound is Reference sound pressure is 1 × 10 in 1500m/s, water-6Pa, imports pressure result file, pulsation is mapped into boundary element mesh.Beginning sound Field and site distribution are calculated, and obtain site sound pressure level frequency spectrum curve in internal field under several groups of cavitation numbers, as shown in fig. 6, during cavitation inception, Noise sound spectral fluctuations narrow scope, main energetic is concentrated at leaf frequency and its high power association frequency, 0~1000Hz frequency range acoustic pressures Level declines slower.With the development of cavitation, the increase of sound pressure level frequency spectrum fluctuation range, the sound pressure level rising nearby of the frequent section of axle is bright Aobvious, leaf frequency and its frequency multiplication slightly have rising, and overall sound pressure level is in rising trend, more meets with experimental result.
Sound field in rotating dipole cavitation Inductive noise is calculated using LMS Virtual.Lab software BEM modules:
Cgns files comprising fluid grid and Fluid pressure information are imported into acoustic boundary element module, by circumferential definition simultaneously Concentrfated load on blade is calculated, segments is set to 10.Boundary element mesh is generated, analysis model is set to by element module on border Direct, interior, element, import boundary element mesh and site grid, and defining trellis-type and carrying out pre-treatment makes net Lattice unit normal direction points to sound field homonymy.Acoustic impedance is 1.5 × 106kg m-2s-1, and the velocity of sound is that sound is referred in 1500m/s, water Press as 1 × 10-6Pa, imports blade surface pressure result file.The definition fan sound source number of blade is 6, and rotating speed is 2950r/min, Time domain ripple is converted into frequency domain pulsation using quick Fourier transformation during noise calculation, sound-filed simulation and site is carried out Calculate, obtain rotating dipole pressure frequency response function in internal field under certain cavitation number.
Step 8) modal calculation is carried out using LMS Virtual.Lab softwares FEM model, BEM/FEM coupling calculating is obtained Cavitation Inductive noise outfield radiated sound field:
The file that pump housing suffix is stp is imported, definition case material is Iron.Finite element grid is divided, sizing grid is set 2mm is set to, absolute tolerance is set to 1mm, the grid defined in finite-element module and structure attribute, the definition of addition containment surfaces is about Beam, and carry out finite element structure modal calculation.Boundary element mesh and pump body structure mode are imported in the element module of border, mould is added State is damped, and 1% mode is added here, sets 360 circumferential directive property sites, and defines each trellis-type, carries out pre-treatment behaviour Make.Stream fluent material density 998kg/m in definition3, velocity of sound 1500m/s in water, acoustic pressure 1 × 10 in water-6Pa, atmospheric density 1.225kg/m3, THE VELOCITY OF SOUND IN AIR 340m/s, reference sound pressure is 2 × 10 in air-5Pa, imports volute pressure pulsation result and goes forward side by side Row data map.Acceleartion boundary condition, sound absorption attribute are defined, coupled sound fields distribution calculating is carried out and site is calculated, obtain outer Radiated sound field distribution map, acoustic power level spectrogram and directive property distribution map, as shown in Fig. 7~Fig. 8.The display radiation of external radiation sound field Noise main sound source position is located near spiral case cut water position, and it is still hydrodynamic noise to illustrate main noise during cavitation.Noise Frequency is focused primarily upon at leaf frequency and its frequency multiplication, and directive property distribution dipolar nature reduction, maximum sensitivity acoustic pressure appears in 120 °, most Small directive property acoustic pressure appears in 330 °.
It should be understood that, although this specification is described according to each embodiment, but not each embodiment only includes one Individual independent technical scheme, this narrating mode of specification is only that for clarity, those skilled in the art will should say Bright book is as an entirety, and the technical solutions in the various embodiments may also be suitably combined, and forming those skilled in the art can be with The other embodiment of understanding.
The a series of detailed description of those listed above illustrating only for the possible embodiments of the present invention, They simultaneously are not used to limit the scope of the invention, all equivalent embodiments made without departing from skill spirit of the present invention or change It should be included in the scope of the protection.

Claims (9)

1. a kind of blade pump cavitation Inductive noise Numerical Prediction Method, it is characterised in that comprise the following steps:
Step 1) three-dimension modeling and mesh generation:Convection current field domain and pump case structure are modeled, and convection current field domain carries out net Lattice are divided, and check mesh quality;
Step 2) non-cavitating numerical computations file is set using CFD software:By step 1) grid import CFD software, carry out side Boundary's condition setting, turbulence model selection, solver parameter setting calculates non-cavitating numerical computations file;
Step 3) turbulence model is modified until reaching precision of prediction:By poor between contrast test and calculated value external characteristics It is different, turbulence model is determined if appropriate for this status condition, if difference is outside allowed band, return to step 2) reselect turbulent flow Model is modified to turbulence model until precision of prediction is reached, until difference is in allowed band, then into step 4);
Step 4) the permanent numerical computations file of cavitation is set:Gas-liquid two-phase physical parameter is set, initial component boundary condition is set, Set inlet-pressure boundary condition, cavitation model selection;
Step 5) cavitation model is modified until reaching precision of prediction:By contrasting cavitation performance parameter and experiment value difference It is different, cavitation model is determined if appropriate for operating mode instantly, the return to step 4 if not meeting) further repaiied by high speed photography test results Positive gasification coefficient, condensation coefficient and density revision function optimization cavitation model are until reach precision of prediction, until meeting calculating essence Degree requires to enter step 6);
Step 6) unsteady cavitation numerical value calculation document is set:Fixed dipole and rotating dipole sound source, which are extracted, to be set, by sky Change permanent original document and iterate to calculate unsteady numerical value, obtain cavitation form;
Step 7) internal field noise is calculated using LMS Virtual.Lab software BEM modules:Fixed and rotating dipole sound source is defeated Enter, boundary element and site mess generation, carry out fluent material attribute definition, boundary condition is set, sound absorption attribute is set, beginning sound Field and site distribution are calculated;
Step 8) modal calculation is carried out using LMS Virtual.Lab softwares FEM model:BEM/FEM coupling calculating obtains outfield Radiated noise, including FEM meshing, modal calculation, damping addition, coupled vibro-acoustic definition and calculating, obtain certain cavitation system Several lower external radiation acoustical power frequency respective functions and sound field point to distribution.
2. a kind of blade pump cavitation Inductive noise Numerical Prediction Method according to claim 1, it is characterised in that the step It is rapid 1) to specifically include following steps:
Step 11) three-dimensional modeling is carried out to pump fluid domain, and three-dimensional modeling is carried out to the pump housing, output suffix is literary for stp model Part;
Step 12) mesh generation is carried out to fluid domain and housing respectively, and mesh quality is checked, to ensure computational accuracy, adjustment Mesh quality is more than 0.2, and minimum grid angle is more than 18 degree, and exports the grid file that suffix is cfx.
3. a kind of blade pump cavitation Inductive noise Numerical Prediction Method according to claim 1, it is characterised in that the step It is rapid 2) to specifically include following steps:
Step 21) by step 1) grid file import in CFD software and carry out computation model setting, select Isothermal Heat transfer model and k-Epsilon turbulence models;Fluid media (medium) is chosen and sets fluid media (medium) density and viscosity;Turnover is set Mouth boundary condition, sets environmental pressure, sets pump impeller rotating speed, and rotation computational fields wall is set to without sliding wall, setting rotation Axle wall relative velocity consistent with impeller direction of rotation is set to 0, uses rotating coordinate systems;Static computational fields wall is set to Without sliding wall, global rest frame is set to;Sound computational fields interface uses MRF Moving references, transformed coordinate system Frozen Rotor are set to, Mesh connection mode is set to GGI;Computational fields interface link model is set to commonly connect silently Connect, Mesh connection mode is set to GGI;Convective term, turbulent flow numerical term, solution step number, convergence criterion residual error, defeated in solver parameter Go out the non-cavitating numerical computations file that suffix is def;
Step 22) non-cavitating file calculates and obtains efficiency lift performance curve, the destination file that output suffix is res.
4. a kind of blade pump cavitation Inductive noise Numerical Prediction Method according to claim 1, it is characterised in that the step It is rapid 4) to specifically include following steps:
Step 41) gaseous medium is added into computational fields, homogeneous phase model is selected, mass transfer is set to Cavitation, And input the saturated vapour pressure of the medium;
Step 42) inlet boundary volume fraction is set, liquid volume fraction is set to 1, gaseous state is set to 0;
Step 43) calculated the res files calculated in the case of non-cavitating as the cavitation under original document one atmospheric pressure of progress, And back pressure is reduced after the completion of calculating, until occurring cavitation phenomenon.
5. a kind of blade pump cavitation Inductive noise Numerical Prediction Method according to claim 1, it is characterised in that the step It is rapid 6) to specifically include following steps:
Step 61) computation model in the permanent calculation document of cavitation is changed to transient state calculating, set transient state to calculate total time and unit Time step;
Step 62) coordinate system is revised as instantaneous rotor stator in the surface model of data boundary;
Step 63) unsteady solver parameter modification, transient state calculates data setting and middle spacer step flow field preserves and set;
Step 64) dipole and the extraction setting of rotating dipole sound source are fixed, pressure fluctuation monitoring point is set, and exports unsteady meter Calculate file;
Step 65) unsteady computation under different cavitation numbers is carried out using each res files of the permanent calculating of cavitation as original document, and Suffix of the output comprising pulsation information is cgns files.
6. a kind of blade pump cavitation Inductive noise Numerical Prediction Method according to claim 1, it is characterised in that the step Rapid 7) middle use LMS Virtual.Lab software BEM modules calculate fixed dipole internal field noise and comprised the following steps:
Step 71) the cgns files comprising fluid grid and Fluid pressure information are imported into acoustic boundary element module calculating pressure arteries and veins Move and preserve;
Step 72) LMS Virtual.Lab Software Create boundary element mesh is used, the interior setting analysis model of BEM modules is Direct, interior, element, import boundary element mesh and site grid, and defining trellis-type and carrying out pre-treatment makes net Lattice unit normal direction points to sound field homonymy;
Step 73) fluent material and attribute in vane pump are defined, pressure result file is imported, pulsation is mapped into boundary element net Lattice;
Step 74) acceleartion boundary condition and sound absorption attribute are defined, start sound field and site distribution is calculated, obtain under certain cavitation number Internal field site sound pressure level and radiated curve.
7. a kind of blade pump cavitation Inductive noise Numerical Prediction Method according to claim 6, it is characterised in that the step Rapid 72) middle comprising the following steps that using LMS Virtual.Lab Software Create boundary element mesh:
Step 721) by file importing LMS Virtual.Lab software of the pump water model suffix for stp, and according to acoustics size Surface grids are divided, acoustic mesh size meets f with calculating frequencymax=c/6L, wherein fmaxFor highest frequency can be calculated, with calculating Time step sets relevant, and c is sound spread speed in fluid media (medium), and L is grid cell length;
Step 722) enter border element module, and choose the entitled bdf of surface grids output suffix grid file;
Step 723) bdf grid files are opened, face will be imported and exported by different characteristic angle, pipeline, impeller and pump chamber is imported and exported It is separated for follow-up data mapping.
8. a kind of blade pump cavitation Inductive noise Numerical Prediction Method according to claim 1, it is characterised in that the step Rapid 7) middle use LMS Virtual.Lab software BEM modules calculating rotating dipole internal field noise comprises the following steps as follows:
Step 75) the cgns files comprising fluid grid and Fluid pressure information are imported into acoustic boundary element module, by circumferential fixed Justice simultaneously calculates concentrfated load on blade;
Step 76) generation boundary element mesh, on border element module by analysis model be set to direct, interior, Element, imports boundary element mesh and site grid, defines trellis-type and carrying out pre-treatment refers to grid cell normal direction To sound field with survey;
Step 77) fluent material and attribute in pump are defined, import blade surface pressure result file;
Step 78) fan sound source and sound absorption attribute are defined, carry out sound-filed simulation and site is calculated, obtain internal field rotation under certain cavitation number Dipole pressure frequency response function.
9. a kind of blade pump cavitation Inductive noise Numerical Prediction Method according to claim 1, it is characterised in that the step Rapid 8) to carry out modal calculation using LMS Virtual.Lab softwares FEM model, BEM/FEM coupling calculating obtains outfield radiation and made an uproar Sound step is as follows:
Step 81) file that pump housing suffix is stp is imported, define case material attribute and preserve;
Step 82) finite element grid, the grid defined in finite-element module and structure attribute are divided, the definition of addition containment surfaces is about Beam, and carry out finite element structure modal calculation;
Step 83) boundary element mesh and pump body structure mode, addition modal damping and directive property are imported in the element module of border Point, and each trellis-type is defined, carry out pre-treatment operation;
Step 84) fluent material attribute and structure attribute are defined, import volute pressure pulsation result and carry out data mapping;
Step 85) acceleartion boundary condition, sound absorption attribute are defined, carry out coupled sound fields distribution calculating and site is calculated, obtain certain External radiation acoustical power frequency respective function and sound field point to distribution under cavitation coefficient, cavitation factor, Toma coefficient.
CN201710325748.2A 2017-05-10 2017-05-10 A kind of blade pump cavitation Inductive noise Numerical Prediction Method Pending CN107273570A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710325748.2A CN107273570A (en) 2017-05-10 2017-05-10 A kind of blade pump cavitation Inductive noise Numerical Prediction Method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710325748.2A CN107273570A (en) 2017-05-10 2017-05-10 A kind of blade pump cavitation Inductive noise Numerical Prediction Method

Publications (1)

Publication Number Publication Date
CN107273570A true CN107273570A (en) 2017-10-20

Family

ID=60074074

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710325748.2A Pending CN107273570A (en) 2017-05-10 2017-05-10 A kind of blade pump cavitation Inductive noise Numerical Prediction Method

Country Status (1)

Country Link
CN (1) CN107273570A (en)

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108036923A (en) * 2018-02-05 2018-05-15 中国石油大学(华东) A kind of convergence-divergence form laminar condition cavitation test device and design method
CN108412805A (en) * 2018-03-06 2018-08-17 清华大学 A kind of cylinder inducer design method, cylindrical inducer and its flow field
CN108629098A (en) * 2018-04-25 2018-10-09 上海法雷奥汽车电器系统有限公司 Method for numerical simulation, device, computing device and the storage medium that generator wind is made an uproar
CN109190166A (en) * 2018-07-31 2019-01-11 江苏大学 A kind of blade pump cavitation determines and state evaluating method and its system
CN109829243A (en) * 2019-02-25 2019-05-31 大连理工大学 A kind of combined type aero-engine high-pressure gear pump modeling and simulating and fault filling method
CN110096812A (en) * 2019-05-05 2019-08-06 湖南凯利特泵业有限公司 A kind of cavitation performance of centrifugal pump automatic simulation method based on CFD platform
CN110206743A (en) * 2019-05-28 2019-09-06 浙江大学 A kind of axial-flow pump cavitation characterization extracting method compareed based on Noise texture and bubble shape
CN110929457A (en) * 2019-11-15 2020-03-27 东南大学 Numerical prediction method for dipole source noise of travelling crane with water surface
CN111104733A (en) * 2019-12-05 2020-05-05 江苏大学 System and method for simulating operating environment of marine pump
CN111400941A (en) * 2019-01-03 2020-07-10 江苏大学 Numerical prediction method for internal reflux and reflux vortex cavitation of vane pump
CN113065292A (en) * 2021-04-25 2021-07-02 中国航发湖南动力机械研究所 Boundary model method for integrated turbine transition section unsteady flow calculation
CN113392601A (en) * 2021-07-12 2021-09-14 东南大学 Geometric feature analysis method of hydrodynamic noise line spectrum of cylinder-like shell
CN114526900A (en) * 2022-02-21 2022-05-24 重庆川仪调节阀有限公司 Cavitation identification method, experiment method and experiment device for flow characteristic experiment of regulating valve
CN114611176A (en) * 2022-03-23 2022-06-10 天津大学 Vinasse-wastewater two-phase flow characteristic simulation method and system in hydraulic pulper
CN117349992A (en) * 2023-12-06 2024-01-05 陕西空天信息技术有限公司 Inducer parameterization design method and device and electronic equipment
CN111199118B (en) * 2019-12-13 2024-05-03 同济大学 Non-flush dual-culvert spout jet noise calculation method, calculation device and storage medium
CN118013894A (en) * 2024-04-09 2024-05-10 蓝深集团股份有限公司 Hydraulic excitation optimization control method based on radial force distribution of pump body
CN111104733B (en) * 2019-12-05 2024-06-07 江苏大学 Marine pump operation environment simulation system and method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103631994A (en) * 2013-11-18 2014-03-12 代元军 Method for predicting acoustic radiation rule value of noise of wind turbine
CN104462652A (en) * 2014-11-07 2015-03-25 中国人民解放军海军工程大学 Design method of hydraulic model of shaftless drive type integrated motor pump jet
CN106050730A (en) * 2016-07-14 2016-10-26 西华大学 Vane pump and impeller blade for vane pump
CN106503323A (en) * 2016-10-17 2017-03-15 江苏大学 A kind of centrifugal multistage pump multiple centrifugal pump flow-induction structural radiation noise numerical prediction method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103631994A (en) * 2013-11-18 2014-03-12 代元军 Method for predicting acoustic radiation rule value of noise of wind turbine
CN104462652A (en) * 2014-11-07 2015-03-25 中国人民解放军海军工程大学 Design method of hydraulic model of shaftless drive type integrated motor pump jet
CN106050730A (en) * 2016-07-14 2016-10-26 西华大学 Vane pump and impeller blade for vane pump
CN106503323A (en) * 2016-10-17 2017-03-15 江苏大学 A kind of centrifugal multistage pump multiple centrifugal pump flow-induction structural radiation noise numerical prediction method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
司乔瑞等: "基于CFD/CA 的离心泵流动诱导噪声数值预测", 《机械工程学报》 *

Cited By (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108036923B (en) * 2018-02-05 2024-02-02 中国石油大学(华东) Design method of convergence-divergence type laminar flow state cavitation test device
CN108036923A (en) * 2018-02-05 2018-05-15 中国石油大学(华东) A kind of convergence-divergence form laminar condition cavitation test device and design method
CN108412805A (en) * 2018-03-06 2018-08-17 清华大学 A kind of cylinder inducer design method, cylindrical inducer and its flow field
CN108412805B (en) * 2018-03-06 2020-01-07 清华大学 Cylindrical inducer design method, cylindrical inducer and flow field of cylindrical inducer
CN108629098A (en) * 2018-04-25 2018-10-09 上海法雷奥汽车电器系统有限公司 Method for numerical simulation, device, computing device and the storage medium that generator wind is made an uproar
CN109190166A (en) * 2018-07-31 2019-01-11 江苏大学 A kind of blade pump cavitation determines and state evaluating method and its system
CN109190166B (en) * 2018-07-31 2023-03-24 江苏大学 Cavitation judgment and state evaluation method and system for vane pump
CN111400941A (en) * 2019-01-03 2020-07-10 江苏大学 Numerical prediction method for internal reflux and reflux vortex cavitation of vane pump
CN109829243B (en) * 2019-02-25 2022-11-18 大连理工大学 Modeling simulation and fault injection method for high-pressure gear pump of combined aero-engine
CN109829243A (en) * 2019-02-25 2019-05-31 大连理工大学 A kind of combined type aero-engine high-pressure gear pump modeling and simulating and fault filling method
CN110096812A (en) * 2019-05-05 2019-08-06 湖南凯利特泵业有限公司 A kind of cavitation performance of centrifugal pump automatic simulation method based on CFD platform
CN110206743A (en) * 2019-05-28 2019-09-06 浙江大学 A kind of axial-flow pump cavitation characterization extracting method compareed based on Noise texture and bubble shape
CN110929457B (en) * 2019-11-15 2023-04-21 东南大学 Numerical prediction method for dipole source noise of traveling crane with waterway surface
CN110929457A (en) * 2019-11-15 2020-03-27 东南大学 Numerical prediction method for dipole source noise of travelling crane with water surface
CN111104733A (en) * 2019-12-05 2020-05-05 江苏大学 System and method for simulating operating environment of marine pump
CN111104733B (en) * 2019-12-05 2024-06-07 江苏大学 Marine pump operation environment simulation system and method
CN111199118B (en) * 2019-12-13 2024-05-03 同济大学 Non-flush dual-culvert spout jet noise calculation method, calculation device and storage medium
CN113065292B (en) * 2021-04-25 2022-07-22 中国航发湖南动力机械研究所 Boundary model method for integrated turbine transition section unsteady flow calculation
CN113065292A (en) * 2021-04-25 2021-07-02 中国航发湖南动力机械研究所 Boundary model method for integrated turbine transition section unsteady flow calculation
CN113392601A (en) * 2021-07-12 2021-09-14 东南大学 Geometric feature analysis method of hydrodynamic noise line spectrum of cylinder-like shell
CN114526900A (en) * 2022-02-21 2022-05-24 重庆川仪调节阀有限公司 Cavitation identification method, experiment method and experiment device for flow characteristic experiment of regulating valve
CN114611176A (en) * 2022-03-23 2022-06-10 天津大学 Vinasse-wastewater two-phase flow characteristic simulation method and system in hydraulic pulper
CN117349992A (en) * 2023-12-06 2024-01-05 陕西空天信息技术有限公司 Inducer parameterization design method and device and electronic equipment
CN117349992B (en) * 2023-12-06 2024-03-08 陕西空天信息技术有限公司 Inducer parameterization design method and device and electronic equipment
CN118013894A (en) * 2024-04-09 2024-05-10 蓝深集团股份有限公司 Hydraulic excitation optimization control method based on radial force distribution of pump body

Similar Documents

Publication Publication Date Title
CN107273570A (en) A kind of blade pump cavitation Inductive noise Numerical Prediction Method
Majidi Numerical study of unsteady flow in a centrifugal pump
Gonza´ lez et al. Numerical simulation of the dynamic effects due to impeller-volute interaction in a centrifugal pump
Luo et al. Tip leakage flow and aeroacoustics analysis of a low-speed axial fan
Zhang et al. Unsteady hydrodynamic forces due to rotor-stator interaction on a diffuser pump with identical number of vanes on the impeller and diffuser
CN103077292A (en) Method for predicting hydraulic noise of centrifugal pump
Torabi et al. The effect of viscosity on performance of a low specific speed centrifugal pump
Rahmati et al. Nonlinear time and frequency domain methods for multirow aeromechanical analysis
Chirkov et al. CFD simulation of pressure and discharge surge in Francis turbine at off-design conditions
Liu et al. Assessment of a turbulence model for numerical predictions of sheet-cavitating flows in centrifugal pumps?
Liu et al. Numerical and experimental studies of hydraulic noise induced by surface dipole sources in a centrifugal pump
Khelladi et al. Flow study in the impeller–diffuser interface of a vaned centrifugal fan
Dai et al. Pressure fluctuation and its influencing factors in circulating water pump
Shojaeefard et al. A parametric study for improving the centrifugal pump impeller for use in viscous fluid pumping
Ji et al. Application of Wray–Agarwal turbulence model in flow simulation of a centrifugal pump with semispiral suction chamber
Si et al. Assessment of cavitation noise in a centrifugal pump using acoustic finite element method and spherical cavity radiation theory
Mao et al. A brief review of numerical solving methods for internal fluid of pumped storage unit
Li et al. Effects of guide vane thickness on pressure pulsation of mixed-flow pump in pumped-storage power station
Pütz Prediction of Rotating Instabilities in Low Pressure Steam Turbines Operating at Low Load
Dai et al. Study on noise characteristics of marine centrifugal pump under different cavitation stages
Lu et al. Study on pressure pulsation induced by cavitation at the tongue of the volute in a centrifugal pump
Ivanov et al. Technique for axial pump characteristics predicting in CFD package Openfoam
Jawad et al. Numerical investigation on the effect of impeller trimming on the performance of a modified compressor
Hamidreza et al. An investigation on turbocharger turbine performance parameters under inlet pulsating flow
Vahdati et al. Computational study of intake duct effects on fan flutter stability

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20171020