Hi, I'm new to FiPy and first just want to say thanks for all the time and effort put into this, it looks like a great learning and research tool!
I also have a comment and question about the simple convection-with-source example at: http://www.ctcms.nist.gov/fipy/examples/convection/generated/examples.convection.source.html The (very minor!) comment is that the equation to be solved is stated in the text of the example as: d(phi)/dx - alpha*phi = 0 but the equation solved is: d(phi)/dx + alpha*phi = 0 which gives the negative gradient of the solution (also in the code the equation is specified by eq = PowerLawConvectionTerm((1,)) + ImplicitSourceTerm(alpha)). My question is about the part of the description reading "The boundary condition at x = L will require the implementation of an outflow boundary condition, which is not currently implemented in FiPy. An ImplicitSourceTerm object will be used to represent this term." I'm not expert in well-posedness and requirements of numerical formulations, but since there is a boundary condition given at x=0 and the equation is first order, why is a boundary condition at x=L required? I tried taking the condition at x=L out and the code still works, but perhaps there is a default Neumann condition implemented anyway? Thanks for your time, Oliver. -- Oliver Maclaren BE(Hons) PhD Candidate Auckland Bioengineering Institute The University of Auckland New Zealand _______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
