COMSOL Multiphysics TutorialNonfrozen Nagumo-Equation (1D)

Denny Otten

Universitt Bielefeld

Fakultt fr Mathematik

Einleitung (Mathematischer Hintergrund)

Das Ziel dieses Tutorials ist die Implementierung der nicht-eingefrorenen und rumlich eindimensionalenNagumo-Gleichung mit Hilfe von COMSOL Multiphysics. Die Namensgebung der Gleichung ist auf denjapanischen Ingenieur Jin-Ichi Nagumo (1926-1999) zurckzufhren. Wir betrachten das folgende Anfangs-Randwertproblem

ut(x, t)u(x, t) = u(x, t) (1 u(x, t)) (u(x, t) a) , |x| < 50, t ]0,[

u

n(x, t) = 0 , |x| = 50, t [0,[

u(x, 0) = u0(x) , |x| < 50, t = 0

Hierbei bezeichnet a ]0, 12[ einen reellen Kontrollparameter, u0 : [50, 50] eine geeignete Anfangs-

funktion und u : [50, 50] [0,[ die gesuchte Lsung.

1

1 Die nicht-eingefrorene Nagumo-Gleichung (1D)

1. Starte COMSOL Multiphysics.

1.1 Modell-Navigator

2. ffne im Modell-Navigator die Registerkarte Neu.

3. Whle dort die Raumdimension 1D.

4. Whle im linken Bildteil Anwendungsmodi COMSOL Multiphysics Partielle Diffe-rentialgleichungen (PDGL) PDGL, Allgemeine Form Zeitabhngige Analyse.

5. Whle unter Abhngige Variablen u, unter Name des Anwendungsmodus PDE und unterElement Lagrange-Linear.

6. Besttige abschlieend auf OK.

1.2 Konstanten festlegen

7. Whle in der Menleiste Optionen Konstanten...

8. Gebe in den Spalten Name, Ausdruck, Wert und Beschreibung die in der folgenden Abbildungenthaltenen Variablen und Werte ein. Die Spalte Wert wird dabei automatisch aus der SpalteAusdruck erzeugt.

Name Ausdruck Wert Beschreibung

a 0.25 0.25 Parameter, 0

1 Die nicht-eingefrorene Nagumo-Gleichung (1D)

16. Whle auf der linken Seite die Registerkarte Gebiete, markiere in der Gebietauswahl das Gebiet1 durch einen Mausklick mit der linken Taste und mache einen Haken bei der Check-Box Aktiv indiesem Bereich.

17. In der Registerkarte Koeffizienten geben wir die folgenden PDGL-Koeffizienten ein

Koeffizient Wert/Ausdruck Beschreibung

-ux FlussvektorF u*(1-u)*(u-a) Quelltermea 0 Massenkoeffizientda 1 Dmpfung/Massen-Koeffizient

18. In der Registerkarte Anfangsbedingung geben wir die folgenden Anfangswerte ein

u(t0) 1/100*(x+50)ut(t0) 0

19. In der Registerkarte Element geben wir in der Element-Auswahl unter Vordefinierte Ele-mente Lagrange-Linear ein (insofern wir dies zu Beginn vergessen haben). Die Eigenschaften sha-pe (Formfunktion), gporder (Integrationsordnung), cporder (Nebenbedingungs-Ordnung) undbnd.gporder (Integrationsgrad fr ultraschwachen Term) bleiben dabei unverndert.

20. In der Registerkarte Schwache Form lassen wir den Inhalt zu Beitrag schwachen Form unver-ndert.

21. Besttige die Eingabe abschlieend mit OK.

Randbedingungen

22. Whle in der Menleiste Physik Randbedingungen (F7)

23. Whle auf der linken Seite die Registerkarte Rnder, markiere in der Randauswahl beide Rnder:Dazu markiere durch einen Mausklick mit der linken Taste den Rand 1, halte die Shift-Taste gedrcktund markiere erneut durch einen Mausklick mit der linken Taste den 2. Rand.

24. In der Registerkarte Koeffizienten geben wir unter Randbedingungen Neumann-Randbedingungein. Fr den Koeffizienten G whlen wir den Wert/Ausdruck 0.

25. Die Registerkarte Schwache Form bleibt unverndert.

26. Besttige die Eingabe abschlieend mit OK.

1.5 Netz-Generierung

27. Markiere im Modell-Verzeichnisbaum den Modell-Namen Geom1 durch einen Mausklick mitder linken Taste.

28. Whle in der Menleiste Netz Freie Netzgenerierungsparameter... (F9)

29. In der Registerkarte Global stellen wir Maximale Elementgre auf 0.1. Der Maximale Skalie-rungsfaktor fr Elementgre wird auf 1 und die Elementvergrerungsrate auf 1.3 gesetzt.

30. Die Registerkarten Gebiet und Rand bleiben unverndert.

31. Klicke nun auf den Button Netz generieren und warte bis das Netz erzeugt wurde.

32. Besttige die Eingabe abschlieend mit OK.

1.6 Modell lsen

1.6.1 Lserparameter

33. Whle in der Menleiste Lsen Lserparameter... (F11)

34. Whle auf der linken Seite unter Lser Zeitabhngig und mache darunter einen Haken bei derCheck-Box Plotten whrend des Lsens.

3

1 Die nicht-eingefrorene Nagumo-Gleichung (1D)

35. In der Registerkarte Allgemein, die sich auf der rechten Seite befindet, geben wir unter Zeitschritt-steuerung die Zeiten 0:0.01:30, die Relative Genauigkeit 0.01 und die Absolute Genauigkeit0.0010 ein. Die Check-Box Komplexe Zahlen erlauben bleibt deaktiviert.

36. Die Registerkarten Zeitschrittsteuerung und Erweitert bleiben in diesem Beispiel unverndert.

37. Besttige die Eingabe mit OK.

1.6.2 Lsermanager

38. Whle in der Menleiste Lsen Lsermanager...

39. In der Registerkarte Anfangswert whlen wir unter Anfangswert Anfangsbedingung und unterWerte fr nicht zu lsende Variablen und Linearisierungspunkt Einstellung aus der An-fangsbedingung verwenden.

40. Besttige die Eingabe abschlieend mit OK.

1.6.3 Modell lsen

41. Whle in der Menleiste Lsen Modell lsen....

42. Beachte: Falls nachtrglich nderungen am Modell (z.B. nderung der Konstanten, Randbedingun-gen, Gebietbedingungen, u.s.w.) vorgenommen werden, muss das Modell zunchst wieder aktualisiertwerden. Dies erfolgt ber die Menleiste: Lsen Modell aktualisieren. Das Modell muss nunerneut gelst werden.

1.7 Postprozessing

1.7.1 Abbildung der Lsung

43. Whle in der Menleiste Postprozessing Plotparameter... (F12)

44. In der Registerkarte Allgemein werden Ploteigenschaften (Linie (Lsungsgraph durch eine Liniedarstellen), Min/Max-Markierung (Markierung von Min/Max des Lsungsgraphen) und Geo-metriebegrenzungen (Hervorhebung der Randpunkte)) angezeigt. Beachte: nderungen in denPlotparametern knnen durch einen Klick mit der linken Maustaste auf Zuweisen direkt (d.h. un-mittelbar) angezeigt werden. Wir aktivieren die Check-Boxen Linie und Geometriebegrenzun-gen. Unter zu verwendende Lsung auf der rechten Seite knnen wir fr Lsung zur Zeit einenbeliebigen Zeitpunkt auswhlen, zu dem wir den Lsungsgraphen betrachten wollen. Wir nehmenhier keine weiteren nderungen vor. Unter Darstellung innerhalb knnen wir auswhlen, ob derLsungsgraph in Comsol (Achsen) oder in einem eigenen Abbildungsfenster (Neue Abbildung)angezeigt werden soll. Wir whlen an dieser Stelle Neue Abbildung. Unter Titel... knnen wir derAbbildung zudem eine individuelle Abbildungsberschrift verleihen. Unter Gltten... knnen wirden Lsungsgraphen zudem glatter machen, insofern dieser unerwartete Kanten in seinem Verlaufaufweist. Wir nehmen auch hier keine weiteren nderungen vor.

45. In der Registerkarte Linie knnen wir unter Liniendarstellung auswhlen, ob der Funktionsgraphangezeigt werden soll. Unter Hhendaten lsst sich bei Vordefinierte Gren einstellen, ob dieLsung selbst oder eine bestimmte zeitliche (bzw. rumliche) Ableitung geplottet werden soll. ImBereich Linienfarbe lsst sich unter Einheitliche Farbe die Farbe des Funktionsgraphen festlegen.An dieser Stelle ndern wir die Farbe des Lsungsgraphen von schwarz auf blau, indem wir auf denButton Farbe... klicken und uns einen entsprechenden Farbton aussuchen.

46. In der Registerkarte Min/Max knnen wir unter Min/Max-Markierung einstellen, ob das Mi-nimum und Maximum der Lsung angezeigt werden sollen. Wir nehmen hier keine nderungenvor.

47. Besttige die Eingabe abschlieend mit OK.

48. Im daraufhin erscheinenden Abbildungsfenster whle zum Speichern der Abbildung das Disketten-symbol in der linken oberen Ecke.

4

1 Die nicht-eingefrorene Nagumo-Gleichung (1D)

49. Im Fenster Bild exportieren knnen die gewnschten Eigenschaften zum Speichern des Bildesvorgenommen werden.

50. Klicke nun auf Exportieren..., whle das gewnschte Verzeichnis sowie den gewnschten Datein-amen und besttige abschlieend auf Export. Das Ergebnis ist in der folgenden Abbildung enthal-ten.

51. Beachte: Htte man in der Registerkarte Allgemein unter Darstellung innerhalb den WertAchsen eingegeben, so muss die Abbildung auf eine andere Art uns Weise gespeichert werden: Dazuwhle in der Menleiste Datei Exportieren Abbildung... und fahre mit den Schritten 78und 79 fort.

1.7.2 Animation des zeitlichen Lsungsverlaufs

52. Whle in der Menleiste Postprozessing Plotparameter... (F12)

53. In der Registerkarte Animation erstellen wir nun ohne groen Aufwand ein kleines Movie. Beachte:Diese lassen sich anschlieend auch problemlos (z.B. mit dem movie15-Paket) in mit Latex erstelltenBeamer-Vortrgen einbinden. Unter Animationseinstellungen whlen wir fr den Ausgabe-TypFilm, Dateityp AVI, Breite (in Pixeln) 640, Hhe (in Pixeln) 480 und unter Bilder pro Se-kunde 20. Unter Erweitert whlen wir als Filmqualitt Beste und besttigen auf OK. Beachte:Falls man die Dateigre des Movies klein halten und zudem einen Qualittsverlust in Kauf nehmenmchte, so sollte man unter Erweitert die Check-Box Animation erstellen aus JPEG-Bildernaktivieren (und eventuell die Filmqualitt reduzieren). Da wir das Beste herausholen wollen, be-lassen wir diese deaktiviert.

54. Besttige die Eingabe abschlieend auf Animation starten und warte bis die Animation beendetist.

55. Im Fenster COMSOL-Animation, in dem die Animation durchgelaufen ist, lsst sich die Anima-tion durch einen Klick mit der linken Maustaste auf das Diskettensymbol speichern.

56. Im daraufhin erscheinenden Fenster Animation speichern als whlen wir das gewnschte Ver-zeichnis sowie den gewnschten Dateinamen und besttigen die Eingaben die einen Klick mit derlinken Maustaste auf den Button Speichern. Wir erhalten das folgende Resultat (Klicke mit derlinken Maustaste auf das folgende Bild).

5

1 Die nicht-eingefrorene Nagumo-Gleichung (1D)

57. Beachte: Um uns in diesem und in den folgenden Abschnitten etwas Zeit zu sparen, haben wir unterLserparameter in der Zeitschrittsteuerung die Zeiten von 0:0.01:30 auf 0:0.1:30 abgendertund das Modell anschlieend erneut gelst. Insofern man die Fehlermeldung Fehler 6092 Speicher-platz berschritten bei der xmesh-Durchfhrung erhlt, sollte man beispielsweise (wie im vorherigenSatz vermerkt) die Zeitschrittweite und somit die Menge an gespeicherten Daten reduzieren.

1.7.3 Modellbereiche darstellen

Zeitliche Vernderung der Wellenfront

58. Whle in der Menleiste Postprozessing Modellbereiche darstellen...

59. In der Registerkarte Allgemein whle unter Plottyp die Linien/Extrusions-Darstellung aus.

60. In der Registerkarte Linie/Extrusion setze Plottyp auf Extrusionsdarstellung. Weiter setze Vor-definierte Gren auf u.

61. Besttige die Eingabe abschlieend mit OK. Das Ergebnis ist in den folgenden zwei Abbildungenthalten.

6

Nonfrozen-Nagumo-Equation-1D-01.aviMedia File (video/avi)

1 Die nicht-eingefrorene Nagumo-Gleichung (1D)

1.8 Ergnzende Einstellungen

62. Whle in der Menleiste Datei Modelleigenschaften....

63. In der Registerkarte Allgemein whlen wir fr Modellname Nonfrozen Nagumo-Equation (1D),fr Autor Vorname Nachname, fr Firma Universitt Bielefeld, fr Abteilung Fakultt fr Ma-thematik, das Feld Referenz lassen wir leer und unter URL geben wir http://www.math.uni-bielefeld.de.

7

1 Die nicht-eingefrorene Nagumo-Gleichung (1D)

64. In der Registerkarte Beschreibung vermerken wir in der Modellbeschreibung Dieses Modellberechnet die Lsung der nicht-eingefrorenen eindimensionalen Nagumo-Gleichung fr das Intervallvon -50 bis 50.

65. In der Registerkarte Bemerkung fhren wir unter Ergebnisbeschreibung auf Die Lsung ist eine(mit der Zeit) wandernde Wellenfront. Die Geschwindigkeit der Front ist bei der nicht-eingefrorenenGleichung unbekannt. Da die Geschwindigkeit ungleich 0 ist, luft die Wellenfront generell aus demBerechnungsbereich hinaus.

66. Die Registerkarte Info wird von COMSOL automatisch ausgefllt.

67. Besttige die Eingabe mit OK.

68. Zuletzt werden wir das Modell speichern. Dazu whle in der Menleiste Datei Speichern un-ter....

69. Wechsle im Fenster Speichern unter in das gewnschte Verzeichnis und gebe den DateinamenNonfrozen-Nagumo-1D.mph ein.

70. Besttige die Eingabe mit Speichern und beende COMSOL Multiphysics.

8

Einleitung (Mathematischer Hintergrund)1 Die nicht-eingefrorene Nagumo-Gleichung (1D)1.1 Modell-Navigator1.2 Konstanten festlegen1.3 Grundgebiet erzeugen1.4 Differentialgleichung definieren1.4.1 PDGL, Allgemeine Form (PDE)

1.5 Netz-Generierung1.6 Modell lsen1.6.1 Lserparameter1.6.2 Lsermanager1.6.3 Modell lsen

1.7 Postprozessing1.7.1 Abbildung der Lsung1.7.2 Animation des zeitlichen Lsungsverlaufs1.7.3 Modellbereiche darstellen

1.8 Ergnzende Einstellungen