HM domain deactivation
Description
Using deactivated_subdomains to disable the H part in HM for a layer in a model, results in enormous displacement fields, consequently wrecking the convergence. That is, if the deactivated layer is only one element in height. If there are more elements between the two hydraulic active layers, it produces nan values for the displacement solution. In THM the same functionality works as expected.
Steps to Reproduce the Problem
- This model without domain deactivation runs fine: HM_test.tar
- Uncomment the domain_deaction, now the model won't converge due to unreasonable displacements
Specifications
- Version: ogs version: 6.4.3-966-g68cd5a0b
- Platform: Ubuntu 22.04.1 LTS (GNU/Linux 5.10.102.1-microsoft-standard-WSL2 x86_64)