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 define a codimention-1 finite element space on the interface

More
4 years 4 months ago #2196 by Guosheng Fu
How to define a codimention-1 finite element space on the interface was created by Guosheng Fu
Hello,

I have a (2D) domain with a closed fitted interface Gamma, and I would like to define a finite element space only on the interface. Then I will use this finite element space to solve the Laplace-Beltrami equation on this interface.

I was thinking of defining a global H1 space on the whole domain, then set false all DOFs away from the interface, and then compress the space. Is this a correct way for the set-up? Or is there more efficient alternatives?

Best,
Guosheng
More
4 years 4 months ago #2197 by mneunteufel
Replied by mneunteufel on topic How to define a codimention-1 finite element space on the interface
Hi Guosheng ,

your approach should work. It is also possible to define the space only on the interface

Code:
fes = H1(mesh, order=2, definedon="", definedonbound="interface")

However, to draw the solution you need to interpolate it into a global H1 space (see attached code)

Best,
Michael

File Attachment:

File Name: defonbnd.py
File Size:1 KB
Attachments:
More
4 years 4 months ago #2205 by Guosheng Fu
Replied by Guosheng Fu on topic How to define a codimention-1 finite element space on the interface
Thanks, Michael.
Time to create page: 0.166 seconds

Login