Journal of Siberian Federal University. Engineering & Technologies 5 (2011 4) 489-503
УДК 532.5.032
Численное и экспериментальное моделирование течения в отсасывающей трубе гидротурбины
С.В. Алексеенкоаб, А.А. Дектеревав, И. В. Литвинов5, А.В. Минакова,в, И.М. Пылевг, А.И. Шандроа, С.И. Шторкаб
а Институт теплофизики им. С. С. Кутателадзе СО РАН Россия 630090, Новосибирск, ак. Лаврентьева, 1 б Новосибирский государственный университет, Россия 630090, Новосибирск-90, ул. Пирогова, 2 в Сибирский федеральный университет, Россия 660041, Красноярск, пр. Свободный, 79 г ОАО «Силовые машины» - ЛМЗ, Россия 195009, Санкт-Петербург, Свердловская наб., 18 1
Received 4.10.2011, received in revised form 11.10.2011, accepted 18.10.2011
В работе проведено экспериментальное и численное изучение структуры закрученного турбулентного потока в модельной отсасывающей трубе гидротурбины. Получено хорошее согласование расчетных и экспериментальных данных. Обнаружена сильная неравномерность потока за поворотным коленом, существенно снижающая пропускную способность отсасывающей трубы.
Ключевые слова: гидротурбина, отсасывающая труба, прецессия вихря, DES-метод.
Введение
В настоящее время на гидроэнергетику приходится 87 % энергии, производимой из возобновляемых источников, на основе которых, в свою очередь, производят одну пятую всей энергии, потребляемой в мире [1]. Повышение устойчивости работы, увеличение КПД гидроэлектростанций невозможны без изучения физических механизмов гидродинамических процессов, имеющих место в проточных частях гидротурбин [2]. Существенную роль среди таких процессов играют нестационарные явления, связанные с образованием крупномасштабных вихревых структур. Одним из механизмов генерации пульсаций потока служит прецессия вихревого жгута [3], образующегося за рабочим колесом гидротурбины в режимах недогрузки или перегрузки работы гидротурбины, когда поток после прохождения через гидротурбину имеет достаточно большую остаточную закрутку. Для прогнозирования резонансных явлений, поис-
* Corresponding author E-mail address: [email protected]
1 © Siberian Federal University. All rights reserved
ка методов подавления неустойчивости необходима детальная информация о характеристиках и структуре потока.
Проведение детальных исследований на натурных ГЭС представляет большую сложность. Выходом может быть применение методов математического (численного) и экспериментального моделирования с использованием лабораторных установок. В контексте сказанного выше целью данной работы стало численное и экспериментальное исследование структуры течения в модельной отсасывающей трубе гидротурбины. Основной упор сделан на выявление особенностей течения, прежде всего, наличия интенсивных пульсаций и неравномерности потока. Экспериментальные исследования проводили с использованием замкнутого гидродинамического контура с максимальным расходом до 100-150 м3/ч. На данном этапе исследования вихревого течения, возникающего в отсасывающей трубе гидротурбины, проводили с использованием рабочего участка, моделирующего геометрию "Ти^ВШЕ-99". Эта геометрия, которую применяют для тестирования численных моделей течения в отсасывающих трубах гидротурбин [4], также хорошо адаптирована к оптическим методам исследования потока. В ходе опытов с использованием цифровой камеры и ЛДА были получены данные по визуализации потока и полям скоростей в различных сечениях отсасывающей трубы. Экспериментальные данные использованы для тестирования результатов численного моделирования. Было продемонстрировано хорошее соответствие измерений и расчетов.
Экспериментальная техника
Стенд представляет собой замкнутый гидравлический контур (рис. 1), состоящий из центробежного химического насоса Х100-80-160Т153 номинальной производительностью 100 м3/ч и напором 3,2 атм. с двигателем 22 кВт, баком емкостью 500 л и собственно рабочего участка (модели отсасывающей трубы). Для снижения потерь давления контур выполнен так, что в нем установлено всего два поворотных колена. Соединения частей контура (трубы, отводы, повороты, краны) выполнены с применением стандартной трубопроводной арматуры из полиэтиле-
Рис.1. Схема экспериментального участка
на (ПЭ80, ПЭ100). Регулирование расхода осуществляется регулированием частоты вращения вала с помощью частотного преобразователя (ЧП) ВЕСПЕР Е2-8300-030Н 22 кВт 380В. Для измерения расхода воды щпЖкШЗт^Нпяпшш с электронным блоком ДО -
800 - 11. Паспортная погрешность измерения |р асхода в рабочем диапазоне для данного стенда состшШЗЯгсЯшщЯЯПЕИ
В качестве модельной геометрии отсасывающей трубы выбрана конфигурация, использованная в проекте TURBINE-99 [4]. Модельный рабочий участок был выполнен целиком из органического стекла, что дает возможность полного оптического доступа ко всему полю течения во всех частях канала и, соответственно, широкого применения оптических бесконтакт-нпщкцЯИЯаЯГДшЯй оящего
иЖцйщщЯЯЯЯщЖГТц На
вВйднЕИшыЗшИшЕиЕашЕдоЕЕВЭ^
пшЯщцйщшЗЗЯщ завих-
рЭёдЕшз!
5=2 3
1 -((1,1 а2 )3 ' (1 -к/ ( 2 )2).
1ап(^), где dl и d2 - меньший и больший диаметры кольцевого канала
завихрителя (рис. 2), ф - угол наклона лопаток. Значение параметра крутки в нашей модели, равное, как и в прототипе TURBINE-99, 0,4, соответствует условиям работы натурных турбин. Приведенные в данной статье результаты были получены при достаточно большом числе Рей-нольдса Re=1,3•Ю5.
Визуализация потока в модели отсасывающей трубы осуществляется с помощью подсветки мелких пузырьков воздуха, подаваемого на вход центробежного насоса. Размеры и концентрация воздушных пузырьков могут регулироваться изменением подачи воздуха, а также варьированием оборотов двигателя насоса и расходом жидкости через байпасную линию. За счет границы межфазного раздела пузыри воздуха хорошо отражают свет и поэтому дают яркий, хорошо видимый трек. В качестве источника света используется 3Вт непрерывный твердо-
Рис. 2. ГеШЯтЩНЯиЩяяштшшЯ Ш1КЩ ЙИУд^ТТЯмЯЖГщШН одится по
часовой стрелке, если смотреть сверху
тельный лазер с диодной накачкой KLM-532/h-3000. Луч лазера преобразовывается с помощью цилиндрической линзы в тонкий световой нож для визуализации течения в выбранном сечении. Регистрация картины течения производилась с помощью цифровой камеры CANON EOS 20D. Для съемки нестационарных процессов использовали скоростные камеры Видео Спринт и PCO 1200hs.
Измеряли локальные скорости с помощью автоматизированной системы ЛДА ЛАД-05 производства ИТ СО РАН. Оптическая система установлена на программно-управляемое 3-координатное перемещающее устройство, позволяющее позиционировать измерительный объем ЛДА по трем координатам с точностью до 12 мкм. Погрешность измерения средней скорости составляет не более 0,5 %.
Расчетная методика
Для сопоставления экспериментальных данных, полученных Hai стенде ИТ СО РАН, про-веденасерия расчетов с целью тестирования и адаптаци и расчетноймодели.
Согласно современным представлениямтрехмерные нестационарные течения вязкой жидкости описываются системойуравнений Навье-Стокса, состоящей из уравнения сохранения массы илиуравнения неразрывное ти:
f + V|» = 0, (1)
и уравнений д вижен ия и ли з акона сохраненияимпульса:
-^ + V(/7V-v) = -VJp + V(t) + F, (2)
где v - вектор скорости жидкости, т - тензор вязких напряжений, F - вектор объемных сил, р - статическоедавление,р - плотность.
Составляющие тензора вязких напряжений^ определяютсякак
(dut du2 du, Л
— +—----—L
дх ( дх, 3 ' dxt
\ j 1 к
гдед- динамическая(молекулярная) вязкость, и, - компоненты вектора скорости, 5У - символ Кронекера.
Считается,чтопри помощи уравнения Навье-Стокса можно описать как ламинарное, так и турбулентное течение, однако насегодня возможности существующей вычислительной техники таковы, что практически все реальные достижения в расчетах турбулентных течений связаны с использованием полуэмпирических моделей турбулентности, использующих подход Рейнольдса. Суть этого подхода заключается в решении осредненных уравнений Навье-Стокса [6,7]:
др\
dt
- + V(pv • v) = -Vp + V(t - pv \') + К,
где V - усредненное по времени скоростей, -р\ • у - тензор рейнольдсовыхнапряжений.
При построении двухпараметрических моделей турбулентности для определения компонент тензора рейнольдсовых напряжений используют гипотезу Буссинеска о изотропной тур-булентнойвязкости:
-р V -V = nt
''du, 9и ^ ^
-+ —L
ЗХ: дХ:
\ 1 '
2
1 дщ
дХгУ
гдед, - турбулентнаявязкостьД - кинетическаяэнергия турбулентных пульсаций.
Использование концепции Буссинеска позволяетпостроить целый рццдвухпараметриче-ских полуэмпирических моделей. В данной работедля замыкания системы RANSypaBHemifi двухзонная дву хпараметрическая модель Ментера SS Т.
Модель Ме нтера записывае тся путем суперпозициимоделейк-еик-ю,основаннойна то м, что моделитипак-е х^чше описываютсвойствасвободных сдвиговых течений,амодели к-ю имеютпреимущество при моделировании пристеночных течений. Плавный переход от1 к-ю модели в пристеночной области к k-емоделивдали от твердыхстенок обеспечивается введе-ниемвесовой эмпирической функции F1. 1Вторая важная деталь модели состоитв изменении стандартной связи между к,аи турбулентнойвязкостью. Модификация состоит
вовведении перехода ic формуле Бредшоу в пристеночной области. Согласно предположению Бредшоунапряжениесдвига впограничном cj^c^e пропорционально энергиитурбулентныс пульсаций.
Ниже пр иведе ныосновныеуравнения M-SST моеел и[7]: ^ + V(pv-k) = V((¡и + ak/л,) -Vk) + P -ß'pmk ,
^ + ) = V((M + aaMl) -V®) + yp—-ßp <o 2 + (1-^;)-(2 pam2-V* -V®)'
ot pt m
В членегенерациитурбулентнойэнергиивведен ограничитель:
P/U ~
P = tti] —- P = min(P, 20 • ß pmk) .
dxs
Весовая фу нкция и ее аргумент:
т = tanh(arg14)
. k , А:1'2 500 w 4ра,Л' ar& = m(n(max(—*-,--); ^ "%)
где положительная часть перекрестных диффузионных членов в уравнении переноса ю:
CDhm 2 тах^/лг^—V2.V®> 10^). т
Выр;шения длявихревой вязкости с у четом гипотезы Бредшоу: ц, = - ^а2с
шах(а1й?; Р20.)
гдевеличина завихренности - □ =20^.
Функцияпереключения ¥2 определяетсяподобно ¥1 :
Р2 = tanh(arg22)
.. к112 500« '
arg2 = таХ(2—-, --)
р ту ра>у
Константы в уравнениях переноса записывают путем суперпозиции констант для модели к-ю модели ^Псох)и константстандартной к-е модели:
^ = ф = {стк,у, ß\.
Набор константдля пристеночного слоя SST модели: ик1 = 0.85 <гт1 = 0.5 Д=0.075
ß' 4ß*
Набор констант для свободных сдвиговых слоев: ак2 = 1.0 ат1 =0.856 Д, = 0.0828
ß\
ß' 4F'
Другие константы, использующиеся в модели:/? == 0.0Î? ат = 0.41 ^=0.31.
Примоделировании течений в гидротурбинах интерес представляетповедение нетолыко осредненных характеристик течения, но и пульсационных (интенсивность пульсаций и частота). Получившее широкоераспр остранение вин женерных расчетах URANS модели турбулентности занижают турбулентные пульсации в закрученном потоке. Чтобы улучшить моделирование нестационарных турбулентных закрученных течений, используют вихреразрешающие методы, как, например, метод моделирования крупных вихрей (Large Eddy Simulation - LES). Однакоего применение требует очень детальнтй сетки, особенно вблизи стенок. В то же ¡время RANS модели достаточно экономичны и хорошоописывают пограничные слои.Чтобы соединить достоинства этихподходов вработе [8], был предложен методмоделирования от-соединенныхвихреfi(DES).ITeрвая Bepc^DES базироваласьнамодели Спяларта-Алмароса. В дальнейшем метод DES стал использоватьсяс другимимоделямитуреулентности,появились различныеего ц:ни.
Метод DES сочетает подходы RANS и LES. В областях потока, где размер вычислительной сетки Д достаточен для разрешения энергонесущих вихрей, т.е. при Д < L, применяется метод LES, а в остальной области - RANS. Перваявеесия DES базироваласьна модели Спаларта-Алмаросa SA, в которойв качестве латактерного линейного масштаба турбулентноетииспонь-зуется ¡расстояние до стенке d. При моделировании отсоединенных вихрейэтот масштаб за-меняетсяфункцией:
d =min {d,CDESis),
где CDES - эмпирическая константа, а Д определяется какмаксимальныйиз трех размеров контрольного объема Дх, Д,, Д2. Таким образом, в пограничных слоях работают RANS модели, а в остальной области - LES.
Такимобразом, для моделирования нестационарного турбулентного течения в данной работе использовала метод DES, основанный на на Zt-coSST модели Meнаера и ограничении турбулентной вязкости: ßt=M,/Fmsа Fdes = шах(4/(CBJlsA),l), Lt = Нъ11\е, CDES = 0,61.
Метод DES позволяет рассмотреть крупномасштабные турбулентные пульсации во всем проточном тракте гидротурбины и оценить их амплитудно-частотные характеристики.
Для аппроксимации конвективных членов уравнений гидродинамики взята схема QUICK -второго порядка точности. Аппроксимация нестационарного члена при URANS и DES моделирования проводилась неявной схемой второго порядка. Величина временного шага при этом удовлетворялаусловиюCFL<2.
Рис. 3. Фрагмент расчетной сетки
ГздшЕтщуЛпммцЗиЛшяш на входе
аций на
входе равнялся 1 %, характерный масштаб турбулентности 0,1 м. На выходе задавались условия фиксированного давления. Закрутка потока в эксперименте обеспечивалась завихрителем, состоящим из 12 лопастей, наклоненных под углом 30 градусов к оси обтекателя.
Для расчета использовалась неструктурированная сетка, состоящая из 420000 полигран-ньлЯТлаЗздЛНЯнЯПШННпК
с т_ЗЯВЯИЕиЕЕ52пЯ5ЭЕ2еВ5щ о описы-
ва ЯЯшЯшкЯВдТИЯЯВД
Щ&щпВГ^ЖщЯуТщщщщ
щВЯщфЯГЯщЯИшЯщ ВиКщСШщшщщш ованы с
пошшГцЯТ£ТуЩЛгащшпшшЩ. ®кшйщ5ЯщВ51у щЖЕЯРЯГЯГЯекЗщЯЯшВЯ скостей
виВЯЯРВк^п^ТЯщЯРшВващВаИ сече-
ницИщщЗщщЗЯЯЯЗЕЕиа^ форми-
рование в поворотной части канала крупномасштабного вихря с горизонтальной осью. Можно увидеть также наличие в конической части выраженной приосевой и пристенных струй, направленных вниз. Особенностью течения выступает также неравномерность течения в диффузоре отсасывающей трубы (за поворотной частью), как видно из визуализации потока в продольных сечениях (рис. 5). Вдоль левой стенки наблюдается выраженный интенсивный проток жидкости в сторону выхода (рис. 5б), в то время как вдоль правой стенки видно хаотичное турбулизированное движение со слабым средним потоком в направлении выхода (рис. 5а). На основе анализа данных скоростной видеосъемки в поперечных сечениях можно также гово-
- 495 -
б)
Рис. 4. Визуализация течения« продольномментральноЛгечении трубы:!а) экспериментальное; б) расчетное поле скорости
рить об отсутствии в конической части отсасывающей трубы выраженного прецессирующего вихревого ядра. Отсутствие выраженной несимметрии мгновенного поля течения в поперечном сечении, характерной для закрученных потоков с прецессирующим вихревым ядром [6], можно видеть на рис. 6.
Количественные данные о поле течения были получены с помощью ЛДА системы. Радиальные профили двух компонент скоростей, измеренные и рассчитанные в поперечном сечении за обтекателем завихрителя, приведены на рис. 7-8. Профиль осевой скорости (рис. 8),
- 496 -
б)
Рис. 5. (а) - визуализация течения в продольном сечении вдоль правой стенки (вид со стороны правой стенки); (б) - визуализация течения в продольном сечении вдоль левой 1стенки (вид ро стороны левой стенки).
направленной вниз, демонстрирует наличие центральной и пристенной кольцевых струй на входе. Разделение потока на две кольцевые струи отражается также в форме профиля окружной компоненты скорости (рис. 8). Профиль радиальной скорости показывает, что центральная кольцевая струя за телом обтекания имеет тенденцию сужения. Внешняя кольцевая струя при этом расширяется, следуя вдоль стенки расширяющегося конуса. Спектральный анализ сигнала ЛДА при измерениях в области конической части канала не выявил определенной
б)
Рис. 6. Визуализация течения в поперечном сечении под обтекателем: а) экспериментальное; б) расчетное поле скорости
дискретной частоты, которая обычно выделяется в потоке с прецессирующим вихрем [9]. Этот факт подтверждает результат, полученный на основе визуальных наблюдений, об отсутствии прецессии вихревого жгута в нашей установке.
Как видно из графиков на рис.7-8, расчетная и экспериментальная формы обоих профилей скорости в данном сечении качественно и количественно хорошо согласуются. Вместе с тем есть и некоторые расхождения, вызванные сложностью экспериментальных измерений в диф-
- 498 -
-1.5 -1-1-1-
-0.1 "0.05 0 0.05 0.1
X
в-в Эксперимент.
Рис. 7. Профиль тангенциальнойскоростиподтелом обтекания
X
в-в Эксперимент,
Рис. 8. Профиль осевой скоростиподтелом обтекания
дДмЯЯ НВВВЯЯтГиИщВшшДшмш льных
профилей скорости. Тем не менее, для такой сложной задачи степень совпадения расчета и эксперимента можно считать достаточно хорошей.
Положение измерительного сечения для поля течения в диффузорной части показаны на рцддИяНВЯВшЯ ЯГлщЖЯйцщвЯжщщЯГО
сйЯИяГЯюЯтЕШИш асчет.
-ВЦ-
1,9 1,6 1.4 1,2 0,91 0,67 0,42 0,17 -0,070
Рис. 9. ИЙвЯцТуйтщЯЯйпвт^цЯЯТЯдГмпга» РЩ! ание
произведено через среднерасходную скорость для сечения Х=430 мм 0.4 м/с: а) экспериментальные изощщт! яя Тр^ЧзхЯЗЗщгам^РщрпВ^РбТЯж^Т
Измерения производили с шагом 10 мм при расходе 13 л/с, соответствующем числу Рейнольдса порядка 105. Как показывает анализ результатов, особенностью течения является выраженная неравномерность потока в горловине выхода из поворотного колена. Имеется интенсивный проток в левом нижнем углу (если смотреть в сторону выхода), а большая часть сечения занята медленным потоком. Положение области минимума примерно совпадает с центром крупномасштабного вихря, видного на картине течения, полученной на основе визуализации. Далее к выходу неравномерность сглаживается, но, тем не менее, остается заметной. Выявленная силь-
- 500 -
ная неравномерность течения в диффузоре, очевидно, приводит к дополнительным потерям энергии, снижающим эффективность работы отсасывающей трубы.
В расчете наблюдается аналогичная структура течения. В этом можно убедиться из сопоставления экспериментальной и расчетной картин течения в различных сечениях, приведенных на рис. 4, 6, 9. Как видно, структура течения в расчете качественно хорошо согласуется с экспериментом. Также наблюдается ярко выраженная неравномерность потока в поворотном колене. Имеется интенсивный проток в левом нижнем углу трубы. Примерно совпадают положения центров крупномасштабного вихря.
Заключение
Таким образом, в данной работе было проведено экспериментальное и численное исследование структуры течения в модельной геометрии отсасывающей трубы гидротурбины (TURBINE-99). Полученные результаты позволяют сделать ряд выводов об особенностях течения в модели данной конфигурации.
Наиболее интересным результатом является то, что в наших исследованиях при данных входных параметрах с завихрителем S=0.4 и при данном расходном режиме, соответствующем турбулентному режиму течения, не образуется выраженное прецессирующее вихревое ядро и не наблюдаются связанные с ним сильные пульсации потока. Этот факт был подтвержден на основе визуальных исследований с использованием скоростной камеры, а также на основе спектрального анализа ЛДА сигнала. Нужно отметить, что в предыдущих исследованиях для аналогичной конфигурации TURBINE-99 и той же степени закрутки потока прецессия вихревого жгута была зафиксирована в расчетах и эксперименте [4, 6]. Отсутствие прецессии вихревого ядра в нашей работе можно объяснить конструкцией использованного завихрите-ля, формирующего специфичный профиль осевой скорости на входе в конус отсасывающей трубы. Прецессия вихря всегда сопровождается зоной возвратного осевого течения под телом обтекания, однако наши входные условия формируют центральную приосевую струю, препятствующую образованию центральной рециркуляционной зоны. В работе [10] показано, что струя воды, испускаемая из отверстия на центральном теле обтекания, может полностью подавить прецессию вихревого ядра. Особенностью нашей работы является то, что устойчивый режим течения без образования прецессии вихря в отсасывающей трубе был получен без использования отдельного канала для подачи воды. Предполагается, что этот результат даст возможность для разработки рекомендаций по конструкции гидротурбины, не генерирующей мощные пульсации потока в отсасывающей трубе.
Еще одной особенностью течения было наличие крупномасштабных вихрей в области выхода из конуса и на выходе из поворотного колена, что может способствовать дополнительной диссипации энергии.
Другой важный результат - обнаруженная сильная неравномерность потока за поворотным коленом. В результате того, что закрутка потока осуществляется по часовой стрелке, а у канала есть выделенное направление, интенсивный поток прижимается к левой стенке. У противоположной стенки, где направление закрутки противоположно направлению течения, поток тормозится, образуя обширную застойную область. Таким образом, проходное сечение в диффузоре отсасывающей трубы используется лишь менее чем на 1/3, что сильно снижает ее
пропускную способность. Этот эффект, очевидно, также будет приводить к дополнительным потерям энергии и требует специального внимания при разработке конструкций практических гидротурбин.
Полученные результаты совместного экспериментального и численного моделирования будут использоваться для оптимизации геометрии отсасывающей трубы для снижения потерь энергии потока и расширения диапазона устойчивой работы гидротурбин.
Работа выполнена при поддержке Министерства образования и науки Российской Федерации в рамках реализации Федеральной целевой программы «Научные и научно-педагогические кадры инновационной России» на 2009 - 2013 годы.
Список литературы
1. World Energy Council, "Survey of Energy Resources", Technical Report, (URL: http://www. worldenergy.org/wec-geis/default.asp), Cited Sep. 2006.
2. Губин М. Ф. Отсасывающие трубы гидроэлектростанций. М.: Энергия, 1970. 270 с.
3. Cassidy J.J., Falvey H.T. Observation of unsteady flow arising after vortex breakdown // Journal of Fluid Mechanics. 41 - 1970. P. 727-736.
4. Cervantes M. J., Engstrom T. F., Gustavsson L. H., Proceedings of the third IAHR/ERCOFTAC Workshop on draft tube flows. Turbine-99 III // Lulea University of Technology, 2005.
5. Гупта А., Лили Д., Сайред Н. Закрученные потоки. М.: Мир, 1987. 590 с.
6. Dekterev A.A., Gavrilov A.A., Minakov A.V., Finnikov K.A. Numerical simulation of unsteady cavitating turbulent flow in water turbine // Proc. of 6th International Symposium On Turbulence, Heat and Mass Transfer, Rome, Italy, 14-18 September. 2009. P. 835-838.
7. Menter F.R. Zonal two equation k-- turbulence models for aerodynamic flows // AIAA Paper. 1993. № 93-2906. 21 p.
8. Spalart P.R., Jou W.-H., Strelets M., Allmaras S.R. Comments on the feasibility of LES for wings and on a hybrid, RANS/LES approach // In Lue, C. and Lue, Z., (eds) Advances in DNS/LES, Proceedings of 1st AFOSR International Conference on DNS/LES, Ruston, LA, August, 4-8, Greyden Press, Columbus, OH. 1997. Pp. 137-147.
9. Shtork S.I., Vieira N.F., Fernandes E.C. On the identification of helical instabilities in a reacting swirling flow // Fuel. 2008. Vol. 87, No. 10-11. P. 2314-2321.
10. Muntean S. , Susan-Resiga R. , Bosioc A. Pressure measurements in a conical diffuser with swirling flow and axial jet control. Center for Advanced Research in Engineering Sciences, National Center for Engineering of Systems with Complex Fluids, "Politehnica" University of Timioara, Romnia. URL:http://www.acadtim.tm.edu.ro/iSMARTflow/pdf/Muntean_4GRoWTH.pdf), 2008.
Numerical and Experimental Simulation fo the Flow in the Water Turbine Draft Tube
Sergey V. Alekseenkoa,b, Aleksander A. Dektereva,c, Ivan V. Litvinovc, Andrey V. Minakovac, Igor M. Pylevd, Andrey I. Shandroa and Sergey I. Shtorka,b
a Institute of Thermophysics SB RAS 1 ak. Lavrentiev. Novosibirsk, 630090 Russia b Novosibirsk State University Pirogov st., Novosibirsk-90, 630090 Russia c Siberian Federal University 79 Svobodny, Krasnoyarsk, 660041 Russia d «Power machines» LMZ 18 Sverdlovskay, Sant- Peterburg, 195009 Russia
We performed experimental and numerical study of the structure of swirling turbulent flow in the model draft tube turbine. Obtained a good agreement between calculated and experimental data. A strong non-uniformity of flow at a turning knee, severely reduces the capacity of the draft tube.
Keywords: water turbine, draft tube, vortex rope, DES method.