[PhilipA]

Hello all,

I have been struggling to run an airfoil, made using NACA 63-212 coordinates, through my first SU2 simulation. I am using the attached config file. (edited from the turbulent NACA 0012 tutorial config file. I made sure to create the mesh following the description on the tutorial page.)

I am coming across a problem where my physical group markers are not showing up in my mesh file. This is the error:

Error in "bool CSU2ASCIIMeshReaderFVM::ReadMetadata(const bool, CConfig *)":
-------------------------------------------------------------------------
Could not find the keyword "NMARK=".
Check the SU2 ASCII file format.
------------------------------ Error Exit -------------------------------

This error occurs when I create the physical groups using a Surface. However, if I use a Curve I get the marker keywords in the mesh file, but it creates a mismatch in NPOIN and number of points in the mesh file (and a segmentation fault if I manually alter the number of points in the mesh file)
-This marker error occurs with save all elements on or off

I tried the solution of this forum post Gmsh does not ouput SU2 markers from physical surfaces, but it did not solve my issue. I have considered editing the airfoil and farfield markers manually into the mesh file, but I struggled to figure out where I would put the markers, as it relies on me knowing exactly which mesh elements are tied to each marker.

I attached a picture of the mesh. The mesh file, geometry script and config file are in the zip file "SU2 error files"

Any help would be much appreciated!

[CleverBoy]

I am not sure but it could be related with defining airfoil and boundaries as a pyhsical surface. Try not to select airfoil as a surface and leave it as it is. Select the boundaries using "physical curves", not surfaces and then try again.

[bigfootedrockmidget]

I reran your geo file and I think the only thing you forgot was that you also have to add the (interior) surface as a physical surface. Not sure if you want to have the interior mesh of the airfoil included as well. In my case, I excluded it, so the last couple of lines of the geo file are :

Plane Surface(2) = {1,2}; 
Coherence; 
//+ 
Physical Surface("interior", 55) = {1}; 
//+ 
Physical Curve("free", 57) = {51, 52, 54, 53}; 
//+ 
Physical Curve("foil", 58) = {15, 16, 28, 27, 26, 25, 24, 23, 22, 21, 20, 19, 18, 17, 29, 14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 42, 50, 49, 48, 47, 46, 45, 44, 43, 3, 41, 40, 39, 38, 37, 36, 35, 34, 33, 32, 31, 30, 2, 1}

[PhilipA]

Hey @CleverBoy and @bigfootedrockmidget. Thank you for taking the time to post your thoughts on my problem.

I have tried to implement both of your proposed solutions but am still running into issues with the mesh file.

I will update you guys if I come across any new information, in the event that I realise I implemented your solutions wrong or if I come across a different solution.

[bigfootedrockmidget]

Just remove your physical markers from the geo file and re-define the foil and the free stream and then define the interior as a surface.

Do not save all elements, this might introduce points that are not used by the mesh and su2 does not like that.This is the NPOIN error that you saw. Below is the mesh that worked for me. Note that it's better to switch to a smooth surface curve.

//+
Point(1) = {1, 0, 0, 1.0};
//+
Point(2) = {0.95012, 0.005660, 0, 1.0};
//+
Point(3) = {0.900250, 0.012240, 0, 1.0};
//+
Point(4) = {0.850350, 0.019470, 0, 1.0};
//+
Point(5) = {0.800420, 0.026980, 0, 1.0};
//+
Point(6) = {0.750430, 0.034510, 0, 1.0};
//+
Point(7) = {0.700430, 0.041820, 0, 1.0};
//+
Point(8) = {0.650380, 0.048700, 0, 1.0};
//+
Point(9) = {0.600290, 0.054910, 0, 1.0};
//+
Point(10) = {0.550160, 0.060300, 0, 1.0};
//+
Point(11) = {0.50, 0.064730, 0, 1.0};
//+
Point(12) = {0.449820, 0.067990, 0, 1.0};
//+
Point(13) = {0.399620, 0.069910, 0, 1.0};
//+
Point(14) = {0.349410, 0.0703, 0, 1.0};
//+
Point(15) = {0.2992, 0.06901, 0, 1.0};
//+
Point(16) = {0.249, 0.06606, 0, 1.0};
//+
Point(17) = {0.19882, 0.06137, 0, 1.0};
//+
Point(18) = {0.14868, 0.0547, 0, 1.0};
//+
Point(19) = {0.098590, 0.04554, 0, 1.0};
//+
Point(20) = {0.073580, 0.039630, 0, 1.0};
//+
Point(21) = {0.048630, 0.03238, 0, 1.0};
//+
Point(22) = {0.02378, 0.02284, 0, 1.0};
//+
Point(23) = {0.01145, 0.01622, 0, 1.0};
//+
Point(24) = {0.00657, 0.0126, 0, 1.0};
//+
Point(25) = {0.00417, 0.01032, 0, 1.0};
//+
Point(26) = {0, 0, 0, 1.0};
//+
Point(27) = {0.00538, -0.00932, 0, 1.0};
//+
Point(28) = {0.00843, -0.0112, 0, 1.0};
//+
Point(29) = {0.01355, -0.01408, 0, 1.0};
//+
Point(30) = {0.02622, -0.01912, 0, 1.0};
//+
Point(31) = {0.05137, -0.02606, 0, 1.0};
//+
Point(32) = {0.07642, -0.03115, 0, 1.0};
//+
Point(33) = {0.10141, -0.0352, 0, 1.0};
//+
Point(34) = {0.15132, -0.04124, 0, 1.0};
//+
Point(35) = {0.20118, -0.04545, 0, 1.0};
//+
Point(36) = {0.251, -0.04816, 0, 1.0};
//+
Point(37) = {0.3008, -0.04957, 0, 1.0};
//+
Point(38) = {0.35059, -0.0497, 0, 1.0};
//+
Point(39) = {0.40038, -0.04849, 0, 1.0};
//+
Point(40) = {0.45018, -0.04609, 0, 1.0};
//+
Point(41) = {0.5, -0.04267, 0, 1.0};
//+
Point(42) = {0.54848, -0.0384, 0, 1.0};
//+
Point(43) = {0.59971, -0.03349, 0, 1.0};
//+
Point(44) = {0.64962, -0.0281, 0, 1.0};
//+
Point(45) = {0.69957, -0.02238, 0, 1.0};
//+
Point(46) = {0.74955, -0.01661, 0, 1.0};
//+
Point(47) = {0.79958, -0.01106, 0, 1.0};
//+
Point(48) = {0.84965, -0.00601, 0, 1.0};
//+
Point(49) = {0.89975, -0.0019, 0, 1.0};
//+
Point(50) = {0.94988, 0.00066, 0, 1.0};
//+
Line(1) = {1, 2};
//+
Line(2) = {2, 3};
//+
Line(3) = {3, 4};
//+
Line(4) = {4, 5};
//+
Line(5) = {5, 6};
//+
Line(6) = {6, 7};
//+
Line(7) = {7, 8};
//+
Line(8) = {8, 9};
//+
Line(9) = {9, 10};
//+
Line(10) = {10, 11};
//+
Line(11) = {11, 12};
//+
Line(12) = {12, 13};
//+
Line(13) = {13, 14};
//+
Line(14) = {14, 15};
//+
Line(15) = {15, 16};
//+
Line(16) = {16, 17};
//+
Line(17) = {17, 18};
//+
Line(18) = {18, 19};
//+
Line(19) = {19, 20};
//+
Line(20) = {20, 21};
//+
Line(21) = {21, 22};
//+
Line(22) = {22, 23};
//+
Line(23) = {23, 24};
//+
Line(24) = {24, 25};
//+
Line(25) = {25, 26};
//+
Line(26) = {26, 27};
//+
Line(27) = {27, 28};
//+
Line(28) = {28, 29};
//+
Line(29) = {29, 30};
//+
Line(30) = {30, 31};
//+
Line(31) = {31, 32};
//+
Line(32) = {32, 33};
//+
Line(33) = {33, 34};
//+
Line(34) = {34, 35};
//+
Line(35) = {35, 36};
//+
Line(36) = {36, 37};
//+
Line(37) = {37, 38};
//+
Line(38) = {38, 39};
//+
Line(39) = {39, 40};
//+
Line(40) = {40, 41};
//+
Line(41) = {41, 42};
//+
Line(42) = {42, 43};
//+
Line(43) = {43, 44};
//+
Line(44) = {44, 45};
//+
Line(45) = {45, 46};
//+
Line(46) = {46, 47};
//+
Line(47) = {47, 48};
//+
Line(48) = {48, 49};
//+
Line(49) = {49, 50};
//+
Line(50) = {50, 1};
//+
Point(51) = {2.3, 4, 0, 1.0};
//+
Point(52) = {2.3, -4, 0, 1.0};
//+
Point(53) = {5.2, 4, 0, 1.0};
//+
Point(54) = {5.2, -4, 0, 1.0};
//+
Point(55) = {2.6, 0, 0, 1.0};
//+
Circle(51) = {51, 55, 52};
//+
Line(52) = {51, 53};
//+
Line(53) = {53, 54};
//+
Line(54) = {54, 52};
//
Curve Loop(1) = {54, -51, 52, 53};
//+
Curve Loop(2) = {13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12};
//+
Plane Surface(1) = {1, 2};
//+
Physical Curve("freestream", 55) = {52, 53, 54, 51};
//+
Physical Curve("airfoil", 56) = {10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 50, 49, 48, 47, 46, 45, 44, 43, 42, 41, 40, 39, 38, 37, 36, 35, 34, 33, 32, 31, 30, 29, 28, 27, 26, 25, 24, 23, 22, 21, 20, 19, 18, 17, 16, 15, 14, 13, 12, 11};
//+
Physical Surface("interior", 57) = {1};

[PhilipA]

@bigfootedrockmidget I just copied your solution and this is the error I get. The same thing occurs if I put save all elements on, but with more points.

------------------- Geometry Preprocessing ( Zone 0 ) -------------------
Two dimensional problem.
76 grid points.
0 volume elements.
2 surface markers.
50 boundary elements in index 0 (Marker = airfoil).
26 boundary elements in index 1 (Marker = farfield).


Error in "void CPhysicalGeometry::DistributeColoring(const CConfig *, CGeometry *)": 
-------------------------------------------------------------------------
Mismatch between NPOIN and number of points listed in mesh file.
Please check the mesh file for correctness.

------------------------------ Error Exit -------------------------------

Method:

• I am exporting the mesh file by first clicking on 2D in the mesh category then exporting into a .su2 file format.
  
• Save all elements is turned off.
  
• The config file is run from the SU2_CFD file location

Extra information

• I am using version 7.4.0 "Blackbird"

• The number of points I input for the geometry is 55, the mesh generates 76. The extra points seem correspond to points where lines in the mesh meet. This makes sense to me, but I am not sure if it is causing a problem.
  • I have tried to delete the extra points manually but it leads to segmentation fault 11

[bigfootedrockmidget]

ok, I see it, it looks like there is a difference in the line loops. The actual additional mesh points and connectivities were not saved. If I use the .geo file below (note: also updated it in the previous post), then the points and connectivities are saved. Gmsh then says:

Info : Writing 586 elements and 331 nodes

//+
Point(1) = {1, 0, 0, 1.0};
//+
Point(2) = {0.95012, 0.005660, 0, 1.0};
//+
Point(3) = {0.900250, 0.012240, 0, 1.0};
//+
Point(4) = {0.850350, 0.019470, 0, 1.0};
//+
Point(5) = {0.800420, 0.026980, 0, 1.0};
//+
Point(6) = {0.750430, 0.034510, 0, 1.0};
//+
Point(7) = {0.700430, 0.041820, 0, 1.0};
//+
Point(8) = {0.650380, 0.048700, 0, 1.0};
//+
Point(9) = {0.600290, 0.054910, 0, 1.0};
//+
Point(10) = {0.550160, 0.060300, 0, 1.0};
//+
Point(11) = {0.50, 0.064730, 0, 1.0};
//+
Point(12) = {0.449820, 0.067990, 0, 1.0};
//+
Point(13) = {0.399620, 0.069910, 0, 1.0};
//+
Point(14) = {0.349410, 0.0703, 0, 1.0};
//+
Point(15) = {0.2992, 0.06901, 0, 1.0};
//+
Point(16) = {0.249, 0.06606, 0, 1.0};
//+
Point(17) = {0.19882, 0.06137, 0, 1.0};
//+
Point(18) = {0.14868, 0.0547, 0, 1.0};
//+
Point(19) = {0.098590, 0.04554, 0, 1.0};
//+
Point(20) = {0.073580, 0.039630, 0, 1.0};
//+
Point(21) = {0.048630, 0.03238, 0, 1.0};
//+
Point(22) = {0.02378, 0.02284, 0, 1.0};
//+
Point(23) = {0.01145, 0.01622, 0, 1.0};
//+
Point(24) = {0.00657, 0.0126, 0, 1.0};
//+
Point(25) = {0.00417, 0.01032, 0, 1.0};
//+
Point(26) = {0, 0, 0, 1.0};
//+
Point(27) = {0.00538, -0.00932, 0, 1.0};
//+
Point(28) = {0.00843, -0.0112, 0, 1.0};
//+
Point(29) = {0.01355, -0.01408, 0, 1.0};
//+
Point(30) = {0.02622, -0.01912, 0, 1.0};
//+
Point(31) = {0.05137, -0.02606, 0, 1.0};
//+
Point(32) = {0.07642, -0.03115, 0, 1.0};
//+
Point(33) = {0.10141, -0.0352, 0, 1.0};
//+
Point(34) = {0.15132, -0.04124, 0, 1.0};
//+
Point(35) = {0.20118, -0.04545, 0, 1.0};
//+
Point(36) = {0.251, -0.04816, 0, 1.0};
//+
Point(37) = {0.3008, -0.04957, 0, 1.0};
//+
Point(38) = {0.35059, -0.0497, 0, 1.0};
//+
Point(39) = {0.40038, -0.04849, 0, 1.0};
//+
Point(40) = {0.45018, -0.04609, 0, 1.0};
//+
Point(41) = {0.5, -0.04267, 0, 1.0};
//+
Point(42) = {0.54848, -0.0384, 0, 1.0};
//+
Point(43) = {0.59971, -0.03349, 0, 1.0};
//+
Point(44) = {0.64962, -0.0281, 0, 1.0};
//+
Point(45) = {0.69957, -0.02238, 0, 1.0};
//+
Point(46) = {0.74955, -0.01661, 0, 1.0};
//+
Point(47) = {0.79958, -0.01106, 0, 1.0};
//+
Point(48) = {0.84965, -0.00601, 0, 1.0};
//+
Point(49) = {0.89975, -0.0019, 0, 1.0};
//+
Point(50) = {0.94988, 0.00066, 0, 1.0};
//+
Line(1) = {1, 2};
//+
Line(2) = {2, 3};
//+
Line(3) = {3, 4};
//+
Line(4) = {4, 5};
//+
Line(5) = {5, 6};
//+
Line(6) = {6, 7};
//+
Line(7) = {7, 8};
//+
Line(8) = {8, 9};
//+
Line(9) = {9, 10};
//+
Line(10) = {10, 11};
//+
Line(11) = {11, 12};
//+
Line(12) = {12, 13};
//+
Line(13) = {13, 14};
//+
Line(14) = {14, 15};
//+
Line(15) = {15, 16};
//+
Line(16) = {16, 17};
//+
Line(17) = {17, 18};
//+
Line(18) = {18, 19};
//+
Line(19) = {19, 20};
//+
Line(20) = {20, 21};
//+
Line(21) = {21, 22};
//+
Line(22) = {22, 23};
//+
Line(23) = {23, 24};
//+
Line(24) = {24, 25};
//+
Line(25) = {25, 26};
//+
Line(26) = {26, 27};
//+
Line(27) = {27, 28};
//+
Line(28) = {28, 29};
//+
Line(29) = {29, 30};
//+
Line(30) = {30, 31};
//+
Line(31) = {31, 32};
//+
Line(32) = {32, 33};
//+
Line(33) = {33, 34};
//+
Line(34) = {34, 35};
//+
Line(35) = {35, 36};
//+
Line(36) = {36, 37};
//+
Line(37) = {37, 38};
//+
Line(38) = {38, 39};
//+
Line(39) = {39, 40};
//+
Line(40) = {40, 41};
//+
Line(41) = {41, 42};
//+
Line(42) = {42, 43};
//+
Line(43) = {43, 44};
//+
Line(44) = {44, 45};
//+
Line(45) = {45, 46};
//+
Line(46) = {46, 47};
//+
Line(47) = {47, 48};
//+
Line(48) = {48, 49};
//+
Line(49) = {49, 50};
//+
Line(50) = {50, 1};
//+
Point(51) = {2.3, 4, 0, 1.0};
//+
Point(52) = {2.3, -4, 0, 1.0};
//+
Point(53) = {5.2, 4, 0, 1.0};
//+
Point(54) = {5.2, -4, 0, 1.0};
//+
Point(55) = {2.6, 0, 0, 1.0};
//+
Circle(51) = {51, 55, 52};
//+
Line(52) = {51, 53};
//+
Line(53) = {53, 54};
//+
Line(54) = {54, 52};//+
//+
Curve Loop(1) = {54, -51, 52, 53};
//+
Curve Loop(2) = {13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12};
//+
Plane Surface(1) = {1, 2};
//+
Physical Curve("freestream", 55) = {52, 53, 54, 51};
//+
Physical Curve("airfoil", 56) = {10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 50, 49, 48, 47, 46, 45, 44, 43, 42, 41, 40, 39, 38, 37, 36, 35, 34, 33, 32, 31, 30, 29, 28, 27, 26, 25, 24, 23, 22, 21, 20, 19, 18, 17, 16, 15, 14, 13, 12, 11};
//+
Physical Surface("interior", 57) = {1};

[PhilipA]

I see, nice spot. Thank you very much for taking the time to help me out with this problem. It is great to see it working now!

Do you have any tips for ensuring that the mesh connectivity and points are saved on the first attempt?

Thanks again