CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Bugs (http://www.cfd-online.com/Forums/openfoam-bugs/)
-   -   Bug in wedgePolyPatchC (http://www.cfd-online.com/Forums/openfoam-bugs/62318-bug-wedgepolypatchc.html)

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.

H


All times are GMT -4. The time now is 17:20.