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 |
| Reporter | albertop | Assigned To | henry | ||
| Priority | normal | Severity | minor | Reproducibility | always |
| Status | resolved | Resolution | fixed | ||
| 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. | ||||
