On Fri, Sep 14, 2012 at 2:41 PM, A. A. Eftekhari <[email protected]> wrote: > Thanks a lot Daniel. I updated the script with more realistic values for > the Buckley Leverett equation. In fact it is not strange that fipy could > not sense the left boundary. We fixed saturation (sw) on the left side > of the domain, i.e., a face value, While the velocity value in the > convection term is calculated using the cell value, which is equal to > zero every where except at the left boundary. I fixed the value of a > cell next to the left boundary and it works now, even without the > diffusion term. We might as well have used the sw.faceValue as the > coefficient of the the convection term, which works but I don't like its > result.
Great. Thanks for getting back with your experience. I'm glad that it all seems to be working. Cheers. -- Daniel Wheeler _______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
