Belousov-Zhabotinsky

Belousov-Zhabotinsky Reaction

A reação de Belousov-Zhabotinsky (BZ) consiste em uma família de reações químicas oscilatórias descobertas inicialmente por Belousov, e posteriormente analisadas por Zhabotinsky. A reação consiste em 3BrO₃⁻ + 5CH₂(CO₂H)₂ + 3H⁺ → 3BrCH(CO₂H)₂ + 4CO₂ + 5H₂O + 2CH₂O₂, e demonstra um comportamento oscilatório não linear até atingir o equilíbrio químico (adicionar imagem da reação). A interação entre a reação e a difusão dos produtos químicos no espaço resultará na auto-organização de ondas viajantes dinâmicas. Seu mecanismo original, foi descrito através de 27 espécies químicas e um total de 80 reações.

Oregonator

Oregonator é um modelo matemático utilizado para descrever de forma mais simples a dinâmica da reação BZ, desenvolvido por Field e Noyes (1974). Foi um modelo não espacial originalmente composto por três variáveis de estado, onde posteriormente, vemos que tornam-se apenas duas. O mecanismo é, inicialmente composto por cinco etapas irreversíveis, onde, A = 3BrO₃ ^-, B = 5CH₂(COOH)₂; 2HCOOH, 3BrCH(COOH)₂ (no geral, estas e demais espécies orgânicas); P = HOBr; X = HBrO₂; Y = Br^-; Z = forma oxidada do catalisador e f = Coeficiente estequiométrico.

		A + Y	$⟶$	X + P
		X + Y	$⟶$	2 P
		A + X	$⟶$	2 X + 2 Z
		2 X	$⟶$	A + P
		B + Z	$⟶$	$\frac{1}{2}f$ Y

Aplicando, então, as equações de taxa, onde v é a taxa da reação e k_i corresponde às constantes de taxa de reação:

v1 = k1 [A][Y]

v2 = k2 [X][Y]

v3 = k3 [A][X]

v4 = k4 [X]2

v5 = k5 [B][Z]

Para construir o modelo Oregonator, é necessário supor que as concentrações de A e B permaneçam constantes (estão associadas às concentrações iniciais dos precursores). Posteriormente, deve-se aplicar as técnicas padrão de cinética química para obter o modelo dinâmico considerando X, Y e Z como variáveis dinâmicas, assumindo que as reações químicas são elementares, ou seja, os coeficientes estequiométricos coincidem com a potência das variáveis dinâmicas. Considerando $\tau$ como o tempo, vemos as seguintes equações de velocidade:

$\frac{d[X]}{d\tau }=k_1[A][Y]-k_2[X][Y]+k_3[A][X]-2k_4[X{]}^2$

$\frac{d[Y]}{d\tau }=-k_1[A][Y]-k_2[X][Y]+\frac{1}{2}f{k}_5[B][Z]$

$\frac{d[Z]}{d\tau }=2k_3[A][X]-k_5[B][Z]$

A análise é simplificada convertendo essas equações em uma forma adimensional:

$x\equiv \frac{2k_4[X]}{k_3[Y]}$

$y\equiv \frac{k_2[X]}{k_3[A]}$

$z\equiv \frac{k_5{k}_4[B][Z]}{(k_3[A]{)}^2}$

$t\equiv k_5[B]\tau$

A partir de operações algébricas com as equações acima, obtemos para x, y e z, o seguinte sistema de equações diferenciais não lineares:

$\frac{dx}{dt}=\frac{qy-xy+x(1-x)}{ϵ}$

$\frac{dy}{dt}=\frac{-qy-xy+fz}{{ϵ}^{\prime }}$

$\frac{dz}{dt}=x-z$

Onde $ϵ\equiv \frac{k_5[B]}{k-3[A]}$, ${ϵ}^{\prime }\equiv \frac{2k_4{k}_5[B]}{k_2{k}_3[A]}$ e $q\equiv \frac{2k_4{k}_1}{k_2{k}_3}$. Como parâmetro $$, é possível considerar a aproximação do estado estacionário da variável y, portanto, $y\equiv \frac{fz}{q+x}$ então as equações são reduzidos para:

$ϵ\frac{dx}{dt}=x(1-x)+f\frac{q-x}{q+x}z$

$\frac{dz}{dt}=x-z$

Dessa forma, o modelo de Oregonator mostra a forma típica de um sistema de feedback químico, ou seja, a variável x, que será reescrito como u, funciona como um ativador, enquanto a variável z, que será reescrita como v, tem o papel de inibidor. Se para as equações termos associados à difusão são adicionados, onde D_u e D_v são os coeficientes de difusão adimensionais, e ${\mathrm{\nabla }}^2$ é o operador Laplaciano, então, o sistema torna-se:

$ϵ\frac{du}{dt}=u(1-u)+f\frac{q-u}{q+u}v+D_u{\mathrm{\nabla }}^{2}u$

$\frac{dv}{dt}=u-v+D_v{\mathrm{\nabla }}^{2}v$

teste

condições iniciais

$u_{i,j}^0=1$ se 0 < 8(0.01i - 0.5) < (0.01j - 0.5) senão = 0

$v_{i,j}^0=1$ se 0 < -(0.01j - 0.5) < 8(0.01i - 0.5) senão = 0

$u_{i,j}^{n+1}=u_{i,j}^n+\frac{(u_{i,j}^n(1-u_{i,j}^n)-f{v}_{i,j}^n(u_{i,j}^n-q)}{(u_{i,j}^n+q)+\mathrm{\nabla }u(\frac{(u_{i+1,j}^n+u_{i-1,j}^n+u_{i,j+1}^n+u_{i,j-1}^n-4u_{i,j}^n)}{D{h}^2})}\frac{dt}{e}$

$v_{i,j}^{n+1}=v_{i,j}^n+(u_{i,j}^n-v_{i,j}^n+\mathrm{\nabla }v(\frac{(v_{i+1,j}^n+v_{i-1,j}^n+v_{i,j+1}^n+v_{i,j-1}^n-4v_{i,j}^n)}{D{h}^2})dt$