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,
Justin

On Tue, Feb 7, 2017 at 4:55 AM, Lawrence Mitchell <lawrence.mitchell@imperial.ac.uk> wrote:

> On 7 Feb 2017, at 09:03, Lawrence Mitchell <lawrence.mitchell@imperial.ac.uk> wrote:
>
>
> PETSc.Options("foo_").getAll()
>
> => {"bar": "-baz_bar"}
>
> I will push a fix for this one.

Please try the branch proposed here:

https://github.com/firedrakeproject/firedrake/pull/1004

Lawrence
_______________________________________________
firedrake mailing list
firedrake@imperial.ac.uk
https://mailman.ic.ac.uk/mailman/listinfo/firedrake