Skip to content
GitLab
Projects Groups Snippets
  • /
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
    • Contribute to GitLab
  • Sign in / Register
  • openfoam openfoam
  • Project information
    • Project information
    • Activity
    • Labels
    • Members
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributors
    • Graph
    • Compare
  • Issues 426
    • Issues 426
    • List
    • Boards
    • Service Desk
    • Milestones
  • Merge requests 9
    • Merge requests 9
  • Deployments
    • Deployments
    • Releases
  • Wiki
    • Wiki
  • Activity
  • Graph
  • Create a new issue
  • Commits
  • Issue Boards
Collapse sidebar
  • Development
  • openfoamopenfoam
  • Issues
  • #2113
Closed
Open
Issue created Jun 04, 2021 by Shannon Leakey@scleakey

cellLimited<cubic> limiter is not differentiable

The cubic polynomial in cubicGradientLimiter.H is

((ar + b)r + 1)r \\
a = \frac{2}{r_t^2} - \frac{2}{r_t^3} \\
b = -\frac32 a r_t

but surely it should be

((ar + b)r + 1)r \\
a = \frac{1}{r_t^2} - \frac{2}{r_t^3} \\
b = -\frac32 a r_t - \frac{1}{2r_t}

to satisfy the constraints from Michalak and Ollivier-Gooch (2008, p.5), specifically that the derivative at the transition point r_t should be zero so it joins up smoothly with 1. This problem is illustrated in the graph below:

cubic

Linking this cfd-online thread here just in case someone replies to it after I create this issue.

Assignee
Assign to
Time tracking