Sistema de reacción-difusión

format_list_bulleted Contenido keyboard_arrow_down

ImprimirCitar

Los sistemas de reacción-difusión son modelos matemáticos que corresponden a varios fenómenos físicos. El más común es el cambio en el espacio y el tiempo de la concentración de una o más sustancias químicas: reacciones químicas locales en las que las sustancias se transforman unas en otras, y difusión que hace que las sustancias se extiendan sobre una superficie en el espacio.

Los sistemas de reacción-difusión se aplican naturalmente en química. Sin embargo, el sistema también puede describir procesos dinámicos de naturaleza no química. Se encuentran ejemplos en biología, geología y física (teoría de la difusión de neutrones) y ecología. Matemáticamente, los sistemas de reacción-difusión toman la forma de ecuaciones diferenciales parciales parabólicas semilineales. Se pueden representar en la forma general $parcial _{t}{boldsymbol {q}}={underline {underline {boldsymbol {D}}}},nabla ^{2}{boldsymbol {q}}+{boldsymbol {R }}({ símbolo de negrita {q}}),$

donde q (x, t) representa la función vectorial desconocida, D es una matriz diagonal de coeficientes de difusión y R representa todas las reacciones locales. Las soluciones de las ecuaciones de reacción-difusión muestran una amplia gama de comportamientos, incluida la formación de ondas viajeras y fenómenos ondulatorios, así como otros patrones autoorganizados como rayas, hexágonos o estructuras más complejas como solitones disipativos. Estos patrones se han denominado "patrones de Turing". Cada función, para la que se cumple una ecuación diferencial de difusión de reacción, representa de hecho una variable de concentración.

Ecuaciones de reacción-difusión de un componente

La ecuación de reacción-difusión más simple está en una dimensión espacial en geometría plana, $parcial _{t}u=Dparcial _{x}^{2}u+R(u),$

también se conoce como la ecuación de Kolmogorov-Petrovsky-Piskunov. Si el término de reacción desaparece, entonces la ecuación representa un proceso de difusión puro. La ecuación correspondiente es la segunda ley de Fick. La elección R (u) = u (1 − u) produce la ecuación de Fisher que se usó originalmente para describir la propagación de poblaciones biológicas, la ecuación de Newell-Whitehead-Segel con R (u) = u (1 − u) para describir Rayleigh –Convección de Bénard, la ecuación más general de Zeldovich–Frank-Kamenetskii con R (u) = u(1 − u)e y 0 < β < ∞ (número de Zeldovich) que surge en la teoría de la combustión, y su caso degenerado particular con R (u) = u − u que a veces también se conoce como la ecuación de Zeldovich.

La dinámica de los sistemas de un componente está sujeta a ciertas restricciones ya que la ecuación de evolución también se puede escribir en la forma variacional $parcial _{t}u=-{frac {delta {mathfrak {L}}}{delta u}}$

y por lo tanto describe una disminución permanente de la "energía libre" ${ matemáticas {L}}$ dada por el funcional ${displaystyle {mathfrak {L}}=int_{-infty}^{infty}left[{tfrac {D}{2}}left(parcial_{x}uright) ^{2}-V(u)derecha],{texto{d}}x}$

con un potencial V (u) tal que R (u) =d V (tu)/eres tú.

En sistemas con más de una solución homogénea estacionaria, una solución típica viene dada por frentes móviles que conectan los estados homogéneos. Estas soluciones se mueven con velocidad constante sin cambiar su forma y son de la forma u (x, t) = û (ξ) con ξ = x − ct, donde c es la velocidad de la onda viajera. Tenga en cuenta que mientras que las ondas viajeras son estructuras genéricamente estables, todas las soluciones estacionarias no monótonas (por ejemplo, dominios localizados compuestos por un par frente-antifrente) son inestables. Para c = 0, hay una prueba simple para esta afirmación:si u ₀ (x) es una solución estacionaria y u = u ₀ (x) + ũ (x, t) es una solución infinitesimalmente perturbada, el análisis de estabilidad lineal produce la ecuación ${displaystyle parcial _{t}{tilde {u}}=Dparcial _{x}^{2}{tilde {u}}-U(x){tilde {u}},qquad U(x)=-R^{primo}(u){Grande |}_{u=u_{0}(x)}.}$

Con el ansatz ũ = ψ (x)exp(− λt) llegamos al problema de valores propios ${sombrero {H}}psi =lambda psi,qquad {sombrero {H}}=-Dparcial _{x}^{2}+U(x),$

de tipo Schrödinger donde los valores propios negativos resultan en la inestabilidad de la solución. Debido a la invariancia traslacional ψ = ∂ _x u ₀ (x) es una función propia neutra con el valor propio λ = 0, y todas las demás funciones propias se pueden ordenar de acuerdo con un número creciente de nodos con la magnitud del valor propio real correspondiente aumenta monótonamente con el número de ceros. La función propia ψ = ∂ _x u ₀ (x) debe tener al menos un cero, y para una solución estacionaria no monótona el valor propio correspondiente λ = 0no puede ser el más bajo, lo que implica inestabilidad.

Para determinar la velocidad c de un frente en movimiento, uno puede ir a un sistema de coordenadas en movimiento y buscar soluciones estacionarias: $Dparcial _{xi }^{2}{sombrero {u}}(xi)+cparcial _{xi }{sombrero {u}}(xi)+R({sombrero { u}}(xi))=0.$

Esta ecuación tiene un buen análogo mecánico como el movimiento de una masa D con posición û en el transcurso del "tiempo" ξ bajo la fuerza R con el coeficiente de amortiguamiento c que permite un acceso bastante ilustrativo a la construcción de diferentes tipos de soluciones. y la determinación de c.

Al pasar de una a más dimensiones espaciales, aún se pueden aplicar una serie de declaraciones de sistemas unidimensionales. Los frentes de onda planos o curvos son estructuras típicas, y surge un nuevo efecto cuando la velocidad local de un frente curvo se vuelve dependiente del radio local de curvatura (esto se puede ver yendo a coordenadas polares). Este fenómeno conduce a la llamada inestabilidad impulsada por la curvatura.

Ecuaciones de reacción-difusión de dos componentes

Los sistemas de dos componentes permiten una gama mucho más amplia de fenómenos posibles que sus contrapartes de un componente. Una idea importante propuesta por primera vez por Alan Turing es que un estado que es estable en el sistema local puede volverse inestable en presencia de difusión.

Sin embargo, un análisis de estabilidad lineal muestra que al linealizar el sistema general de dos componentes ${displaystyle {begin{pmatrix}parcial _{t}u\parcial _{t}vend{pmatrix}}={begin{pmatrix}D_{u}&0\0&D_{v} end{pmatrix}}{begin{pmatrix}partial _{xx}u\partial _{xx}vend{pmatrix}}+{begin{pmatrix}F(u,v)\G(u,v)end{matrix}}}$

una perturbación de onda plana ${tilde {boldsymbol {q}}}_{boldsymbol {k}}({boldsymbol {x}},t)={begin{pmatrix}{tilde {u}}(t)\{ tilde {v}}(t)end{pmatrix}}e^{i{boldsymbol {k}}cdot {boldsymbol {x}}}$

de la solución homogénea estacionaria satisfará ${begin{pmatrix}parcial _{t}{tilde {u}}_{boldsymbol {k}}(t)\parcial _{t}{tilde {v}}_{boldsymbol { k}}(t)end{pmatrix}}=-k^{2}{begin{pmatrix}D_{u}{tilde {u}}_{boldsymbol {k}}(t)\D_ {v}{tilde {v}}_{boldsymbol {k}}(t)end{pmatrix}}+{boldsymbol {R}}^{prime }{begin{pmatrix}{tilde { u}}_{boldsymbol {k}}(t)\{tilde {v}}_{boldsymbol {k}}(t)end{pmatrix}}.$

La idea de Turing solo se puede realizar en cuatro clases de equivalencia de sistemas caracterizados por los signos del jacobiano R ′ de la función de reacción. En particular, si se supone que un vector de onda finito k es el más inestable, el jacobiano debe tener los signos ${begin{pmatrix}+&-\+&-end{pmatrix}},quad {begin{pmatrix}+&+\-&-end{pmatrix}},quad {begin{ pmatrix}-&+\-&+end{pmatrix}},quad {begin{pmatrix}-&-\+&+end{pmatrix}}.$

Esta clase de sistemas se denomina sistema activador-inhibidor por su primer representante: cerca del estado fundamental, un componente estimula la producción de ambos componentes mientras que el otro inhibe su crecimiento. Su representante más destacado es la ecuación de FitzHugh-Nagumo ${begin{alineado}parcial_{t}u&=d_{u}^{2},nabla ^{2}u+f(u)-sigma v,\tau parcial_{t }v&=d_{v}^{2},nabla ^{2}v+uvend{alineado}}$

con f (u) = λu − u − κ que describe cómo viaja un potencial de acción a través de un nervio. Aquí, d _u, d _v, τ, σ y λ son constantes positivas.

Cuando un sistema activador-inhibidor sufre un cambio de parámetros, se puede pasar de condiciones en las que un estado fundamental homogéneo es estable a condiciones en las que es linealmente inestable. La bifurcación correspondiente puede ser una bifurcación de Hopf a un estado homogéneo de oscilación global con un número de onda dominante k = 0 o una bifurcación de Turing a un estado de patrón global con un número de onda finito dominante. Este último en dos dimensiones espaciales típicamente conduce a patrones de rayas o hexagonales.

Bifurcación de Turing subcrítica: formación de un patrón hexagonal a partir de condiciones iniciales ruidosas en el sistema de reacción-difusión de dos componentes anterior del tipo Fitzhugh-Nagumo.
Condiciones iniciales ruidosas en t = 0.
Estado del sistema en t = 10.
Estado casi convergente en t = 100.

Para el ejemplo de Fitzhugh-Nagumo, las curvas de estabilidad neutra que marcan el límite de la región linealmente estable para la bifurcación de Turing y Hopf están dadas por ${begin{alineado}q_{text{n}}^{H}(k):&{}quad {frac {1}{tau }}+left(d_{u}^{2} +{frac {1}{tau }}d_{v}^{2}right)k^{2}&=f^{prime }(u_{h}),\[6pt]q_{ text{n}}^{T}(k):&{}quad {frac {kappa }{1+d_{v}^{2}k^{2}}}+d_{u}^ {2}k^{2}&=f^{primo}(u_{h}).end{alineado}}$

Si la bifurcación es subcrítica, a menudo se pueden observar estructuras localizadas (solitones disipativos) en la región histérica donde el patrón coexiste con el estado fundamental. Otras estructuras que se encuentran con frecuencia comprenden trenes de pulsos (también conocidos como ondas viajeras periódicas), ondas espirales y patrones objetivo. Estos tres tipos de solución también son características genéricas de ecuaciones de reacción-difusión de dos (o más) componentes en las que la dinámica local tiene un ciclo límite estable

Otros patrones encontrados en el sistema de reacción-difusión de dos componentes anterior del tipo Fitzhugh-Nagumo.
Espiral giratoria.
Patrón objetivo.
Pulso estacionario localizado (solitón disipativo).

Ecuaciones de reacción-difusión de tres o más componentes

Para una variedad de sistemas, se han propuesto ecuaciones de reacción-difusión con más de dos componentes, por ejemplo, la reacción de Belousov-Zhabotinsky, para coagulación sanguínea, ondas de fisión o sistemas de descarga plana de gas.

Se sabe que los sistemas con más componentes permiten una variedad de fenómenos que no son posibles en sistemas con uno o dos componentes (por ejemplo, pulsos de funcionamiento estables en más de una dimensión espacial sin retroalimentación global). Se proporciona una introducción y una descripción sistemática de los posibles fenómenos en función de las propiedades del sistema subyacente.

Aplicaciones y universalidad

En tiempos recientes, los sistemas de reacción-difusión han atraído mucho interés como modelo prototipo para la formación de patrones. Los patrones mencionados anteriormente (frentes, espirales, objetivos, hexágonos, rayas y solitones disipativos) se pueden encontrar en varios tipos de sistemas de reacción-difusión a pesar de las grandes discrepancias, por ejemplo, en los términos de reacción locales. También se ha argumentado que los procesos de reacción-difusión son una base esencial para los procesos relacionados con la morfogénesis en biología e incluso pueden estar relacionados con el pelaje de los animales y la pigmentación de la piel. Otras aplicaciones de las ecuaciones de reacción-difusión incluyen invasiones ecológicas, propagación de epidemias, crecimiento tumoral, dinámica de ondas de fisión y cicatrización de heridas.Otra razón del interés en los sistemas de reacción-difusión es que, aunque son ecuaciones diferenciales parciales no lineales, a menudo hay posibilidades de un tratamiento analítico.

Experimentos

Los experimentos bien controlables en sistemas de reacción-difusión química se han realizado hasta ahora de tres maneras. En primer lugar, se pueden utilizar reactores de gel o tubos capilares llenos. En segundo lugar, se han investigado los pulsos de temperatura en las superficies catalíticas. En tercer lugar, la propagación de pulsos nerviosos en ejecución se modela utilizando sistemas de reacción-difusión.

Aparte de estos ejemplos genéricos, ha resultado que, en circunstancias apropiadas, los sistemas de transporte eléctrico como plasmas o semiconductores pueden describirse en un enfoque de reacción-difusión. Para estos sistemas se han llevado a cabo varios experimentos de formación de patrones.

Tratamientos numéricos

Un sistema de reacción-difusión se puede resolver usando métodos de matemáticas numéricas. Existen varios tratamientos numéricos en la literatura de investigación. También para geometrías complejas se proponen métodos de solución numérica. Los sistemas de reacción-difusión se describen con el mayor grado de detalle con herramientas de simulación basadas en partículas como SRSim o ReadDy que emplean, por ejemplo, dinámicas de reacción reversibles de partículas que interactúan.

Contenido relacionado

Más resultados...