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
