hello,
can someone help me with figuring out the unreasonable capacity loss of long life modeling.
I am doing long cycle life modeling on my own cell with relatively large capacity. The parameters were based on Mohtat2020, with some modifications in accordance with measuring and fitting results.
The code is like this:
model = pybamm.lithium_ion.DFN({
"SEI": "reaction limited",
"SEI porosity change": "true",
"lithium plating": "partially reversible",
"lithium plating porosity change": "true",
"particle mechanics": "swelling and cracking",
"SEI on cracks": "true",
"loss of active material": "stress-driven",
})
var_pts = {
"x_n": 5, # negative electrode
"x_s": 5, # separator
"x_p": 5, # positive electrode
"r_n": 10, # negative particle
"r_p": 10, # positive particle
}
sim = pybamm.Simulation(model, var_pts=var_pts)
param = pybamm.ParameterValues("Mohtat2020")
experiment = pybamm.Experiment(
[
"Discharge at 1 C until 2.8V (10 minute period)",
"Rest for 0.5 hour",
"Charge at 2 C until 4.4 V (10 minute period)",
"Rest for 0.5 hour"
] * 2000
)
The discharge capacity simulation result matches very well with the experiment value after several rounds of parameters fitting. In details, the discharge_capacity decreased from 115Ah to nearly 104Ah. However, the LLI loss is unseasonably low, less than 0.2Ah, as following picture.
Why is that?
And during simulation, I faced several kinds of errors as following:
1)t = 528.534 and h = 8.75955e-15, the corrector convergence failed repeatedly or with |h| = hmin.
2)538.078, , mxstep steps taken before reaching tout.
3)t = 75.9756, repeated recoverable residual errors.
4)[ERROR] solution.check_ys_are_not_too_large(397): Solution for ‘y[1:2]’ exceeds the maximum allowed value of `100000.0. This could be due to incorrect scaling, model formulation, or parameter values. The maximum allowed value is set by 'pybammm.settings.max_y_value
By the way, my pybamm version is 25.1.1