Forum Message

 

 

We have moved the forum to https://forum.ngsolve.org . This is an archived version of the topics until 05/05/23. All the topics were moved to the new forum and conversations can be continued there. This forum is just kept as legacy to not invalidate old links. If you want to continue a conversation just look for the topic in the new forum.

Notice

The forum is in read only mode.

How to set up boundary condition for VectorL2FESpace

More
3 years 11 months ago - 3 years 11 months ago #2693 by dong
I tried to set up u=g on the boundary for Stokes problem. I defined the DG FEspace as follows.
Code:
V = VectorL2(mesh,order=2, dirichlet ='left|right|top|bottom',dgjumps=True) Q = L2(mesh,order=1) X = FESpace([V,Q],dgjumps=True)

When I run the code, there is a error at the following line

gfu.components[0].Set(g, BND, definedon=mesh.Boundaries("left|right|top|bottom"))

VectorL2FESpace does not have an evaluator for BND!


The code was attached.
Could you please show me how to set up the boundary condition for VectorL2FESpace or fix this error? Thank you so much.
Last edit: 3 years 11 months ago by dong.
More
3 years 11 months ago #2694 by Guosheng Fu
l2 space doesnot have dofs on the bdry. you need to impose dirichlet bc weakly
More
3 years 11 months ago #2695 by dong
I'm a newbie, I don't understand what you've said. Could you please tell me what 'weakly' meant here? Do you know any examples about imposing dirichlet bc weakly with Ngsolve? Thank you.
Time to create page: 0.132 seconds