Н.А. Слёзкин - Динамика вязкой несжимаемой жидкости (1159534), страница 17
Текст из файла (страница 17)
В дзнном случае все елиничные векторы изменяют своь направление при переходе нз Рнс. 20. одной точки в лружую, если этот перехол связан с изменением двух координат 0 и 2. Вектор угловой ско- рости поворота с изменением угла 0 будет направлен по касатель- ной к коорлипатной линии у, т, е. Вектор жс угловой скорости поворота с изменением угла ф будет направлен по оси Е, и поэтому булем иметь: ы = — (соэ0г — э!п оюя). дт Подставляя в уравнение (3.1) выражения (3.2), (3.4) и (3.6) и собирая коэффициенты слева и справа при одинаковых единичных векторах, получим слелующие лифференциальиые уравнения лвиження сплошной среды в цилинлрических координатах в коипоненгах напряжений: 82 лиеевгвнцилльные твлвнения лвижвння вязкой жидкости ьгл.
и Учитывая эти равенства и равенства для производных от единичных векторов по времени, получим: дЮ„ дв =ЮзХ(з= — '=(соз Ою,— з(п 01з) Хю, = з)п 0)з, дЮ, з дт — ю', — ' = (соз Мю — зюп 01) Х1я= соя 0(з, дю, т О, д - — — (соз ОЯ,— з)п ОЩ(ю~= — соз 0)з — ~1пОКы двз в дЮ дз = юзХ~г = з (3.9) дюз дз =ююХ"ю= На основании этих равенств булем иметь: до' д д— ,„= дз (ол(в+оввз+ отюз) =- дол доь до„ дз з+ дз з+ дз з+ из дР д д д (олью+пью' ~ о Ч— Чз т т дол дов . до = — г + — ю(+- — тю(+ дт ' дт з ду + ов 5!и юююз+ов соя 01з — о, соз юьюз о зьп Оюь д д д— „,Жив) = дз 1Ю зььь 0(Рвлюв+ Рььюэ+ Рвь(з)1 = (3ПО) = ьь соз 0 (Рввюз+ Рввюз+ ЮзвтЮз) + дль дльз . дль, .
+)~"" 0(дз ' + дз ' + дз'"юз+Р ' — Р4) д д дй ()зь)тьРз) = д. М(Ртл)ю+Рььюю+Р юЛ = Юдоьв . ОРьь дР +Р всозОЮж — Юэ сов бюз — Р„„эьп Ою ), Подставляя в уравнение (3.1) выражения (3,2), (3.8) и (3.10) и приравнивая коэффициенты при одинаковых елиничных векторах слева ь и справа, получим следующие дифференциальные уравнения вижеь ия сплошной среды в сферических кооодинатах в компонентах л гвлвнкнив пеевносл полной знеегии напряжений: дерк дек е д'рк г дпк ее -(- ез дт дк( А' дз РерпО дт р — к+од — + — ' — + — т — — — т = 1 гдркв 1 дРьк 1 'др в = рл+ — ~ — + — — + — — ь Г !( дв0 я, дз рзрпО д 1 + ~ (2рлк+ 0!К Орьн — р„— р де, деь е де„е де+в е ц~ дт дЛ' ' А' д0 )серп О дэ Аь др — +ел — -'- — ="+ — т — —,+ — ' — — ес)сО = 1 МИко 1 дкм 1 д.е, Рь+ 1 + + — — те+ Р д0 10е1п 0 дт 1 + —,(рмссаб+Зр,„— р с)ЕО)1, де еь дет е'т де е ев — '+ — -+ — — + — — + + — ' дг дЙ д дв де1пО дт рз лр 1 гдркр ! дрьр 1 ь ( д)0 )О дв ' 10 зря О дт + — (2р 01к О+ Зр л)~, (зп !) й 4.
Уравненне переноса полной энергии Переходя к выводу уравнения изменения энергии в фиксированном элементарном объеме, заметим, что в термодинамике под внутренней энергией системы полразумевается та часть полной энергии, которая зависит от температуры, объдма и химического состава системы, при этом, если пренебрегать энергией взаимодействия частиц системы друг с другом, то внутренняя энергия обладает свойством аддитивности, н поэтому можно ввести понятие удельной внугвренней энергии е, представрряюшей внутреннюю энергию единицы массы. Если мы будем рассматривать фиксированный объем, не изменяюрцийся во времени, то полная удельная энергия елиницы массы будет 1 состоять из кинетической энергии — (ря и внутренней энергии е. 2 Изменение полной энергии массы в фиксированном малом объеме за малый промежуток времени Ы будет составляться из о~дельных изменений за счет: !) входа и выхода масс через границы обьвма, 2) элементарной работы объвмной силы Р, 3) элементарной работы векторов напряжений р,, ре, рз и 4) притока тепла благодаря теплопроводности, г(ругие источники изменения полной энергии (излученне и пр.) мы учитывать не будем.
Подсчитаем отдельные изменения полйой энергии в фиксированном параллелепипеде с ребраии Оэо еэа, Оэв. 34 диФФеРенциАльные УРАвнениЯ дВижениЯ ВЯзкОЙ инакости Ггл. и В момент Г масса, содержащаяся В фиксированном объвме Ьгеазйзз, бУлем иметь полпУю энеРгию, РавнУю [( —, + а) р] Н,Н Нлйд, йд ед,, а в момент Г-+ЬГ булет иметь: р( —,+з)] Н Н Няйд,йдаьда= = ~ [р( — +-а)] + — [р ( — +а)] бт+... ) Н НЕН Зд,йдаЗда. [Гох ( 2 + В)НВНА] едэ еда йт. Через противоположную грань, проходящую через точку с коорлинатами д +эд,, д, дз, выходящая масса вынесет количество полной энергии, равное [ГО,НАНА ( — + В)] едя ьдз — ( [зп, ( — + е) Н.,НА] + + — [Ро ( — + е) НЯНЯ]йд,+ ... ~едэйдзМ.
Следовательно, внутри параллелепипеда задержится следующее коли- чество полной энергии: — — [ро ( — +з) Н. Н ~йд едчедаМà — .. г4. 2) Повторяя такие же рассуждения по отношению к граням, перпе~гдикулярным к касзтельным к координатным линиям дя и д, получим: — — ) РВ, )--+ а~ Н НЯ1 йд 'дяйдя бт —... 44.3) Складывая выражения (4.'2) и 44.3), получим приращение полной энергии в фиксированном объеме за счет Входа и выхода масс через Следовательно, приращение полной энергии в фиксированном объеие за промежуток времени ОГ представится в виде Л вЂ” )Р~ — +яд ЬГНЕНВНаедгьдяйдз г Через грань, перпенликулярную к касательной к координатной линни д, и прохолящую через точку О с координатами ди дю дз входящая масса ро,НЯНзодяйдт внесет с собой в параллелепипел следующее количество полной энергии; 86 диеевтянцилльныв ягавняння лвижяния низкой жидкости (гл.
и Через перелнюю грань, перпендикулярную к касательной к координатной, линии д,, за промежуток времени Ы будет передано по закону Фурье следующее количество тепла: — (х — Нзйз) йд, йда Ьд дТ Над, з в з а Через противоположную грань за тот же промежуток времени будет передано количество тепла, равное — (х — — — а) ОдзВдадт = / дТ НаНа~ ддг Нг )я,чля, Следовательно, внутри параллелепипела задержится следующее коли чество тепла: д I дТ НаНя' — (х — а ' ) йд„йд. бдяща+... (4.9) вдг'х дд, Н, Т Повторяя такие же рассуждения по отношению к остальныи граням параллелепипела, получим: д ~ дТНН' — — 3 1)3д Од од, ба+..., 1 (х — ' а) йд,йд йдзбт+... д дТ НчНа ада ада Н, (4.10) Складывая выражения (4.9) н (4.10) н деля на термический эквивалент А, получим то приращение полной энергии в фиксированном объйме, которое обусловлено процессом теплопроволностн: дда Х дда На ) Других источников изиепения полной энергии в рассматриваемом объйме яет, поэтому приращение энергии, представленное выражением (4.1), мы должны приравнять сумме отдельных приращений, представленных формулами (4.4), (4.5), (4.8) и (4.11).
Обе части получешюго равенства разделим на НгНлНдбд,одайдабт и перейддм к пределу, стягивая параллелепипед в точку и уменьшая промежуток времени бс до нуля. В результате получим 9 41 УРАВНЕНИИ ПЕРЕНОСА ПОЛНОЙ ЗНЕРГИИ следующее уравнение изменения полной энергии в фиксированной точке области, занятой средой: + — '['(-;-'+2) ° .|+[Роз(2-+ ) Ч.)+ г Н2НВНА'Где х ' е г ' двг (4.12) Уравнение (4.12) можно также называть уравнением переноса полной энергии.
Оно в своей простейшей форме было ввелено впервые в рассмотрение Н. А. Умовым ') в 18?3 г. Группируя слагаемые в правой части (4,12), пояучим: =Рр' )г ННРН ( д? [(Рох 2 +рог — Рх ' 1' — АН до ) НВНА)+ +дог[(р 2 2 +р г ре АН д ) г 21+ дГ/ Уг х дТА + — [(ро —;+р *- —,о. У вЂ” — — ) НН1). (4,13) дог)1 2 2 2 Ангдов) Вырагкение в фигурной скобке в правой части (4.!3) представляет собой ливергенцию особого вектора, который можно назвать вектором плотности потока переноса полной энергии. Обозначая этот вектор через Е, для его компонент будем иметь следующие выражения; (4.14) В Умов Н.
А„ Избранные сочинения, Гостехнздат, 1950. ! )гг "="' (-+з)— г Уг "вг2 ! ) Ее — — Роз ( 2 + ) х дТ У вЂ” — —, АН, двх ' А дТ АН2 д?2 ' дт АНгдвг При этих обозначениях уравнение переноса полной энергии (4.13) представится в виде д ~Р( — + )] =-РЕ. — ~ — (ЕзНгНг)+ В декартовых координатах уравнение переноса полной энергии (4,15) будет иметь вид д дГ дбз дЕ, дг(Р(2 + )) Р (дх+ду+д ) где нооекции Е , Е„ и Е вектора плотности потока полной энергии равны /Уз х хдТ Е =Ри( — +з) — р У вЂ” — —, (,2 ) з А дх' 7У .ду Е =Резь-+зз — -р У вЂ” — ' —, з ~2 *) з А ду' грг Х гду Е = Рсз( — + з)1 -- тз У вЂ” — —. 1,2,) " "А дз' (4.17) ф 5. Уравнение изменения внутренней энергии Преобразуем полученное уравнение (4.12) переноса полной энер- гии, Так как — = — У ° У, 2 то уравнение (4,!2) можно предстанить в виде ~дУ и, дУ егдУ оздУ „~ з ( г ) з 1 дг Нзддз Нгддз Нз ддз + — (Р.ННг)))+Р~-'+ — — + — — + — =) = ддз з ' з ) ) 'ьдг Нзддз Нзддг ' Нзддз! дУ дУ, др 1 (д / дХНгН„'~ Нзддз з Нзддз ' з Нзсдз ' АНзНзНз(ддз( ддз Нз 7 „д ( ду НаНз) д ~ ду 77,Нз)~ 88 диеевввнцилльныв ттлвнвпия движения вязкой жидкости (гл.