# Making a cylindrical hole in a cubic domain (wellbore) in GMsh

 Register Blogs Members List Search Today's Posts Mark Forums Read

 July 29, 2015, 18:36 Making a cylindrical hole in a cubic domain (wellbore) in GMsh #1 New Member   Igor Shovkun Join Date: Jul 2015 Posts: 2 Rep Power: 0 Hello everyone! I am trying to create a 3D mesh to simulate geomechanics in a cubic domain with a deviated wellbore. I have build all the geometry, surface (2D) mesh generation works fine. When I try to create a 3D mesh for a volume, GMsh fills the hollow space as well. Could anybody help me with that? There follows my code: ft = 0.3048; xSize = 20*ft; ySize = 20*ft; zSize = 20*ft; esize = 1*ft; deg = 2*Pi/360; alpha = 85*deg; // angle between hor plane and well beta = 0*deg; // angle between x and well r = 2*ft; // well radius // parameters of the horizontal cross-section e = Cos(alpha); a = r/Sin(alpha); b = a*Sqrt(1-e^2); // creating ellipse Point(1) = {a*Cos(beta),a*Sin(beta),0,esize}; // first point of the arc Point(2) = {0,0,0,esize}; // center Point(3) = {a/2*Cos(beta),a/2*Sin(beta),0,esize}; // point in major axis Point(4) = {b*Sin(beta),b*Cos(beta),0,esize}; // end point of the arc Point(5) = {-a*Cos(beta),-a*Sin(beta),0,esize}; // first point of the arc Point(6) = {-b*Sin(beta),-b*Cos(beta),0,esize}; // end point of the arc Ellipse(1) = {1,2,3,4}; Ellipse(2) = {4,2,3,5}; Ellipse(3) = {5,2,3,6}; Ellipse(4) = {6,2,3,1}; // Line Loop(5) = {1,2,3,4}; Line Loop(5) = {-4,-3,-2,-1}; out[] = Extrude{ zSize*Cos(alpha)*Cos(beta), zSize*Cos(alpha)*Sin(beta), zSize} { // Surface{6}; // Loop{5}; Line{-4,-3,-2,-1}; }; // // Points // Bottom Point(31) = {xSize/2,ySize/2,0,esize}; Point(32) = {-xSize/2,ySize/2,0,esize}; Point(33) = {-xSize/2,-ySize/2,0,esize}; Point(34) = {xSize/2,-ySize/2,0,esize}; // Top Point(35) = {xSize/2,ySize/2,zSize,esize}; Point(36) = {-xSize/2,ySize/2,zSize,esize}; Point(37) = {-xSize/2,-ySize/2,zSize,esize}; Point(38) = {xSize/2,-ySize/2,zSize,esize}; // Lines // Bottom Line(31) = {31,32}; Line(32) = {32,33}; Line(33) = {33,34}; Line(34) = {34,31}; Line Loop(35) = {31,32,33,34}; surf2 = news; Plane Surface(surf2) = {35,5}; // bottom surf // // Top Line(36) = {35,36}; Line(37) = {36,37}; Line(38) = {37,38}; Line(39) = {38,35}; Line Loop(40) = {36,37,38,39}; Line Loop(41) = {out[0],out[5],out[10],out[15]}; // upper ellipse // Plane Surface(666) = {41}; surf3 = news; Plane Surface(surf3) = {40,41}; //top surf // Front Line(42) = {34,38}; Line(43) = {38,37}; Line(44) = {37,33}; Line(45) = {33,34}; Line Loop(46) = {42,43,44,45}; surf4 = news; Plane Surface(surf4) = {46}; // Left Line(47) = {33,37}; Line(48) = {37,36}; Line(49) = {36,32}; Line(50) = {32,33}; Line Loop(51) = {47,48,49,50}; surf5 = news; Plane Surface(surf5) = {51}; // Back Line(52) = {32,36}; Line(53) = {36,35}; Line(54) = {35,31}; Line(55) = {31,32}; Line Loop(56) = {52,53,54,55}; surf6 = news; Plane Surface(surf6) = {56}; // Right Line(57) = {31,35}; Line(58) = {35,38}; Line(59) = {38,34}; Line(60) = {34,31}; Line Loop(61) = {57,58,59,60}; surf7 = news; Plane Surface(surf7) = {61}; surf8 = news; // Surface Loop(surf8) = {surf2,surf5,surf4,surf3,surf6,surf7,-9,-13,-17,-21} ; Surface Loop(surf8) = {surf2,-9,-13,-17,-21,surf3,surf5,surf4,surf7,surf6} ; Volume(1) = {surf8} ; //

 July 30, 2015, 10:44 Solution #2 New Member   Igor Shovkun Join Date: Jul 2015 Posts: 2 Rep Power: 0 In case anyone needs the solution, I got this! Just add Coherence; As the last line, and It will work as a charm.

 Tags cylindrical hole, wellbore

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post gocarts OpenFOAM Announcements from Other Sources 36 March 29, 2017 08:05 GGerber OpenFOAM Installation 54 April 24, 2015 16:02 allenzhao OpenFOAM Installation 127 January 30, 2009 20:08 hjasak OpenFOAM Post-Processing 69 April 24, 2008 01:24 jens_klostermann OpenFOAM Bugs 11 June 28, 2007 17:51

All times are GMT -4. The time now is 10:55.