On Fri, Jul 1, 2016 at 2:44 PM, Gopalakrishnan, Krishnakumar <k.gopalakrishna...@imperial.ac.uk> wrote: > > Now, my question is, for the explicit term in the BC, can't I just define in > a regular Neumann BC style of definition ? , i.e. > phi.facegrad.constrain([alpha/ (1 - beta*dx/2)], mesh.facesRight]) . i.e. > will this ordinary Neumann correctly add up with the Implicit Boundary > condition term , implemented using the ImplicitSourceTerm method ? > > Or, should I treat this explicit BC term as part of the PDE itself, (like > the trick done for the implicit BC), ie. Do we need something along the lines > of fp.SourceTerm((mesh.faceNormals * explicitCoeff* > mesh.facesRight).divergence)) ? In either case, the diffusion > coefficient along the problematic boundary has been forced to be zero.
I recommend trying both methods and implementing some 1D test cases to double check. I always do this no matter how much I trust the code. -- Daniel Wheeler _______________________________________________ fipy mailing list fipy@nist.gov http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]