Слёзкин Н.А. Динамика вязкой несжимаемой жидкости (1123892), страница 17
Текст из файла (страница 17)
Уравненне переноса полной энергии Переходя к выводу уравнения изменения энергии в фиксированном элементарном объеме, заметим, что в термодинамике под внутренней энергией системы полразумевается та часть полной энергии, которая зависит от температуры, объдма и химического состава системы, при этом, если пренебрегать энергией взаимодействия частиц системы друг с другом, то внутренняя энергия обладает свойством аддитивности, н поэтому можно ввести понятие удельной внугвренней энергии е, представрряюшей внутреннюю энергию единицы массы. Если мы будем рассматривать фиксированный объем, не изменяюрцийся во времени, то полная удельная энергия елиницы массы будет 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 диеевввнцилльныв ттлвнвпия движения вязкой жидкости (гл.
и 89 й 5! гелвнение изменания внттгеннвй энаггии Первая фигурная скобка в силу уравнения (3.!), вторая — на основании уравнения (1,8) обращаются в нуль. В результате, получим уравнение др дУ ! др 1 Гд /дТН»Н»! Ндд +' е Н.деа ' Рэ Ндаз + ЛИНН Где ( до П )+ Тзк как выражение в скобке в левой части (5.1) представляет собой индивидуальную произнодную от внутренней энергии фиксированной частицы, то полученное уравнение есть уривнение изменения внутренней энергии фиксированной чистииы с аостоянной массой. В декартовых координатах уравнение (5.1) изменения внутренней энергии представится в зиле Гд» де д» д»т д У д У д У р( —.+и — +о — +ш — )=р .— +Р . — +р .— + (,дс дх ду дх,) х дх в ду ' де +.-~д( -)-.-(Я.-)+-.
(х-.-)1 Раскрывая скалярные произведения трах первых слагаемых в праной части (5.2), получим; дУ дУ, ддр ди да дш дх ! е ду+' " д» Р™дх+Р»едх+~х» дх ди, до дш ди да дж + ду~~Р»ед +Ре д +! хд +~ш! +~»» Учитывая соотношения взаимности касательных напряжений (10.17) главы 1 и обозначения компонент скоростей деформаций (5.5) главы 1, будем иметь: др дУ, дУ Рх ' дх+Р» ' ду ~ Р» д. =р .л +ревене+р е, +2р „е, + 2р»»„»+2р е, . (5.3) С помощью соотношений (11.!), (11.16) главы 1, представляющих обобщенную гипотезу Ньютона лля вязкой жидкости, равенство (5.3) можно записать; дУ дУ дУ Р .
— +.а — +Р» ° — — = дх ' в ду д" 2и = — Рй — '(к' — — ~) бз+ 2!» (=-'з +ее +»е -1-2(ее +»з +»'„)). (5.4) 90 дифеееенцилльныв еелвнения движения вязкой жидкости (гл. и Таким образом, уравнение (5.2) изменения внутренней энергии фиксировинной частицы вязкой ясидкости представится в виле /де де де дез р( — +и — +о — +ш — )= 'Хдг дх ду дгг' +(х' — --")У+29(ег +ез +ге +2(гг -+ез +ег )). (5.5) Уравнение (5.5) можно рассматривать как уравнение притока внутренней энергии зо единицу времени в фиксированной частице вязкой жидкости. Источниками изменения внутренней энергии частицы вязкой жидкости, таким образом, будут: 1) теплота, поступающая благодаря процессу теплопроводности, 2) работа сил давлений, связанная с изменением плотности частиц, и 3) некоторая часть работы вязких напряжений.
Для случая так называемого совершенного (идеального) газа внутренняя энергия единицы массы равна с зт е А ' (5.6) Принимая теплоемкость с„ постоянной и подставляя значение е в уравнение (5.5), получим следующее уравнение притока теала для совершенного вязкого газа: + — (л — )~ — рй+()У вЂ” ~~) ба+ 2р ~( — ) +( — ) +( — ) + й 6. Дифференциальные уравнения движения вязкой несжимаемой жидкости В 6 3 были установлены лнфференциальные урзвнения лвижения жидкости в напряжениях. Чтобы написать эти уравнения через проекции вектора скорости, необходимо воспользоваться соотношениями, представляющими компоненты тензора напряжения через компоненты тензора скоростей деформации.