[shuji_type97]

Hi Foamers,

Thanks for your reading.

I am still at the groping stage of using OpenFoam. I want to access the data "KD" calculated in the InterfaceCompositionPhaseChangePhaseSystem.C of MultiphaseEulerFoam because I want to get the Kla calculated in OpenFOAM and save it in the time folder.

May I ask whether is possible to access it by using mesh().lookupObject method? Actually, I tried that. Cause I know the Saturated Mass fraction, I want to calculate in ControlDict functions it by dividing mass transfer by rho like below:

code
       #{
            //Look up the dmdt
            const volScalarField& msp = mesh().lookupObject<volScalarField>("interfaceCompositionPhaseChange:dmidtfSp.O2.gasAndLiquid");
            
            //Look up the rho liquid
            const volScalarField& rho = mesh().lookupObject<volScalarField>("thermo:rho.liquid");
            
            volScalarField Kla
            (
               IOobject
               (
                  "Kla",
                  mesh().time().timeName(),
                  mesh(),
                  IOobject::NO_READ,
                  IOobject::AUTO_WRITE
               ),
               mesh(),
               dimensionedScalar("Kla", dimless/dimTime, 0.0)
            );
          
           Kla = msp/rho;
           
           if(mesh().time().outputTime())
           {
              Kla().write();
           }
           
       #};

No error was reported, but also no output in the time folder. I found in that Head file, dmidtfSp is private. Is this the reason why I cannot access it? I see my method is quite indirect. Is there any way to get the "KD" directly?

Any BIG BOSS could give me some suggestions about doing this?

Best wishes

[shuji_type97]

Hi Foamers,

Sorry for my stupidity. I make a fool mistake in my code.

in "functionsobjective" of controlDict, there is no block for "code". Now I use codeWrite, and it works.

Best wishes