[red.yxg]

Dear all:

Now I am using interFoam to simulate two-phase flow. And I want to write out private data in interfaceProperties.C like curvature. Because in createFields.H of interFoam there is an object ''interface'' of class ''interfaceProperties''. So I just add this in createFields.H :

volScalarField 	K
    (
        IOobject
        (
            "K",
            runTime.timeName(),
            mesh,
            IOobject::NO_READ,
            IOobject::AUTO_WRITE
        ),
        interface.K()
    );

After recompile interFoam and run a case, the value of K in different time step is the same. It seems that K is just calculated from alpha distribution of the first time step instead of the updated alpha. Is this because createFields.H is included before the time loop ?
So how to solve this problem ? Any hint is appreciated.

[alexeym]

Hi,

you should update your K after interface.correct() call to copy newly calculated values to your variable. You can do it with simple K = interface.K().

Also you can try to get K from object registry and modify its wirteOpts and path, but it'll be more complicated.

[red.yxg]

Hi, Alexey:

Thank you for your reply. It works! Yes, this is simple to implement than writing another .C like writeCellCenters.C.

[Mahmoud_aboukhedr]

Dear All,

I tried to follow what you said to plot K, but it failed to compile.
What I did is ,,
1) Added this piece of code to the creatPatch.h file in interfoam
Info<< "Reading field K\n" << endl;
volScalarField K
(
IOobject
(
"K",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
K = interface.K()
);

2) Then compiled

I got this error each time ..

createFields.H: In function ‘int main(int, char**)’:
createFields.H:40:12: error: ‘interface’ was not declared in this scope
In file included from myinterFoam.C:61:0:
/home/ofoam/OpenFOAM/OpenFOAM-2.3.x/src/finiteVolume/lnInclude/readTimeControls.H:38:8: warning: unused variable ‘maxDeltaT’ [-Wunused-variable]
make: *** [Make/linux64GccDPOpt/myinterFoam.o] Error 1

I tried to change the code to
Info<< "Reading field K\n" << endl;
volScalarField K
(
IOobject
(
"K",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh
);

it complie but nothing happens in the K file after the run


Any help will do ..

Thanks mahmoud