New question #140910 on DOLFIN:
https://answers.launchpad.net/dolfin/+question/140910
setting a mixed space (velocity field v, pressure field p) as
V = VectorFunctionSpace(mesh, 'CG', 1)
Q = FunctionSpace(mesh, 'DG', 0)
ME = MixedFunctionSpace([V,Q])
and by describing
null = Constant(0.0)
class TopDrag(Expression):
def eval(self, vec, x):
vec[0] = 10.0
vec[1] = 0.0
vec[2] = 0.0
def dim(self):
return 3
vOnTop = TopDrag()
wanted to reach:
a given velocity on a subdomain (not here but defined in a proper way with
compile_subdomains() )
DirichletBC(ME.sub(0), vOnTop, top)
zero velocity in one direction
DirichletBC(ME.sub(0).sub(1), null, front)
and a given pressure on a subdomain
DirichletBC(ME.sub(1), pIn, left)
which I couldn't :)
Any sharp eyes where I am not thinking as FEniCS wants?
--
You received this question notification because you are a member of
DOLFIN Team, which is an answer contact for DOLFIN.
_______________________________________________
Mailing list: https://launchpad.net/~dolfin
Post to : [email protected]
Unsubscribe : https://launchpad.net/~dolfin
More help : https://help.launchpad.net/ListHelp
