[Fairy]

Hello Foamers,
I hope you are all doing well.

I’m encountering a complex issue in OpenFOAM 5.x and would greatly appreciate your assistance. I have a class that calculates the temperature at a specific point in each cell, which is called every timestep.

The calculation of temperature at that specific point is conducted using cell, cellPoint, cellPintFace interpolation schemes.
The class for calling the interpolation header files is as follows:

Temp.C

interpolationCellPoint<scalar> voidfractionInterpolator_(voidfraction_);
interpolationCellPoint<vector> UInterpolator_(U_);
interpolationCellPoint<scalar> TInterpolator_(tempField_);

 for(int index = 0;index < particleCloud_.numberOfParticles(); ++index)
    {
            cellI = particleCloud_.cellIDs()[index][0];
            if(cellI >= 0)
            {
                    vector position = particleCloud_.position(index);
                    voidfraction = voidfractionInterpolator_.interpolate(position,cellI);
                    Ufluid = UInterpolator_.interpolate(position,cellI);
                    Tfluid = TInterpolator_.interpolate(position,cellI);
            }
   }

However, I’ve noticed that the line for:

interpolationCellPoint<scalar> TInterpolator_(tempField_);

suddenly gets stuck after a large number of iterations without any error. This behaviour occurs with the schemes cellPoint and cellPointFace, but not with cell, which is not accurate enough.



I’m wondering if this could be related to a memory leak for the TInterpolator_ because of its enourmous usage during simulation. The interpolationCellPoint.C file suggests a possible cause for this scenario:

template<class Type>
Foam::interpolationCellPoint<Type>::interpolationCellPoint
(
    const GeometricField<Type, fvPatchField, volMesh>& psi
)
:
    interpolation<Type>(psi),
    psip_
    (
        volPointInterpolation::New(psi.mesh()).interpolate
        (
            psi,
            "volPointInterpolate(" + psi.name() + ')',
            true        // use cache
        )
    )
{
    // Uses cellPointWeight to do interpolation which needs tet decomposition
    (void)psi.mesh().tetBasePtIs();
}

If it is indeed a memory issue, how can I reset the memory allocation for TInterpolator_ every timestep to prevent this problem?


Please guide me how to solve this issue!

[Fairy]

Nobody knows how I can reset the memory for his interpolation variable every time step?
I really got stuck into this issue.