Изменения

исходная смесь при темп T0 (при которой скорость реакц очень маленькая). начинаем стеночку прогреватьПредставим банку, заполненную хим акт жидкостью, тепло передается близлеж слоям смесейа может твердое тело, нагреваютсяперемешаны хим реагенты, так как скорость реакц увеличивается которые способны взаимод друг с темп, в них начинает происходить реак, как только пошла, начинается выделение собственного тепла реакциидругом. Это тепло передается след слоям, они тоже прогремваются и тдОдну из стенок банки начинают прогревать (поддерживая температуру стенки <tex>T_w</tex>). При некоторых условий формируется тепловоцй фронт химической реакциинагревании происходит хим реакция. (РИСУНОК распределение температуры, ось z, Эта хим волна реакции).реакция удовлетворяет 2 свойствам:Tm = T0 + QC (кажется, я ошиблась в формуле)Tm температура адиаьатического прохождения # скорость реакции, когда все тепло реакции на нагрев смеси"сильно" увеличивается с температуройQ - тепловой эффект хим # происходит "сильное" выделение тепла этой реакцииС - теплоемкость

Исходная смесь имеет температуру <tex>T_0</tex> (при которой скорость реакции очень маленькая). Начинаем стеночку прогревать, тепло передается близлежащим слоям смесей, нагреваются, так как скорость реакции увеличивается с температурой, в них начинает происходить реакция. Как только пошла, начинается выделение собственного тепла реакции. Это тепло передается след слоям, они тоже прогреваются и тд. При некоторых условиях формируется тепловой фронт химической реакции. (РИСУНОК распределение температуры, ось z, хим волна реакции).<br> <tex>Tm = T_0 + \frac{Q}{C}</tex><br>, где <br><tex>T_m</tex> температура адиабатического прохождения реакции, когда все тепло реакции на нагрев смеси<br><tex>Q</tex> - тепловой эффект хим реакции<br><tex>C</tex> - теплоемкость * Как ведет концентрация реагинтовреагентов?

(в чем мер конц <b>концентрация</b> = отношение плотности вещ вещества к полной плотности смеси <tex>x = \frac {\rho_A} {\rho_A + \rho_B} </tex> , <tex>0 \leq x \leq 1</tex>) Перед фронтом когда один реагент концентрация = 1. После фронта асимтотически выходит на 0.

перед фронтом когда один геагент концентраци я = 1* Как ведет себя скорость?Дает очень узкий пик в какой-то малой зоне. Перед зоной скорость реакции мала, так как температура мала, а после фронта асимтотически выходит на 0мала, так как реагент скушался. Расчеты показывают, что это очень узкий пик.

Например, как ведет себя скоростьинициировать фронт?дает оче узкий пик в кокой-то малой зоне. перед зоной скорость реакц малаДопустим, так как температура малаустанавливаем температуру стенки <tex>T_w</tex>; если <tex>T_w = T_m</tex>, а после малато волна без проблем идет, так как реагент скушался. расчеты показываеютесли меньше, что это оче узкий пикто существует критическое значение <tex>T^*</tex> для инициирования волны. Тогда

Наример как фронт инициировать? Допустим устанавливаем температуру стенки Tw, есди эта Tw = Tm то волна без проблем идет, если меьшн, то есть критической значения * <tex> T_w \lesssim T^* для инициирования волны\le T_m </tex> - нет "поджига" (т.е.волна не начинается, инертный прогрев)

* <tex>T_0 T^* \lesssim T* T_w \le TmT_m </tex> - нет "поджиг" (наверное, тут вместо T* Tw)с задержкой

<tex>T* \lesssim Tw \le Tm</tex> - "поджиг" с задержкой<tex> Tm T_m \le TwT_w </tex> - быстрый "поджиг" (мнгновенно)

Когда волна отходит , она забывает об начальном условии. Влияние другой стенки. и тп.При определенном соотношении параметров, которые характеризуют эту волну, она может терять устойчивость. Что происходит посде после потери. ? Если терят теряет в одномерной моде (?) то есть сохраняет свою плоскую структуру, то формиру.тся формирются другие устойчивые режимы, например колеебтельныеколебательные, тоесть то есть волна движется, то ускор ускоряясь, то замедляясь, ; дальше мождет поизойти может произойти бифуркация, и воникнуть 2х периодические колебание, то есть делает таие такие колебания с большим периодом и маленьким.И при определенном наборе параметров возникает хаотичское хаотическое поведение, волна , сохраняя плоскую форму , распространяется колебательно,но вообще не периодичсекипериодически, поведение похоже на чаотичскоехаотическое. пример динам Пример динамического хаоса. : поведение похоже на хаосЮ хаос, но описывается детерминир детерминированной закономерностью.

Если 3д , то очаги(2шт) по спирали двигаются в одну сторону.Могут распасться на несколько очагов - спиновая воллнаволна. Всякие чудеса Совершенно детерминир система - такое сложное поведение =)  ''' КАК МОДЕЛИРОВАТЬ '''

Совершенно детерминир система В одномерном случае ситема описывается 2мя функциями:# <tex>x(t, z)</tex> - такое сложное поведениеконцентрация# <tex>T(t, z)</tex> - температура

КАК МОДЕЛИРОВАТЬв одномеррном случае ситема опис 2мя фунцик<tex>\left\{ \begin{matrix} \frac {\partial x} {\partial t} - D \frac{\partial^2 x}{\partial z^2} = W(tx, zT) \\ \rho C \frac {\partial T}{\partial t} - кончентрация , температура \lambda \frac{\partial^2 T}{\partial z^2} = - \rho Q W (tx, zT)\end{matrix} \right.</tex><br>, где <tex>D</tex> - коэффициет диффузии

Первое - уравнение диффузии. Справа - скорость хим реакции:<br><tex>\left\{ \begin{matrix} \frac {\partial x} {\partial t} - D \frac{\partial^2 x}{\partial t^2} = W(x, T) \\ C = - K x^a exp ( - \frac {\partial TE}{\partial t} - \lambda \frac{\partial^2 R T}{\partial z^2} = - \rho Q W (x,T)\end{matrix} \right.</tex>D <br>, где K - константа скорости реакции. К, а - коэффициет диффузиипорядок реакции, Е - энергия активациии - константы

первое - уравнение диффузии. справа скорость хим реакции<tex>WЧто такое переход из вещества А в В (xРИСУНОК енергия связи, T) = - K x^a exp ( - \frac{E}{R T}барьер.)</tex>K - константа скорости реакцииТо есть чтобы произошла реакция необходимо преодолеть молекулярный барьер. КЭкспонента формуле показывает, а - порядок реакции, Е - енергия активациии - константыкакая часть модекул больше барьера.

Что такое переход из вещ А в В (РИСУНКО енергия связи, барьерНадо решить ту систему уравнений.) То есть чтобы проихощла рекция необходимо преодолеть молек барьер. Экспонента формуле показывает, какая часть модекул больше барьера

Надо решуть ту систему уравнений.* Граничные условия.:

<tex>T|_{z = 0} = T_w</tex> - темпер температура стенки <tex>\frac {\partial x} {\partial z} |_{z = l} = 0, \frac {\partial T} {\partial z} |_{z = l} = 0</tex>. На самом деле, все это не важно условия на дальнем конце, пока фронт не подойдет к ней.

<tex>\frac {\partial x} {\partial z} |_{z = l} = 0, \frac {\partial T} {\partial z} |_{z = l} = 0</tex> На самом деле все это не важно услоивя на дальнем конце, пока фронт не подойдет к ней.* Начальные условия:

с вер * ЗамечаниеС вероятность 99 рпоцентов процентов не получится, ; надо представлять структуру того, что происходит. То есть нельзя формально применять методы, должен быт быть предварительный физ анализ. Поэтому нужны оценки.

<tex>U \sim [\frac{2 K \lambda}{Q \rho \triangle T} (\frac{R T m^2}{E})^2 e^{-\frac{E}{R T m}}] ^ {Характерная величина скорости фронта для случая когда, порядок реакции а = 1/2}</tex>

<tex>U \sim [\frac{2 K \lambda}{Q \rho \triangle T} (\frac{R T_{m^2}}{E})^2 e^{-\frac{E}{R T_m}}] ^ {1/2}</tex><br>, где<br>К - конст реакции, <br><tex>\triangle T</tex> - насколько среда прогревается, <br><tex>\lambda</tex> - коэффициент теплопроводности<br>Q - топловой тепловой эффект реакции

<br><tex>\triangle T = T_m - T_0 = \frac{Q}{C}</tex><br> T_m - ьемпература температура адиабатического прохожденя реакции, то есть насколько прогрелась

есть сравнительно широкая зона подогрева <tex>\delta_t</tex> и сравнительно узкая зна зона реакции <tex>\delta_r</tex>. То есть температура увелич в сравнительно широкой облачти, а реакция контертруется (?) в более узкой зоне.

<tex>\delta_T \sim \frac {\varkappa}{U} = \frac{\lambda}{p c \rho С U}</tex>, <tex>\varkappa</tex>- коэфф темепературопроводности

<tex>\delta_D \sim \frac {D/} {U } </tex> , где D - коэфф диффузии

<tex>\delta_r \sim \delta_T \beta</tex> ??<br>, где <tex>\beta =\frac{R T_m}{E} \ll 1</tex> - условние "сильной " зависимости скор реакц от температуры

<tex>\beta gamma =\frac{R T_m^2}{E\triangle T} = \frac{R T_m^2}{E (T_m - T_0)} = \frac{R T_m^2 C}{E Q} \ll 1</tex> - условние условие "сильной " зависимости скор реакц от темпертурыэкзотермичности реакии

* Как подбирать шаги по времени?Должны разрешить наименьший физ масштаб. нужно чтобы # на <tex>\gamma = delta_r</tex> укладывалось хотя юы несколько пространственных шагов , # <tex> \frac{R T_m^2}{E triangle z\triangle T} = lesssim \frac{R T_m^2}{E (T_m - T_0)} = delta_r</tex>,# <tex>\frac{R T_m^2 c}{E Q} delta_T \ll 1l </tex> l - условие "сильной" экзотермичности реакииразсер области, то есть чтобы фронт поместился.

Предже всего , получить обычный фронт, потом варьируя параметры залезть за критичсекие режимы. Что способствует переходу за крит режимы: D↓, K↑, и одновременно (K↑, Е↑ таким образом что <tex>K e^{-\frac{E}{l t m}} = const </tex>- может привести к релаксационным колебаниям)

(*)Для желающих 2мерную задачу.

<tex>K = 1,.6 * \cdot 10^6 </tex> 1 /c константа скорости реакции

<tex>E = 8 * \cdot 10^4 </tex> Дж/Моль энергия активации

<tex>R = 8,.314 </tex> Дж/(Моль * К)унив газовая постоянная

<tex>Q =7 * \cdot 10^5 </tex> Дж/кг тепловой эффект реакции

<tex>\lambda = 1,0.13 </tex> Дж/(м * с * К) теплопроводность <ref>У меня немного по-другому <tex>\lambda = 0.13 </tex> Дж/(м * с * К)</ref>

<tex>D \sim 8 * \cdot 10^{-12}</tex> м^2/c коэффиц диффуздиффузии. Диффузия в жидк и твердых телах очень маленькая. для Для начала не реальную юрать брать D, а звять взять не физ значение а такое, что число Льюиса <tex>L_e = \frac{D}{\varkappa} = \frac{D \rho C} {\lambda} = 1</tex>. Это даст ситуацию подобия уравнений переноса тепла и переноса массы.

"Препроцессинг" - интерактивный ввод параметров физических и вычислительных(шаги колво кол-во шагов...)