View Issue Details
|ID||Project||Category||View Status||Date Submitted||Last Update|
|0003684||OpenFOAM||Bug||public||2021-06-10 07:32||2021-06-10 08:44|
|Fixed in Version||dev|
|Summary||0003684: Incorrect coefficients of polynomial in limiter|
|Description||The 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.
|Tags||No tags attached.|