Python BCs for mixed FE (as e.g. in the HM process)
Description
For problems with equal order shape functions (e.g. the pure M process) complex BCs can be nicely prescribed using the python interface. However, this is not possible for mixed finite elements (as e.g. required for the hydro-mechanics process).
Steps to Reproduce the Problem
- Compile the latest OGS-Version with the flag OGS_USE_PYTHON = ON.
- Run simHM_Bense_transient.prj
using the meshes:
basinmesh_domain.vtu
basinmesh_boundary_top.vtu
basinmesh_boundary_bottom.vtu
basinmesh_boundary_left.vtu
basinmesh_boundary_right.vtu
and the python files for the BCs: glacierBCsOGS.py glacierclass.py
Expected behavior: The simulation should run.
Actual behavior: The simulation starts, but then quits with:
error: No d.o.f. found for (node=54280, var=0, comp=0). That might be due to the use of mixed FEM ansatz functions, which is currently not supported by the implementation of Python BCs. That excludes, e.g., the HM process.
Can we extend the python BCs? That would be really helpful!
Specifications
- Version: 6.4
- Platform: Ubuntu 20.04.1 LTS