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@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@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@imperial.ac.uk
https://mailman.ic.ac.uk/mailman/listinfo/firedrake




--