Ксе к — различия между версиями

Материал из Викиконспекты
Перейти к: навигация, поиск
м (rollbackEdits.php mass rollback)
 
(не показано 14 промежуточных версий 5 участников)
Строка 1: Строка 1:
 +
''' О ЗАДАЧЕ '''
 +
 
Представим банку, заполненную хим акт жидкостью, а может твердое тело, перемешаны хим реагенты, которые способны взаимод друг с другом. Одну из стенок банки начинают прогревать (поддерживая температуру стенки <tex>T_w</tex>). При нагревании происходит хим реакция. Эта хим реакция удовлетворяет 2 свойствам:
 
Представим банку, заполненную хим акт жидкостью, а может твердое тело, перемешаны хим реагенты, которые способны взаимод друг с другом. Одну из стенок банки начинают прогревать (поддерживая температуру стенки <tex>T_w</tex>). При нагревании происходит хим реакция. Эта хим реакция удовлетворяет 2 свойствам:
 
# скорость реакции "сильно" увеличивается с температурой
 
# скорость реакции "сильно" увеличивается с температурой
Строка 45: Строка 47:
 
# <tex>T(t, z)</tex> - температура  
 
# <tex>T(t, z)</tex> - температура  
  
<tex>\left\{ \begin{matrix} \frac {\partial x} {\partial t} - D \frac{\partial^2 x}{\partial t^2} = W(x, T) \\  
+
<tex>\left\{ \begin{matrix} \frac {\partial x} {\partial t} - D \frac{\partial^2 x}{\partial z^2} = W(x, T) \\  
C \frac {\partial T}{\partial t} - \lambda \frac{\partial^2 T}{\partial z^2} = - \rho Q W (x,T)\end{matrix} \right.</tex>
+
\rho C \frac {\partial T}{\partial t} - \lambda \frac{\partial^2 T}{\partial z^2} = - \rho Q W (x,T)\end{matrix} \right.</tex>
<ref>
 
У меня немного по-другому 2-ое уравнение: <tex> \rho C \frac {\partial T}{\partial t} - \lambda \frac{\partial^2 T}{\partial z^2} = - \rho Q W (x,T) </tex>
 
</ref>
 
 
<br>, где <tex>D</tex> - коэффициет диффузии
 
<br>, где <tex>D</tex> - коэффициет диффузии
  
Строка 56: Строка 55:
 
<br>, где K - константа скорости реакции. К, а - порядок реакции, Е - энергия активациии - константы
 
<br>, где K - константа скорости реакции. К, а - порядок реакции, Е - энергия активациии - константы
  
Что такое переход из вещ А в В (РИСУНКО енергия связи, барьер.) То есть чтобы проихощла рекция необходимо преодолеть молек барьер. Экспонента формуле показывает, какая часть модекул больше барьера
+
Что такое переход из вещества  А в В (РИСУНОК енергия связи, барьер.) То есть чтобы произошла реакция необходимо преодолеть молекулярный барьер. Экспонента формуле показывает, какая часть модекул больше барьера.
  
Надо решуть ту систему уравнений.
+
Надо решить ту систему уравнений.
Граничные условия.
+
 
 +
* Граничные условия:
  
 
<tex>x|_{z = 0} = 0</tex>
 
<tex>x|_{z = 0} = 0</tex>
  
<tex>T|_{z = 0} = T_w</tex> - темпер стенки
+
<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>. На самом деле, все это не важно условия на дальнем конце, пока фронт не подойдет к ней.
  
Начальные условия
+
* Начальные условия:
  
 
<tex>x|_{t = 0} = \left\{ \begin{matrix} 1, z \ne 0 \\  
 
<tex>x|_{t = 0} = \left\{ \begin{matrix} 1, z \ne 0 \\  
Строка 76: Строка 76:
 
T_w, z = 0\end{matrix} \right.</tex>
 
T_w, z = 0\end{matrix} \right.</tex>
  
с вер 99 рпоцентов не получится, надо представлять структуру того, что происходит. То есть нельзя формально применять методы, должен быт предварительный физ анализ. Поэтому нужны оценки
+
* Замечание
 +
С вероятность 99 процентов не получится; надо представлять структуру того, что происходит. То есть нельзя формально применять методы, должен быть предварительный физ анализ. Поэтому нужны оценки.
  
Лценки:
 
  
Характерная величинаа скорости фронта для случая когда, порядок реакции а = 1
+
''' Оценки: '''
  
<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>
+
Характерная величина скорости фронта для случая когда, порядок реакции а = 1
  
К - конст реакции, <tex>\triangle T</tex> - насколько среда прогревается, <tex>\lambda</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>
Q - топловой эффект реакции
+
<br>, где
 +
<br>К - конст реакции,  
 +
<br><tex>\triangle T</tex> - насколько среда прогревается
 +
<br><tex>\lambda</tex> - коэффициент теплопроводности
 +
<br>Q - тепловой эффект реакции
  
<tex>\triangle T = T_m - T_0 = \frac{Q}{C}</tex>
+
<br><tex>\triangle T = T_m - T_0 = \frac{Q}{C}</tex>
 
+
<br> T_m - температура адиабатического прохожденя реакции, то есть насколько прогрелась
T_m - ьемпература адиабатического прохожденя реакции, то есть насколько прогрелась
 
  
 
По структуре фронта (ГРАФИКИ структура фронта)
 
По структуре фронта (ГРАФИКИ структура фронта)
есть сравнительно широкая зона подогрева <tex>\delta_t</tex> и сравнительно узкая зна реакции <tex>\delta_r</tex>. То есть температура увелич в сравнительно широкой облачти, а реакция контертруется в более узкой зоне.
+
есть сравнительно широкая зона подогрева <tex>\delta_t</tex> и сравнительно узкая зона реакции <tex>\delta_r</tex>. То есть температура увелич в сравнительно широкой облачти, а реакция контертруется (?) в более узкой зоне.
  
<tex>\delta_T \sim \frac {\varkappa}{U} = \frac{\lambda}{p c  U}</tex>, <tex>\varkappa</tex>- коэфф темепературопроводности
+
<tex>\delta_T \sim \frac {\varkappa}{U} = \frac{\lambda}{\rho С U}</tex>, <tex>\varkappa</tex>- коэфф темепературопроводности
  
 
диффузионный масштаб (может не совпадать с тепловым)
 
диффузионный масштаб (может не совпадать с тепловым)
<tex>\delta_D \sim D/U </tex> D - коэфф диффузии
+
<tex>\delta_D \sim \frac {D} {U} </tex>, где D - коэфф диффузии
  
<tex>\delta_r \sim \delta_T \beta</tex> ??
+
<tex>\delta_r \sim \delta_T \beta</tex>
 +
<br>, где <tex>\beta =\frac{R T_m}{E} \ll 1</tex> - условние "сильной " зависимости скор реакц от температуры
  
<tex>\beta =\frac{R T_m}{E} \ll 1</tex> - условние "сильной " зависимости скор реакц от темпертуры
+
<tex>\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 = \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>\delta_r</tex> укладывалось хотя юы несколько пространственных шагов  ,  
 
# на <tex>\delta_r</tex> укладывалось хотя юы несколько пространственных шагов  ,  
 
# <tex> \triangle z\lesssim \delta_r</tex>,
 
# <tex> \triangle z\lesssim \delta_r</tex>,
# <tex>\delta_T \ll l </tex> l - разсер области, то еть чтоб фрон поместился.
+
# <tex>\delta_T \ll l </tex> l - разсер области, то есть чтобы фронт поместился.
 +
 
 +
''' Задача '''
  
Предже всего получить обычный фронт, потом варьируя параметры залезть за критичсекие режимы. Что способствует переходу за крит режимы: D↓, K↑, и одновременно (K↑, Е↑ таким образом что <tex>K e^{-\frac{E}{l t m}} = const </tex>)
+
Предже всего, получить обычный фронт, потом варьируя параметры залезть за критичсекие режимы. Что способствует переходу за крит режимы: D↓, K↑, и одновременно (K↑, Е↑ таким образом что <tex>K e^{-\frac{E}{l t m}} = const </tex> - может привести к релаксационным колебаниям)
  
 
(*)Для желающих 2мерную задачу.  
 
(*)Для желающих 2мерную задачу.  
Строка 120: Строка 125:
 
<tex>E = 8 \cdot 10^4 </tex>  Дж/Моль энергия активации
 
<tex>E = 8 \cdot 10^4 </tex>  Дж/Моль энергия активации
  
<tex>R = 8.314 </tex> Дж/(Моль * К)
+
<tex>R = 8.314 </tex> Дж/(Моль * К) унив газовая постоянная
  
 
<tex>a = 0..2</tex> - порядок реакции. лучше начинать с 1
 
<tex>a = 0..2</tex> - порядок реакции. лучше начинать с 1
Строка 132: Строка 137:
 
<tex>C = 1980</tex> Дж/(кг * K) теплоемкость
 
<tex>C = 1980</tex> Дж/(кг * K) теплоемкость
  
<tex>\lambda = 1.13 </tex> Дж/(м * с * К) теплопроводность  
+
<tex>\lambda = 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>. Это даст ситуацию подобия уравнений переноса тепла и переноса массы.  
+
<tex>D \sim 8 \cdot 10^{-12}</tex> м^2/c коэффиц диффузии.
 +
Диффузия в жидк и твердых телах очень маленькая. Для начала не реальную брать D, а взять не физ значение а такое, что число Льюиса <tex>L_e = \frac{D}{\varkappa} = \frac{D \rho C} {\lambda} = 1</tex>. Это даст ситуацию подобия уравнений переноса тепла и переноса массы.  
  
"Препроцессинг" - интерактивный ввод параметров физических и вычислительных(шаги колво шагов...)
+
"Препроцессинг" - интерактивный ввод параметров физических и вычислительных(шаги кол-во шагов...)
  
 
"Процессор" - солвер
 
"Процессор" - солвер
Строка 142: Строка 151:
 
"Постпроцессор" - визуализация. Температура, концентрация, W скорость реакции. (интересно - анимация, прям волна бежит)
 
"Постпроцессор" - визуализация. Температура, концентрация, W скорость реакции. (интересно - анимация, прям волна бежит)
  
   
+
 
=== Возможные альтернативные варианты формул: ===
+
== Решения ==
 +
=== Неявный метод ===
 +
 
 +
Общий вид:
 +
 
 +
<tex> \frac{T^{n+1}-T^{n}}{\Delta t}=L_{n}T^{n+1} </tex>
 +
 
 +
 
 +
В нашем случае (вычисляем <tex>(n+1)</tex>-ый слой):
 +
 
 +
<tex> \frac {X_{i}^{n+1} - X_{i}^{n}} {\Delta t}  -  D \frac {X_{i-1}^{n+1} - 2X_{i}^{n+1} + X_{i+1}^{n+1}} {\Delta z^{2}}  = W(X_{i}^{n},\ T_{i}^{n}) </tex>
 +
<ref> В слагаемом с <tex> D </tex> у нас была производная по времени, но тогда мало понятно как решать, а ещё решение которое нашёл Вова содержит в этом месте производную по координате. </ref>
 +
 
 +
<tex> \rho C \frac {T^{n+1} - T^{n}} {\Delta t} - \lambda \frac {T_{i-1}^{n+1} - 2T_{i}^{n+1} + T_{i+1}^{n+1}} {\Delta z^{2}}  = -\rho Q W(X_{i}^{n},\ T_{i}^{n}) </tex>
 +
 
 +
Решается методом прогонки ( + внутренние итерации)
 +
 
 +
== Возможные альтернативные варианты формул: ==
 
   
 
   
 
<references/>
 
<references/>

Текущая версия на 19:40, 4 сентября 2022

О ЗАДАЧЕ

Представим банку, заполненную хим акт жидкостью, а может твердое тело, перемешаны хим реагенты, которые способны взаимод друг с другом. Одну из стенок банки начинают прогревать (поддерживая температуру стенки [math]T_w[/math]). При нагревании происходит хим реакция. Эта хим реакция удовлетворяет 2 свойствам:

  1. скорость реакции "сильно" увеличивается с температурой
  2. происходит "сильное" выделение тепла этой реакции

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

  • Как ведет концентрация реагентов?

начальный реагент A -> B (в продукт B) (в чем мер концентрация = отношение плотности вещества к полной плотности смеси [math]x = \frac {\rho_A} {\rho_A + \rho_B} [/math] , [math]0 \leq x \leq 1[/math])

Перед фронтом когда один реагент концентрация = 1. После фронта асимтотически выходит на 0.

  • Как ведет себя скорость?

Дает очень узкий пик в какой-то малой зоне. Перед зоной скорость реакции мала, так как температура мала, а после мала, так как реагент скушался. Расчеты показывают, что это очень узкий пик.

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

  • [math] T_w \lesssim T^* \le T_m [/math] - нет "поджига" (т.е. волна не начинается, инертный прогрев)
  • [math] T^* \lesssim T_w \le T_m [/math] - "поджиг" с задержкой
  • [math] T_m \le T_w [/math] - быстрый "поджиг" (мнгновенно)

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

Совершенно детерминир система - такое сложное поведение =)


КАК МОДЕЛИРОВАТЬ

В одномерном случае ситема описывается 2мя функциями:

  1. [math]x(t, z)[/math] - концентрация
  2. [math]T(t, z)[/math] - температура

[math]\left\{ \begin{matrix} \frac {\partial x} {\partial t} - D \frac{\partial^2 x}{\partial z^2} = W(x, T) \\ \rho C \frac {\partial T}{\partial t} - \lambda \frac{\partial^2 T}{\partial z^2} = - \rho Q W (x,T)\end{matrix} \right.[/math]
, где [math]D[/math] - коэффициет диффузии

Первое - уравнение диффузии. Справа - скорость хим реакции:
[math]W(x, T) = - K x^a exp ( - \frac{E}{R T})[/math]
, где K - константа скорости реакции. К, а - порядок реакции, Е - энергия активациии - константы

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

Надо решить ту систему уравнений.

  • Граничные условия:

[math]x|_{z = 0} = 0[/math]

[math]T|_{z = 0} = T_w[/math] - температура стенки

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

  • Начальные условия:

[math]x|_{t = 0} = \left\{ \begin{matrix} 1, z \ne 0 \\ 0, z = 0\end{matrix} \right.[/math]


[math]T|_{t = 0} = \left\{ \begin{matrix} T_0, z \ne 0 \\ T_w, z = 0\end{matrix} \right.[/math]

  • Замечание

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


Оценки:

Характерная величина скорости фронта для случая когда, порядок реакции а = 1

[math]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}[/math]
, где
К - конст реакции,
[math]\triangle T[/math] - насколько среда прогревается
[math]\lambda[/math] - коэффициент теплопроводности
Q - тепловой эффект реакции


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

По структуре фронта (ГРАФИКИ структура фронта) есть сравнительно широкая зона подогрева [math]\delta_t[/math] и сравнительно узкая зона реакции [math]\delta_r[/math]. То есть температура увелич в сравнительно широкой облачти, а реакция контертруется (?) в более узкой зоне.

[math]\delta_T \sim \frac {\varkappa}{U} = \frac{\lambda}{\rho С U}[/math], [math]\varkappa[/math]- коэфф темепературопроводности

диффузионный масштаб (может не совпадать с тепловым) [math]\delta_D \sim \frac {D} {U} [/math], где D - коэфф диффузии

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

[math]\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[/math] - условие "сильной" экзотермичности реакии

  • Как подбирать шаги по времени?

Должны разрешить наименьший физ масштаб. нужно чтобы

  1. на [math]\delta_r[/math] укладывалось хотя юы несколько пространственных шагов ,
  2. [math] \triangle z\lesssim \delta_r[/math],
  3. [math]\delta_T \ll l [/math] l - разсер области, то есть чтобы фронт поместился.

Задача

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

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

Параметры:

[math]K = 1.6 \cdot 10^6 [/math] 1 /c константа скорости реакции

[math]E = 8 \cdot 10^4 [/math] Дж/Моль энергия активации

[math]R = 8.314 [/math] Дж/(Моль * К) унив газовая постоянная

[math]a = 0..2[/math] - порядок реакции. лучше начинать с 1

[math]Q =7 \cdot 10^5 [/math] Дж/кг тепловой эффект реакции

[math]\rho = 830 [/math] кг / м^3

[math]T_0 = 293[/math] K

[math]C = 1980[/math] Дж/(кг * K) теплоемкость

[math]\lambda = 0.13 [/math] Дж/(м * с * К) теплопроводность [1]

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

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

"Процессор" - солвер

"Постпроцессор" - визуализация. Температура, концентрация, W скорость реакции. (интересно - анимация, прям волна бежит)


Решения

Неявный метод

Общий вид:

[math] \frac{T^{n+1}-T^{n}}{\Delta t}=L_{n}T^{n+1} [/math]


В нашем случае (вычисляем [math](n+1)[/math]-ый слой):

[math] \frac {X_{i}^{n+1} - X_{i}^{n}} {\Delta t} - D \frac {X_{i-1}^{n+1} - 2X_{i}^{n+1} + X_{i+1}^{n+1}} {\Delta z^{2}} = W(X_{i}^{n},\ T_{i}^{n}) [/math] [2]

[math] \rho C \frac {T^{n+1} - T^{n}} {\Delta t} - \lambda \frac {T_{i-1}^{n+1} - 2T_{i}^{n+1} + T_{i+1}^{n+1}} {\Delta z^{2}} = -\rho Q W(X_{i}^{n},\ T_{i}^{n}) [/math]

Решается методом прогонки ( + внутренние итерации)

Возможные альтернативные варианты формул:

  1. У меня немного по-другому [math]\lambda = 0.13 [/math] Дж/(м * с * К)
  2. В слагаемом с [math] D [/math] у нас была производная по времени, но тогда мало понятно как решать, а ещё решение которое нашёл Вова содержит в этом месте производную по координате.