View Issue Details

IDProjectCategoryView StatusLast Update
0000794OpenFOAMBugpublic2013-03-21 14:56
ReporterLieven Assigned Tohenry  
PriorityhighSeveritymajorReproducibilityalways
Status resolvedResolutionfixed 
Summary0000794: The coeff_ in anisotropicFilter appears to be calculated wrong.
DescriptionIf I understand correctly
coeff_.internalField().replace
(
    d,
    (2.0/widthCoeff_)*mesh.V()
    /fvc::surfaceSum(mag(mesh.Sf().component(d)))().internalField()
);

returns the average length of the cell in the d direction divided by widthCoeff_. To calculate the coeff_, seems to me that you need the square of this length. Hence the code should read

coeff_.internalField().replace
(
    d,
    1/widthCoeff_*
    sqr(
        2.0*mesh.V()
        /fvc::surfaceSum(mag(mesh.Sf().component(d)))().internalField()
    )
);
Additional InformationOther versions of OF show the same error (also 2.2.x)
TagsNo tags attached.

Activities

henry

2013-03-21 14:56

manager   ~0002032

Thanks for the bug-report. Your are quite correct and your suggested change is included in commit 6f9b0212f3bbae61303b1eb3e9d70fd55c847a25 in OpenFOAM-2.2.x

Issue History

Date Modified Username Field Change
2013-03-20 16:12 Lieven New Issue
2013-03-21 14:56 henry Note Added: 0002032
2013-03-21 14:56 henry Status new => resolved
2013-03-21 14:56 henry Resolution open => fixed
2013-03-21 14:56 henry Assigned To => henry