Расчетно-экспериментальное исследование кавитационного обтекания гидрокрыла NACA0015

Abstract

Одним из основных возобновляемых источников энергии являются гидроэнергетические ресурсы, используемые для выработки электрической энергии на гидроэлектростанциях. Существенную научную и техническую проблему для проектирования оборудования гидроэлектростанций представляет собой кавитация, обусловленная большими скоростями течения воды вблизи лопастей рабочего колеса. Появление кавитации характерно и для других турбомашин, насосов, гидравлического оборудования, движителей судов и т. д. Кавитация приводит к снижению энергетической эффективности и повышенному износу оборудования. Потребности моделирования турбулентного обтекания лопастных систем, работа которых сопровождается кавитацией, требуют развития современных численных методов, способных с приемлемой точностью предсказывать возникновение кавитации и описывать динамику течения. Цель работы: исследование кавитационного течения в окрестности модельного гидропрофиля NACA0015 с помощью экспериментальных и численных методов, сопоставление результатов моделирования с данными измерений, анализ влияния модели турбулентности на результаты расчета течения в пограничном слое со стороны разрежения гидрокрыла. Методы исследования. Пространственная структура и динамика кавитационных каверн изучалась с помощью высокоскоростной визуализации, по данным которой также производилась оценка интегральных параметров каверн. Двумерные распределения средней скорости и турбулентных характеристик вблизи гидрокрыла были получены методом анемометрии по изображениям частиц (PIV) как в одно- (некавитирующих), так и в двухфазных (кавитирующих) течениях. При численном моделировании течений использовались методы вычислительной гидродинамики, основанные на решении уравнений Рейнольдса для турбулентного течения с помощью метода контрольного объема на трехмерной сетке из гексаэдральных ячеек. Учет дисперсной фазы (кавитационных пузырей и их скоплений) производился путем решения уравнения переноса доли пара. Турбулентность описывалась с помощью двухпараметрической модели k-[омега] SST, дифференциальной модели переноса рейнольдсовых напряжений (RSM) и метода моделирования отсоединенных вихрей (DES). Результаты. В экспериментах был получен полный набор данных, позволяющих проводить прямой анализ результатов измерений и численного моделирования, а также сравнение расчетных моделей между собой. Выполненные расчеты показали, что даже в случае небольшой паровой каверны за передней кромкой гидропрофиля для корректного описания присоединения потока к поверхности ниже по течению требуется привлечение метода замыкания второго порядка - дифференциальной модели переноса рейнольдсовых напряжений. Когда размер каверны становится большим, нестационарность потока является определяющим фактором турбулентного переноса. Моделирование с использованием вихреразрешающих методов позволило выявить периодическую динамику паровой каверны на низких частотах. Присоединение потока за каверной в нестационарном случае также лучше предсказывается моделью переноса рейнольдсовых напряжений.

Hydropower resources, used for power generation in hidraulic power plants, are one of the major renewable energy sources. A significant scientific and technical challenge for designing the equipment of hydroelectric power plants is a cavitation caused by high velocity of water flow near the runner blades. Cavitation is also typical for other turbomachinery, pumps, hydraulic equipment, ship propulsors, etc. Cavitation leads to decrease in energy efficiency and increased wear of equipment. The needs in modeling turbulent flow around blade assemblies, the operation of which is accompanied by cavitation, require the development of modern numerical approaches capable of forecasting the cavitation occurrence and describing flow dynamics with satisfactory accuracy. The main aim of the research is to investigate cavitation flow in a vicinity of NACA0015 hydrofoil by means of experimental and numerical methods, to compare the simulation and measurement results, to analyze the effect of turbulence model on calculation of the boundary layer flow over the suction side of a hydrofoil. The methods used in the study. The spatial structure and dynamics of partial gas-vapor cavities were studied by high-speed imaging which was also used to estimate integral parameters of the cavities. Two-dimensional distributions of mean velocity and turbulent characteristics around the hydrofoil were measured by Particle Image Velocimetry (PIV) in both single- (non-cavitating) and two-phase (cavitating) flows. In numerical simulation of flows the authors have applied the methods of computational fluid dynamics based on solving the Reynolds equations for turbulent flow, using the control volume approach on a three-dimensional mesh consisting of hexahedral cells. The dispersed phase (cavitation bubbles and their clouds) was calculated by solving the equation of vapor fraction transfer. Turbulence was described by means of k-[omega] SST two-equation model, differential model of the Reynolds stress transfer (RSM) and the detached eddy simulation method (DES). The results. The authors obtained full set of data in the experiments. The data allow direct analysis of the results of numerical modeling and measurements and comparing the computational models. The carried out calculations showed that, even in the case of a small vapor cavity behind the leading edge of the hydrofoil, the use of the second-order closure method - differential model of the Reynolds stress transfer - was required to describe correctly the flow reattachment to the wall downstream. When the cavity became long, the flow unsteadiness was a determinant factor of turbulent transfer. Modeling with the eddy-resolving approaches allowed identifying the periodic dynamics of vapor cavities at low frequencies. In an unsteady case, the flow reattachment downstream of the cavity was better predicted by the model of the Reynolds stress transfer.