
- •«Череповецкий государственный университет» Инженерно-технический институт
- •Технология математического моделирования металлургических процессов курс лекций
- •1. Методология разработки математических моделей металлургических процессов
- •1.1. Роль моделей при разработке и освоении новых металлургических технологий и агрегатов
- •1.2. Методология разработки математических детерминированных моделей металлургических процессов
- •Элементы системного подхода
- •Р Рис. 1.2. Структурная схема процесса разработки модели и моделирования объекта ис. 1.3. Условное представление объекта управления
- •2. От объекта до формулировки задач исследования
- •2.1. Основные понятия
- •2.2. Примеры проблемных ситуаций
- •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. Решение задач исследования
- •Уравнение теплопереноса
- •Уравнение массопереноса
- •Уравнения гидродинамики
- •Уравнение сплошности
- •Точное решение задачи теплопроводности в охлаждаемой пластине при граничных условиях третьего рода
- •Точное решение задачи теплопроводности в пластине, нагреваемой излучением от окружающей среды
- •Точное решение задачи теплопроводности в охлаждаемом цилиндре при граничных условиях третьего рода
- •Точное решение задачи теплопроводности в цилиндре, нагреваемой излучением от окружающей среды
- •Литература
Тестирование алгоритма решения задачи теплопроводности при теплообмене излучением с окружающей средой
Вывод критериев подобия для граничных условий излучением
Выведем критерии подобия для процесса нагрева пластины, охлаждаемой с обеих поверхностей одинаково путем излучения. Температура в начале процесса равномерно распределена по толщине пластины. Теплофизические свойства материала не зависят от температуры. Примем, что начало координат х = 0 находится в плоскости симметрии пластины. С учетом симметрии охлаждения пластины математическое описание процесса нагрева включает
уравнение теплопроводности:
, (7.11)
которое решается в области: , , (7.12)
начальное условие при t = 0: , (7.13)
граничные условия при
х = 0: , (7.14)
x
= S:
,
(7.15)
где
t
– время, Т
– температура, х
– координата по толщине пластины,
отсчитываемая от плоскости симметрии,
– температуропроводность, ,
с,
– теплопроводность, теплоемкость и
плотность материала пластины,
– начальное значение температуры, 2S
– толщина пластины, tк
– конечное время нагрева,
– коэффициент теплоотдачи, л
= 0
- коэффициент излучения серого тела, 0
– коэффициент излучения абсолютно
черного тела,
– степень черноты серого тела (поверхности
пластины), Тср
– температура греющей среды.
Преобразуем
приведенную систему уравнений
(7.11)-(7.15) таким образом, чтобы получить
безразмерные комплексы, включающие
следующие величины: температура, время,
координата и коэффициент излучения.
Для этого все уравнения системы разделим
на величину Тср.
Умножим обе части (11) на
.
Умножим первое неравенство (7.12) на
,
а второе - на
.
Числитель и знаменатель правой части
условия (7.15) умножим на
.
Обе части уравнений (7.14) и (7.15) умножим
на
.
В результате система (7.11)-(7.15) преобразуется
в систему (7.11)-(7.15):
,
(7.11)
, , (7.12)
:
,
(7.13)
,
(7.14)
(7.15)
Обозначим полученные безразмерные комплексы в уравнениях (7.11)-(7.15) следующим образом:
1.
Безразмерная температура -
,
(7.16)
2. Безразмерное время (число Фурье) - , (7.17)
3. Безразмерная координата - , (7.18)
4.
Безразмерная начальная температура -
,
(7.19)
5.
Критерий Старка -
,
(7.20)
Полученные критерий Старка имеет смысл, близкий к критерию Био, и его можно представить, как отношение термического сопротивления тела к термическому сопротивлению границы тела:
,
где
- коэффициент теплоотдачи излучением.
С учетом обозначений система (7.11)-(7.15) примет вид:
(7.11)
(7.12)
(7.13)
(7.14)
(7.15)
Таким образом, решение системы (7.11)-(7.15) зависит от четырех комплексов и имеет вид:
(7.21)
Точные
решения системы (7.11)-(7.15)
представлены в виде таблиц
для различных значений Sk
и
в работе [17] и используются в качестве
тестов.