[firedrake] Enforcing slip boundary condition strongly
David Ham
David.Ham at imperial.ac.uk
Fri Jul 10 10:25:21 BST 2015
Hi Colin,
That does not appear to be what Justin is asking for. He is asking how to
set a Dirichlet condition on one component of a VectorFunctionSpace. I'm
not sure we can do that right now, but I think this should be fixable. I
want to check what the Dolfin syntax for this is (so we maintain
compatibility) and then try to fix it. (Of course someone might want to
pleasantly surprise me and tell me we can do this already!)
David
On Fri, 10 Jul 2015 at 10:14 Cotter, Colin J <colin.cotter at imperial.ac.uk>
wrote:
> Hi Justin,
> Just to be clear, are you talking about enforcing normal component of
> the vector RT0 quantity to be zero?
>
> all the best
> --cjc
>
> On 9 July 2015 at 21:28, Justin Chang <jychang48 at gmail.com> wrote:
>
>> Hi everyone,
>>
>> I am solving a mixed poisson (aka: Darcy) problem using the VMS formalism
>> (Masud and Hughes, CMAME, 2002) which enables one to use equal order
>> lagrangian interpolation. In this formulation, the flux/velocity is
>> prescribed strongly and the pressure is prescribed weakly, like how it's
>> done with RT0 elements.
>>
>> The main difference is that velocity is now a VectorFunctionSpace as
>> opposed to a FunctionSpace. My question is, how would I implement the flux
>> boundary condition in a Dirichlet sense? Say I am working with a unit
>> square for now. I want to enforce v_x = some value(s) for x=0/x=1, and v_y
>> - some values(s) for y = 0/y = 1. How would this be done in Firedrake?
>>
>> Thanks,
>> Justin
>>
>> _______________________________________________
>> firedrake mailing list
>> firedrake at imperial.ac.uk
>> https://mailman.ic.ac.uk/mailman/listinfo/firedrake
>>
>>
>
>
> --
> http://www.imperial.ac.uk/people/colin.cotter
>
> www.cambridge.org/9781107663916
>
>
>
-------------- next part --------------
HTML attachment scrubbed and removed
More information about the firedrake
mailing list