УДК 532.543
О ФОРМАЛЬНОМ СОВПАДЕНИИ УРАВНЕНИЙ ДВИЖЕНИЯ ДВУХМЕРНЫХ В ПЛАНЕ ОТКРЫТЫХ ВОДНЫХ ПОТОКОВ И УРАВНЕНИЙ ДВИЖЕНИЯ ИДЕАЛЬНОГО ГАЗА
© 2008 г. В.Н. Коханенко, И.В. Папченко, Н.Г. Папченко
Донской государственный аграрный университет Donskoy State Agrarian University
В работе приведена аналогия между уравнением движения двухмерных в плане открытых водных потоков и уравнением движения идеального газа. Этот факт даст широкие перспективы для теории и практического развития аналитических решений граничных задач по течению двухмерных в плане открытых водных потоков. Он позволит хорошо разработанные методы моделирования в газовой динамике заимствовать для решения сопоставимых задач по течению двухмерных в плане открытых водных потоков
Ключевые слова: двухмерный открытый водный поток, потенциальность потока, интеграл Бернулли, плоскость годографа скорости, идеальный газ, адиабатный процесс.
The research work represents the analogy between the moving bimeasured equations for water stream and moving ideal gas equations. This fact gives a good opportunity for theory and practice of analytical decisions of border problems in connection with the moving binary measured equations for open water streams. It let adopt well worked out simulation methods at gas dynamic to decide confronted questions in moving binary measured equations for water stream.
Keywords: binary measured open water flow, potential of the flow, integral bernoulli, plane godograph velocities, ideal gas, adiabatic process.
Развитие возможностей вычислительной техники и появление мощных математических пакетов программ создало возможности для быстрого построения и апробирования различной сложности математических моделей по течению водных потоков. При этом наибольший интерес представляют собой наиболее простые модели, дающие, однако, удовлетворительные по степени адекватности результаты для практики. Известно, что построение математической модели какого-либо физического процесса начинается с выбора системы уравнений, описывающей сам процесс. Часто математически формально система уравнений исследуемого процесса по внешнему виду сводится к уже известной системе, решение краевых задач на базе которой также известно и его можно заимствовать для решения новой исследуемой задачи.
Целью настоящей работы является приведение известных уравнений движения идеального газа и уравнений движения, двухмерных в плане открытых водных потоков, при определенных ограничениях на поток к сопоставимому виду и демонстрации их совпадения.
Рассмотрим вначале уравнения движения двухмерных в плане открытых водных потоков. Для этого воспользуемся в качестве исходной системы известной системой двухмерных в плане уравнений движения открытого водного потока в случае горизонтального плоского отводящего русла, в которой не учитываются силы трения сопротивления потоку [1]:
du du dh u--hu--h g — = 0;
dx dy dx
du du dh u--hu--h g— = 0;
dx dy dy
dd — ( hu )+ —( hu) = 0,
dx dy
(1)
где и, и - осредненные по глубине потока проекции местной скорости; h - местная глубина потока; Оху -плоская прямоугольная декартова система координат (в плоскости течения потока); Ох - продольная ось симметрии потока; Оу - дополняет ось Ох до правой системы координат; g - ускорение силы тяжести.
Система уравнений в частных производных (1) является, в соответствии с общепринятой в математической литературе [2] классификацией, существенно нелинейной системой, и поиск ее регулярных решений представляет определенные трудности. Поэтому для ее упрощения примем дополнительные условия, а именно потенциальность потока, т.е. существование функции ф = ф(х,у) такой, что выполнялись бы равенства:
d< dx
d< dy
= u.
(2)
u
Линия ф = ф(х,y) = const называется эквипотенциальной. Тогда с учетом условий (2) система уравнений (1) трансформируется к виду:
V2
--ъ h = H = const;
2 g
д д — ( hu )+ —(ho) = 0,
дх ' дУ
(3)
где V2 = и2 + и2; а постоянная «Н» определяется, если известны модуль вектора скорости - У0 и глубина потока - h0 в некоторой характерной его точке
у 2
по формуле: Н = Н 0 = —— + h 0.
2 Я
Первое уравнение системы (3) представляет собой интеграл Бернулли для двухмерных в плане водных потоков. Второе уравнение системы (3) является уравнением неразрывности потока и определяет существование функции тока у = у(х,у), такой, что выполняются равенства:
h дш —u =—'
дУ
h дш
h дх
Линия
ш(Х У ) = C
(4)
(5)
называется линией тока, которая обладает свойством непроницаемости потока вдоль нее. Из условия (5) следует:
, ду ду d у = —— dx +—— dy = дх ду
= — — иdx + — udy = — (-udx + udy) = 0,
При совпадении направления векторов dS и V . dy и
выполняется условие tg6 = — = —, из которого сле-
dx и
дует условие (6).
Линии тока перпендикулярны эквипотенциалам в физической плоскости. Это следует из уравнения:
^ф ^^г+^7 +*,- ) =
uhu uhu - +-= 0.
h
h
Итак, систему уравнений (2) - (4) будем считать базовой для двухмерных в плане открытых водных потоков. Перепишем ее еще раз для лучшего восприятия по тексту работы:
u=
дф дх
u =
дф
~дду'
V2 V02
h = H0 - V H0 = h0 +Vg;
(7)
h дш —u = -
дУ
h дш
h дх
Введем понятие квадрата скоростного коэффици-
V2 2
ента х = -—г— = X и преобразуем систему (7) к виду:
2 gH с
дф дх
= Но (1 -т); h =
дф
ду'
h о (1 -т).
-т
(8)
h дш
—u = —; h
ду
h дш
-u = --
h
дх
, dx dy
т.е. udx = udy, — = — .
u u
(6)
Условие (6) показывает, что касательная к линии тока совпадает с направлением вектора скорости жидкой частицы в рассматриваемой точке потока (рисунок).
Из системы (8) получим зависимость между функциями ф и у , характеризуемую формулами:
дш 1 - т дф
ду 1 - т 0 дх
дш 1 - т дф
дх 1 -т0 ду
(9)
Преобразуем далее систему (9), переходя в плоскость годографа скорости, т.е. к переменным (х,6) [3]. В результате имеем систему
Схема к пояснению понятия линии тока
дф 9т
hr
1- 3т дш
[ 0 т
(1 -т)2
(10)
дф = 2 К__т дш
д9 Нп 1 -т дт .
h
0
0
u
и =
0
0
0
X
Для удобства введем новую функцию тока
у = —— у, отличающуюся от у(х, у) постоянным н о
множителем h0/Н0 . В результате система (10) преобразуется к виду
Зф Зт
1 1 - 3т
2
Зу
(1 -т)2 зе
т Зу
(11)
5ф = 2 зе 1 -т Зт
Или, убирая волну над у в системе (11), окончательно получим систему уравнений двухмерных в плане потоков в плоскости годографа скорости:
Зф Зт
1 1 - 3т
2
Зу
(1 -т)2 Зе
т Зу
(12)
— = 2 Зе 1-т Зт
Системе (12) эквивалентно уравнение в частных производных второго порядка:
З [ 2т Зу 1 1 - 3т З 2у
-2—2- = 0.
2т(1 -т)2 Зе2
Зт 11 -т Зт
(13)
З 2у т—i- + -
1 Зу 1 - 3т З 2у
= 0.
(14)
Зт2 1 -т Зт 4т(1 -т) Зе2 Из сравнения уравнения (14) с уравнением общего
[2] видно, что
,З 2у „„З 2у ^З 2у „ вида A—т- + 2В—- + C—2 = 0
Зт2
A = т, В = 0
ЗтЗе
с =
Зе2 1 - 3т 4т (1 -т)
Д = В2 -AC =
-т(1 - 3т) 3т-1 п 1
= —^-г- = —-- < 0 , при 0 < т < -, и уравнение
4т (1 -т) 4 (1 — т) 3
(14) относится к эллиптическому типу.
Далее рассмотрим идеальный газ. Потенциальное движение газа согласно [4] требует для адиабатного процесса выполнение равенства:
p = kpy, т.е. p = p(p).
(15)
где p - давление газа; k - коэффициент пропорциональности; у - показатель адиабаты; p - плотность газа.
А.С. Чаплыгин рассматривал движение газа с сохранением тепла частицами газа ввиду малой теплопроводности и малой лучеиспускательной способности газа.
Для потенциала скоростей имеем функцию ф = ф( У)
Зф Зф u = — , и = — .
Зх Зу
Уравнение неразрывности записывается в виде
фи + фи = 0
Зх Зу
(16)
Для спокойных потоков 0 < т < —; для бурных потоков — < т < 1 .
3
Уравнения (12), (13) решают вопрос о течении жидкости, если известна область переменных (т, 6),
соответствующая этому течению, если даны значения «у» на граничных линиях тока; если везде, внутри области, функция «у» вместе со своими первыми производными конечна, однозначна и непрерывна, а количество «т» не превышает 1/3 и обращается в нуль лишь в некоторых точках контура для спокойных потоков и превышает 1/3 и обращается в единицу лишь в некоторых точках контура для бурных потоков. Для спокойных потоков уравнение (13) относится к эллиптическому типу. Оно переписывается в раскрытом виде:
Закон Бернулли с использованием уравнения (15) приводит к соотношению
Р = Р о
( V 2 ^
1--
2а
(17)
ТЛ2 2 2 kур 0
где V = u + и2, а =
у-1
у-1 _ 1
Р =-7 , Р 0 = const.
у-Г
При этом р 0 - плотность газа в точке с V = 0 .
Р
V 2
Полагая т = ^— [4] из (17) получим: р = р 0 (1 -т)
В литературе [5] эта же формула приводится в виде
__1_
р 0 I у -1 2 | Т-'
—0 = 11--X I . Из уравнения (16) следует
Р I У +1 )
существование функции у , определяемой равенствами:
'-Р-u = Р 0 Зу
Зу '
[Р 0 Зу Зх
Далее, переходя в плоскость годографа вектора скорости, получим уравнения С.А. Чаплыгина для идеального газа [4] в виде:
и
Зф
"Зё
2т Зу
(1 -т)
ß Зт
Зф 1 -(2ß + 1)т Зу ;
2т (1 -т)Р+: 30 ^
(18)
2т Зу | + 1 -(2ß + 1)т З 2у _ 0
-I-
Зт[(1 -т)Р Зт
2т(1 -т)
ß+1
З92
(19)
С.А. Чаплыгин в своих исследованиях рассматривал случай 0 < х < ^ 1 ^ дозвукового течения совершенного газа и разработал методы решения ряда практических задач по течению идеального газа при его дозвуковых режимах.
Как видно из сравнения уравнений (12), (13) и (18), (19), они совпадают при р = 1. Следовательно, в плоскости годографа вектора скорости при одних и тех же граничных условиях решения соответствующих задач для газа при р = 1 и двухмерного в плане спокойного потока будут совпадать. Этот факт дает широкие перспективы для теории и практического развития аналитических решений граничных задач по течению двухмерных в плане открытых водных потоков. Он позволит хорошо разработанные методы моделирования в газовой динамике заимствовать для
решения сопоставимых задач по течению двухмерных в плане открытых водных потоков.
В заключение работы приведем для интересующихся исследователей формулу, связывающую параметры потока в плоскости годографа скорости жидких частиц с параметрами потока в физической плоскости течения потока в дифференциальном виде:
dz = | dф + i — dу |—е'6 , ^ h ^ V '
где z = х + 'у , позволяющую, имея решение граничной задачи в плоскости (х, 6), перейти в плоскость (х, у), т.е. полностью решить соответствующую задачу.
Литература
1. Емцев Б.Т. Двухмерные бурные потоки. М., 1967.
2. Справочник по математике для научных работников и инженеров / Под общ. ред. Г. Корн, Т. Корн. М., 1970.
3. Коханенко В.Н., Волосухин Я.В., Ширяев В.В., Коха-ненко Н.В. Моделирование одномерных и двухмерных открытых водных потоков / Под общей ред. В.Н. Коханенко. Ростов н/Д, 2007.
4. Чаплыгин С.А. Избр. тр. Механика жидкости и газа. Математика. Общая механика. М., 1976.
5. Лойцянский Л.Г. Механика жидкости и газа: 5-е изд. М.,
1978.
Поступила в редакцию 9 июня 2008 г.
Коханенко Виктор Николаевич - д-р техн. наук, профессор, зав. кафедрой механики и гидравлики Донского государственного аграрного университета, пос. Персиановский. Тел.: 3-55-21.
Папченко Игорь Владимирович - ассистент кафедры механики и гидравлики Донского государственного аграрного университета, пос. Персиановский.
Папченко Наталья Геннадиевна - ассистент кафедры механики и гидравлики Донского государственного аграрного университета, пос. Персиановский.