<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style id="owaParaStyle" type="text/css">
<!--
p
{margin-top:0;
margin-bottom:0}
p
{margin-top:0;
margin-bottom:0}
p
{margin-top:0;
margin-bottom:0}
p
{margin-top:0;
margin-bottom:0}
p
{margin-top:0;
margin-bottom:0}
p
{margin-top:0;
margin-bottom:0}
p
{margin-top:0;
margin-bottom:0}
p
{margin-top:0;
margin-bottom:0}
-->
P {margin-top:0;margin-bottom:0;}</style>
</head>
<body ocsi="0" fpstyle="1">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Hi Andreas,<br>
<br>
yes, thanks, now it is perfectly clear.<br>
I will put the interval just a big bigger than the actual boundary, so the comparison is safe.<br>
<br>
Best,<br>
<div>Andrea<br>
<div style="font-family:Tahoma; font-size:13px">__________________________________________________________<br>
<br>
Andrea Sacconi<br>
PhD student, Applied Mathematics<br>
AMMP Section, Department of Mathematics, Imperial College London,<br>
London SW7 2AZ, UK<br>
a.sacconi11@imperial.ac.uk<br>
</div>
</div>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div style="direction: ltr;" id="divRpF430852"><font face="Tahoma" color="#000000" size="2"><b>From:</b> Dedner, Andreas [A.S.Dedner@warwick.ac.uk]<br>
<b>Sent:</b> 16 May 2013 10:56<br>
<b>To:</b> Sacconi, Andrea; dune@dune-project.org<br>
<b>Subject:</b> RE: Boundary IDs of the faces of a cube within the DGF file<br>
</font><br>
</div>
<div></div>
<div>
<div style="direction:ltr; font-family:Tahoma; color:#000000; font-size:10pt">Okay, if one boundary of your cube is a=(0,0,0) to b=(10,10,0) then you could<br>
use that in the boundarydomain and it would probably work but it would have to <br>
involve a double comparison of the point generated form the interval block to see<br>
if they lie inside the interval [a,b]. So the domainblock does not generate any boundary<br>
one just provides an interval [a,b] and for each boundary segement generated by the<br>
interval block (or any other process for that matter) a test is performed to see if all<br>
the vertices of the segment fall within this interval. Using the exact boundaries to define<br>
this domain is a bit dangerous since the double comparison can go wrong. <br>
That's why in the example I extend the interval - thats why you have the -1 and 11 (-10, 100 or -0.1 and 10.1 would work just as well). In fact one should probably use<br>
(-1, -1, -1) to (11, 11, 0.0001)...<br>
Does that explain the process?<br>
Best<br>
Andreas<br>
<br>
<div style="font-family:Times New Roman; color:#000000; font-size:16px">
<hr tabindex="-1">
<div id="divRpF85323" style="direction:ltr"><font face="Tahoma" color="#000000" size="2"><b>From:</b> Sacconi, Andrea [a.sacconi11@imperial.ac.uk]<br>
<b>Sent:</b> 16 May 2013 10:45<br>
<b>To:</b> Dedner, Andreas; dune@dune-project.org<br>
<b>Subject:</b> RE: Boundary IDs of the faces of a cube within the DGF file<br>
</font><br>
</div>
<div></div>
<div>
<div style="direction:ltr; font-family:Tahoma; color:#000000; font-size:10pt">
<div style="direction:ltr; font-family:Tahoma; color:#000000; font-size:10pt">
<div style="direction:ltr; font-family:Tahoma; color:#000000; font-size:10pt">Hi Andreas, hi all,<br>
<br>
yes, I had a look at the documentation, and since I didn't find what I was looking for I posted my question.<br>
<br>
Let's have a look at the example you pointed: examplegrid6.dgf<br>
<div>The grid is a cube, with dimension 10, precisely from the origin (0, 0, 0) to (10, 10, 10).<br>
The comment says that you want the bottom boundary to have ID 2, and 1 for all the others.<br>
The bottom boundary is a square with South-West corner (0, 0, 0) and North-East corner (10, 10, 0).<br>
<br>
Why do you use the interval from (-1, -1, -1) to (11, 11, 0) ?? Why -1 and 11?<br>
Does this command work exactly as the block "Interval" used to generate the grid? I.e., we generate a cuboid, not a segment with extremes (A, B, C) and (D, E, F). If so, we would have a cuboid, with dimensions 12, 12, 1.<br>
Would it be the same if I used an interval from (0, 0, 0) to (10, 10, 0)?
<div style="font-family:Tahoma; font-size:13px"><br>
Hope that now my questions are more clear.<br>
Thanks again for your help!<br>
<br>
Andrea<br>
__________________________________________________________<br>
<br>
Andrea Sacconi<br>
PhD student, Applied Mathematics<br>
AMMP Section, Department of Mathematics, Imperial College London,<br>
London SW7 2AZ, UK<br>
a.sacconi11@imperial.ac.uk<br>
</div>
</div>
<div style="font-family:Times New Roman; color:#000000; font-size:16px">
<hr tabindex="-1">
<div id="divRpF883834" style="direction:ltr"><font face="Tahoma" color="#000000" size="2"><b>From:</b> Dedner, Andreas [A.S.Dedner@warwick.ac.uk]<br>
<b>Sent:</b> 16 May 2013 09:25<br>
<b>To:</b> Sacconi, Andrea; dune@dune-project.org<br>
<b>Subject:</b> RE: Boundary IDs of the faces of a cube within the DGF file<br>
</font><br>
</div>
<div></div>
<div>
<div style="direction:ltr; font-family:Tahoma; color:#000000; font-size:10pt">Hi Andrea.<br>
Have you had a look at the documentation (e.g.<span class="preprocessor"> examplegrid6.dgf</span>)<br>
<a href="http://www.dune-project.org/doc/doxygen/dune-grid-html/group___dune_grid_format_parser.html" target="_blank">http://www.dune-project.org/doc/doxygen/dune-grid-html/group___dune_grid_format_parser.html</a><br>
If you find that unclear then please describe in more detail what is unclear to you<br>
so I can adjust the examples and the documentation.<br>
Best<br>
Andreas<br>
<div style="font-family:Times New Roman; color:#000000; font-size:16px">
<hr tabindex="-1">
<div id="divRpF144098" style="direction:ltr"><font face="Tahoma" color="#000000" size="2"><b>From:</b> Sacconi, Andrea [a.sacconi11@imperial.ac.uk]<br>
<b>Sent:</b> 14 May 2013 13:50<br>
<b>To:</b> Dedner, Andreas; dune@dune-project.org<br>
<b>Subject:</b> RE: Boundary IDs of the faces of a cube within the DGF file<br>
</font><br>
</div>
<div></div>
<div>
<div style="direction:ltr; font-family:Tahoma; color:#000000; font-size:10pt">Hi all,<br>
<br>
thanks Andreas for your answer.<br>
I had a look at the boundarydomain, and the documentation says: "Each line consists of an integer greater than zero and two vectors describing an interval in dimworld space."<br>
So, and I am not sure if it is directly possible; if you have a 3D cuboid tessellated with tetrahedra, and you want all the triangles lying on one of the faces of the cuboid to be identified by the same boundary ID, how do you do that, without enumerating all
the boundary faces of the elements involved?<br>
<br>
Thanks again!<br>
<div>Andrea<br>
<div style="font-family:Tahoma; font-size:13px">__________________________________________________________<br>
<br>
Andrea Sacconi<br>
PhD student, Applied Mathematics<br>
AMMP Section, Department of Mathematics, Imperial College London,<br>
London SW7 2AZ, UK<br>
a.sacconi11@imperial.ac.uk<br>
</div>
</div>
<div style="font-family:Times New Roman; color:#000000; font-size:16px">
<hr tabindex="-1">
<div id="divRpF328008" style="direction:ltr"><font face="Tahoma" color="#000000" size="2"><b>From:</b> Dedner, Andreas [A.S.Dedner@warwick.ac.uk]<br>
<b>Sent:</b> 13 May 2013 14:41<br>
<b>To:</b> Sacconi, Andrea; dune@dune-project.org<br>
<b>Subject:</b> RE: Boundary IDs of the faces of a cube within the DGF file<br>
</font><br>
</div>
<div></div>
<div>
<div style="direction:ltr; font-family:Tahoma; color:#000000; font-size:10pt">Hi.<br>
You can use the interval block to define your domain (like in 2d) but that has nothing to do with the<br>
boundary ids. Have a look at the boundarydomain block which should do exactly what you need<br>
(it also uses an "interval" to define a region with the same boundary id).<br>
Best<br>
Andreas<br>
<div style="font-family:Times New Roman; color:#000000; font-size:16px">
<hr tabindex="-1">
<div id="divRpF85542" style="direction:ltr"><font face="Tahoma" color="#000000" size="2"><b>From:</b> dune-bounces+a.s.dedner=warwick.ac.uk@dune-project.org [dune-bounces+a.s.dedner=warwick.ac.uk@dune-project.org] on behalf of Sacconi, Andrea [a.sacconi11@imperial.ac.uk]<br>
<b>Sent:</b> 13 May 2013 14:37<br>
<b>To:</b> dune@dune-project.org<br>
<b>Subject:</b> [Dune] Boundary IDs of the faces of a cube within the DGF file<br>
</font><br>
</div>
<div></div>
<div>
<div style="direction:ltr; font-family:Tahoma; color:#000000; font-size:10pt">Hi everybody,<br>
<br>
I would like to ask you a question about how to prescribe boundary IDs within a DGF file.<br>
I just had a closer look at the DGF format, but I couldn't find how to do what I need.<br>
<br>
Let's say that the domain for your computations is a cube (or a rectangular cuboid, it doesn't matter), and you want to impose different boundary conditions on the six faces.<br>
In order to do that, I would introduce different boundary IDs to be read in the code, to guarantee that the faces are treated accordingly.<br>
<br>
If you want to assign the same ID to all the faces of elements lying on a certain big face of the cube, is it possible to do it directly, without enumerating all those faces?<br>
I mean, in 2D there is the excellent command "Interval".<br>
What happens if I use the keyword "Interval" (in 3D), passing the four segments forming the perimeter of the (big) face?<br>
<br>
Thanks in advance! <br>
<br>
<div>Andrea<br>
<div style="font-family:Tahoma; font-size:13px">__________________________________________________________<br>
<br>
Andrea Sacconi<br>
PhD student, Applied Mathematics<br>
AMMP Section, Department of Mathematics, Imperial College London,<br>
London SW7 2AZ, UK<br>
a.sacconi11@imperial.ac.uk<br>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>