Unrecognized electrode nodes
We used gmsh
to produce an ERT model. In this model, a fissure was set as a cube in the cylinder, in which three little cylinders are across it. The electrodes, a series of points, were arrayed along the edge of the cube. However, there existed some problems during the performance:
- The added
point
s (electrodes) could not be recognized by the software; - When making the
point
s in thevolume
, thecylinder
could not be recognized bygmsh
.
I wonder what the mistake is and how to deal with it?
The detailed code are listed as follows:
# Q1
> for i in range(1, 29):
> p = gmsh.model.occ.addPoint(R, 0, i * 0.1 - Hf/2, 0, 1000+i)
> mesh.findCell([1, 0, 0])
Reading hb1.msh2...
Nodes: 2174
Entries: 6269
Points: 28
Lines: 0
Triangles: 0
Quads: 0
Tetrahedra: 6241
Creating mesh object...
Dimension: 3-D
WARNING: No boundary conditions found. Setting Neumann on the outer edges by default.
Regions: 4 (1, 2, 3, 4)
Marked nodes: 28 (99,)
Done.
Mesh: Nodes: 2174 Cells: 6241 Boundaries: 14333
Node: Node: 0x1e462c24670 id: 62 1 0 0 marker: -99
Traceback (most recent call last):
File ".\hb2.py", line 12, in <module>
mesh.findCell([1, 0, 0])
RuntimeError: C:/msys64/home/halbm/gimli/gimli/core/src/mesh.cpp:864 GIMLI::Cell* GIMLI::Mesh::findCell(const RVector3&, size_t&, bool) const no cells for this node. This is a corrupt mesh
# 2
> gmsh.model.occ.addCylinder(0, 0, -H/2, 0, 0, H, R, 1)
> gmsh.model.mesh.embed(0, [p], 3, 1)
C:\ProgramData\Anaconda3\lib\site-packages\numpy\__init__.py:148: UserWarning: mkl-service package failed to import, therefore Intel(R) MKL initialization ensuring its correct out-of-the box operation under condition when Gnu OpenMP had already been loaded by Python process is not assured. Please install mkl-service package, see http://github.com/IntelPython/mkl-service
from . import _distributor_init
Error : Volume 1 does not exist
Traceback (most recent call last):
File "c:\Users\Hsueh\Desktop\hb\hb1.py", line 52, in <module>
gmsh.model.mesh.embed(0, [p], 3, 1)
File "C:\ProgramData\Anaconda3\lib\site-packages\gmsh.py", line 3505, in embed
raise Exception(logger.getLastError())
Exception: Volume 1 does not exist
[Done] exited with code=1 in 0.362 seconds
[Running] python -u "c:\Users\Hsueh\Desktop\hb\hb1.py"
C:\ProgramData\Anaconda3\lib\site-packages\numpy\__init__.py:148: UserWarning: mkl-service package failed to import, therefore Intel(R) MKL initialization ensuring its correct out-of-the box operation under condition when Gnu OpenMP had already been loaded by Python process is not assured. Please install mkl-service package, see http://github.com/IntelPython/mkl-service
from . import _distributor_init
Error : Volume 1 does not exist
Traceback (most recent call last):
File "c:\Users\Hsueh\Desktop\hb\hb1.py", line 52, in <module>
gmsh.model.mesh.embed(0, [p], 3, 1)
File "C:\ProgramData\Anaconda3\lib\site-packages\gmsh.py", line 3505, in embed
raise Exception(logger.getLastError())
Exception: Volume 1 does not exist
and the source code of python
is attached as follows: hb1.py