ramc2struc.geo

//Gmsh script for ramc2struc.geo
//+
//Wall points
Point(1) = {0, 0, 0, 1.0};
Point(2) = {0.12856, 0.15052, 0, 1.0};
Point(3) = {1.2954, 0.33533, 0, 1.0};
//Domain points
Point(4) = {-0.3048, 0, 0, 1.0};
Point(5) = {-0.053514, 0.49318, 0, 1.0};
Point(6) = {1.2954, 1.4732, 0, 1.0};
//Circular arc centers
Point(7) = {0.1524, 0, 0, 1.0};
Point(8) = {0.3048, 0, 0, 1.0};
//Lines
Circle(1) = {1, 7, 2};
Circle(2) = {4, 8, 5};
Line(3) = {1, 4};
Line(4) = {5, 6};
Line(5) = {3, 6};
Line(6) = {2, 3};
//+
Line(7) = {2, 5};
//+
//Loops
Curve Loop(1) ={1, 7, -2, -3};
Curve Loop(2) ={6, 5, -4, -7};
//+
//Surfaces
Plane Surface(1) = {1};
Plane Surface(2) = {2};
// Mesh Parameters
n_nose = 400;
r_nose = 1;
n_cone = 800;
r_cone = 0.4;
n_w_d = 1600;
r_w_d = 1.008;
//+
Transfinite Curve {1, 2} = n_nose Using Progression r_nose;
Transfinite Curve {6, 4} = n_cone Using Bump r_cone;
Transfinite Curve {3, 7, 5} = n_w_d Using Progression r_w_d;
//+
Transfinite Surface {1};
Transfinite Surface {2};
//+
Recombine Surface {1, 2};
//+
Extrude {0, 0, 0.2} {
	Surface{1}; Surface{2}; Layers {1}; Recombine;
}
//+
//Boundary Conditions
Physical Volume("FLUID", 52) = {1, 2};
Physical Surface("INLET", 53) = {24};
Physical Surface("FARFIELD", 54) = {38};
Physical Surface("WALL", 55) = {16, 46};
Physical Surface("OUTLET", 56) = {42};
Physical Surface("SYMMETRY", 57) = {28, 1, 2, 29, 51};
//+