Abstract
This paper investigates how OpenFOAM matrix solver works based on mesh information of cell and face indices for each cell. The smoother code of Gauss-Seidel in OpenFOAM uses minimum information of start- and end-indices of owner face for each cell to make it efficient in computational time, but can give rise to incorrect calculation of matrix solver in a grid construction. To complement the weakness of Gauss-Seidel smoother, OpenFOAM provides diagonal incomplete LU (DILU) decomposition smoother. The DILU smoother is robust to the face-ordering but spends a lot of computational time to solve the matrix due to an inefficient usage of owner cell and neighbour cell for all faces. This paper proposes a robust modification of Gauss-Seidel smoother in OpenFOAM maintaining its computational efficiency. The present essence is to get perfect information of adjacent cell indices regardless of face ordering.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.