dolfin team mailing list archive

["B. Emek Abali" <question153567@xxxxxxxxxxxxxxxxxxxxx>]

Question #153567 on DOLFIN changed:
https://answers.launchpad.net/dolfin/+question/153567

    Status: Open => Answered

B. Emek Abali proposed the following answer:
try to change the order of the boundary apply,

a better solution is prescribe left only horizontal velocities
bc1 = DirichletBC(system.sub(0).sub(0), Constant(2.0), left)

and leave the verticals as zero Neumann conditions,

or define it over the pressure as traction in the horizontal which is
actually the same as above, on the corner you will have some dispersion
which you wont see if the channel is long enough.

-- 
You received this question notification because you are a member of
DOLFIN Team, which is an answer contact for DOLFIN.