[firedrake] Prescribe pressure at a corner for Darcy's equation

[Justin Chang]

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 at imperial.ac.uk> wrote:

>
> > On 7 Feb 2017, at 09:03, Lawrence Mitchell <lawrence.mitchell at 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 at imperial.ac.uk
> https://mailman.ic.ac.uk/mailman/listinfo/firedrake
>
-------------- next part --------------
HTML attachment scrubbed and removed
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 2D_Darcy_Homo_Iso.py
Type: text/x-python
Size: 7948 bytes
Desc: not available
URL: <http://mailman.ic.ac.uk/pipermail/firedrake/attachments/20170213/c1270df5/attachment.py>