Maybe some bug in leastSquares
Maybe some bug in leastSquaresVectors.c
some error in calculation of dd in leastsquare vector
void Foam::leastSquaresVectors::makeLeastSquaresVectors () const
const vector& d = pd[patchFacei];
dd maybe done by the following:
also a question:
for the couple boundary,dd is calculated only for P cell(owner),why not for N cell(neighbour)?
Yes you are right this is a bu
Yes you are right this is a bug, the second use of pw[patchFacei] is as (1 - pw[patchFacei]) but the first is inconsistent and incorrect. Here is the corrected http://www.cfd-online.com/OpenFOAM_D...hment_icon.gif leastSquaresVectors.C to replace the file in OpenFOAM-1.4.1/src/finiteVolume/finiteVolume/gradSchemes/leastSquaresGrad
On the coupled boundaries only the owner contribution needs to be calculated, the neighbour contribution being calculated by the "other-side" of the couple.
Thanks for the bug report and proposed solution
|All times are GMT -4. The time now is 18:23.|