region tet/mesh command

Submitted by alvriver on Wed, 04/26/2017 - 15:11

Hi!

I trying to use a particular STL volume like a particle generation region using the command "tet/mesh", but I'm having problems. The "tet/mesh" command ask me to use a "ASCII VTK file", so, I'm using PV to convert my STL volume to the VTK file, but this apparently doesn't works.

region bc mesh/tet file meshes/EmisorVTK2.vtk scale 0.001 side in units box

When I compile this, the program tell me:

region bc mesh/tet file meshes/EmisorVTK2.vtk scale 0.001 side in units box
ERROR: Illegal region mesh/tet command (../region_mesh_tet.cpp:70)
xk��� (../region_mesh_tet.cpp:70)

Help me please!

Álvaro Rivera Torres.

Daniel Queteschiner | Fri, 04/28/2017 - 16:07

STL files describe surface meshes, region mesh/tet requires a volume mesh made up of tetrahedrons. I doubt that paraview can create a tetrahedral volume from an STL surface.

alvriver | Wed, 05/03/2017 - 12:57

I solve my problem, but now I have this problem:

region bc mesh/tet file meshes/Emisorsolo01.1.vtk scale 0.001 move 0. 0. 0. rotate 0. 0. 0. units box side in
Note: Skipping empty line in VTK mesh file
Note: Skipping empty line in VTK mesh file
Note: Skipping empty line in VTK mesh file
Note: Skipping empty line in VTK mesh file
WARNING: Region mesh/tet: too many node neighbors, mesh is of bad quality; 'all_in yes' might not work correctly (../region_mesh_tet.cpp:460)
ERROR: Mesh (null): The minumum angle of mesh element 15 (line -1) is lower than the specified curvature. Increase mesh quality or decrease curvature (currently 0.256235°)
ERROR on proc 1: Fix mesh: Bad mesh, cannot continue. You can try setting 'curvature' to 1e-5 or lower or use 'curvature_tolerant yes' (../surface_mesh_I.h:647)
(../surface_mesh_I.h:647)
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 1 in communicator MPI_COMM_WORLD
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------
ERROR: Mesh (null): The minumum angle of mesh element 15 (line -1) is lower than the specified curvature. Increase mesh quality or decrease curvature (currently 0.256235°)
ERROR on proc 3: Fix mesh: Bad mesh, cannot continue. You can try setting 'curvature' to 1e-5 or lower or use 'curvature_tolerant yes' (../surface_mesh_I.h:647)
(../surface_mesh_I.h:647)
ERROR: Mesh (null): The minumum angle of mesh element 15 (line -1) is lower than the specified curvature. Increase mesh quality or decrease curvature (currently 0.256235°)
ERROR on proc 0: Fix mesh: Bad mesh, cannot continue. You can try setting 'curvature' to 1e-5 or lower or use 'curvature_tolerant yes' (../surface_mesh_I.h:647)
(../surface_mesh_I.h:647)
ERROR: Mesh (null): The minumum angle of mesh element 15 (line -1) is lower than the specified curvature. Increase mesh quality or decrease curvature (currently 0.256235°)
ERROR on proc 2: Fix mesh: Bad mesh, cannot continue. You can try setting 'curvature' to 1e-5 or lower or use 'curvature_tolerant yes' (../surface_mesh_I.h:647)
(../surface_mesh_I.h:647)
[ubuntu:05570] 3 more processes have sent help message help-mpi-api.txt / mpi-abort
[ubuntu:05570] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages

Help please!

Álvaro Rivera Torres.

tjleps | Mon, 02/26/2018 - 05:20

Have you managed to solve this problem, I can only find a curvature_tolerant command in mesh/surface.

richti83's picture

richti83 | Mon, 02/26/2018 - 15:23

As there is no parser interface between region mesh tet and surface_mesh the only way I know is to change the source code and recompile liggghts.
in file surface_mesh_I.h set line 59 curvature_tolerant_(true) and maybe minAngle_hardLimit_(0.9999995), to a higher value ....

Be warned, that to omit warnings can lead to unexpected behavior and that there is no support when something goes wrong.

original file: https://github.com/CFDEMproject/LIGGGHTS-PUBLIC/blob/master/src/surface_...

I'm not an associate of DCS GmbH and not a core developer of LIGGGHTS®
ResearchGate | Contact