Erraten, was eine einfache partielle Differentialgleichung physikalisch beschreibt

Gibt es eine einfache Möglichkeit, eine teilweise andere Gleichung zu betrachten und ein Gefühl dafür zu bekommen, welche Art von Phänomen sie physikalisch beschreibt? Ich habe eine Gleichung, die so aussieht:

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

Nach der Untersuchung und einigen Suchen schien es, als würde es versuchen, ein Phänomen zu beschreiben, bei dem Advektion und Diffusion gleichzeitig auftreten. Ist dies eine korrekte Interpretation?

Beachten Sie, dass C 3 > 0 , C 2 < 0 , C 1 < 0 , C 0 > 0 Und v ( X , T ) ist Geschwindigkeit. Ich bin mir nicht sicher, aber es scheint vernünftig anzunehmen, dass die Geschwindigkeit entweder negativ oder positiv sein kann. Auch, A ( X , T ) sollte positiv sein.

Beachten Sie, dass ich nicht nach einer Lösung für die Gleichung suche, sondern nur nach der physikalischen Interpretation, welche Art von Physik beteiligt wäre.

Konnten Sie nachsehen, ob es in Polyanins Büchern steht ?
Was ist mit dem Vorzeichen der Konstanten? A , v , C ich ?
Im Allgemeinen sollten Fragen nicht quer gepostet werden, es sei denn, Sie haben versucht und es nicht geschafft, auf einer Seite nützliche Beiträge zu erhalten. Also werden wir wahrscheinlich entweder diese oder die auf der Matheseite schließen.
@David Zaslavsky: Danke, dass du mich informiert hast. Ich verstehe. Crossposting werde ich in Zukunft vermeiden.

Antworten (4)

Versuchen wir, die Gleichung in ungefährer Form endlicher Differenzen umzuschreiben :

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
Wo Δ T -- ist ein Zeitschritt, und H - Leerschritt.
Der Ausdruck wird am Ende zu Ihrem PDE Δ T 0 , H 0 .

Die linke Seite beschreibt, wie hoch die Menge ist A Änderungen während eines Zeitschritts an einem bestimmten Punkt (ich hoffe, das ist offensichtlich). Lassen Sie uns Term für Term sehen, was auf der rechten Seite steht:

  • C 3 : Wenn A ist an "benachbarten" Punkten insgesamt größer als im Punkt selbst: A ( X + H , T ) + A ( X H , T ) > 2 A ( X , T ) Dann A wird bei zunehmen X . Sonst wird es weniger. Also der Begriff "Kräfte" A zu einer Art lokalem Gleichgewicht.
    Dies ist ein Standardbegriff zur Beschreibung verschiedener Diffusions- oder Wärmeverteilungsprozesse , die Ihr System zu einem gewissen Gleichgewicht zwingen.
  • C 2 : Hier muss man mit Schildern aufpassen. Lassen Sie uns anrufen X + H "der Punkt nach rechts" und X H wird "der Punkt nach links". Wenn v > 0 dann "fließt es nach rechts" und wenn v < 0 es "fließt nach links". Jetzt überprüfen Sie, ob Ihr Fluss in Ihren Punkt geht X (sagen wir, fließt nach rechts vom Punkt nach links), dann steigt es an A bei X . Und nimmt ab, wenn der Fluss von Ihrem Punkt weggeht.
    Das ist ein Standardbegriff, der Dinge beschreibt, die mit dem Strom getragen werden v und normalerweise in Form von wesentlichen Derivaten abgeleitet
  • C 1 Und C 0 : Diese beiden sind trivial, weil sie lokal sind. C 0 gibt nur einen konstanten Beitrag zum Wachstum von A , Und C 1 ändert (in Ihrem Fall hemmt) die Wachstumsrate proportional zum Wert von A .
    Sie können sie auch als lokale ODE für einen bestimmten Punkt verstehen: A ˙ = C 1 A + C 0

Um es zusammenzufassen: Lassen Sie uns das sagen A ( X , T ) beschreibt die Dichte von beispielsweise Bakterien in einem Röhrchen. Dann C 3 beschreibt, wie sie sich ausbreiten, C 2 beschreibt, wie sie von einem Strom getragen werden v ( X , T ) Flüssigkeit in der Röhre, C 0 -- ist eine Wachstumsrate neuer Bakterien und schließlich C 1 -- ist dafür verantwortlich, dass sich das Wachstum durch die Überbevölkerung verlangsamt.

Dies ist eine unglaublich gute und nützliche Antwort. Ich wünschte, ich hätte noch eine Stimme.
In der Bakterieninterpretation C 0 ist wirklich eine Einwanderungsrate, keine Wachstumsrate: Neue Bakterien erscheinen, auch wenn anfänglich keine vorhanden waren. C 1 , kann aber durchaus als dichteunabhängige Sterblichkeitsrate (bzw. Auswanderungsrate) pro Kopf interpretiert werden.

Als Alternative zur Wärmeinterpretation von Christian Blatter A könnte die Konzentration von Partikeln beschreiben, die auf einer eindimensionalen Substratoberfläche (oder einer zweidimensionalen, bei der wir eine der Dimensionen ignorieren) adsorbiert sind.

  • Neue Partikel werden mit einer Rate adsorbiert C 0 pro Längeneinheit.
  • Adsorbierte Partikel lösen sich schnell von der Oberfläche C 1 pro Teilchen.
  • Die Partikel bewegen sich entlang der Oberfläche (oder die Oberfläche selbst bewegt sich) mit mittlerer Nettogeschwindigkeit C 2 v .
  • Während der Bewegung diffundieren die Partikel auch über die Oberfläche mit Diffusionskoeffizient C 3 .

In jedem Fall beschreibt diese Gleichung die Dynamik einer Größe A Diffusion und Advektion über einen eindimensionalen Raum erfahren, während sie auch einer konstanten (dh nullten Ordnung) lokalen Akkumulation und einem Zerfall erster Ordnung unterliegen.

(Als mathematischer Ökologe war mein erster Gedanke, es als räumliches Populationsmodell zu interpretieren, aber es passt nicht wirklich gut zu dieser Interpretation: Es gibt kein A 2 Begriff, der die lokale Dichteregulierung beschreiben könnte.)

Das ist eine faszinierende Interpretation. Wo finde ich mehr zu dieser Betrachtungsweise?

Sie haben ein dünnes zylindrisches Rohr entlang der X -Achse, die mit etwas Dichtegas gefüllt ist ρ ( X , T ) . Die Temperatur des Gases ist A ( X , T ) , und das Gas bewegt sich zusammen mit dem "Massenstrom" M ( X , T ) := ρ ( X , T ) v ( X , T ) , Wo v bezeichnet die tatsächliche Geschwindigkeit einzelner Teilchen. (Der ρ fehlt in deiner Gleichung). Wärme wird durch Wärmeleitung und mittels Konvektion transportiert. Zusätzlich X -Achse ist ein elektrischer Draht, der Wärme mit konstanter Rate erzeugt, und an der Oberfläche der Röhre haben wir einen Wärmeverlust in Richtung des Weltraums, wobei letzterer eine Temperatur ist 0 .

Ihre Gleichung beschreibt die zeitliche Änderungsrate der Temperatur in einem „Längenelement“ an X zum Zeitpunkt T . Die einzelnen Terme auf der rechten Seite berücksichtigen die Beiträge von Leitung, Konvektion, Oberflächenverlust zum Weltraum und Erwärmung.

Das Vernachlässigen von Termen und das Lösen der Gleichung für idealisierte Anfangsbedingungen ist eine Möglichkeit, die Bedeutung der einzelnen Terme zu untersuchen.

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

Legen Sie beispielsweise alle außer fest C 0 auf Null und erhalten A = C 0 T + A ( X , 0 ) . Der Begriff mit C 0 liefert die Menge dargestellt durch A mit konstanter Rate. Es ist ein Quellbegriff.

Stellen Sie alle aber ein C 1 auf Null und erhalten A = A ( X , 0 ) e C 1 T . Dies ist eine Art nichtlinearer Quellterm -- A wird in einer positiven Rückkopplungsrate proportional zu seiner Menge zugeführt.

Stellen Sie alle aber ein C 2 auf Null und Sie finden A = A ( X , 0 ) δ ( X C 2 v T ) vorausgesetzt v ist konstant. Sie können für verallgemeinern v x-abhängig. Dies ist ein Advektionsbegriff, der die Bewegung beschreibt A bei Geschwindigkeit C 2 v .

Stellen Sie alle aber ein C 3 auf Null und Sie finden die klassische Diffusionsgleichung. Hier ist die Lösung für eine Punktquellen-Anfangsbedingung etwa so A ( X , T ) = 1 C 3 T e X 2 / ( C 3 T ) obwohl ich einige Konstanten falsch habe. Dies beschreibt eine Tendenz zum "Ausbreiten" oder Verschmieren aufgrund unterschiedlicher Transporteigenschaften zwischen Komponenten von A .

Um es kurz zu machen, Sie haben Diffusion mit Diffusivität C 3 , Advektion mit Geschwindigkeit C 2 v , nichtlineare Versorgung mit Rate C 1 , und konstante Versorgung mit Rate C 0 .