I recently noticed that the faceValue of a variable is incorrect when a Neumann (constant gradient) boundary condition is applied. For example, let's assume the following:
dx = 1 nx = 5 dTr/dx = 4 for right face Tl = 300 for left face x = [0.5, 1.5, 2.5, 3.5, 4.5] The steady solution to this problem would be: T = [302, 306, 310, 314, 318] T.faceValue = [300, 304, 308, 312, 316, 320] However, in FiPy I get the following for the faceValues T.faceValue = [300, 304, 308, 312, 316, 318] I'm assuming this is a bug in the faceValue code. Would this be correct? - - - - - - - - - - - - - - - - - L.Bryce Whitson Jr. [email protected]
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
