Question regarding the F2F particle tracking algorithm implementation.
I have read Niklas Nordin's thesis  about LPT. There he describes the face to face tracking algorithm (Chapter 4.2). Looking at the code in src/lagrangian/basic (in OF 1.5) I could find the F2F article again. He also mentions, that the mesh needs to be converted into a tetrahedral mesh. Otherwise calculating, whether the particle is in a cell or not, may fails. The solution to this problem is simply to decompose the hexahedral mesh into a tetrahedral mesh.
I played a bit with icoLagrangianFoam  and the cavityLagrangian case, where a tetrahedral mesh is generated using blockMesh. However when I look at the of code to calculate the lambdas (ParticleI.H ), which tells you whether a particle lies in a cell or not, I just find a reference to polyMesh. I cannot figure out where the mesh is decomposed into a hexahedral mesh.
It has been rewritten since then, you should look here instead.
Particle tracking in unstructured, arbitrary polyhedral meshes for CFD
and molecular dynamics,
G.B. Macpherson, N.Nordin and H.G Weller
Communications In Numerical Methods In Engineering (2008)
thx for the hint!
|All times are GMT -4. The time now is 18:16.|