Subsequent meshing does not remesh the boundary
Context
I'm trying to create a sequence of meshes with increasingly fine mesh size near certain points. The goal is to use these in a FEM simulation where I know that there is a singularity near these mesh points.
To do this I use a field to create a mesh with locally a very small mesh size. Then in a loop I decrease the values in the field and create a new mesh, thus creating a set of finer and finer meshes.
Problem
When calling Mesh 2
a second time, only the interior of the domain is remeshed but not the boundary. Consequently I get a mesh that is very fine in the interior but very coarse on the boundary. See the example below (corresponding example minimal.geo)
Expectation
A subsequent calls to Mesh
would remesh all parts. Alternatively, there should be a clear way to reset the state of the meshing algorithm so that the whole model gets remeshed, e.g. by being able to delete the current mesh.
Tried workarounds
I tried calling Mesh 1
before calling Mesh 2
. This works in this simple example. For a more complicated script it does not. sqrt2-mesh2-diamond-local.geo.
Info
- os: Arch linux
- install source: https://aur.archlinux.org/packages/gmsh-bin
- gmsh version 4.11.1