Numerische Methoden stellen einen wesentlichen Bestandteil des  Auslegungsprozesses moderner Turbomaschinen dar und tragen zu einer erheblichen Steigerung der Zeit- und Ressourceneffizienz bei. Eine Vielzahl der hierzu verwendeten Verfahren beruht auf der Annahme eines kalorisch perfekten Gases, welche im Kontext konventioneller Arbeitsmedien wie beispielsweise Luft eine geeignete Näherung der realen Fluideigenschaften ermöglicht. Im Zuge der Erschließung alternativer Energiequellen erlangen zunehmend solche Prozesse eine besondere Relevanz, für welche eine auf den Anwendungsfall abgestimmte Wahl des Arbeitsmediums möglich ist. Die in der Folge aufgrund ihrer besonderen thermophysikalischen Eigenschaften eingesetzten Fluide lassen sich nur unzureichend durch das Modell des kalorisch perfekten Gases beschreiben, was ihre Klassifikation als nicht-ideal bedingt. Da eine Berücksichtigung dieser Nicht-Idealität sowie eventueller Phasenwechselvorgänge anhand etablierter numerischer Methoden häufig nicht möglich ist, bedarf die Auslegung der zugehörigen Turbomaschinen zunächst einer Entwicklung geeigneter Verfahren. Der im Rahmen dieser Arbeit vorgestellte Ansatz zur Berechnung numerischer Flussterme leistet hierzu einen Beitrag und kann sowohl auf ein- als auch zweiphasige Strömungen nicht-idealer Fluide angewendet werden. Er zeichnet sich, neben einer hohen numerischen Robustheit gegenüber Diskontinuitäten, insbesondere durch seine Anwendbarkeit auf beliebige Zustandsgleichungen aus. Durch eine Vielzahl an Modellierungsarten ist eine Beschreibung der im Zuge eines Phasenwechsels entstehenden dispersen Phase mit unterschiedlichem Detailgrad möglich. So kann die Tropfengrößenverteilung entweder als mono- oder polydispers modelliert werden, wobei sie im letztgenannten Fall anhand ihrer statistischen Momente beschrieben wird. Auch ermöglicht das entwickelte Verfahren eine Berücksichtigung von Geschwindigkeitsdifferenzen zwischen den Phasen, wodurch die Basis für eine detaillierte Untersuchung der Bewegung der dispersen Phase im Strömungsfeld geschaffen wird. Anhand einer Tabellierung der thermophysikalischen Größen kann zudem eine im Kontext der Auslegung relevante Steigerung der Rechenzeiteffizienz erreicht werden. Die hiermit verbundene Einbuße in Bezug auf die Beschreibungsgenauigkeit wird in Folge der Verwendung eines Taylorreihen-Ansatzes zur Interpolation minimiert. Durch eine Verifizierung und Validierung anhand einer repräsentativen Auswahl an Testfällen kann die Anwendbarkeit des entwickelten Verfahrens auf ein- und zweiphasige Strömungen kompressibler nicht-idealer Fluide unterschiedlicher molekularer Komplexität nachgewiesen werden. Dabei wird insbesondere der Phasenwechsel in Form der homogenen Nicht-Gleichgewichtskondensation in hoher Übereinstimmung mit den Ergebnissen experimenteller Untersuchungen beschrieben.