If the angle between the boundary face normal Â and the normalized vector from the cell centroid to the boundary face centroid êb (skewness angle) is larger 90 • , the transfer coefficient κ turns to negative values resulting in numerical problems and a nonphysical solution. Therefore, the skewness angle is limited to 88 • , corresponding to Â • êb = 0.035, in the implemented calculation routine of the transfer coefficients. This enables simulations based on computational meshes containing boundary cells with skewness angle close to 90 • or larger. Nevertheless, for an accurate solution, it is recommended to prevent skewness angles larger than 88 • in the applied mesh. All meshes used in the presented study follow this recommendation.