Электронный научный журнал "Математическое моделирование, компьютерный и натурный эксперимент в естественных науках" http://mathmod.esrae.ru/ URL статьи: mathmod.esrae.ru/15-50 Ссылка для цитирования этой статьи:
Агеев Р.В., Евдокимова Е.В., Ковалева И.А., Могилевич Л.И. Динамика осесимметричного течения вязкой несжимаемой жидкости в упругой трубе кругового и кольцевого сечений // Математическое моделирование, компьютерный и натурный эксперимент в естественных науках. 2017. №3.
Выполнено при поддержке гранта РФФИ № 16-01-00175a._
УДК 681.03.06:531.383:532.516
ДИНАМИКА ОСЕСИММЕТРИЧНОГО ТЕЧЕНИЯ ВЯЗКОЙ НЕСЖИМАЕМОЙ ЖИДКОСТИ В УПРУГОЙ ТРУБЕ КРУГОВОГО И
КОЛЬЦЕВОГО СЕЧЕНИЙ
Агеев Р.В.1, Евдокимова Е.В.2, Ковалева И.А.3, Могилевич Л.И.4
Саратовский государственный технический университет имени Гагарина Ю.А.,
Россия, Саратов, r255@mail.ru Саратовский государственный технический университет имени Гагарина Ю.А.,
Россия, Саратов, eev2106@mail.com Саратовский государственный технический университет имени Гагарина Ю.А.,
Россия, Саратов, irinakovaleva1406@gmail.com 4Саратовский государственный технический университет имени Гагарина Ю.А.,
Россия, Саратов, mogilevich@sgu.ru
THE DYNAMICS OF VISCOUS INCOMPRESSIBLE LIQUID AXISYMMETRIC FLOW IN ELASTIC PIPE OF CIRCULAR AND
ANNULAR CROSS-SECTIONS
Ageev R.V.1, Evdokimova E.V2, Kovaleva I.A.3, Mogilevich L.I.4
eev2106@mail.com 3Yuri Gagarin State Technical University of Saratov, Russia, Saratov,
irinakovaleva1406@gmail. com 4Yuri Gagarin State Technical University of Saratov, Russia, Saratov,
mogilevich@sgu.ru
Аннотация. Рассматриваются задачи о движении вязкой несжимаемой жидкости в упругой цилиндрической оболочке под воздействием некоторых факторов, как с круговым так и кольцевым сечением. Во всех этих задачах имеются источники движения систем «жидкость- упругое тело».Для выявления динамических характеристик систем необходимо решать связанные задачи динамики упругих оболочек и вязкой несжимаемой жидкости.
Упрощающими задачу факторами являются малость радиуса оболочки по сравнению с её длиной, как в теории смазки. В получающихся задачах гидроупругости вводится ещё предположение о малости прогибов оболочки по сравнению c радиусом оболочки или толщиной слоя жидкости, но не с толщиной оболочки. В этих условиях уравнения гидродинамики линеаризуются. Уравнения динамики упругих элементов остаются нелинейными, если таковыми были изначально. При отсутствии возможности найти точное решение, применяется метод итерации для решения уравнений гидродинамической теории смазки.
Abstract. The problems of viscous incompressible liquid movement in an elastic cylinder shell as of circular, as well as of annular cross-sections under the impact of some factors are considered. All these problems contain the systems "liquid - elastic body" movement sources. To reveal dynamic characteristic of these systems, it is necessary to solve the connected problems of an elastic shell and viscous incompressible liquid dynamics. The factor, making the problem simpler, is a small size of a shell radius in comparison with its length, as in lubrication theory. The problems under consideration contain this supposition about a small size of the shell sagging in comparison with a shell radius, or liquid layer thickness, but not with a shell thickness. Under these conditions hydrodynamics equations become linear. The elastic elements dynamics equations remain non-liner, if they were such initially. If it is not possible to find the exact solution, the iteration method for solving the lubrication theory hydrodynamics equations applied.
В случае линейной теории оболочек и в предположении гармонических по времени воздействиях на систему, решение задачи может определяться в виде гармонической зависимости всех параметров задачи по времени. Это возможно в силу влияния вязкости жидкости, которая приводит к демпфированию колебаний и быстрому затуханию по времени решения из-за начальных условий и выходу на стационарные колебания.
Эти условия позволяют найти точные решения линеаризованных уравнений механики жидкости, записанных в виде уравнений теории смазки, но с учётом локального члена инерции для любых значений колебательного (смазочного) числа Рейнольдса оставаясь в рамках ламинарного движения.
В указанных случаях нормальное напряжение жидкости (давление) на поверхности упругого тела значительно больше касательного напряжения и последним обычно пренебрегают[1-15].
Рассмотрим возможность получения приближённого значения решения уравнений динамики жидкости методом итерации, пренебрегая на первом шаге локальным членов инерции и учитывая его на втором шаге итерации. При этом можно доказать, что метод итерации сходится при условии, что колебательное число Рейнольдса меньше единицы [1]. Этот подход позволяет отказаться от требования гармонического закона по времени всех параметров жидкости и упругих элементов. Появляется возможность решать нелинейные уравнения динамики оболочек или при негармонических законах изменения по времени источников движения.
На рисунках 1,2 представлены системы «жидкость - упругое тело». На рис. 1 представлена упругая оболочка с жидкостью 3 внутри [5]. На рис 2 представлены две соосные оболочки с жидкостью между ними [6,7,14].
Для упругой оболочки (рис. 1) внутренний радиус оболочки 2 значительно меньше её длины, то есть R << I и вводится малый параметр щ = ^ << 1. Прогиб упругих элементов 1,2 значительно меньше толщины слоя
жидкости 3, то есть Ж «Я, Ж « 8п, и — = 0(Л), — = О(Л), порядка Я.
Рис. 1
Рис. 2
Рассмотрим ламинарное течение в упругой цилиндрической оболочке. Уравнения динамики жидкости в осесимметричном случае имеют вид [16]
V+V+ух V+1др = у
д г дг х дх р дг
д Х +1 дУг
дг г дг г
Хг.+дК
дК дКх дК + + V
—- +--— = V
дх р дх
гд 2К
- +—
дг г дг дх
V дК + Ук- +—х = 0 дг г дх
Здесь Кх, Кг - проекции вектора скорости жидкости на продольную ось Ох и нормаль к этой оси в цилиндрической системе координат (пг, i), х-продольная координата вдоль оси оболочки трубы, г-расстояние от оси оболочки трубы, с внутренним радиусом - R1, длиной I.
Граничные условия - условия прилипания жидкости на внутренней оболочке, условия ограниченности скорости на оси трубы и условия на концах трубы для давления
дЖ т. ди „
V =--, К = — при г = Я- Ж
г дл дt
дК дК
г -К = 0, г-^ = 0при г=0
Р = РоС )при X = 1/2
и х = = - I/
и первое условие для давления можно заменить условием симметрии давления
= 0 при х=0.
Здесь и - продольное упругое перемещение оболочки; Ж - прогиб, положительный к центру кривизны; R1 - внутренний радиус оболочки толщины ^ и длины 1.
Поверхностная нагрузка определяется напряжением со стороны жидкости
Чгг = Чп =- Р + 2Р
д¥г (д¥г д¥х
-г-, Чх = Чх =Р -г- + ^
дг I дх дг
Введем безразмерные переменные
С = , 1= —, Т = М , Ж = wmUз, и = ияи„
= Wm®Uz , Ух = Р = Ро +
Р™я®
<< 1, Re = Я2®/ V
Ш = /I << 1, Л = я
Учитывая наличие малых параметров у и X получим уравнения динамики жидкости (теория смазки с учетом локального члена инерции) и граничные условия
дР = о, Re = 1 , 1^)+^= 0,
д£ дт д, С дс{ дс) с дСУ а д,
диз тт иш ди „ .. „ и, = —, и, = яТ при £ = 1, £
ди, ди, —£ = С—1 дС дС
= о при £ = 0.
Р = 0 при , = -1, — = 0 при , = 0.
Полагая гармоническую зависимость от времени давления, скоростей жидкости и упругих перемещений оболочки находим
Р = 1-1 ^ ди
-а. Jо
= Яе(а-1)
V2 ^я дт2
•»о дт2
+16/ Г
d1 + 8/
дт ^ дт
При этом
Чгг = Чп =- Р =
д ШUmU1 3^ + 16/д7[ WmU
Чгх = Чх =р
-л2 -л
Re(a-1) -2- + 8/ — дт дт
здесь а-1 =
Re p0 - 2VReq0
Re p0 - 4VReq0 + 4s0 Re p0 - ^VReq0 + 4s0
Re p0 - ^VReq0 + 4s,
p0 = ber2VRe + bei2 VRe, q0 = berVRe + bei&VRe - ber& VRebeiVRe,
r0 = berVReber& VRe + beiVRebei& VRe, s0 = ber&2 VRe + bei&2 VRe,
ber м bei функции Кельвина нулевого порядка, штрих означает производную.
Отметим, что
а = >3, Y = 1 при Re ^ 0
а=1, а-1 = V2/Re , y = 2VRe/8, при Re > 20
Отказываясь от гармоничности по времени параметров течения, что означает отсутствие гармоничности p0 (t) и, считая а = у. , где
отсутствие гармоничности р0 ^)
t0 - характерное время, применим метод итерации к задаче (6).
На первом шаге итерации отбросим локальный член инерции и получим
i f 1 итщ dU1 ri dU3
Ui=1Z -1>6|
Urn¥ dU1
На втором шаге итерации будем иметь
р = £ j 2Re
J° дт
^ f + 8
Можно снова ввести частотозависимые коэффициенты а(ю) и у(ю) и считать, что Re - любое, при этом получаем
Р = £ j 2Re
= Re(а -1)
+16/
+ 8у
что формально совпадает с точным решением (7)для гармонических законов изменения параметров среды, но позволяет использовать эти формулы как для негармонического изменения со временем р0, так и для нелинейных уравнений динамики оболочки.
Следует отметить, что в [1] доказана сходимость метода итерации. Рассмотрим динамику тонкого слоя вязкой несжимаемой жидкости, заключённого между двумя соосными упругими оболочками (см. рис. 2).
Внешняя оболочка 1 имеет внутренний радиус R1, толщину ^((1) и радиус срединной поверхности R(1) = Л + 2. Внутренняя оболочка 2 имеет
внешний радиус R2, толщину ^>2) и радиус срединной поверхности
R (2) = Л + h0;2Y2. Ширина зазора 8, заполненного вязкой несжимаемой жидкостью значительно меньше Л15 Л2. Характерная длина 1 (длина волны или оболочки) значительно больше радиусов Я1, Л2.
Уравнения движения вязкой несжимаемой жидкости имеют вид (1). На границах оболочек и жидкости выполняются условия прилипания ди(г) — дШ(г)
Vx = V = —^^ при г = Л — Ш(г) (1=1,2) (11)
здесь и() - продольное упругое перемещение оболочек 1,2; Ш()- прогиб, положительный к центру кривизны оболочек 1,2; индекс 1=1 относится к внешней, а - 1=2 к внутренней оболочкам.
Напряжение со стороны жидкости определяются формулами (4). Принимая за характерную длину величину I (длина волны или половины длины оболочки) переходим к безразмерным переменным
г* = г-^ х* = Х ^ = ^t,Ш) = ^ ы^и(г) = и и/0,
с»& 7& 7& т 3 & т 1 &
V, = ^ ог V = WmC-88 Ох, Р = р0 + - =
Ро(1 —/О
у = «1 Л = 8 «V-г = £<а&-¡Г = °(£)>^Т = °(£1/2), 7 = °(£ (12)
Л2 8 I Л2 I I
Учитывая наличие малых параметров у и Л, получим уравнение динамики
жидкости и граничные условия
= оЯе ^ + ^ = ^Ёи = 0 (13)
у-ч * Э Л .* -Л * -Л *2 & -Л * -Л * V &
дг дt дх дг дг дх
ди, & ^ * ^ * ~ 8 8со
иг =--V, их = о при г = 1, г = 0;Яе = —.—1
Для гармонических по времени решений из (13), в которых вместо с0/1 положим с, получим
г г ди (2) ди (1) ~ д2и (2) д2и (1)
р = И ™ -иг - ^ + аКе( - Чг ^х-д —щн^—(14)
дг &г =1 4 дt * дt 2 4 дt 2 дt
дг * г=0 л дt * дt * & 2 V дt 2 д1 *2 &
где 282 = 81®/у = Яе - число Рейнольдса
, ч 1 _g3(shg— s1n g)_
у(с) = т - ,-:—--—---(15)
а(с) =
е[е(сЬе + cos £) — (^Ы — sin £)]
£ (ск£ + COS £) — 2£^£ + sin £) + (ск£ — COS £)
Отметим, что
а =1,2 и у=1 при£^0; а = 1, а-1 = 1/s, y_s/6 при s > 4 (16)
Решаем (13) методом итерации [1]. На первом шаге итерации найдём:
я=.2! [ I (zf-UW,
^ *2 *чг 2U3(2) д2U3(1&& * —x _ 6(r - r ) I (----—)dx ,
at ^ J at2 at2
Подставим во второе уравнение (13), на втором шаге итерации
находим
я=!! [12( Ul-Ц!)+dU--dUr-na-jx
JJ v dt * at * 5 v dt2 at2 n
duk|.1=! [6(£U^-dU^)+±Re(dU--dUV (17)
dr * lr_ JL v dt dt * & 10 V dt2 dt2 & V 7
U. o=! [6( Щ1 -duC)+—Re( dU--dUV
dr &r =0 J V dt dt * & 10 V dt *2 dt 2 }
Эти формулы позволяют найти нормальное и касательное напряжение на оболочках
>-1„ _( П г-Л _Г „ PV0Wm DV П i- •
ИГ qn _ (-1)г-1 qrr _ [-Р0 ^я](-1)
q (1) _ q (1) _ PVC0lWm | . (18)
_ _ §2 &r*(18)
q (2) _ q (2) _РуФт_ U |
qx qrx §2 dr * lr-_0
Как видно из (18) эти величины имеют противоположные значения на верхней и нижней оболочках.
Отметим, что характерное время t0 = I / с0 может быть заменено на частоту 1/а>в случае колебательного процесса, а не волнового.
Если на втором шаге итерации ввести поправочные частотозависимые коэффициенты а перед локальным членом и у - перед силами трения, то формулы (14) совпадает с (18). Это совпадение является формальным, но (14) при этом можно использовать и для негармонических волновых процессов при с^ /1 = ш. Этот случай включает в себя и нелинейные уравнения оболочек.
Надо только учесть что
у = 1, а = —,-= — при а << 1.
Работа выполнена при финансовой поддержке гранта РФФИ № 16-01-00175a.