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
[email protected]
http://www.ctcms.nist.gov/fipy
[ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]