BESS-Containerdeflagration: Abschätzung der Druckentlastungsfläche

Die Berechnung setzt das Niederdruckverfahren für Gase nach NFPA 68:2018, §§7.2.1, 7.2.6 und 7.4, bei einem anfänglichen Absolutdruck von 100 kPa um. Sie verwendet den allgemeinen Gemischkoeffizienten und nicht die vereinfachte Kohlenwasserstoffbeziehung $C=0.0223lambda S_u$.

Für die Stoffmengenanteile $x_i$ beträgt der Sauerstoffbedarf $nu_{O_2}=0.5(x_{H_2}+x_{CO})+2x_{CH_4}+3x_{C_2H_4}$ und es gilt $x_{st}=[1+nu_{O_2}/0.2095]^{-1}$. Die ideale Gasdichte und die Schallgeschwindigkeit werden für das Gas-Luft-Gemisch beim Äquivalenzverhältnis der gemessenen maximalen Brenngeschwindigkeit bestimmt. Die Massenstromdichte bei Schallgeschwindigkeit wird mit idealer isentroper Strömung berechnet.

Mit $p_0=100,mathrm{kPa}$ und $p_{max}$ als Überdruck lautet die vollständige Gleichung für geringe Trägheit

Der Solver aktualisiert den hydraulischen Durchmesser der quadratischen Entlastungselemente $D_v=sqrt{A_{v0}/n}$, die Reynolds-Zahl-abhängige Flammenverstärkung durch Flamme und Entlastungsöffnung, den Hindernisfaktor und die Korrektur des Längen-Durchmesser-Verhältnisses, bis die relative Flächenänderung unter $10^{-11}$ liegt (höchstens 32 Iterationen). Anschließend werden die Trägheitsschwelle der Entlastungselemente und die Korrekturen für translatorische bzw. schwenkbare Bewegung bestimmt. Derselbe Bemessungskern wird für die Bemessung zum Zieldruck, die Druckvariationen und das durchgerechnete NIST-Referenzbeispiel verwendet. Die Druckrückrechnung durchsucht jedes zusammenhängende gültige Intervall separat mit 45 Bisektionsschritten. Eine Nullstelle im stetigen Bereich muss die Beziehung für die vorgeschlagene Fläche erfüllen; bei einem übersprungenen ungültigen Intervall aufgrund der Ansprechdruckgrenzen oder einem abwärts gerichteten Sprung der Trägheitskorrektur wird der erste zulässige Druck als ausdrücklich gekennzeichnete obere Schranke ausgegeben, einschließlich des tatsächlichen Flächenresiduums ungleich null.

Henriksen et al. (2021), Tabellen 2, 5 und 7 (attached s6 pp.2,5), liefern die beiden zusammengehörigen Zusammensetzungen und Messdatensätze. Die Bemessung verwendet die tabellierten Maximalwerte; $K_G=(dp/dt)_{max}V_{test}^{1/3}$ mit $V_{test}=0.020,mathrm{m^3}$. Die PCHIP-Interpolation stellt die gemessenen Brenngeschwindigkeitskurven dar, ohne ein kinetisches Mischungsgesetz zu unterstellen. Molkov–Bragin (2015) gilt speziell für Wasserstoff und wird zur Bewertung der Anwendbarkeit herangezogen, nicht zur Übertragung von Wasserstoffkoeffizienten auf gemischte Batteriegase.

Die Verifikation umfasst NIST TN 2183 (2021), Anhang A3: Der im Rechner eingesetzte Solver liefert 0,921 m² gegenüber dem dort gerundeten Ergebnis von 0,9 m² für schwenkbare Entlastungselemente, bei einem wirksamen Trägheitsfaktor von 1,588. Auch die eingebaute Wand mit 1 m² hält bei Anwendung der im Rechner eingesetzten Druckrückrechnung die publizierte Druckgrenze von 6,9 kPa ein. Im Anhang nicht angegebene nominelle Transporteigenschaften werden im Referenzfall ausdrücklich festgelegt; geprüft wird das publizierte gerundete Bemessungsbeispiel, kein messtechnisch erfasstes Ereignis in einem Batteriecontainer.

• Ein einzelnes, nicht mit anderen Räumen verbundenes rechteckiges Gehäuse, anfangs bei 100 kPa Absolutdruck; die publizierten Gasfälle gelten bei 300 K. Benutzerdefinierte Messwerte müssen zur eingegebenen Anfangstemperatur und zu diesem Druck passen.
• Das gesamte Gehäusevolumen ist beteiligt; keine Minderung durch das von Racks verdrängte Volumen, einen begrenzten Gasvorrat, Verdünnung durch die Lüftungsanlage oder eine Teilwolke.
• Gleich große quadratische freie Öffnungen im gewählten Dach oder in einer einzelnen Stirnwand, ungehinderte Ausströmung ins Freie, $C_d=0.70$, keine Entlastungskanäle oder Flammendurchschlagsicherungen und keine Behinderung der Bewegung der Entlastungselemente. Eine Prüfung der Anordnung in Reihen berücksichtigt den eingegebenen Randabstand zwischen den Öffnungen und am Umfang der Montagefläche; halten Sie mit diesem Abstand Fläche für Rahmen, Rippen und andere nicht nutzbare Bereiche frei. Dies ist eine vorläufige Prüfung des Platzbedarfs, keine Detailauslegung der Entlastungselemente.
• Die Bewertung des Flammenwegs erfolgt getrennt von der Montageart. Die gesamte Gehäuselänge ist auch bei Dachmontage der konservative Standardansatz; ein kürzerer zugeordneter Entlastungsabschnitt erfordert eine dokumentierte Bewertung der Verteilung der Entlastungsöffnungen. Die schematische Darstellung kennzeichnet die gewählte Montagefläche, ohne einzelne Positionen vorzugeben.
• Die anfängliche Luftgeschwindigkeit beträgt höchstens 5 m/s; keine Zündung durch einen Flammenstrahl aus einem anderen Gehäuse. Die Hindernisoberfläche dient als Ersatzgröße zur Vorabschätzung und bildet die Rackanordnung nicht räumlich aufgelöst ab.
• Isentropenexponenten und Viskosität sind ausdrücklich eingegebene Werte zur Vorabschätzung; ersetzen Sie diese durch Transport- und thermochemische Daten des Gemischs. Die molaren Massen beruhen auf Standardatomgewichten. Reaktionschemie der Wärmefreisetzung und laminare Brenngeschwindigkeit werden nicht durch lineare Mischung abgeleitet.

Verfahrensgrundlage geprüft am 2. November 2026. Nur eine Abschätzung auf Grundlage des genannten Verfahrens von 2018; daraus wird keine Konformität mit einer späteren Ausgabe abgeleitet. Prüfen Sie anhand der geltenden Norm. Beschränken Sie die Anwendung auf $L/Dle5$, $S_u<3,mathrm{m/s}$, einen Überdruck im geschlossenen Prüfgefäß unter 10 bar, eine Masse der Entlastungselemente je Öffnungsfläche von höchstens 40 kg/m² und ein Gehäusevolumen über 1 m³. Dieses Werkzeug zur Vorabschätzung begrenzt die Hindernisoberfläche außerdem auf das Doppelte der Gehäuseinnenfläche; dichtere Einbauten erfordern eine validierte Analyse.

Die Ansprechdruckgrenzen werden bei jedem reduzierten Explosionsdruck angewendet: Bei höchstens 10 kPa muss der Ansprechdruck mindestens 2,4 kPa niedriger liegen; oberhalb von 10 kPa muss er strikt unter 75 % des reduzierten Explosionsdrucks liegen. Dadurch können voneinander getrennte gültige Intervalle entstehen. Für ungültige Zielwerte wird kein Flächenergebnis ausgegeben; Druckvariationen lassen ungültige Drücke aus. Die obere Grenze beträgt 50 kPa. Außerhalb des gültigen Bereichs sind der Druck und das zugehörige Rückrechnungsresiduum nicht verfügbar.

Die empirische Trägheitsschwelle der Entlastungselemente erzeugt eine Unstetigkeit. Liegt die vorgeschlagene Fläche innerhalb dieses Sprungs, existiert keine exakte Rückrechnung: Die Druckausgabe ist eine konservative obere Schranke am ersten zulässigen Ast, und die Bewertung benennt die Trägheitsschwelle. Eine durch die Ansprechdruckgrenzen entstehende Lücke liefert entsprechend eine gekennzeichnete obere Schranke auf dem gültigen oberen Ast. Das Flächenresiduum einer Schranke bleibt erhalten und muss nicht null sein. Innerhalb einer Lücke wird weder eine Interpolations- noch eine Bisektionsnullstelle ausgewiesen; diese Schranken beschreiben die empirische Bemessungsbeziehung, keinen simulierten Öffnungsverlauf der Entlastungselemente.

Die publizierten $P_{ex}$-Tabellen kennzeichnen den Druckbezug nicht ausdrücklich. Die Werte werden hier konservativ als Überdrücke behandelt; falls es sich um Absolutdrücke handelt, erhöht diese Interpretation den Schätzwert. Bestätigen Sie den Druckbezug für einen konkreten Gasprüfbericht. Auch der Fließtext der Quelle weicht bei den Maxima der Brenngeschwindigkeit und der Druckanstiegsgeschwindigkeit von den Tabellen ab; die Tabellenwerte werden beibehalten. Bei der späteren Erhöhung aufgrund von Hindernissen im Beispiel aus NFPA-Anhang I scheint der angegebene Term −0,2 zu fehlen; die Prüfung der Übernahme dieses Beispiels verwendet daher dessen konsistente erste Iteration. Die Verifikation des Endergebnisses verwendet stattdessen das unabhängige NIST-Beispiel mit schwenkbaren Entlastungselementen.

Für die Druckvariation werden je zusammenhängendem Intervall bis zur eingegebenen Gehäusedruckfestigkeit 121 Punkte berechnet oder jeweils 81 Punkte auf jeder Seite einer Trägheitsschwelle. Die gezeichneten Linien werden an Ansprechdrucklücken und Trägheitssprüngen unterbrochen. Die Flammenkurven verwenden 141 PCHIP-Punkte über die gemessenen Brennstoff-Luft-Verhältnisse. Damit sind weder CFD, transiente Öffnungsdynamik der Entlastungselemente, eine äußere Druckspitze, Detonation, eine Berechnung der Tragfähigkeit noch ein Teilvolumenmodell abgedeckt. Für BESS-Anlagen mit dichten Einbauten, das Strukturverhalten, die sichere Ausströmung, Rückhaltevorrichtungen und die Systemintegration sind qualifizierte ingenieurtechnische Planung und repräsentative Prüfungen erforderlich.