Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.

Working of Non-Faradaic Reactions Node under Porous Electrode Node in Lithium-Ion Battery Interface

Please login with a confirmed email address before reporting spam

Hi,

I am currently working on an End-of-Life (EOL) battery model and trying to model transition metal dissolution from the electrodes as one of the degradation mechanisms. The reaction association with this mechanism is given in the attached "Mn_Dissolution.png" discussed in the paper (Lee, Yoonkoo et al. “A Comprehensive Experimental and Modeling Study on Dissolution in Li-Ion Batteries.” Journal of The Electrochemical Society (2019): n. pag.).

The rate of Mn dissolution is modeled with the help of above reaction as,

where is the rate constant and is the concentration of hydrogen ions obtained by solving the species balance given in the atttached "Species_Balance.png".

As seen in the "Mn_Dissolution.png", Li+ ions are generated during this reaction which needs to be coupled back to the Li+ species balance which is solved internally in the Lithium-Ion Battery Interface. So, to do this, I used the Non-Faradaic Reactions Node under Porous Electrode Node where we can add this source as , which is coupled back to the Li+ species balance as given in the COMSOL equations "COMSOL_Eqs.png" of the Porous Electrode Node.

Some concerns arose while implementing the above mechanism:

  1. Given that the term in "COMSOL_Eqs.png" is with a '+' sign, signifying that, a positive value (R_Mndiss * a_pos * 2/y_Mn * t_factor) entered in the tab in "NonFaradaic_Rxn.png" should result in the increasing concentration of Li+, however a decreasing trend was observed. In another run, a negative value (-R_Mndiss * a_pos * 2/y_Mn * t_factor) was entered which gave an increasing trend for Li+ concentration. What could be the possible reason behind this behaviour?
  2. Secondly, what is the role of "Reaction rate (mol/(m^3s))" value under the "Reaction rate for dissolving-depositing species" (See "NonFaradaic_Rxn.png"). Should it be kept zero or equal to the magnitude of term. One of the COMSOL examples employing the Non-Faradaic Reactions node Adsorption-Desorption Voltammetry used the "Reaction rate (mol/(m^2s))" value as equal to entered for the source term (See "NonFaradaic_Rxn2.png"). Please clarify regarding the term's magnitude as well as sign.


0 Replies Last Post Jan 2, 2024, 5:44 a.m. EST
COMSOL Moderator

Hello Shantanu Sontakke

Your Discussion has gone 30 days without a reply. If you still need help with COMSOL and have an on-subscription license, please visit our Support Center for help.

If you do not hold an on-subscription license, you may find an answer in another Discussion or in the Knowledge Base.

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.