CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Bugs (
-   -   Bug in wedgePolyPatchC (

su_junwei February 12, 2009 08:09

Dear Foamers? There seems
Dear Foamers?

There seems to be a bug in the implementation of initTransforms() in wedgePolyPatch.C(line 44)

If the face[0] in this patch is in xy plane, patchNormal_ will equal {0 0 1} (after normal treatment) .
centerNormal_ is {0 0 0.5} before normal treatment and {0 0 1} after normal treatment.

While evaluating the axis_ using
axis_ = centreNormal_ ^ patchNormal_;
we get a zero vector for axis_.

In the following normal treatment for axis_, we will get a zero division exception.

It seems that the method for evaluation of centreNormal_ can't always give us a feasible result.

Is it right?

Junwei su

henry February 12, 2009 08:13

> If the face in this patch is
> If the face[0] in this patch is in xy plane

This should not be the case. Take a look at the documentation on how to setup wedge cases.


All times are GMT -4. The time now is 15:30.