Skip to content

faPatch::makeCorrectionVectors fails, depending on the mesh decomposition

Summary

Running the tutorial under $FOAM_TUTORIALS/incompressible/adjointOptimisationFoam/sensitivityMaps/motorBike might crash in faPatch::makeCorrectionVectors with a division by zero, depending on the mesh decomposition.

In specific, running the case in 88 processors (attached) causes the issue while running the one under the tutorials folder (20 processors) does not.

Steps to reproduce

Download and run the attached case (apologies for the relatively large core count but it's the only available setup that reproduces the issue). The issue appears when solving the sensitivity smoothing PDEs on the faMesh. To speed up testing, I have set the primal and adjoint iterations to a very small number. If the default laplaceScheme is set to Gauss uncorrected; in faSchemes, the problem does not manifest.

motorBike_88.tar.gz

Relevant logs and/or images

See attached the last bit of the log file including the error.

err

Environment information

  • OpenFOAM version : v2312
  • Operating system : CentOS 7.9
  • Compiler : gcc 9.3

/label bug