Сравнение методов расчета неустановившегося движения жидкости в круговых цилиндрических оболочках Текст научной статьи по специальности «Механика и машиностроение»

Механика жидкости и газа

СРАВНЕНИЕ МЕТОДОВ РАСЧЕТА НЕУСТАНОВИВШЕГОСЯ ДВИЖЕНИЯ ЖИДКОСТИ В КРУГОВЫХ ЦИЛИНДРИЧЕСКИХ ОБОЛОЧКАХ

Ф.В. РЕКАЧ, канд. техн. наук, доцент Российский университет дружбы народов, 117198, Москва, ГСП - 6, ул. Миклухо-Маклая, 6

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

КЛЮЧЕВЫЕ СЛОВА: стабилизатор давления, неустановившееся движение жидкости, цилиндрическая оболочка.

В качестве одного из эффективных средств гашения волновых процессов в трубопроводных системах рассматриваются стабилизаторы давления (СД). Принцип их работы основан на распределенном по длине трубопровода диссипативном и упругодемпфирующем воздействии на пульсирующий поток перекачиваемой среды. Наибольший эффект гашения достигается при диссипации энергии пульсаций на перфорационных отверстиях, равномерно распределенных по длине стабилизатора, а также вследствие демпфирования, обусловленного податливостью упругих элементов стабилизатора. Отличительная особенность стабилизаторов давления заключается в том, что они не нарушают форму трубопровода и имеют минимальные гидравлические сопротивления.

Стабилизатор работает следующим образом. При распространении волны повышенного давления (рис. 1) происходит перетекание транспортируемой среды через отверстия 2 из трубопровода 1 в демпфирующую камеру 3, упругий элемент при этом сжимается. При понижении давления (рис. 2) упругий элемент увеличивается в объеме, заполняя свободное пространство демпфирующей камеры, часть заполняющей ее среды вновь перетекает в трубопровод 1, обеспечивая сглаживание провала давления. Гашение пульсаций осуществляется также за счет дросселирования среды через перфорационные отверстия 2. При большой суммарной податливости упругих камер и оптимальном варианте суммарной площади перфорации можно добиться максимального уменьшения

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

Рис. 1 Рис. 2

Среди точных методов наиболее рациональными с точки зрения простоты, точности и легкости задания исходных данных являются методы ДЛАламбера и Лапласа.

Метод Д'Аламбера. Неустановившееся движение несжимаемой жидкости в трубопроводе без учета трения (после преобразований) описывается уравнениями движения и неразрывности следующего вида [1]:

д2р(х, г) _ д2p(x, г)

ах2 " С2 дг2 ' д 2Q (х, г) _ ^ддОхА

дх2 " С2 дг2 ,

где р(х, г) - давление в трубопроводе [н/м2], () - расход жидкости в магистрали [м3/с], F- площадь поперечного сечения магистрали [м2], с - скорость распространения волн давления [м/с], р- плотность жидкости [кг/м3].

Идея метода Д'Аламбера состоит в представлении неустановившегося движения жидкости на -ом участке в виде прямой и обратной волн, при этом общее решение системы (1) для участка х1 _ //-1 ^ I 1 записывается в виде:

Q (х, t )=--• р

Pt (х, t) = -c •

f

f

t - Xl

+g i

f \ t +

V

с

i J

( \ ( \

хг хг

t —- - gi t + —

с с

V у V у

(2)

Затем учитываются граничные условия и производится «сшивка» участков . Этим методом решались задачи на вынужденные колебания со стабилизатором давления [4].

Метод Лапласа. Система уравнений, описывающая движение жидкости в трубопроводе со стабилизатором, считая жидкость несжимаемой средой р = const, имеет вид [3]:

dpi(Хt)

дх

dPi(Xt) dt

= -Р-

dvI (х, t)

dt

+

2avI (х, t)

--рс^ dvI (х,t) дх

(3)

с

dPii fc t ) _ р

-_ - Р

dx

dvii fc t).

dt

2avII (x, t )

cPii (xt) _ Р 2 dvn (x, t) -_ ~рст -

dt

dx

dPn

dx

_~Р-

dvi,

dt

4

dv

"T-v„ep +_ 0, (8)

аш dx

dp

ПК

fnK

dv

-v„„„ --

fK _ 0, (10)

dx

ÊPa. _-pc2 .^YA.

dx dPK

_-Р-

dv

ПК

dt

dt

dx

dx (12)

_-Р'~

dv

dt

PlII - PnK _ A Vnep + B Vnep ' \Vnep

(5)

(6) (7) (9)

(11)

(13),

магистра, cK - ско-

где pI, vI, pII, vII, pIII, vIII - давление и скорость на участках ли I, II, III; a - коэффициент трения в магистрали трубопровода; с. рость звука в трубе, камере с жидкостью; vnep - скорость перетекания потока из трубопровода в предкамеру; pnK, vnK, pK, vK - давление и ско-

^ i 2 2 \

рость в предкамере и камере стабилизатора; fnK _ — (d ПК - dUI ) - площадь

поперечного сечения предкамеры кольцевой формы; dIII, dnK - средний диаметр трубы магистрали и внешний диаметр кольца предкамеры; A и B -постоянные коэффициенты, получаемые экспериментально.

Уравнения (3) - (13) представляют собой систему относительно одиннадцати неизвестных функций Pi , vI, Pii , vu, рш, Уш, Рпк, vnK, Рк , vK, vnep, в

которой (3), (5), (7), (9) и (11) являются уравнениями равновесия элемента потока; (4), (6), (8), (10) и (12) - уравнениями неразрывности и уравнение (13) выражает условие перетекания жидкости через перфорацию.

Длина волны, которую гасит стабилизатор, как правило, значительно превышает его длину 1СТ . Поэтому принимается гипотеза о том, что давление в различных элементах стабилизатора не меняется по координате x.

Система (3) - (13) сводится к четырем уравнениям, к которым вместе с граничными условиями (в граничных сечениях задается известная величина давления или расхода) применяется преобразование Лапласа по времени. Преобразованная система дифференциальных уравнений уже не является системой в частных производных и имеет точное решение. Переходя обратно во временную область с помощью формулы обращения Меллина, методом Д'Аламбера находятся оригиналы функций pI(5,x), vI(5,x), pII(5,x), vII(5,x).

Сравнение методов.

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

2) Методом Лапласа легче решаются задачи на гидроудар (в граничном сечении с постоянным расходом мгновенно закрывается задвижка) [2]. Также исследовались задачи с нелинейным законом перетекания жидкости через перфорацию, что в методе Д'Аламбера сделать трудно (если вообще возможно). Задачи с наличием трения в магистральном трубопроводе можно решать обоими методами.

3) При решении задач на вынужденные колебания метод Д'Аламбера более удобен, так как легче задавать гармоники вынужденных колебаний.

4) Часть результатов решения по рассматриваемым методам сравнивались и разница амплитуд давления не превосходила 10%.

Выводы:

1) При решении методом Д'Аламбера значительно легче рассчитывать гидравлические системы, включающие несколько разных конструктивных участков, однако труднее учитывать нелинейности задач.

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

3) Решение задач точными методами дает правдоподобную физическую интерпретацию, что может быть использовано при поверке решений численными методами.

Л и т е р а т у р а

1. Чарный И.А. Неустановившееся движение реальной жидкости в трубах. 2 -е изд // М.: Недра, 1975 - 296 с.

2. Ганиев Р.Ф., Низамов Х.Н., Дербуков Е.И. Волновая стабилизация и предупреждение аварий в трубопроводах. - М.: Изд-во МГТУ им. Н.Э. Баумана, 199 с.

3. Рекач Ф.В. Расчет стабилизатора давления диссипативного принципа действия методом Лапласа с учетом трения в магистральном трубопроводе и нелинейной зависимости скорости перетекания жидкости через перфорации// Строительная механика инженерных конструкций и сооружений. - 2007. - №4. - С. 56-60.

4. Рекач Ф.В. Исследование вынужденных колебаний в круговых цилиндрических оболочках методом Д'Аламбера со стабилизатором давления диссипативного типа// Строительная механика инженерных конструкций и сооружений. -2007.-№2.- С. 47-52.

5. Young W.R. and Wolfe C.L. Generation of surface waves by shear-flow instability// J. Fluid Mechanics. - 2014. - Vol. 739. - P. 276 - 307.

Reference

1. Œarny, I.A. (1975). Unsteady Motion of Real Liquid in Tubes, Moscow: Nedra, 296 p.

2. Ganiev, R.F, Nizamov, H.N., Derbukov, E.I. (1996). Wave Stabilization and Prevention of Accidents in Pipe-Lines, Moscow: Izd-vo MGTU, 260 p.

3. Rekach, F.V. (2007). Analysis of dissipative principle pressure stabilizer by Laplace method taking friction in main tube and nonlinear dependence of flouting liquid speed through perforation in consideration, Structural Mechanics of Engineering Constructions and Buildings, № 4, р. 56-60.

4. Rekach, F. V. (2007). Analysis of forced oscillations in circular cylindrical shells by D'Alamber method with dissipative principle pressure stabilizer, Structural Mechanics of Engineering Constructions and Buildings, № 2, р.47-52.

5. Young, W.R. and Wolfe, C.L. (2014) Generation of surface waves by shear-flow instability, J. Fluid Mechanics, Vol. 739, p. 276 - 307.

METHODS OF COMPARISON IN ANALYSIS OF UNSTEADY LIQUID MOTION IN CIRCULAR CYLINDRICAL SHELLS

F.V. Rekach

Peoples Friendship University of Russia, Moscow

A short essay of pressure pipelines analysis on water hammer and pressure oscillations by exact methods is described in this paper.

KEY WORDS: oscillations, pressure stabilizer, unsteady flow, cylindrical shell.