[firedrake] Enforcing slip boundary condition strongly

[Justin Chang]

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
-------------- next part --------------
HTML attachment scrubbed and removed