Adivinando lo que una simple ecuación diferencial parcial describe físicamente

¿Hay una manera fácil de mirar una ecuación parcial diferente y tener una idea de qué tipo de fenómenos está describiendo físicamente? Tengo una ecuación que se ve así:

t A = C 3 X X 2 A + C 2 X ( v A ) + C 1 A + C 0 .

Al inspeccionarlo y hacer algunas búsquedas, parecía que podría estar tratando de describir un fenómeno que tenía advección y difusión al mismo tiempo. ¿Es esta una interpretación correcta?

Tenga en cuenta que C 3 > 0 , C 2 < 0 , C 1 < 0 , C 0 > 0 y v ( X , t ) es la velocidad. No estoy seguro, pero parece razonable suponer que la velocidad puede ser negativa o positiva. También, A ( X , t ) debería ser positivo.

Tenga en cuenta que no estoy buscando una solución a la ecuación, solo la interpretación física de qué tipo de física estaría involucrada.

¿ Has conseguido comprobar si está en los libros de Polyanin ?
¿Qué pasa con el signo de las constantes? A , v , C i ?
En general, no se supone que las preguntas se publiquen de forma cruzada, a menos que haya intentado y no haya podido obtener información útil en un sitio. Así que probablemente terminaremos cerrando este o el del sitio de matemáticas.
@David Zaslavsky: Gracias por avisarme. Entiendo. Evitaré las publicaciones cruzadas en el futuro.

Respuestas (4)

Intentemos reescribir la ecuación en forma aproximada de diferencias finitas :

A ( X , t + Δ t ) A ( X , t ) Δ t = C 3 A ( X + h , t ) + A ( X h , t ) 2 A ( X , t ) h 2 +
+ C 2 v ( X + h , t ) A ( X + h , t ) v ( X h , t ) A ( X h , t ) 2 h + C 1 A ( X , t ) + C 0
Dónde Δ t -- es un paso de tiempo, y h -- paso espacial.
La expresión se convierte en tu PDE, en el límite Δ t 0 , h 0 .

El lado izquierdo describe cuánto la cantidad A cambia durante un paso de tiempo en un punto dado (espero que sea obvio). Veamos qué hay en el lado derecho término por término:

  • C 3 : Si A es mayor en total en los puntos "vecinos" que en el punto mismo: A ( X + h , t ) + A ( X h , t ) > 2 A ( X , t ) entonces A aumentará en X . De lo contrario, disminuirá. Así que el término "fuerzas" A a algún tipo de equilibrio local.
    Este es un término estándar para la descripción de varios procesos de difusión o distribución de calor que obligan a su sistema a alcanzar cierto equilibrio.
  • C 2 : Aquí hay que tener cuidado con las señales. llamemos X + h "el punto a la derecha" y X h será "el punto a la izquierda". Si v > 0 entonces "fluye hacia la derecha" y si v < 0 "fluye hacia la izquierda". Ahora compruebas que si tu flujo entra en tu punto X (digamos, fluye hacia la derecha desde el punto hacia la izquierda) luego aumenta A en X . Y disminuye si los flujos se alejan de su punto.
    Ese es un término estándar que describe cosas llevadas con la corriente. v y generalmente derivado en términos de derivado sustancial
  • C 1 y C 0 : Estos dos son triviales, porque son locales. C 0 solo da una contribución constante al crecimiento de A , y C 1 cambia (inhibe en su caso) la tasa de crecimiento proporcional al valor de A .
    También puede entenderlos en términos de EDO local para un punto dado: a ˙ = C 1 a + C 0

Para resumir: digamos que A ( X , t ) describe la densidad de, digamos, bacterias en un tubo. Entonces C 3 describe cómo se difunden, C 2 describe cómo son llevados por una corriente v ( X , t ) de líquido en el tubo, C 0 -- es una tasa de crecimiento de nuevas bacterias y, finalmente, C 1 -- es responsable de que el crecimiento se ralentice debido a la sobrepoblación.

Esta es una respuesta increíblemente buena y útil. Ojalá tuviera otro voto.
En la interpretación de las bacterias, C 0 es realmente una tasa de inmigración, no una tasa de crecimiento: aparecen nuevas bacterias incluso si no había ninguna presente inicialmente. C 1 Sin embargo, sí puede interpretarse como una tasa de mortalidad (y/o emigración) per cápita independiente de la densidad.

Como alternativa a la interpretación del calor de Christian Blatter, A podría describir la concentración de partículas adsorbidas en una superficie de sustrato unidimensional (o bidimensional, donde ignoramos una de las dimensiones).

  • Las partículas nuevas se adsorben a una velocidad C 0 por unidad de longitud.
  • Las partículas adsorbidas se separan de la superficie a una velocidad C 1 por partícula.
  • Las partículas se mueven a lo largo de la superficie (o la superficie misma se mueve) con una velocidad neta media C 2 v .
  • Mientras se mueven, las partículas también se difunden sobre la superficie con coeficiente de difusión C 3 .

En cualquier caso, esta ecuación describe la dinámica de alguna cantidad A experimentando difusión y advección sobre un espacio unidimensional, mientras que también experimenta acumulación local constante (es decir, de orden cero) y decaimiento de primer orden.

(Como ecólogo matemático, mi primer pensamiento fue interpretarlo como un modelo de población espacial, pero en realidad no se ajusta muy bien a esa interpretación: no hay A 2 término que podría describir la regulación de la densidad local.)

Esta es una interpretación intrigante. ¿Dónde puedo encontrar más sobre esta forma de verlo?

Tienes un tubo cilíndrico delgado a lo largo del X -eje que está lleno de algún gas de densidad ρ ( X , t ) . La temperatura del gas es A ( X , t ) , y el gas se mueve junto con el "flujo de masa" metro ( X , t ) := ρ ( X , t ) v ( X , t ) , dónde v denota la velocidad real de las partículas individuales. (El ρ falta en su ecuación). El calor se transporta por conducción de calor y por convección. Además, el X -eje es un cable eléctrico que produce calor a un ritmo constante, y en la superficie del tubo tenemos una pérdida de calor hacia el espacio exterior, siendo este último una temperatura 0 .

Su ecuación describe la tasa de cambio temporal de la temperatura en un "elemento de longitud" en X en el momento t . Los términos individuales del lado derecho representan las contribuciones de conducción, convección, pérdida de superficie al espacio exterior y calentamiento.

Despreciar términos y resolver la ecuación para condiciones iniciales idealizadas es una forma de estudiar lo que significa cada término.

t A = C 3 X 2 A + C 2 X ( v A ) + C 1 A + C 0

Por ejemplo, configure todos menos C 0 a cero y obtener A = C 0 t + A ( X , 0 ) . El término con C 0 suministra la cantidad representada por A a un ritmo constante. Es un término fuente.

Establecer todo menos C 1 a cero y obtener A = A ( X , 0 ) mi C 1 t . Este es un tipo de término fuente no lineal: A se suministra a un ritmo proporcional a su cantidad en una realimentación positiva.

Establecer todo menos C 2 a cero y encuentras A = A ( X , 0 ) d ( X C 2 v t ) asumiendo v es constante Puedes generalizar para v dependiente de x. Este es un término de advección que describe el movimiento de A a velocidad C 2 v .

Establecer todo menos C 3 a cero y encuentras la clásica ecuación de difusión. Aquí la solución para una condición inicial de fuente puntual es algo así como A ( X , t ) = 1 C 3 t mi X 2 / ( C 3 t ) aunque tengo algunas constantes mal. Esto describe una tendencia a la "extensión" o manchado debido a las características de transporte diferencial entre los componentes de A .

Para resumir, tienes difusión con difusividad C 3 , advección con velocidad C 2 v , suministro no lineal a tasa C 1 , y suministro constante a tasa C 0 .