URANS COMPUTATION OF CAVITATING FLOWS AROUND SKEWED PROPELLERS

Cavitating flows around skewed propellers are investigated numerically by means of the unsteady Reynolds Averaged Navier-Stokes (RANS) Equation method. The standard k-ε turbulence and the modified Z-G-B cavitation models are employed. A measured nominal wake is used for the inlet velocity boundary condition. Predicted cavitating evolution processes and tip cavity patterns are compared with experimental observations. In addition, the influence of the skew angles on the cavitation and unsteadiness performances of propellers operating in a non-uniform wake is also studied. Results show that the modified Z-G-B cavitation model performs better to simulate the cavitating flow cases studied in this paper. Comparisons demonstrate that the skewed propeller with a skew angle of 20 o is the best choice for a given stern wake with a assigned thrust and the minimum force fluctuations.     

UNSTEADY FLOW ANALYSIS AND EXPERIMENTAL INVESTIGATION OF AXIAL-FLOW PUMP

The three-dimensional unsteady turbulent flow in axial-flow pumps was simulated based on Navier-Stoke solver embedded with k-ε RNG turbulence model and SIMPLEC algorithm.Numerical results show that the unsteady prediction results are more accurate than the steady results,and the maximal error of unsteady prediction is only 4.54%.The time-domain spectrums show that the static pressure fluctuation curves at the inlet and outlet of the rotor and the outlet of the stator are periodic,and all have four peaks and four valleys.The pressure fluctuation amplitude increases from the hub to the tip at the inlet and outlet of the rotor,but decreases at the outlet of the stator.The pressure fluctuation amplitude is the greatest at the inlet of the rotor,and the average amplitude decreases sharply from the inlet to the outlet.The frequency spectrums obtained by Fast Fourier Transform(FFT) show that the dominant frequency is approximately equal to the blade passing frequency.The static pressure on the pressure side of hydrofoil on different stream surfaces remains almost consistent,and increases gradually from the blade inlet to the exit on the suction side at different time steps.The axial velocity distribution is periodic and is affected by the stator blade number at the rotor exit.The experimental results show that the flow is almost axial and the pre-rotation is very small at the rotor inlet under the conditions of 0.8 QN -1.2 QN.Due to the clearance leakage,the pressure,circulation and meridional velocity at the rotor outlet all decrease near the hub leakage and tip clearance regions.     

NUMERICAL SIMULATIONS OF CAVITATING FLOWS

BRIEFINTRODUCTIONOFTHEPAPER :　Thisthesismainlyresearchesonthehydrodynamiccharacteristicsandmechanismfor 3Dcavitatingflowsaroundaxisymmetricbodies ,aswellasthe2Dcavitatingflowsaroundhydrofoils .Anewcavi tatingflowmodel,whichinvolvesviscousandmul ti phaseeffects,isestablishedintwo phaseflowcategory .Accordingtothelocalizedvariationofdensitywithinpredominantlyincompressiblewatermediumandthecharacteristicsofsoundspeedinwater vapormixture ,arelationbetweendensityandpressureisassumed ,in…     

Shedding frequency of sheet cavitation around axisymmetric body at small angles of attack

Cavity shedding of cavitating flows around an axisymmetric body belongs to the unsteady cavitating flows in the condition of steady incoming current.The periodic characteristics of unsteady cavitating flows around an axisymmetric body at small angles of attack are investigated experimentally and numerically.The evolution and shedding process of the three-dimensional sheet cavitation are computed numerically by the Reynolds averaged Navier-Stokes equations and the RNG k-?model.The modification approach for eddy viscosity coefficient in the transition area of the two-phase flow is adopted to reproduce the shedding process of cavitating flows.The computed frequency of the cavity shedding coincides with the experimental data for the cases of unsteady cavitating flows around axisymmetric bodies with four headforms.Given the cavitation number,the shedding process of the cavitating flow depends heavily on the headform of the axisymmetric body.If the angle of attack of the axisymmetric body is greater than a critical value,the violent shedding of the sheet cavitation seems to be depressed.     

NUMERICAL SIMULATIONS OF 2D PERIODIC UNSTEADY CAVITATING FLOWS

1. INTRODUCTION Cavitation is a natural phenomenon especially existing in liquids. A cavitating flow generally involves a large number of vapor structures such as bubbles or vortices which are convecting downstream. When they reach high pressure zones, th…     

低比转速离心泵空化性能的数值模拟

为了对低比转速离心泵空化发生时的内部气液两相流动进行分析,应用计算流体力学(CFD)软件对比转速为66的离心泵空化性能进行气液两相流场的数值研究。对不同有效空化余量时叶轮内部气泡分布的研究表明,泵在进口压力较高时就已经在叶片的进口背面产生空化初生,临界空化余量和许用空化余量时,气泡在叶片表面和流道内部均有分布,占据了部分叶轮流道,影响叶轮内部能量交换,并可能对泵造成空蚀。对该实型泵进行试验研究,验证了数值计算的精确度。     

NUMERICAL PREDICTION OF BLADE FREQUENCY NOISE OF CAVITA- TING PROPELLER

The blade frequency noise of a cavitating propeller in a uniform flow is analyzed in the time domain. The unsteady loading (of a dipole source) and the sheet cavity volume (of a monopole source) on the propeller surface are calculated by a potential-based surface panel method. Then the time-dependent pressure and the cavity volume data are used as the input for the Fowcs Williams-Hawkings formulation to predict the acoustics pressure. The integration of the noise source is performed over the true blade surface rather than the ideal blade surface without thickness. The noise characteristics of the cavitating propeller are discussed. With the sheet cavitation, the thickness (cavitation) noise is larger than the loading noise and is the dominant noise source.The noise directivity is not as clear as that of the noise under a non-cavitation condition. The cavitation noise is attenuated mores lowly than the non-cavitation noise.     

THE WALL EFFECT ON VENTILATED CAVITATING FLOWS IN CLOSED CAVITATION TUNNELS

For ventilated cavitating flows in a closed water tunnel, the wall effect may exert an important influence on cavity shape and hydrodynamics, An isotropic mixture multiphase model was established to study the wall effect based on the RANS equations, coupled with a natural cavitation model and the RNG k-ε turbulent model. The governing equations were discretized using the finite volume method and solved by the Gauss-Seidel linear equation solver on the basis of a segregation algorithm. The algebraic multigrid approach was carried through to accelerate the convergence of solution. The steady ventilated cavitating flows in water tunnels of different diameter were simulated for a conceptual underwater vehicle model which had a disk cavitator. It is found that the choked cavitation number derived is close to the approximate solution of natural cavitating flow for a 3-D disk. The critical ventilation rate falls with decreasing diameter of the water tunnel. However, the cavity size and drag coeflicient are rising with the decrease in tunnel diameter for the same ventilation rate, and the cavity size will be much different in water tunnels of different diameter even for the same ventilated cavitation number.     

On the numerical simulations of vortical cavitating flows around various hydrofoils

This paper reviews the numerical models of various cavitating flows around hydrofoils. Numerical models relating to cavitation flows, including mass transfer models and turbulence models, are summarized at first. Then numerical results and analysis of flow characteristics for the cavitating flows around twisted hydrofoils, truncated hydrofoil and tip leakage are discussed respectively. For mean flow fields, Reynolds averaged Navier-Stokes(RANS) simulation associated with a kind of nonlinear turbulence model is found to be an economic and robust numerical approach for different kinds of cavitating flows including cloud cavitation, tip cavitation and tip leakage cavitation. To predict the fluctuations of pressure and velocity, large eddy simulation(LES) and detached eddy simulation(DES) are two effective approaches. Finally, a few open questions are proposed for future research.     

An integral calculation approach for numerical simulation of cavitating flow around a marine propeller behind the ship hull

In this paper, the unsteady cavitating turbulent flow around a marine propeller is simulated based on the unsteady Reynolds averaged Navier-Stokes(URANS) with emphasis on the hull-propeller interaction by an integral calculation approach, which means the propeller and hull are treated as a whole when the cavitating flow is calculated. The whole calculational domain is split to an inner rotating domain containing a propeller and an outer domain containing a hull. And the two split sections are connected together in ANSYS CFX by using the GGI interfaces and the transient rotor stator frame change/mixing model. The alternate rotation model is employed for the advection term in the momentum equations in order to reduce the numerical error. Comparison of predictions with measurements shows that the propeller thrust coefficient can be predicted satisfactorily. The unsteady cavitating flow around the propeller behind the ship hull wake shows quasi-periodic features including cavity inception, growth and shrinking. These features are effectively reproduced in the simulations which compare well to available experimental data. In addition, significant pressure fluctuations on the ship hull surface induced by the unsteady propeller cavitation are compared with experimental data at monitoring points on the hull surface. The predicted amplitudes of the first components corresponding to the first blade passing frequencies match well with the experimental data. The maximum error between the predictions and the experimental data for the pressure pulsations is around 8%, which is acceptable in most engineering applications.     

ASSESSMENT OF A CENTRAL DIFFERENCE FINITE VOLUME SCHEME FOR MODELING OF CAVITATING FLOWS USING PRECONDITIONED MULTIPHASE EULER EQUATIONS

A numerical treatment for the prediction of cavitating flows is presented and assessed.The algorithm uses the preconditioned multiphase Euler equations with appropriate mass transfer terms.A central difference finite volume scheme with suitable dissipation terms to account for density jumps across the cavity interface is shown to yield an effective method for solving the multiphase Euler equations.The Euler equations are utilized herein for the cavitation modeling,because some certain characteristics of cav...     

绕扭曲翼型三维非定常空泡脱落结构的数值分析

该文以一种三维扭曲翼型为研究对象,采用均衡流空化模型和基于滤波器的湍流修改模式模拟了绕翼型非定常空化流动。计算结果显示,在一定条件下绕翼型的空泡流动存在两种脱落结构,即主体脱落和二次脱落。进一步的流场分析表明回射流和侧面射流是这两种空泡脱落结构产生的主要原因。计算结果与试验的比较说明,该文采用的数值模拟方法可成功模拟实验观测的空泡脱落现象。     

AN APPROACH IN MODELING TWO-DIMENSIONAL PARTIALLY CAVITATING FLOW

1. INTRODUCTIONCavitation is a phenomenon that widely occurs in liquids. When the flow leads to the local pressure below the saturation vapor pressure of water (the fluid medium, in our research), usually incipient cavitation firstly occurs on the solid b…     

Large eddy simulation of turbulent attached cavitating flow with special emphasis on large scale structures of the hydrofoil wake and turbulence-cavitation interactions

In this paper, the turbulent attached cavitating flow around a Clark-Y hydrofoil is investigated by the large eddy simulation(LES) method coupled with a homogeneous cavitation model. The predicted lift coefficient and the cavity volume show a distinctly quasi-periodic process with cavitation shedding and the results agree fairly well with the available experimental data. The present simulation accurately captures the main features of the unsteady cavitation transient behavior including the attached cavity growth, the sheet/cloud cavitation transition and the cloud cavitation collapse. The vortex shedding structure from a hydrofoil cavitating wake is identified by the Q- criterion, which implies that the large scale structures might slide and roll down along the suction side of the hydrofoil while being further developed at the downstream. Further analysis demonstrates that the turbulence level of the flow is clearly related to the cavitation and the turbulence velocity fluctuation is much influenced by the cavity shedding.     

Determination of frequencies of oscillations of cloud cavitation on a 2-D hydrofoil from high-speed camera observations

A method is presented to determine significant frequencies of oscillations of cavitation structures from high-speed camera recordings of a flow around a 2-D hydrofoil. The top view of the suction side of an NACA 2412 hydrofoil is studied in a transparent test section of a cavitation tunnel for selected cloud cavitation regimes with strong oscillations induced by the leading-edge cavity shedding. The ability of the method to accurately determine the dominant oscillation frequencies is confirmed by pressure measurements. The method can resolve subtle flow characteristics that are not visible to the naked eye. The method can be used for noninvasive experimental studies of oscillations in cavitating flows with adequate visual access when pressure measurements are not available or when such measurements would disturb the flow.     

Application of modified κ-ω model to predicting cavitating flow in centrifugal pump

Considering the compressibility of the cavity in the cavitating flow, this paper presents a modified κ-ω model for predicting the cavitating flow in a centrifugal pump, in which the modified κ-ω model and Schnerr-Sauer cavitation model were combined with ANSYS CFX. To evaluate the modified and standard κ-ω models, numerical simulations were performed with these two models, respectively, and the calculation results were compared with the experimental data. Numerical simulations were executed with three different values of the flow coefficient, and the simulation results of the modified κ-ω model showed agreement with most of the experimental data. The cavitating flow in the centrifugal pump obtained by the modified κ-ω model at the design flow coefficient of 0.102, was analyzed. When the cavitation number decreases, the cavity initially generates on the suction side of the blade near the leading edge and then expands to the outlet of the impeller, and the decrease of the total pressure coefficient mainly occurs upstream of the impeller passage, while the downstream remains almost unaffected by the development of cavitation.     

INTERFERENCE OF SIDE STRUT WITH THE NATURAL CAVITATING FLOWS AROUND A SUBMERGED VEHICLE IN WATER TUNNEL EXPERIMENTS

To apply the measurements of model experiment in water tunnel to the actual sailing condition,it is necessary to know accurately the strut effect and its rule.In the present work,the corresponding interferences of one-side strut and two-side strut on the natural cavitating flows around a submerged vehicle in water tunnel were investigated numerically,using the homogeneous equilibrium two-phase model coupled with a natural cavitation model.The numerical simulation results show that the strut types have distinct effects on the hydrodynamic properties.For the same given upstream velocity and downstream pressure,the existence of the strut leads to an increment of natural cavitation number,reduces the low-pressure region and depresses the pressure on the vehicle surface near the sides of strut.In the case of given cavitaiton number,the influences of the two-side strut on the drag and lift coefficients are both enhanced along with the increment of attack angle,however the influence of the one-side strut gradually gets stronger on the drag coefficient but weaker on the lift coefficient contrarily.In addition,based on the present numerical results,a correction method by introducing the sigmoidal logistic function is proposed to eliminate the interference from the foil-shaped strut.     

NUMERICAL SIMULATION OF VENTILATED CAVITATING FLOW AROUND A 2D FOIL

By using a pressure-based method and the finite volume method in the framework of the time dependent Reynolds-averaged Navier-Stokes equations, the authors studied the unsteady process of ventilated cavities generated forcing air through an orifice in a 2D hydrofoil without natural cavitation physically. The computation was carried out with the Volume Of Fluid （VOF） technique to track the gas-liquid two-phase interface. The results of simulation indicate that the ventilation rate is an important parameter in determining the morphology of cavity. There exists a critical value to convert sheet cavity into supereavity. A high ventilation rate can induce a two phase interface fluctuation and enable the ventilated eavitating flow to present a characteristic of periodicity.     

双吸泵小流量工况下的空化特性研究

为研究双吸泵在小流量工况下叶轮内部空化特性,同时进一步说明小流量工况相比于设计工况时的空化特性差异,结合均质两相流模型和SST k-ω湍流模型,对双吸泵小流量工况和设计工况下的全流道空化流场进行数值模拟,以分析不同流量工况下空化分布与发展情况,以及空化对各叶片载荷造成的影响。研究结果表明:适当减小双吸泵进口流量,有助于改善双吸泵的空化性能;在小流量工况下空化首先发生于叶片吸力面头部靠后盖板附近,而且此处的空泡体积分数最大,这一空化特征同设计工况有所差异;随着NPSH的降低,叶轮内空化不断加强,但是小流量工况下的空化强度始终不及设计工况;不同空化状态会导致叶片吸力面压力的变化,从而表现为叶片表面载荷分布的变化。     

Numerical analysis of cavitation within slanted axial-flow pump

In this paper, the cavitating flow within a slanted axial-flow pump is numerically researched. The hydraulic and cavitation performance of the slanted axial-flow pump under different operation conditio...