Bug in patchInternalField()
I believe that there is a problem in implementation of the function patchInternalField().
The distribution is OpenFoam-1.7.1.
The problem was discovered while working on updateCoeffs() function of derivedFvPatchField (derived from mixedFvPatchField):
Instead, some cells contain zero value.
New boundary condition is applied only on the inner patch.
Anomalies appear only on the faces that are next to the edges, selective, not all. Usually several first ones (like in the example above).
If the "inner" patch is divided into 6 ones, with applied new derivedBC, anomalies are not observed.
Presented usage of the function is quite basic- this drives to conclusion that there is something wrong deeper in the code. Please this bug into consideration.
I would also be very grateful for any other possible explanation of the described problem.
Did this get resolved? Maybe you should report it directly to http://www.openfoam.com/bugs/
did you create a bug report for this? I have not been able to find one...
At the moment, I am facing the same problem (using OF 1.6-ext), as patch values calculated by patchInternalField() are different from internal field values in the patch near cells.
Or did you find a solution for this in the meantime?
Edit: Could someone check if the following is right or wrong:
Ok, having a deeper look at my results, I can so far state that the "bug" in my case only resulted from the interpolation/display scheme in ParaView:
Values from U.boundaryField()[patchi].patchInternalField() were different from U internalField values at regarded probe locations, because ParaView normally seems to use (mesh) point values instead of cell values for display. Using a CellDatatoPointData filter before setting the probe location in contrast results in identical values.
i am facing exactly the same Problem as Pawel does. I also get four wrong values as the first entries in my Field list.
Is there a Solution to this? Is it even a Bug or did i do something wrong?
this issue disappeared when I moved to OF-2.0.0. In the end I did not find out what caused it, but used a "walk around".
You can create a function and use it instead of patchInternalField
In the .H file:
Hope it helps a bit.
|All times are GMT -4. The time now is 06:53.|