Hi again,
So this may or may not be related to the original issue. I was able to define corner IDs needed for my pressure BCs but I cannot seem to enforce both no-flow and atmospheric pressure simultaneously. Enforcing no-flow seems to negate the weakly prescribed pressure flux and my pressure solution is complete wack. On the otherhand, not enforcing no flow gives me a pressure profile I expect, but I no longer have a no flow domain.
Attached is my code, the custom boundary IDs begin at line 39, the weak formulations begin at line 164. Run the code as:
python 2D_Darcy_Homo_Iso.py 600 150 2
The code right now will generate a nonsensical pressure solution because I have no-flow enforced on all boundaries (including the corners).
Any ideas on what to do here?
Thanks,