Homec4science
Diffusion ConSol (master) Browse

ConSol/ConSol/Examples/HybridTest/ShockVortex/shockVortex8QUAD.geo
a31395923f2dmaster

Last Change
Raw File

shockVortex8QUAD.geo
View Options

// Variables
boxdim = 1;
N = 2*28*2+1;
N0 = 2;
gridsize = boxdim/(N-1);
// Points
Point(1) = {0,0,0, gridsize};
Point(2) = {boxdim,0,0, gridsize*3};
Point(3) = {boxdim, boxdim,0, gridsize*3};
Point(4) = {0,boxdim,0, gridsize};
Point(5) = {boxdim/3-N0*gridsize,0,0, gridsize};
Point(6) = {boxdim/3,0,0, gridsize};
Point(7) = {boxdim*2/3,0,0, gridsize};
Point(8) = {boxdim*2/3+N0*gridsize,0,0, gridsize};
Point(9) = {boxdim, boxdim/3-N0*gridsize,0, gridsize};
Point(10) = {boxdim, boxdim/3,0, gridsize};
Point(11) = {boxdim, boxdim*2/3,0, gridsize};
Point(12) = {boxdim, boxdim*2/3+N0*gridsize,0, gridsize};
Point(13) = {boxdim*2/3+N0*gridsize,boxdim,0, gridsize};
Point(14) = {boxdim*2/3,boxdim,0, gridsize};
Point(15) = {boxdim/3,boxdim,0, gridsize};
Point(16) = {boxdim/3-N0*gridsize,boxdim,0, gridsize};
Point(17) = {0, boxdim*2/3+N0*gridsize,0, gridsize};
Point(18) = {0, boxdim*2/3,0, gridsize};
Point(19) = {0, boxdim/3,0, gridsize};
Point(20) = {0, boxdim/3-N0*gridsize,0, gridsize};
Point(21) = {boxdim/3-N0*gridsize,boxdim/3-N0*gridsize,0, gridsize};
Point(22) = {boxdim/3,boxdim/3-N0*gridsize,0, gridsize};
Point(23) = {boxdim/3,boxdim/3,0, gridsize};
Point(24) = {boxdim/3-N0*gridsize,boxdim/3,0, gridsize};
Point(25) = {boxdim*2/3,boxdim/3-N0*gridsize,0, gridsize};
Point(26) = {boxdim*2/3+N0*gridsize,boxdim/3-N0*gridsize,0, gridsize};
Point(27) = {boxdim*2/3+N0*gridsize,boxdim/3,0, gridsize};
Point(28) = {boxdim*2/3,boxdim/3,0, gridsize};
Point(29) = {boxdim*2/3,boxdim*2/3,0, gridsize};
Point(30) = {boxdim*2/3+N0*gridsize,boxdim*2/3,0, gridsize};
Point(31) = {boxdim*2/3+N0*gridsize,boxdim*2/3+N0*gridsize,0, gridsize};
Point(32) = {boxdim*2/3,boxdim*2/3+N0*gridsize,0, gridsize};
Point(33) = {boxdim/3-N0*gridsize,boxdim*2/3,0, gridsize};
Point(34) = {boxdim/3,boxdim*2/3,0, gridsize};
Point(35) = {boxdim/3,boxdim*2/3+N0*gridsize,0, gridsize};
Point(36) = {boxdim/3-N0*gridsize,boxdim*2/3+N0*gridsize,0, gridsize};
Line(1) = {1,5};
Line(2) = {5,6};
Line(3) = {6,7};
Line(4) = {7,8};
Line(5) = {8,2};
Line(6) = {2,9};
Line(7) = {9,10};
Line(8) = {10,11};
Line(9) = {11,12};
Line(10) = {12,3};
Line(11) = {13,3};
Line(12) = {14,13};
Line(13) = {15,14};
Line(14) = {16,15};
Line(15) = {4,16};
Line(16) = {17,4};
Line(17) = {18,17};
Line(18) = {19,18};
Line(19) = {20,19};
Line(20) = {1,20};
Line(21) = {20,21};
Line(22) = {21,22};
Line(23) = {22,25};
Line(24) = {25,26};
Line(25) = {26,9};
Line(26) = {19,24};
Line(27) = {24,23};
Line(28) = {23,28};
Line(29) = {28,27};
Line(30) = {27,10};
Line(31) = {18,33};
Line(32) = {33,34};
Line(33) = {34,29};
Line(34) = {29,30};
Line(35) = {30,11};
Line(36) = {17,36};
Line(37) = {36,35};
Line(38) = {35,32};
Line(39) = {32,31};
Line(40) = {31,12};
Line(41) = {5,21};
Line(42) = {21,24};
Line(43) = {24,33};
Line(44) = {33,36};
Line(45) = {36,16};
Line(46) = {6,22};
Line(47) = {22,23};
Line(48) = {23,34};
Line(49) = {34,35};
Line(50) = {35,15};
Line(51) = {7,25};
Line(52) = {25,28};
Line(53) = {28,29};
Line(54) = {29,32};
Line(55) = {32,14};
Line(56) = {8,26};
Line(57) = {26,27};
Line(58) = {27,30};
Line(59) = {30,31};
Line(60) = {31,13};
Line Loop(1) = {1,41,-21,-20};
Line Loop(2) = {2,46,-22,-41};
Line Loop(3) = {3,51,-23,-46};
Line Loop(4) = {4,56,-24,-51};
Line Loop(5) = {5,6,-25,-56};
Line Loop(6) = {21,42,-26,-19};
Line Loop(7) = {22,47,-27,-42};
Line Loop(8) = {23,52,-28,-47};
Line Loop(9) = {24,57,-29,-52};
Line Loop(10) = {25,7,-30,-57};
Line Loop(11) = {26,43,-31,-18};
Line Loop(12) = {27,48,-32,-43};
Line Loop(13) = {28,53,-33,-48};
Line Loop(14) = {29,58,-34,-53};
Line Loop(15) = {30,8,-35,-58};
Line Loop(16) = {31,44,-36,-17};
Line Loop(17) = {32,49,-37,-44};
Line Loop(18) = {33,54,-38,-49};
Line Loop(19) = {34,59,-39,-54};
Line Loop(20) = {35,9,-40,-59};
Line Loop(21) = {36,45,-15,-16};
Line Loop(22) = {37,50,-14,-45};
Line Loop(23) = {38,55,-13,-50};
Line Loop(24) = {39,60,-12,-55};
Line Loop(25) = {40,10,-11,-60};
Plane Surface(1) = {1};
Plane Surface(2) = {2};
Plane Surface(3) = {3};
Plane Surface(4) = {4};
Plane Surface(5) = {5};
Plane Surface(6) = {6};
Plane Surface(7) = {7};
Plane Surface(8) = {8};
Plane Surface(9) = {9};
Plane Surface(10) = {10};
Plane Surface(11) = {11};
Plane Surface(12) = {12};
Plane Surface(13) = {13};
Plane Surface(14) = {14};
Plane Surface(15) = {15};
Plane Surface(16) = {16};
Plane Surface(17) = {17};
Plane Surface(18) = {18};
Plane Surface(19) = {19};
Plane Surface(20) = {20};
Plane Surface(21) = {21};
Plane Surface(22) = {22};
Plane Surface(23) = {23};
Plane Surface(24) = {24};
Plane Surface(25) = {25};
Physical Surface("tri_0") = {13};
Physical Surface("qdr_0") = {1};
Physical Surface("qdr_1") = {3};
Physical Surface("qdr_2") = {5};
Physical Surface("qdr_3") = {11};
Physical Surface("qdr_4") = {15};
Physical Surface("qdr_5") = {21};
Physical Surface("qdr_6") = {23};
Physical Surface("qdr_7") = {25};
Physical Surface("mid_0") = {8};
Physical Surface("mid_1") = {12};
Physical Surface("mid_2") = {14};
Physical Surface("mid_3") = {18};
Physical Surface("q2q_0") = {2};
Physical Surface("q2q_1") = {4};
Physical Surface("q2q_2") = {6};
Physical Surface("q2q_3") = {10};
Physical Surface("q2q_4") = {16};
Physical Surface("q2q_5") = {20};
Physical Surface("q2q_6") = {22};
Physical Surface("q2q_7") = {24};
Physical Surface("rst_1") = {7};
Physical Surface("rst_2") = {9};
Physical Surface("rst_3") = {17};
Physical Surface("rst_4") = {19};
Physical Line("inflow") = {20,19,18,17,16};
Physical Line("outflow") = {6,7,8,9,10};
Physical Line("slipWall") = {1,2,3,4,5,15,14,13,12,11};
Physical Line("q2i") = {23,58,38,43,22,47,27,42,24,57,29,52,34,59,39,54,32,49,37,44};
Physical Line("t2i") = {28,53,33,48};
Physical Line("q2q") = {41,46,51,56,45,50,55,60,21,26,31,36,25,30,35,40};
Transfinite Curve{1,21,26,31,36,15,5,25,30,35,40,11,20,41,46,51,56,6,16,45,50,55,60,10} = (N-1)/3+1 - N0;
Transfinite Curve{2,22,27,32,37,14,4,24,29,34,39,12,19,42,47,52,57,7,17,44,49,54,59,9} = N0+1;
Transfinite Curve{3,23,28,33,38,13,18,43,48,53,58,8} = ((N-1)/3) + 1;
Transfinite Surface{1,2,3,4,5,6,7,8,9,10,11,12,14,15,16,17,18,19,20,21,22,23,24,25};
Recombine Surface{1,2,3,4,5,6,7,8,9,10,11,12,14,15,16,17,18,19,20,21,22,23,24,25};
c4science · Help