Hi All,
I have a problem where I would want to specify fixed value boundary
conditions for times t<t0 and then switch to fixed gradient boundary
conditions for t>=t0. The way I've implemented this is:
a.constrain(a_0,(mesh.exteriorFaces))
while t<tottime:
if t > = t0:
a.faceGrad.constrain(0.,mesh.exteriorFaces)
eq.solve(dt=timeStepDuration)
t+=dt
But what it seems to be happening is that the fixed value bc is not
being overwritten by the fixed gradient one. How can I do to remove the
fixed value bc? Is there a way to check the boundary conditions on a
variable?
Thanks,
Adrian.
_______________________________________________
fipy mailing list
[email protected]
http://www.ctcms.nist.gov/fipy
[ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
