I got the error :

D:\Nanokappa-main>python nanokappa.py -ff parameters_test.txt
Traceback (most recent call last):
File "D:\Nanokappa-main\nanokappa.py", line 2, in
from classes.Geometry import Geometry
File "D:\Nanokappa-main\classes\Geometry.py", line 11, in
import trimesh as tm
ModuleNotFoundError: No module named 'trimesh'

Trimesh should be install, can yo clarify which version?

Opening this issue to keep it on the record.

• Problem appearing on Laurent's macbook.

• Triangle and pyembree were installed via pip. The rest was installed from conda.

• Invalid division happening in this part of the code:

    offset_dist_norm = - v_norm*self.offset/np.sum(velocities*normals, axis = 1, keepdims = True) # ||d|| = ||v|| h / v . n
    offset_dist      = - velocities*offset_dist_norm/v_norm                                       # d = ||d|| v/||v||
  

• Does not happen on Windows nor Linux.

The volume triangulation routine currently uses scipy.delaunay.triangulate to divide the mesh into tetrahedra. Despite Scipy's routine being applied only to convex geometries, additional checks in Nanokappa are made to detect holes and corners, as well as new points are added to avoid holes inside the mesh that should not be there. Nevertheless sometimes these are not enough, and we still find particles outside of the mesh being detected as being inside, and vice versa.

There is a bounds check using the bounding box of the geometry, eventually particles that are outside will be moved inside. However this may take a lot of time if their velocity is small or the geometry is large. Also this does not mean it solves the triangulation issue inside the mesh.

Example geometry that reproduces this problem:

--geometry         freewire
--dimensions       100 200 100 0 300 200 150 0 250 200 100 0 150 200 50 0 100 200 10 6    

The way these calculations are being made seem inconsistent.

These may be a consequence of a combination of:

1. the lack of precision of T estimation throughout the volume and
2. the boundary scattering model.

The heat flux calculation actually seems to be faulty, not changing properly with area and Delta T.