29 add funel example with polyhedra (#32)

* add funnel example

* It was not the good box

* push example

* better perfs

---------

Co-authored-by: rprat-pro <rp269144@pleiades1217.intra.cea.fr>

example/polyhedra/funnel/alpha3.shp

@@ -0,0 +1,25 @@
+
+<
+name alpha3
+radius 0.237770037899799
+preCompDone y
+nv 6
+0.475540075799599 0 0
+-0.475540075799599 0 0
+0 0.475540075799599 0
+0 -0.475540075799599 0
+0 0 0.475540075799599
+0 0 -0.475540075799599
+ne 3
+0 1
+2 3
+4 5
+nf 0
+obb.extent 0.713310113699398 0.713310113699398 0.713310113699398
+obb.e1 1 0 0
+obb.e2 0 1 0
+obb.e3 0 0 1
+obb.center 0 0 0
+volume 0.523598775598299
+I/m 0.13266 0.13266 0.13266
+>

example/polyhedra/funnel/box.stl

@@ -0,0 +1,72 @@
+solid Visualization Toolkit generated SLA File
+ facet normal -0 -1 -0
+  outer loop
+   vertex 258.57806396484375 20.84356689453125 -70
+   vertex 258.57806396484375 20.84356689453125 -37.99993896484375
+   vertex 19.180915832519531 20.84356689453125 -37.99993896484375
+  endloop
+ endfacet
+ facet normal 0 -1 0
+  outer loop
+   vertex 258.57806396484375 20.84356689453125 -70
+   vertex 19.180915832519531 20.84356689453125 -37.99993896484375
+   vertex 19.180915832519531 20.84356689453125 -70
+  endloop
+ endfacet
+ facet normal -1 0 0
+  outer loop
+   vertex 19.180915832519531 20.84356689453125 -70
+   vertex 19.180915832519531 20.84356689453125 -37.99993896484375
+   vertex 19.180915832519531 260.337158203125 -37.99993896484375
+  endloop
+ endfacet
+ facet normal -1 0 -0
+  outer loop
+   vertex 19.180915832519531 20.84356689453125 -70
+   vertex 19.180915832519531 260.337158203125 -37.99993896484375
+   vertex 19.180915832519531 260.337158203125 -70
+  endloop
+ endfacet
+ facet normal -0 -0 -1
+  outer loop
+   vertex 19.180915832519531 260.337158203125 -70
+   vertex 258.57806396484375 260.337158203125 -70
+   vertex 258.57806396484375 20.84356689453125 -70
+  endloop
+ endfacet
+ facet normal 0 0 -1
+  outer loop
+   vertex 19.180915832519531 260.337158203125 -70
+   vertex 258.57806396484375 20.84356689453125 -70
+   vertex 19.180915832519531 20.84356689453125 -70
+  endloop
+ endfacet
+ facet normal 1 0 0
+  outer loop
+   vertex 258.57806396484375 260.337158203125 -70
+   vertex 258.57806396484375 260.337158203125 -37.99993896484375
+   vertex 258.57806396484375 20.84356689453125 -37.99993896484375
+  endloop
+ endfacet
+ facet normal 1 0 0
+  outer loop
+   vertex 258.57806396484375 260.337158203125 -70
+   vertex 258.57806396484375 20.84356689453125 -37.99993896484375
+   vertex 258.57806396484375 20.84356689453125 -70
+  endloop
+ endfacet
+ facet normal -0 1 0
+  outer loop
+   vertex 19.180915832519531 260.337158203125 -70
+   vertex 19.180915832519531 260.337158203125 -37.99993896484375
+   vertex 258.57806396484375 260.337158203125 -37.99993896484375
+  endloop
+ endfacet
+ facet normal 0 1 0
+  outer loop
+   vertex 19.180915832519531 260.337158203125 -70
+   vertex 258.57806396484375 260.337158203125 -37.99993896484375
+   vertex 258.57806396484375 260.337158203125 -70
+  endloop
+ endfacet
+endsolid

example/polyhedra/funnel/funnel.msp

@@ -0,0 +1,68 @@
+grid_flavor: grid_flavor_dem
+
+includes:
+  - config_polyhedra.msp
+
+input_data:
+  - lattice:
+      structure: SC
+      types: [ 0 ]
+      size: [ 1.5 , 1.5 , 1.5 ]
+      repeats: [ 200 , 200 , 150 ]
+      enlarge_bounds: 0.0 m
+      region: FUNNEL and BZ
+  - read_shape_file:
+     filename: alpha3.shp
+  - radius_from_shape
+  - density_from_shape
+  - set_rand_velocity:
+     var: 0.1
+     mean: [0.0,0.0,0.0]
+  - set_rand_vrot_arot
+  - set_quaternion
+  - inertia_from_shape
+
+
+setup_drivers:
+  - add_stl_mesh:
+     id: 0
+     filename: box.stl
+     minskowski: 0.01
+  - add_stl_mesh:
+     id: 1
+     filename: funnel.stl
+     minskowski: 0.01
+  - update_grid_stl_mesh
+
+particle_regions:
+   - FUNNEL:
+      quadric:
+        shape: cylz
+        transform:
+          - scale: [120, 120, 1]
+          - translate: [150,150, 0]
+   - BZ:
+      bounds: [[0, 0, 110],[300,300, 210]]
+
+compute_force:
+  - gravity_force
+  - compute_hooke_interaction:
+     config: { rcut: 0.0 m , dncut: 0.0 m, kn: 10000, kt: 8000, kr: 0.0, fc: 0.0, mu: 0.1, damp_rate: 0.999} 
+     config_driver: { rcut: 0.0 m , dncut: 0.0 m, kn: 10000, kt: 8000, kr: 0.0, fc: 0.0, mu: 0.1, damp_rate: 0.999} 
+
+domain:
+  cell_size: 4 m
+  periodic: [false,false,false]
+
+write_vtklegacy:
+  ascii: true
+  ghost: false
+
+global:
+  simulation_dump_frequency: 100000
+  simulation_end_iteration: 200000
+  simulation_log_frequency: 100
+  simulation_paraview_frequency: 1000
+  dt: 0.0001 s 
+  rcut_inc: 0.4 m
+  enable_stl_mesh: true