Uso ideal da divisão de Strang (para equação de difusão da reação)

Fiz uma observação estranha ao calcular a solução para uma simples equação de difusão da reação 1D:

$\frac{\partial}{\partial t}a=\frac{\partial^2}{\partial x^2}a-ab$

$\frac{\partial}{\partial t}b=-ab$

$\frac{\partial}{\partial t}c = a$

O valor inicial de é uma constante ( ), e só estou interessado na integral sobre de a ( ). O propósito de e a equação é apenas para avaliar esta integral.$b$$b(0,x)=b_0$$a$$0$$1$$\int_0^1a(t,x)dt$$c$$\frac{\partial}{\partial t}c = a$

Usei um esquema de divisão de Strang para o acoplamento entre difusão e reação (uma reação de meio passo, depois uma difusão de passo completo e, novamente, uma reação de meio passo), um esquema de Crank Nicholson para a difusão e uma solução analítica para a reação ( incluindo a equação ).$\frac{\partial}{\partial t}c = a$

Como uma etapa da solução analítica era mais que um fator 3 mais lenta que uma etapa do esquema Crank Nicholson, tentei fazer mais de uma etapa Crank Nicholson para cada etapa da reação. Eu esperava seguir com menos etapas do esquema de divisão do Strang, para que eu fosse mais rápido no geral.

No entanto, o efeito oposto pode ser observado, a saber, que muito mais etapas para o esquema de divisão de Strang são necessárias se mais de uma etapa da Crank Nicholson for usada. (Estou preocupado apenas com a precisão da integral sobre , que parece convergir mais rápido que própria.) Depois de pensar por algum tempo, notei que o mesmo efeito também ocorre para , e que eu até entendo o porquê desse caso. O ponto é que, se eu fizer exatamente uma etapa do Crank Nicholson, o esquema geral se transformará em uma regra trapezoidal (se ).$a$$a$$b(t,x)=b_0=0$$b(t)=0$

Portanto, se eu tratasse como parte da etapa de difusão, aumentar o número de etapas da Crank Nicholson (provavelmente) não levaria a uma precisão geral reduzida (como observado). Mas isso parece anular o propósito de usar uma solução analítica para a parte de reação (não linear e potencialmente muito rígida) do sistema.$\frac{\partial}{\partial t}c = a$

Então, eis a minha pergunta: existe uma maneira melhor de tratar no contexto de uma divisão de Strang, do que tratá-lo como parte da etapa da reação ou tratar como parte da etapa de difusão. Quero evitar ser "forçado" a usar exatamente uma etapa do Crank Nicholson para a difusão. (Por exemplo, em 3D, eu preferiria resolver a difusão analiticamente por uma FFT em vez de usar Crank Nicholson. É claro que também posso combinar a FFT com a Crank Nicholson, para que não seja tão importante.)$\frac{\partial}{\partial t}c = a$

Vou escrever isso como uma resposta, embora não responda diretamente à pergunta.

Ao conectar a segunda equação e a terceira equação na primeira e conectar a terceira na segunda, dê-se: Reorganizando essas duas opções: Agora, podemos integrar os dois uma vez em , deixando para a primeira equação:

$$\begin{align} \frac{\partial^2 c}{\partial t^2} &= \frac{\partial^2}{\partial x^2}\frac{\partial c}{\partial t} + \frac{\partial b}{\partial t} \\ \frac{\partial b}{\partial t} &= - \left(\frac{\partial c}{\partial t}\right) b \end{align}$$ $$\begin{align} \frac{\partial}{\partial t}\left(\frac{\partial c}{\partial t} - \frac{\partial^2 c}{\partial x^2} - b \right) &= 0 \\ \frac{1}{b}\left(\frac{\partial b}{\partial t}\right) &= -\frac{\partial c}{\partial t} \end{align}$$ $t$ $$\frac{\partial c}{\partial t} - \frac{\partial^2 c}{\partial x^2} = b + A(x)$$ Usando a terceira equação, podemos expressar a "constante" da integração como . A segunda equação é um pouco mais complicada. Reescrevendo um pouco, temos: Isso leva à solução Ou exponenciamento fornece: E, finalmente, conectando-o ao PDE para dá $A(x)=a_0-\frac{\partial^2 c_0}{\partial x^2}-b_0$ $$\int_0^t \frac{1}{b(x,t')}\left(\frac{\partial b(x,t')}{\partial t'}\right) dt' = -\int_0^t \frac{\partial c(x,t')}{\partial t'} dt'$$ $$\ln b(x,t) - \ln b_0(x) = - c(x) + c_0(x)$$ $$\ln \frac{b}{b_0} = -c + c_0$$ $$b = b_0 e^{c_0-c}$$ $c$ $$\frac{\partial c}{\partial t} - \frac{\partial^2 c}{\partial x^2} = b_0 e^{c_0-c}+A(x)$$

Substituindo por , ou equivalente usando as condições iniciais , esta equação simplifica para Agora, você deve conseguir literatura considerável sobre a melhor maneira de resolver essa equação. Não sei se Crank-Nicholson é uma boa escolha para o termo exponencial, mas parece plausível. Alguns cuidados provavelmente devem ser tomados para garantir que em lugar ou a solução possa explodir rapidamente.$c$$c-c_0$$c_0(x)=0$

$$\frac{\partial c}{\partial t} = \frac{\partial^2 c}{\partial x^2} + a_0 - (1-e^{-c}) b_0$$ $c > c_0$

Eu só passei por essa derivação duas vezes, então pode haver um erro ou dois, mas tem a sensação certa para mim. Se todos os lugares, então esta é claramente a solução certa e, de outra forma, tem um cheiro de plausibilidade.$b_0 = 0$

Resolver isso até e avaliar deve fornecer a resposta que você está procurando.$t=1$$c(x,1)$