CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > SU2

SU2 issue with number of mesh points

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 12, 2023, 05:20
Default SU2 issue with number of mesh points
  #1
New Member
 
giulio romano
Join Date: Dec 2023
Posts: 3
Rep Power: 2
giulio romano is on a distinguished road
Hello to everyone, I'm new to this forum and I would ask some information about the mesh usage in SU2. I'm performing some calculation on the fuselage of an airplage. To generate the hybrid mesh I'm using GMSH but when I try to run the simulation in SU2 I obtain this error:

Error in "void CPhysicalGeometry:istributeColoring(const CConfig*, CGeometry*)":

-------------------------------------------------------------------------
Mismatch between NPOIN and number of points listed in mesh file.
Please check the mesh file for correctness.

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

application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0
[unset]: write_line error; fd=-1 buf=:cmd=abort exitcode=1
:
system msg for write_line failure : Bad file descriptor

Does someone know what is this issue about and if it is a problem in the mesh generation or in the reading of the mesh from SU2? There is a way to solve it while keep on using GMSH ? I attach the cfg file and the python script for the mesh to this post, so you can better understand the issue.
Thank you in advance
Attached Images
File Type: jpg mesh_image.jpg (115.6 KB, 11 views)
Attached Files
File Type: txt fuselage_cfg.txt (5.9 KB, 4 views)
File Type: txt fuselage_mesh_py.txt (2.7 KB, 5 views)
giulio romano is offline   Reply With Quote

Old   December 13, 2023, 15:07
Default
  #2
New Member
 
A Roy
Join Date: Sep 2022
Posts: 15
Rep Power: 3
A.Roy is on a distinguished road
Are you able to post either the mesh file or the iges geometry referenced in the gmsh python script to generate the mesh?
A.Roy is offline   Reply With Quote

Old   December 13, 2023, 15:36
Default
  #3
Senior Member
 
bigfoot
Join Date: Dec 2011
Location: Netherlands
Posts: 529
Rep Power: 17
bigfootedrockmidget is on a distinguished road
Common problem with such imported geometries is that it usually only imports the surface and the elements/cells were not generated in gmsh. Can you check in gmsh if you have any internal elements? You can also check in the su2 file if there are any elements (NELEM=0) in the file.
bigfootedrockmidget is offline   Reply With Quote

Old   December 14, 2023, 02:33
Default
  #4
New Member
 
giulio romano
Join Date: Dec 2023
Posts: 3
Rep Power: 2
giulio romano is on a distinguished road
Quote:
Originally Posted by A.Roy View Post
Are you able to post either the mesh file or the iges geometry referenced in the gmsh python script to generate the mesh?
Hello thank you for the reply here it is the iges file.
Attached Files
File Type: txt fuselage_bli_2D.iges.txt (34.0 KB, 1 views)
giulio romano is offline   Reply With Quote

Old   December 14, 2023, 02:44
Default
  #5
New Member
 
giulio romano
Join Date: Dec 2023
Posts: 3
Rep Power: 2
giulio romano is on a distinguished road
Quote:
Originally Posted by bigfootedrockmidget View Post
Common problem with such imported geometries is that it usually only imports the surface and the elements/cells were not generated in gmsh. Can you check in gmsh if you have any internal elements? You can also check in the su2 file if there are any elements (NELEM=0) in the file.
Hello thank you for the reply. I checked the SU2 file and the NELEMENT is not 0. But by looking at the SU2 mesh file I realized that the NPOIN = 50342 in my case but there are only 50341 points written below, i tryed to change in the mesh file NPOIN=50341 but SU2 gave me a different error:

"/usr/include/c++/8/bits/stl_vector.h:932: std::vector<_Tp, _Alloc>::reference std::vector<_Tp, _Alloc>:perator[](std::vector<_Tp, _Alloc>::size_type) [with _Tp = long unsigned int; _Alloc = std::allocator<long unsigned int>; std::vector<_Tp, _Alloc>::reference = long unsigned int&; std::vector<_Tp, _Alloc>::size_type = long unsigned int]: Assertion '__builtin_expect(__n < this->size(), true)' failed.
[1] 667514 abort (core dumped) SU2_CFD case15_BLI_fuselage_test3_GMSH.cfg"

How can I check in GMSH if I have any internal elements ?

I'm sorry that I can't attach to this post the mesh in SU2 and msh format so you can visualize them, because they are too large.
giulio romano is offline   Reply With Quote

Old   December 14, 2023, 05:07
Default
  #6
Senior Member
 
bigfoot
Join Date: Dec 2011
Location: Netherlands
Posts: 529
Rep Power: 17
bigfootedrockmidget is on a distinguished road
gmsh also uses helper points to construct your mesh, for instance the center of a circle is a helper point. Usually, these points are ignored if you deselect the option 'save all elements'. If this option is already deselected, then you could go through the su2 mesh points and check if there is a point that is not used in the mesh. I posted a python script for this a while ago on the forum. But the best thing is to solve it in gmsh.
bigfootedrockmidget is offline   Reply With Quote

Reply

Tags
gmsh, mismatch_npoint, su2

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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 Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
flow over cylinder in openFoam saeed jamshidi OpenFOAM Pre-Processing 3 August 11, 2023 15:16
Decomposing meshes Tobi OpenFOAM Pre-Processing 22 February 24, 2023 09:23
decomposePar problem: Cell 0contains face labels out of range vaina74 OpenFOAM Pre-Processing 37 July 20, 2020 05:38
foam-extend_3.1 decompose and pyfoam warning shipman OpenFOAM 3 July 24, 2014 08:14
[Commercial meshers] Trimmed cell and embedded refinement mesh conversion issues michele OpenFOAM Meshing & Mesh Conversion 2 July 15, 2005 04:15


All times are GMT -4. The time now is 18:50.