我正在使用fipy来模拟线性化的
Poisson-Boltzmann equation,这本质上是
我假设我可以将f(x)建模为边界条件.如果epsilon(x)是常量,fipy可以处理这个:
phi = CellVariable(mesh)
dielectric_solvent = 80.0
dielectric_inner = 4.0
LHS = (DiffusionTerm(coeff = dielectric_solvent))
RHS = phi
eq = LHS == RHS
dr = np.linalg.norm(mesh.faceCenters, axis=0)
mask = (dr<.5 mesh.exteriorfaces>
phi.constrain(1, mask)
mask = (dr>.5) * mesh.exteriorFaces
phi.constrain(0, mask)
sol = eq.solve(var=phi)
赠送:
我想做的是让epsilon(x)在空间中作为函数变化,但DiffusionTerm只能取常数.如何实现空间变化的电介质项?