CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM (https://www.cfd-online.com/Forums/openfoam/)
-   -   LPT incompressible fluid (https://www.cfd-online.com/Forums/openfoam/202974-lpt-incompressible-fluid.html)

Lighto June 13, 2018 07:42

LPT incompressible fluid
 
Hi everyone,

I'm a new member of the community.
I'm struggling with OpenFOAM to implement a lagrangian particle tracking on an incompressible fluid (i.e. blood cells in blood). I've followed many tutorials online but any of them worked for me. How can I implement LPT in icoFoam? And how can I post-process particles' properties such as their velocity and inducted shear stress?
I'm using OpenFOAM 5.0 on OpenSuse 42.2.

Thanks a lot,
Lighto

guin June 13, 2018 18:53

I think that the closest to what you are looking for is the solver DPMFoam. Indeed it is more a 2-way (kinematic) coupling between pimpleFoam + lagrangian particles with option for particle collisions.


The post-processing is quite easy:
1. paraFoam (load all objects)
2. use extractBlock filter to select your particles
3. use glyph filter e.g as spheres, scaled by the diameter and map onto them whatever property available from the list


Regarding the shear stress, I don't think there is such a thing for the particles in OF-5.0. May you rather be looking for something as the "immersed boundary" method? You can find an implementation in the foam-extend version, but keep in mind that it will not work efficiently if you want to simulate many particles simultaneously.

Lighto June 14, 2018 06:57

Thank you for the advice. Basically DMPFoam can solve pressure and velocity fields for water and simultanaously inject particles that follow the fluid velocity, is it right? So what I have to do is to specify the properties of the fluid inside trasportProperties and turbulenceProperties.fluid and than replace rho.air, U.air ecc with rho.fluid, U.fluid ecc inside the directory 0/ and the file system/fvSolution?


All times are GMT -4. The time now is 11:48.