Question regarding the F2F particle tracking algorithm implementation.

andrewryan · October 21, 2009, 17:47

Hi all

I have read Niklas Nordin's thesis [1] 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 [2] 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 [3]), 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.

Any ideas?

[1] http://powerlab.fsb.hr/ped/kturbo/Op...sNordinPhD.pdf
[2] http://openfoamwiki.net/index.php/Co...LagrangianFoam
[3] http://freefoam.sourceforge.net/doc/...8H_source.html

niklas · October 21, 2009, 18:18

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)

andrewryan · October 22, 2009, 02:55

thx for the hint!

October 21, 2009, 17:47	Question regarding the F2F particle tracking algorithm implementation.	#1
andrewryan Member Andrew Ryan Join Date: Mar 2009 Posts: 47 Rep Power: 17	Hi all I have read Niklas Nordin's thesis [1] 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 [2] 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 [3]), 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. Any ideas? [1] http://powerlab.fsb.hr/ped/kturbo/Op...sNordinPhD.pdf [2] http://openfoamwiki.net/index.php/Co...LagrangianFoam [3] http://freefoam.sourceforge.net/doc/...8H_source.html

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
[OpenFOAM] How to do particle tracking animation with Paraview?	John.I	ParaView	3	June 24, 2015 10:45
Steady state particle tracking?	John	CFX	5	August 11, 2013 19:34
restarting lagrange (particle tracking) simulation	dbdias	CFX	0	September 22, 2007 20:26
particle tracking collection tank	Neser	CFX	1	June 1, 2006 11:06
Particle tracking in s.s and u.s.s fluid field.	yueroo	FLUENT	0	April 6, 2001 23:02

October 21, 2009, 18:18		#2
niklas Super Moderator Niklas Nordin Join Date: Mar 2009 Location: Stockholm, Sweden Posts: 693 Rep Power: 29	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)

October 22, 2009, 02:55		#3
andrewryan Member Andrew Ryan Join Date: Mar 2009 Posts: 47 Rep Power: 17	thx for the hint!