Skip to content

Latest commit

 

History

History
134 lines (110 loc) · 4 KB

README.md

File metadata and controls

134 lines (110 loc) · 4 KB

Lithium-ion Battery Realisation Algorithms (LiiBRA)

Build Status SciML Code Style ColPrac: Contributor's Guide on Collaborative Practices for Community Packages DOI:10.1016/j.est.2022.105637

Create and Simulate Reduced Order Lithium-Ion Battery Models

LiBRA provides an open-source implementation of realisation algorithms used for generating reduced-order state-space models. This work aims to develop real-time capable physics-informed models for deployment onto embedded hardware. LiiBRA provides capabilities in offline and online model creation, as well as a framework for lithium-ion degradation predictions. For more information on LiiBRA, and the computationally-informed discrete realisation algorithm (CI-DRA), please refer to the publication above.

Further examples are located in the "examples" directory. Please open an issue if you have requests or ideas for additional examples.

Install (Julia 1.7 and later)

(v1.7) pkg> add LiiBRA

(Type ] to enter package mode.)

Example Usage

using LiiBRA, Plots

Setup:

Sₑ = 4
Sₛ = 2
Cell = Construct("LG M50")
Spatial!(Cell, Sₑ, Sₛ)
Ŝ = collect(1.0:-0.25:0.0)
SOC = 0.75
Cell.Const.T = 298.15

Realisation:

A,B,C,D = Realise(Cell,Ŝ);

HPPC Simulation:

Results = HPPC(Cell,Ŝ,SOC,4.0,-3.0,A,B,C,D);

Plotting Results:

plotly()
plot(Results.t, Results.Cell_V;
     legend=:topright,
     color=:blue,
     bottom_margin=5Plots.mm,
     left_margin = 5Plots.mm,
     right_margin = 15Plots.mm,
     ylabel = "Terminal Voltage (V)",
     xlabel = "Time (s)",
     title="HPPC Voltage",
     label="Voltage",
     size=(1280,720)
    )

plot(Results.t, Results.Ce;
     legend=:topright,
     bottom_margin=5Plots.mm, 
     left_margin = 5Plots.mm, 
     right_margin = 15Plots.mm, 
     ylabel = "Electrolyte Concen. (mol/m³)", 
     xlabel = "Time (s)",
     title="Electrolyte Concentration",
     label=["Neg. Separator Interface" "Neg. Current Collector" "Pos. Current Collector" "Pos. Separator Interface"], 
     size=(1280,720)
    )

plot(Results.t, Results.Cseₚ;
     legend=:topright,
     bottom_margin=5Plots.mm, 
     left_margin = 5Plots.mm, 
     right_margin = 15Plots.mm, 
     ylabel = "Concentration (mol/m³)", 
     xlabel = "Time (s)",
     title="Positive Electrode Concentration",
     label=["Current Collector" "Separator Interface"], 
     size=(1280,720)
    )

plot(Results.t, Results.Cseₙ;
     legend=:topright,
     bottom_margin=5Plots.mm, 
     left_margin = 5Plots.mm, 
     right_margin = 15Plots.mm, 
     ylabel = "Concentration (mol/m³)", 
     xlabel = "Time [s]", 
     title="Negative Electrode Concentration",
     label=["Current Collector" "Separator Interface"],
     size=(1280,720)
    )

Bug Tracking

Please report any issues using the Github issue tracker. All feedback is welcome.