[ssa_cfd]

Hi all,


I am trying to create a volScalarField for storing reaction rates kf and kr in my calculation. In my simulation case, the kf and kr variables are constant throughout the simulation.



I created a volScalarFields for the reaction rates as follows.

PtrList<volScalarField> kf(N);
Info << "Creating fields for kf" << endl;
for(int i = 0; i < N; i++)
{
    word kfi ("kf-un-named-reaction-" + name(i));
    kf.set
    (
        i,
        new volScalarField 
    (
            IOobject
            (
                kfi,
                runTime.timeName(),
                mesh,
                IOobject::NO_READ,
                IOobject::NO_WRITE
            ),
            mesh
,
           dimensionedScalar(kfi, dimless, 0.0)
        )
    );
}

PtrList<volScalarField> kr(N);
Info << "Creating fields for kr" << endl;
for(int i = 0; i < N; i++)
{
    word kri ("kr-un-named-reaction-" + name(i));
    kr.set
    (
        i,
        new volScalarField 
    (
            IOobject
            (
                kri,
                runTime.timeName(),
                mesh,
                IOobject::NO_READ,
                IOobject::NO_WRITE
            ),
            mesh
,
           dimensionedScalar(kri, dimless, 0.0)
        )
    );
    
}

But when I check, it creates volScalarField also for the Previous Iteration. It creates 4 volScalarFields for each reaction. (kf, kfPrevIter, kr, krPrevIter). Since the number of reactions is large, the simulation takes too much memory and slows down the simulation.



I want to make it constant such that it creates only two volScalarField (kf and kr).



Does anyone know how to do it.?


Thanks in advance,
ssa.

[adhiraj]

Would it help if you do not register the variables?

kf.set
(
    i,
    new volScalarField 
    (
        IOobject
        (
            kfi,
            runTime.timeName(),
            mesh,
            IOobject::NO_READ,
            IOobject::NO_WRITE,
            false
         ),
         mesh,
         dimensionedScalar(kfi, dimless, 0.0)
    )
);

[ssa_cfd]

If I do not register the fields, it deletes all the four (kf, kfPrevIter, kr, krPrevIter) volScalarFields.



I couldn't access kf and kr in the calculations.

[ssa_cfd]

I found a function called checkout() in the regIOobject.c ,


but it also deletes all the four fields.



Does someone know how to delete only the prevIter fields??