View Issue Details

IDProjectCategoryView StatusLast Update
0003684OpenFOAMBugpublic2021-06-10 08:44
ReporteralbertopAssigned Tohenry 
PrioritynormalSeverityminorReproducibilityalways
Status resolvedResolutionfixed 
Product Version 
Fixed in Versiondev 
Summary0003684: Incorrect coefficients of polynomial in limiter
DescriptionThe coefficents a_ and b_ in cubicGradientLimiter.H do not satisfy some of the conditions on the derivatives from the original reference.

Expanding the polynomial that defines r leads to a r^3 + b^2 + r, which is zero and has derivative equal to 1 for r = 0. However, the behavior for r= rt is incorrect (the value of the polynomial should be 1 and its derivative should be zero).

The correct values of the coefficients are:

a = -2/rt^3 + 1/rt^2
b = -1/(2*rt) - 3*a*rt/2

Originally found by Shannon Leakey ( https://www.cfd-online.com/Forums/openfoam-programming-development/236503-polynomial-celllimited-cubic-seems-dodgy.html ). I have verified the results reported in the CFD-online post and are correct.
TagsNo tags attached.

Activities

henry

2021-06-10 08:44

manager   ~0012053

Resolved by commit eea13a52b2e02d32f5cccf82e59bff53b38c18d3

Issue History

Date Modified Username Field Change
2021-06-10 07:32 albertop New Issue
2021-06-10 08:44 henry Assigned To => henry
2021-06-10 08:44 henry Status new => resolved
2021-06-10 08:44 henry Resolution open => fixed
2021-06-10 08:44 henry Fixed in Version => dev
2021-06-10 08:44 henry Note Added: 0012053