Wechselstrom-Selbsterwärmung von Batterien bei −30 °C: durch Plating begrenzte Aufheizung

Der Frequenzbereichsansatz folgt Zhang et al. (2015), dem Optimierungskonzept mit konstanter Polarisation von Ruan et al. (2016) und der Elektrodenpotenzialbegrenzung von Ge et al. (2016). Liu et al. (2023), pp. 4–7, enthält die hier verwendeten, frei zugänglichen Gleichungen für Heizleistung und Anodenpotenzial bei sinusförmiger Anregung.

Für jede Elektrode gilt $Z_k=[R_{ct,k}+Z_{D,k}]/[1+Y_{dl,k}(R_{ct,k}+Z_{D,k})]$. Die planare Diffusion mit reflektierender Randbedingung wird wie in impedance.py durch $Z_D=R_D/[stanh(s)]$, $s=sqrt{jomegatau_D}$ beschrieben. Die CPE-Admittanz lautet $Y=C_*omega_* (jomega/omega_*)^a$ mit $omega_*=1,mathrm{s}^{-1}$; der Eingabewert in F ist eine äquivalente Kapazität bei dieser Kreisfrequenz. Ein separater Deckschichtzweig aus $R_{SEI}parallel$CPE sowie $R_0$ und $jomega L$ bilden $Z_{cell}$.

Die Arrhenius-Skalierung lautet $g(E,T)=exp[(E/R)(1/T-1/T_{ref})]$ mit absoluten Temperaturen und $T_{ref}=253.15$ K. Durchtritts- und Deckschichtwiderstände verwenden jeweils eigene Aktivierungsenergien. Sowohl der Diffusionswiderstandsparameter als auch die Diffusionszeit werden mit dem Diffusionsfaktor skaliert; dadurch bleibt die chemische Kapazität im Niederfrequenzbereich unverändert.

Der Scheitelstrom ist der kleinste Wert aus $I_{max}$, $(U_n-U_{reserve})/|Z_n|$, $V_{pol}/|Z_n+Z_p+Z_{SEI}|$ und $V_{AC}/|Z_{cell}|$. Damit bleibt das im Modell berechnete minimale Anodenpotenzial jeder Periode mindestens auf Höhe der Reserve. Der Anodenzweig enthält die Diffusion; seine Spannung umfasst daher sowohl kinetische als auch Konzentrationspolarisation. Die Spannungsreserve an den Klemmen muss relativ zur OCV angegeben werden.

Es gilt $P=I_{peak}^2operatorname{Re}(Z_{cell})/2=I_{rms}^2operatorname{Re}(Z_{cell})$ und $C_{th}dT/dt=P-H(T-T_{amb})$. Der Frequenzdurchlauf integriert $dt/dT=C_{th}/[P-H(T-T_{amb})]$ mit dem logarithmischen Integral jedes lokal linearen Abschnitts der Nettoheizleistung und verwirft Pfade mit Aufheizstillstand. Die logarithmische Form erfasst den starken Anstieg der Aufheizdauer bei gegen null gehender Nettoheizleistung, ohne die Endpunktsingularität einer Trapezintegration. Die Frequenzoptimierung innerhalb der vorgegebenen Grenzen verwendet eine adaptive lokale Verfeinerung desselben Integrals, ausgewertet am tatsächlichen temperaturabhängigen Arbeitspunkt.

Das gewählte manuelle, feste oder adaptive Steuerprofil wird mit Intervallhalbierung integriert, bis die lokale Differenz zwischen grober und feiner Diskretisierung eine relative Toleranz von $10^{-5}$ einhält. Ein zweiter Durchlauf beginnt mit doppelt so vielen Intervallen und einer auf ein Viertel verringerten Toleranz. Nicht hinreichend aufgelöste Fehler werden gemeldet. Bei adaptiver Steuerung wird der interpolierte Frequenzkandidat an jeder ausgewerteten Temperatur mit der zeitoptimalen Festfrequenz verglichen; verwendet wird die Variante mit der höheren Heizleistung.

Eine gemeinsame Arbeitspunktfunktion liefert die Werte für die thermische DGL sowie die dargestellten Verläufe von Heizleistung, Strom und Anodenpotenzial, auch bei Temperaturen unterhalb des Startrasters. Mit dem Dormand–Prince-Verfahren werden Temperatur, kumulierte erzeugte Wärme und kumulierte Wärmeabgabe an die Umgebung gemeinsam integriert. Bei erreichbarer Zieltemperatur wird die gesamte Aufheizung erfasst, unabhängig von der gewählten Position der Animation.

Real- und Imaginärteile werden analytisch mit skalierten hyperbolischen Ausdrücken ausgewertet, um Überläufe und wiederholte komplexe Rechnungen zu vermeiden. Eine direkte komplexe Impedanzberechnung bleibt als rechnerische Gegenprüfung erhalten.

• Eine einzelne Zelle mit Graphitanode, symmetrischer Sinusstrom, kein Gleichstromanteil, nahezu konstanter SOC, gleichmäßige Zelltemperatur und konstante spezifische Wärmekapazität. Die zugrunde liegende Zelle ist eine Graphit/NMC-Zelle bei 50 % SOC.
• Liu et al., Table 1 (p. 7): Masse 56 g und spezifische Wärmekapazität 1,24 J/(g·K). Anoden-OCP: 0,1354 V (p. 5). Aktivierungsenergien des Ladungsdurchtritts: 76,8 und 81,9 kJ/mol (p. 9).
• Bei −20 °C näherungsweise aus Liu et al., Fig. 4 (p. 10), abgelesen: Anode $R_{ct}=0.58$ Ω, CPE-Skalierungsparameter 0,72 F und Exponent 0,73; Kathode 0,14 Ω, 4,5 F und 0,85; ohmscher Zellwiderstand 0,012 Ω und Deckschichtwiderstand 0,0075 Ω, CPE-Skalierungsparameter der Deckschicht 0,2 F und Exponent 0,87. Die Werte sind gerundete Diagrammablesungen und stammen nicht aus einer veröffentlichten Parametertabelle. CPE-Skalierungsparameter und -Exponenten werden als temperaturunabhängig angenommen.
• Messtechnisch zu bestätigende Beispielannahmen: Diffusionswiderstandsparameter 0,10/0,05 Ω und Diffusionszeiten 100/100 s; Aktivierungsenergie der Diffusion 30 kJ/mol; Aktivierungsenergien des ohmschen Widerstands und der Deckschicht 10/20 kJ/mol; Induktivität 1 µH; Wärmeleitwert zur Umgebung 0,05 W/K; Strom- und Spannungsgrenzen. Damit liegt keine Charakterisierung einer kommerziellen Zelle vor.
• Die reversible entropische Wärme hat den Mittelwert null und hebt sich über vollständige Perioden auf. Umrichterverluste und die einer separaten Stromversorgung entnommene Energie werden nicht berücksichtigt; die erzeugte Wärme entspricht der elektrischen Wirkverlustenergie der Zelle.

Diese Berechnung bildet eine modellbasierte Anodenpotenzialbegrenzung ab und belegt nicht, dass in einer nicht charakterisierten Zelle kein Lithium-Plating auftritt. Elektrodenpotenzial, Impedanz, Reserve und Spannungsreserve an den Klemmen müssen für die vorgesehene Zellchemie, den SOC, den Alterungszustand und die Temperatur gemessen werden. Bei großer Anregungsamplitude kann die lineare EIS ihre Gültigkeit verlieren; räumliche Gradienten im Elektrolyten und in den Partikeln werden nicht aufgelöst. Beginnen Sie mit einer Rampe, halten Sie den Gleichstromanteil bei null und validieren Sie das Profil unter Überwachung von Temperatur und Spannung auf Zellebene.

Das veröffentlichte EIS-Beispiel umfasst −20 bis 25 °C und 0,01 Hz bis 10 kHz. Der Betrieb bei −30 °C beruht auf einer Arrhenius-Extrapolation. Die Erweiterung um planare Diffusion mit reflektierender Randbedingung und die angenommenen Diffusionsparameter bilden das angepasste Modell von Liu nicht nach; die experimentellen Aufheizzeiten der Veröffentlichung werden nicht als Validierung herangezogen. Die beiden Diffusionsreferenztests prüfen unabhängige veröffentlichte Zahlenwerte; die Wärmekapazitätsprüfung kontrolliert Quelldaten und Einheiten.

Die Frequenzsuche verwendet 501 logarithmisch verteilte Frequenzen und 151 Temperaturstützstellen. Anschließend wird das beste Festfrequenzintervall mit adaptiver Zeitquadratur verfeinert. Der adaptive Frequenzkandidat wird zwischen den berechneten Stützstellen interpoliert. Im voreingestellten Band über fünf Dekaden verursacht der Rasterabstand nahe einem Optimum eine Frequenzunsicherheit von bis zu etwa 1,2 %; breitere Bänder vergröbern das Raster. Die Heizleistung wird bei jeder Quadratur- und DGL-Auswertung aus der Impedanz neu berechnet und niemals auf ihren Anfangswert festgesetzt. Außerhalb des Temperaturbereichs des Steuerprofils wird die Frequenz auf dem jeweiligen Endpunktwert gehalten; bei adaptivem Betrieb wird weiterhin mit der Festfrequenz als Rückfalloption verglichen.

Die Zeitintegration halbiert Intervalle adaptiv und vergleicht das logarithmische Integral eines Abschnitts mit dem seiner beiden Hälften. Zulässig sind höchstens 24 Verfeinerungsdurchläufe und 4097 gespeicherte Temperaturstützstellen. Die Berechnung warnt, wenn die Fehlerschätzung oder die strengere Prüfung der gewählten Strategie nicht hinreichend aufgelöst bleibt. Nahe einem Aufheizstillstand reagiert das Ergebnis zudem sehr empfindlich auf Zellparameter und Wärmeverluste. Frequenzen, bei denen das Ziel unerreichbar ist, erscheinen bewusst als Lücken. Erreicht keine Festfrequenz die Zieltemperatur, bleibt das Optimum leer; die Animation verwendet dann die Rasterfrequenz mit der höchsten anfänglichen Heizleistung.

Bei erreichbarer Zieltemperatur zeigen Diagramme und Animation die vollständige Aufheizung, auch bei einer Dauer über 200000 s. Erzeugte und abgegebene Wärme umfassen dieses gesamte Intervall. Fälle mit Aufheizstillstand zeigen einen einstündigen Verlauf; ihre Energieergebnisse beziehen sich auf diese Stunde. Die eingegebene verstrichene Zeit wählt lediglich eine Position innerhalb der zwischengespeicherten Lösung und verkürzt niemals die Energieintegration. Die Animation interpoliert die zwischengespeicherten thermischen und elektrischen Verläufe. Sie veranschaulicht die langsame Aufheizung und zeigt keinen Wellenverlauf einzelner Wechselstromperioden. Ein Modell mit räumlich einheitlicher Temperatur kann lokale Überhitzungen im Zellinneren übersehen, insbesondere bei dicken Zellen oder schneller Aufheizung. Temperaturgradienten im Batteriepack, Verbindungen, Balancing, Schaltharmonische, nichtlineares Elektrodenverhalten und die Kinetik des Lithium-Platings erfordern weiterführende Modelle.