<html><head><style type='text/css'>p { margin: 0; }</style></head><body><div style='font-family: arial,helvetica,sans-serif; font-size: 12pt; color: #000000'>Hi Andreas, thank you for your reply.<br>I'm using GeometryGrid and ALU together. I provide ALU with a DiscreteDeformationFunction (my coarse element parametrization). Any time a new point appears in the mesh after a refinement I move that point using the DiscreteDeformationFunction.<br>The problem appears for local refinement. On the other hand, if I globally refine the mesh I have no problem about the divergence of the cells: it is zero everywhere.<br>However, if I refine only a bunch of cells, then the coarse ones, sharing intersection with the fine ones, have non-zero divergence.<br>I'm doing all this using the following methods:<br>intersection centerOuterUnitNormal<br>and codimension 1 volume method (choosing the element face corresponding to the intersection)<br>I hope it is clearer.<br>Thanks,<br><br>Marco<br><br><br><hr id="zwchr"><div style="color:#000;font-weight:normal;font-style:normal;text-decoration:none;font-family:Helvetica,Arial,sans-serif;font-size:12pt;"><b>Da: </b>"Andreas Dedner" <a.s.dedner@warwick.ac.uk><br><b>A: </b>dune@dune-project.org<br><b>Inviato: </b>Martedì, 24 febbraio 2015 12:14:33<br><b>Oggetto: </b>Re: [Dune] parametrized AMR and cell divergence<br><br>
I'm not quite sure I understand the question. Are you seeing
problems<br>
when refining the grid using ALU's standard globalRefine method or
do the problems<br>
only appear when you use your local parametrization and after moving
nodes around<br>
(I'm not quite clear how you do that with ALU....).<br>
Andreas<br>
<br>
<br>
<div class="moz-cite-prefix">On 24/02/15 11:00, Marco Cisternino
wrote:<br>
</div>
<blockquote cite="mid:637362868.4442210.1424775619676.JavaMail.root@optimad.it">
<style>p { margin: 0; }</style>
<div style="font-family: arial,helvetica,sans-serif; font-size:
12pt; color: #000000">Good morning,<br>
I would like to ask duners about deformed AMR.<br>
I'm working with body fitted mesh and ALUGrid.<br>
</div>
</blockquote>
<blockquote cite="mid:637362868.4442210.1424775619676.JavaMail.root@optimad.it">
<div style="font-family: arial,helvetica,sans-serif; font-size:
12pt; color: #000000">In order to refine my mesh I compute a
local parametrization for every coarse element, then I refine
the mesh moving new nodes using my parametrization.<br>
Therefore I compute the cell divergence summing over the cell
intersection<br>
normal . flux * area<br>
where intersection normals and areas are computed by Dune
methods centerUnitNormal and volume(this for a face), while flux
is (1,1,1) everywhere.<br>
I expect 0 over the entire mesh and this is true on the coarse
grid, but for a refined grid this is not true for coarse cells
sharing intersections with finer cells.<br>
Does anyone see something like this? And is there a cure? Am I
doing something wrong?<br>
<br>
Thanks for any hint.<br>
Bests,<br>
<br>
Marco<br>
<br>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre>_______________________________________________
Dune mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Dune@dune-project.org" target="_blank">Dune@dune-project.org</a>
<a class="moz-txt-link-freetext" href="http://lists.dune-project.org/mailman/listinfo/dune" target="_blank">http://lists.dune-project.org/mailman/listinfo/dune</a>
</pre>
</blockquote>
<br>
</div><br></div></body></html>