[cjc96]

Hello everyone,

I'm currently working on a custom (PISO-based) Eulerian-Lagrangian Solver in OpenFOAM 7, primarily based on this presentation and existing work carried out by my department in OpenFOAM 2.4.

I'm trying to extract the data for post-processing in MATLAB, OpenFOAMs built-in 'particleTracks' function doesn't seem overly versatile, unless I am missing something here also.

My problem is reading the position data from:

*/lagrangian/kinematicCloud/positions

My data is of the format:

(0.08930289 0.6773368 0.1431805 0.09017981) 1135715 3437655 1

As you can see, there are seemingly four 'position' values

Whereas data I have from the work in 2.4 (albeit from a different problem) is formatted as so:

(-0.189055 -0.169927 0.157027) 22739688

Which has three 'position' values, as one would expect

Taking a rather crude measurement (using the ruler) of the corresponding particle position in ParaView shows an actual position of:

expectedPosition = [1.50022 0.031163 0.56261]

I've tried to manipulate the numbers based on an initial injection point of:

injectorPosition = [1.01, -0.15, 0.414]

But even then I can't get the numbers to add up.

If anyone has any experience in this regard I would be extremely grateful!