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

Actuator disk mesh generation problem with GMSH

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 30, 2023, 08:37
Unhappy Actuator disk mesh generation problem with GMSH
  #1
New Member
 
Tony Govoni
Join Date: Mar 2023
Posts: 1
Rep Power: 0
tgovoni is on a distinguished road
Hello,

I am trying to generate a SU2 simulation of a simple actuator disk as the one proposed in the testcase propeller.cfg (https://github.com/su2code/SU2/tree/...rans/propeller)

I use GMSH for the generation of the geometry and the mesh and in order to get started and see if everything was okay I made a first very simple test case with a disk in a box.

I am using the exact same .cfg file of the propeller test case, except that I dont have the nacelle surfaces so I commented them.

Before changing anything related to the simulation, it is impossible for me to get this first case running since I always have some errors concerning the mesh.

I tried a lot of different things but I get always the two same errors :

either:

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.


or if I activate the option SaveAll of GMSH.:

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.


I cannot understand where this error is coming from, my mesh file is "simillar" to the one used for the propeller test case.

I tried to switch from different version of GMSH and SU2 but nothing solve this issue.

I am currently using GMSH 4.9.5 (as the latest version generate additional errors with the physical groups of volumes elements) and the latest SU2 8.0.0

Here is the python code generating the mesh (python 3.11.4) :

Code:
import gmsh

gmsh.initialize()

gmsh.model.add("mesh actuator disk")

# can redefine a namespace:
model = gmsh.model
model_kernel = model.occ
model_mesh = model.mesh

# define a volume by a box of 6 surface with one interal surface that is a disk in the center of the box

# box : center is in 0,0,0 and size is 1,1,1
box = model_kernel.addBox(-0.5, -0.5, -0.5, 1, 1, 1)
# disk : center is in 0,0,0 and radius is 0.1
disk = model_kernel.addDisk(0, 0, 0, 0.1, 0.1)
model_kernel.synchronize()

# emmebed the disk in the box for the mesh
model_mesh.embed(2, [disk], 3, box)

# define physical groups
# 1 : the box walls will be named FF
# 2 : the disk will be named DISK

# get the box walls
tag_box = [1, 2, 3, 4, 5, 6]
tag_disk = [7]

# define physical groups
model.addPhysicalGroup(2, tag_box, 1)
model.setPhysicalName(2, 1, "FARFIELD")
model.addPhysicalGroup(2, tag_disk, 2)
model.setPhysicalName(2, 2, "ACTDISK")
model.addPhysicalGroup(3, [1], 3)
model.setPhysicalName(3, 3, "fluid")
# mesh
model_mesh.generate(3)


# save all the mesh
gmsh.option.setNumber("Mesh.SaveAll", 1)
# save the mesh .su2 format
gmsh.write("ActuatorDisk.su2")
Someone as any Idea of whats going on here ?

Many thanks
tgovoni is offline   Reply With Quote

Reply

Tags
gmsh, npoinerror, su2 8.0.0

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
Mesh Generation problem Kelysamex STAR-CCM+ 2 August 6, 2021 02:50
[snappyHexMesh] SHM problem : KVLCC2 with appendage mesh sc.park OpenFOAM Meshing & Mesh Conversion 1 March 13, 2016 13:28
[Gmsh] Scripted version of "2D Mesh Generation Tutorial for GMSH" laubeg OpenFOAM Meshing & Mesh Conversion 1 April 14, 2013 08:32
[snappyHexMesh] snappyHexMesh won't work - zeros everywhere! sc298 OpenFOAM Meshing & Mesh Conversion 2 March 27, 2011 21:11
unstructured vs. structured grids Frank Muldoon Main CFD Forum 1 January 5, 1999 10:09


All times are GMT -4. The time now is 00:37.