Лойцянский Л. Г. - Механика жидкости и газа (1067432), страница 14
Текст из файла (страница 14)
1д1~ 1 /ды дих = 0 + тилеф (О! 1'е! 0) !11= 2 ( — + — ) «я~й; 2 (, д«дг )о в) д л я т о ч к и Мо. хо = 0 -+ илеф(0, О, ля)а!г= — 1 — + — ) гзо!Г 1 /ди, ди!л 2, дг дх)о 1 /ды дот «о=о + Флеф(0! О! гь)Ж= — ( + — ) гоЖ, 2 (,ду дг)о го =аз+ шлеф(0 О, гз) геле!= ха+ (д ) гог1е'=го ~1 +( — ) Ю). Составим теперь, например, сноросогь относительного удлинения отрезка МоМ„' это будет, с точностью до малых высших порядков: М М ММ )гУ х (!+( — ) дг1 +...— х! МоМ! хе де (д е о аналогично будет и для других диагональных колепонент. В начальный момент 1 угол между осями Ох и Оу был равен 90', косинус его — нулю.
Обозначим через т „уменьшение (скошение) этого угла к моменту 1+И; тогда получим, в силу бесконечной малости угла у — е., со$ (МоМ! МоМл) — соя ( ухи) = $1п уо!в = или, используя известную связь косинуса угла между двумя напРавлениями с их направляющими косинусами, в данном случае с точностью до малых высшего порядка, равными: х! =- хл + илеф (хн у, = 0 + ол,ф (х„ г! = 0 + гольф (х1 О, О)дг=х,+( д") х,ж=х,~1+®1 дГЗ, О, О)д1= — ( — -+ ) .,,11, 1 где ди' 2 (дх ду), О, О)дс= — '(ди ( д ) х,д1,. дх о ЗЛЕМЕНТЫ ТЕОРИИ ПОЛЯ. КИНЕМАТИКА СРЕДЫ (гл. ! а) для направления М Мг, х, у, 1 сди дих г, 1 тди дэ1 хг ' х! 2 дх ду)о ' х1 2 ~де дх)о г б) для направления М,Мг! х., 1 сдо ди~ у., уг 2 дх ду)о ' уг получим = ! ' — ( — + х-) с!г+ — ( — + †) с!! 1 + беск.
мал. 2-го пор. 1 lди дих 1 сдо дих 2 (,дх ду)е 2 (,дх ду)е 1 дге да Окончательно найдем (опуская индекс нуль) для скорости т сношения угла лОу: и аналогичные формулы для других направлений. Деформационная часть поля скоростей характеризуется еще одной важной величиной — сноростью относительного обяежного расширения в данной точке, которую можно определить как предел 1 Ы 1пп — — (Ьт), АА где дт — малый объем, в которол~ взята точка.
Эта физическая скалярная величина носит наименование дивергенцаи (расходимости) скоростного поля и обозначается символом б!чУ, так что можно на- писать б!ч У = !!ш — — (~Ь). 1 ьт.+ еач "г (58) Чтобы не смешивать бесконечно малые приращения в пространстве с приращениями во времени, условимся в тех случаях, когда зто может повести к недоразумениям, обозначать пространственные дифференциалы символом 8, временные — ь!. Тогда (58) дает б1ч У = — — (3т). 1 Ея И (59) В дальнейшем часто придется находить производную по времени от элементарного объема жидкости.
По (59) имеем: й иг — (дт) = б!ч У ° 8т. (59') 9 11. Скорость объемного расширения жидкости. Интегральные представления дифференциальных операторов поля. Основные интегральные формулы $11! скОРОсть Овъемного РАсшиРенйя жидкости Для определения величины 61РЧ воспользуемся приемом, идея которого восходит еше к Эйлеру.
Возьмем в данный момент времени 1 элементарную трубку тока и двумя произвольно наклонными сечениЯми (Рис. 9) ~Го, и г1ао выделим некотоРый объем АВСА1 = 8т. За время Ж объем сместится в положение А'В'С'Ю' и, вообще говоря, изменится, причем, как легко сообразить, полное изменение объема АВС0 в этом случае будет равно: г1 (8т) = объем А'В'С'В' — объем АВС0 = = объем ИЭ'С'С вЂ” объем АЛ'В'В, так как объем А'В'СО является общей частью объема трубки в начальный и следующий моменты. Проведем внешние по отношению к объему АВСЮ нормали п, и по и внутреннюю нормаль пг, а также отметим векторы скоростей Ч, и Уя в сечениях Ио, и 1ао Тогда будем иметь; объем АА'В'В =с1а, ° л, = = — Фо, ° (г, г1г ° соз (Ч„п1) =- = — 1г,соз(Чп и,) г11 ° Иа„ объем ВО'С'С = ~Ья ло = 1' соя(Чо, п~ Ж ° газ и, следовательно, о (Ьт) — = Ч„,Но, + 1гя„гЬЕ.
(60) Возьмем теперь в поле скоростей любой конечный объем;, разобьем его поверхностями трубок тока на бесчисленное множество элементарных объемов 8т; при этом входные и выходные сечения г1а заполнят всю поверх- Рнс. 9. ность а, ограничивающую объем Просуммируем равенства, подобные (60), по всем трубкам тока; тогда, очевидно, получим обшую формулу для любого конечного объема: ич (' р ч. — — ) Уч~й= ~ 1'соз(Ч, п)а~о= ~ п ° Чда. (61) Ф о Согласно (58), получим теперь следующее интегральное представление дивергенца скорости: 61Р Ч = !пп — ~ 1г„~й = !Нп — ~ п ° Ч гЬ, (62) а..+ оа' ., ао.о о~о . а.
ао ЭЛРМЕНТЫ ТЕОРИИ ПОЛЯ. КИНЕМАТИКА СРЕДЫ (гл. ! где г!е — поверхность, ограничиваюп!ая малый объем Ьт, заключающий в себе точку, в которой определяется ц!ч Ч; при стремлении !!т к нулю поверхность дч стягивается в эту точку. Замечая, что выражение !Т„бч представляет секундный объемный расход жидкости сквозь замкнутую поверхность Ьч, ограничивающую объем Ьт. содержащий внутри себя точку, в которой определится дивергенция, можем еще определить величину б!ч Ч как предел отношении секундного обеемного расхода вкидкости сквозь замкнутую поверхность к обеему, ограниченному этой поверхностью, когда поверхность стягивается к точке, в которой определяется дивергенция.
Как видно из хода доказательства, объем Ьт совершенно произволен по Рис. 10. форме. Выберем за ит элементарный декартов координатный параллелепипед (рис. 10); тогда, составляя непосредственно поверхностный интеграл в правой части (62) от значения Уч по всем шести граням (эти значения показаны на рисунке), получим: б!ч Ч = !Нп — — — р и + — Ьх) Ьуйе — идущее+ ~~~ ( + !о + — Ьу) Ьх Ьг — о Ьх Ьи + до ду дм +(го+ — Ьг) Ахи — шйх Ау+ б. м.
выс. Пор.~. откуда найдем искомое выражение дивергенции скорости в прямоугольных декартовых координатах: ди до дге б!ТЧ= — + — + —. дх ду дг ' По заданным уравнениям поля скоростей й!ч Ч в данный момент легко вычисляется. Для облегчения запоминания формулы (63) приведем интггвальныв ьоюжлы ее символический зид б!ч Ч -:= Р 11. (64) Раскрывая правую часть по правилу скалярного произведения и вспоминая (24), получим формулу (63).
Заметим, что дивергенция может рассматриваться как некоторая общая дифференциальная операция, совершаемая над любой векторной функцией и определяемая формулами (62) и (63), куда вместо вектора Ч надлежит вставить дифференцируемый вектор и. В этом случае уже нельзя говорить о скорости обьемного расширения, а вы- ражение ) а„ь!е= ( и ° ас!о, Ю а б!та= !пп — з! а„с!в= 1!гп — ~ п ° аь1в. (62') т.+ь дт ь;-ьь ич ьа Аа Выражение дивергенции вектора а в декартовых прямоугольных координатах будет, аналогично (63), иметь вид: 61ча х+ У+ ь дан да да дх ду дв ' (63') Приведенный ранее вывод формулы (63) почти буквально можно повторить для элельента объема в любой системе криволинейных координат (полярных, цилиндрических, сферических и др.) и получить, таким образом, выражение днвергенцин вектора-функции в криволинейной системе координат; это будет сделано далее в гл. 711.
Из формулы (62') легко выводится важная для дальнейшего интегральная формула, впервые указанная в 1834 г. знаменитым русским академиком М. В. Остроградским (1801 †18). Разобьем любой конечный объем т на большое число малых объемов Ьт; обозначим поверхность, ограничивающую т, через а, а амтв через Ьв. б з . пмь л. г. л ь аь. где интеграл берется по некоторой (вообще говоря, не замкнутой) поверхности е, называют потока,и вектора а через поверхность е, Если вектор а представляет скорость жидкости, то поток вектора а совпадает с объемом ькидкости, протекающим через поверхность е в единицу времени, т. е.
с секундным объемным расходом жидкости сквозь сечение а, что приводит к ранее данному определению дивергенции скорости. В общем случае дивергенция вектора определяется как предел отношения потока вектора сквозь замкнутую поверхность к обьему, ограниченному этой поверхностью, когда поверхность стягивается к точке, в которой определяется дивергенция, т. е.
элвмвнты твоРии поля. кинвмАтикА сРвды (гл. ! Тогда, согласно (62'), будем иметь для элементарного объема Д;; 66 п ° а е!ч = б!ча ° Де+а Дтч Г (65) Ьа где е — малая величина, идущая к нулю с уменьшением Дг. Просуммируем обе части равенства (65) по всем объемам Дтч образующим конечный объем; получим: ~~.", ~ п .
а е!е = ~~" 6!ч а Дт+ ~чр ~е . Дт. Ьа В сумме, стоящей слева, взаимно сократятся все элементы интегралов, взятые по общим границам двух соседних малых объемов, так как сама вектор-функция а в силу непрерывности имеет одинаковое значение на границе со стороны какого объема не совершался бы подход к граничащей поверхности, в то же время внешняя нормаль к поверхности, ограничивающей один из малых объемов, является внутренней нормалью к той же поверхности лля смежного малого объема; поэтому в рассматриваемой сумке часть слагаемых, равных между собою по вели!ине и противоположных по знаку, сократится.
Останутся лишь элементы интеграла, распространенные по внешним частям поверхности е, окружавощей объев! т, т. е. поверхностный интеграл по поверхности е. С правой стороны, если устремить к нулю малые объеввы Дт, останется объемный интеграл от б!Ра, взятый по объему е, так как второй член справа, как сумма малых четвертого порядка, обратится в нуль. Таким образом, получим интегральную Формулу. ~ и ае!ч=~ 6!Раач (66) или ) а„а!е = ) д!Раавт.