Berechnung kompressibler Freistrahlen mit polynärer Diffusion und Auftrieb