| ... | ... | @@ -4,15 +4,15 @@ It is important to note that within the granular flows field there are many diff |
|
|
|
|
|
|
|
One of the greatest advantage of the MigFlow Software is its compatibility with the [LMGC90](https://git-xen.lmgc.univ-montp2.fr/lmgc90/lmgc90_user/wikis/home) Software (for installation we recommand to follow the instructions on their website). A contact solver, named _scontact_ is provided with MigFlow to compute the velocities and trajectories of the grains, forbidding interpenetration, in a lagrangian way. It uses a time-stepping method that corrects iteratively the free velocities of the grains, obtained by applying external forces, to compute a set of velocities giving after displacement of the grains a stationary state without intersection. But the _scontact_ module can be disconnected from the fluid part and replaced by LMGC90 to take advantage of the huge variety of contact laws implemented in it. Some test cases are given to present how to call LMGC90 from the Users Interface:
|
|
|
|
|
|
|
|
- [melangeur.py](https://git.immc.ucl.ac.be/fluidparticles/migflow/blob/oneFluid/testcases/couette-2d/betonlmgc.py) presents the mixing of grains with different radii in a rotating drum using LMGC90 to take into account the impact of the friction on the mixing efficiency.
|
|
|
|
- [depot.py](https://git.immc.ucl.ac.be/fluidparticles/migflow/blob/oneFluid/testcases/smallDepot-3d/depot.py) presents a tridimensionoal deposit in a box filled with water.
|
|
|
|
- [melangeur.py](https://git.immc.ucl.ac.be/fluidparticles/migflow/blob/master/testcases/couette-2d/betonlmgc.py) presents the mixing of grains with different radii in a rotating drum using LMGC90 to take into account the impact of the friction on the mixing efficiency.
|
|
|
|
- [depot.py](https://git.immc.ucl.ac.be/fluidparticles/migflow/blob/master/testcases/smallDepot-3d/depot.py) presents a tridimensional deposit in a box filled with water.
|
|
|
|
|
|
|
|
### Adaptive Mesh
|
|
|
|
|
|
|
|
MigFlow is based on an unresolved Computational Fluid Dynamics-Discrete Element Method model. The discretisation method we choose for the CFD part is the Finite Element Method. On the one hand, unstructured grids offer a great flexibility to represent complex geometries with which grains interact. On the other hand, it is easy to adapt the mesh in order to refine it around the area of interest to accurately capture flow fields during the grains motion.
|
|
|
|
|
|
|
|
- [drop.py](https://git.immc.ucl.ac.be/fluidparticles/migflow/blob/oneFluid/testcases/drop-2d/SimpleDrop/drop.py) presents the fall of a cloud made of grains in a viscous fluid and makes use of this feature to refine the mesh around the cloud area during the fall. This way, it is possible to highly fasten the computation because the area of interest (i.e. the cloud) is small compared with the entire computational domain.
|
|
|
|
- [avalanch1fluidnofriction.py](https://git.immc.ucl.ac.be/fluidparticles/migflow/blob/oneFluid/testcases/avalanch/avalanch1fluid/avalanch1fluidnofriction.py) presents the fall of a granular column immersed in soft water.
|
|
|
|
- [drop.py](https://git.immc.ucl.ac.be/fluidparticles/migflow/blob/master/testcases/drop-2d/SimpleDrop/drop.py) presents the fall of a cloud made of grains in a viscous fluid and makes use of this feature to refine the mesh around the cloud area during the fall. This way, it is possible to highly fasten the computation because the area of interest (i.e. the cloud) is small compared with the entire computational domain.
|
|
|
|
- [avalanch1fluidnofriction.py](https://git.immc.ucl.ac.be/fluidparticles/migflow/blob/master/testcases/avalanch/avalanch1fluid/avalanch1fluidnofriction.py) presents the fall of a granular column immersed in soft water.
|
|
|
|
|
|
|
|
```python
|
|
|
|
fluid.adapt_mesh(Max. element size, Min. element size, Total number of elements, a, b)
|
| ... | ... | @@ -24,4 +24,4 @@ fluid.adapt_mesh(Max. element size, Min. element size, Total number of elements, |
|
|
|
|
|
|
|
To consider two continuous phases, an advective equation for the concentration is added to the standard Navier-Stokes equations. See the following example to run test cases with two different continuous phases.
|
|
|
|
|
|
|
|
- [avalanch2fluids.py](https://git.immc.ucl.ac.be/fluidparticles/migflow/blob/oneFluid/testcases/avalanch/avalanch2fluids/avalanch2fluids.py) presents the fall of a granular column immersed in sea water in a 2D domain filled with soft water. |
|
|
|
- [avalanch2fluids.py](https://git.immc.ucl.ac.be/fluidparticles/migflow/blob/master/testcases/avalanch/avalanch2fluids/avalanch2fluids.py) presents the fall of a granular column immersed in sea water in a 2D domain filled with soft water. |
