No label defined (Q3838017)

Zur numerischen Lösung von Zweiphasenströmungen in einem heterogenen, porösen Medium benutzen wir den impes-Ansatz. Am mathematischen Modell wird deutlich, daß drei wichtige Fälle bei der numerischen Modellierung zu unterscheiden sind. Es handelt sich um den konvektions- und den diffusionsdominanten sowie den stark nichtlinearen Fall. Bei den Vergleichsrechnungen wurden die Beispiele in diesen drei Klassen eingeteilt, und so konnte überprüft werden, wie sich das entwickelte Verfahren im Vergleich zu anderen Methoden aus der neueren Literatur in typischen Situationen verhält. Im konvektionsdominanten Fall stellt die vorgestellte Methode ein effektives Verfahren dar. Insbesondere werden Unstetigkeiten ohne Verschmierungseffekte scharf aufgelöst, und es treten keine Oszillationen in der numerischen Lösung auf. Die Feinstrukturen werden sehr genau approximiert, weil die konvektiven Anteile mit einem expliziten Verfahren und kurzen Zeitschritten diskretisiert werden. Dadurch werden die hohen Frequenzen nicht so stark wie bei der Verwendung von impliziten Verfahren gedämpft. Da Sättigungs- und Druckgleichung entkoppelt sind und die Sättigungsgleichung voll explizit gelöst wird, ist das benutzte Verfahren trotz der einschneidenden CFL-Bedingung schneller als die Verfahren, die den gekoppelten Ansatz als Grundlage haben. Die entwickelte Methode ist auch im stark nichtlinearen Fall den Verfahren überlegen, die den gekoppelten Ansatz benutzen. Während in anderen Methoden numerische Probleme mit dem Newton-Verfahren auftauchen, verhält sich unser Verfahren sehr stabil. Hierbei wirken sich die bei unserem Verfahren benutzten relativ kleinen Zeitschriften positiv aus. Das Zeitschrittverfahren ist so variabel, daß es auch für schwach kompressible Fluide einsetzbar ist. Zudem kann die Ordnung des Zeitschrittverfahrens erhöht werden. Bei den Beispielen wurde jedoch wegen der schon geringen Zeitschrittweite darauf verzichtet.