- •Швыдкий Владимир Серафимович дифференциальные уравнения и численные методы
- •Екатеринбург
- •Введение
- •Методы исследования
- •Физический смысл дифференциальных уравнений
- •Глава 1. Дифференциальные уравнения в частных производных
- •Основные положения
- •1.1.1. Характеристики дучп
- •Интерпретация на физической основе
- •1.2. Гиперболические дифференциальные уравнения в частных производных
- •1.3. Параболические дифференциальные уравнения в частных производных
- •Интерпретация с помощью характеристик
- •Интерпретация на физической основе
- •1.4. Эллиптические дифференциальные уравнения в частных производных
- •Глава 2. Математическое описание физических процессов
- •2.1. Общий вид законов сохранения
- •2.2. Характеристическая макроскопическая скорость и диффузионный поток
- •2.3. Баланс массы
- •2.4. Уравнение сохранения количества движения
- •2.5. Уравнение баланса энергии
- •2.6. Обобщённое дифференциальное уравнение
- •2.7. Выбор координат
- •Односторонние и двухсторонние координаты
- •Глава 3. Методы дискретизации
- •Сущность численных методов
- •Р ис. 3.2. Изменение относительной температуры по длине трубы
- •Методы получения дискретных аналогов
- •Метод формальной замены производной конечно-разностными отношениями. Этот метод основан на разложении в ряд Тейлора достаточно гладких функций. Рассмотрим узловые точки, показанные на рис. 3.3.
- •Запишем, например, следующее разложение по формуле Тейлора
- •Анализ методов аппроксимации. Вернёмся к модельной задаче, сформулированной выше, и выявим связь точности численного решения со способом аппроксимации.
- •Основные правила построения дискретных аналогов
- •Заключение
- •Глава 4. Численное моделирование процессов теплопроводности
- •4.1. Стационарная одномерная теплопроводность
- •Нестационарная одномерная теплопроводность
- •Двух- и трёхмерные задачи
- •4.4. Некоторые геометрические соображения
- •4.5. Заключение
- •Глава 5. Конвекция и диффузия
- •5.1. Рассматриваемая задача
- •5.2. Установившиеся одномерные конвекция и диффузия
- •5 .3. Дискретный аналог для двухмерных задач
- •5.4. Дискретный аналог для трёхмерных задач
- •5.5. Односторонняя пространственная координата
- •5.6. Схемная искусственная диффузия
- •5.7. Заключение
- •Глава 6. Расчёт поля течения
- •6.1. Обоснование необходимости в специальной методике
- •6.2. Трудности расчёта поля давления
- •6.3. Переменные "завихренность – функция тока"
- •6.4. Шахматная сетка
- •6.5. Поправки скорости и давления
- •6.6. Уравнение для поправки давления
- •6.7. Модифицированный алгоритм simpler
- •6.8. Заключение
- •Оглавление
- •Глава 1. Дифференциальные уравнения в частных производных...................................12
- •Глава 2. Математическое описание физических процессов............................................28
- •Глава 3. Методы дискретизации.........................................................................................49
- •Глава 4. Численное моделирование процессов теплопроводности.................................72
- •Глава 5. Конвекция и диффузия..........................................................................................109
- •Глава 6. Расчёт поля течения................................................................................................134
- •Библиографический список
2.5. Уравнение баланса энергии
Мы не будем здесь подробно выводить уравнение баланса энергии потока, поскольку это потребовало бы от нас "внедрения" в области физической химии, термодинамики необратимых процессов и т.д. Процедура такого вывода достаточно подробно рассмотрена в нашем учебнике (В.С. Швыдкий, Н.А. Спирин и др. Элементы теории систем и численные методы моделирования процессов тепломассопереноса: Учебник для вузов. – М.: "Интермет Инжиниринг", 1999. – 520 с.), поэтому отметим лишь два момента.
Для
условий работы большинства
теплотехнологических агрегатов градиент
давления в рабочем пространстве
пренебрежимо мал и, следовательно,
термоэффектом, обусловленным давлением,
можно пренебречь. Незначительны также
работа сил давления и аэродинамического
нагрева, перенос теплоты за счёт диффузии
под действием внешних сил (сил тяжести),
а также диффузионный термоэффект. Учёт
этих обстоятельств позволяет существенно
упростить уравнение баланса энергии.
Однако здесь есть и новые аспекты. Прежде
всего, это касается коэффициента
теплопроводности среды ,
который следует заменить его эффективным
значением эф,
включающим турбулентный перенос. Далее
необходимо учесть энергию турбулентных
пульсаций или интенсивность
турбулентности
Ет
= (1/2)
для
которой в турбулентном потоке также
должен выполняться своеобразный закон
сохранения (уравнение переноса), и,
наконец, необходимо определить тепловые
эффекты физико-химических превращений.
Все эти особенности турбулентных потоков излагаются в курсах механики жидкостей и газов и тепломассопереноса. Учитывая, что нашей целью является установление общих свойств различных уравнений, не будем останавливаться на деталях промежуточных преобразований, а сразу запишем в окончательном виде уравнение энергии для двумерного потока, включающее общее уравнение неразрывности:
(2.52)
где
(2.53)
Здесь Hj – тепловой эффект j – той реакции; Rj – эффективная скорость этой реакции; j – коэффициент переноса (диффузии) кинетической энергии турбулентности [ = ( + т)/(Е)]; Е – эмпирический коэффициент модели турбулентности; Qj* - теплота переноса, обусловленная движением среды под действием поля внешних сил; Dэф j – коэффициент диффузии j – того компонента; Jj* - плотность диффузионного потока (в широком смысле слова) этого компонента.
В векторной форме записи уравнение (2.53) принимает вид
(2.54)
В широко распространённой в настоящее время модели турбулентности, включающей два уравнения, в качестве одного из них входит уравнение переноса кинетической энергии турбулентности Ет
(2.55)
где - скорость диссипации (рассеивания) энергии турбулентности, а G – скорость генерации этой энергии. Аналогичное дифференциальное уравнение записывается для переменной .
Нетрудно заметить, что все записанные ранее уравнения "укладываются" в рамки единого представления при некоторых ограничениях. Об этих ограничениях мы поговорим несколько подробнее ниже, а сейчас сформулируем некоторое рабочее правило: всегда при разработке математической модели объекта необходимо выявлять общие тенденции в развитии различных процессов.
