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 #...
To upload designs, you'll need to enable LFS and have an admin enable hashed storage. More information