[Ake]

Hi,

i am creating new fields in createFields.H dict which are the copy of a certain given field. For example the "wallShearStress" and a "new field".

volVectorField wallShearStress
 (
    IOobject
    (
       "wallShearStress",
        runTime.timeName(),
        mesh,
        IOobject::READ_IF_PRESENT,
        IOobject::AUTO_WRITE
    ),
    mesh,
 dimensionedVector("wallShearStress",dimensionSet( 0, 2, -2, 0, 0, 0, 0),Zero)
 );
volVectorField newField
(
    IOobject
    (
       "newField",
        runTime.timeName(),
        mesh,
        IOobject::READ_IF_PRESENT,
        IOobject::AUTO_WRITE
    ),
    mesh,
 dimensionedVector("newField",dimensionSet( 0, 2, -2, 0, 0, 0, 0),Zero)
);

In the controlDict i added the wallShearStress function, which calculates the shear stress in every time step.

    wallShearStress
    {
        type            wallShearStress;
        libs            ("libfieldFunctionObjects.so");
        executeControl  timeStep;
        writeControl    writeTime;
    }

In my solver.C file i added

newField=wallShearStress;

so the field "newField" is exactly a copy of the field "wallShearStress" in every time step. This is just an example, i actually do this to have/create fields which can be easily modified by me. This might be a "beginners solution" because i am not very good in OpenFoam programming but it works properly for me.


I tried to do exactly the same with a field which is created by a lagrangian post-processing cloud function "OwnPatchCollisionDensity" ("OwnPatchCollisionDensity" function is just a modified version of the standard "patchCollisionDensity" function). This function is added in the kinematicCloudProperties and it writes a boundary field which is called "collision" in my stored time folders.

cloudFunctions
{

    myCollisionFunction
    {
        type            OwnPatchCollisionDensity;
    }
  }

The function generally works properly. But, if i do the same procedure as described above (rename the field "wallShearStress" to "collision"), both fields, "collision" and "newField", remain empty. So the cloudFunction do not overwrite my created field as it does in the wallShearStress case.





1. Somebody has an idea why the fields "collision" and "newField" remain empty ? Maybe it is because the field "collision" is only a boundaryField? But isnt wallShearStress also only a boundaryField?


2. Maybe somebody has an other idea how to get access to the field of a lagrangian post-processing cloudFunction?




Thanks in advance!