CC BY
32
УДК 681.121.842+532.5
И. А. Яценко, А. Н. Сабирзянов, В. А. Фафурин, В. Б. Явкин
КОЭФФИЦИЕНТ ИСТЕЧЕНИЯ ТОРЦЕВОГО СОПЛА БЕЗ ВХОДНОГО УЧАСТКА
Ключевые слова: расход газа, торцевое сопло, эталон массового расхода газа, коэффициент истечения, неопределенность измерений, число Рейнольдса, относительный диаметр, математическое моделирование, модели турбулентности,
стандартное сопло.
Проведены экспериментальные и численные исследования коэффициента истечения торцевых сопел без входного участка. Эксперименты проводились на государственном эталоне массового расхода газа ФГУП ВНИИР. Получены экспериментальные зависимости коэффициента истечения торцевого сопла от числа Рейнольдса и относительного диаметра. Численные исследования проведены с использованием современных средств вычислительной гидродинамики. Представлены результаты моделирования коэффициента истечения сопла ИСА 1932 в условиях стандартного монтажа и показана возможность рассчитывать его в пределах неопределенности, установленной стандартом. Проведено сопоставление результатов численного моделирования коэффициента истечения торцевых сопел с полученными экспериментальными данными. Расхождение экспериментальных и расчетных значений коэффициента истечения торцевых сопел, полученных по моделям турбулентности со стандартным набором модельных констант, не превышает 1,5 % в диапазоне чисел Рейнольдса выше 105. Более удовлетворительное совпадение экспериментальных и расчетных данных получено при модификации констант в моделях турбулентности.
Keywords: gas flowrate, end-holed nozzle, mass flowrate standard, discharge coefficient, measurement uncertainty, Reynolds number, relative diameter, mathematic modeling, turbulence model, standard nozzle.
The article describes experimental and numerical investigation of discharge coefficient of end-holed nozzles. Experiments have been held on a national standard of mass flow rate at Federal State Unitary Enterprise "All-Russian Research Institute ofFlow Metering". As a result a dependence of discharge coefficient of end-holed nozzle and orifice on Reynolds number and relative diameter has been acquired. Numerical investigations have been held using contemporary CFD methods. Article contains results of mathematical modeling of discharge coefficient of ISA 1932 nozzle for a standard mounting thus proving the ability of its calculation within uncertainty of an appropriate standard. Also article gives the comparison of the result of numerical modeling and experimental investigation. Deviation of the calculated vales on and experimental values for a turbulence models with a standard set of constants for a Reynolds number more than 105 does not exceed 1,5 %. Better A satisfactory agreement between the experimental and calculated data has been obtained by modifying constants in turbulence models.
Введение
Для измерения расходов жидкостей и газов при заборе из большого объема или при сбросе в большой объем могут использоваться торцевые сужающие устройства (СУ), которые монтируются, соответственно, на входе или на выходе измерительного трубопровода (ИТ). Первый случай измерений удобен при испытаниях различных компрессоров, горелочных устройств и т.п. Второй - при решении экологических и др. задач. Для широкого применения торцевых СУ необходимы знания их метрологических характеристик в широком диапазоне режимных и конструктивных параметров. В настоящий момент применение торцевых СУ ограничено отсутствием официальных Правил и методик измерения расхода на их основе.
В классе торцевых СУ для измерения расхода на входе в ИТ наибольшее распространение имеют стандартные диафрагмы. Результаты первых исследований коэффициента расхода торцевых диафрагм были опубликованы в 1934 году Штахом [1]. Независимо от относительного диаметра получено, что коэффициент расхода при Ре > 5,5-104 составляет величину 0,6. Другие исследования [2-4] подтвердили данное значение. Отсутствие зависимости коэффициента расхода от числа Рейнольдса и относительного диаметра свидетельствует о широком и универсальном характере использования торцевой диафрагмы.
В условиях стандартного монтажа сопло по сравнению с диафрагмой обладает более стабильными характеристиками при длительной эксплуатации, имеет меньшие потери давления и в меньшей степени реагирует на турбулентные пульсации потока, но неопределенность коэффициента истечения больше, чем у диафрагмы [5]. Во многом достоинства сопла определяются достаточно полным сжатием струи, что обеспечивает профиль сопла, и меньшим вихреобразованием за соплом [6].
Для стандартного сопла ИСА 1932 существует зависимость коэффициента истечения от числа Рейнольдса и относительного диаметра в = СС/й [7], где СС, й - диаметры отверстия СУ и ИТ, соответственно. Предположение об отсутствии качественных отличий в структуре потока за СУ в условиях стандартного монтажа и при торцевом исполнении дает основание для существования зависимостей коэффициента истечения торцевого сопла от Ре и от в.
В отношении метрологических
характеристик торцевых сопел данные, представленные в литературе, весьма скудные. Так, Штах для торцевых сопел без входного участка с ввых = 0,4^0,7 (индекс «вых» означает, что относительный диаметр рассчитывался по диаметру ИТ на выходе СУ, т.к. на входе й ^ ж и в ^ 0) при Ре > 5,5-104 независимо от относительного
диаметра получил значение коэффициента расхода 0,99. В более поздних работах [8] с погрешностью 0,6 % было получено значение коэффициента расхода 0,992.
Анализ современного состояния исследований метрологических характеристик торцевых сопел без входного участка позволяет сделать следующие выводы:
- коэффициент расхода практически не зависит от числа Рейнольдса, что вызывает сомнения при сопоставлении зависимостей коэффициента истечения стандартного сопла;
- отсутствует зависимость коэффициента расхода от геометрических параметров сопла.
Цель работы
Изложенное выше определило цель работы - исследовать метрологические характеристики торцевого сопла. Достижение данной цели предполагает выполнение следующих задач:
- проведение экспериментальных и численных исследований коэффициента истечения торцевых сопел с забором из безграничного объема (отсутствует входной трубопровод);
- численное моделирование коэффициента истечения сопла ИСА 1932 в условиях стандартного монтажа [7] в целях верификации расчетных моделей.
Экспериментальная часть
Экспериментальное определение
коэффициента истечения торцевого сопла проводилось на государственной эталонной установке Эу-2 Государственного первичного эталона массового и объемного расходов газа ГЭТ-118-2006 ФГУП ВНИИР с расширенной стандартной неопределенностью < 0,094 %.
В качестве торцевого сопла исследовалось сопло ИСА 1932 для ИТ диаметром 200 мм. Для проведения экспериментов были изготовлены 2 торцевых сопла согласно требованиям ГОСТ 8.586.3-2005 [7] с ^вых= 0,6 и 0,75. Исследуемые торцевые сопла с угловым отбором давления через кольцевую камеру осреднения монтировались на трубопроводе диаметром 200 мм и длиной 4 м, что превышало минимально допустимую длину выходного участка после СУ [7]. Подвод воздуха к торцевому СУ соответствовал естественным нестесненным условиям, что обеспечивалось отсутствием препятствий, нарушающих движение воздуха на расстояниях не менее 20с1 перед СУ по направлению оси и не менее 10С перпендикулярно оси от СУ [1].
В процессе экспериментальных
исследований набором калиброванных критических сопел задавался эталонный расход газа Цтэт, измерялись давление р, температура Т, влажность поверочной среды в невозмущенных условиях помещения и перепад давления Др на исследуемом торцевом сопле. По результатам измерений давления, температуры и влажности косвенным методом определялась плотность среды р.
Уравнение для определения коэффициента истечения имеет вид:
С = Чтэт/П С 2 , (1)
где е - коэффициент расширения, учитывающий сжимаемость среды, который рассчитывался по формуле
£ =
'кт 2 к
1 - т
(к-1)1 к \
к -1
1-т
;
здесь т = 1 - Др/ р .
Неопределенность коэффициента истечения зависит от точности измерений.
Для обеспечения точности задаваемого расхода дт,эт проводилась калибровка критических сопел по методике ВНИИРа [9].
Перепад давлений на торцевом сопле измерялся преобразователем перепада давлений «Yocogawa» Б1Л 110А с пределом приведенной погрешности + 0,15%. Верхний предел измерений перепада давлений 10 кПа. С уменьшением перепада давления пределы допускаемой основной приведенной погрешности увеличиваются. В нижнем диапазоне 130 - 150 Па относительная погрешность преобразователя достигает 5%.
Для повышения точности измерений преобразователь перепада давления подвергался индивидуальной калибровке с помощью калибратора давления Метран-517. Для проведения калибровки использовался эталонный модуль Метран-518 с диапазоном измерения давлений 0 -10 кПа и предел допускаемой приведенной погрешности + 0,02%. Результаты калибровки показали, что погрешность преобразователя перепада давления не выходит за пределы погрешности калибратора. В расчетах погрешности преобразователя перепада учитывались как случайные, так и инструментальные погрешности.
Точность измерения температуры и влажности обусловлена параметрами
преобразователя ИВТМ-7: в диапазоне от -20°С до 60°С абсолютная погрешность измерения температуры ДТ = 0,1 °С; абсолютной погрешность измерения влажности в диапазоне от 2 до 98 % составляла Дф = 1 %. Методическая погрешность определения давления насыщенного водяного пара не превышала 0,05 %.
Неопределенность коэффициента истечения согласно ГОСТ 8.586.5 рассчитывалась как неопределенность измерений величин, связанных с коэффициентом истечения функциональной зависимостью (1), по формуле
ис =
2 Ы )
где
и', - коэффициент влияния
неопределенность параметра х,, входящего в уравнение (1). В свою очередь, коэффициенты влияния рассчитывались следующим образом
,=1
и
дС х/
' дх1 С
Расчеты показали, что в диапазоне от 180 до 3500 Па по перепаду давления на сопле инструментальная погрешность не превышает 0,61 %.
Результаты измерений и проведенные расчеты показали, что случайная составляющая для различных диапазонов перепада давления изменялась от 0,02% (при больших перепадах) до 0,1% (при перепадах около 120 Па). Эти значения составляют менее 1/3 от величины инструментальной погрешности преобразователя перепада.
Рассчитанные значения неопределенности коэффициента истечения не превышают неопределенность, определяемую стандартом.
Результаты экспериментальных исследований
Исследования коэффициента истечения торцевых сопел проведены в диапазоне 3-104 < Рей < 4,5-105, рассчитанного относительно диаметра ИТ.
При измерениях в низком диапазоне чисел Рейнольдса возникали трудности с точностью измерений малых перепадов давления. В свою очередь ограничения по наибольшему значению
определялись возможностями подачи эталонного расхода поверяемой среды.
Экспериментальные зависимости С= f (Рей) приведены на рис.1. Полученные значения коэффициента истечения торцевого сопла в диапазоне Рей > 105 соответствуют
экспериментальным данным [8]. Показано, что коэффициент истечения торцевого сопла увеличивается с ростом Ре и зависит от ввых. Неопределенность коэффициента истечения для
сопла с ввых= 0,75 выше, чем для сопла с j
х= 0,6.
Рис. 1 - Экспериментальные зависимости коэффициента истечения торцевых сопел: х - А,ых= 0,6; о - ввых= 0,75
При перепадах Др > 200 Па,
соответствующих Рей > 105, расширенная неопределенность < 1 % и зависимостям С = f (Рей) свойственно монотонное увеличение, а при Др < 100 Па (Рей < 105) неопределенность
коэффициента истечения > 1 % и монотонный характер зависимостей С= f (ReD) отсутствует.
При ReD > 105 максимальные отклонения предельных экспериментальных значений коэффициента истечения относительно их среднего значения < 1 %, что хорошо согласуется с расчетами неопределенности. С ростом числа Рейнольдса разброс экспериментальных значений коэффициента истечения, как и его неопределенность, монотонно убывает.
Численные исследования
Экспериментальное определение
коэффициента истечения в широком диапазоне режимных и конструктивных параметров торцевых сопел требует существенных материальных затрат. Особенно проблемно получить экспериментальные данные для торцевых сопел при высоких числах Рейнольдса. Численные методы, позволяющие корректно рассчитывать коэффициент истечения, могут позволить сократить объем
экспериментальных работ и способствовать получению инженерных зависимостей для расчета коэффициента истечения.
Применение численных методов для расчета коэффициента истечения стандартной диафрагмы, как показано в [10,11], обеспечивает точность определения C в пределах систематической погрешности.
Апробация применения численных методов к расчету коэффициента истечения торцевых сопел проведена на примере расчета течения на участке ИТ с соплом ИСА 1932 в условиях стандартного монтажа [7]. Рассматривался ИТ диаметром D =200 мм и сопло ИСА 1932 с относительным диаметром в диапазоне 0,4 < в < 0,6. Расчеты выполнялись с применением программного обеспечения ANSYS Fluent.
Моделирование проводилось в
осесимметричной постановке в рамках изотермического несжимаемого течения. Геометрия расчетной области отражала все конструктивные особенности монтажа и обеспечивала отсутствие влияния граничных условий на результаты расчета. Сетка создавалась с выделением пограничных слоев, что обеспечивало возможность сеточного разрешения вязкого подслоя. Расчеты проводились по двум моделям турбулентности - k-s RNG и k-a SST со стандартным набором эмпирических констант. Коэффициент истечения сопла ИСА 1932 при стандартных условиях монтажа рассчитывался по формуле
' п
С = цт - С2 (У V 1-в4 Jл/2Дpp .
Результаты моделирования показали возможность рассчитывать коэффициент истечения в пределах неопределенности, установленной
стандартом (УС = 0,8 % для в < 0,6). Независимо
от в при Рей > 3-104 относительные отклонения расчетных значений коэффициента истечения от данных [7] не превышают 0,7 %. В качестве примера
на рис.2. приведены результаты расчета для сопла с в = 0,6.
Следует отметить, что, например, для в =0,4 модели турбулентности к-а 88Т и к-е Я^в практически одинаково прогнозируют коэффициент истечения (отличия не превышают 0,5 %). С увеличением в появляются определенные отличия в структуре потока по рассматриваемым моделям турбулентности и соответственно растут различия в коэффициенте истечения. Для в =0,6 при Рво < 2-104 расчетные значения коэффициента истечения выходят за рамки допустимой стандартом неопределенности.
С
0,96 0.95 0,94 0,93 0,92
---
104
105
Rer
Рис. 2 - Сопоставление расчетных зависимостей С= f (Reo) сопла ИСА 1932 при ß =0,6 в условиях стандартных монтажа с данными стандарта
(_) и допускаемой неопределенностью
(_______) [7]: ...... - модель турбулентности к-ю
SST;__- k-sRNG
Зависимости С = f (Reo) торцевого и стандартного сопел отличаются между собой. Очевидно, что зависимости коэффициента истечения стандартного сопла ИСА 1932 не пригодны для оценки значений коэффициента истечения торцевого сопла.
Целью численных исследований течения в ИТ с торцевым соплом был выбор модели турбулентности, обеспечивающей лучшее согласование результатов расчета с экспериментальными данными.
Для организации численных исследований использовалась постановка задачи, аналогичная расчету сопла в стандартных условиях. Расчетная область строилась в соответствии с условиями эксперимента. Граница расчетной области на входе располагалась по радиусу на расстояние »10d от оси сопла, что исключало ее влияние на результаты расчета. Сетка строилась с необходимым разрешением вязкого подслоя. Для проведения расчетов использовались те же 2 модели турбулентности - k-sRNG и k-ю SST.
Расчеты со стандартным набором констант в выбранных моделях турбулентности показали, что существует расхождение расчетных и экспериментальных данных (рис. 3). Качественный характер полученных численных зависимостей С = f (Reo) соответствует экспериментальному.
Модели турбулентности k-s RNG и k-ю SST практически одинаково прогнозируют коэффициент истечения торцевого сопла. Существенные отличия в результатах расчета между моделями турбулентности k-s RNG и k-ю SST наблюдаются для сопла с ßвых =0,6 в диапазоне Reo < 3-104. В других случаях отличия не превышают 0,4 %.
Наибольшие расхождения расчетных и экспериментальных данных существует в диапазоне Reo <105 и составляет > 2 %. Следуют отметить, что при Reo <105 экспериментальные данные получены с неопределенностью, превышающую допустимую, о чем свидетельствует разброс результатов измерений. С увеличением Re расхождение уменьшается, но сохраняется до 1,5 %.
С
0,995 0,990 0,985 0,9 S0 0,975 0,970 0,965 0,960
-
•
****** * • • • • *
> .7
/ г / ß =с В Ы >: .6
/ ; /
104
105
Rer
Рис. 3 - Сопоставление экспериментальных и расчетных зависимостей С= f (Reo) торцевых
сопел:.........- модель турбулентности к-ю SST;
_- модель k-s RNG ; - модель k-s
RNG с модифицированным набором констант
Полученное расхождение расчетных и экспериментальных данных может быть вызвано многими причинами, одна из которых, неравномерность профиля скорости по соплу. Выравнивание профиля скорости в математической модели можно обеспечить за счет повышения интенсивности турбулентности, что возможно за счет корректировки значений констант в моделях турбулентности.
Модель k-s RNG является более привлекательной для исследований влияния значений модельных констант на величину коэффициента истечения, т.к. без учета эффектов вращения имеет всего три модельных константы, а модель k-ю SST - пять. Кроме этого, при Reo ^ 105 рассматриваемые модели турбулентности
практически одинаково прогнозируют коэффициент истечения.
В модели Я^в используются константы Си , С1е и С2е . Константа Си необходима для расчета турбулентной вязкости / ¡, которая в пределе Ре ^ ж имеет значение 0,0845, а в низкорейнольдсовой области определяется расчетным путем. Константы С1е и С2е соответственно используются в членах, отвечающих за генерацию и диссипацию кинетической энергии турбулентности (стандартные значения констант: С1е=1,42; С2е=1,68). При заданном перепаде давления последовательно исследовалось влияние той или иной константы при постоянных значениях других. Влияние констант на коэффициент истечения исследовалось для ИТ с различными ввых торцевых сопел в зависимости от Ре.
Исследования показали, что максимальные значения коэффициента истечения торцевого сопла получаются при значениях констант См = 0,16, С\Е= 1 и С2е= 2. Модификация констант приводит к снижению расхождения между расчетом и экспериментом. После корректировки констант в диапазоне Рей > 105 расхождение результатов расчета с экспериментальными данными не превышает 0,5 % для сопла с ввых = 0,6 и 0,9 % для сопла с ввых = 0,75. Сопоставления расчетных зависимостей С= f (Рей) с откорректированными значениями констант и со стандартным набором модельных констант отражены на рис.3.
Заключение
По результатам проведенных исследований можно сделать следующие выводы.
1. Получены экспериментальные зависимости коэффициента истечения торцевых сопел от числа Рейнольдса и от относительного диаметра.
2. Численные исследования коэффициента истечения стандартного сопла ИСА 1932 показали адекватную возможность расчета коэффициента истечения в пределах неопределенности, установленной стандартом.
3. Расчетные значения коэффициента истечения торцевого сопла по моделям турбулентности со стандартным набором модельных констант в диапазоне Рей > 105 дает
расхождение с экспериментальными данными до 1,5 %.
4. Модификация констант в моделях турбулентности приводит к уменьшению расхождения между результатами численного и экспериментального исследований.
Литература
1. П.П. Кремлевский, Расходомеры и счетчики количества: Справочник. Л.: Машиностроение, 1989.701с.
2. В.Ф. Рис, Н.А. Широков, Энергомашиностроение, №10, 13-14 (1984)
3. ISO/DIS 5221. Draft International standard. Guide to methods of measuring air flow rate in an handling duct. 1980.
4. Р.С. Гареев, И.А. Быков, В.А. Фафурин, В.Б. Явкин, И.Ф. Яценко, Вестник Казан. технол. ун-та, т.16, №17,.197-200 (2013)
5. ГОСТ 8.586.1-2005 (ИСО 5167-1:2003). Измерение расхода и количества жидкостей и газов с помощью стандартных сужающих устройств. Часть 1. Принцип метода измерений и общие требования.
6. В.П. Преображенский, Теплотехнические измерения и приборы: Учебник для вузов по специальности «Автоматизация теплотехнических процессов». М.: «Энергия», 1978. 704 с.
7. ГОСТ 8.586.3-2005 (ИСО 5167-3:2003). Измерение расхода и количества жидкостей и газов с помощью стандартных сужающих устройств. Часть 3. Сопла и сопла Вентури. Технические требования.
8. E. Eujen, O. Burcke, BWK, Bd.12, №10, 63-64 (1960)
9. Методика калибровки эталонных критических сопел при атмосферном давлении на исходной эталонной установке Государственного первичного эталона единиц объемного и массового расходов газа ГЭТ 118-2006. МК 2567988-03-2008. Федеральное государственное унитарное предприятие «Всероссийский научно-исследовательский институт расходометрии». Казань, 2008 г., 19 с.
10. Р.И. Ганиев, Н.А. Николаев, А.Н. Сабирзянов, В.А. Фафурин, В.Б. Явкин, Изв. Вузов. «Авиационная техника», №4, 21-24 (2008)
11. Р. А. Тырышкин, А.Н. Сабирзянов, В. А. Фафурин, В.В. Фефелов, В.Б. Явкин, Вестник Удмуртского университета. Математика. Механика. Компьютерные науки. Вып.2, 109-116 (2010)
© И. А. Яценко - технич. дир., Метрологический центр ООО "СТП"; office@ooostp.ru; А. Н. Сабирзянов - к.т.н., доц. каф. реактивных двигателей и энергетических установок КНИТУ им. А.Н.Туполева, ANSabirzyanov@kai.ru; В. А. Фафурин -д.т.н., проф., зав. каф. систем автоматизации и управления технологическими процессами КНИТУ, fafurin_va@kstu.ru; В. Б. Явкин - к.т.н., доц. каф. реактивных двигателей и энергетических установок КНИТУ им. А.Н.Туполева, vd.yavkin@gmail.com.
