- •«Череповецкий государственный университет» Инженерно-технический институт
- •Технология математического моделирования металлургических процессов курс лекций
- •1. Методология разработки математических моделей металлургических процессов
- •1.1. Роль моделей при разработке и освоении новых металлургических технологий и агрегатов
- •1.2. Методология разработки математических детерминированных моделей металлургических процессов
- •Элементы системного подхода
- •Р Рис. 1.2. Структурная схема процесса разработки модели и моделирования объекта ис. 1.3. Условное представление объекта управления
- •2.2. Примеры проблемных ситуаций
- •2. От объекта до формулировки задач исследования
- •2.1. Основные понятия
- •2.3. Пример выполнения этапа разработки «От объекта до формулировки задач исследования»
- •3.Физическое описание объекта
- •Источники информации и пример физического описания
- •Применение системного подхода к формализации физического описания
- •Р ис.3.2. Структурная схема системы процесса вакуумной обработки
- •4. Математическое описание процесса тепломассопереноса в металлургических агрегатах
- •4.1. Виды балансовых уравнений
- •Некоторые упрощения уравнения теплопереноса
- •Выбор начального и граничных условий
- •4.3.1. Выбор начального условия
- •4.3.2. Выбор граничных условий
- •4.4. Учет выделения тепла кристаллизации при затвердевании двойных сплавов и чистых металлов
- •5. Математическое описание гидродинамических процессов в объеме жидкой стали металлургических агрегатов.
- •Вынужденная конвекция жидкости в ограниченном объеме
- •5.2 Свободная конвекция жидкости в ограниченном объеме
- •6. Разработка численной модели
- •6.1. Метод конечных разностей и пример применения
- •Второй производной по координате
- •6.2. Алгоритмы решения конечно-разностной задачи теплопроводности в цилиндре
- •7. Тестирование компьютерной модели тепловых процессов
- •7.1. Введение в тестирование
- •Тестирование алгоритма решения задачи теплопроводности при граничных условиях третьего рода
- •Вывод критериев подобия для граничных условий третьего рода
- •Тестирование алгоритма решения
- •Тестирование алгоритма решения задачи теплопроводности при теплообмене излучением с окружающей средой
- •Вывод критериев подобия для граничных условий излучением
- •Тестирование алгоритмов решения
- •Дополнительные методы настройки алгоритмов численного решения
- •8. Проверка адекватности и адаптация модели
- •Р ис. 8.1. Кривые охлаждения оси (1) и поверхностного слоя (2):
- •9. Выбор исходных данных для моделирования
- •10. Исследование объекта с помощью модели
- •Р ис. 10.3. Схема к расчету количества вариантов моделирования:
- •11. Решение задач исследования
- •Уравнение теплопереноса
- •Уравнение массопереноса
- •Уравнения гидродинамики
- •Уравнение сплошности
- •Точное решение задачи теплопроводности в охлаждаемой пластине при граничных условиях третьего рода
- •Точное решение задачи теплопроводности в пластине, нагреваемой излучением от окружающей среды
- •Точное решение задачи теплопроводности в охлаждаемом цилиндре при граничных условиях третьего рода
- •Точное решение задачи теплопроводности в цилиндре, нагреваемой излучением от окружающей среды
- •Литература
4. Математическое описание процесса тепломассопереноса в металлургических агрегатах
4.1. Виды балансовых уравнений
В основе детерминированных моделей лежат уравнения, описывающие законы сохранения вещества и энергии, состояния среды и ее движения, равновесия, химические реакции и т.п. Сюда могут входить уравнения алгебраические, дифференциальные, интегральные и т. п., которыми можно описать металлургические процессы. В данном пособии ограничимся рассмотрением только уравнением сохранения энергии и вещества, которые составляют основу всех остальных металлургических процессов. Уравнения сохранения энергии и вещества представляют собой уравнения баланса. Балансы бывают материальные и энергетические, мгновенные и за промежуток времени, локальные и по ограниченному объему пространства.
Рассмотрим примеры балансов:
Материальный баланс углерода в ванне печи за время плавки:
,
(4.1)
где
– начальная и конечная концентрация
углерода (%) в ванне;
– количество углерода, выгоревшего за
время плавки;
- масса i-го компонента шихты с содержанием
углерода Сi;
-
коэффициент использования кислорода;
-
количество кислорода, поступившего в
ванну в % от садки печи.
Материальный баланс жидкости в проточной ванне за единицу времени.
(4.2)
где
и
- количество поступающей и уходящей
жидкости, м3/с;
- скорость изменения объема жидкости
в ванне, м3/с.
Мгновенный баланс тепла в проточной ванне:
(4.3)
где
и
– теплосодержание массы поступающей
и уходящей жидкости в единицу времени,
Вт;
- скорость изменения теплосодержания
жидкости в ванне, Вт;
- теплоемкость жидкости при постоянном
давлении Р, Дж·кг/0С;
–
плотность жидкости, кг/м3;
φ – источник (сток) тепла в ванне, Вт; Т1
и Т2
– температура поступающей и уходящей
жидкости, 0С.
Мгновенный баланс тепла в неподвижном единичном элементарном объеме, характеризуемом точкой с координатами (x, y, z), через который протекает жидкость со скоростью
описывается
уравнением:
(4.4)
где
с и ρ - теплоемкость и плотность материала;
f
– источник или сток тепла в точке
(x,y,z);
- поток тепла к элементарному объему,
;
i,
j,
k
– орты декартовой системы координат;
λ – коэффициент теплопроводности,
– субстанциональная производная;
,
– компоненты вектора скорости
.
Вывод уравнения (4.4) приведен в приложении 1.
Подставляя
выражение для
в уравнение (4.4), при
и используя обозначение
,
получим (4.4) в форме:
(4.5)
Мгновенный баланс примеси в элементарном объеме с координатами (x,y,z) ненеподвижной системе координат, через который протекает жидкость со скоростью описывается уравнением:
(4.6)
где
fc
– источник или сток примеси,
- поток примеси к элементу объема путем
диффузии, С – концентрация примеси в %
от массы объема (может в объемных
процентах и молярных долях),
,
D
– коэффициент диффузии примеси в
жидкости (м2/с),
– субстанциональная производная,
– изменение концентрации в элементе
за счет процесса диффузии,
– конвективные члены, отвечающие за
перенос примеси за счет конвекции
(движения) жидкости. Вывод уравнения
(4.6) приведен в приложении 9.
