[cojua8]

Hello!,

I'm working on a lagrangian project in which I need to add a List member to kinematicParcel. This list is of type List<tuple2<scalar, vector>> :

KinematicParcel.H

...
// Basset
            //- storage of dt [s] and Ur [m/s]
            List<Tuple2<scalar, vector>> BassetList_;
...

In KinematicParcel.C, I have:

...
       if
        (
            p.age()
         <= td.cloud().constProps().dtB()*(td.cloud().constProps().NB() + 1)
        )
        {
            p.BassetList().resize(p.BassetList().size() + 1);
        }
        // rota la lista: el mas reciente se ubicara al principio
        inplaceRotateList(p.BassetList(), 1);

        p.BassetList().first().first() = dt;
        p.BassetList().first().second() = p.Uc() - p.U();
...

with NB = 2000 and dtB = 0.0001, BassetList has a max size of about 200000.

In file KinematicParcelIO.C, i have:

...
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //

template<class ParcelType>
Foam::string Foam::KinematicParcel<ParcelType>::propertyList_ =
    Foam::KinematicParcel<ParcelType>::propertyList();

template<class ParcelType>
const std::size_t Foam::KinematicParcel<ParcelType>::sizeofFields_
(
    offsetof(KinematicParcel<ParcelType>, BassetList_)
  - offsetof(KinematicParcel<ParcelType>, active_)
);


// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //

template<class ParcelType>
Foam::KinematicParcel<ParcelType>::KinematicParcel
(
    const polyMesh& mesh,
    Istream& is,
    bool readFields
)
:
    ParcelType(mesh, is, readFields),
    active_(false),
    typeId_(0),
    nParticle_(0.0),
    d_(0.0),
    dTarget_(0.0),
    U_(Zero),
    rho_(0.0),
    age_(0.0),
    tTurb_(0.0),
    UTurb_(Zero),
    BassetList_(),
    rhoc_(0.0),
    Uc_(Zero),
    muc_(0.0),
    curlUc_(Zero)
...

(I do not write BassetList to file, since it writes big sized files (200k lines per particle, for thousands of particles).

The problem is that when I run my case with just 50 particles, I get a memory error (running on cluster) after a second of simulation