[oswald]

Dear Users,

I am trying to implement some forces in the lagrangian library. For this purpose I need the curl of the carrier velocity field at the particles' positions. So I included a interpolator similar for this in kinematicParcel, similar to the already existing interpolators.

    curlUInterp_
    (
        interpolation<vector>::New
        (
            cloud.solution().interpolationSchemes(),
            fvc::curl(cloud.U())
        )
    ),

This is in principle working. But after some time of running a simulation, I get an output like the following:

0] #0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #1 Foam::sigSegv::sigHandler(int) in "/opt/openfoam/openfoam211/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #2 in "/lib/libc.so.6"
[0] #3 Foam::interpolationCellPoint<Foam::Vector<double> >::interpolate(Foam::Vector<double> const&, Foam::tetIndices const&, int) const in "/opt/openfoam/openfoam211/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[0] #4
[0] at ~/OpenFOAM/user-2.1.1/src/lagrangian/incompressible/intermediate/lnInclude/KinematicParcel.C:67
...

Line 67 in KinematicParcel.C calls this new interpolator:

    curlUc_ = td.curlUInterp().interpolate(this->position(), tetIs);

All the other interpolators work fine. Does anyone have an idea what the error could be? And maybe how to solve it?

[dradenkovic]

Hello oswald,

Have you found the error?

Regards,
Darko

[oswald]

Hi Darko,

I don't remember exactly what the solution was, but at the moment I call the interpolator in KinematicParcel.C like this:

    curlUc_ = td.curlUcInterp().interpolate(this->position(), cellI);

and I changed the Cloud so that it has a field called curlUc and I can use this in KinematicParcelTrackingData.H:

    curlUcInterp_
    (
        interpolation<vector>::New
        (
            cloud.solution().interpolationSchemes(),
            cloud.curlUc()
        )
    ),

Hope this helps
oswald

[dradenkovic]

Thank you.

Regards,
Darko