Strenge zukünftige Limitierungen des Schadstoffausstoßes werden große Probleme bei dem Betrieb von Gasturbinen- oder Flugtriebwerksbrennkammern bereiten. Um den steigenden Anforderungen an diese Systeme Rechnung zu tragen ist es notwendig, den involvierten turbulenten Verbrennungsprozess genauer zu untersuchen. Dafür eignen sich detaillierte numerische Untersuchungen. Obwohl die direkte numerische Simulation (DNS) turbulente Strömungen ohne Modellierung abbilden kann, ist sie dennoch aufgrund des hohen Bedarfs an Speicherplatz und Rechenzeit beschränkt auf einfache Konfigurationen. Weniger kritisch ist dieser Punkt in Bezug auf die Grobstruktursimulation (LES, Large Eddy Simulation), welche sich in jüngster Zeit als ein vielversprechendes Werkzeug erwiesen hat, um turbulente Verbrennungsprozesse abzubilden. Im Vergleich zu statistischen Ansätzen (RANS, Reynolds Averaged Navier Stokes) bietet die LES den entscheidenden Vorteil, zeitaufgelöste Informationen zu liefern und trägt damit der inhärent instationären Natur der Strömung Rechnung. Daher ist diese Methode Gegenstand dieser Arbeit, welche aus den folgenen zwei Teilen besteht. Zunächst werden verschiedene Ansätze zur Modellierung der Feinstruktur (SGS, Subgrid-Scale) untersucht. Das Ziel ist dabei, ein besseres Verständnis für die Sensitivitäten der Feinstrukturmodellierung zu erlangen und ein geeignetes, möglichst einfaches Modell für die nachfolgenden Untersuchungen auszuwählen. Hierfür wird der Effekt von verschiedenen Formulierungen des Modells von Smagorinsky mit jeweils verschiedenen Filterweiten untersucht sowie weiterhin die Robustheit der dynamischen Prozedur, das Verhalten verschiedener dynamischer Modelle und kombinierter Modelle und schließlich die Performanz eines Eingleichungsmodells. Die bisher genannten Modelle sind dennoch nicht geeignet, die anisotrope Natur der Strömung wiederzugeben. Daher wird in dieser Arbeit ein anisotropes Eingleichungsmodell basierend auf der kinetischen Energie der Feinstruktur (SGS) entwickelt. Zugrunde liegt eine anisotrope Formulierung der Wirbelviskosität (Eddy Viscosity). Als generischer Testfall für die vorgestellten Modelle dient die Simulation einer turbulenten Kanalströmung. Schließlich werden einige der Modelle noch für die Simulation der Strömung in einem rotierenden Kanal benutzt. Es kann gezeigt werden, dass das entwickelte Modell gute Übereinstimmung mit experimentellen Daten erzielt. Weiterhin wird noch ein dynamischer Lagrange-Ansatz verwendet. Der zweite Teil beschäftigt sich mit der Grobstruktursimulation turbulenter, vorgemischter Verbrennung. Im speziellen werden das Flame Surface Density Modell sowie der ATF-Ansatz (Artificially Thickened Flame) behandelt. Zu diesem Zweck wird ein einfacher Ansatz zur Dichtekopplung vorgeschlagen, um den inkompressiblen Löser zu erweitern. Dazu wird die Divergenz des Geschwindigkeitsfeldes berechnet und die Dichte wird damit korrigiert. Die Divergenz wird auch für die konvektive Ausstromrandbedingung verwendet. Weiterhin werden eine Behandlung der Flamme-Wand-Interaktion sowie verschiedene Modelle zur Flammfrontfaltung vorgestellt. Als Testfall für die Simulationen dient eine durch einen Staukörper stabilisierte Flamme. Der Vergleich von experimentellen Messdaten und numerischen Ergebnissen zeigt, das die implementierte numerische Methode eine vielversprechende Übereinstimmung liefert.