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.

2D Geometry with partly shared edge

More
4 years 3 weeks ago #2518 by creativeworker
2D Geometry with partly shared edge was created by creativeworker
Hello,

we would like to build a 2D Geometry (for example two nested rectangles) that represent one half of a symmetric problem, so the inner rectangle is completely shifted to the edge of the sourounding rectangle. At that moment we cannot mesh the geometry.
How could we implement this?
Code:
from netgen.geom2d import * from ngsolve import * geo = SplineGeometry() geo.AddRectangle ( (0,0), (2,2), bcs=["b","r","t","l"], leftdomain=1) geo.AddRectangle ( (0,0.5), (1,1.5), bcs=["b2","r2","t2","l2"], leftdomain=2, rightdomain=1) geo.SetMaterial (1, "outer") geo.SetMaterial (2, "inner") mesh = Mesh(geo.GenerateMesh(maxh=0.1))
More
4 years 3 weeks ago #2523 by mneunteufel
Replied by mneunteufel on topic 2D Geometry with partly shared edge
Hi creativeworker,

to construct such a mesh you need to specify the points and lines by hand
Code:
from netgen.geom2d import * from ngsolve import * geo = SplineGeometry() pnts = [ (0,0), (2,0), (2,2), (0,2), (0,1.5), (0,0.5), (1,0.5), (1,1.5)] pind = [ geo.AppendPoint(*pnt) for pnt in pnts ] geo.Append(['line',pind[0],pind[1]],leftdomain=1,rightdomain=0,bc="b") geo.Append(['line',pind[1],pind[2]],leftdomain=1,rightdomain=0,bc="r") geo.Append(['line',pind[2],pind[3]],leftdomain=1,rightdomain=0,bc="t") geo.Append(['line',pind[3],pind[4]],leftdomain=1,rightdomain=0,bc="l") geo.Append(['line',pind[4],pind[5]],leftdomain=2,rightdomain=0,bc="l") geo.Append(['line',pind[5],pind[0]],leftdomain=1,rightdomain=0,bc="l") geo.Append(['line',pind[5],pind[6]],leftdomain=2,rightdomain=1,bc="b2") geo.Append(['line',pind[6],pind[7]],leftdomain=2,rightdomain=1,bc="r2") geo.Append(['line',pind[7],pind[4]],leftdomain=2,rightdomain=1,bc="t2") geo.SetMaterial (1, "outer") geo.SetMaterial (2, "inner") mesh = Mesh(geo.GenerateMesh(maxh=0.1))

Best,
Michael
The following user(s) said Thank You: creativeworker
Time to create page: 0.137 seconds

Login