Hello friends,
I am trying to use the mesh/tet command and so far I had no success.
I used 2 meshes in the study. The first mesh is a STL mesh. The second is a VTK mesh. Both meshes are in the same scale (mm). I scaled the meshes by 1e-3, because the unit system of the study is SI. If I run a test without the VTK mesh, everthing run ok. But if I run the study with the VTK mesh, liggghts return this error:
ERROR: VTK mesh file is incompatible with simulation box: One or more vertices outside simulation box (input_mesh_tet.cpp:298)
If I scale the VTK mesh by 1e-4, the study run without problem.
The question is why liggghts return this problem. Both meshes are in the same dimensions.
This figure shows the meshes:
https://mega.nz/#!89FSSZgZ!ehGD7NSK7GeDwnFi0dt_0Y5oj0QbPPSczEpuB6dmhHM
The files of the study:
cad0.stl: https://mega.nz/#!819DkKQR!54Mu9KAzTsDzEyrQRpxY7rMrZEQUtWzT0PasNo195SI
cad1.vtk: https://mega.nz/#!JssXmBxL!LK7znSVf4Bw0ONNnrRAc0sI2NnI5FOQkUTkPKDgnntE
in.study: https://mega.nz/#!9kUBXaTQ!Z4UMPnt6CNdhtd_gJJ8Ov-7VJdamZovPe_fpdYJlZOo
data.head: https://mega.nz/#!d9U0DbSD!f07VbVcum7UhnNqMc3x10Qk3ca64vnsLl8whY2XxUak
The VTK mesh was generate by Hypermesh 14 with this template:
VTK Hypermesh template: https://mega.nz/#!p10QAAAR!syO5o4mQc9azQiotkO7LSfa26xm6y5tUyV3-RwfOTe0
I thank you in advance.
Daniel Oliveira
| Attachment | Size |
|---|---|
 selection_006.png | 55.29 KB |
arnom | Wed, 07/05/2017 - 15:09
You are reading the tet file
You are reading the tet file before creating the domain. This results in default values for the domain being taken and because of this you get the error.
DanielBR | Tue, 08/22/2017 - 00:11
Thank you
Thank you arnom. I modified my study file and I can run the study now.
But, I have faced another problem.
My study run until the timestep 5778 and suddenly liggghts crash with this message.
INFO: Particle insertion insert0: inserted 138799 particle templates (mass 4.545524e+04) at step 1
- a total of 138799 particle templates (mass 4.545524e+04) inserted so far.
321 0.00015582891 138799 0.050021079 0.0028846309 18.511525 3449.311
963 0.00015582891 138799 0.15006324 0.0024352799 15.627908 3758.8042
1605 0.00015582891 138799 0.25010539 0.0027843542 17.868452 3632.1326
2247 0.00015582891 138799 0.35014755 0.0058024432 37.235987 3189.0578
2889 0.00015582891 138799 0.45018971 0.0038323181 24.593114 2951.1745
3531 0.00015582891 138799 0.55023187 0.0026175719 16.79773 3020.2099
4173 0.00015582891 138799 0.65027403 0.0023755511 15.244611 3110.5615
4815 0.00015582891 138799 0.75031618 0.0020498718 13.154631 3201.6214
5457 0.00015582891 138799 0.85035834 0.0027661627 17.751281 3330.2649
5778 0.00015582891 138799 0.90037942 0.0026882973 17.251596 3316.1053
ERROR on proc 0: This feature is not available for tet mesh regions (region_mesh_tet.cpp:201)
ERROR on proc 1: This feature is not available for tet mesh regions (region_mesh_tet.cpp:201)
ERROR on proc 4: This feature is not available for tet mesh regions (region_mesh_tet.cpp:201)
ERROR on proc 5: This feature is not available for tet mesh regions (region_mesh_tet.cpp:201)
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 4
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 5
I suspect that there are a problem with "Insert/Rate/Region" command:
fix insert0 nve_group insert/rate/region seed 49979687 distributiontemplate partic01 insert_every ${every} overlapcheck yes vel constant 0. 0. 0. mass ${TotalMass} massrate ${FactMassrate} region inReg ntry_mc 10000
This is my new study script: https://mega.nz/#!RhMS3YrB!KJeF3lKCswokHsb--B-1x1Ur0dhhEXIw8rYT35_WNJs
I tried many things but I could not solve or understand this.
Can you help with this?
Thank you very much.
arnom | Tue, 08/22/2017 - 09:10
I'm not sure why this error
I'm not sure why this error appears after a certain time. One guess is the following: Is it possible that particles that already have been inserted are located close to or inside the insertion region?
tjleps | Wed, 02/28/2018 - 21:45
close mesh/tet
Is it not possible to insert particles in two mesh/tet regions that are close together? I'm trying to do this currently and getting the same error.
arnom | Tue, 08/22/2017 - 09:16
I can't run your new script
I can't run your new script with the meshes above. If the problem persists, please put all your data in one zip file and upload that to mega so that I can reproduce it here.
DanielBR | Tue, 08/29/2017 - 06:32
You did make the question
You did make the question "One guess is the following: Is it possible that particles that already have been inserted are located close to or inside the insertion region?" and I suspect that this is the problem. But, I insert mass in the same frequency when I use a region/block command and the study run without problem.
The compressed file with all files of the study: https://mega.nz/#!Q98xlRYB!f3prxbGMJdW1S3IV8VzQK0QE8l8Gu9PUH1iqwRG6IUs
Thank you very much.