[mokard]

Hey Guys,
i actually have a problem extracting the lagrangian data out of an simulation with the sprayFoam solver.
For other i already use swak4Foam 0.4.4 and pyFoam 6.8.1 at Openfoam 5.x ........
Sampling lagrangian data
and
for sprayCloud
i wanted to do get the size distribution of the particles on a slice (e.g. i would say cuttingplane is the method here).
Also the velocity and so on is also from interest.
I checked out the swak4foam examples for lagrangian data extraction, but badly there is a problem with the type of sprayCloud which is used by the sprayFoam solver.
the controlDict with the functionobject looks like this with the libs.

libs (
"libsimpleSwakFunctionObjects.so"
"libswakLagrangianParser.so"
);

cloudPosition {
        type swakExpression;
        outputControlMode timeStep;
        expression "pos()";
        accumulations (
            min
            max
            average
            weightedAverage
        );
        verbose true;

        valueType cloud;
        cloudName sprayCloud;
    }
    cloudWeights {
        $cloudPosition;
        expression "weight()";
        //   alternateCloudType Cloud<basicKinematicParcel;
    }
         cloudVelocityDifference {
        $cloudPosition;
        expression "U-fluidPhase(U)";
        interpolationSchemes {
            U cellPoint;
        }

and the error message from the terminal comes here

swak4Foam: Setting default mesh
--> FOAM Warning : 
    From function bool Foam::functionObjectList::read()
    in file db/functionObjects/functionObjectList/functionObjectList.C at line 675
    Caught FatalError 
--> FOAM FATAL ERROR: 
Unknown  cloud type Cloud<basicSprayParcel> (cloud: Cloud<basicSprayParcel>) 
Consider setting the parameter 'alternateCloudType' in "/home/op50-ub14045-amd64/OpenFOAM/op50-ub14045-amd64-5.0/run/posttest_calculations_func_2_7/system/controlDict.functions.cloudPosition" to a valid type

Valid cloudType are :

9
(
Cloud<ReaderParticle>
Cloud<basicKinematicParcel>
Cloud<basicReactingMultiphaseParcel>
Cloud<basicReactingParcel>
Cloud<basicThermoParcel>
kinematicCloud<basicKinematicParcel>
kinematicCloud<basicReactingMultiphaseParcel>
kinematicCloud<basicReactingParcel>
kinematicCloud<basicThermoParcel>
)

The main problem is the cloud type which the sprayFoam solver uses.
This is the "basicSprayParcel" but this one can't be used with swak4Foam.
Does anyone know a possibility to extract the data from the cloud .... a solution without paraview seems to be better.
Another way of doing that is also welcome.

Thanks for an answer.
Greetings
Sebastian

[foamiste]

Hello, I am trying to do the same actually (extracting data concerning the lagrangian sprayCloud) but I dom't know how can I do it.

I tried to install swak4Foam for it doesn't compile correctly, how did you install it? knowing that I am using the development version of OpenFoam OF-dev.

I added your code to my controlDict, the simulation runs correctly without giving any errors and neither any results concerning what i am demanding, It acts like not accessing to the code added in controlDict.

Please if you find something concerning the data extraction make me know.

Best regards
Asmaa

[mokard]

Hey foamiste,

the functionobject code you see there is from
Bernhard Gschaiders examples in the example dir in swak4foam source.
Another user had a similar problem .... you can see the post here
Sampling lagrangian data
Badly all those automatic postprocessing was not successfull, but you can extract the particle data manually from paraview.
To extract the particel data you have to do this
Sampling lagrangian data

because i did not want to do something on the sourcecode i exported the lagrangian data manually from paraview.
After that i processed all the data by some scilab scripts i wrote in the recent days.
For example the diameter distribution in a plane and the average velocity of these particles and so on .....
As i could see you wanted to have all particles in an annulus -> a scilab script did this job for me.
Scilab is a free alternative to Matlab.
so i did a workaround in this case.
Greetings
Sebastian