МЕТОДИЧЕСКОЕ И ПРОГРАММНОЕ ОБЕСПЕЧЕНИЕ ОПРЕДЕЛЕНИЯ ТЕПЛОПРОВОДНОСТИ ГРАНУЛИРОВАННОГО МАТЕРИАЛА
(Работа выполнена при поддержке гранта РФФИ № 08-08-00980) П.А. Трубаев, д.т.н.
(Белгородский государственный технологический университет им. В.Г. Шухова, trubaev@mail.ru)
Ключевые слова: теплообмен, метод определения теплопроводности, программное обеспечение численного моделирования температурного поля гранулы.
Процесс теплопроводности в обрабатываемом материале оказывает значительное влияние на теплообмен в теплотехнологических установках производства стройматериалов, в частности, обжига цементного клинкера. Поэтому моделирование их тепловой работы требует использования отсутствующих зависимостей коэффициента теплопроводности силикатных материалов от их состава и температуры. Стационарные методы исследования теплопроводности разработаны наиболее полно, но они характеризуются низкими рабочими температурами, сложностью экспериментальных установок и значительным временем проведения эксперимента [1]. Нестационарные методы и методы неравновесной термодинамики основаны на исследовании меняющихся во времени температурных полей [2]. Основная трудность их реализации состоит в приближении используемой математической модели к реальным экспериментальным условиям.
В работе предложен метод определения зависимости коэффициента теплопроводности гранулированного материала от температуры, основанный на решении нестационарной обратной задачи теплопроводности. Принцип метода заключается в экспериментальном определении температур в двух точках гранулы с диаметром 5-25 мм при ее нагреве или охлаждении и моделировании температурного поля гранулы численным решением дифференциального уравнения нестационарной теплопроводности в сферических координатах [3] с приближением модели к экспериментальным данным. По изменению поля температуры производится расчет теплового потока, который равен изменению теплосодержания гранулы за определенный промежуток времени, и расчет коэффициента теплопроводности исследуемого материала.
Необходимо отметить отличие предлагаемого метода от классической обратной задачи. Приведение в соответствие результатов численного моделирования и экспериментальных данных путем изменения коэффициента теплопроводности для рассматриваемой задачи невозможно, так как найденные коэффициенты будут содержать значительную погрешность, связанную с неточностью определения параметров модели (радиуса и условий теплоотдачи на поверхности) из-за неправильной формы гранулы. В предлагаемой методике моделирование служит для точного расчета количества теплоты, теряемой гранулой, по которо-
му и определяется коэффициент теплопроводности.
Исследование теплопроводности по предлагаемой методике состоит из следующих этапов.
1. Исследование материала гранулы, определение истинной и кажущейся плотности, эквивалентного диаметра.
2. Экспериментальное определение температур в двух точках гранулы 10 и измеряемых с помощью двух термопар (рис. 1), при ее нагреве в муфельной печи до 1 000 °С и при охлаждении в естественных условиях от 1 000 до 200.. .300 °С.
Внутренняя часть
Термопары
ч \ э \ . -\ \
У ✓
Гк Я
Рис. 1. Закладка термопар при экспериментальном измерении температуры гранулы
3. Приближение температур, определенных в эксперименте 10 и и рассчитанных по модели теплопроводности 1М и 1М (рис. 2), выполняемое путем изменения граничных условий модели и других параметров, измерение которых связано с погрешностью.
4. Расчет коэффициента теплопроводности для каждого момента замера температуры, выполняемый по тепловому потоку Q, проходящему через сферу радиусом гк, ограничивающую внутреннюю часть гранулы, площади поверхности сферы и по разнице температуры в малом объеме у ее поверхности (рис. 2).
5. Определение по полученной табличной зависимости коэффициента теплопроводности X от температуры 1 в виде регрессионного уравнения Х=Г(1) и итерационное уточнение зависимости.
Согласно изложенной методике разработано программное обеспечение для численного моделирования температурного поля гранулы, обработки результатов эксперимента и определения
1
¿к 1 |\ 11 \ 11 \ ■ 1 \ 11 \ 11 \
/ / / ' м [ .' 'о \ 1 1 хА \ N 1 » I 1 1
1 к \ < 0 V* Ш >1
а б
Рис. 2. Расчетное температурное поле гранулы: а - нагрев; б - охлаждение
коэффициентов зависимости коэффициента теплопроводности от температуры (рис. 3).
Алгоритм работы программы следующий.
1. Считывание из файла экспериментальных данных: радиуса гранулы И, плотности гранулы р, времени от начала нагрева (охлаждения) до начала замеров Дт0, температур 10 т и ^ т, т=1...№э, где N'3 - число всех экспериментальных точек (как правило, 50-80).
2. Удаление экспериментальных точек, не входящих в заданный диапазон температур, сглаживание экспериментальных данных полиномом заданной степени. Результат - температуры 10 m и 1к т, т=1...№,, где N.3 - число экспериментальных точек.
3. Численное моделирование температурного поля гранулы 1" ^ !=0...^ ]=0..^м, где N - максимальный номер точки расчетной сетки по радиусу гранулы; N - число итераций по времени, определяемое из условия заданного шага по времени Дтм (как правило, Дтм=0,1 сек.).
4. Минимизация квадратичного отклонения 8 экспериментальных и расчетных температур путем изменения заданных параметров методом покоординатного спуска с дроблением шага (итерационный повтор п. 3).
5. Расчет коэффициента теплопроводности ЯшОт), т=1...№,, и коэффициентов регрессионного уравнения А=(1), где Щ) - квадратичное а0+ +а11+а212 или кубическое а^а^+а^+а^3 уравнение с оцениваемыми коэффициентами а0,а1;а2, а3.
6. Итерационный расчет пп. 3-5 до стабилизации коэффициентов регрессионного уравнения
В программе реализованы следующие возможности:
- выбор количества точек сетки N и шага по времени в модели Дтм;
- выбор вида регрессионного уравнения А=(1), используемого в модели теплопроводности,
которым может быть полином 2-й или 3-й степени;
- выбор наличия или отсутствия сглаживания исходных экспериментальных данных и степени сглаживающего полинома (от 2-й до 4-й);
- задание диапазона температур экспериментальных данных для исключения вносящих заметную погрешность начального участка, на который оказывает влияние процесс извлечения гранулы из муфельной печи или помещения гранулы в нее, и конечного участка, когда температуры мало изменяются по времени и по радиусу;
- выбор схемы дискретизации: Кранка-Никол-сона, неявной и совмещенной, предложенной в работе [3];
- выбор способа приближения расчетных и экспериментальных данных путем изменения следующих параметров: коэффициента ка для корректировки коэффициента теплоотдачи на поверхности; времени до начала замеров Дт0; начальной температуры гранулы при охлаждении или температуры в муфельной печи при нагреве; диаметра И. Для перечисленных способов возможно установление их активности или неактивности и очередности применения, которая задается перемещением способа в списке на экране;
- выбор между использованием для вычисления коэффициента теплопроводности сглаженных экспериментальных значений температур или температурного поля модели;
- расчет температурного напора Д1 по двум соседним узлам сетки с номерами к и к-1 или по двум точкам, с номером к и лежащей рядом с ней на заданном малом расстоянии Дг<<гк-гк-1, температура во второй точке при этом определяется аппроксимацией полиномом Лагранжа температурного поля т;
- задание начальных значений для коэффициентов, используемых при приближении;
- наглядное представление хода расчета в виде демонстрации в реальном времени в графическом виде экспериментальных точек, температур, полученных в результате моделирования, полу-
ченной зависимости для коэффициента теплопроводности.
Графический интерфейс программы позволяет визуально оценивать точность аппроксимации экспериментальных данных полиномом, степень приближения температур модели и эксперимента, вид получаемой зависимости для коэффициента теплопроводности, дает возможность изменением исходных данных подбирать наиболее удачный вариант расчета. Использование графического диалогового интерфейса и наглядного представления хода и результатов расчета оказало также большую помощь при разработке и отладке алгоритма программы.
Для тестирования программы в качестве исходных данных задавалось распределение температур, полученное аналитическим решением уравнения нестационарной теплопроводности при постоянных теплофизических параметрах. Результаты расчета коэффициента теплопроводности при различных вариантах параметров процесса (плотности, теплоемкости, коэффициента теплоотдачи) не отличались от задаваемого при аналитическом решении более чем на 0,5 %.
Предложенный метод использован для определения теплопроводности заводских гранул АО
«Белгородский цемент», АО «Осколцемент» и гранул, полученных в лабораторных условиях обжигом в высокотемпературной печи. Статистическая проверка показала однородность полученных данных с результатами измерения теплопроводности в установках стационарного плоского и цилиндрического слоя, в которых свойства определялись при температурах 200-400 °С [4].
Можно заключить, что разработанный метод позволяет с достаточной степенью точности определять коэффициент теплопроводности гранул при температурах до 1 000 °С и их зависимость от температуры, не требуя дорогостоящего специализированного оборудования.
Литература
1. Осипова В.А. Экспериментальное исследование процессов теплообмена. - М.: Энергия, 1979. - 320 с.
2. Бек Дж., Блакуэл Б., Сент-Клэр Ч. Некорректные обратные задачи теплопроводности. / Пер. с англ. - М.: Мир, 1989. -312 с.
3. Трубаев П.А., Беседин П.В. Численное моделирование процесса охлаждения клинкерных гранул в колосниковом холодильнике. // Изв. вузов. Строительство. - 2004. - № 6. -С.120-125.
4. Беседин П.В., Трубаев П.А. Исследование и оптимизация процессов в технологии цементного клинкера. - Белгород: Изд-во БГТУ им. В.Г. Шухова, БИЭИ, 2004. - 420 с.
ПЕРСПЕКТИВЫ ИСПОЛЬЗОВАНИЯ МЕТОДОВ ДИНАМИЧЕСКОЙ ОПТИМИЗАЦИИ ПОРТФЕЛЯ АКТИВОВ И ПАССИВОВ КОММЕРЧЕСКОГО БАНКА
Т.М. Вандышева
(Технологический институт Южного федерального университета, г. Таганрог, wandischewatm@list.ru)
Ключевые слова: условия обеспечения устойчивого развития кредитной организации, методы оптимизации, управление активами и пассивами, многоэтапное стохастическое программирование.
Повышение значимости банковской системы в российской экономике позволяет определить как одну из приоритетных задачу повышения устойчивости и эффективности функционирования кредитных организаций. Данное положение было закреплено в Стратегии развития банковского сектора Российской Федерации на период до 2008 г. в качестве основной цели на среднесрочную перспективу.
Одним из основных условий обеспечения устойчивого развития кредитной организации является эффективное управление активами и пассивами предприятия, которое возможно с помощью методов и моделей портфельной оптимизации. Данный инструментарий позволяет решить задачу формирования эффективного портфеля кредитной организации посредством определения таких пропорций распределения средств между допустимыми активами, которые обеспечат минимальный
риск портфеля при определенном (заданном) уровне доходности.
Исторически первым математическую модель формирования оптимального портфеля в статическом случае как задачу исследования операций и теории игр сформулировал Г. Марковиц [1]. Публикации ученого вызвали большой поток научных работ по финансовой теории и теории исследования операций (например [2,3]).
Интерес к использованию методов оптимизации при формировании эффективного портфеля огромен, при этом наиболее перспективными представляются исследования, связанные с динамическими моделями управления активами и пассивами. Целью данной публикации является аналитический обзор проблем и перспектив практического использования моделей данного вида в деятельности предприятий банковского сектора.