HPRefinement

Netgen
1

Hi everyone!
I’m experimenting with the HPRefinement function in Netgen.
In particular, I’m trying to construct a Powell-Sabin split. To achieve this I’ve constructed a new refinement reference triangle, i.e.

int reftrig_Powell_splitedges[][3] =
{
  { 1, 2, 4 },
  { 2, 3, 5 },
  { 3, 1, 6 },
  { 0, 0, 0 },
};
int reftrig_Powell_splitfaces[][4] =
{
  { 1, 2, 3, 7 },
  { 0, 0, 0, 0 }
};

HPREF_ELEMENT_TYPE reftrig_Powell_newelstypes[] =
{
  HP_TRIG, HP_TRIG, HP_TRIG, HP_TRIG, HP_TRIG, HP_TRIG,
  HP_NONE, 
};
int reftrig_Powell_newels[][8] =
{
  { 1, 4, 7 }, 
  { 4, 2, 7 }, 
  { 2, 5, 7 },
  { 5, 3, 7 },
  { 3, 6, 7 },
  { 6, 1, 7 },
};
HPRef_Struct reftrig_Powell =
{
  HP_TRIG,
  reftrig_Powell_splitedges, 
  reftrig_Powell_splitfaces,
  0,
  reftrig_Powell_newelstypes,
  reftrig_Powell_newels
};

Now the structure I get out is more or less correct, but I don’t understand why the edge splits don’t match across two elements.
Any ideas?
image

2

Hi Umberto,

It looks like it is due to that line, which we should certainly remove:

Joachim

3

So removing L679 seem to create some issues with face splits, for example, it broke the Alfeld split. I’ve introduced a similar quantity known as fac2 to solve this issue. Now the Powell-Sabin refinement works ! I’ve opened a merge request in the TU Gitlab :slight_smile:

4

thanks Umberto, is merged and will be in the coming nightly,
Joachim

5

I tested the Powell-Sabin mesh and the mesh that comes looks great! Yet I think something is off with the new vertex label, in fact, if I try to solve the Poisson equation I get something absolutely wrong. I don’t understand what is going wrong tho because with the Alfeld split coded in the same way everything works. Any idea of what I could have done wrong?

image
image758×448 19.2 KB

6

can you share the example?

7

Sure !
test_powellsabin.py (875 Bytes)

8

fixed and committed, Joachim

first_problem_fixed
first_problem_fixed1446×1004 48.4 KB

1 Like
9

I think you posted the wrong image, it actually looks good now :wink:

image
image764×513 18 KB

1 Like
10

yes sure, thanks, I just found the “first_problem_solved” picture more interesting, and it clearly tells what was the second problem :slight_smile:

11

May I ask what were the first and second issues ?

12

the Powell-Sabin split was implemented via hp-Refinement (by you).

New vertices are defined by a weighted partitioning of edges. Shared vertices between two elements are identified by indices of the end-points. While the ordering of the end-points coincide by the “traditional” use of the hp-refinement, this was not the case here. Solution: When the new point is in the middle of the edge, identify them independently of ordering of old vertices.

Problem 2: Boundary elements (segments) were not refined. Thus, the Dirichlet-flag was set only for original vertices (as we can see in the first picture).

Now we are looking forward to see use-cases of the splitted mesh

Joachim

1 Like