Die frueher entwickelte dreidimensionale einphasige Version des Rechenprogrammes BACCHUS zur Beschreibung der Thermohydraulik eines Kuehlmittels in hexagonaler Buendelgeometrie lieferte die Grundlage fuer die Entwicklung der zweiphasigen Version, die in diesem Bericht dokumentiert wird. Das Zwei-Phasen-Schlupfmodell mit dem Spezialfall des homogenen Gleichgewichtsmodells wird detailliert beschrieben. Im Vergleich zu den ueblichen Modellen, die in der Literatur beschrieben werden, enthaelt es mehrere Verbesserungen sowohl in der physikalischen Modellierung als auch in der numerischen Behandlung der Grundgleichungen. Das Modell der getrennten Phasen fuer hoehere Dampfgehalte wird in allen analytischen Einzelheiten beschrieben, die zum vollen Verstaendnis seiner Implementierung im Rechenprogramm noetig sind. Probleme im Zusammenhang mit dem Uebergang zwischen den beiden Modellen in einer integrierten Programmversion werden anschliessend diskutiert. Das Rechenprogramm enthaelt eine zusaetzliche Option zur Modellierung von aktiven und passiven, durchlaessigen und undurchlaessigen Blockaden. Diese Option wird zusaetzlich dokumentiert. Neue numerische Methoden zur Loesung des algebraischen Gleichungssystems, das durch die Linearisierung der Grundgleichungen entsteht, haben die Methoden, die frueher fuer die einphasige Version des Rechenprogramms benutzt wurde, vollstaendig ersetzt und werden in diesem Bericht detailliert beschrieben. Anschliessend ist ein Abschnitt einem Ueberblick ueber die Verifikation des Programms gewidmet, die sich ueber mehrere Jahre erstreckt und von stationaeren einphasigen unbeheizten Buendelexperimenten bis zu schnellen transienten zweiphasigen Experimenten in elektrisch beheizten 37-Stabbuendeln geht. (orig.)