Hello everyone,
I am creating this mesh in order to use it as a size function. When creating the mesh from a geometry everything seems to work fine, however the problem occurs when importing the mesh again.
Before exporting the mesh, I have 4 volumes whose IDs are ordered from 1 to 4. However, when exporting and re-importing the mesh, the IDs of these volumes change, and it happens that for volumes 3 and 10 the same IDs are shared and it causes a lot a troubles. Does anyone knows how to solve it?
This is the way I am creating the mesh
# -*- Python -*- (syntax highlighting)
# ----------------------------------------------------------------------
reset
${Units("si")} # --------------------------------------------
# Import slab, name the surfaces and body for future reference
# --------------------------------------------
import Acis 'geometry.sat'
# ----------------------------------------------------------------------
# Create tet4 mesh at 4.0 km resolution.
# ----------------------------------------------------------------------
${dx=40.0*km}
volume all size {dx}
volume all scheme tetmesh
# ----------------------------------------------------------------------
# Generate the mesh
# ----------------------------------------------------------------------
mesh surface all
mesh volume all
# ----------------------------------------------------------------------
# Export exodus file
# ----------------------------------------------------------------------
set large exodus off
# Create one block so all sizing function information is in one block.
block 1 volume all
export mesh "mesh_cellsize.exo" dimension 3 overwrite
