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

De Física Computacional
Ir para navegação Ir para pesquisar
Sem resumo de edição
Sem resumo de edição
Linha 21: Linha 21:


== Análise de estabilidade ==
== Análise de estabilidade ==
Em geral, a análise de estabilidade linear de sistemas reativos-difusivos se constitui em duas etapas. Inicialmente se determinam os estados de equilíbrio homogêneo estáveis (igualando os coeficientes de difusão a zero). A seguir se analisa se a inclusão dos coeficientes de difusão gera instabilidade naqueles estados.
A difusão normalmente tende a homogeneizar o sistema. O surgimento de padrões complexos e não homogêneos ocorre nos casos de exceção a essa regra, quando a difusão torna os estados de equilíbrio instáveis. <ref name=Biologia>http://mcb111.org/w13/w13-lecture.html#the-gray-scott-model</ref>


'''Nota''': A análise em toda esta seção pressupõe sempre que os parâmetros e coeficientes de difusão são positivos.
'''Nota''': A análise em toda esta seção pressupõe sempre que os parâmetros e coeficientes de difusão são positivos.
Linha 60: Linha 56:
Ambas desigualdades são imediatamente satisfeitas para quaisquer valores de <math>F, k, D_{u}</math>, e <math>D_{v}</math>. Portanto, o estado de equilíbrio <math>(u^{*}, v^{*}) = (1, 0)</math> permanece estável no modelo de Gray-Scott mesmo após a inclusão dos coeficientes de difusão, sejam quais forem os valores desses coeficientes (lembrando que estamos nos restringindo a valores positivos dos parâmetros e coeficientes).
Ambas desigualdades são imediatamente satisfeitas para quaisquer valores de <math>F, k, D_{u}</math>, e <math>D_{v}</math>. Portanto, o estado de equilíbrio <math>(u^{*}, v^{*}) = (1, 0)</math> permanece estável no modelo de Gray-Scott mesmo após a inclusão dos coeficientes de difusão, sejam quais forem os valores desses coeficientes (lembrando que estamos nos restringindo a valores positivos dos parâmetros e coeficientes).


Esse é um resultado à primeira vista surpreendente. Em geral, o surgimento de padrões complexos e não homogêneos em sistemas reativos-difusivos está relacionado à desestabilização de um ou mais estados de equilíbrio homogêneo causada pela introdução dos coeficientes de difusão (conhecida como instabilidade de Turing)<ref name=Biologia>http://mcb111.org/w13/w13-lecture.html#the-gray-scott-model</ref>.


=== Estados de Equilíbrio Não Triviais ===
Entretanto, no caso do modelo de Gray-Scott, o surgimento de padrões complexos e não homogêneos '''não''' decorre da instabilidade de Turing, uma vez que o surgimento de padrões não triviais nesse modelo ocorre mesmo quando apenas o estado de equilíbrio trivial <math>(u^{*}, v^{*}) = (1, 0)</math> está presente <ref name=Gros>C. Gros, "Complex and Adaptive Dynamical Systems". Springer-Verlag, Berlim, 2015.</ref>.




Como o modelo claramente exibe padrões complexos não homogêneos 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 estados de equilíbrio não triviais, soluções do sistema de equações (1).<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>
=== Estados de Equilíbrio Não Triviais ===




Os outros estados de equilíbrio do modelo de Gray-Scott, com a restrição <math>F\geq4(F+k)^2</math>, são <math>(u_{+},v_{-})</math> e <math>(u_{-},v_{+})</math>, com  
outros dois estados de equilíbrio que são soluções não triviais do sistema de equações (1). Desde que seja obedecida a condição <math>F\geq4(F+k)^2</math>, esses estados são <math>(u_{+},v_{-})</math> e <math>(u_{-},v_{+})</math>, com<ref name=Gros></ref>





Edição das 21h19min de 20 de fevereiro de 2022

Introdução

Descrição do Modelo

The reaction-diffusion system described here involves two generic chemical species U and V, whose concentration at a given point in space is referred to by variables u and v. As the term implies, they react with each other, and they diffuse through the medium. Therefore the concentration of U and V at any given location changes with time and can differ from that at other locations.

There are two reactions which occur at different rates throughout the space according to the relative concentrations at each point:

U + 2V → 3V

V → P

P is an inert product. It is assumed for simplicity that the reverse reactions do not occur (this is a useful simplification when a constant supply of reagents8 prevents the attainment of equilibrium). Because V appears on both sides of the first reaction, it acts as a catalyst for its own production.

The overall behavior of the system is described by the following formula, two equations which describe three sources of increase and decrease for each of the two chemicals:


Análise de estabilidade

Nota: A análise em toda esta seção pressupõe sempre que os parâmetros e coeficientes de difusão são positivos.


Estado de Equilíbrio Trivial

O modelo de Gray-Scott depende dos parâmetros e dos coeficientes de difusão . É fácil mostrar que, ignorando os termos de difusão, o sistema possui estado de equilíbrio estável em para quaisquer valores dos parâmetros.

Demonstração. O sistema de equações do modelo, com e , fazendo , é dado por



Logo, é trivial que o sistema acima é satisfeito quando . Esse estado de equilíbrio é estável porque a matriz jacobiana possui traço negativo e determinante positivo[1].

Se agora incluímos os termos de difusão e , deve-se levar em consideração a matriz . Aqui, é a matriz jacobiana dos termos de reação, é 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 :



Para que o estado de equilíbrio 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



Ambas desigualdades são imediatamente satisfeitas para quaisquer valores de , e . Portanto, o estado de equilíbrio permanece estável no modelo de Gray-Scott mesmo após a inclusão dos coeficientes de difusão, sejam quais forem os valores desses coeficientes (lembrando que estamos nos restringindo a valores positivos dos parâmetros e coeficientes).

Esse é um resultado à primeira vista surpreendente. Em geral, o surgimento de padrões complexos e não homogêneos em sistemas reativos-difusivos está relacionado à desestabilização de um ou mais estados de equilíbrio homogêneo causada pela introdução dos coeficientes de difusão (conhecida como instabilidade de Turing)[2].

Entretanto, no caso do modelo de Gray-Scott, o surgimento de padrões complexos e não homogêneos não decorre da instabilidade de Turing, uma vez que o surgimento de padrões não triviais nesse modelo ocorre mesmo quando apenas o estado de equilíbrio trivial está presente [3].


Estados de Equilíbrio Não Triviais

Há outros dois estados de equilíbrio que são soluções não triviais do sistema de equações (1). Desde que seja obedecida a condição , esses estados são e , com[3]



Referências

  1. 1,0 1,1 H. Sayama, "Introduction to the Modeling and Analysis of Complex Systems". Open SUNY Textbooks, Geneseo, NY, 2015.
  2. http://mcb111.org/w13/w13-lecture.html#the-gray-scott-model
  3. 3,0 3,1 C. Gros, "Complex and Adaptive Dynamical Systems". Springer-Verlag, Berlim, 2015.