Modelo de Gray-Scott: mudanças entre as edições

De Física Computacional
Ir para navegação Ir para pesquisar
Wallec (discussão | contribs)
Wallec (discussão | contribs)
Linha 32: Linha 32:
<math>(F+D_{u}\omega^2)(F+k+D_{v}\omega^2) > 0</math>
<math>(F+D_{u}\omega^2)(F+k+D_{v}\omega^2) > 0</math>


<math>-2F - (D_{u}+D_{v}\omega^2) - k < 0</math>
<math>-2F - (D_{u}+D_{v})\omega^2 - k < 0</math>




Ambas desigualdades são imediatamente satisfeitas para quaisquer valores positivos de <math>F, k, D_{u}</math>, e <math>D_{v}</math>. Portanto, o estado de equilíbrio <math>(u^{*}, v^{*}) = (1, 0)</math> é estável no modelo de Gray-Scott nessas condições.
Ambas desigualdades são imediatamente satisfeitas para quaisquer valores positivos de <math>F, k, D_{u}</math>, e <math>D_{v}</math>. Portanto, o estado de equilíbrio <math>(u^{*}, v^{*}) = (1, 0)</math> é estável no modelo de Gray-Scott nessas condições.


Como o modelo claramente exibe padrões complexos para diversos valores positivos de <math>F, k, D_{u}</math>, e <math>D_{v}</math>, devem então existir outros estados de equilíbrio. De fato, há outros dois.
Como o modelo claramente exibe padrões complexos para diversos valores positivos de <math>F, k, D_{u}</math>, e <math>D_{v}</math>, devem então existir outros estados de equilíbrio e pelo menos um deles deve ser instável para determinados valores dos parâmetros. De fato, há outros dois estados de equilíbrio.<ref name=Wang>Tingting Wang et al., "Fractional Gray–Scott model: Well-posedness, discretization, and simulations, Computer Methods in Applied Mechanics and Engineering", Volume 347, 2019, pp. 1030-1049, ISSN 0045-7825, https://doi.org/10.1016/j.cma.2019.01.002.</ref>
 
Os outros estados de equilíbrio do modelo de Gray-Scott, ditos não triviais, são <math>(u_{+},v_{-})</math> e <math>(u_{-},v_{+})</math>, com
 
 
<math>u_{\pm} = \frac{1}{2}(1 \pm \sqrt{1 − 4\gamma^2F})</math>
 
 
 
<Reference\>


==Teste==
==Teste==

Edição das 17h09min de 19 de fevereiro de 2022

Análise de estabilidade

O modelo de Gray-Scott depende dos parâmetros F,k e dos coeficientes de difusão Du,Dv. É fácil mostrar que, ignorando os termos de difusão, o sistema possui estado de equilíbrio em (u*,v*)=(1,0).

Afirmação: O modelo de Gray-Scott possui estado de equilíbrio homogêneo em (u*,v*)=(1,0).

Demonstração: Resolvendo o sistema de equações do modelo com ut=0 e vt=0 e fazendo Du=Dv=0, temos


0=uv2+F(1u)

0=uv2(F+k)v


É então trivial que o sistema acima é satisfeito quando (u*,v*)=(1,0).

Esse estado de equilíbrio é estável no modelo de Gray-Scott para quaisquer valores positivos de F,k,Du, e Dv. É possível mostrar isso a partir da determinação dos autovalores da matriz


(JDω2)|f=feq


Aqui, J é a matriz jacobiana dos termos de reação, D é a matriz diagonal dos termos de difusão e ω é o parâmetro que determina a frequência espacial das perturbações. A demonstração da validade desse método pode ser encontrada na referência[1]. Aplicando ao modelo de Gray-Scott em (u*,v*)=(1,0)


((Fv22uvv2Fk+2uv)(Du00Dv)ω2)|(u,v)=(1,0)=(FDuω200FkDvω2)


Para que o estado de equilíbrio (u*,v*)=(1,0) seja estável é necessário que o determinante da matriz acima seja positivo e o seu traço seja negativo. Obtém-se então


(F+Duω2)(F+k+Dvω2)>0

2F(Du+Dv)ω2k<0


Ambas desigualdades são imediatamente satisfeitas para quaisquer valores positivos de F,k,Du, e Dv. Portanto, o estado de equilíbrio (u*,v*)=(1,0) é estável no modelo de Gray-Scott nessas condições.

Como o modelo claramente exibe padrões complexos para diversos valores positivos de F,k,Du, e Dv, devem então existir outros estados de equilíbrio e pelo menos um deles deve ser instável para determinados valores dos parâmetros. De fato, há outros dois estados de equilíbrio.[2]

Os outros estados de equilíbrio do modelo de Gray-Scott, ditos não triviais, são (u+,v) e (u,v+), com



<Reference\>

Teste

  1. H. Sayama, "Introduction to the Modeling and Analysis of Complex Systems", p. 287. Open SUNY Textbooks, Geneseo, NY, 2015.
  2. Tingting Wang et al., "Fractional Gray–Scott model: Well-posedness, discretization, and simulations, Computer Methods in Applied Mechanics and Engineering", Volume 347, 2019, pp. 1030-1049, ISSN 0045-7825, https://doi.org/10.1016/j.cma.2019.01.002.