Finite Area Method: ZeroGradient Boundary condition changed
I use the zeroGradient boundary condition for outlets in thin film simulations (with faSavageHutterFoam from the avalanche module). This strategy worked with v1812 but does not after an update to v1912.
Steps to reproduce
Outlets, modeled as grad(Us) = 0, grad(h) = 0 in shallow water/thin film simulations will show distortions.
Run the appended case with OpenFOAM-v1912. It relies on the avalanche module.
This example is taken from section 6.2 in https://www.researchgate.net/publication/323184565_A_finite_area_scheme_for_shallow_granular_flows_on_three-dimensional_surfaces
Result of the example with OpenFOAM-v1812
Result of the example with OpenFOAM-v1912
Note the high peak of the height at the outlet.
What is the current bug behaviour?
The height grows monotonically near the outlet.
What is the expected correct behavior?
The height should reach a steady state near the outlet similar to the rest of the domain.
Relevant logs and/or images
- OpenFOAM version : v1912
- Operating system : Ubuntu (Linux Mint)
- Hardware info :
- Compiler : gcc
I expect the change in the zeroGradient boundary condition.