Skip to content

GitLab

  • Projects
  • Groups
  • Snippets
  • Help
    • Loading...
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
    • Contribute to GitLab
  • Sign in / Register
D
dorie
  • Project overview
    • Project overview
    • Details
    • Activity
    • Releases
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributors
    • Graph
    • Compare
  • Issues 31
    • Issues 31
    • List
    • Boards
    • Labels
    • Service Desk
    • Milestones
  • Merge Requests 9
    • Merge Requests 9
  • Operations
    • Operations
    • Incidents
    • Environments
  • Analytics
    • Analytics
    • Repository
    • Value Stream
  • Wiki
    • Wiki
  • Members
    • Members
  • Activity
  • Graph
  • Create a new issue
  • Commits
  • Issue Boards
Collapse sidebar
  • dorie
  • dorie
  • Issues
  • #136

Closed
Open
Opened Jan 07, 2019 by Santiago Ospina De Los Ríos@sospinarMaintainer

Neumman BC only applies in the gravity direction

Description

Local operator for Richards applies Neumman BC only in the gravity direction. My question is why? Yeah, most of our study cases are in the vertical direction, but who says that there cannot be pumps or lateral fluxes under certain conditions? Am I missing something?

// update residual according to byType flag
if (bcType == BCType::Neumann)
{
  // flux is given parallel to gravity vector
  normal_flux *= std::abs( ig.centerUnitOuterNormal() * param.gravity() );

  // update residual (flux term)
  for (unsigned int i = 0; i<lfsv_s.size(); i++)
    r_s.accumulate(lfsv_s,i, normal_flux * phi_s[i] * factor);

  continue;
}

Related issues

See #...

Edited Jan 07, 2019 by Santiago Ospina De Los Ríos
To upload designs, you'll need to enable LFS and have an admin enable hashed storage. More information
Assignee
Assign to
None
Milestone
None
Assign milestone
Time tracking
None
Due date
None
Reference: dorie/dorie#136