Periodische Randbedingung

aus Wikipedia, der freien Enzyklopädie
(Weitergeleitet von Periodische Randbedingungen)
Wechseln zu: Navigation, Suche

Periodische Randbedingungen werden in analytischen oder numerischen Modellrechnungen gewählt, um eine gesonderte Behandlung von Rändern zu vermeiden oder das Gebiet, über den sich die Rechnung erstreckt, verkleinern zu können. Periodische Randbedingungen in d Dimensionen können als Kompaktifizierung des Raumes zu einem Torus in d+1 Dimensionen aufgefasst werden.[1]

Anwendungsgebiete sind in der Festkörperphysik kristalliner Materialien, Molekulardynamik, Monte-Carlo-Simulationen sowie bei Simulationen auf Gittern wie Gittereichtheorien.

Periodische Randbedingungen in zwei Dimensionen. Der Rand der Simulationsbox ist schwarz dargestellt.

Bei kontinuierlichen Teilchensimulationen mit periodischen Randbedingungen treten Teilchen an einem Rand der Simulationsbox aus und treten am gegenüberliegenden Rand wieder ein.

Mathematische Definition[Bearbeiten]

Eine partielle Differentialgleichung ist eine Gleichung, in der Ableitungen einer unbekannten Funktion

u:\Omega\to\R

für eine offene Teilmenge \Omega\subset\R^d, d\ge 2 vorkommen. Man spricht von periodischen Randbedingungen, wenn \Omega von der Form

\Omega=(a_1,b_1)\times\ldots\times (a_d,b_d)\subset\R^d

ist und man verlangt, dass es eine stetige Fortsetzung u:\overline{\Omega}\to\R mit

u(x_1,\ldots,x_{i-1},a_i,x_{i+1},\ldots,x_d)=u(x_1,\ldots,x_{i-1},b_i,x_{i+1},\ldots,x_d)\ \forall (x_1,\ldots,x_d)\in\Omega, i=1,\ldots,d

geben soll.[2]

\Omega muss nicht quaderförmig sein, aber eine vollständige periodische Abdeckung des Raumes erlauben[3] (siehe Raumfüllung und Parkettierung).

Kompaktifizierung[Bearbeiten]

Eine rechteckige Fläche wird zu einem Torus zusammengelegt; ihre Ränder verschwinden.
Die für die 3D-Ansicht nötige Verzerrung tritt in der Anwendung nicht auf.

Periodische Randbedingungen in allen d Dimensionen entsprechen einem flachen Torus in d+1 Dimensionen, d. h. nicht gekrümmter Oberfläche. Die Oberfläche ist deshalb nicht als Teilmenge des d+1-dimensionalen Raumes beschreibbar, sondern als kartesisches Produkt von d Kreisen.

Erhaltungsgrößen[Bearbeiten]

Periodische Randbedingungen erlauben Energie- und Impulserhaltung, verletzen aber die Drehimpulserhaltung. Formal: Drehimpulserhaltung ist Folge der hier nicht gegebenen Invarianz der physikalischen Gesetze gegenüber Drehungen des Bezugssystems. Anschaulich: Ein lokal gestarteter Wirbel wächst durch Impulsdiffusion bis auf die Längenskala der Simulationsbox an und wird dann systematisch vernichtet.

Beispiele[Bearbeiten]

Periodische Lösung für ein periodisches Problem[Bearbeiten]

Im ersten Beispiel werden periodische Randbedingungen mit verschiedenen nichtperiodischen Randbedingungen kombiniert. Es geht um die Simulation der Umströmung von Turbinenschaufeln in einem Ringspalt. Es gelte die Annahme (oder Näherung), dass die Umströmung aller N Schaufeln gleich sei. Das Problem hat damit eine N-zählige Drehsymmetrie, die ausgenutzt werden soll, indem als Simulationsgebiet ein Sektor gewählt wird, der lediglich eine der Schaufeln enthält. Am Ende soll die für den Sektor gefundene Lösung N-fach vervielfältigt zu dem Ring zusammengesetzt werden. Der Rechenaufwand beschränkt sich damit auf einen Bruchteil des Gesamtgebiets und die Gesamtlösung hat automatisch die erwartete Symmetrie, ist periodisch. Damit die Lösung aber physikalisch sinnvoll sein kann, muss sie an den Nahtstellen zusammenpassen, so glatt wie es von der Lösung auch im Inneren des Simulationsgebiets verlangt wird. Diese Forderung macht die periodischen Randbedingungen aus. Zudem muss in diesem speziellen Beispiel die Lösung an den festen Innen- und Außenflächen entsprechende Randbedingungen erfüllen (z. B. Strömungsgeschwindigkeit = Wandgeschwindigkeit, falls die Grenzschicht im Detail modelliert wird) und stromauf- und -abwärts mit den Lösungen für weitere Turbinenstufen passen (es gilt, eine gemeinsame Lösung zu finden).

Wie man zu Lösungen kommt, die die periodischen Randbedingungen erfüllen, hängt davon ab, wie überhaupt die physikalischen Größen des Problems (hier Druck, Temperatur und Geschwindigkeitskomponenten) repräsentiert werden sollen. Eine Möglichkeit wäre, als Summe von Basisfunktionen, die einzeln jeweils periodisch und glatt sind, also z. B. sinusförmig vom Winkel um die Drehachse abhängen, siehe Fourier-Reihe. Ein Rand tritt dann nicht auf. Eine andere, hier geeignetere Klasse von Verfahren nutzt zahlreiche Stützstellen (Gitterpunkte), siehe z. B. Finite-Volumen-Verfahren. Dabei werden Gitterpunkte an einer Seite des Simulationsgebiets mit solchen auf der anderen Seite als benachbart definiert, um der Randbedingung Genüge zu tun.

Ignorieren von Oberflächenwellen[Bearbeiten]

Hauptartikel: Born-von-Kármán-Modell

In einem quaderförmigen Einkristall sind Phononen Anregungen mechanischer stehender Wellen. Ignoriert man, dass die Atome in Randlage abweichende Nachbarschaftsverhältnisse haben, vereinfacht sich die Wellengleichung auf die für das unendlich ausgedehnte Gitter. Die Lösungen sind Überlagerungen ebener Wellen der Art

u(x,y,t)=u_0 e^{\mathrm i (k_x x + k_y y - \omega t)}\,.

Darin sind die Komponenten von u die der Auslenkungen aller Atome der kristallographischen Elementarzelle aus ihren Ruhelagen, x und y Ortskoordinaten (der Einfachheit halber in nur zwei Dimensionen), t die Zeit und \omega die Frequenz, die abhängig ist vom Wellenvektor k und von der im Amplitudenfaktor u_0 enthaltenen Polarisation.

Der einzige berücksichtigte Oberflächeneffekt ist die Beschränkung auf diskrete Wellenvektoren:

k_x = \frac{\pi\,n_x}{N_x a_x}
k_y = \frac{\pi\,n_y}{N_y a_y}

mit natürlichen Zahlen n und der Größe N des Kristalls in Einheiten der Größe a der Elementarzelle des Kristalls. Mit diesen k erfüllen die Lösungen die Periodizitätsbedingung

u(x,y,t) = u(x+2 N_x a_x,y,t)=u(x,y+2 N_y a_y, t).

Im Born-von-Kármán-Modell werden periodische Randbedingungen auch Born-von Kármán-Randbedingungen genannt.[4]

Kompression ohne Schraubstock[Bearbeiten]

Das Bild des Schraubstocks steht für große äußere Kräfte. Würde man simulierte Schraubstockbacken auf die Atome einer kleinen Simulationsbox einwirken lassen, wären die Ergebnisse unbrauchbar. Besonders große Kräfte sind nötig, um Materialien bei nicht im Labor realisierbaren Bedingungen, wie jenen im Erdkern, zu untersuchen (z. B. ihren Elastizitätstensor zu bestimmen). Durch die Verwendung periodischer Randbedingungen wird dieses Problem umgangen, indem nicht Kräfte (mechanische Spannungen) vorgegeben und Verformungen beobachtet werden, sondern umgekehrt. Die Simulationsbox enthält dazu eine einzige Elementarzelle des Kristalls und wird willkürlich verformt (Kantenlängen und Winkel). Für jede Geometrie der Elementarzelle wird die Lage der Atome variiert und jeweils die elektronische Energie berechnet (siehe Born-Oppenheimer-Näherung), wobei die Dichtefunktionaltheorie mit periodischen Basisfunktionen zum Einsatz kommt. Aus der Abhängigkeit der elektronischen Energie von der Geometrie der Elementarzelle bei relaxierter Anordnung der Atome ergeben sich die äußeren mechanischen Spannungen als Rechenergebnis.

Anwendungsgebiete mit nichtperiodischer Problemstellung[Bearbeiten]

Der Vorteil der Wandlosigkeit durch periodische Randbedingungen lässt sich auch bei der Simulation von Systemen nutzen, die eigentlich nicht periodisch sind. Die Größe der Simulationsbox wird dann willkürlich, aber nicht beliebig gesetzt: Sie muss größer sein als die Entfernungen, über die Korrelationen auftreten. Das kann ein großes Molekül in Lösung sein, eine Materialverformung oder eine Dichteschwankung nahe eines Phasenübergangs.

Für kurzreichweitige Wechselwirkungen können Cutoff-Radien eingeführt werden, ab denen keine explizite Wechselwirkung zwischen Teilchen mehr berechnet wird, allerdings können weiterhin analytisch erhaltene Zusatzterme für die abgeschnittenen Wechselwirkungen berücksichtigt werden[5]. In dem Fall der Verwendung dieser Cutoff-Radien kann man ein Kriterium für eine nicht zu unterschreitende Größe der Simulationsbox angeben, wobei häufig noch andere Kriterien zu beachten sind, die deutlich größere Boxgröße erzwingen: Der kleinste Durchmesser der Simulationsbox sollte mindestens doppelt so groß sein, wie der größte derartig verwendete Cutoff-Radius, ansonsten sieht ein Teilchen in der Simulationsbox (Central Box) eine Kopie seiner selbst in einer benachbarten Box. Weiterhin ist dies das Kriterium, damit man die Minimum Image Convention für diese Wechselwirkungen verwenden darf. Diese Konvention besagt, dass man nur Wechselwirkungen mit Teilchen in den nächsten benachbarten Boxen der Central Box beachten muss.

Der folgende Absatz ist nicht hinreichend mit Belegen (beispielsweise Einzelnachweisen) ausgestattet. Die fraglichen Angaben werden daher möglicherweise demnächst entfernt. Bitte hilf der Wikipedia, indem du die Angaben recherchierst und gute Belege einfügst. Näheres ist eventuell auf der Diskussionsseite oder in der Versionsgeschichte angegeben. Bitte entferne zuletzt diese Warnmarkierung.
Siehe Diskussionsseite.

Thermodynamischer Limes und Kontinuumslimes[Bearbeiten]

Bei vielen Molekulardynamiksimulationen wird der thermodynamische Limes, d. h. N \to \infty, V \to \infty bei N/V =\text{const}, ermittelt. Dazu werden die Messgrößen für verschiedene Größen der Simulationsbox zu unendlichen Volumen approximiert. In der Regel verschwinden die Effekte der periodischen Randbedingungen in diesem Grenzfall.

Analog wird bei Gittersimulationen der Kontinuumslimes, d. h. der Grenzfall für einen Gitterabstand a\rightarrow0, ermittelt.

Einzelnachweise[Bearbeiten]

  1. Hugo Reinhardt: Quantenmechanik 1: Pfadintegralformulierung und Operatorformalismus. ISBN 348671516X, S. 312 (Google Books).
  2. Jürgen Jost: Partial differential equations (= Graduate Texts in Mathematics. 214) Third edition, Springer, New York 2013, ISBN 978-1-4614-4808-2<--Seite?-->.
  3. Oleg Jardetzky, Michael D. Finucane (Hrsg.): Dynamics, Structure, and Function of Biological Macromolecules. IOS Press, 2001, ISBN 1586030329, S. 3 (Google Books).
  4. Mireille Defranceschi, Claude Le Bris: Mathematical Models and Methods for Ab Initio Quantum Chemistry. (= Lectures in Chemistry. 74), Springer, 2000 ISBN 3540676317, S. 96 (Google Books).
  5.  M. P. Allen, D. J. Tildesley: Computer Simulation of Liquids. Oxford University Press, 1989, ISBN 0198556454, S. 24 (eingeschränkte Vorschau in der Google-Buchsuche).