CC BY
87
ПРОЦЕССЫ И АППАРАТЫ ХИМИЧЕСКОЙ ТЕХНОЛОГИИ
УДК 681.121.842+532.5
Р. С. Гареев, И. А. Быков, В. А. Фафурин,
В. Б. Явкин, И. А. Яценко
ЭКСПЕРИМЕНТАЛЬНОЕ И РАСЧЕТНОЕ ИССЛЕДОВАНИЕ ТОРЦЕВЫХ ДИАФРАГМ
Ключевые слова: Торцевая диафрагма, коэффициент истечения, эталонная установка, математическое моделирование, модели турбулентности, коэффициент сужения струи, коэффициент гидравлических потерь.
Экспериментальным и расчетным путем получены зависимости коэффициента истечения трех типов торцевых диафрагм в диапазоне чисел Рейнольдса Re=2,5x104 -4х105 и относительного диаметра диафрагмы /3=0,4 - 0,75. Эксперименты были проведены на эталонной установке ФГУП «ВНИИР». На основе численного решения уравнений для стационарного турбулентного течения несжимаемой жидкости проведены расчеты коэффициента истечения торцевой диафрагмы. Для получения большей точности при расчете с использованием моделей турбулентности k-sRNG и k-m SST была проведена модификация стандартного набора констант моделей. Проведен анализ уравнения для коэффициента истечения, показано, что определяющее влияние на коэффициент истечения оказывают только коэффициенты сужения струи и гидравлических потерь.
Keywords: End diaphragm, discharge coefficient, reference installation, mathematical modeling, turbulence model, the coefficient of
narrowing jet, the coefficient of hydraulic losses.
Experimental and calculated dependences of the coefficient expiration of the three types of end diaphragms were obtained in the range of Reynolds numbers Re=2,5x104 - 4x105 and the relative diameter diaphragm /3=0,4 - 0,75. The experiments were performed on a reference installation ofFGUP VNIIR. The discharge coefficient of end diaphragm was calculated on the basis of numerical solutions of equations for stationary turbulent flow of an incompressible fluid.
To achieve greater accuracy in the calculation with the use of models of turbulence k-sRNG and k-m SST constant modification of a standard set of models was held. The analysis of the equation for the discharge coefficient was held and it is shown that the determining influence on the discharge coefficient have only narrowing jet and hydraulic losses
coefficients.
В практике измерения расходов приходится иметь дело с нестандартными применениями приборов переменного перепада давления. К такому случаю относится использование торцевых диафрагм, отличающихся от стандартной диафрагмы отсутствием входного, выходного или входного и выходного трубопроводов. Торцевые диафрагмы могут применяться в условиях ограничений на длины трубопроводов. Их применение имеет определенные преимущества при заборе из большого объема или сбросе в большой объем. Для возможности практического применения этих средств измерений необходимо знать их метрологические характеристики в широком диапазоне режимных и конструктивных параметров.
В отличие от стандартных диафрагм, в которых обязательно присутствие длинного входного участка, необходимого для формирования профиля скорости полностью развитого турбулентного потока и выходного участка, в торцевых диафрагмах отсутствует выходной или входной участки либо оба участка. Будем обозначать в дальнейшем диафрагму без входного участка Д1, без выходного участка Д2, без входного и выходного участков Д3. В литературе сведения об использовании и метрологических характеристиках торцевых диафрагм представлены крайне скудно. Исследования проводились только для диафрагм без входного участка. Согласно данным, представленным в [1] коэффициент расхода торцевых диафрагм близок к 0,6 при
Яе>5,5х10 и относительном диаметре /=0,5 - 0,8 и не зависит от / Результаты исследований, представленных в [2] подтверждают эти данные. В [3] приведено эмпирическое уравнение для определения коэффициента истечения диафрагмы без входного участка
ч0,7
(1)
Это уравнение гарантирует определение С с неопределенностью не превышающей 1% и справедливо при условии Б>2ё, где Б - диаметр трубопровода выходного участка, ё - диаметр диафрагмы. Недостатком уравнения является отсутствие зависимости коэффициента истечения от отношения ё/Б. При установке торцевой диафрагмы должно быть выполнено условие удаленности стесняющих поток поверхностей, которое по данным [1] составляет не менее 10ё и не менее 8ё по данным [3]. Отсюда следует, что коэффициент истечения торцевой диафрагмы должен совпадать с коэффициентом истечения стандартной диафрагмы при /<0,125.
Для диафрагм без выходного участка в [3] предлагается применять эмпирическое уравнение для коэффициента истечения стандартной диафрагмы (см. [5]), увеличивая неопределенность до 1%. В работе [3] нет формулы для расчета коэффициента истечения для Д3.
Получение более точных уравнений, учитывающих изменение геометрических и режимных
параметров в более широких пределах, потребует проведения большого объема дорогостоящих экспериментальных исследований. Однако, с учетом развития в настоящее время программных комплексов по моделированию течений жидкостей и газов, значительный объем работ может быть проведен путем численного моделирования, с последующей проверкой в относительно небольшом по объему экспериментальном исследовании.
В данной работе представлены результаты экспериментального и расчетного исследования коэффициента истечения торцевой диафрагмы с забором из безграничного объема (отсутствует входной трубопровод), полученные в диапазоне изменения
4 5
числа Рейнольдса Яе=2,5х10 - 4х10 и относительного диаметра диафрагмы /=0,4 - 0,75. Экспериментальное определение коэффициента расхода проводилось на государственной эталонной установке ЭУ-2 Государственного первичного эталона единиц объёмного и массового расходов газа ГЭТ 118-2006, с расширенной стандартной неопределенностью при 0,092% (при к=2) [6][7].
Результаты экспериментов представлены на рис.1 в виде зависимости коэффициента истечения от числа Рейнольдса. На этом же рисунке нанесены зависимости коэффициента истечения стандартной диафрагмы, полученные расчетом по эмпирическому уравнению [5] при /=0,1 и 0,75. С ростом числа Рейнольдса у всех испытанных диафрагм коэффициента истечения монотонно уменьшается. По характеру зависимости от числа Рейнольдса полученные результаты качественно соответствуют зависимостям коэффициента истечения стандартных диафрагм, однако количественное расхождение достаточно велико и в некоторых случаях превышает 1,5%. Присутствует расслоение экспериментальных данных по относительному диаметру диафрагмы /=ё/Б. Расчетная зависимость коэффициента истечения стандартной диафрагмы, полученная при /=0,1, дает заниженные на 1 - 1,5% по сравнению с экспериментом значения С при Яе<2х10 и на 0,9 -1% при Яе>2х10 ; расчет, выполненный для /=0,75, наоборот, дает завышенные значения С при малых числах Яе и неплохо согласуется с экспериментальными для диафрагм с /=0,6 и 0,75 при Яе>1х10 . Уравнение (1) хорошо согласуется (в пределах 0,5%) с результатами экспериментов для /=0,4, а также для /=0,6 при Яе<2х10 . С ростом числа Рейнольдса расхождение формулы (1) и данных для /=0,6 возрастает, а при /=0,75 расхождение формулы (1) с данными эксперимента достигает 1% и более.
Из представленных результатов можно сделать вывод, что эмпирическое уравнения для коэффициента истечения стандартной диафрагмы [5] не пригодно для определения коэффициента истечения торцевой диафрагмы, а уравнение (1) хорошо согласуется с экспериментальными данными при малых / и числах Рейнольдса Яе<2х10 .
Для выяснения возможности применения численных методов к расчету коэффициента истечения торцевой диафрагмы были проведены расчеты, выполненные в осесимметричной постановке.
Геометрия расчетной области была выбрана такой, чтобы передать все конструктивные особенности монтажа диафрагмы (фланцы, крепление, отбор давления). Расчетная сетка содержала 250 000 прямоугольных ячеек, ориентированных вдоль стенок, с выделением пограничных слоев и плавным увеличением размера ячейки к ядру потока. Сгущения сетки в пристеночных областях были такими, что величина у+^у/ V < 1,5 для 99% пристеночных
* I-------
ячеек, где u = л тw/Р , - напряжения трения на
стенке, у - расстояние от центра граничной ячейки до стенки. Это обеспечивало возможность сеточного разрешения ламинарного подслоя. В окрестности угловой точки диафрагмы величина у+ превышала указанное выше значение и могла достигать 5.
Red*1 0-6
О /=0,75 О /=0,6
Д /=0,4 данные [3]
Стандарт, Ь=0,1 Стандарт, Ь=0,75
Рис. 1 - Зависимость коэффициента истечения торцевой диафрагмы от числа Рейнольдса
Основной целью расчетов был выбор модели турбулентности и модификация констант модели для получения лучшего согласования с экспериментальными данными. Основываясь на результатах работы [8] для проведения расчетов были выбраны 2 модели турбулентности - к-е и к-ю 88Т.
Расчеты, проведенные со стандартным набором констант в модели турбулентности к-ю 88Т, показали, что полученная зависимость коэффициента истечения от числа Рейнольдса качественно согласуется с экспериментальными данными, однако имеется значительное постоянное смещение расчетных данных, около 1%, практически одинаковое для всех / (рис.2). Такой же результат был получен и для модели турбулентности к-е со стандарт-
ным набором констант. Эта модель также дает значительное смещение относительно экспериментальных точек. Устранение смещения было достигнуто корректировкой констант модели турбулентности, как для модели к-е Я^в, так и для модели к-ю 88Т.
Так как смещение, полученное при расчетах приводит к заниженным значениям С, вызванным, в том числе, большей неравномерностью профиля скорости перед диафрагмой, то корректировка была выполнена таким образом, чтобы обеспечить большие значения интенсивности турбулентности вблизи входа в диафрагму, что должно привести к выравниванию профиля скорости. После соответст-
вующего изменения констант в моделях турбулентности постоянное смещение экспериментальных и расчетных данных было сокращено до 0,5%. Из представленных на рис.2 графиков видно, что после модификации констант расчетные и экспериментальные данные согласуются между собой значительно лучше, чем при стандартном наборе констант. Отметим также, что результаты расчетов по модели к-ю 88Т с модифицированными константами хорошо согласуются с расчетами по формуле (1).
Red *10-6
О Эксперимент, b=0,6 — — Расчет, k-e RNG стандартные константы
Расчет, k-e RNG модифицированные константы
Ке/Ю'6
Д Эксперимент, Ь=0,4
— - Расчет, к-м' БвТ стандартные константы
Расчет, к-м ББТ модифицированные константы
— Расчет по формуле [3]
Рис. 2 - Результаты расчетов С численным методом
Для оценки влияния типа диафрагмы на коэффициент истечения были проведены численные расчеты течения в диафрагмах трех типов, а также в стандартной диафрагме с использованием модели турбулентности кю-ББТ со стандартным набором констант. По результатам расчетов были построены зависимости коэффициента истечения, С от числа Рейнольдса, представленные на рис.3.
Из результатов расчетов следует, что характер изменения С в зависимости от числа Рейнольдса одинаков для всех типов диафрагм - с ростом Яе коэффициент истечения монотонно убывает, но величина С зависит от типа диафрагмы и относитель-
ного диаметра (за исключением Д3), при этом, различие может достигать 2%.
ЯеЛО*
—без выходного участка :
-А- без входного участка ;
без входного и выходного участка :
А стандартная диафрагма.
Рис. 3 - Зависимости коэффициента истечения торцевых диафрагм от числа Рейнольдса по результатам численных расчетов
Для объяснения полученных результатов рассмотрим уравнение коэффициента истечения, полученное на основе одномерного подхода [4]
с- , (2)
л/аЬ-фЬ + ^вЬ-К-^8 )^4
где /и - коэффициент сужения струи, а - коэффициент Кориолиса, учитывающий неравномерность профиля скорости, у - коэффициент отбора давления, учитывающий разность давлений в точке отбора давления и в горле струи, £аЬ - коэффициент гидравлических потерь, индекс а обозначает сечение в невозмущенном потоке перед диафрагмой, индекс Ь
- сечение в горле струи за диафрагмой.
Эта формула получена для стандартной диафрагмы, но справедлива для всех типов диафрагм. Для диафрагм без входного и без входного и выходного участков коэффициенты аа=0, уа=0, т.к. сечение а—а располагается на большом удалении от диафрагмы, где скорость течения пренебрежимо мала. Такой же результат получается и формально, если принять, что поток забирается из бесконечно большого объема, вследствие чего /=0, т.е. (аа-уа)и /=0. Для диафрагмы без выходного участка влияние коэффициентов аа=0, уа=0 сильно ослабляется, т.к. произведение и /<0,05 при /<0,6.
Наиболее сильно влияющим параметром является коэффициент сужения струи, и Величина коэффициента сужения определялась по положению «граничной» линии тока, которая для каждой диафрагмы подбиралась, как линия тока наиболее близко подходящая к угловой точке диафрагмы и являлась, тем самым, условной границей струи. На этой линии тока определялось положение минимальной радиальной координаты, по которому рассчитывал-
ся коэффициент сужения струи. Для диафрагмы без выходного участка (Д2) получено наибольшее значение коэффициента сужения струи ц=0,677. У двух других типов диафрагм коэффициенты сужения струи одинаковы, ц=0,63. Большие значения и для диафрагмы Д2 обусловлены тем, что поток на входе в диафрагму в большей степени ориентирован вдоль оси диафрагмы, в то время как у диафрагм Д1 и Д3 картина линий тока близка к стоку. Таким образом, большая величина коэффициента истечения диафрагмы без выходного участка объясняется большей величиной коэффициента сужения струи.
Однако, несмотря на то, что у диафрагм Д1 и Д3 коэффициент сужения струи одинаков, коэффициент истечения немного различается. Более того, приблизительно такое же различие отмечается и в коэффициентах истечения стандартной диафрагмы и Д2. Для объяснения этих различий рассмотрим влияние коэффициентов аЬ, £аЬ, щ.
Сравнение коэффициентов аЬ и щЬ, выполненное на основе результатов расчетов, показало, что различия между этими коэффициентами для разных типов диафрагм достаточно малы и не оказывают заметного влияния на коэффициент истечения. Так, значения аЬ изменяются от 1,028 до 1,032, а щЬ от 0,072 до 0,075, а влияние этих изменений не превышает 0,2%.
Более значительное влияние оказывает коэффициент гидравлических потерь, характеризующий потери полного давления в струе между сечениями а и Ь. Как отмечается в [1] коэффициент гидравлических потерь определяется потерями в отрывной зоне перед диафрагмой, потерями на трение и потерями при обтекании диафрагмы; при этом отрывная зона за диафрагмой не влияет на потери давления. Из сравнения зависимостей коэффициента истечения различных типов диафрагм, приведенного на рис.3 следует, что большие значения С получаются для диафрагм, у которых отсутствует выходной участок, т.е. отсутствует отрывная зона за диафрагмой. Таким образом, наименьшее значение коэффициента гидравлических потерь соответствует диафрагмам, у которых отсутствует выходной участок, следовательно, влияние отрывной зоны перед диафрагмой пренебрежимо мало, а потери в отрывной зоне за диафрагмой имеют большее значение.
По результатам проведенного исследования могут быть сделаны следующие выводы:
1. Коэффициент истечения торцевой диафрагмы зависит от типа и относительного диаметра диафрагмы (при наличии входного или выходного участков).
2. Наиболее сильное влияние на коэффициент истечения оказывает коэффициент сужения струи.
3. При одинаковых значениях коэффициента сужения струи разные значения коэффициента истечения получается вследствие различия в коэффициентах гидравлических потерь.
Литература
1. П. П. Кремлевский, В кн. Расходомеры и счетчики количества веществ. Справочник: Кн. 1.- 5е изд. перераб. и доп., Политехника, Санкт-Петербург, 2002, С. 409 с.
2. В. Ф. Рис, Н. А. Широков, В Кн. Коэффициенты расхода торцевых диафрагм при нестесненном входе воздуха, Энергомашиностроение, 10, C. 13-14.
3. ISO Technical Report. ISO/TR 15377:2007 (E). Measurement of fluid flow by means of pressure differential devices
- Guidelines for the specification of orifice plates, nozzles and Venturi tubes beyond the scope of ISO 5167.
4. ГОСТ 8.586.1-2005 (ИСО 5167-1:2003). Измерение расхода и количества жидкостей и газов с помощью стандартных сужающих устройств. Часть 1. Принцип метода измерений и общие требования.
5. ГОСТ 8.586.2-2005 (ИСО 5167-2:2003). Измерение расхода и количества жидкостей и газов с помощью стандартных сужающих устройств. Часть 2. Диафрагмы Технические требования.
6. Р. С. Гареев, А. В. Красавин, В. А. Фафурин, В. Б. Яв-кин, И. А. Яценко, Экспериментальное и расчетное определение коэффициента истечения торцевой диафрагмы, Законодательная и прикладная метрология, 6, С. 1317, (2012).
7. Красавин В.М., Раинчик С.В., Ившин В.П., Красавин А.В., Быков И.А., Лобарев С.В. Исходная эталонная установка государственного первичного эталона единиц объёмного и массового расходов газа ГЭТ 11-2006. Процедура калибровки критических сопел// Вестник Казан. технол. ун-та, №10, c.285-287 (2013).
8. Ганиев Р. И., Николаев Н. А., Сабирзянов А. Н., Фафурин В. А., Явкин В. Б., Выбор сетки и модели турбулентности для расчета коэффициента расхода стандартной диафрагмы, Изв.Вузов. «Авиационная техника», 4, С. 21-24 (2008).
© Р. С. Гареев - зам. нач. отдела разработки методик (методов) измерений, Метрологический центр ООО “СТП”, [email protected]; И. А. Быков - науч. сотр. ФГУП «ВНИИР», [email protected]; В. А. Фафурин - д-р техн. наук, проф., зав. каф. систем автоматизации и управления технологическими процессами КНИТУ; В. Б. Явкин, канд. техн. наук, доц. КНИТУ им. А.Н.Туполева - КАИ (КНИТУ-КАИ); И. А. Яценко - техн. дир., Метрологический центр ООО “СТП”, [email protected].
