Die vorliegende Arbeit befaßt sich mit schwachbeladenen polydispersen turbulenten Mehrphasenströmungen. Das Verfahren, das zur Beschreibung der Zweiphasenströmung verwendet wird, basiert auf der Euler-Lagrange Methode. Die Turbulenz in der kontinuierlichen Phase wird mittels Grobstruktursimulation (Large Eddy Simulation, LES) modelliert. Die meisten Untersuchungen mittels Direkten Numerischen Simulation (DNS) neueren Datums, die einen Einblick in die Wechselwirkungsmechanismen zwischen Partikeln und umgebendem Fluid gewähren, beschreiben einfache Strömungen, sehr oft mit homogener Turbulenz und sehr niedrigen Reynolds Zahlen. Im Gegensatz dazu eröffnet LES die Möglichkeit diese Studien auf Strömungen mit wesentlich höheren Reynolds Zahlen und komplexeren Geometrien durchzuführen. Des Weiteren ist das Verfahren weniger restriktiv gegenüber Partikeln, die größer als die Kolmogorov Längenskalen sind. Die wesentliche Zielsetzung dieser Arbeit ist es, den Stand des Wissens bezüglich Zweiphasen-LES einer kritischen Prüfung zu unterziehen und hierbei Schwachpunkte zu identifizieren. Es werden Vorschläge für die weitere Modellentwicklung gemacht, um die Eignung der Methode zur Vorhersage von mehrphasigen Strömungsphänomenen zu steigern. Die Zuverlässigkeit der Übertragung von Erkenntnissen von einer Konfiguration auf andere wird abgeschätzt. Die Schlussfolgerungen sind auf einer systematischen Parameterstudie der relevanten Strömungskenngrößen, wie der Reynolds Zahl oder der Stokes Zahl, begründet, und sichern somit einen großen Anwendungsbereich. Hierzu werden einige partikelbeladene Strömungskonfigurationen untersucht. Diese sind: zwei ebene Kanalströmungen, ein Freistrahl und ein Spray mit Verdampfung bei einer vergleichsweise niedrigen Temperatur. Besondere Beachtung findet die Vorhersage der bevorzugten Akkumulation der Partikel in einer turbulenten Umgebung. Diese ist von besonderer Wichtigkeit bei Simulationen von Mischungs- und Verbrennungsvorgängen in turbulenten Strömungen. Mittels einer systematischen Parameterstudie der Partikel Stokes Zahl (hier definiert als das Verhältnis zwischen einer charakteristischen Partikelzeitskala und der Kolmogorov Zeitskala) und der Reynolds Zahl der Strömung, wird aufgezeigt, dass die maximale bevorzugte Akkumulation immer bei einer konstanten Stokes Zahl auftritt und dass diese Stokes Zahl nicht von der Reynolds abhängt. Die Intensität der Akkumulation wiederum ist eine Funktion der Reynolds Zahl. Der Einfluss eines stochastischen Kleinskalen-Dispersionsmodells auf die Partikelakkumulation ist bewertet und als unwesentlich für Partikel eingestuft, deren charakteristische Zeitskalen in derselben Größenordnung wie die Kolmogorov Zeitskala liegen. Eine zusätzliche Schwierigkeit für die Zweiphasen-LES ergibt sich aus der Tatsache, dass die Standard Euler-Lagrange Methode unphysikalisch hohe Partikelanzahldichten in den Bereichen bevorzugter Konzentration zuläßt, da hier die Partikel-Partikel Kollision nicht berücksichtigt wird. Die Verwendung eines Modells das diese beschreibt, scheint unumgänglich. Der Einfluss der dispersen Phase auf die Fluidturbulenz wird in einer geschlossenen und in einer freien Scherströmung betrachtet. Die ungenaue Modellierung der Partikel-Wand Interaktion in wandgeschlossenen Strömungen kann zu unterschiedlichen Lösungen für das Partikelgeschwindigkeitsfeld führen, und daher auch zu unterschiedlichen Intensitäten der Turbulenzmodifikation. Es sollte berücksichtigt werden, dass diese unterschiedlichen Lösungen keineswegs unphysikalisch sein müssen, wie auch experimentelle Untersuchungen zeigen. Im Gegensatz dazu wird die Turbulenzmodifikation in der freien Scherströmung sehr genau vorhergesagt. Darüber hinaus wird gezeigt, dass diese Vorhersage auch für unterschiedliche Massenbeladungen der dispersen Phase zuverlässig ist. Im Rahmen dieser Arbeit wird ein effizientes nicht iteratives Fractional Step Verfahren für die Druckgeschwindigkeitskopplung des umgebenden Fluids verwendet. Die Kopplung dieser Methode mit der Zwei-Wege gekoppelten Euler-Lagrange Methode ist erfolgreich für Konfigurationen isothermer Strömungen und solche mit kleiner Dichteänderung validiert, wobei sowohl experimentelle als auch DNS Daten als Referenz verwenden werden.