Berechnungsschema von Speichern
Die Simulation von Speichern erfolgt mit einem neuentwickelten Baustein zur Berechnung von Speichern, deren Prozessfunktionen bereichsweise linear abzubilden sind. Der Baustein ist eine Weiterentwicklung des Ansatzes von /OSTROWSKI, 1991/ und wird detailliert bei /MEHLER, 1995/ beschrieben. Er ermöglicht die simultane Lösung der Kontinuitätsgleichung für mehrere Prozesse ohne aufwendige Iterationen und wird nachfolgend kurz erläutert.
Für einen Speicher, dessen Inhalt von mehreren Zu- bzw. Ablaufprozessen abhängig ist, kann die Kontinuitätsgleichung wie folgt dargestellt werden:
Die Entnahmeterme sind in der Regel nichtlineare Funktionen des Speicherinhaltes (z.B. die Entnahme aus dem Bodenspeicher mit den Prozessfunktionen) Diese Funktionen werden bereichsweise linearisiert.
Für jede Entnahmefunktion kann nach Lineariserung eine Geradengleichung aufgestellt werden, die nur noch vom Speicherinhalt abhängig ist. Die Steigung "m" der Geraden ändert sich von Stützstelle zu Stützstelle. Somit gibt es für jede vom Speicherinhalt abhängige Funktion einen bereichsweise linearisierten Verlauf entlang der Speicherfüllung. Die Funktion selbst kann mit einem für jeden Zeitschritt konstanten Faktor (A) skaliert werden, der alle weiteren Abhängigkeiten als Produkt zusammenfasst.
Die Kontinuitätsgleichung kann nun umformuliert werden zu:
Nach Ausmultiplizieren der Klammer wird die Kontinuitätsgleichung zu:
Diese Gleichung ist eine inhomogene lineare Differentialgleichung erster Ordnung und besitzt folgende Lösung:
Damit steht die Speicherfüllung zu jedem Zeitpunkt fest. Tritt innerhalb eines Zeitintervalls eine Bereichsüberschreitung ein, sind die Größen C1 und C2 mit den jeweils aktuellen Steigungen und Achsenabschnittswerten der bereichsweise linearisierten Funktionen neu zu berechnen. Die simultane Berechnung der Abgabenfunktionen wird durch Einsetzen der Speicherinhaltsgleichung in die jeweilige Geradengleichung erreicht. Allgemein ausgedrückt gilt für die mittlere Intensität aller Abgaben:
Mit diesem Berechnungsschema können alle Speicher, deren Prozesse bereichsweise linear zu beschreiben sind, berechnet werden. In Talsim-NG werden mit diesem Baustein die Bodenprozesse, die Speicher sowie die Transportstrecken berechnet.