Hi Adrian,
It's a little bit of a hack, but I got something to work, see below.
Unfortunately, I think you need to explicitly delete the fixed value
constraint from the CellVariable's faceConstraints list, which
requires some knowledge of FiPy's internals and remember the order in
which the
Hi All,
I have a problem where I would want to specify fixed value boundary
conditions for times t=t0. The way I've implemented this is:
a.constrain(a_0,(mesh.exteriorFaces))
while t = t0:
a.faceGrad.constrain(0.,mesh.exteriorFaces)
eq.solve(dt=timeStepDuration)
t+=dt
