Страница 43 из 150
В общем случае система располагает n степенями свободы, а интенсиалы изменяются вдоль всех трех координат х , у и z одновременно; такое поле интенсиалов именуется трехмерным. Для вывода простейших уравнений нестационарного переноса используются второе и третье начала ОТ, а также третье частное уравнение пятого начала. В системе мысленно выделяется элементарный объем dV . Количество данного вещества, вошедшего в этот объем за время dt , сопоставляется с количеством вещества, вышедшего из этого объема за то же время. Разница между этими количествами идет на изменение интенсиалов рассматриваемого объема. В результате получается дифференциальное уравнение нестационарного переноса вещества [12, с.303; 14, с.348; 16, с.41; 17, с.104; 18, с.414; 21, с.195]. Здесь для простоты мы ограничимся случаем, когда система располагает всего двумя степенями свободы (n = 2), а ее интенсиалы изменяются только вдоль одной координаты х (одномерное поле интенсиалов). В этих условиях дифференциальное уравнение нестационарного переноса приобретает вид
U1 = L11Z1 + L12Z2 (157)
U2 = L21Z1 + L22Z2
где
U1 = ??P11(?P1/?t) ; U2 = ??P22(?P2/?t) ;
Z1 = ?2P1/?x2 ; Z2 = ?2P2/?x2 ;
?P11 = KP11/m ; ?P22 = KP22/m ;
? - плотность вещества системы, кг/м3; ? - удельная массовая емкость системы по отношению к данному веществу; m - масса системы, кг.
Для гипотетического частного случая, когда n = 1 и поле интенсиала одномерное, находим
U = LZ
или
?P/?t = D(?2P/?x2) (158)
где D - диффузивность:
D = L/(??) (159)
Из выражения (158) в частном случае получаются известные дифференциальные уравнения теплопроводности Фурье, второго закона Фика и т.д. Методы решения дифференциальных уравнений типа (157) разрабатывались Н.А. Буткевичюсом [6] [ТРП, стр.160-161].
14. Особенности применения нестационарного уравнения.
По поводу дифференциального уравнения нестационарного переноса типа (157) требуется сделать несколько замечаний. Прежде всего надо сказать, что границы применимости этого уравнения неодинаковы для различных форм явлений. Эти границы определяются конкретной спецификой явлений и степенью отклонения системы от состояния равновесия.
Если система находится вблизи состояния равновесия, когда перенос осуществляется под действием малых разностей интенсиалов, то уравнение (157) справедливо для любых явлений. С увеличением степени неравновесности результаты рассмотрения отдельных явлений с помощью уравнения (157) заметно искажаются, так как возникают дополнительные степени свободы, начинает заметно сказываться неучтенная специфика распространения и взаимодействия соответствующих веществ и т.д. Например, вблизи равновесия механическая степень свободы, определяемая равенством (43), ничем не осложняется. С увеличением разности давлений появляется скорость перемещения объектов, заметно отличающаяся от нуля, а с нею и новая кинетическая (метрическая) степень свободы. Неучет этой новой степени может привести к существенным ошибкам. Другой пример: при малой скорости жидкость движется ламинарно, при большой движение становится турбулентным, вихревым, то есть появляется дополнительная вращательная степень свободы. Третий пример: распространение электрического заряда вблизи состояния равновесия не влечет за собой никаких неприятностей. С возрастанием разности электрических потенциалов движение заряда сопровождается возникновением кинетической степени свободы и магнитного поля, которыми уже невозможно пренебречь.
В противоположность этому для некоторых других явлений уравнение (157) оказывается справедливым при очень больших отклонениях системы от состояния равновесия. К числу таких явлений относятся вермические (термические), диффузионные и некоторые другие.
Очевидно, что с целью избежания ошибок надо заранее учесть в уравнениях необходимые специфику и дополнительные степени свободы, то есть должны быть заранее выведены более общие и полные уравнения. Тогда при любом отклонении системы от состояния равновесия будут получены правильные результаты. Вблизи состояния равновесия эти общие уравнения должны приводить к более простым частным уравнениям типа (157). Все эти вопросы подробнее затрагиваются при выводе уравнений Максвелла [21] [ТРП, стр.162].
Глава ХII. Шестое начало ОТ.
1. Вывод уравнения.
Согласно пятому началу ОТ, распространение любого данного вещества сопровождается увлечением всех остальных, входящих в переносимый ансамбль. Эффект увлечения одних веществ ансамбля другими определяется перекрестными проводимостями, или коэффициентами увлечения, причем указанному эффекту присущи многие интересные особенности. Чтобы установить эти особенности с количественной стороны, надо вывести соответствующие дифференциальные уравнения. Целесообразно это сделать в самом общем виде, введя группу особых, важнейших для структуры и ее симметрии и, вообще, для термодинамики характеристик А , которые являются функциями главных независимых переменных, входящих в качестве аргументов в основные уравнения ОТ, и измеряются в единицах работы или энергии (в джоулях). Смысл этих характеристик зависит от конкретных значений аргументов и конкретных условий взаимодействия системы и окружающей среды. Для простоты рассуждений ограничимся системой с двумя степенями свободы (n = 2). В термодинамике применительно к термомеханической системе величины А принято именовать характеристическими функциями, или термодинамическими потенциалами.
Выше было показано, что аргументами уравнений могут служить не только экстенсоры, но и интенсиалы. Следовательно, при двух степенях свободы число независимых переменных у каждой из функций должно быть равно двум, а общее число экстенсоров и интенсиалов - четырем. Поэтому количество возможных вариантов аргументов, а значит, и искомых функций А должно соответствовать числу сочетаний из четырех по два, то есть шести. Получается следующий набор аргументов:
(Е1 ; Е2) , (Р1 ; Р2) , (Е1 ; Р2) , (160)
(Е2 ; Р1) , (Е1 ; Р1) , (Е2 ; Р2) .
Третий и четвертый, а также пятый и шестой аргументы дают попарно тождественные результаты, если поменять местами индексы 1 и 2. Применим эти аргументы для определения функций А и вывода на их основе соответствующих законов симметрии структуры.
Нетрудно сообразить, что первый аргумент (?1 ; Е2) приводит к первой характеристической функции ?1 , которая представляет собой не что иное, как энергию U , то есть
А1 = U = F1(?1 ; Е2) Дж (161)
dА1 = dU = Р1d?1 + Р2dЕ2 Дж (162)
Это соответствует прежним уравнениям (30) и (35). Далее автоматически следуют законы структуры (73) и ее симметрии (85) и т.д. Равенство (85) служит исходным звеном в первой цепочке законов симметрии, фактически являющейся следствием применения первого аргумента перечня (160). Кстати, такого типа равенства получили название дифференциальных соотношений, или тождеств, термодинамики, или соотношений Максвелла.
Первое дифференциальное тождество термодинамики (85) мы выводили, когда исходная характеристическая функция ?1 (энергия U ) была уже известна из чисто физических соображений. В отличие от этого при использовании второго аргумента (Р1 ; Р2) нам предстоит найти не только второе тождество, но также и саму исходную функцию А2 . Общий вид второй характеристической функции следующий: