Error exporting mesh from SALOME to SU2

Hi,

I’m trying to export a 3D mesh from SALOME to use in SU2 for CFD, but I’m running into errors.

  1. When I export directly to SU2 format, SALOME shows: “Export Failed: TypeError: Cannot unpack non-iterable CellBlock object.”
  2. If I export in CGNS format instead, it completes without errors, but SU2 fails to read it, showing: "Parsing config file zone 0, Assertion failed: head && *head file ../externals/cgns/hdf5/H5CX.c, line 3579.”

Is there a fix or workaround for this issue?

Regards,
N469

I had that kind of errors also when I tried to use directly the meshio suite, I think it is most likely a problem there. If I may suggest, I’d rather use the cfdmsh.py to export natively from SALOME to SU2.

I have been using the same script, but I often get the NPOIN mismatch error, especially while using viscous layers.

Regards,
N469

Hi N469, I don’t recall having that issue, can you be more specific ? NPOIN error, I understand you are having that on the SU2 side, do you check in SALOME whether there are double nodes, prior to the export ? Are you considering a 3D or 2D case?

Hi Giovanni,

I’m working on a full aircraft CFD simulation in SU2 and have attached the relevant scripts and error logs for reference. When I run the 3D CFD case without viscous layers, it completes successfully in both serial and parallel modes.

However, once I add viscous layers, the meshing process finishes without any issues and exports correctly using cfdmsh.py, but SU2 then produces different errors depending on whether the run is in serial or parallel. I’ve already tried refining the mesh and testing various approaches, but the problem consistently appears as soon as viscous layers are included.

When running in parallel, I encounter the following error:

------------------- Geometry Preprocessing ( Zone 0 ) -------------------
Three dimensional problem.
507036 grid points before partitioning.
961980 volume elements before partitioning.
3 surface markers.
19732 boundary elements in index 0 (Marker = SYMM).
230 boundary elements in index 1 (Marker = FARFIELD).
32066 boundary elements in index 2 (Marker = BODY).
Executing the partitioning functions.
Building the graph adjacency structure.
Calling ParMETIS… graph partitioning complete (8056 edge cuts).
Distributing ParMETIS coloring.

job aborted:
[ranks] message

[0] process exited without calling finalize

[1-3] terminated

-— error analysis -----

[0] on LAPTOP-3RJTS8M9
D:\CRM\salomefiles\P30\SU2_CFD ended prematurely and may have crashed. exit code 0xc0000005

-— error analysis -----
Traceback (most recent call last):
File “D:\CRM\salomefiles\P30\parallel_computation3.py”, line 126, in
main()
File “D:\CRM\salomefiles\P30\parallel_computation3.py”, line 72, in main
parallel_computation(options.filename, options.partitions, options.compute)
File “D:\CRM\salomefiles\P30\parallel_computation3.py”, line 103, in parallel_computation
info = SU3.run.CFD(config)
File “D:\CRM\salomefiles\P30\SU3\run\interface.py”, line 118, in CFD
run_command(the_Command)
File “D:\CRM\salomefiles\P30\SU3\run\interface.py”, line 295, in run_command
raise exception(message)
RuntimeError: Path = D:\CRM\salomefiles\P30,
Command = mpiexec -n 4 “D:\CRM\salomefiles\P30\SU2_CFD” “D:\CRM\salomefiles\P30\config_CFD.cfg”
SU2 process returned error ‘4294967293’

Whereas, the same mesh and configuration file in serial modes gives:

------------------- Geometry Preprocessing ( Zone 0 ) -------------------
Three dimensional problem.
507036 grid points.
961980 volume elements.
3 surface markers.
19732 boundary elements in index 0 (Marker = SYMM).
230 boundary elements in index 1 (Marker = FARFIELD).
32066 boundary elements in index 2 (Marker = BODY).

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 -------------------------------

Please let me know if you have come across a similar issue or have any suggestions.

NP - No Prism
P20 - 20 Layers
P30 - 30 Layers

Regards,
N469

GEOM.zip (4.6 MB)

Mesh_1_NP.py (13.2 KB)

Mesh_1_P20.py (13.2 KB)

Mesh_1_P20.py (13.2 KB)

cfdmsh.py (391.4 KB)