Barnes-Hut-Algorithmus

Der Barnes-Hut-Algorithmus ist ein 1986 von Josh Barnes und Piet Hut veröffentlichtes Näherungsverfahren, das eine effektive Berechnung der Kräfte in einem N-Körperproblem ermöglicht. Im Gegensatz zur direkten Aufsummierung der Kräfte, deren Rechenaufwand mit $ O(N^{2}) $ ansteigt, reduziert sich der Aufwand beim Barnes-Hut-Algorithmus auf $ O(N\log N) $.

Motivation

Die Simulation von N-Körperproblemen ist eine Standardaufgabe der Astronomie. Bei einem solchen Problem bewegt sich eine Anzahl (N) von Körpern unter dem Einfluss einer Kraft, die wiederum von den Positionen aller andern Körper abhängt. Als Beispiel sei hier die Bewegung von Sternen im Gravitationsfeld einer Galaxie verwiesen. Die direkte Integration ist mit zunehmender Anzahl an Körpern sehr aufwendig, da die auf einen Körper wirkende Kraft die Summe aller Kräfte ist die von den anderen Körpern auf ihn wirken. Berechnungen auf Basis der direkten Kräftesummation werden daher mit steigender Körperanzahl (N>10000) schnell uneffektiv, denn die Gesamtzahl der Kraftberechnungen ist:

$ {\frac {1}{2}}N(N-1)\approx O(N^{2}) $

Diesem Trend lässt sich zwar durch Verwendung hochparallelisierter Computerhardware (GPU) entgegenwirken, allerdings verschiebt sich das Problem damit lediglich zu größeren Teilchenanzahlen. Für effektive Simulationen mit mehreren Millionen oder gar Milliarden Partikeln sind Algorithmen notwendig, deren Berechnungsaufwand nicht quadratisch mit der Teilchenzahl ansteigt. Einer dieser Algorithmen ist der von Barnes und Hut.

Partikelverteilung die zwei benachbarten Galaxien nachempfunden ist.

Vollständiger Barnes-Hut-Baum. Das Gebiet wird solange rekursiv unterteilt, bis sich jedes Teilchen in einer eigenen Zelle befindet.

Die Zellen des Barnes-Hut Baumes, die zur Berechnung der Kraft auf ein Teilchen im Koordinatenursprung herangezogen werden.

Die Berechnungen für Wechselwirkende Galaxien ist ein typischer Anwendungsfall für Barnes-Hut-Simulationen.

Funktionsweise

Der Barnes-Hut-Algorithmus verringert die Anzahl zu berechnenden Kräfte durch geeignetes Zusammenfassen von Teilchengruppen zu Pseudoteilchen. Die Grundidee dabei ist, dass die von einer Partikelgruppe auf ein Einzelpartikel ausgeübte Kraft in sehr guter Näherung durch die Wirkung einer einzelnen Masse im Massenschwerpunkt der Partikelgruppe approximiert werden kann. So kann man beispielsweise die Kraft, die von der Andromeda-Galaxie auf die Sonne ausgeübt wird, sehr gut durch eine Punktmasse nähern, die sich im Zentrum der Andromeda-Galaxie befindet und die deren Masse hat. Diese Näherung ist allerdings nur zulässig, wenn der Abstand $ r $ der Gruppe vom Einzelteilchen groß und der Gruppendurchmesser $ d $ im Verhältnis zum Abstand klein ist. Das Verhältnis von Gruppendurchmesser zu Gruppenentfernung wird als Multipol-Akzeptanz Kriterium (engl.: MAC; Multipole-Acceptance-Criterion) bezeichnet:

$ \theta ={\frac {d}{r}} $

Je kleiner $ \theta $, umso besser die Approximation. Überschreitet $ \theta $ einen bestimmten Schwellwert, so sollte die Näherung nicht mehr angewendet werden, um größere Fehler zu vermeiden. Der Algorithmus setzt dieses Prinzip durch rekursives Unterteilen des Simulationsbereiches in Quadranten (2D) bzw. Oktanten (3D) um. Die Partikel werden in den Knoten der so entstandenen Baumstruktur gespeichert. Ist die Entfernung eines Knotens von einem Einzelteilchen groß genug, dann erfolgt die Kraftberechnung nicht mehr direkt zwischen den Teilchen, sondern zwischen dem Einzelteilchen und dem Massenschwerpunkt des Knotens.

Die von einer Partikelverteilung auf ein Einzelpartikel ausgeübte Kraft kann durch die von einer Punktmasse ausgehenden Kraft angenähert werden, wenn die Entfernung der Partikelgruppe zum Einzelteilchen groß im Verhältnis zur Gruppengröße ist.

Die gravitative Wirkung von Sternhaufen und Stern B auf Stern A kann infolge der großen Entfernung als Punktmasse approximiert werden. Doch auch innerhalb des Sternhaufens kann die Gravitationskraft des Sternhaufens auf Stern B durch eine Punktmasse angenähert werden, da Stern B weit genug vom Sternhaufen entfernt ist.

Siehe auch

• Dreikörperproblem

Weblinks

Deutsch

• Detaillierte Beschreibung der Implementierung des Barnes-Hut-Algorithmus

Englisch

• J. Barnes und P. Hut: A hierarchical O(N log N) force-calculation algorithm in Nature, 324(4), Dezember 1986
• Fast Hierarchical Methods for the N-body Problem, Part 1, CS267, Lecture 24, 11. April 1996
• Video einer Partikelsimulation mit 8,7 Millionen Partikeln unter Verwendung des Barnes-Hut-Algorithmus gerechnet auf mehreren GPUs (youtube.com)
• PEPC – The Pretty Efficient Parallel Coulomb Solver, eine parallele open-source Implementation des Barnes-Hut Algorithmus mit austauschbarem Wechselwirkungskern für vielseitige Anwendungen