Skip to content
GitLab
Projects Groups Snippets
  • /
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
    • Contribute to GitLab
  • Register
  • Sign in
  • dumux dumux
  • Project information
    • Project information
    • Activity
    • Labels
    • Members
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributors
    • Graph
    • Compare
  • Issues 102
    • Issues 102
    • List
    • Boards
    • Service Desk
    • Milestones
  • Merge requests 91
    • Merge requests 91
  • CI/CD
    • CI/CD
    • Pipelines
    • Jobs
    • Schedules
  • Deployments
    • Deployments
    • Environments
    • Releases
  • Packages and registries
    • Packages and registries
    • Container Registry
  • Monitor
    • Monitor
    • Incidents
  • Analytics
    • Analytics
    • Value stream
    • CI/CD
    • Repository
  • Wiki
    • Wiki
  • Activity
  • Graph
  • Create a new issue
  • Jobs
  • Commits
  • Issue Boards
Collapse sidebar
  • dumux-repositoriesdumux-repositories
  • dumuxdumux
  • Issues
  • #934
Closed
Open
Issue created Sep 29, 2020 by Kilian Weishaupt@kweisMaintainer

[periodic] Check if current implementation for box is correct

For box, we consider a "primary" and "secondary" DOFs, i.e. the secondary DOF is periodically corresponding to the primary one.

We add the value of the secondary DOF residual to the one of the primary one

 res[m.first] += res[m.second];

and enforce

jac[m.first][it.index()] += (*it)

This makes sense, I guess.

For the secondary DOF, we adapt the matrix entries such that the is a -1 for the periodic DOF and a 1 for the own DOF. All other entries are zero. I guess the idea is to set the secondary DOF values equal to the primary ones.

I think we should also adapt the residual of the secondary DOF such that

res[m.second] = curSol[m.second] - curSol[m.first];
Edited Sep 29, 2020 by Timo Koch
Assignee
Assign to
Time tracking