Schnellladen ohne Lithium-Plating: sichere C-Rate mit SPMe

Das vollständige kanonische SPMe, Marquis et al. (2019), dimensionsbehaftete Gl. (48), beigefügte S. 32–33, wird mit den angepassten LG-M50-Parametern aus Chen et al. (2020), Tabelle IX, beigefügte S. 20, gelöst, wie sie in PyBaMM 25.10.2 als Chen2020 implementiert sind. Die gemessenen Diffusionskoeffizienten und Konzentrationen aus Tabelle VII werden nicht anstelle dieses angepassten Parametersatzes verwendet.

Für jedes Partikel gilt $partial_t c_s=r^{-2}partial_r(r^2D_spartial_r c_s)$, mit verschwindendem Fluss im Zentrum und vorgegebenem faradayschem Fluss an der Oberfläche. Über die drei Zellbereiche gilt $epsilonpartial_t c_e=partial_x(epsilon^{1.5}D_{e,0}partial_xc_e)+(1-t^+)a j/F$, mit stetiger Konzentration und stetigem Fluss an den Übergängen sowie undurchlässigen Randbedingungen an den Stromableitern. Ladeströme sind in der Benutzeroberfläche positiv; der im Artikel bei Entladung positive Strom hat das entgegengesetzte Vorzeichen.

Kugelschalenvolumina und harmonisch gemittelte Diffusionswiderstände an den Grenzflächen bilden konservative Finite-Volumen-Operatoren: 28 radiale Zellen je Partikel und 18/6/18 Zellen in Dickenrichtung. Eine massengewichtete symmetrische Eigenwertzerlegung berücksichtigt sämtliche diskreten Diffusionsmoden. Jedes Konstantstromsegment wird mit dem exakten Matrixexponential fortgeschrieben, einschließlich der Nullmode für die Erhaltungsgröße. Damit werden die vollständigen räumlich diskretisierten instationären Gleichungen gelöst, ohne Näherungsansatz für das Partikelprofil und ohne Ersatz durch einen stationären Elektrolytzustand.

Die Austauschstromdichte ist $j_0=m(T)sqrt{c_s(c_{s,max}-c_s)}langlesqrt{c_e}rangle$; die symmetrische Butler–Volmer-Gleichung wird analytisch zu $eta=2RT/F,mathrm{asinh}(j/(2j_0))$ umgestellt. Chens Aktivierungsenergien des Ladungstransfers von 35,0 und 17,8 kJ/mol liefern $m(T)=m_{ref}exp[(E_a/R)(1/T_{ref}-1/T)]$ (beigefügte S. 17). Der kanonische Elektrolytdiffusionskoeffizient und die kanonische Elektrolytleitfähigkeit werden gemäß Gl. (48) bei 1000 mol/m³ ausgewertet.

Die über die Elektroden gemittelte Spannung enthält Beiträge aus Reaktion, Elektrolytkonzentration sowie ohmschen Verlusten im Elektrolyten und in der Festphase. Die lokalen Festphasen- und Elektrolytpotentiale der Anode werden im Zentrum jeder Anodenzelle und an beiden Rändern rekonstruiert; das Minimum von $phi_s-phi_e$ wird mit 0 V gegen Li/Li⁺ und der eingegebenen Potentialreserve verglichen, nicht das gemittelte Anodenpotential oder das Potential des willkürlich auf Masse gelegten Stromableiters.

Eine Bisektion mit 12 Iterationen prüft für jedes Kandidatensegment 121 ungleichmäßig verteilte Zeitpunkte. Die abschließenden CC-Kurven verwenden 221 Abtastpunkte, jedes der vier SOC-Intervalle der Stufenladung jeweils 161. Die Suche für die Stufenladung übernimmt beide Partikelkonzentrationsfelder und das Elektrolytkonzentrationsfeld zwischen den Intervallen und verhindert Stromerhöhungen.

• Eine neue, anfangs relaxierte LG-M50-Zelle mit homogenen Partikelkonzentrationen und einer Elektrolytkonzentration von 1000 mol/m³; konstante, räumlich homogene Zelltemperatur.
• Der nominelle SOC wird ausgehend von den angepassten Konzentrationen bei Vollladung (29866 und 17038 mol/m³) auf eine Kapazität von 5 Ah bezogen. Die daraus resultierenden Stöchiometrien am entladenen Ende folgen aus der Lithiumerhaltung auf dieser Kapazitätsskala; sie entsprechen nicht den gemessenen Endwerten aus Tabelle VII.
• Konstante angepasste Feststoffdiffusionskoeffizienten, feste Partikelradien und Porositäten, ideale Elektrolytaktivität, Überführungszahl 0,2594 und Bruggeman-Exponent 1,5. Das Graphit–SiOₓ-Gemisch der Anode wird wie in diesem Parametersatz durch ein effektives Partikel dargestellt.
• Die Temperatur verändert die Ladungstransferkinetik und die Temperaturspannung. Feststoffdiffusion und Elektrolyttransport behalten die Referenzwerte bei 25 °C, da dieser Quellensatz keine Aktivierungsenergien für deren Temperaturabhängigkeit liefert.
• Keine Nebenreaktionen, kein SEI-Schichtwiderstand, keine Hysterese, keine thermische Rückkopplung, keine mechanischen Spannungen, kein Lithiumverlust und kein Verlust an Aktivmaterial.

„Zulässig“ bedeutet, dass der aufgelöste SPMe-Verlauf oberhalb der geforderten lokalen Potentialreserve der Anode und unterhalb der eingegebenen Klemmenspannungsgrenze bleibt sowie einen Lithiierungsgrad der Partikel von 0,0001–0,9999 und eine Elektrolytkonzentration von 50–2500 mol/m³ einhält. Diese Konzentrationsgrenzen sind hier gewählte rechnerische Betriebsgrenzen, keine experimentell nachgewiesenen Gültigkeitsgrenzen. Die Laderatensuche reicht von 0,005C bis zur eingegebenen Suchobergrenze; die Breite des Suchintervalls beträgt höchstens (Suchobergrenze − 0,005)/4096 C. Kein Ergebnis beansprucht ein Maximum oberhalb dieser Suchobergrenze.

Das kanonische SPMe ist ein asymptotisches Modell poröser Elektroden. Große Abweichungen der Elektrolytkonzentration von 1000 mol/m³, steile Bereiche des Ruhepotentials (OCP), hohe Ströme oder stark ungleichmäßige Reaktionen verringern seine Genauigkeit. Das linearisierte Elektrolytpotential und die Transportkoeffizienten bei Referenzbedingungen gehören zum zitierten kanonischen Modell. Die Temperaturextrapolation unterhalb von Chens kinetischen Messungen bei 25–60 °C ist besonders unsicher; die vernachlässigte Temperaturabhängigkeit des Transports kann die Leistungsfähigkeit kalter Zellen erheblich überschätzen. Eine positive Potentialreserve berücksichtigt weder Parameterunsicherheit noch lokale Defekte oder die Keimbildung beim Plating.

Das endliche Rechengitter und die begrenzte Anzahl zeitlicher Prüfpunkte können dünne anfängliche Diffusionsschichten oder kurze Potentialeinbrüche übersehen. Sämtliche Finite-Volumen-Moden werden berücksichtigt; begrenzt sind lediglich die räumliche Auflösung und die zeitliche Abtastung. Der Lithiierungsgrad an der Partikeloberfläche wird quadratisch extrapoliert. Die radiale Darstellung verwendet PCHIP-Interpolation durch die berechneten radialen Zellen, die Elektrolytdarstellung PCHIP-Interpolation zwischen den 42 berechneten Punkten in Dickenrichtung; die Animation wertet den gespeicherten Segmentzustand zum gewählten Zeitpunkt aus, statt die Laderatensuche erneut abzuspielen.

Unabhängige PyBaMM-Referenzrechnungen mit 60 Zellen je Partikel und 60/20/60 Zellen im Elektrolyten konfigurieren den kanonischen Transport nach Gl. (48) und die Kinetik mit gemittelter Austauschstromdichte ausdrücklich; es sind keine unveränderten Läufe des PyBaMM-Composite-SPMe. Diese Prüfungen verifizieren die Implementierung, nicht die Genauigkeit des Schnellladens realer Zellen. Das Profil mit vier SOC-Intervallen ist stufenweise zulässig und muss nicht global zeitoptimal sein. Prüfen Sie die tatsächliche Zelltemperatur, Alterung, Kühlung und die Ladegrenzen des Herstellers gesondert. Der Testverlauf wird bewusst über die Spannungsgrenze hinaus fortgesetzt, um das vorgeschlagene CC-Fenster zu untersuchen; verwerfen Sie seine Vorhersagen, sobald eine Konzentrationsgrenze verletzt wird.