ПРИКЛАДНАЯ МАТЕМАТИКА И МЕТОДЫ МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ
УДК 544.45+519.2
ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ ТЕПЛОВОЙ СТАБИЛЬНОСТИ ГРАНУЛЫ КАТАЛИЗАТОРА С ВНУТРЕННИМ ТЕПЛОВЫДЕЛЕНИЕМ В СЛУЧАЙНОМ ПОЛЕ ТЕМПЕРАТУРЫ СРЕДЫ
И.В. Деревич, Д.Д. Галдина
МГТУ им. Н.Э. Баумана, Москва, Российская Федерация e-mail: [email protected]; [email protected]
Предложен метод имитационного моделирования температуры гранулы с внутренним тепловыделением, погруженной в среду со случайными флуктуациями температуры. Метод основан на решении системы обыкновенных дифференциальных стохастических уравнений, представляющей флуктуации температуры среды случайным процессом с конечным временем затухания автокорреляционной функции. Проведено тестирование метода на основе сопоставления с точными аналитическими результатами. Показано, что поведение температуры гранулы с внутренним тепловыделением в случайной среде качественно отличается от результатов, получаемых в рамках детерминированного подхода. Рассчитаны времена ожидания теплового взрыва в среде с флуктуациями температуры при различных режимных параметрах.
Ключевые слова: стохастическое обыкновенное дифференциальное уравнение, автокорреляционная функция, тепловой взрыв, диаграмма Семенова, флуктуации температуры.
NUMERICAL STUDY OF THERMAL STABILITY OF A CATALYST GRANULE WITH INTERNAL HEAT RELEASE IN THE RANDOM FIELD OF MEDIUM TEMPERATURE
I.V. Derevich, D.D. Galdina
Bauman Moscow State Technical University, Moscow, Russian Federation e-mail: [email protected]; [email protected]
A method is proposed for simulating the temperature of a granule with internal heat release that is immersed in a medium with random temperature fluctuations. The method is based on solving the system of stochastic ordinary differential equations that presents the medium temperature fluctuations as a random process with the finite time of the autocorrelation function decay. The method testing is conducted on the basis of comparing the test data with exact analytical results. It is shown that the behavior of temperature of the granule with internal heat release in a random medium differs qualitatively from the results obtained within the deterministic approach. The time that elapses before the heat explosion in a medium with temperature fluctuations is calculated for different regime parameters.
Keywords: stochastic ordinary differential equation, autocorrelation function, heat explosion, Semenov’s diagram, temperature fluctuations.
Введение. Постановка задачи. Как правило, каталитические процессы синтеза сопровождаются выделением теплоты. Практический
ISSN 1812-3368. Вестник МГТУ им. Н.Э. Баумана. Сер. “Естественные науки”. 2014. №2
3
пример — синтез тяжелых углеводородов в процессе синтеза Фишера-Тропша [1]. Синтез искусственной нефти позволяет решить экологические проблемы, связанные с утилизацией попутного жирного газа на нефтяных месторождениях, вовлечь в экономическую деятельность труднодоступные месторождения природного газа, обеспечить экономическую независимость стран от импорта углеводородного сырья.
В Технологическом институте сверхтвердых и новых углеродных материалов (ФГБНУ ТИСНУМ, Троицк) разрабатывается промышленный реактор на 5 000 м3/ч синтез-газа с проектной мощностью до 500кг/ч стабилизированных жидких углеводородов. В реакторе с неподвижным слоем гранулированного катализатора синтез происходит внутри пористых гранул и сопровождается выделением большого количества теплоты. Теплота экзотермической реакции отводится из объема гранулы к границе, с которой отводится за счет теплоотдачи в жидкие продукты синтеза. Превышение мощности тепловыделения над мощностью теплоотвода с поверхности приводит к неконтролируемому увеличению температуры — тепловому взрыву. Потеря тепловой стабильности гранул катализатора может вызвать тепловой взрыв всего реактора. В связи с этим исследование границ и условий начала теплового взрыва — актуальная задача, определяющая стабильные технологические режимы работы установок синтеза.
Условия, приводящие к потере тепловой стабильности в детерминированных условиях, достаточно хорошо изучены [1-6]. Существует критическая температура, превышение которой вызывает существенное увеличение температуры гранулы. Ситуация качественно меняется, когда температура среды является случайной величиной. В этом случае всегда с ненулевой вероятностью существует флуктуация температуры, значение которой превысит критическое значение, что может вызвать потерю тепловой стабильности. Исследованию влияния случайного шума на поведение систем с резким изменением состояний посвящены, например, работы [7-11]. Результаты настоящего исследования также могут быть использованы при расчетах условий воспламенения дисперсного топлива в авиационных и ракетных двигателях, в энергетике, для оценок вероятности возникновения взрыва при хранении и транспортировке дисперсных горючих материалов.
Влияние шума на взрывное поведение систем рассмотрено в работах [9-16] в рамках аппарата функции плотности вероятности. Это требует использования современных методов случайных процессов и функционального анализа и позволяет получить результаты, имеющие практическое значение. Однако метод функции плотности вероятности абстрагируется от некоторых важных деталей химической кинетики,
4
ISSN 1812-3368. Вестник МГТУ им. Н.Э. Баумана. Сер. “Естественные науки”. 2014. № 2
например, сложно включить детальную кинетическую схему процесса синтеза, а также внутреннюю структуру каталитической гранулы.
В этой ситуации следует использовать методы прямого численного моделирования динамики изменения температуры на основе решения стохастических обыкновенных дифференциальных уравнений [1-19]. Классический подход, развиваемый в методах имитационного моделирования, основан на уравнении Ланжевена [20], где случайным источником является процесс Гаусса с дельта-коррелированной во времени автокорреляционной функцией (белый шум).
Известно, что белый шум — абстракция, энергия этого процесса бесконечна и процесс недифференцируем [21]. Для генерации случайного процесса с конечным временем затухания автокорреляционной функции используется прием, основанный на решении системы стохастических дифференциальных уравнений [17-19]. С помощью такой технологии можно моделировать случайные процессы с широкой гаммой автокорреляционных функций.
В настоящей работе предложен метод численного моделирования случайной температуры гранулы с внутренним тепловыделением в среде с флуктуациями температуры, имеющими конечное время затухания автокорреляционной функции. Этот метод позволяет изучить особенности процесса потери тепловой устойчивости и может быть применен в дальнейшем для моделирования стохастического поведения как температуры, так и концентраций реагентов внутри гранулы с учетом детальной кинетики синтеза. Проведено тестирование предложенного метода на примере динамики случайных систем, осредненные параметры которых имеют аналитическое описание. Проиллюстрированы различные сценарии поведения температуры частицы в стохастической среде. Представлены результаты расчетов среднего времени ожидания теплового взрыва. Сформулированы задачи дальнейших исследований.
Уравнения для температуры частицы с внутренним тепловыделением. Диаграмма Семенова. Рассмотрим сферическую частицу радиусом dp, м, находящуюся в жидкой среде с температурой 0/, K. Внутри частицы протекают экзотермические химические реакции с тепловым эффектом Q, Дж/моль. Скорость химической реакции задаем по закону Аррениуса с энергией активации E, Дж/моль. Теплота с поверхности частицы отводится в жидкость за счет теплоотдачи, по закону Ньютона, коэффициент теплоотдачи а, Вт/(м2-К). Уравнение для средней температуры частицы 0Р имеет вид
d0 _ E
тРсР-00~ = aSP (0/ - 0Р) + VpQAe RO&p ,
ISSN 1812-3368. Вестник МГТУ им. Н.Э. Баумана. Сер. “Естественные науки”. 2014. №2
5
где mp — масса частицы, кг; Sp = ndp — площадь поверхности частицы, м2; Vp = nd^/6 — объем сферической частицы, м3; A — предэкспоненциальный множитель, моль/(м^с); R° — универсальная газовая постоянная, Дж/(моль^К); cp — теплоемкость материала частицы.
Уравнение для температуры частицы можно переписать в релаксационном виде
d0p = 0/ 0p + ^Ае-R°0V
dt тв PpCp
(1)
Здесь то = dpcp/ap — время тепловой релаксации частицы, т.е. характерное время, за которое температура частицы достигает постоянной температуры жидкой среды.
Уравнение (1) приводим к безразмерному виду. Выбираем соответствующие масштабы. Поскольку далее исследуется влияние флуктуаций температуры несущей фазы, температура жидкой фазы задается как
0/ (t) = (0/ > + df (t),
где осредненная температура (0/> не зависит от времени и однородна в пространстве; угловые скобки обозначают осреднение по множеству случайных “траекторий” в пространстве температур — осреднение по ансамблю. Флуктуации температуры среды в/ заданы так, что (в/> = 0. В качестве временного масштаба выбран интегральный временной масштаб Те затухания автокорреляционной функции температуры среды. Определение этого масштаба будет приведено далее. Масштаб температуры — осредненная температура несущей среды. Для выбранных масштабов уравнение (1) принимает вид
dr
(1 + 7/) - Г
E*
+ Q*e Гр ,
(2)
dT Q0
где т = t/TE — безразмерное время; rp = 0p/ (0/> — безразмерная температура частицы; Q0 = т0/Те — параметр тепловой инерции частицы; у/ = в// (0/> — безразмерные флуктуации температуры несущей среды; E* = E/R° (0/> — безразмерная энергия активации; Q* = QA/ (ppcp (0/>) — безразмерная мощность тепловыделений.
Согласно уравнению (2), температура частицы изменяется в результате двух механизмов: 1) теплоотдачи с поверхности (первое слагаемое в правой части уравнения); 2) тепловыделения (второе слагаемое в правой части уравнения).
На основе анализа диаграммы Семенова покажем существование критической температуры, бесконечно малое превышение, которой приводит к неконтролируемому увеличению температуры частицы —
6
ISSN 1812-3368. Вестник МГТУ им. Н.Э. Баумана. Сер. “Естественные науки”. 2014. № 2
тепловому взрыву. Анализ диаграммы Семенова выполним для стационарной температуры несущей среды yf = 0. Находим стационарные температуры частицы
= 0 Гр 1
dr = 0 О©
Q* e гр .
Вводим безразмерные мощность теплоотдачи W* = (Гр — 1) /О© и мощность тепловыделения WQ = Q* exp (—E*/Гр). На рис. 1 представлены функции, описывающие безразмерные тепловыделения и теплоотдачу. Видно, что существует область с тремя стационарными температурами частицы, ограниченная касательными, положение которых определяется значениями параметра тепловой релаксации частицы.
На касательных c и b (см. рис. 1) гранула возвращается в стационарное состояние с низкой или высокой температурой соответственно. В целях изучения типа стационарных температур было проведено численное интегрирование нелинейного уравнения (2) без учета флуктуаций температуры среды. Динамику изменения температуры гранулы, если начальная температура близка к значению второго корня на диаграмме Семенова (см. рис. 1), иллюстрирует рис. 2. Видно, что бесконечно малое превышение температуры ГII приводит к потере тепловой стабильности, температура гранулы стремится к высокой температуре Гщ. Если начальная температура меньше значения Гц на малую величину, то температура гранулы переходит в область низких температур Гь близких к температуре среды. С позиции прогнозирования начала теплового взрыва температуру ГII можно полагать критической:
Гсг = ГН.
Рис. 1. Диаграмма Семенова
Рис. 2. Диаграмма изменения температуры гранулы с начальной температурой, близкой к значению второго корня, взятого с диаграммы Семенова:
1 - Г = 1,02885; 2 - Гп = 2,39015; 3 -Гш = 18,98075
ISSN 1812-3368. Вестник МГТУ им. Н.Э. Баумана. Сер. “Естественные науки”. 2014. №2
7
Автокорреляционные функции флуктуаций температуры. Точные результаты. В этом разделе будут получены точные результаты, с которыми сопоставим данные прямого численного моделирования стохастического процесса флуктуаций температуры частицы. Точные решения существуют для линейных уравнений. Рассмотрим уравнение для флуктуаций температуры гранулы (1) без источника химической реакции
d°P (t) = °f(t)- °Р (t) (3)
dt t©
Здесь предполагается, что средние значения флуктуаций температур гранулы и среды (dp (t)) = {Of (t)) = 0. Уравнение (3) — уравнение релаксационного типа. Для t ^ т© начальное значение температуры гранулы “забудется” и не будет вносить вклад в динамику изменения температуры. Флуктуации температуры гранулы будут определяться только параметрами среды, поэтому для сокращения записи начальное значение температуры гранулы примем равным нулю. Случайный процесс Of (t) — статистически стационарный случайный процесс Гаусса с корреляционной функцией
(Of (to Of (t")) = (02) ¥f (t - t"). (4)
Используя метод спектрального анализа случайных процессов [20, 21], находим из уравнений (3) и (4) корреляции и интенсивность флуктуаций температуры гранулы для различных видов автокорреляционной функции флуктуаций температуры среды Фf (t).
Случайные флуктуации температуры среды представим в виде спектрального разложения по частотам
Of (t)
1
2п
e-iwtda (ш),
(5)
где da (ш) — случайная мера в пространстве частот.
Подставив выражение для спектрального разложения (5) в формулу
(4), получим
(Of Ю Of Ю)
1
(2Л7
e«yt'W't" (da (ш') da (ш''))
(6)
Из условия статистической стационарности случайного процесса
(4) вытекает функциональный вид коррелятора случайной меры в пространстве частот:
(da (ш') da (ш'')) = 2п (Of) ^f (ш') 8 (ш' + ш"). (7)
Здесь 8 (ш) — дельта-функция Дирака; ^f (ш) — спектр случайного процесса Of (t).
8
ISSN 1812-3368. Вестник МГТУ им. Н.Э. Баумана. Сер. “Естественные науки”. 2014. № 2
Подставив формулу (7) в выражение (6), запишем связь между автокорреляционной функцией и ее спектром
(в, («') в, (t")> = (в,) ф, (f -t'') =
в2.
2п
<1ш''ф, (ш'') dJеш t+^tS (ш' + ш'') =
2п
ц, (j) dj.
Ф, (t) = 2n e~-utф, (ш) dj.
(8)
При t' = t'' из формулы (8) определим условие нормировки спектра
2п Iф,(ш) dj = 1.
(9)
Процедура нахождения спектра по автокорреляционной функции следует из формулы (8):
Ф, (ш) = егш8Ф , (s) ds.
(10)
-ГО
На основе спектрального представления (5) получим условие дифференцируемости случайного процесса: процесс дифференцируем, если существует величина ((dO,/dt)2) [21]. Рассмотрим более общую корреляцию
dO, (t') doO, (t'') dt dt
1
(2П)2
f ff iw't' t j ( f\ j ( H\\
ш ш e + (da (ш ) da (ш )).
Для раскрытия корреляции случайных мер используем функциональный вид (7), тогда
dO, (t') doO, (t'') dt dt
— f dш"ш"ф, (ш'') f dш'ш'еь
8(ш' + ш'') =
(O2\ р
= , e-iw(t'-t")ш2ф, (ш) с1ш. (11)
2п J
— га
ISSN 1812-3368. Вестник МГТУ им. Н.Э.Баумана. Сер. “Естественные науки”. 2014. №2
9
В соответствии с формулой (11) существование производной случайного процесса связано со сходимостью интеграла
ОО
и)2ф/ (w) dw.
-га
Из формул (8) и (11) также следует представление для среднего квадрата производной случайного процесса через автокорреляционную функцию
d2*/ (t)
dt2
t=о
Решение уравнения (3) с нулевым начальным условием имеет вид
t
1 f _t—
вр (t) = — e т- 9/ (s) ds. T© J о
Двухвременная корреляция флуктуаций температуры гранулы равна
t' t"
1 f t'-s' 1 Г t''-s''
(9p(t')9p(t")) = - ds'e- r- _ ds"e~ - (9/(s')9/(s'')). (12)
t©J t© J
оо
Корреляцию флуктуаций температуры среды в формуле (12) представляем в спектральном виде (8):
(вр (t') вр (t'O) =
t'
ds'e
о
t'-s' 1 те —
T©
ds''e
о
dwe-Ms'-s''V/ (ш). (13)
Интегралы в формуле (13) вычисляются в предположении t ^ t©:
t'
e т-
т©
t
s_
e т-
(1-гшге)
ds'
о
e т-т©
s_
e т-
(1+гшг-)
ds''
о
1
1 — iwT©
■ , il
e-,wt' - e т-
1
1 + iwT©
e^ - e
t-
e-iut!
1 — iwT© ’
e%ut"
1 + iwT©
Подставим эти формулы в выражение (13) и найдем связь между спектром флуктуаций температуры среды и корреляционной функцией флуктуаций температуры гранулы
(9р (t') 9р (t''))
e-,"C'-‘") ^ (ш) 2dw.
1 + (wT©)2
(14)
10
ISSN 1812-3368. Вестник МГТУ им. Н.Э. Баумана. Сер. “Естественные науки”. 2014. № 2
Согласно формуле (14), для времени большего, чем время тепловой релаксации гранулы, флуктуации температуры также являются статистически стационарным случайным процессом. Из формулы (14) можно получить выражения для автокорреляционной функции и квадрата дисперсии флуктуаций температуры гранулы:
<«,2) t («) <«2>
, м сю
(°/) [ е-гшг Фг (ш)
2 п J 1 + (шт©)2
-[
<f> [ Ф(ш) 2 dw.
2п J 1 + (шт©)2
-[
dw;
(15)
Рассмотрим два частных случая автокорреляционной функции флуктуаций температуры среды. Случайный процесс Of (t) — дельтакоррелированный во времени случайный процесс. Автокорреляционная функция (4) имеет вид
(Of (С) Of (?)) = <0?> tf (t - t") = <0?> 2то8 (t - t"), (16)
где t0 — временной микромасштаб.
Интегральный временной масштаб автокорреляционной функции (16) равен
tf (s) ds = 2t0 J 8 (s) ds = t0.
0 0
Спектр дельта-коррелированного случайного процесса получим из формул (10) и (16)
фf (ш) = 2то eiws8 (s) ds = 2t0. (17)
-[
В соответствии с формулами (9) и (17) энергия белого шума бесконечна. Подставив спектр дельта-коррелированного процесса (17) в формулу (15), получим выражение для корреляции флуктуаций температуры гранулы
<op> t (t)
е~ш1
2то
1 + (шт©)2
dw
0?
т t
—е т® т©
(18)
Интеграл в формуле (18) вычисляется методами теории вычетов. Из формулы (18) следуют выражения для квадрата дисперсии и для автокорреляционной функции флуктуаций температуры гранулы
0?> = 7
tp (t) = е те .
(19)
t
ISSN 1812-3368. Вестник МГТУ им. Н.Э.Баумана. Сер. “Естественные науки”. 2014. №2
11
Согласно формуле (19), дисперсия флуктуаций температуры гранулы меньше, чем среды. Интегральный временной масштаб автокорреляционной функции Фр (t) совпадает со временем тепловой релаксации частицы. Таким образом, уравнение (3), когда источником флуктуаций является белый шум, привело к статистически стационарному случайному процессу с экспоненциально спадающей автокорреляционной функции с конечным временем затухания.
Рассмотрим случай, когда $f (t) — случайный процесс с экспоненциально спадающей автокорреляционной функцией
t
Фf (t) = e-Te . (20)
Спектр этого случайного процесса
Ю
^f (ш) = J e*wt--TEdt
— OO
2 TE
1 + (шТЕ )2.
(21)
Этот спектр (см. (11)) соответствует случайному процессу с конечной энергией, но недифференцируемому. Из формул (15) и (21) получим выражение для корреляции флуктуаций температуры гранулы
Ю Фр (t)
/ о\ Ю
(f 0 2 TEe~<“‘
2п J [1 + (шт©)'2] [1 + (шТе)2] '
— Ю
Отметим, что этот случайный процесс является уже дифференцируемым. Интеграл вычисляется методами теории вычетов
tt
$Р> Фр (t) = <$
e Te - (t©/Te) e Te
(22)
1 - (t©/Te)2
При t = 0 из выражения (22) следует формула для интенсивности флуктуаций температуры гранулы
$р> =
1 + (t©/Te )
(23)
Видно, что можно выделить два типа гранул. Первый тип — мелкие гранулы, время тепловой релаксации которых существенно меньше интегрального временного масштаба флуктуаций температуры среды: t© ^ TE .В этом случае квадраты дисперсий флуктуаций температуры гранулы и среды близки: ($р) « ($2). Второй тип — крупные гранулы, время тепловой релаксации которых существенно больше, чем интегральный временной масштаб: t© ^ TE. Интенсивность флуктуаций температуры гранулы мала: ($р) « (TE/т©) ($2). С учетом формул (22) и (23) следует выражение для автокорреляционной функции
12
ISSN 1812-3368. Вестник МГТУ им. Н.Э. Баумана. Сер. “Естественные науки”. 2014. № 2
флуктуаций температуры гранулы
% (t) =
e Te - (tq/Te) e Te 1 — (tq /Te )
(24)
Интегральный временной масштаб флуктуаций температуры гранулы составляет
Ю
Tq = J Фр (t) dt = Te + tq. (25)
о
В соответствии с формулами (24) и (25) гранулы с малой инерцией tq ^ Т хорошо вовлекаются во флуктуации среды и их осредненные параметры близки к характеристикам среды:
t
Фр (t) w e - Т; Tq w Te.
Гранулы с большой тепловой инерцией tq ^ ТЕ образуют случайный процесс с автокорреляционной функцией и временньш масштабом, отличающимся от параметров среды,
_±_
Фр (t) w e -T® ; Tq w tq. (26)
Полученные точные результаты используются для тестирования численного алгоритма имитационного моделирования температуры гранулы в случайном поле.
Система стохастических уравнений. Тестирование алгоритма.
С учетом изложенного выше, моделирование случайного процесса с конечным временем затухания автокорреляционной функции возможно только на основе численного решения системы стохастических обыкновенных дифференциальных уравнений.
Рассмотрим систему уравнений, описывающих генерацию флуктуаций температуры в среде с конечным временем затухания автокорреляционной функции и в грануле с учетом тепловыделения:
dYf(т^ /ч , ч drp(t)_(1 + Yf(t)) — rp(t)
E*
= n(T )—y f (t );
fQ*e Mt), (27)
dT dT Qq
где n (t) — начальный источник флуктуаций, являющийся случайным процессом Гаусса с автокорреляционной функцией
<n Ю П (Л) = 2t0 (,2> S (t' — t").
При записи первого уравнения системы (27) учитывается, что временной масштаб совпадает с временным интегральным масштабом ТЕ. Интегрирование системы уравнений (27) осуществляется явным методом Эйлера [18]:
ISSN 1812-3368. Вестник МГТУ им. Н.Э. Баумана. Сер. “Естественные науки”. 2014. №2
13
Yf
(n+l)
Y/n+1) + Дr{n(n) - 7<n>} ;
P(n + 1) =
-L p
P
I /О* r(n
+ Q e lp
E *
n(n)
(28)
Здесь n — номер шага по времени; случайное приращение источ-
лого шума) с нулевым средним значением и единичной дисперсией. Шаг по времени в уравнениях (28) выбирается достаточно малым:
Согласно уравнениям (27) и (28), при отсутствии тепловыделения результаты аналитического исследования и результаты, полученные путем осреднения численных реализаций системы (28), должны быть близки.
Влияние тепловой инерции гранул на характер флуктуаций температуры без тепловыделения иллюстрирует рис. 3. Увеличение тепловой инерции приводит к снижению амплитуды флуктуаций и к их более гладкому поведению.
На рис. 4 показано сопоставление результатов по относительной дисперсии флуктуаций температуры гранулы, полученных путем осреднения данных имитационного моделирования системы (28) и по аналитической формуле (23). Увеличение тепловой инерции понижает интенсивность флуктуаций температуры гранул.
Автокорреляционные функции флуктуаций температуры частицы представлены на рис. 5. Результаты численного моделирования удовлетворительно согласуются с точными данными, определенными
УрУр
Рис. 3. Пример случайного поведения температуры среды и частиц при = 1
(я) и 3 (б) для Yp (-0 и Yf (2)
ника флуктуаций моделируется как Дтц(n) = £(n) л/2РПУт^Ат, где £(n) — случайная реализация нормированного процесса Гаусса (бе-
Дт = 10 2 min(1, О©). Сбор информации для получения осреднен-ных параметров начинается со времени т > 102 max (1, О©).
о
з
6
9
12 15 18 г
14
ISSN 1812-3368. Вестник МГТУ им. Н.Э. Баумана. Сер. “Естественные науки”. 2014. №2
методом спектрального анализа случайных процессов. Возрастание тепловой инерции гранулы увеличивает область затухания автокорреляционной функции.
Сопоставление точных аналитических результатов и данных численного эксперимента позволяет утверждать об адекватности предложенного в работе численного алгоритма моделирования, основанного на решении системы стохастических уравнений (28).
Моделирование теплового взрыва. Среднее время ожидания взрыва. В этом разделе приведены результаты, иллюстрирующие различные сценарии поведения температуры частицы с внутренним тепловыделением во флуктуирующей среде. Изменение актуальной температуры гранулы с тепловыделением показано на рис. 6, где rf = 1 + Yf (t) — актуальная температура среды. Флуктуации скорости химической реакции вносят заметный вклад в случайную температуру гранулы. Начальная температура гранулы меньше, чем критическое значение, соответствующее второму корню Гп, взятому с диаграммы Семенова. Случайный процесс с ненулевой вероятностью может выйти за любой уровень. По истечении некоторого случайного времени актуальная температура частицы превысит критическое значение Гсг = Гп и произойдет
Рис. 4. Отношение дисперсий флуктуаций температуры гранулы и среды по результатам моделирования (1) и полученные по формуле (23) (2)
Рис. 5. Автокорреляционные функции флуктуаций температуры гранулы (1) и среды (2) по результатам численного моделирования (о, д, V, □ ), и полученные по формулам (24) и (20) (кривые); штриховой линией показана корреляционная функция (26)
ISSN 1812-3368. Вестник МГТУ им. Н.Э. Баумана. Сер. “Естественные науки”. 2014. №2
15
Рис. 6. Изменения актуальных температур среды Гf (1) и гранулы Гр (2) до
теплового взрыва ((Yp)° — квадрат дисперсии флуктуаций температуры среды без учета тепловыделения) при ^© = 2
потеря тепловой устойчивости. Этот сценарий проиллюстрирован на рис. 7.
Время ожидания теплового взрыва определяем как среднее время первого пересечения случайными флуктуациями температуры гранулы критического уровня Гсг, среднее время ожидания теплового взрыва тсг в зависимости от начальной температуры гранулы Гр, с которой она подается во флуктуирующую среду, показано на рис. 8. По мере приближения начальной температуры гранулы к критической среднее время потери тепловой устойчивости резко сокращается. Критическая температура зависит от параметра тепловой инерции гранулы.
Заключение. Разработана процедура численного моделирования случайной температуры гранулы с внутренним тепловыделением
Рис. 7. Пример возникновения теплового взрыва для yf (1), YP (2) и Гд (3) при ^© = 4
16
ISSN 1812-3368. Вестник МГТУ им. Н.Э. Баумана. Сер. “Естественные науки”. 2014. №2
Рис. 8. Зависимость среднего времени ожидания теплового взрыва при = 1 (1), 1,5 (2) и 2 (3)
в жидкой среде с флуктуациями температуры. Интенсивность тепловыделения описывается законом Аррениуса.
Для флуктуаций температуры среды предложена процедура численной генерации случайного процесса Гаусса с экспоненциально затухающей автокорреляционной функцией. Автокорреляционные функции и дисперсии флуктуаций температуры гранулы без тепловыделения, полученные методом имитационного моделирования, сопоставляются с точными формулами, найденными путем спектрального анализа случайных процессов.
Представлен анализ влияния флуктуаций температуры среды на процесс возникновения теплового взрыва. Проиллюстрирована динамика теплового взрыва для различных значений параметра тепловой инерции начальной температуры частицы, уровня дисперсий температуры среды, мощности тепловыделений и энергии активации. На основе осреднения данных имитационного моделирования установлена зависимость среднего времени ожидания начала теплового взрыва. Обнаружен эффект стохастического дрейфа температуры частицы с внутренним тепловыделением к ее критическому значению.
Направление дальнейших исследований в области численного моделирования связано с привлечением в разработанный метод реальных кинетических схем процесса синтеза Фишера-Тропша на кобальтовых катализаторах, а также исследования, посвященного учету перемежаемости, характеризуемой логарифмически нормальным законом распределения случайной температуры среды.
Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (проект 11-08-00645-а).
ISSN 1812-3368. Вестник МГТУ им. Н.Э. Баумана. Сер. “Естественные науки”. 2014. №2
17
ЛИТЕРАТУРА
1. Steynberg A.P., Dry M.E., Davis B.H., Breman B.B. Fischer — Tropsch Reactors // Studies in Surface Science and Catalysis. Vol. 152. Fischer-Tropsch Technology. Amsterdam: Elsevier, 2004. P. 64-195.
2. Зельдович Я.Б., Баренблатт Г.И., Либрович В.Б., Махвиладзе Г.М. Математическая теория горения и взрыва. М.: Наука, 1980. 478 с.
3. Франк-Каменецкий Д.А. Диффузия и теплопередача в химической кинетике. М.: Наука, 1987. 502 с.
4. Мержанов А.Г., Руманов Э.ДНелинейные эффекты в макроскопической кинетике // Успехи физических наук. 1987. Т 151. № 4. С. 553-593.
5. Худяев С.И.Пороговые явления в нелинейных уравнениях. М.: Физматлит, 2003. 272 с.
6. Варнатц Ю., Маас У., Диббл Р. Горение. Физические и химические аспекты, моделирование, эксперименты, образование загрязняющих веществ. М.: Физматлит, 2003. 352 с.
7. Хорстхемке В., Лефевр Р. Индуцированные шумом переходы: Теория и применение в физике, химии, биологии. М.: Мир, 1987. 400 с.
8. Derevich I.V., Gromadskaya R.S. Rate of chemical reactions with regard to temperature fluctuations // Theoretical Foundations of Chemical Engineering. 1997. Vol. 31. No. 4. P. 392-397.
9. Федотов С.П., Третьяков М.В. Стационарные режимы гетерогенной химической реакции при наличии внешних шумов // Химическая физика. 1988. Т. 7. № 11. С. 1533-1537.
10. Федотов С.П., Третьяков М.В. О стохастическом воспламенении частицы // Химическая физика. 1991. Т. 10. № 2. С. 238-241.
11. Медведев В.Г., Телегин В.Г., Телегин Г.Г. Статистический анализ кинетики адиабатического теплового взрыва // Физика горения и взрыва. 2009. Т. 45. № 3. С. 44-48.
12. Деревич И.В., Зайчик Л.И. Уравнение для плотности вероятности скорости и температуры частиц в турбулентном потоке, моделируемом гауссовым случайным полем // Прикладная математика и механика. 1990. Т. 54. С. 767-774.
13. Derevich I.V. Effect of temperature fluctuations of fluid on thermal stability of particles with exothermic chemical reaction // Int. J. Heat Mass Transfer. 2010. V. 53. P. 5920-5932.
14. Деревич И.В. Влияние флуктуаций температуры среды на тепловой взрыв одиночной частицы // Физика горения и взрыва. 2011. Т. 47. № 5. С. 1-12.
15. Derevich I.V. Temperature oscillation in a catalytic particle of Fischer — Tropsch synthesis // Int. J. Heat and Mass Transfer. 2010. Vol. 53. P. 135-153.
16. Кляцкин В.И. Стохастические уравнения глазами физика (основные положения, точные результаты и асимптотические приближения). М.: Физматлит, 2001. 528 c.
17. Liang G.Y., Cao L., Wu D./.Approximate Fokker-Planck equation of system driven by multiplicative colored noises with colored cross-correlation // Physica A. 2004. Vol. 335. P. 371-384.
18. Gillespie D.T. Exact numerical simulation of the Ornstein-Uhlenbeck process and its integral // Physical Review E. 1996. Vol. 54. No. 2. P. 2084-2091.
19. Ilie S., Teslya A. An adaptive stepsize method for the chemical Langevin equation // J. Chem. Phys. 2012. Vol. 136. P. 184101 (14).
20. Гардинер К.В. Стохастические методы в естественных науках. М.: Мир, 1986. 527 c.
21. Леви П. Стохастические процессы и броуновское движение. М.: Наука, 1972. 376 с.
18
ISSN 1812-3368. Вестник МГТУ им. Н.Э. Баумана. Сер. “Естественные науки”. 2014. № 2
REFERENCES
[1] Steynberg A.P., Dry M.E., Davis B.H., Breman B.B. Fischer-Tropsch Reactors. Fischer-Tropsch Technology. Studies in Surface Science and Catalysis. Amsterdam: Elsevier, 2004. vol. 152, pp. 64-195.
[2] Zel’dovich Ya.B., Barenblatt G.I., Librovich V.B., Makhviladze G.M. Matematicheskaya teoriya goreniya i vzryva [Mathematical theory of combustion and explosion]. Moscow, Nauka Publ., 1974. 432 p.
[3] Frank-Kamenetskiy D.A. Diffuziya i teploperedacha v khimicheskoy Kinetike [Diffusion and heat transfer in chemical kinetics]. Moscow, Nauka Publ., 1987. 502 p.
[4] Merzhanov A.G., Rumanov E.N. Nonlinear effects in macroscopic kinetics. Usp. Fiz. Nauk [Sov. Phys.-Usp.], 1987, vol. 151, no. 4, pp. 553-593.
[5] Khudyaev S.I. Porogovye yavleniya v nelineynykh uravneniyakh [Threshold Phenomena in Nonlinear Equations]. Moscow, Fizmatlit Publ., 2003. 272 p.
[6] Varnatts Yu., Maas U., Dibbl R. Gorenie. (Russ. ed.: Varnatts Yu., Maas U., Dibbl R. Fizicheskie i khimicheskie aspekty, modelirovanie, eksperimenty, obrazovanie zagryaznyayushchikh veshchestv. Moscow, Fizmatlit Publ., 2003. 352 p.).
[7] Xofrstemaxe V., Air in R. Noise-tnduced transitions: Theory and applications in physics, chemistry, biology. New York: Wiley, 1987. 400 p. (Russ. ed.: Moscow, Mir Publ., 1987. 400 p.).
[8] Derevich I.V., Gromadskaya R.S. Rate of chemical reactions with regard to temperature fluctuations. Theoretical Foundations of Chemical Engineering, 1997, vol. 31, no. 4, pp. 392-397.
[9] Fedotov S. P, Tret’yakov M.V. Stationary regimes of heterogeneous chemical reaction in the presence of ambient noise. Khimicheskaya fizika [Chemical Phys.]. 1988, vol. 7, no. 11, pp. 1533-1537.
[10] Fedotov S.P., Tret’yakov M.V. Stochastic ignition particles. Khimicheskaya fizika [Chemical Physics], 1991, vol. 10, no. 2, pp. 238-241 (in Russ.).
[11] Medvedev V.G., Telegin V.G., Telegin G.G. Statistical analysis kinetics of adiabatic thermal explosion. Fizika goreniya i vzryva [Combustion, Explosion, and Shock Waves], 2009, vol. 45, no. 3, pp. 44-48 (in Russ.).
[12] Derevich I.V., Zaychik L.I. The equation for the probability density of the velocity and temperature of the particles in a turbulent flow, the simulated Gaussian random field. Prikladnaya matematika i mekhanika, 1990, vol. 54, pp. 767-774.
[13] Derevich I.V. Effect of temperature fluctuations of fluid on thermal stability of particles with exothermic chemical reaction. Int. J. Heat Mass Transfer, 2010, vol. 53, pp. 5920-5932.
[14] Derevich I.V. Influence of temperature fluctuations on the thermal explosion of single particle. Fizika goreniya i vzryva [Physics of combustion and explosion], 2011, vol. 47, no. 5, pp. 1-12.
[15] Derevich I.V. Temperature oscillation in a catalytic particle of Fischer-Tropsch synthesis. Int. J. Heat and Mass Transfer, 2010, vol. 53, pp. 135-153.
[16] Klyatskin EI.Stokhasticheskie uravneniya glazami fizika (osnovnye po- lozheniya, tochnye rezul’taty i asimptoticheskie priblizheniya) [Stochastic equations by eyes physics (basic framework, exact results and asymptotic approximations)]. Moscow, Fizmatlit Publ., 2001. 528 p.
[17] Liang G.Y., Cao L., Wu .J. Approximate Fokker-Planck equation of system driven by multiplicative colored noises with colored cross-correlation. Physica A. 2004, vol. 335. pp. 371-384.
[18] Gillespie D.T. Exact numerical simulation of the Ornstein - Uhlenbeck process and its integral. Physical Review E. 1996, vol. 54, no. 2, pp. 2084-2091.
[19] Ilie S., Teslya A. An adaptive stepsize method for the chemical Langevin equation. J.Chem. Phys., 2012, vol. 136, pp. 184101 (14).
ISSN 1812-3368. Вестник МГТУ им. Н.Э. Баумана. Сер. “Естественные науки”. 2014. №2
19
[20] Gardiner K.V. Stochastic methods in the natural sciences. New York: Wiley, (Russ. ed.: Stokhasticheskie metody v estestvennykh naukakh. Mir Publ., 1986. 527 p.).
[21] Levi P. Stokhasticheskie protsessy i brounovskoe dvizhenie [Stochastic processes and Brownian motion]. (Russ. ed.: Moscow, Nauka Publ., 1972. 376 p.
Статья поступила в редакцию 02.07.2013
Игорь Владимирович Деревич — д-р техн. наук, профессор кафедры “Прикладная математика” МГТУ им. Н.Э. Баумана. Автор 80 научных работ в области математического моделирования физических процессов и технических систем, термодинамики, турбулентных двухфазных потоков, теории случайных процессов, химической кинетики.
МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.
I.V. Derevich — Dr. Sci. (Eng.), professor of “Applied Mathematics” department of the Bauman Moscow State Technical University. Author of 80 publications in the field of mathematical simulation of physical processes and technical systems, thermodynamics, turbulent two-phase flows, theory of random processes, chemical kinetics.
Bauman Moscow State Technical University, Vtoraya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.
Дарья Денисовна Галдина — аспирантка кафедры “Прикладная математика” МГТУ им. Н.Э. Баумана. Автор семи научных работ в области математической статистики, реакторов Фишера-Тропша, тепловой стабильности, теплового взрыва.
МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.
D.D. Galdina — post-graduate of “Applied Mathematics” department of the Bauman Moscow State Technical University. Author of seven publications in the field of mathematical statistics, Fischer-Tropsch reactors, thermal stability, heat explosion. Bauman Moscow State Technical University, Vtoraya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.
20
ISSN 1812-3368. Вестник МГТУ им. Н.Э. Баумана. Сер. “Естественные науки”. 2014. № 2