On Fri, Jun 22, 2012 at 4:29 PM, Allawala, Altan <[email protected]>wrote:
> Hi Daniel, > > Eureka. That certainly fixed it. In fact, I applied your changes to the > other approach too (ie: my original code where I expressed the triple > derivative by placing a convection term inside a diffusion term) and it's > now giving the same correct answer as well. It is interesting that changing > the Upwind to the CentralDifference convection term didn't make a > difference in this case (probably because the second equation in this setup > had only a convection term, no diffusion) but the actual change occurred > when the left constraint was abolished in favor of a dirac-delta function > approach. The problem in FiPy is being able to easily write the boundary conditions. We've gone round in circles with this. For me, writing the boundary condition as a source allows me to understand what I'm doing a little better. Anyway, good luck with your work. > The FiPy works in mysterious ways! > :-) -- Daniel Wheeler
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
