Source code for mpet.electrode.materials.LiC6

import numpy as np




[docs]
def LiC6(self, y, ybar, muR_ref):
    """ Ferguson and Bazant 2014 """
    muRtheta = -self.eokT*0.12
    muR1homog, muR2homog = self.graphite_2param_homog(
        y, self.get_trode_param("Omega_a"), self.get_trode_param("Omega_b"),
        self.get_trode_param("Omega_c"), self.get_trode_param("EvdW"))
    muR1nonHomog, muR2nonHomog = self.general_non_homog(y, ybar)
    muR1 = muR1homog + muR1nonHomog
    muR2 = muR2homog + muR2nonHomog
    actR1 = np.exp(muR1/self.T)
    actR2 = np.exp(muR2/self.T)
    muR1 += muRtheta + muR_ref
    muR2 += muRtheta + muR_ref
    return (muR1, muR2), (actR1, actR2)