СПОСОБ ОПРЕДЕЛЕНИЯ РАСХОДНОЙ ХАРАКТЕРИСТИКИ ГИДРАВЛИЧЕСКОГО ТРАКТА В ПЕРЕХОДНОЙ ОБЛАСТИ ИСТЕЧЕНИЯ Российский патент 2022 года по МПК G01M10/00 

Описание патента на изобретение RU2765801C1

Изобретение относится к области измерительной технике и может быть использовано для расчета пропускной способности проектируемых гидравлических трактов транспортных и дозирующих систем в химической, нефтехимической, авиационной, текстильной, лакокрасочной и других отраслях промышленности, в частности, - для определения (прогнозирования) расходных характеристик систем автоматизированного порционного дозирования натурных жидкостей с поддержанием постоянного перепада давлений на линии налива, с отмериванием дозы по времени.

Основной задачей, которую приходится решать при проектировании тракта транспортирования жидких продуктов, является определение его пропускной способности (объемного расхода) Q при различных значениях действующего (располагаемого) напора Н, т.е. его расходной характеристики - функции Q=F(H). Предлагаемый способ определения расходной характеристики основан на анализе уравнения Бернулли для установившегося течения несжимаемой жидкости. Это уравнение для гидравлического тракта, состоящего из нескольких местных сопротивлений, соединенных участками трубопровода длиной Li с внутренним диаметром di, с концевой раздачей через запорное устройство с цилиндрическим сливным наконечником с внутренним диаметром du имеет вид:

где - скорость истечения жидкости из сливного наконечника (ƒH=πdH2/4 - площадь проходного сечения сливного наконечника), g - ускорение силы тяжести, α - коэффициент Кориолиса потока на выходе сливного наконечника, λi - коэффициент сопротивления трения на i-м участке трубопровода, ξ - суммарный коэффициент потерь напора на местных сопротивлениях в тракте, который определяется конструктивными параметрами местных сопротивлений и зависит от числа Рейнольдса потока (Re) на выходе сливного наконечника.

При заданных значениях параметров νH, Li, di и dH и известной кинематической вязкости натурной жидкости - ν легко рассчитать входящие в уравнение (1) коэффициенты трения λi. Для течения с числами Рейнольдса Rei<104 на отдельных участках трубопровода для определения λi используется формула ламинарного режима течения где - скорость течения на i-м участке трубопровода. Таким образом, для определения интересующей нас функции Q=F(H) необходимо определить входящие в (1) заранее неизвестный коэффициент потерь напора на местных сопротивлениях ξ и зависящий от характера распределения скоростей на выходе сливного наконечника коэффициент Кориолиса α.

Известен способ определения коэффициента полного гидравлического сопротивления тракта - (RU 2240525 С1, 20.11.2004), который реализуется путем проливки натурного тракта на переходном режиме натурной жидкостью, используемой в технологическом процессе.

Недостатками данного способа является необходимость испытаний тракта на натурной жидкости, что требует дополнительного контроля изменений ее вязкости в процессе измерений, а также ограниченность способа применения на натурных трактах большой протяженности, что требует их размещения на участке испытаний. Недостатками способа являются также невозможность прогнозирования расчетной расходной характеристики тракта с варьируемыми величинами его конструктивных параметров и высокая трудоемкость.

Известен также способ определения расходной характеристики гидравлического тракта и устройство для его осуществления (RU 2709034 С1, 13.12.2019.), принятый за прототип.

Недостатком способа-прототипа является недостаточная точность и достоверность определения и прогнозирования расходной характеристики проектируемого тракта для области перехода от турбулентного к ламинарному режиму истечения из-за использования при расчетах суммарного коэффициента потерь напора на местных сопротивлениях ξ значения коэффициента Кориолиса α=1,05 в уравнении (1), характерного для турбулентного режима истечения.

Задачей настоящего изобретения является повышение точности и достоверности определения расходной характеристики проектируемого тракта при работе на жидкостях различной вязкости (натурных жидкостях), которое достигается за счет учета изменения числовых значений коэффициента Кориолиса α в области перехода от турбулентного к ламинарному режиму истечения (в переходной области), и возможность ее прогнозирования при варьировании заданного диапазона изменения расхода и длины Li образующих тракт участков трубопроводов.

Техническим результатом способа является повышение точности и достоверности определения расчетной расходной характеристики гидравлического тракта в переходной области при транспортировании по нему жидких продуктов различной вязкости и возможность ее прогнозирования при варьировании заданного диапазона изменения расхода и длины образующих тракт участков трубопроводов.

Технический результат предлагаемого способа определения расходной характеристики гидравлического тракта в переходной области истечения достигается тем, что модель тракта с проектными конструктивами местных сопротивлений и участков соединительных трубопроводов длиной с внутренним диаметром di с концевой раздачей через запорное устройство с цилиндрическим сливным наконечником с внутренним диаметром dH проливают водопроводной водой посредством ее многократного порционного дозирования в мерную тару при различных, поддерживаемых на заданном постоянном уровне, величинах действующего напора Hj (где j=1,2…, n=31 - число проливок в каждом цикле проливок), ступенчато уменьшаемых с шагом 10 см от максимального см до промежуточного значения см - и далее, с шагом 2 см, - до нижней границы диапазона измерений - см; с помощью электронных таблиц программы «Microsoft Excel» последовательно вычисляют: объемные расходы соответствующие напорам Hj, где объем дозы Vj регистрируют с помощью взвешивания на электронных весах, а время дозирования tj задают и отмеривают с помощью электронного таймера; скорости истечения из сливного наконечника потери напора - и числа Рейнольдса потока - на его выходе, где коэффициент Кориолиса α первоначально принимают равным 1,05 как при турбулентном режиме истечения, а кинематическую вязкость воды νB рассчитывают по ее температуре в данном цикле проливок; скорости течения - числа Рейнольдса - коэффициенты трения - и потери напора на трение - на отдельных участках трубопровода; суммарные потери напора на трение - на местных сопротивлениях в тракте - и суммарные коэффициенты потерь напора на местных сопротивлениях - полученную в табличном виде экспериментальную расходную характеристику тракта Qj=F(Hj) делят на зоны условно турбулентного, условно переходного и условно ламинарного режимов истечения посредством последовательного сравнительного анализа числовых значений стандартного отклонения ее точек и выбора количеств точек (р≥6 - в верхней области диапазона измерений) и (q≥6 - в нижней области диапазона) с минимальными значениями стандартного отклонения от графика соответствующей этим точкам степенной аппроксимирующей функции; при этом полученные для вышеуказанных зон значения табличных функций (где индексы s=1,2,3 относятся к условно турбулентной, условно переходной и условно ламинарной зонам, соответственно) корректируют с учетом влияния коэффициента Кориолиса α, который увеличивают от значения α=1,05 при турбулентном режиме истечения до α=2 при ламинарном режиме в соответствии с принятой эмпирической функциональной зависимостью α=ϕ(νH) для вычисления скорректированных значений скоростных напоров на выходе сливного наконечника - и соответствующих скорректированных значений табличных функций для каждой из выделенных зон строят графики функций и путем их регрессионного анализа получают уравнения аппроксимирующих функций - определяют точки пересечения графиков: Re1 и Re2 фиксирующих, соответственно, верхнюю и нижнюю границы фактической переходной зоны, для прогнозируемой расходной характеристики; формируют ряд произвольно, с небольшим шагом уменьшаемых значений расхода Qk и - соответствующий ряд скоростей истечения прогнозируемой расходной характеристики натурной жидкости; для каждого значения k и для каждой из выделенных фактических зон режимов истечения натурной жидкости вычисляют соответствующие этим зонам числа Рейнольдса на выходе сливного наконечника - (где ν - кинематическая вязкость натурной жидкости), и по полученным уравнениям аппроксимирующих функций вычисляют суммарные коэффициенты потерь напора на местных сопротивлениях - для прогнозируемой расходной характеристики натурной жидкости вычисляют: скорости течения - числа Рейнольдса - коэффициенты трения - потери напора на трение на отдельных участках трубопровода - и суммарные потери напора на трение - потери напора на выходе сливного наконечника - и на местных сопротивлениях в тракте - расчетные значения располагаемого напора - по полученным в табличном виде для каждой из выделенных фактических зон режимов истечения функциям путем их аппроксимации степенными функциями получают соответствующие графики и уравнения зависимости расхода от располагаемого напора, т.е. определяют для каждой из трех фактических зон режимов истечения искомые расходные характеристики натурной жидкости.

Увеличение коэффициента а от 1,05 до 2 начинается для полученной в цикле проливок экспериментальной расходной характеристики модели гидравлического тракта от верхней границы условно переходной зоны истечения, а для прогнозируемой расходной характеристики - от верхней границы фактической переходной зоны; увеличение коэффициента α до значения α=2 заканчивается на нижней границе ламинарной зоны экспериментальной расходной характеристики; эмпирическая функциональная зависимость α=ϕ(νH) имеет вид гиперболической функции α=aH+b (где а и b - коэффициенты аппроксимации), график которой строится по двум точкам в программе «Advanced Grapher» - [νH1-1,05] и [νH2-2], где νH1 и νH2 - скорости истечения, соответствующие верхней границе условной переходной и нижней границе условной ламинарной зон экспериментальной расходной характеристики; для вычисления значений скоростных напоров в зоне условного турбулентного режима истечения для экспериментальной расходной характеристики, а также значений скоростных напоров - в зоне фактического турбулентного режима истечения для прогнозируемой расходной характеристики натурной жидкости принимается значение α=1,05.

На чертеже представлены: 1 - дозирующее устройство (ДУ); 2 - устройство управления (УУ); 3 - расходная емкость с водопроводной водой; 4 - расходный резервуар (РР) из нержавеющей стали полезным объемом 30 л (диаметр резервуара - 20 см, высота - 100 см) с герметичной крышкой; 5 - уровнемерная трубка; 6 - барботажная трубка; 7 - питающая трубка; 8 - трубопровод линии пополнения РР; 9 - запорный вентиль на линии пополнения РР; 10 - мерная емкость (ME); 11 - запорное устройство - пневмоклапан (шаровой кран с двухсторонним пневмоприводом); 12 - сливной наконечник; 13 - гибкий трубопровод на линии налива воды в ME; 14 - струйный эжектор; 15 - 17 - запорные вентили; 18 - регулятор давления П; 19 - пневмотумблер «Контроль ПЗ/Контроль П»; 20 - пневмотумблер «Давление в РР»; 21 - задатчик давления ПЗ; 22, 23 - пневмореле; 24, 25 - регулируемые дроссели; 26 - пятимембранный элемент сравнения; 27 - блок усилителей мощности; 28 - сборный коллектор; 29 - образцовый манометр; 30 - пневмоэлектронное временное устройство; 31 - электронный таймер; 32 - электропневмопреобразователь; 33 - логический элемент «НЕ».

Испытуемый гидравлический тракт состоит из стальной питающей трубки (ПТ) 7, гибкого трубопровода 13 из ПВХ, запорного устройства 11 - двухпозиционного пневмоклапана в виде шарового крана с двухсторонним пневмоприводом, снабженного сливным наконечником (СН) 12.

РХ САД представляет собой функциональную зависимость объемного расхода Q на выходе СН 12 от величины располагаемо напора где р - избыточное давление в РР 4, γB=1 г/см3 - удельный вес воды, Н1 - текущее положение уровня жидкости в РР, H2=66 см - высота установки СН относительно нижней торцевой кромки ПТ 7 (высота истечения), р+γBH1=П - полное давление в РР, измеряемое БТ 6.

В состав УУ 2 входят пневматический пропорциональный регулятор давления П (РД) 18 и пневмоэлектронное временное устройство (ВУ) 30.

РД 18 содержит пневмотумблер 19 «КонтрольПЗ/КонтрольП» и пневмореле 20 для контроля по манометру 29 заданного (ПЗ) (при выключенном положении тумблера 19) или фактического (П) (при включенном положении тумблера 19) значений полного давления в РР, пневмотумблер 20 «Давление в РР» для включения подачи сжатого воздуха в РР, задатчик 23 давления ПЗ, пневмореле 22 и 23, регулируемые дроссели 24 и 25, пятимембранный элемент сравнения 26 и блок усилителей мощности 27. Пневмотумблеры 19 и 20 и задатчик 21 функционально образуют пульт управления. Функции собственно регулятора давления выполняет пятимембранный элемент сравнения 26.

РД 18 обеспечивает автоматическое поддержание заданной постоянной величины давления в каждом j-ом цикле проливки (дозирования) путем подкачки сжатого воздуха в РР от усилителей мощности 27 через сборный коллектор 28 и запорный вентиль 16, обеспечивая постоянство величины располагаемого напора а следовательно, и расхода Qj при открытом положении пневмоклапана 11.

ВУ 30 содержит электронный таймер 31 с цифровой индикацией времени дозирования, устанавливаемого посредством набора кнопками «SET», и электропневмо-преобразователь (ЭПП) 32 и логический элемент «НЕ» 33. Запуск операции дозирования производится по команде от сенсорной кнопки «СТАРТ», при касании которой вырабатывается управляющий сигнал (24 в) на переключение пневмоконтактов ЭПП 32 и перевод пневмоклапана 11 в открытое положение на время дозирования воды в ME 10.

Для выявления закономерностей зависимости параметра ξ от величин конструктивных параметров гидравлического тракта ДУ измерения расхода воды проводились с использованием в качестве трубопровода 13 гладкостенных гибких трубок из ПВХ трех типоразмеров внутреннего диаметра: d1=0,6; 0,8 и 1,0 см и длины L1=328; 433 и 277 см. Типоразмеры питающей трубки 7 (d2=1,5 см и L2=112 см) во всех циклах измерений оставались неизменными. В качестве СН 12 использовались металлические трубки общей длиной 20 мм с зенковкой входной кромки внутренним диаметром выходного тонкостенного канала dH=0,2; 0,3; 0,4; 0,5 и 0,6 см.

Рабочий цикл измерений представлял собой последовательные замеры расхода воды в гидравлическом тракте ДУ 1 для каждой (из 3-х) выбранной пары типоразмеров [d1_L1] гибкого трубопровода и сливного наконечника (dH) весовым методом. При этом каждый из 15 (3×5) рабочих циклов измерений заключался в последовательных проливках тракта водопроводной водой при различных постоянных значениях располагаемого напора где ПЗj - заданная величина полного давления Пj (j=1, 2…31), H2=66 см - высота истечения; при этом H1=Hmax=184 см, H31=Hmin=4 см. Значения давлений устанавливали задатчиком давления 21 по образцовому манометру 29 (пределы измерения давления - 0-1 кг/см2; 1 дел. шкалы - 2 см вод. ст.; класс точности - 0,16) с дискретностью 10 см вод. ст. при вод. ст. и далее (при ПЗj=100-70 см вод. ст.) - с дискретностью 2 см вод. ст.

Для оперативной обработки, представления и анализа полученных в ходе измерений массивов экспериментальных и расчетных данных использовались компьютерные программы «Microsoft Excel» и ((Advanced Grapher».

В результате проведения операций проливки для каждой пары [d1_L1]_dH типоразмеров конструктивных параметров гидравлического тракта ДУ и для каждого значения располагаемого напора Hj, в соответствии с заявляемым способом, последовательно вычислялись: соответствующий напору Н, объемный расход воды где - объем дозы, Gj[г] - вес дозы, регистрируемый с помощью взвешивания на электронных весах, γВ=1 г/см3 - удельный вес воды, tj[c] - время дозирования; скорость истечения воды из СН - и соответствующее этой скорости число Рейнольдса скорости течения в гибкой - и питающей - трубках; соответствующие этим скоростям числа Рейнольдса и коэффициенты трения и и потери напора на трение в гибкой - и в питающей - трубках; суммарные потери напора на трение - скоростной напор на выходе СН - (коэффициент Кориолиса α первоначально принят равным 1,05 как при турбулентном режиме истечения); суммарные потери на местных сопротивлениях - суммарный коэффициент потерь напора на местных сопротивлениях

На последующем этапе обработки полученных экспериментальных и расчетных данных с помощью программы «Advanced Grapher» полученные в табличном виде экспериментальные РХ тракта, Qj=F(Hj), по процедуре, подробно описанной в патенте-прототипе (и, кратко, - в п. 1 Формулы заявляемого изобретения) делили на зоны турбулентного, условно переходного и ламинарного режимов истечения.

Далее в качестве иллюстрации методики и вычислительных процедур, используемых для осуществления заявляемого способа, а также для его экспериментальной проверки (вычисления и анализа погрешностей прогнозирования РХ) циклы проливки гидравлического тракта ДУ (фиг.1) и измерения расхода воды проводились с использованием указанных выше трех типоразмеров внутреннего диаметра (d1) гибкого трубопровода и пяти типоразмеров (dH) внутреннего диаметра СН при различной условно «высокой» (УВ) и условно «низкой» (УН) сезонной температуре воды, т.е. при ее различной вязкости. Наибольшая разность температуры воды, при этом, составила ≈14°С. В результате этого этапа было получено 15 (3×5) массивов экспериментальных и расчетных данных параметров течения с УВ и столько же массивов данных с УН температурой воды.

Методика экспериментальной проверки предлагаемого способа заключалась в определении расчетных РХ, которые получились бы при УН температуре воды, по экспериментальным и расчетным данным, полученным при ее УВ температуре. О достоверности получаемых результатов прогнозирования РХ можно было судить по результатам вычисления и анализа относительных погрешностей расчетных значений расходов QPk (где k=j=1,2,…31) при их сравнении с соответствующими экспериментальными значениями QЭk, полученными в циклах проливки гидравлического тракта ДУ водой при ее УН температуре.

Вычисление погрешностей прогнозирования РХ САД воды при ее УН температуре проводилось в двух вариантах выбора значений коэффициента Кориолиса α - для α=1,05=const (как при турбулентном режиме истечения) и для переменных значений α, устанавливаемых в соответствии с представленным в п. 2. Формулы изобретения алгоритмом определения эмпирической зависимости α=a/νH+b, где а и b - коэффициенты аппроксимации.

Для организации компьютерных вычислительных процедур в программе ((Microsoft Excefo для α=1,05 были сформированы три пары папок (по числу типоразмеров гибкого трубопровода 13 на фиг.1), каждая из которых содержала парные листы «d1+- d1-» с полученными на измерительной установке данными параметров течения с УВ и с УН температурой воды, соответственно. Аналогично, для α=a/νH+b были сформированы три пары папок с листами «d1+A - d1-A», содержащими скорректированные (в соответствии с рассматриваемой далее процедурой) расчетные данные параметров течения.

Таким образом, процедура прогнозирования РХ САД воды при УН температуре по экспериментальным и расчетным данным, полученным при УВ температуре, в соответствии с уравнением Бернулли - описывающим течение в гидравлическом тракте ДУ измерительной установки, и в соответствии с п. 1 Формулы изобретения имела следующую последовательность операций.

1. Деления на листе «d1+А» полученной в цикле «проливок» табличной функции Qj=F(Hj) на условные зоны режимов истечения и корректировка значений потерь напора на выходе СН для условно переходной и ламинарной зон - и, соответственно, табличных значений функции согласно предложенному алгоритму определения α.

2. Проведение регрессионного анализа скорректированных значений табличных функций построение соответствующих выделенным зонам графиков аппроксимирующих функций и определение их уравнений.

3 Определение точек пересечения графиков: Re1 и Re2 фиксирующих, соответственно, верхнюю и нижнюю границы фактической переходной зоны для прогнозируемой РХ воды при УН температуре.

4. Выделение на листе «d1-А», содержащем ряд экспериментальных значений расхода воды при УН температуре - и - соответствующий ряд скоростей истечения прогнозируемой РХ, фактических зон режимов истечения путем образования пустых разделительных строк. При этом экспериментальные и расчетные данные параметров течения для чисел относились к фактической зоне турбулентного режима истечения, для чисел - к фактической зоне ламинарного режима, а для чисел - к фактической переходной зоне прогнозируемой РХ. Если графики функций и не пересекались или точка Re1>Rek, то прогнозируемая РХ не имела зону турбулентного режима истечения.

5. Корректировка на листе «d1-А» значений потерь напора на выходе СН для фактических переходной и ламинарной зон с учетом уравнения α=a/νH+b для коэффициента Кориолиса.

6. Формирование третьего дополнительного листа в первый столбец которого заносились числа Rek с листа «d1-А», а во второй столбец - вычисленные по полученным на 2-й операции уравнениям расчетные значения суммарных коэффициентов потерь напора на местных сопротивлениях - для прогнозируемой РХ.

7. Формирование четвертого дополнительного листа «ПРОГНОЗ», в первый, второй и третий столбцы которого заносились вычисленные на листе «d1-A» значения суммарных потерь на трение - скорректированные (в зависимости от а) значения скоростных напоров на выходе СН - и экспериментальные значения скоростей истечения νHk, соответственно. В четвертый и пятый столбцы листа «ПРОГНОЗ» копировались парные значения с листа полученные в предыдущей операции, в шестой столбец - вычисленные по формуле значения суммарных потерь напора на местных сопротивлениях и в седьмой столбец - вычисленные по формуле расчетные значения располагаемого напора. Для удобства построения точек прогнозируемой РХ в последний 8-й столбец копировались экспериментальные значения QЭk с листа «d1-А».

8. По полученным для каждой из выделенных зон табличным функциям путем их регрессионного анализа в программе «Advanced Grapher» и аппроксимации степенными функциями были получены уравнения расчетной РХ воды при УН температуре для фактических зон турбулентного (s=1), переходного (s=2) и ламинарного (s=3) режимов истечения.

Следующий этап проверки предложенного способа определения РХ касался процедуры вычисления относительных погрешностей расчета расхода прогнозируемой РХ САД воды при УН температуре. Отметим, что для всех 15-ти массивов экспериментальных и расчетных данных, полученных для УН температуры воды, и обоих вариантах вычисления коэффициента Кориолиса α (а=1,05 и α=a/νH+b) эти погрешности имеют максимальные значения в верхней и нижней точках диапазона изменения располагаемого напора, т.е. при Н=Hmax=184 см и Н=Hmin=4 см. Вычисление погрешностей проводилось по формулам:

где и - экспериментальные значения расхода при H=184 см и H=4 см, соответственно, а и - расчетные значения расхода при тех же значениях H.

Параметры и определялись посредством использования операции «трассировки» в программе «Advanced Grapher» графиков аппроксимирующих степенных функций, полученных в результате регрессионного анализа табличных функций и на листе «d1-А», соответствующих фактическим турбулентной и ламинарной зонам истечения. Соответственно, параметры и определялись посредством использования операции «трассировки» графиков аппроксимирующих степенных функций, полученных в результате регрессионного анализа табличных функций и на листе «ПРОГНОЗ».

Ниже приведена сводная таблица расчетных значений относительных погрешностей расхода для верхней и нижней границ диапазона изменения располагаемого напора.

Как следует из приведенной таблицы, максимальные значения погрешностей имеют место на нижней границе диапазона измерений (при Н=4 см). При вычислении α по эмпирической формуле - α=а/νH+b эти погрешности уменьшаются, что свидетельствует о правомерности предложенного в п. 2 Формулы изобретения алгоритма учета влияния на прогнозируемую РХ коэффициента Кориолиса.

Полученные значения погрешностей прогнозирования РХ, максимальная из которых получена для трубки с внутренним диаметром d1=0,8 см и СН с внутренним диаметром dH=0,5 см, имеют достаточно малую величину и вполне приемлемы для практического применения предложенного способа для расчета РХ проектируемой системы дозирования натурной жидкости.

Таким образом, предлагаемый способ позволяет в лабораторных условиях оперативно и с высокой степенью достоверности прогнозировать расходные характеристики проектируемого тракта с варьируемыми физическими параметрами натурной жидкости (удельным весом γ и кинематической вязкостью ν), диапазонами изменения расхода и различной длиной образующих тракт участков трубопроводов.

Похожие патенты RU2765801C1

название год авторы номер документа
СПОСОБ ОПРЕДЕЛЕНИЯ РАСХОДНОЙ ХАРАКТЕРИСТИКИ ГИДРАВЛИЧЕСКОГО ТРАКТА ДЛЯ ОБЛАСТИ ПЕРЕХОДА ОТ ТУРБУЛЕНТНОГО К ЛАМИНАРНОМУ РЕЖИМУ ИСТЕЧЕНИЯ 2019
  • Безменов Василий Серафимович
RU2709034C1
СПОСОБ ОПРЕДЕЛЕНИЯ РАСХОДНОЙ ХАРАКТЕРИСТИКИ ГИДРАВЛИЧЕСКОГО ТРАКТА ПРИ ТУРБУЛЕНТНОМ РЕЖИМЕ ИСТЕЧЕНИЯ 2018
  • Безменов Василий Серафимович
RU2680987C1
СПОСОБ ОПРЕДЕЛЕНИЯ РАСХОДНОЙ ХАРАКТЕРИСТИКИ ГИДРАВЛИЧЕСКОГО ТРАКТА И УСТРОЙСТВО ДЛЯ ЕГО ОСУЩЕСТВЛЕНИЯ 2015
  • Безменов Василий Сирафимович
RU2582486C1
Способ управления электроприводом башенной насосной установки 1987
  • Кудрявцев Владимир Иванович
SU1472881A1
УЛУЧШАЮЩАЯ ЭКСПЛУАТАЦИОННЫЕ ЭНЕРГЕТИЧЕСКИЕ ХАРАКТЕРИСТИКИ МАШИН ЖИДКАЯ ПРИСАДКА, УГЛЕВОДОРОДНАЯ ЖИДКОСТЬ НА ОСНОВЕ НЕФТЕПРОДУКТОВ, ИСПОЛЬЗУЕМАЯ В МАШИНАХ, И ЖИДКОЕ УГЛЕВОДОРОДНОЕ ГОРЮЧЕЕ 2006
  • Гапонов Валерий Дмитриевич
  • Стернин Леонид Евгеньевич
  • Фатуев Игорь Юрьевич
  • Чванов Владимир Константинович
RU2343187C2
СПОСОБ ОПРЕДЕЛЕНИЯ КОЭФФИЦИЕНТА ГИДРАВЛИЧЕСКОГО СОПРОТИВЛЕНИЯ ТРАКТА И УСТРОЙСТВО ДЛЯ ЕГО ОСУЩЕСТВЛЕНИЯ 2003
  • Пинке И.М.
RU2240525C1
СПОСОБ ОПРЕДЕЛЕНИЯ ХАРАКТЕРИСТИК БЕЗРАЗМЕРНЫХ ПАРАМЕТРОВ ТЕЧЕНИЯ ПОТОКОВ РАБОЧЕЙ ЖИДКОСТИ В ДРОССЕЛЬНЫХ ОКНАХ ЗОЛОТНИКОВЫХ ГИДРОРАСПРЕДЕЛИТЕЛЕЙ 2004
  • Белоногов Олег Борисович
RU2282065C2
Способ определения расходных характеристик (РХ) струйных датчиков расхода (СДР) 2021
  • Вологодский Николай Витальевич
  • Канунников Юрий Александрович
  • Сторожев Валерий Васильевич
  • Цацорина Лидия Усмановна
RU2770512C1
Способ формирования управляющего сигнала по углу крена модели гиперзвукового летательного аппарата (ГЛА) для контроля аэродинамической идентичности по числам Рейнольдса траекторий полёта модели и натурного изделия при проведении опережающих лётных исследований аэродинамических характеристик 2017
  • Ловицкий Лаврентий Лаврентьевич
  • Садртдинов Владислав Диясович
  • Бадретдинова Айгуль Булатовна
RU2650331C1
СПОСОБ ПРОЕКТИРОВАНИЯ ЦЕНТРОБЕЖНЫХ НАСОСОВ 2013
  • Кушнарев Владимир Иванович
  • Кушнарев Иван Владимирович
  • Обозный Юрий Сергеевич
RU2542160C1

Иллюстрации к изобретению RU 2 765 801 C1

Реферат патента 2022 года СПОСОБ ОПРЕДЕЛЕНИЯ РАСХОДНОЙ ХАРАКТЕРИСТИКИ ГИДРАВЛИЧЕСКОГО ТРАКТА В ПЕРЕХОДНОЙ ОБЛАСТИ ИСТЕЧЕНИЯ

Изобретение относится к области измерительной технике и может быть использовано для расчета пропускной способности проектируемых гидравлических трактов транспортных и дозирующих систем в химической, нефтехимической, авиационной, текстильной, лакокрасочной и других отраслях промышленности, в частности, - для определения (прогнозирования) расходных характеристик систем автоматизированного порционного дозирования натурных жидкостей с поддержанием постоянного перепада давлений на линии налива, с отмериванием дозы по времени. Предлагаемый способ характеризуется тем, что модель тракта с проектными конструктивами местных сопротивлений и участков соединительных трубопроводов длиной с внутренним диаметром di с концевой раздачей через запорное устройство с цилиндрическим сливным наконечником с внутренним диаметром dH проливают водопроводной водой посредством ее многократного порционного дозирования в мерную тару при различных, поддерживаемых на заданном постоянном уровне, величинах действующего напора Hj (где j = 1, 2 …, n=31 - число проливок в каждом цикле проливок), ступенчато уменьшаемых с шагом 10 см от максимального Hjmax=Hj=1=184 см до промежуточного значения Hj=Hj=18=34 см - и далее, с шагом 2 см, - до нижней границы диапазона измерений - Hjmin=H31=4 см. При этом с помощью электронных таблиц программы «Microsoft Excel» последовательно вычисляют: объемные расходы Qj [см3/c]=Vj/tj, соответствующие напорам Hj, где объем дозы Vj регистрируют с помощью взвешивания на электронных весах, а время дозирования tj задают и отмеривают с помощью электронного таймера; скорости истечения νHj[см/с]=QjH из сливного наконечника (где ƒн=πdH2/4), потери напора - hHj[см]=α⋅νHj2/2g и числа Рейнольдса потока - RejHjdHB на его выходе, где коэффициент Кориолиса α первоначально принимают равным 1,05 как при турбулентном режиме истечения, а кинематическую вязкость воды νB рассчитывают по ее температуре в данном цикле проливок; скорости течения - νij[см/с]=νHj(dH/di)2, числа Рейнольдса - ReijijdiB, коэффициенты трения - λij=64/Reij, и потери напора на трение - на отдельных участках трубопровода; суммарные потери напора на трение - hTj[смМ]=Σi[hij], на местных сопротивлениях в тракте - hMj[см]=Hj-hHj-hTj и суммарные коэффициенты потерь напора на местных сопротивлениях - ξj=2ghMjHj2. Полученную в табличном виде экспериментальную расходную характеристику тракта Qj=F(Hj) делят на зоны турбулентного, условно переходного и условно ламинарного режимов истечения посредством последовательного сравнительного анализа числовых значений стандартного отклонения ее точек и выбора количеств точек Qj=1,2,…p (р≥6 - в верхней области диапазона измерений) и Qj=34,33,…q (q≥6 - в нижней области диапазона) с минимальными значениями стандартного отклонения от графика соответствующей этим точкам степенной аппроксимирующей функции. Технический результат заключается в повышении точности и достоверности определения расчетной расходной характеристики гидравлического тракта в переходной области при транспортировании по нему жидких продуктов различной вязкости и возможности ее прогнозирования при варьировании заданного диапазона изменения расхода и длины образующих тракт участков трубопроводов. 1 з.п. ф-лы, 1 ил., 1 табл.

Формула изобретения RU 2 765 801 C1

1. Способ определения расходной характеристики гидравлического тракта в переходной области истечения, характеризующийся тем, что модель тракта с проектными конструктивами местных сопротивлений и участков соединительных трубопроводов длиной с внутренним диаметром di с концевой раздачей через запорное устройство с цилиндрическим сливным наконечником с внутренним диаметром dH проливают водопроводной водой посредством ее многократного порционного дозирования в мерную тару при различных, поддерживаемых на заданном постоянном уровне, величинах действующего напора Hj, где j=1,2…, n=31 - число проливок в каждом цикле проливок, ступенчато уменьшаемых с шагом 10 см от максимального см до промежуточного значения см - и далее, с шагом 2 см, - до нижней границы диапазона измерений - см; с помощью электронных таблиц программы «Microsoft Excel» последовательно вычисляют: объемные расходы соответствующие напорам Hj, где объем дозы Vj регистрируют с помощью взвешивания на электронных весах, а время дозирования tj задают и отмеривают с помощью электронного таймера; скорости истечения из сливного наконечника, где ƒH=πdH2/4, потери напора - и числа Рейнольдса потока - на его выходе, где коэффициент Кориолиса α первоначально принимают равным 1,05 как при турбулентном режиме истечения, а кинематическую вязкость воды νB рассчитывают по ее температуре в данном цикле проливок; скорости течения - числа Рейнольдса - коэффициенты трения - и потери напора на трение - на отдельных участках трубопровода; суммарные потери напора на трение - на местных сопротивлениях в тракте - и суммарные коэффициенты потерь напора на местных сопротивлениях - полученную в табличном виде экспериментальную расходную характеристику тракта Qj=F(Hj) делят на зоны турбулентного, условно переходного и условно ламинарного режимов истечения посредством последовательного сравнительного анализа числовых значений стандартного отклонения ее точек и выбора количеств точек , р ≥ 6 - в верхней области диапазона измерений и , q ≥ 6 - в нижней области диапазона с минимальными значениями стандартного отклонения от графика соответствующей этим точкам степенной аппроксимирующей функции; при этом полученные для вышеуказанных зон значения табличных функций , где индексы s=1, 2, 3 относятся к условно турбулентной, условно переходной и условно ламинарной зонам, соответственно, корректируют с учетом влияния коэффициента Кориолиса α, который увеличивают от значения α=1,05 при турбулентном режиме истечения до α=2 при ламинарном режиме в соответствии с принятой эмпирической функциональной зависимостью α=ϕ(νH) для вычисления скорректированных значений скоростных напоров на выходе сливного наконечника - и соответствующих скорректированных значений табличных функций для каждой из выделенных зон строят графики функций и путем их регрессионного анализа получают уравнения аппроксимирующих функций - определяют точки пересечения графиков: и фиксирующих, соответственно, верхнюю и нижнюю границы фактической переходной зоны, для прогнозируемой расходной характеристики; формируют ряд произвольно, с небольшим шагом уменьшаемых значений расхода и - соответствующий ряд скоростей истечения прогнозируемой расходной характеристики натурной жидкости; для каждого значения k и для каждой из выделенных фактических зон режимов истечения натурной жидкости вычисляют соответствующие этим зонам числа Рейнольдса на выходе сливного наконечника - , где ν - кинематическая вязкость натурной жидкости, и по полученным уравнениям аппроксимирующих функций вычисляют суммарные коэффициенты потерь напора на местных сопротивлениях - для прогнозируемой расходной характеристики натурной жидкости вычисляют: скорости течения - числа Рейнольдса - коэффициенты трения - потери напора на трение на отдельных участках трубопровода - и суммарные потери напора на трение - потери напора на выходе сливного наконечника - и на местных сопротивлениях в тракте - расчетные значения располагаемого напора - по полученным в табличном виде для каждой из выделенных фактических зон режимов истечения функциям Qk=F(HPk) путем их аппроксимации степенными функциями получают соответствующие графики и уравнения зависимости расхода от располагаемого напора, т.е. определяют для каждой из трех фактических зон режимов истечения искомые расходные характеристики натурной жидкости.

2. Способ по п. 1, отличающийся тем, что увеличение коэффициента α от 1,05 до 2 начинается для полученной в цикле проливок экспериментальной расходной характеристики модели гидравлического тракта от верхней границы условно переходной зоны истечения, а для прогнозируемой расходной характеристики - от верхней границы фактической переходной зоны; увеличение коэффициента α до значения α=2 заканчивается на нижней границе ламинарной зоны экспериментальной расходной характеристики; эмпирическая функциональная зависимость α=ϕ(νH) имеет вид гиперболической функции α=a/νH+b, где а и b - коэффициенты аппроксимации, график которой строится по двум точкам в программе «Advanced Grapher» - [νH1-1,05] и [νH2-2], где νH1 и νH2 - скорости истечения, соответствующие верхней границе условной переходной и нижней границе условной ламинарной зон экспериментальной расходной характеристики; для вычисления значений скоростных напоров в зоне условного турбулентного режима истечения для экспериментальной расходной характеристики, а также значений скоростных напоров - в зоне фактического турбулентного режима истечения для прогнозируемой расходной характеристики натурной жидкости принимается значение α=1,05.

Документы, цитированные в отчете о поиске Патент 2022 года RU2765801C1

СПОСОБ ОПРЕДЕЛЕНИЯ РАСХОДНОЙ ХАРАКТЕРИСТИКИ ГИДРАВЛИЧЕСКОГО ТРАКТА ДЛЯ ОБЛАСТИ ПЕРЕХОДА ОТ ТУРБУЛЕНТНОГО К ЛАМИНАРНОМУ РЕЖИМУ ИСТЕЧЕНИЯ 2019
  • Безменов Василий Серафимович
RU2709034C1
СПОСОБ ОПРЕДЕЛЕНИЯ РАСХОДНОЙ ХАРАКТЕРИСТИКИ ГИДРАВЛИЧЕСКОГО ТРАКТА ПРИ ТУРБУЛЕНТНОМ РЕЖИМЕ ИСТЕЧЕНИЯ 2018
  • Безменов Василий Серафимович
RU2680987C1
СПОСОБ ОПРЕДЕЛЕНИЯ РАСХОДНОЙ ХАРАКТЕРИСТИКИ ГИДРАВЛИЧЕСКОГО ТРАКТА И УСТРОЙСТВО ДЛЯ ЕГО ОСУЩЕСТВЛЕНИЯ 2015
  • Безменов Василий Сирафимович
RU2582486C1
US 7636640 B2, 22.12.2009.

RU 2 765 801 C1

Авторы

Безменов Василий Серафимович

Даты

2022-02-03Публикация

2020-09-08Подача