Olivier Coutier-Delgosha

Evaluation of the Turbulence Model Influence on the Numerical Simulations of Unsteady Cavitation

Unsteady cavitation in a Venturi-type section was simulated by two-dimensional computations of viscous, compressible, and turbulent cavitating flows. The numerical model used an implicit finite volume scheme (based on the SIMPLE algorithm) to solve Reynolds-averaged Navier-Stokes equations, associated with a barotropic vapor/liquid state law that strongly links the density variations to the pressure evolution. To simulate turbulence effects on cavitating flows, four different models were implemented (standard k−e RNG; modified k−e RNG; k−ω with and without compressibility effects), and numerical results obtained were compared to experimental ones. The standard models k−e RNG and k−ω without compressibility effects lead to a poor description of the self-oscillation behavior of the cavitating flow. To improve numerical simulations by taking into account the influence of the compressibility of the two-phase medium on turbulence, two other models were implemented in the numerical code: a modified k−e model and the k−ω model including compressibility effects. Results obtained concerning void ratio, velocity fields, and cavitation unsteady behavior were found in good agreement with experimental ones. The role of the compressibility effects on turbulent two-phase flow modeling was analyzed, and it seemed to be of primary importance in numerical simulations.

A joint experimental and numerical study of mechanisms associated to instability of partial cavitation on two-dimensional hydrofoil

The present work was carried out in the scope of a numerical-experimental collaborative research program, whose main objective is to understand the mechanisms of instabilities in partial cavitating flow. Experiments were conducted in the configuration of a rectangular foil located in a cavitation tunnel. Partial cavitation was investigated by multipoint wall-pressure measurements together with lift and drag measurements and numerical videos. The computations were conducted on two-dimensional hydrofoil section and are based on a single fluid model of cavitation: the liquid/vapor mixture is considered as a homogeneous fluid whose composition is regulated by a barotropic state law. The algorithm of resolution is derived from the SIMPLE approach, modified to take into account the high compressibility of the medium. Several physical features were pointed out by this joint approach. Particularly two distinct cavity self-oscillation dynamics characterized by two different frequencies (dynamics 1 and dynamics 2) ...

Experimental and Numerical Studies in a Centrifugal Pump With Two-Dimensional Curved Blades in Cavitating Condition

In the presented study a special test pump with two-dimensional curvature blade geometry was investigated in cavitating and noncavitating conditions using different experimental techniques and a three-dimensional numerical model implemented to study cavitating flows. Experimental and numerical results concerning pump characteristics and performance breakdown were compared at different flow conditions. Appearing types of cavitation and the spatial distribution of vapor structures within the impeller were also analyzed. These results show the ability of the model to simulate the complex three-dimensional development of cavitation in a rotating machinery, and the associated effects on the performance.

Influence of the Cavitation Model on the Simulation of Cloud Cavitation on 2D Foil Section

For numerical simulations of cavitating flows, many physical models are currently used. One approach is the void fraction transport equation-based model including source terms for vaporization and condensation processes. Various source terms have been proposed by different researchers. However, they have been tested only in different flow configurations, which make direct comparisons between the results difficult. A comparative study, based on the expression of the source terms as a function of the pressure, is presented in the present paper. This analytical approach demonstrates a large resemblance between the models, and it also clarifies the influence of the model parameters on the vaporization and condensation terms and, therefore, on the cavity shape and behavior. Some of the models were also tested using a 2D CFD code in configurations of cavitation on two-dimensional foil sections. Void fraction distributions and frequency of the cavity oscillations were compared to existing experimental measurements. These numerical results confirm the analytical study.

Simulation of unsteady cavitation with a two-equation turbulence model including compressibility effects

Unsteady effects associated with cavitation were investigated by numerical simulations in three configurations. The simplest one was a Venturi-type section in which the cavitation sheet oscillates periodically with vapour cloud shedding. The second one was a hydrofoil whose unsteady cavitating behaviour depends on the angle of attack, and the most complex one was a cascade of three hydrofoils. In this last configuration, in addition to the unsteadiness associated with each cavity, a coupling between the three channels was also observed. These cavitating flows were simulated by 2D computations. Resolution of Reynolds-averaged Navier-Stokes equations was based on a finite-volume discretization associated with a pressure correction algorithm. Cavitation was simulated by using a barotropic vapour/liquid state law that links the fluid density evolution to the pressure variations. As standard k–ϵ RNG or k–ω turbulence models were found to be weakly efficient to simulate unsteady cavitation, influence of the com...

Internal structure and dynamics of sheet cavitation

The internal structure and the dynamics of two-dimensional (2D) sheet cavitation on the suction side of a 2D foil section were investigated experimentally. Experiments were conducted in a cavitation tunnel and situations ranging from steady sheet cavitation to unsteady cloud cavitation were obtained by varying the foil incidence and the cavitation number. Using a novel endoscopic technique, coupled with x-ray attenuation measurements, the two-phase morphology and the void fraction within the sheet cavitation were investigated. Supplemental information on the instantaneous shape of the sheet cavity and its instability frequency were also obtained by visualization and pressure measurements, respectively. The investigations focused on (a) the void fraction distribution and (b) the frequency of the cavity oscillations. It was found that the void ratio can reach as much as 50% depending on the conditions of operation, and the Strouhal numbers are around 0.25 in the case of partial cavity instability and 0.12 i...

Numerical Prediction of Cavitating Flow on a Two-Dimensional Symmetrical Hydrofoil and Comparison to Experiments

This paper presents comparisons between two-dimensional (2D) CFD simulations and experimental investigations of the cavitating flow around a symmetrical 2D hydrofoil. This configuration was proposed as a test case in the “Workshop on physical models and CFD tools for computation of cavitating flows” at the 5th International Symposium on cavitation, which was held in Osaka in November 2003. The calculations were carried out in the ENSTA laboratory (Palaiseau, France), and the experimental visualizations and measurements were performed in the IRENav cavitation tunnel (Brest, France). The calculations are based on a single-fluid approach of the cavitating flow: the liquid/vapor mixture is treated as a homogeneous fluid whose density is controlled by a barotropic state law. Results presented in the paper focus on cavitation inception, the shape and the general behavior of the sheet cavity, lift and drag forces without and with cavitation, wall pressure signals around the foil, and the frequency of the oscillations in the case of unsteady sheet cavitation. The ability of the numerical model to predict successively the noncavitating flow field, nearly steady sheet cavitation, unsteady cloud cavitation, and finally nearly supercavitating flow is discussed. It is shown that the unsteady features of the flow are correctly predicted by the model, while some subtle arrangements of the two-phase flow during the condensation process are not reproduced. A comparison between the peer numerical results obtained by several authors in the same flow configuration is also performed. Not only the cavitation model and the turbulence model, but also the numerical treatment of the equations, are found to have a strong influence on the results.

Numerical Model to Predict Unsteady Cavitating Flow Behavior in Inducer Blade Cascades

The cavitation behavior of a four-blade rocket engine turbopump inducer is simulated. A two-dimensional numerical model of unsteady cavitation was applied to a blade cascade drawn from an inducer geometry. The physical model is based on a homogeneous approach of cavitation, coupled with a barotropic state law for the liquid/vapor mixture. The numerical resolution uses a pressure-correction method derived from the SIMPLE algorithm and a finite volume discretization. Unsteady behavior of sheet cavities attached to the blade suction side depends on the flow rate and cavitation number. Two different unstable configurations of cavitation are identified. The mechanisms that are responsible for these unstable behaviors are discussed, and the stress fluctuations induced on the blade by cavitation instabilities are estimated.

Numerical simulation of the unsteady cavitation behavior of an inducer blade cascade

One source of unsteadiness in turbopump inducers consists in a rotating cavitation behavior, characterized by different cavity shapes on the different blades, which leads to super- or subsynchronous disturbances. This phenomenon is simulated for the case of a simple two-dimensional blade cascade corresponding to a typical four-blade inducer. A numerical model of unsteady cavitating flows was adapted to take into account nonmatching connections and periodicity conditions. Single-channel and four-channel computations were performed, and in the latter case, nonsymetrical unstable flow patterns were obtained. Limits of stability according to the mass flow rate and the cavitation number are presented. Qualitative comparisons with experiments, instability criterion, and the mechanisms of instabilities are also investigated

Experimental Study of a Cavitating Centrifugal Pump During Fast Startups

The start-up of rocket engine turbopumps is generally performed only in a few seconds. It implies that these pumps reach their nominal operating conditions after only a few rotations. During these first rotations of the blades, the flow evolution in the pump is governed by transient phenomena, based mainly on the flow rate and rotation speed evolution. These phenomena progressively become negligible when the steady behavior is reached. The pump transient behavior induces significant pressure fluctuations which may result in partial flow vaporization, i.e. cavitation. An existing experimental test rig has been updated in the LML laboratory (Lille, France) for the start-ups of a centrifugal pump. The study focuses on cavitation induced during the pump start-up. Instantaneous measurement of torque, flow rate, inlet and outlet unsteady pressures, and pump rotation velocity enable to characterize the pump behavior during rapid starting periods. Three different types of fast start-up behaviors have been identified. According to the final operating point, the start-up is characterized either by a single drop of the delivery static pressure, by several low-frequency drops, or by a water hammer phenomenon that can be observed both a the inlet and outlet of the pump. A physical analysis is proposed to explain these three different types of transient flow behavior.