CC BY
34
УДК 535,004.942
К. Ю. Литвинцев
ОСОБЕННОСТИ ИСПОЛЬЗОВАНИЯ КОНЕЧНО-ОБЪЕМНОГО, ДИСКРЕТНО-ОРДИНАТНОГО И ДИФФУЗИОННОГО ПРИБЛИЖЕНИЯ ДЛЯ УРАВНЕНИЯ РАДИАЦИОННОГО ТЕПЛОПЕРЕНОСА
Приведено сравнение методов решения уравнения радиационного теплопереноса: конечно-объемного, диск-ретно-ординатного и диффузионного приближения - на примере решения прикладных и тестовых задач. Рассмотрена целесообразность применения этих методов для различных задач.
Ключевые слова: излучение, КОМ, ДОМ, диффузионное приближение.
Моделирование процесса радиационного теплообмена является очень сложной и ресурсоемкой задачей, так как в этом случае в отличие от остальных процессов теп-ломассопереноса каждый элементарный объем среды находится в непосредственном взаимодействии со всеми другими элементарными объемами и решение интегрально-дифференциальных уравнений, описывающих данное явление, очень трудоемко. Размерность уравнения, описывающего радиационный теплоперенос, может достигать шести измерений: три пространственных координаты х, у, г, две угловых координаты описывающие направление распространения излучения, и частота излучения.
Проблема выбора модели радиационного теплопереноса для расчета прикладных задач заключается в поиске баланса между временем счета и точностью, так как временные затраты на решение только уравнения радиационного теплопереноса (УРТ)
<11 (г,5) . . . . . . . .
—------ = _Р(^)^(^%^) + к(г)/4 (^) +
ст(г)
471
(1)
быть существенно выше Р1 -приближения. Отличие КОМ от ДОМ заключается в том, что при переходе к разностному аналогу (1), в ДОМ происходит интегрирование только по объему:
,еИ‘ к,еИ‘
(X-----+ | ------+ ту------ШУ =
скх с1у ск)
= | (-РI1 + Б1 )ёуап,
(2)
а в КОМ - еще и по угловому пространству [2]:
г1Т^
| |—с1УсЮ.= | | (-р 1,+Б,)ауап, (3)
/7е
А01АУ и* АО1 АV
5' (г,?) = к(г)/Дг,?) + о (г)
4%
(4)
могут существенно превосходить затраты на решение всех остальных процессов. Здесь/-интенсивность, Вт/(ср • м2); /4 - интенсивность абсолютно черного тела; ст - коэффициент рассеивания,^1; к -коэффициентпоглощения,^1; Р - коэффициент затухания (Р = к + ст), м-1; О - телесный угол, рад; г - радиус-вектор, м; 5 - угловое направление, м; Ф - функция рассеивания.
Краткий обзор методов решения уравнения радиационного теплопереноса. Ниже будут рассмотрены три подхода к решению (1): диффузионное (Р1) приближение, конечно-объемный метод (КОМ), дискретно-ординатный метод (ДОМ).
Для решения прикладных задач с радиационным теп-лопереносом наиболее распространено использование Р1-приближения. Оно основывается на предположении от изотропности поля излучения. Основные его плюсы -это простота и относительно малые затраты времени на одну итерацию. Ограничение, накладываемое на применение этого метода, вытекает из предположения, лежащего в его основе, сравнительно слабая анизотропия поля излучения [1].
Методы ДОМ и КОМ являются лучевыми. Они не имеют ограничений по применению, их точность зависит прежде всего от дискретизации углового пространства, однако временные затраты на одну итерацию могут
где м, х, Ь - направляющие косинусы; /-1-е угловое направление для ДОМ или 1-й дискретный телесный угол для КОМ.
Сравнение диффузионного приближения КОМ и ДОМ при решении различных задач. Область применения диффузионного приближения - это оптически толстые среды. Примером такой задачи может служить моделирование газовой топки.
Как правило, в газовых котлах критерий Буггера » 1 и Р1-приближение достаточно хорошо описывает распределение поля излучения. Но в данной статье приведен пример нештатной работы котла, когда через две горелки поступает только воздух без топлива. Эго способствует возникновению сильного градиента температуры в области расположения данных горелок. Моделирование проводилось на основе программы «вРкж», разработанной совместно сотрудниками Красноярского филиала Институт теплофизики Сибирского отделения Российской академии наук, ООО «ТОРИНС» и Сибирского федерального университета. Использование методов Р1, КОМ и ДОМ для решения (1) приводит к формированию отличной друг от друга структуры поля течения (рис. 1). Так, Р1 и КОМ с 8 направлениями имеют структуру поля температур, отличную от остальных (рис. 1). При использовании ДОМ с направлениями поле температур выглядит схожим с полем ДОМ с 24 направлениями и КОМ с 3 2 направлениями, но применение в ДОМ такой малой дискретизации приводит к большой ошибке по тепловым потокам на стенках топки (табл. 1).
Таким образом, использование Р1-приближения в оптически толстых средах дает малое отклонение по
тепловому потоку от лучевых методов (в данном случае ~3 %), но может привести формированию другой структуры поля течения. Что касается лучевых методов, то использование малой угловой дискретизации в этом случае приводит к существенным ошибкам по тепловым потокам на стенки топки для ДОМ и по распределению поля температур для КОМ, поэтому для лучевых методов целесообразнее использовать большее число направлений.
от эталонного решения составила порядка 1 %. но по тепловым потокам на стенку - уже около 13 %. Лучевые же методы показали решение, очень близкое к эталонному. Т=1200 К, е = 0.85
“■ I 1 Г“
0 0.4 0.8
Рис. 2. Идеализированная печь: а - постановка задачи; б - тепловой поток с горячей стенки для разных методов решения УРТ: Р1-приближения: ДОМ с 48 дискретными направлениями ВО Б&, РШМЖ2 - КОМ с 32 дискретными направлениям; НОТТЕЬ - эталонного решения зональным методом [4]
Рис. 1. Распределение поля радиационной температуры в котле (шкала 20-1 800 °С): а - Р1; б - КОМ (8 направлений); в - ДОМ (8 направлений); г - КОМ (32 направления); д - ДОМ (24 направления)
При снижении критерия Буггера меньше единицы в Р1 -приближении, не только падает точность расчета, но и существенно возрастает количество итераций, требующихся для сходимости решения. Это может привести к тому, что основным критерием окончания расчета будет сходимость радиационной задачи, что в целом увеличивает время решения всей задачи радиационного теплопереноса.
В работе [3], где рассмотрено решение для идеализированной печи (рис. 2), критерий Буггера равен 0,8. Результаты расчетов сравнивались с эталонным решением задачи методом Хоттеля. Ошибка по температуре для Р1
Лучевые методы по сравнению с диффузионным приближением сходятся быстрее, особенно в оптически тонких средах, но время, требуемое ими на одну итерацию возрастает примерно в;? раз, где п - число угловых направлений. Сократить временные затраты можно путем замораживания радиационных источниковых членов в уравнении энергии в течение нескольких итераций Один из недостатков лучевых методов (в большей степени КОМ) связан с их использованием для геометрически сложных объектов. В этом случае может наблюдаться угловое перекрытие, т. е. пересечение грани контрольного объема и дискретного телесного угла (рис. 3), и при недостаточном угловом разбиении в сильноанизотроп-ных средах возникает серьезная ошибка. Для разрешения этой проблемы в КОМ можно использовать дополнительное разбиение телесного угла на два угла плоскостью грани, пересекающей телесный угол.
Таблица 1
Сравнение методов по температуре и величине теплового потока на стенки котла
Мегод решения уравнения радиационного теплопереноса Величина тепловой потока, Вт/м2
КОМ (8 направлений) 1,32Е + 08
КОМ (32 направления) 1,30Е + 08
Р1 1,26Е.+ 08
ДОМ (8 направлений) 1,55В+ 08
ДОМ (24 направления) 1,29Е + 08
Представим расчет поля радиационной температуры в цилиндре радиусом 0,5 м и высотой 4 м, к = ОД м1 (рис. 4), стенки которого имеют различные температуру и степень черноты. При малом количестве угловых направлений наблюдается сильный дисбаланс по энергии, если не учитывать угловое разбиение (табл. 2). В этой задаче среда является оптически тонкой и Р1 -приближение показывает сильное расхождение с лучевыми методами.
Рис. 3. Угловое перекрытие
х < 0,5 м коэффициент рассеивания б принимается равным 2 м1, в остальной - равным нулю. На нижней стенке задается поток излучения д. Остальные стенки считаются абсолютно черными и неизлучающими. Размерность пространственной сетки составляет 22 х 22 х 22 узла.
Рис. 5. Схематическое отображение угловой многоблочно сти 0,5 м і 0,5 м <------------->\<-------->
а = 2 м
а = 0 м
і
а б в г д
Рис. 4. Сравнение методов на сетке с неортогональными ячейками (для Р1 шкала 600...630 °С, для лучевых -500...670 °С): а -сетка; б-Р1-приближение; в-ДОМ (48 направлений); г - КОМ (8 направлений с учетом углового перекрытия); д - КОМ (288 направлений без учета углового перекрытия)
Использование угловой многоблочности в КОМ. Ко-
нечно-объемный метод обладает дополнительным инструментом, позволяющим сократить время расчета. Этот инструмент-угловая многоблочность [4], когда различные расчетные области задачи могут иметь различную угловую дискретизацию (рис. 5).
Такой подход позволяет сокращать время расчета, если примерно известны оптические свойства среды. Приведем пример использования угловой многоблочности (рис. 6). Расчетная область имеет кубическую форму. Среда принимается неизлучающей, т. е. к = 0. Во всей области
X Поток излучения q Рис. 6. Использование угловой многоблочности
Представим результаты расчетов трех вариантов (рис. 7): первый вариант - 288 телесных углов, второй - 32 телесных угла, третий - 32 телесных угла в области а = 2 м1 и 288 телесных углов, в области а = 0 м1. При использовании малого количества угловых разбиений (2-й вариант) наблюдался нефизичный провал. Для получения точного решения было существенно увеличено разбиение углового пространства (1-й вариант). Это привело к существенному росту времени счета (табл. 3). Использование угловой многоблочности (3-й вариант) позволило получить решение, достаточно близкое к 1-му варианту, при этом общее время счета уменьшилось в три раза. Возросло только количество итераций, что связано с дополнительной обработкой границы угловых блоков.
Интенсивный рост производительности вычислительных систем делает использование лучевых методов доступным для решения прикладных задач, но в геометрически сложных объектах в условиях сильной анизотропии излучения требуется достаточно сильная дискретизация этих методов по угловому пространству, что может привести к серьезному увеличению ресурсоемкости задачи. Диффузионное приближение в области прикладного моделирования до сих пор остается более привлекательным по соотношению времени счета к точности ре-
Таблица 2
Пример влияния углового перекрытия
Угловое разбиение (кол-во углов) Потери энергии за счет схемы (\Епад - ЕИСХ\/ Ев сх), %
КОМ без учета углового перекрытия КОМ с учетом углового перекрытия ДОМ
8 207 12 00
24 57 5 77
48 10 - 4
288 4 - -
шения, однаш еш использование в оптически тонких сре- Библиографический список
дах и при наличии сильных градиентов поля температуры приводит, к существенным ошибкам по распределе- * • Четверушкин, Б. Н. Математическое моделированию поля излучения и радиационным потокам на повер- ние задач динамики излучающего газа / Б. Н. Черверуш-
кин. М.: Наука, 1988.
2. Chai, J. С. Finitc-volumc method for radiation heat transfer, to Appear in Advances in Numerical Heat Transfer / J. C. Chai, S. V Patankar. Editors: Minkowycz and Sparrow, Vol. 2.2000.
3. Файвлэнд, В. А. Решение уравнения радиационного тепло переноса методом дискретных ординат /В. А. Файвлэнд//Аэрокосмическая техника. 1989. №9.
4. Chai, J. С. Angular-multiblock procedure for radiation heat transfer / J. C. Chai, J. P. Moder // Intern. Conf. on Computational Heat and Mass Transfer. Famagusta, Cyprus, 1999.
вдоль линии г = 1 м,у = 0,5 м
Таблица 3
Сравнение сходимости и времени расчета для вариантов без использования угловой многоблочности и с ее использованием
Угловое разбиение Число итераций Общее время счета, с Затраты времени на одну итерацию, с
32 6 6 1
288 6 219 36,5
32 + 288 7 73 10,4
K. Yu. Litvitsev
FEATURES OF THE FINITE - VOLUME METHOD, DISCRETE - ORDINATES METHOD AND PI APPROXIMATION APPLICATION FOR RADIATION HEAT TRANSFER EQUATION
The article presents a comparison of the finite - volume method and PI approximation for radiation heat transfer equation. The expediency of these methods application for various problems are considered.
Keywords: radiation, FVM, DOM, PI approximation.
УДК519.68
Е.С. Семенкин, К. А. Токмин
КОЭВОЛЮЦИОННЫЙ ПОДХОД ДЛЯ РЕШЕНИЯ ЗАДАЧ УСЛОВНОЙ МНОГОКРИТЕРИАЛЬНОЙ ОПТИМИЗАЦИИ
Предложен коэволюционный подход к решению задач многокритериальной оптимизации с ограничениями на основе метода обобщенного лагранжиана и эволюционных алгоритмов. Проведено исследование эффективности разработанного алгоритма на тестовых и практических задачах.
Ключевые слова: задача многокритериальной условной оптимизации, коэволюционный алгоритм.
Практические задачи оптимизации не всегда могут быть удобных для оптимизации свойств или, по крайней мере,
эффективно решены средствами классической теории оп- информации о таких свойствах и т. д. К тому ж характерной
тимизации [ 1; 2], так как эти задачи обладают свойствами, чертой, объединяющей классические подходы к решению
существенно затрудняющими их решение: дискретными или многокритериальных задач, является то, что в кавдом из них
смешанными переменными, алгоритмически заданными задача многокритериальной оптимизации сводится к одной
целевыми функциями и (или) ограничениями, отсутствием или нескольким задачам однокритериальной оптимизации.
хности исследуемого объекта.
оооооооооо X. М
Ф 1-ый вариант А, 2-ой вариант-^^З-ой вариант
Рис. 7. Распределение падающего теплового потока
