Source code for mpet.electrode.materials.LiFePO4

import numpy as np



[docs]
def LiFePO4(self, y, ybar, muR_ref):
    """ Bai, Cogswell, Bazant 2011 """
    muRtheta = -self.eokT*3.422
    muRhomog = self.reg_sln(y, self.get_trode_param("Omega_a"))
    muRnonHomog = self.general_non_homog(y, ybar)
    muR = muRhomog + muRnonHomog
    actR = np.exp(muR/self.T)
    muR += muRtheta + muR_ref
    return muR, actR