EBSILON®Professional Online Dokumentation
|
Leitungsanschlüsse |
|
|
|
1 |
Systemgrenze mit Wertvorgabe |
|
Zur Spezifikation der Stoffzusammensetzung und Stoffeigenschaften (Heizwert, Dichte, ...) siehe ”Stoffwerte”.
Allgemeines Stoffwerte Gleichungen Bauteilformen Beispiel
Bauteil 1 wird verwendet, um eintretende Ströme an den Schnittstellen des Systems zu definieren. Beispiele sind:
Der Mediumzustand ist durch Festlegung der folgenden Werte gekennzeichnet:
Bei einer Vorgabe der Werte (P,H, M) außerhalb des Gültigkeitsbereichs der Start- und Grenzwerte (ModellEinstellungen --> Simulation ) gibt es eine Fehlermeldung.
Wenn H gegeben ist, wird die Temperatur berechnet. Wenn auch T gegeben ist, zeigt Ebsilon eine Fehlermeldung an.
Drei der fünf Werte reichen aus, um den Mediumszustand eindeutig zu beschreiben. Üblicherweise sind P, T und M vorgegeben. H kann mit P und T berechnet werden. Q ergibt sich dann aus H und M.
Wenn der Energiestrom Q gegeben ist und der Massenstrom M gegeben ist, wird die Enthalpie H berechnet.
Mitunter ist die Vorgabe von allen drei Zuständen mit Bauteil 1 nicht notwendig, da beispielsweise einzelne Werte bekannt sind und bereits an anderen Stellen des Systems festgelegt sind. Würde ein solcher Wert nochmals mit Bauteil 1 eingestellt, wäre eine Uneindeutigkeit die Folge.
Mit dem Vorgabewert LOAD wird der im Vorgabewert M eingetragene Wert multipliziert, um den Massenstrom auf einer Leitung zu setzen. Dies ist hilfreich, wenn man in Unterprofilen verschiedene Lastfälle durchrechnen möchte. Man erspart sich dann die Berechnung des Massenstroms im jeweiligen Profil und braucht beim Kopieren oder Ändern der Vorgabe nur noch M im Design-Profil zu ändern.
Der Standardwert für LOAD ist 1.0, und solange man in allen Profilen LOAD auf 1.0 lässt, hat er keine Auswirkungen.
ACHTUNG:
Für Dampf im 2-Phasen Gebiet muss H eingegeben werden, da P und T nicht ausreichend sind, um H eindeutig festzulegen.
Grundlagen:
Bisher wurde in Ebsilon davon ausgegangen, dass alle Fluide langsam genug strömen, so dass die kinetische Energie vernachlässigbar ist und totale und statische Zustandsgrößen gleichgesetzt werden können. Es wurde deshalb nicht zwischen totalen und statischen Werten unterschieden.
Einzige Ausnahme war die Dampfturbine (Bauteil 122). Da hier entsprechend hohe Strömungsgeschwindigkeiten auftreten, wurden bei der Berechnung der Dampfturbine die kinetischen Anteile berücksichtigt. Dabei wurde angenommen, dass die auf den Leitungen gespeicherten Größen grundsätzlich totale Größen sind.
In Release 15 besteht nun die Möglichkeit, überall zwischen totalen und statischen Größen zu unterscheiden. Dabei wird die bisher bei der Dampfturbine getroffene Annahme, dass die Leitungswerte totale Größen sind, auf alle Leitungen angewandt.
Die Umrechnung zwischen totalen und statischen Größen erfordert die Kenntnis der Strömungsgeschwindigkeit vel (da die Abkürzung v in Ebsilon bereits für das spezifische Volumen verwendet wird, steht v für die Geschwindigkeit nicht zur Verfügung).
In den Ebsilon-Standardeinheiten (vel in m/s, H in kJ/kg) gilt:
Htot = Hstat + Hkin, mit
Hkin = 0.0005 * vel²
Die Strömungsgeschwindigkeit kann dabei wahlweise
Die Vorgabe erfolgt beim Randwert (Bauteil 1) oder Startwert (Bauteil 33) durch einen der Vorgabewerte VEL_SET, A_SET oder D_SET.
Als Ergebnis werden dann alle drei Größen VEL, A und D angezeigt.
Zusätzlich kann in Release 15 auch die potentielle Energie berücksichtigt werden. Für diese gilt in Ebsilon-Standardeinheiten (z in m, H in kJ/kg):
Hpot = 0.001 * g * z,
wobei g = 9.81 m²/s die Fallbeschleunigung ist und z die Höhe, bezogen auf ein gewähltes Nullniveau. Die Höhe wird in Bauteil 1 bzw. 33 durch einen Vorgabewert Z_SET spezifiziert und als Ergebniswert Z angezeigt.
Insgesamt ergibt sich
Htot = Hstat + Hkin + Hpot
Die Unterscheidung zwischen totalen und statischen Größen erfolgt nur innerhalb des Bauteils 1 bzw. 33. Es ist nicht sinnvoll, die Strömungsgeschwindigkeit als Attribut der Leitung zu betrachten, da sich der Leitungsquerschnitt und somit die Strömungsgeschwindigkeit von einer Stelle der Leitung zur anderen ändern kann. Die Vorgabe und Berechnung muss deshalb stets gemeinsam mit der Vorgabe der Geschwindigkeit (bzw. Querschnitt oder Durchmesser) erfolgen.
Um unterschiedliche Strömungsgeschwindigkeiten an verschiedenen Stellen der Leitung betrachten zu können, besteht die Möglichkeit, mehrere Bauteile 33 auf einer Leitung zu setzen. Eine Vorgabe statischer Größen ist allerdings nur in einem Bauteil möglich.
Wenn die Strömungsgeschwindigkeit bzw. die Höhe bekannt ist, werden in Bauteil 1 folgende Ergebniswerte angezeigt:
Bauteil 1 bzw. 33 ermöglicht nicht nur die Anzeige, sondern auch die Vorgabe statischer Zustandsgrößen. Dadurch können in Wärmeschaltbildern angegebene oder auch gemessene statische Größen jetzt so vorgegeben werden.
Hierzu muss – außer im letzten Fall - die Strömungsgeschwindigkeit wie oben beschrieben entweder direkt oder über Querschnitt oder Durchmesser vorgegeben sei.
Folgende Möglichkeiten sind implementiert:
Dieser Fall ist insbesondere interessant, wenn mit diversen Druckmessgeräten (Prandlt-Staurohr, Pitot-Rohr) statische und totale Drücke gemessen werden, um dadurch den Durchsatz zu bestimme
Vorgabe von Zusammensetzungen für Luft / Rauchgas als Medium:
Mit Bauteil 1 oder 33 ist es möglich, die Zusammensetzung eines Fluids zusätzlich zu den P/T/H/M/Q Werten vorzugeben
Die Spezifikation der Werte der Analyse und damit in Verbindung stehender Attribute (Heizwert, Dichte, ...) ist im Abschnitt "Vorgabe von Stoffeigenschaften" beschrieben.
Vorgabe von Zusammensetzungen :
Zusammensetzungen können über Rand- bzw. Startwerte (Bauteil 1 bzw. 33) oder über Messwerte (Bauteil 46) vorgegeben werden. Die Vorgabe über Messwerte ist dann erforderlich, wenn eine Validierung der Zusammensetzungen vorgenommen werden soll.
Das Zusammenspiel zwischen den beiden Vorgabemöglichkeiten war relativ kompliziert. Da die Summe aller Substanzen 1 ergeben muss, konnten nicht alle Substanzen über Messwerte vorgegeben werden, sondern man brauchte noch einen Freiheitsgrad für die Normierung, wozu in der Regel der größte Anteil im Bauteil Rand- bzw. Startwert genommen wurde.
Ab Release 11 kann man die vollständige Zusammensetzung entweder über einen Rand- bzw. Startwert oder über Messwerte vorgeben:
Bei der Vorgabe der vollständigen Zusammensetzung über einen Rand- bzw. Startwert, gibt es keine Änderungen in der Spezifizierung. In diesem Bauteil kann unmittelbar geprüft werden, ob die Summe alle Substanzen 1 ergibt. Falls dies nicht der Fall ist, wird eine Fehlermeldung generiert.
Bei der Vorgabe über Messwerte ist für jede vorkommende Substanz ein Messwert auf die Leitung zu platzieren. Auf derselben Leitung wird außerdem ein Rand- bzw. Startwert benötigt, um die Zusatzinformationen abzusetzen, die nicht über Messwerte definiert werden können (zum Beispiel der Kohletyp oder cp-Koeffizienten beim benutzerdefinierten Fluid). In diesem Startwert ist außerdem anzugeben, dass bei Integration der Materialgleichungen in die Gleichungsmatrix für die nicht spezifizierten Substanzen Gleichungen mit dem Wert 0 abgesetzt werden (andernfalls müsste man für jede Substanz, die nicht enthalten ist, einen Messwert mit dem Wert 0 setzen).
Falls die Summe aller Messwerte nicht 1 ergibt, wird eine Fehlermeldung ausgegeben.
Darstellung der Stoffanteile für die Bauteil 1 und 33 auf dem Blatt "Vorgabewerte" :
Bei den Bauteilen zur Vorgabe einer Zusammensetzung wurde die Möglichkeit geschaffen, diese (und sonstige Koeffizienten) nicht nur in den Dialogen auf dem Blatt „Stoff-Anteile“ einzutragen, sondern alle Werte auf dem Hauptblatt (Vorgabewerte) darzustellen. Dies bietet zwar weniger Komfort, ist aber unter Umständen hilfreich insbesondere in Zusammenhang mit der Excel-Schnittstelle.
Die Aktivierung dieser Option erfolgt in der Auswahlbox „Ansichtseinstellungen“ unterhalb der Komponentengrafik durch den Eintrag „Zusammensetzung auch als Vorgabewerte anzeigen“.
Beispiel : Vorgabewerte (mit Stoffanteilen)
- Erste Eigenschaftenseite "Vorgabewerte"
- Letzte Eigenschaftenseite "Vorgabewerte"
Hinweis: Zugriff auf XH2OG und XH2OL :
Durch die Zusammenlegung von XH2OG und XH2OL zu XH2O (entsprechend auch XNH3L und XNH3G zu XNH3 sowie XCO2 und XCO2L zu XCO2, siehe Kapitel 1.2.5.1) sind Zuweisungen auf die alten Variablen bei den Vorgabewerten (Bauteil 1 und 33) nicht mehr möglich. Bei den Ergebniswerten (Leitungsergebnisse) stehen die alten Bezeichnungen jedoch weiterhin zur Verfügung.
Berücksichtigung von nicht gasförmigen Bestandteilen beim spezifischen Volumen von Gasen
Für die Berechnung des spezifischen Volumens (und damit auch der Dichte) von Gasen (Luft, Rauchgas, Gas, Rohgas) wurden nur die gasförmigen Bestandteile berücksichtigt, da der Anteil der flüssigen und festen Bestandteile zum spezifischen Volumen wegen der höheren Dichte in der Regel unerheblich ist.
Das spezifische Volumen dieser Bestandteile kann in der Regel nicht berechnet werden, da üblicherweise hierfür nur die elementare Zusammensetzung bzw. die allgemeine Spezifikation „Asche“ vorgegeben wird.
Es wurde die Möglichkeit geschaffen, für diesen Anteil eine Dichte im Vorgabewert „Dichte für Elementarzusammensetzung“ (RHOELEM) vorzugeben. Bis Release 10 wurde die Dichtevorgabe (Vorgabewert RHO) nur für Öl verwendet. Zur Vereinheitlichung ist auch bei Öl RHOELEM zu verwenden. RHO steht als Vorgabewert nicht mehr zur Verfügung.
Als Ergebniswerte auf den Leitungen gibt es sowohl RHO (mittlere Dichte des Gesamtstroms) als auch RHOELEM (Dichte des durch die Elementaranalyse gegebenen Anteils).
Wird bei der Dichte (RHOELEM) ein Wert von 0 eingetragen, wird der Anteil der als Elementaranalyse gegebenen Stoffe bei der Ermittlung des spezifischen Volumens vernachlässigt.
Bei nicht gasförmigen Bestandteilen, deren chemische Zusammensetzung bekannt ist, wird das spezifische Volumen aus den entsprechenden Stoffdaten ermittelt. Dies betrifft flüssiges H2O, NH3 und CO2, für die Bibliotheken in Ebsilon integriert sind, sowie die neuen Substanzen für die Direktentschwefelung, für die folgende Konstanten verwendet werden:
• CaSO4 2960 kg/m³
• CaCO3 2730 kg/m³
• CaO 3370 kg/m³
• Ca(OH)2 2240 kg/m³
• MgCO3 2960 kg/m³
• MgO 3580 kg/m³
Falls in der Kohle gebundenes Wasser (H2OB) vorhanden ist, wird dieses als Bestandteil der Kohle angesehen, d.h. es wird davon ausgegangen, dass dieser Anteil bereits in RHOELEM enthalten ist. H2O auf einer Kohleleitung (beispielsweise Regenwasser zwischen den Kohlestückchen) wird dagegen separat mit den Stoffdaten für H2O berechnet.
Die Bauteile 1 bzw. 33 können auch für die Vorgabe von binären Mischungen und für das Universalfluid verwendet werden.
Bei binären Mischungen besteht die Auswahl zwischen Ammoniak / Wasser oder Wasser / Lithiumbromid als Arbeitsmedium. Außerdem ist der Massenanteil XI des Kältemittels (Ammoniak bei Ammoniak / Wasser, Wasser bei Wasser / Lithiumbromid) vorzugeben.
Beim Universalfluid sind in der Vorgabetabelle zunächst die zu verwendenden Stoffdaten-Bibliotheken vorzugeben. In der Spalte „Spezifikation“ ist dann die jeweils gewünschte Zusammensetzung und weitere Parameter einzugeben. Die hier einzutragende Zusammensetzung bezieht sich auf die Massenanteile innerhalb des Teilstroms, der von der jeweiligen Bibliothek berechnet wird.
Welche Stoffe zur Auswahl stehen und welche zusätzlichen Eingabeparameter benötigt werden, hängt von der jeweiligen Bibliothek ab.
Bei FDBR ist als zusätzlicher Eingabeparameter ein Berechungstyp (cp-type) einzugeben. Dies ersetzt die Auswahl des Leitungstyps bei den klassischen Ebsilon-Leitungen und (bei Kohleleitungen) des Kohletyps. Es stehen folgende Werte für den cp-type zur Verfügung:
In Abhängigkeit von der hier gewählten Einstellung werden die übrigen Parameter verwendet. Es ist zu beachten, dass beim benutzerdefinierten Fluid für das cp-Polynom nicht die globalen Parameter aus den Modelleinstellungen, sondern die hier eingetragenen Werte verwendet werden. Auf diese Weise können auch unterschiedliche benutzerdefinierte Fluide innerhalb einer Schaltung verwendet werden.
Bei Verwendung der LibHuAir_xiw für feuchte Luft ist als Parameter xi der Wassergehalt (Massenanteil Wasser pro Gesamtmasse vorzugeben.
Hinweis: Die ebenfalls von der Hochschule Zittau/Görlitz entwickelte LibHuAir verwendet im Gegensatz zur LibHUAir_Xiw als Eingangsparameter die absolute Feuchte (Massenanteil
Wasser pro Masse trockener Luft) und liefert als Ergebnis die Enthalpie und Entropie bezogen pro Masse trockener Luft. Diese Bibliothek ist jedoch nicht in Ebsilon eingebunden, da die LibHUAir_Xiw besser in die Philosophie von Ebsilon passt.
Vorgabe der relativen Luftfeuchte:
Dies ist im Randwert (Bauteil 1) und im Startwert (Bauteil 33) möglich.
Dabei ist zu beachten, dass
Bei Verwendung der LibAmWa für Ammoniak / Wasser - Gemische ist als Parameter xi der Ammoniakgehalt (Massenanteil Ammoniak pro Gesamtmasse) vorzugeben.
Bei Verwendung der LibWaLi für Wasser / Lithiumbromid - Gemische ist als Parameter xi der Wassergehalt (Massenanteil Wasser pro Gesamtmasse) vorzugeben.
Bei Verwendung der LibH2 für Wasserstoff ist zu beachten, dass die thermischen Eigenschaften des Wasserstoffs bei tiefen Temperaturen von der relativen Ausrichtung der Kernspins im H2-Molekül abhängen. Beim Parawasserstoff sind die Spins entgegengesetzt ausgerichtet, beim Orthowasserstoff parallel. Normalerweise liegt eine Mischung beider Varianten („Normalwasserstoff“) vor. Parawasserstoff ist insbesondere für die Kernspintomatographie interessant. Über ein Flag kann wahlweise Normal- oder Parawasserstoff eingestellt werden.
Bei Verwendung der Refprop-Bibliothek gibt es Umschaltmöglichkeiten für den Referenzzustand. Standardmäßig verwendet Ebsilon den von Refprop als „Default“ vorgegebenen Referenzzustand. Alternativ werden angeboten:
Bei der Vorgabe der Nullpunkte ist darauf zu achten, dass diese in den Standardeinheiten der REFPROP erfolgen muss:
Eine Einheitenumrechnung ist an dieser Stelle nicht möglich.
Im Feld "Berechnungsmodell" kann man unter verschiedenen Berechnungsmodellen auswählen, die die REFPROP-Bibliothek bietet :
Der GERG-2008 Mechanismus wird für Erdgas-Gemische standardmäßig verwendet, nicht allerdings für Reinstoffe. Für die Reinstoffe sind die in REFPROP umgesetzten Gleichungen genauer, da sie nicht (wie in der GERG-2008) auf Grund von Rechengeschwindigkeit abgekürzt wurden. Die Unterschiede sind jedoch äußerst gering (solange sich kein Wasser im Gemisch befindet). Die GERG-2008 (programmiert in Bochum von Wagner und Kunz) weist aber in jedem Fall höhere Konsistenz und Genauigkeiten im Phasengleichgewicht auf.
Um REFPROP-Fluide einzufügen, klickt man in das Feld "Zum Einfügen anklicken..." in der Spalte Material der Tabelle Zusammensetzung rechts im Fenster "Universalfluid-Zusammensetzung-Refprop". Es erscheint dann eine Liste aller verfügbare REFPROP-Fluide, aus der man einen Eintrag auswählen kann. In den folgenden Zeilen kann man dann weitere Komponenten auswählen. Die so gewählten Komponenten werden von der REFPROP-Bibliothek als Gemisch gerechnet. Man sollte sich dabei auf wenige Komponenten beschränken, da die Rechenzeit mit der Anzahl der Komponenten im Gemisch stark ansteigt. Leider sind auch nicht alle Kombinationen möglich, man erhält dann eine entsprechende Warnmeldung (z.B. "Binäre Wechselwirkungsparameter für diese Mischung nicht verfügbar, ideales Lösungsverhalten angenommen").
Als Alternative zu den selbst zusammengestellten Mischungen mit frei einstellbaren Stoffanteilen kann man auch fertig konfigurierte Mischungen der REFPROP-Bibliothek verwenden. Einige dieser Gemische werden als "pseudo pure fluid" in derselben Liste wie die reinen Stoffe angezeigt. Diese werden auch wie reine Stoffe berechnet, mit effektiven (gemittelten) Stoffparametern. Um auf andere vordefinierte Gemische zugreifen zu können, muss man mit der rechten Maustaste in die Spalte "Material" klicken und im Kontextmenü den Eintrag "Standard-Zusammensetzung einfügen..." auswählen. Dann erhält man die fertig konfigurierten Gemische zur Auswahl.
Bei allen Substanzen ist zusätzlich die Angabe eines cp-Korrekturfaktors möglich. Dies ist ein konstanter Faktor, mit dem die spezifische Wärmekapazität cp und damit auch die Enthalpie H und die Entropie S multipliziert werden. Dieser kann für den Fall genutzt werden, dass die genauen Stoffdaten für einen bestimmten Stoff unbekannt sind, man jedoch einen ähnlichen Stoff kennt, dessen thermodynamische Eigenschaften durch Anpassung der spezifischen Wärmekapazität in einem Punkt genau genug beschrieben werden können.
Man beachte, dass dieser Korrekturfaktor nichts mit dem Vorgabewert CPCORR zu tun hat. CPCORR bezieht sich ausschließlich auf den Asche-Anteil und wird bei FDBR-Fluiden verwendet (zum Beispiel auf Kohle- und Rauchgas-Leitungen). Der hier beim Universalfluid verfügbare cp-Korrekturfaktor bezieht sich dagegen auf die komplette Zeile in der Stoffanteil-Tabelle.
Beim Fluidtyp "Thermoflüssigkeit" stehen diverse Thermoöle und Salzschmelzen zur Auswahl. Es besteht auch die Möglichkeit, Stoffwertkoeffizienten selbst vorzugehen. Die Bedeutung der Koeffizienten ist im Abschnitt Thermoflüssigkeiten im Kapitel "Erläuterungen zur Berechnung - Stoffwerte" beschrieben.
Für die Berechnung des Heizwertes (NCV/GCV) eines brennbaren Fluids wurde ein Schalter (FNCV) eingebaut, mit dem man festlegen kann, was bei einer Änderung der Zusammensetzung passieren soll:
In früheren Versionen war es erforderlich, nach einer Änderung der Stoffwertzusammensetzung den Knopf „Berechneten Wert übernehmen“ zu betätigen, um den Heizwert anzupassen. Da dies leicht vergessen werden konnte, wurde dadurch eine potentielle Fehlerquelle vermieden. Die (häufig lästige) Warnung einer Abweichung zwischen vorgegebenem und berechnetem Heizwert ist deshalb zum Kommentar heruntergestuft worden.
Ebsilon bietet neue Möglichkeiten zur Vorgabe von Heizwerten.
Es besteht die Möglichkeit, auch den oberen Heizwert (GCV) vorzugeben bzw. sich anzeigen zu lassen (siehe dazu auch Kapitel "Vorgabe von Stoffeigenschaften"). Die Vorgabe erfolgt im Bauteil 1 (Randwert), 33 (Startwert) oder Bauteil 46 (Messwert), für die Anzeige steht ein Leitungswert GCV, aber auch das Bauteil 45 (Wertanzeige) zur Verfügung.
Dazu wurde das Bauteil 1 (Randwert) bzw. Bauteil 33 (Startwert) um folgende Vorgabewerte erweitert
(siehe dazu Beispiel : Vorgabewerte (mit Stoffanteilen - Erste Eigenschaftenseite "Vorgabewerte" ):
Diese Erweiterung der Vorgabewerte wird sichtbar, wenn die Darstellung der Stoffanteile für die Bauteile 1 und 33 auf dem Blatt "Vorgabewerte" gewählt wird.
Hinweis:
Der Vorgabewert NCV (unterer Heizwert) wird nur bei entsprechender Einstellung des Schalters FNCVSRC verwendet.
Wenn man mit der rechten Maustaste auf einen Zahlenwert der Zusammensetzung klickt, gibt es zu den bisherigen Skalierungsmöglichkeiten zwei weitere Varianten:
Beim benutzerdefinierten Fluid erfolgt die Vorgabe der cp-Koeffizienten (ab Release 8) beim Bauteil 1 bzw. 33 und nicht mehr bei den Modelleinstellungen. Damit können auch verschiedene benutzerdefinierte Fluide in einer Schaltung verwendet werden. Wenn man solche Fluide allerdings mischen will, muss das Universalfluid verwendet werden. Andernfalls werden die cp-Werte der ersten Leitung verwendet und eine Warnung ausgegeben.
Heizwert, Zusammensetzung und ergänzende Attribute (wie Kohletyp, Z-Faktor beim Öl) konnten nur gemeinsam gesetzt oder nicht gesetzt werden. Damit war es zum Beispiel nicht möglich, mit einem Bauteil 33 nur den Heizwert vorzugeben und nicht gleichzeitig die Zusammensetzung. Es gibt eine Auswahlbox, mit der eingestellt wird, welche dieser Größen gesetzt werden soll. Dabei können auch zwei oder alle drei gesetzt werden.
Zur Steigerung der Genauigkeit kann eine Realgaskorrektur durchgeführt werden. Im Vergleich zur Berechnung als ideales Gas ergibt sich dadurch allerdings eine deutliche Erhöhung
der Rechenzeit.
Häufig befinden sich in der Schaltung jedoch nur einige wenige Leitungen, für die die Realgaskorrektur signifikant ist. Die Anwendung auf sämtliche Luft- und Rauchgasleitungen mit
Drücken im Bereich des Atmosphärendrucks ist in der Regel überflüssig. Aus diesem Grunde wurde die Möglichkeit geschaffen, die Realgaskorrektur leitungsspezifisch zu definieren.
Die Vorgabe der anzuwendenden Realgaskorrektur
erfolgt im Bauteil 1 bzw. 33 (Rand- bzw. Startwert) auf dem Blatt „Stoffanteile“.
Die Definition gilt dann für die jeweilige Leitung und wird dann entlang des Hauptstroms weitergereicht. Welche Korrektur verwendet wurde, ist auf der Leitung im Ergebniswert
FREALGC ersichtlich. Wenn bei einer Zusammenführung zwei unterschiedliche Realgaskorrekturen zusammentreffen, wird eine Warnung ausgegeben.
Die Konvergenzgenauigkeit ist in Ebsilon eine modellweite Einstellung. Sie ist eine obere Schranke für die zulässige relative Änderung einer Variablen (Massenstrom, Druck, Enthalpie) von einem Iterationsschritt zum nächsten. Erst wenn bei allen Variablen die relative Änderung betragsmäßig kleiner als diese Schranke ist, wird das Iterationsverfahren erfolgreich beendet.
Normalerweise wird diese relative Änderung auf den Wert der Variable bezogen. Wenn sich also beispielsweise ein Massenstrom von 50 auf 51 kg/s ändert, ist die relative Änderung 1%. Bei einer Konvergenzgenauigkeit von 10-7 (dies ist der Standardwert) darf sich dieser Massenstrom von 50 nur noch um 0.000005 kg/s (also 5 mg/s) ändern, damit der Wert als konvergent angesehen wird. Bei einem Massenstrom von 0.01 kg/s wäre die zulässige Änderung dann nur noch 10-9 kg/s. Solch kleine Änderungen sind aber in der Praxis nicht signifikant und würden die Iteration unnötigerweise in die Länge ziehen. In vielen Fällen könnte aufgrund des „numerischen Rauschens“ auch gar keine Konvergenz mehr erzielt werden.
Aus diesem Grunde wurde in Ebsilon eine minimale Bezugsgröße für die Berechnung der relativen Änderung definiert. Wenn der Wert der Variable kleiner als die Bezugsgröße war, wurde die relative Änderung nicht in Bezug auf die Variable, sondern in Bezug auf die Bezugsgröße berechnet. Dadurch wurden (bei einer Konvergenzgenauigkeit von 10-7) Massenstromschwankungen von weniger als 0.000002 kg/s, Druckschwankungen von weniger als 0.0000002 bar und Enthalpieschwankungen von weniger als 0.00006 kJ/kg nicht mehr als Konvergenzhindernis angesehen.
Bis Release 11 waren diese Mindestbezugsgrößen fest im Code verankert. Ab Release 12 besteht die Möglichkeit, sie individuell für einzelne Leitungen vorzugeben. Dadurch können weniger interessante Bereiche der Schaltung mit einer geringeren Genauigkeit gerechnet werden und dadurch die Rechenzeit möglicherweise verkürzt werden.
Die Vorgabe dieser Mindestbezugsgrößen erfolgte bis Release 12 im Bauteil 1 bzw. 33 (Rand- bzw. Startwert) in den Einträgen
MINREFITP
MINREFITH und
MINREFITM
Wenn man diese Werte leer lässt, verwendet Ebsilon die Standardwerte. Die Bezugsgrößen gelten dann für die jeweilige Leitung und werden dann entlang des Hauptstroms weitergereicht.
Bei einer Zusammenführung werden also die Bezugsgrößen vom Nebenanschluss (Anschluss 3) ignoriert. Dadurch wird es erleichtert, im Nebenzweig mit einer geringeren Genauigkeit zu rechnen, während man hinter der Zusammenführung mit dem Hauptstrom wieder die höhere Genauigkeit hat.
Hinweis:
Zusätzlich wurde die, bei Bauteil 1 und 33, implementierte Funktionalität zur individuellen Vorgabe der Konvergenzgenauigkeit in das Bauteil 147 (Begrenzer) verlagert. Aus Kompatibilitätsgründen bleibt sie in Bauteil 1 und 33 jedoch weiterhin verfügbar. Die dafür benötigten Eingabefelder MINREFITP, MINREFITH und MINREFITM wurden lediglich ausgeblendet, können bei Bedarf jedoch sichtbar geschaltet werden.
Auf Elektroleitungen können mit einem Bauteil 1
Hinweis: Da aus historischen Gründen in Ebsilon auf Elektroleitungen die Leistung als Grundvariable verwendet wird, wird bei einer Vorgabe der Stromstärke die Leistung berechnet und auf die Leitung geschrieben. Aus diesem Grunde ist eine gleichzeitige Vorgabe von Stromstärke und Leistung nicht möglich. Da zur Berechnung der Leistung die Spannung und die Phase benötigt werden, muss bei Vorgabe der Stromstärke auch Spannung und Phase im selben Bauteil 1 vorgegeben werden.
|
Wenn P gegeben, dann Wenn H gegeben, dann Wenn T gegeben, dann Wenn M gegeben, dann Wenn Q gegeben, dann |
Die Belegung der fünf Werte mit den Rand-Vorgabewerten erfolgt nur beim 2. Iterationsdurchgang.
falls H1<=0, , dann H1 = f(P1,T1) mit Vorgabe von P1 und T1
falls H1 > 0, , dann T1 = f(P1,H1) mit Vorgabe von P1 und H1
falls Q1 <= 0, , dann Q1 = M1 * H1
falls Q1 > 0, , dann H1 = Q1 / M1 falls M1 > 0
oder M1 = Q1 / H1 falls H1 > 0
 |
Bauteilform 1 |
Klicken Sie hier >> Bauteil 1 Demo << um ein Beispiel zu laden.
