pythonfinite-element-analysisfipy

Variable dependent flux boundary condition in fipy

I have a case where reaction flux is present at a boundary and flux is proportional to the concentration variable c. Will the following correctly apply the boundary condition:

c.faceGrad.constrain(k*c.faceValue, mesh.physicalFaces["ReactingFace"]);
Solution

  • It sounds like you have a Robin Condition. Please revise your question if that documentation doesn't help.