The one place for your designs
To enable design management, you'll need to meet the requirements. If you need help, reach out to our support team for assistance.
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 #...