FiPy has issues with non-orthogonal meshes (i.e. when the vector between the cell center isn't parallel to the face normal). We did make an attempt at fixing this, which resulted in two diffusion terms,
https://github.com/usnistgov/fipy/blob/develop/fipy/terms/diffusionTermCorrection.py and https://github.com/usnistgov/fipy/blob/develop/fipy/terms/diffusionTermNoCorrection.py As I recall, the "corrected" diffusion term was quite unstable for moderate to high non-orthogonality. I'd like to address this issue again and I was wondering if anyone had opinions regarding the latest / best / most tractable approach for dealing with non-orthogonality in CC-FV. I haven't followed the literature on this for quite some time. -- Daniel Wheeler _______________________________________________ fipy mailing list fipy@nist.gov http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]