Lithium Plating Submodel if statement #2625
-
|
Hello, I am working on a submodel for lithium plating from the O'Kane 2022 submodel and I wanted to use a Butler-Volmer equation that turns on and off depending on the overpotential like this:
But when I try to use an if..then statement in the submodel with eta_stripping I get the error: File ~\env\lib\site-packages\pybamm\expression_tree\symbol.py:660, in Symbol.bool(self) NotImplementedError: Boolean operator not defined for Symbols. What can I do to get around this/what am I doing wrong? |
Beta Was this translation helpful? Give feedback.
Replies: 1 comment 2 replies
-
|
Instead of the if statement, try multiplying each of the statement with |
Beta Was this translation helpful? Give feedback.
-
|
Hello, I am Hanyu and I also have a question related to this issue: Are there any specific reasons why the conditions of lithium plating and SEI growth (eta_stripping should be less than zero ) are not added to the current model? Like issues with the solver? Also, the following error occurs when adding the following code to the lithium plating function: But it seems like the bug only occurs at certain cycles where the value of the function changes from zero to one. psetup failed: .../casadi/interfaces/sundials/idas_interface.cpp:849: Calculating Jacobian failed |
Beta Was this translation helpful? Give feedback.
-
|
It's hard to tell where the problem is just with the log. Can you share a minimum working example? SEI growth and lithium plating can be enabled through model options, you can find an example here: https://docs.pybamm.org/en/latest/source/examples/notebooks/models/coupled-degradation.html |
Beta Was this translation helpful? Give feedback.
Instead of the if statement, try multiplying each of the statement with
(eta_stripping >= 0)and(eta_stripping < 0)which should define a heaviside function that does what you want.