Skip to content

Fix bug at the domain's edges caused by incorrect neighbors

There was a bug at the edges of the fluid domain in positive coordinate direction:

2022_08_17_BugAnKante1_t2460

For certain gird sizes the fluid nodes on the edges had invalid neighbors. To make them visible we added a function which writes the neighbors to a vtk files: NeighborDebugWriter::writeNeighborLinkLinesDebug(para.get()); There is a file for each of the 27 directions. The file contains lines which connect each node to the neighbor in the respective direction. The image below shows the neighbors in PPO-direction: 2022_08_22_Neighbor_PP0_NE_bug

To fix this bug void GridImp::setStopperNeighborCoords(uint index) was modified:

vf::Math::lessEqual(y + delta, endY) was changed to vf::Math::lessEqual(y + delta, endY + (0.5 * delta))

The image below visualizes the new version after the bugfix:

2022_09_02_Neighbor_of_Stopper_Test

Edited by Anna Wellmann

Merge request reports