def evalNa(self): listNa = [] for xp in self.coords: a = xp / (2.0 * sqrt(self.De * self.time)) Na = self.C0 * erfc(a) listNa.append(Na) return listNa
def evalNa(self): listNa=[] for xp in self.coords: a = xp / (2.0 * sqrt(self.De * self.time)) Na = self.C0 * erfc(a) listNa.append(Na) return listNa
def evalPH(self): listPh=[] for xp in self.coords: a = xp / (2.0 * sqrt(self.De * self.time)) Na = self.C0 * erfc(a) d = (Na*Na)+(4.0*(self.K_Na+(self.K_Silice * self.K_Silicique))) d = pow(d,0.5) pH = -1.*(log10(0.5*(d-Na))) #pH = -1.*(log(0.5*(d-Na)))/(log(10.)) listPh.append(pH) return listPh
def evalPH(self): listPh = [] for xp in self.coords: a = xp / (2.0 * sqrt(self.De * self.time)) Na = self.C0 * erfc(a) d = (Na * Na) + (4.0 * (self.K_Na + (self.K_Silice * self.K_Silicique))) d = pow(d, 0.5) pH = -1. * (log10(0.5 * (d - Na))) #pH = -1.*(log(0.5*(d-Na)))/(log(10.)) listPh.append(pH) return listPh
def chargeFunction(listOfPoints, time, h0, k0, s0, d0): h = [] k = k0 / s0 for point in listOfPoints: h.append(H0 * erfc(point / (2. * sqrt(k0 * time / (s0 * d0))))) return h
def chargeFunction(listOfPoints,time,h0,k0,s0): h = [] k = k0/s0 for point in listOfPoints: h.append( H0*erfc ( point/( 2.*sqrt(k0*time/s0) ) ) ) return h