A hexagonal collision model for the numerical solution of the Boltzmann equation In kompakter Form werden die Hauptergebnisse der diskreten Boltzmann Gleichung basierend auf hexagonalen Elementen vorgestellt. Zwecks Lösung dieses Problems mittels des hexagonalen diskreten Geschwindigkeitsmodells werden im R2̂ automatisch beliebig große Sechseckgitter generiert. Zur Identifikation jeder Sechseckstruktur wird gezeigt dass der Mittelpunkt eines beliebigen regulären Hexagons entweder in die Mitte eines Basishexagons fällt oder ein Gitterknoten ist. Wir beweisen dass bei Zugrundelegung des binären Stoßgesetzes der globale Stoßoperator in einem beschränkten Sechseckgitter in R2̂ nur eine künstliche Invariante besitzt die auch aufgezeigt wird. Wir formulieren ein N Schicht Modell zum Aufstellen von generellen Formeln für alle möglichen regulären Hexagons auf dem Gitter G_N dieser Schicht und beweisen damit ihre Existenz. Dazu bestimmen wir den numerischen Aufwand (flops) zur Auswertung des Boltzmann Stoßoperators im N Schicht Modell. Weiterhin entwickeln wir die kinetische Theorie der diskreten Boltzmann Gleichung für eine hexagonale Diskretisierung in R3̂. Das hexagonale Stoßmodell in R3̂ wird vorgestellt und das lokale Stoßmodell dazu ist ein 12 Geschwindigkeitsmodell entsprechend den 12 Ecken eines kubischen Oktahedron (hexagonaler Kubus bzw. h Kubus). Die Berücksichtigung nur des binären Stoßgesetzes in dem lokalen Stoßmodell führt hier auf drei künstliche Invarianten. Aber bei Einbeziehung des Drei Teilchen Stoßgesetzes wird das Auftreten dieser künstlichen Invarianten vermieden. Wir beweisen dass das 3D hexagonale Modell die grundlegenden Eigenschaften der klassischen kinetischen Theorie erfüllt. Schließlich zeigen wir noch dass dieses 3D Modell der genannten Theorie auch mit dem 2 Teilchen Stoßgesetz ab dem 216 Geschwindigkeitsmodell genügt. Wir präsentieren die Konstruktionen der Gleichgewichtsverteilung für das allgemeine 2D N Schicht Modell und für das 3D Modell wobei die Gleichgewichtsverteilung sich auf die Parameter von Masse Momenten und kinetischen Energie beziehen. Wir geben dazu numerische Ergebnisse für das 2D als auch 3D hexagonale Modell an.