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"]);