[danny123]

Hello,

I want to change the implementation of residual control in OpenFoam somewhat (interDymFoam is used). Right now, if you add additional PISO loops or Nonorthogonal corrections, the abstol is taken from the last iteration. I want to change this that it takes always the value of corrPISO == 1 (for 0 there is no pressure equation) and NonOrthogonalCorr == 0.

In order to do this, I thought to add an object in pimpleControlI.H:

inline bool Foam::pimpleControl::firstPISONOC() const
{
    return ((corrPISO_ == 1) && (corrNonOrtho_ == 0));
}

Now, this object can be called in PimpleControl.C similar to "storeIni". My problem is, that the object needs to be stored somewhere (abstol is local). I could do something similar like defining a new object like

residualControl_[fieldI].initialabstol

But where do I have to do this? I found in the C++ guide an object in solutionControl.H:

public:

    struct fieldData
    {
        wordRe name;
        scalar absTol;
        scalar relTol;
        scalar initialResidual;
    };

Can I just add the object to this?

Thanks for the support.

Regards,

Daniel