<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<meta name="Generator" content="Microsoft Exchange Server">
<!-- converted from text --><style><!-- .EmailQuote { margin-left: 1pt; padding-left: 4pt; border-left: #800000 2px solid; } --></style>
</head>
<body>
<meta content="text/html; charset=UTF-8">
<style type="text/css" style="">
<!--
p
{margin-top:0;
margin-bottom:0}
-->
</style>
<div dir="ltr">
<div id="x_divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:#000000; font-family:'Times New Roman',Times,serif">
<p>Hello Ruth Vazquez ,</p>
<p><br>
</p>
<p>Thank you very much. If I use the same procedure (boolean fragments) to enclose a concentric cylinders within </p>
<p>a cube, it results in 2D meshes. Can you please let me know what I am doing wrong here ?</p>
<p><br>
</p>
<p></p>
<div>SetFactory("OpenCASCADE");</div>
<div>Mesh.Algorithm3D = 4; // 3D mesh algorithm (1=Delaunay, 4=Frontal, 5=Frontal Delaunay, 6=Frontal Hex, 7=MMG3D, 9=R-tree)</div>
<div>R = 0.01;</div>
<div>R1 = 0.05;</div>
<div>w = R;</div>
<div>Box(1) = {-R1,-R1,-R1, 2*R1,2*R1,2*R1};</div>
<div><br>
</div>
<div>Cylinder(2) = {0, 0, -0.9*R1, 0, 0, 1.6*R1, 0.2*R, 2*Pi}; // inner cylinder</div>
<div>Cylinder(3) = {0, 0, -0.9*R1, 0, 0, 1.6*R1, 0.25*R, 2*Pi}; // inner cylinder</div>
<div>Cylinder(4) = {0, 0, -0.9*R1, 0, 0, 1.6*R1, 0.9*R, 2*Pi}; // inner cylinder</div>
<div>Cylinder(5) = {0, 0, -0.9*R1, 0, 0, 1.6*R1, 1*R, 2*Pi}; //outer cylinder</div>
<div><br>
</div>
<div>aa() = BooleanFragments{ Volume{1}; Delete; }{ Volume{2,3,4,5}; Delete;};</div>
<div><br>
</div>
<div>Physical Volume(600) = {aa(0)};</div>
<div>Physical Volume(601) = {aa(1)};</div>
<div>Physical Volume(602) = {aa(2)};</div>
<div>Physical Volume(603) = {aa(3)};</div>
<div>Physical Volume(604) = {aa(4)};</div>
<div><br>
</div>
<div>Mesh.CharacteristicLengthMax = 0.001;</div>
<br>
<p></p>
<p><br>
</p>
<div id="x_Signature">
<div name="x_divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<div name="x_divtagdefaultwrapper" style="margin:0px"><font face="Times New Roman, Times, serif" color="#0000ff"><br>
</font></div>
<div name="x_divtagdefaultwrapper" style="margin:0px">
<hr>
</div>
<blockquote style="margin:0 0 0 40px; border:none; padding:0px">
<div name="x_divtagdefaultwrapper" style="margin:0px"><font face="Times New Roman, Times, serif">Best Regards,</font></div>
<div name="x_divtagdefaultwrapper" style="margin:0px"><font face="Times New Roman, Times, serif"><br>
</font></div>
<div name="x_divtagdefaultwrapper" style="margin:0px"><span style="font-family:"Times New Roman",Times,serif">Sathyanarayan Rao, PhD student</span></div>
<div name="x_divtagdefaultwrapper" style="margin:0px"><span style="font-family:"Times New Roman",Times,serif; background-color:rgb(255,255,255); white-space:pre-wrap">Earth and Life Institute/Environmental Sciences (ELI-e)</span></div>
<div name="x_divtagdefaultwrapper" style="margin:0px"><span style="background-color:rgb(255,255,255); font-family:"Times New Roman",Times,serif; white-space:pre-wrap">Université catholique de Louvain</span></div>
<div name="x_divtagdefaultwrapper" style="margin:0px"><span style="font-family:"Times New Roman",Times,serif; white-space:pre-wrap; background-color:rgb(255,255,255)">c.037, Croix du Sud 2, Louvain-la-Neuve, Belgium
</span></div>
<div name="x_divtagdefaultwrapper" style="margin:0px"><span style="background-color:rgb(255,255,255); white-space:pre-wrap">Phone: 010473827 ( intercom 73827)</span></div>
</blockquote>
<div name="x_divtagdefaultwrapper" style="margin:0px">
<div>
<pre cols="72" style="white-space:pre-wrap; background-color:rgb(255,255,255)"><hr></pre>
</div>
</div>
</div>
</div>
</div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> gmsh <gmsh-bounces@ace20.montefiore.ulg.ac.be> on behalf of gmsh-request@ace20.montefiore.ulg.ac.be <gmsh-request@ace20.montefiore.ulg.ac.be><br>
<b>Sent:</b> 09 February 2018 14:39:01<br>
<b>To:</b> gmsh@onelab.info<br>
<b>Subject:</b> gmsh Digest, Vol 181, Issue 20</font>
<div> </div>
</div>
</div>
<font size="2"><span style="font-size:10pt;">
<div class="PlainText">Send gmsh mailing list submissions to<br>
gmsh@onelab.info<br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="http://onelab.info/mailman/listinfo/gmsh">http://onelab.info/mailman/listinfo/gmsh</a><br>
or, via email, send a message with subject or body 'help' to<br>
gmsh-request@onelab.info<br>
<br>
You can reach the person managing the list at<br>
gmsh-owner@onelab.info<br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of gmsh digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
1. Re: Create a 3D cube with each face having a physical surface<br>
(Ruth Vazquez Sabariego)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Fri, 9 Feb 2018 13:34:20 +0000<br>
From: Ruth Vazquez Sabariego <ruth.sabariego@kuleuven.be><br>
To: Sathyanarayan Rao <sathyanarayan.rao@uclouvain.be><br>
Cc: "gmsh@onelab.info" <gmsh@onelab.info><br>
Subject: Re: [Gmsh] Create a 3D cube with each face having a physical<br>
surface<br>
Message-ID: <00C42B21-C06B-4419-AED8-9752CD5A868B@kuleuven.be><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Using your second option, you may use BooleanFragments.<br>
Try at the end of your file:<br>
<br>
// Surface{7,8,9,10,11,12} In Volume{1}; // Comment out this line<br>
Physical Volume(100) = {1};<br>
<br>
aux() = BooleanFragments{ Volume{1}; Delete; }{ Surface{7:12}; Delete;};<br>
Printf("", aux()); // For checking the numbers of the generated geometrical entities<br>
<br>
// For visibility<br>
Recursive Color Red{Volume{1};}<br>
Recursive Color Yellow{Surface{7:12};}<br>
<br>
With regard to your first option, you can only create rectangles in the XY plane, but you can translate and rotate to place them where you want.<br>
<br>
HTH<br>
Ruth<br>
<br>
?<br>
Prof. Ruth V. Sabariego<br>
KU Leuven<br>
Dept. Electrical Engineering ESAT/Electa, EnergyVille<br>
<a href="http://www.esat.kuleuven.be/electa">http://www.esat.kuleuven.be/electa</a><br>
<a href="http://www.energyville.be">http://www.energyville.be</a><br>
<br>
Free software: <a href="http://gmsh.info">http://gmsh.info</a> | <a href="http://getdp.info">
http://getdp.info</a> | <a href="http://onelab.info">http://onelab.info</a><br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
On 8 Feb 2018, at 22:13, Sathyanarayan Rao <sathyanarayan.rao@uclouvain.be<mailto:sathyanarayan.rao@uclouvain.be>> wrote:<br>
<br>
Dear GMSH team,<br>
<br>
I am trying to generate a cube with each face having a rectangular surface that will be later identified as a sensor.<br>
<br>
1) In case I want to use rectangle command of opencascade is it possible to create rectangle other than in XY plane ?<br>
<br>
<br>
<br>
```<br>
SetFactory("OpenCASCADE");<br>
Mesh.Algorithm3D = 4; // 3D mesh algorithm (1=Delaunay, 4=Frontal, 5=Frontal Delaunay, 6=Frontal Hex, 7=MMG3D, 9=R-tree)<br>
R = 0.05;<br>
w = R;<br>
Box(1) = {-R,-R,-R, 2*R,2*R,2*R};<br>
Rectangle(10) = {-R/2,-R/2,-R, R,R};<br>
Rectangle(11) = {-R/2,-R/2,R, R,R};<br>
```<br>
<br>
2) When I manually create surfaces in each face of the cube, they come as separate regions in 3D mesh.<br>
<br>
```<br>
SetFactory("OpenCASCADE");<br>
Mesh.Algorithm3D = 4; // 3D mesh algorithm (1=Delaunay, 4=Frontal, 5=Frontal Delaunay, 6=Frontal Hex, 7=MMG3D, 9=R-tree)<br>
R = 0.05;<br>
R1 = R;<br>
w = R;<br>
Box(1) = {-R,-R,-R, 2*R,2*R,2*R};<br>
q = 0.005;<br>
<br>
/* Right electrode */<br>
RE = 300;<br>
<br>
Point(RE+1) = { R1, R1/2, R1/2, q};<br>
Point(RE+2) = { R1, -R1/2, R1/2, q};<br>
Point(RE+3) = { R1, -R1/2, -R1/2, q};<br>
Point(RE+4) = { R1, R1/2, -R1/2, q};<br>
Line(RE+13) = {RE+4, RE+1};<br>
Line(RE+14) = {RE+1, RE+2};<br>
Line(RE+15) = {RE+2, RE+3};<br>
Line(RE+16) = {RE+3, RE+4};<br>
Line Loop(RE+7) = {RE+16, RE+13, RE+14, RE+15};<br>
Plane Surface(7) = {RE+7};<br>
Physical Surface(10000) = {7};<br>
<br>
/* Left electrode */<br>
<br>
RE = 400;<br>
<br>
Point(RE+1) = { -R1, R1/2, R1/2, q};<br>
Point(RE+2) = { -R1, -R1/2, R1/2, q};<br>
Point(RE+3) = { -R1, -R1/2, -R1/2, q};<br>
Point(RE+4) = { -R1, R1/2, -R1/2, q};<br>
Line(RE+13) = {RE+4, RE+1};<br>
Line(RE+14) = {RE+1, RE+2};<br>
Line(RE+15) = {RE+2, RE+3};<br>
Line(RE+16) = {RE+3, RE+4};<br>
Line Loop(RE+7) = {RE+16, RE+13, RE+14, RE+15};<br>
Plane Surface(8) = {RE+7};<br>
Physical Surface(10001) = {8};<br>
<br>
<br>
// create plate electrodes<br>
/* back electrode */<br>
<br>
RE = 500;<br>
<br>
Point(RE+1) = { R1/2, -R1, R1/2, q};<br>
Point(RE+2) = { -R1/2, -R1, R1/2, q};<br>
Point(RE+3) = { -R1/2, -R1,-R1/2, q};<br>
Point(RE+4) = { R1/2, -R1, -R1/2, q};<br>
Line(RE+13) = {RE+4, RE+1};<br>
Line(RE+14) = {RE+1, RE+2};<br>
Line(RE+15) = {RE+2, RE+3};<br>
Line(RE+16) = {RE+3, RE+4};<br>
Line Loop(RE+7) = {RE+16, RE+13, RE+14, RE+15};<br>
Plane Surface(9) = {RE+7};<br>
Physical Surface(10002) = {9};<br>
<br>
<br>
/* front electrode */<br>
<br>
RE = 600;<br>
<br>
Point(RE+1) = { R1/2, R1, R1/2, q};<br>
Point(RE+2) = { -R1/2, R1, R1/2, q};<br>
Point(RE+3) = { -R1/2, R1,-R1/2, q};<br>
Point(RE+4) = { R1/2, R1, -R1/2, q};<br>
Line(RE+13) = {RE+4, RE+1};<br>
Line(RE+14) = {RE+1, RE+2};<br>
Line(RE+15) = {RE+2, RE+3};<br>
Line(RE+16) = {RE+3, RE+4};<br>
Line Loop(RE+7) = {RE+16, RE+13, RE+14, RE+15};<br>
Plane Surface(10) = {RE+7};<br>
Physical Surface(10003) = {10};<br>
<br>
/* top electrode */<br>
<br>
RE = 700;<br>
<br>
Point(RE+1) = { R1/2,R1/2,R1, q};<br>
Point(RE+2) = { -R1/2,R1/2, R1,q};<br>
Point(RE+3) = { -R1/2, -R1/2,R1, q};<br>
Point(RE+4) = { R1/2, -R1/2,R1, q};<br>
Line(RE+13) = {RE+4, RE+1};<br>
Line(RE+14) = {RE+1, RE+2};<br>
Line(RE+15) = {RE+2, RE+3};<br>
Line(RE+16) = {RE+3, RE+4};<br>
Line Loop(RE+7) = {RE+16, RE+13, RE+14, RE+15};<br>
Plane Surface(11) = {RE+7};<br>
Physical Surface(10004) = {11};<br>
<br>
<br>
/* bottom electrode */<br>
<br>
RE = 800;<br>
<br>
Point(RE+1) = { R1/2,R1/2,-R1, q};<br>
Point(RE+2) = { -R1/2,R1/2, -R1,q};<br>
Point(RE+3) = { -R1/2, -R1/2,-R1, q};<br>
Point(RE+4) = { R1/2, -R1/2,-R1, q};<br>
Line(RE+13) = {RE+4, RE+1};<br>
Line(RE+14) = {RE+1, RE+2};<br>
Line(RE+15) = {RE+2, RE+3};<br>
Line(RE+16) = {RE+3, RE+4};<br>
Line Loop(RE+7) = {RE+16, RE+13, RE+14, RE+15};<br>
Plane Surface(12) = {RE+7};<br>
Physical Surface(10005) = {12};<br>
<br>
Surface{7,8,9,10,11,12} In Volume{1};<br>
Physical Volume(100) = {1};<br>
```<br>
<br>
Somehow the command "Surface{} In Volume" seems to not care for what I want. Please let me know<br>
how can I get a cube with each face having a rectangular surface with its own physical number.<br>
<br>
<br>
<br>
<br>
<br>
________________________________<br>
Best Regards,<br>
<br>
Sathyanarayan Rao, PhD student<br>
Earth and Life Institute/Environmental Sciences (ELI-e)<br>
Universit? catholique de Louvain<br>
c.037, Croix du Sud 2, Louvain-la-Neuve, Belgium<br>
Phone: 010473827 ( intercom 73827)<br>
<br>
________________________________<br>
<br>
_______________________________________________<br>
gmsh mailing list<br>
gmsh@onelab.info<<a href="mailto:gmsh@onelab.info">mailto:gmsh@onelab.info</a>><br>
<a href="http://onelab.info/mailman/listinfo/gmsh">http://onelab.info/mailman/listinfo/gmsh</a><br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://onelab.info/pipermail/gmsh/attachments/20180209/e0eb17ae/attachment.html">http://onelab.info/pipermail/gmsh/attachments/20180209/e0eb17ae/attachment.html</a>><br>
<br>
------------------------------<br>
<br>
Subject: Digest Footer<br>
<br>
_______________________________________________<br>
gmsh mailing list<br>
gmsh@onelab.info<br>
<a href="http://onelab.info/mailman/listinfo/gmsh">http://onelab.info/mailman/listinfo/gmsh</a><br>
<br>
<br>
------------------------------<br>
<br>
End of gmsh Digest, Vol 181, Issue 20<br>
*************************************<br>
</div>
</span></font>
</body>
</html>