[Gmsh] Sphere problem

mhengs at gmx.at mhengs at gmx.at
Thu Oct 28 14:16:59 CEST 2004


Hi,

i am working with 2D structured meshed spheres. Everything works fine, but
i have problem with the small elements on the top of the spheres. My solver
cannot work with this elements.

Maybe you have a solution how to cope with this problem.

Markus

Here is an exmaple of 2x2x2 sphereblock:

/********************************************************************* 
 *
 *  
 *  source file: example.igs
 *  translater version: jIGStoGEO 0.2  
 *  
 *  output files:
 *  example.geo 
 *  example.msh
 *
 *********************************************************************/

 
Function Sphere 
 
  p1 = newp; Point(p1) = {x,  y,  z,  1.0} ;
  p2 = newp; Point(p2) = {x+r,y,  z,   1.0} ;
  p3 = newp; Point(p3) = {x,  y+r,z,   1.0} ;
  p4 = newp; Point(p4) = {x,  y,  z+r, 1.0} ;
  p5 = newp; Point(p5) = {x-r,y,  z,   1.0} ;
  p6 = newp; Point(p6) = {x,  y-r,z,   1.0} ;
  p7 = newp; Point(p7) = {x,  y,  z-r, 1.0} ;
 
  c1 = newreg; Circle(c1) = {p2,p1,p7};
  c2 = newreg; Circle(c2) = {p7,p1,p5};
  c3 = newreg; Circle(c3) = {p5,p1,p4};
  c4 = newreg; Circle(c4) = {p4,p1,p2};
  c5 = newreg; Circle(c5) = {p2,p1,p3};
  c6 = newreg; Circle(c6) = {p3,p1,p5};
  c7 = newreg; Circle(c7) = {p5,p1,p6};
  c8 = newreg; Circle(c8) = {p6,p1,p2};
  c9 = newreg; Circle(c9) = {p7,p1,p3};
  c10 = newreg; Circle(c10) = {p3,p1,p4};
  c11 = newreg; Circle(c11) = {p4,p1,p6};
  c12 = newreg; Circle(c12) = {p6,p1,p7};
 
 
  l1 = newreg; Line Loop(l1) = {c5,c10,c4};
  l2 = newreg; Line Loop(l2) = {c9,-c5,c1};
  l3 = newreg; Line Loop(l3) = {c12,-c8,-c1};
  l4 = newreg; Line Loop(l4) = {c8,-c4,c11};
  l5 = newreg; Line Loop(l5) = {-c10,c6,c3};
  l6 = newreg; Line Loop(l6) = {-c11,-c3,c7};
  l7 = newreg; Line Loop(l7) = {-c2,-c7,-c12};
  l8 = newreg; Line Loop(l8) = {-c6,-c9,c2};
 
 
  s1 = newreg; Ruled Surface(newreg) = {l1};
  s2 = newreg; Ruled Surface(newreg) = {l2};
  s3 = newreg; Ruled Surface(newreg) = {l3};
  s4 = newreg; Ruled Surface(newreg) = {l4};
  s5 = newreg; Ruled Surface(newreg) = {l5};
  s6 = newreg; Ruled Surface(newreg) = {l6};
  s7 = newreg; Ruled Surface(newreg) = {l7};
  s8 = newreg; Ruled Surface(newreg) = {l8};
 
 Transfinite Line {c5,c10,c4} = 15 Using Progression 1.0;
 Transfinite Line {c9,-c5,c1} = 15 Using Progression 1.0;
 Transfinite Line {c12,-c8,-c1} = 15 Using Progression 1.0;
 Transfinite Line {c8,-c4,c11} = 15 Using Progression 1.0;
 Transfinite Line {-c10,c6,c3} = 15 Using Progression 1.0;
 Transfinite Line {-c11,-c3,c7} = 15 Using Progression 1.0;
 Transfinite Line {-c2,-c7,-c12} = 15 Using Progression 1.0;
 Transfinite Line {-c6,-c9,c2} = 15 Using Progression 1.0;
 
 Transfinite Surface {s1} = {p4,p3,p2};
 Transfinite Surface {s2} = {p7,p3,p2};
 Transfinite Surface {s3} = {p7,p6,p2};
 Transfinite Surface {s4} = {p4,p6,p2};
 Transfinite Surface {s5} = {p4,p5,p3};
 Transfinite Surface {s6} = {p4,p6,p5};
 Transfinite Surface {s7} = {p7,p6,p5};
 Transfinite Surface {s8} = {p7,p3,p5};
 
 Recombine Surface {s1};
 Recombine Surface {s2};
 Recombine Surface {s3};
 Recombine Surface {s4};
 Recombine Surface {s5};
 Recombine Surface {s6};
 Recombine Surface {s7};
 Recombine Surface {s8};
 
Return
 
x = 0;
y = 0;
z = 0;
r = 0.5;
 
 
For w In {1:2}
For u In {1:2}
For t In {1:2}
 
 
 x += 1.0;
 
  Call Sphere ;
 
 
EndFor
x = 0;
y += 1.0;
 
EndFor
z += 1.0;
x = 0;
y = 0;
 
 
EndFor
 

-- 
NEU +++ DSL Komplett von GMX +++ http://www.gmx.net/de/go/dsl
GMX DSL-Netzanschluss + Tarif zum supergünstigen Komplett-Preis!