ZeroDivisionError in running code in Windows

[Akila1993]

Hi All,

I am trying to run the following peace of code and gives me the ZeroDivisionError. Can anyone comment about it. This code was perfectly working with MacOS.

%pip install "pybamm[plot,cite]" -q   
import pybamm
import matplotlib.pyplot as plt


model = pybamm.lithium_ion.DFN(
    {   
        "thermal": "lumped",
        "surface form": "differential", 
        "SEI": "solvent-diffusion limited",
        "SEI porosity change": "true",
        "lithium plating": "partially reversible",
        "lithium plating porosity change": "true",  
        "particle mechanics":"swelling only",
        "SEI on cracks": "true",
        "loss of active material": "stress-driven",
        
    }
)

---------------------------------------------------------------------------
ZeroDivisionError                         Traceback (most recent call last)
Cell In[3], line 1
----> 1 model = pybamm.lithium_ion.DFN(
      2     {   
      3         "thermal": "lumped",
      4         "surface form": "differential", 
      5         "SEI": "solvent-diffusion limited",
      6         "SEI porosity change": "true",
      7         "lithium plating": "partially reversible",
      8         "lithium plating porosity change": "true",  
      9         "particle mechanics":"swelling only",
     10         "SEI on cracks": "true",
     11         "loss of active material": "stress-driven",
     12         
     13     }
     14 )

File ~\env\lib\site-packages\pybamm\models\full_battery_models\lithium_ion\dfn.py:32, in DFN.__init__(self, options, name, build)
     29 self.x_average = False
     30 super().__init__(options, name)
---> 32 self.set_submodels(build)
     34 pybamm.citations.register("Doyle1993")

File ~\env\lib\site-packages\pybamm\models\full_battery_models\lithium_ion\base_lithium_ion_model.py:59, in BaseModel.set_submodels(self, build)
     56 self.set_total_interface_submodel()
     58 if build:
---> 59     self.build_model()

File ~\env\lib\site-packages\pybamm\models\full_battery_models\base_battery_model.py:907, in BaseBatteryModel.build_model(self)
    905 def build_model(self):
    906     # Build model variables and equations
--> 907     self._build_model()
    909     # Set battery specific variables
    910     pybamm.logger.debug("Setting voltage variables ({})".format(self.name))

File ~\env\lib\site-packages\pybamm\models\base_model.py:552, in BaseModel._build_model(self)
    548     self.build_fundamental()
    550 self.build_coupled_variables()
--> 552 self.build_model_equations()

File ~\env\lib\site-packages\pybamm\models\full_battery_models\base_battery_model.py:877, in BaseBatteryModel.build_model_equations(self)
    872 for submodel_name, submodel in self.submodels.items():
    873     pybamm.logger.verbose(
    874         "Setting rhs for {} submodel ({})".format(submodel_name, self.name)
    875     )
--> 877     submodel.set_rhs(self.variables)
    878     pybamm.logger.verbose(
    879         "Setting algebraic for {} submodel ({})".format(
    880             submodel_name, self.name
    881         )
    882     )
    884     submodel.set_algebraic(self.variables)

File ~\env\lib\site-packages\pybamm\models\submodels\interface\sei\sei_growth.py:225, in SEIGrowth.set_rhs(self, variables)
    222         l_cr = variables["Negative particle crack length [m]"]
    223         dl_cr = variables["Negative particle cracking rate [m.s-1]"]
    224     spreading_outer = (
--> 225         dl_cr / l_cr * (self.phase_param.L_outer_crack_0 - L_outer)
    226     )
    227     spreading_inner = (
    228         dl_cr / l_cr * (self.phase_param.L_inner_crack_0 - L_inner)
    229     )
    230 else:

File ~\env\lib\site-packages\pybamm\expression_tree\symbol.py:575, in Symbol.__truediv__(self, other)
    573 def __truediv__(self, other):
    574     """return a :class:`Division` object."""
--> 575     return pybamm.simplified_division(self, other)

File ~\env\lib\site-packages\pybamm\expression_tree\binary_operators.py:1067, in simplified_division(left, right)
   1064     return (1 / right) * left
   1066 # Check for Concatenations and Broadcasts
-> 1067 out = _simplified_binary_broadcast_concatenation(left, right, simplified_division)
   1068 if out is not None:
   1069     return out

File ~\env\lib\site-packages\pybamm\expression_tree\binary_operators.py:682, in _simplified_binary_broadcast_concatenation(left, right, operator)
    680     return left._unary_new_copy(operator(left.orphans[0], right))
    681 elif isinstance(right, pybamm.Broadcast) and left.domain == []:
--> 682     return right._unary_new_copy(operator(left, right.orphans[0]))
    684 # Concatenation commutes with elementwise operators
    685 # If one of the sides is constant then commute concatenation with the operator
    686 # Don't do this for ConcatenationVariable objects as these will
    687 # be simplified differently later on
    688 if isinstance(left, pybamm.Concatenation) and not isinstance(
    689     left, pybamm.ConcatenationVariable
    690 ):

File ~\env\lib\site-packages\pybamm\expression_tree\binary_operators.py:1059, in simplified_division(left, right)
   1057 # anything divided by zero raises error
   1058 if pybamm.is_scalar_zero(right):
-> 1059     raise ZeroDivisionError
   1061 # Move constant to always be on the left
   1062 # For a division, this means (var / constant) becomes (1/constant * var)
   1063 if right.is_constant() and not left.is_constant():

ZeroDivisionError: 

[brosaplanella]

Which PyBaMM version are you using, and which one were you using in Mac?

[Akila1993]

In Mac: : Version: 22.10
In windows : Version: 23.5

Please also let me know the commands for both upgrade and degrade of the version for MacOS and windows please.

I think since I need to also work with PyBaMM-EIS the version should be compatible with that which already is fine with In Mac: : Version: 22.10.

[brosaplanella]

I think it is complaining because the initial crack length is zero, so that might be the cause.