[vishwesh]

I was following the work done by Aurelia on VOF and LPT conversion and in the move function described in the solidParticleCloud.C file, the interpolationCellPoint function is not working in parallel. It doesn't give any error but just doesn't proceed after that and gets stuck there. But when I try to run the code in serial, it works just fine. Can anyone help me out in finding a solution for this?

The move function:

void Foam::solidParticleCloud::move(const dimensionedVector& g, bool INJECT, vector posP, vector UP, scalar dP)
{
    Info<<"Function called"<<endl;
    const volScalarField& rho = mesh_.lookupObject<const volScalarField>("rho");
    Info<<"rho: "<<rho<<endl;
    const volVectorField& U = mesh_.lookupObject<const volVectorField>("U");
    const volScalarField& nu = mesh_.lookupObject<const volScalarField>("nu");
    const volScalarField& alphaW = mesh_.lookupObject<const volScalarField>("alpha.water");
    const volScalarField& p = mesh_.lookupObject<const volScalarField>("p");

    Info<<"lookup object worked"<<endl;

    Info<<"rho interpolation"<<endl;
    interpolationCellPoint<scalar> rhoInterp(rho);
//    Info<<"rho: "<<rhoInterp<<endl;
    Info<<"U interpolation"<<endl;
    interpolationCellPoint<vector> UInterp(U);
//    Info<<"U: "<<UInterp<<endl;
    Info<<"nu interpolation"<<endl;
    interpolationCellPoint<scalar> nuInterp(nu);
//    Info<<"nu: "<<nuInterp<<endl;
    Info<<"alpha.water interpolation"<<endl;
    interpolationCell<scalar> alphaWInterp(alphaW);
//    Info<<"alpha: "<<alphaWInterp<<endl;
    Info<<"p interpolation"<<endl;
    interpolationCellPoint<scalar> pInterp(p);
//    Info<<"p: "<<pInterp<<endl;

    Info<<"Interpolation complete"<<endl;

    smom_ = vector::zero;
    correctalphaW_=0;
    correctU_=vector::zero;

    Info<<"smom part complete"<<endl;
    solidParticle::trackingData
        td(*this, rhoInterp, UInterp, nuInterp, alphaWInterp, pInterp, g.value());
//      solidParticle::trackingData
//        td(*this,rho,U,nu,alphaW,p,g.value());
    Info<<"td initialized"<<endl;
    Info<<"Second move called"<<endl;
    Cloud<solidParticle>::move(td, mesh_.time().deltaTValue());
    Info<<"Second move completed"<<endl;
    // Inject
    if (INJECT){
    //    if(mesh_.time().value()> td.cloud().tInjStart_ && mesh_.time().value()< td.cloud().tInjEnd_)
      //  {
            Info<<"Injection called"<<endl;
            this->inject(td,posP,UP,dP);
       // }
    }
}

Result:

Function called
rho: dimensions      [1 -3 0 0 0 0 0];

internalField   uniform 1;

boundaryField
{
    inlet
    {
        type            symmetry;
    }
    outlet
    {
        type            symmetry;
    }
    walls
    {
        type            symmetry;
    }
    procBoundary0to1
    {
        type            processor;
        value           uniform 1;
    }
    procBoundary0to8
    {
        type            processor;
        value           uniform 1;
    }
    procBoundary0to16
    {
        type            processor;
        value           uniform 1;
    }
    procBoundary0to17
    {
        type            processor;
        value           uniform 1;
    }
    procBoundary0to24
    {
        type            processor;
        value           uniform 1;
    }
    procBoundary0to25
    {
        type            processor;
        value           uniform 1;
    }
}

lookup object worked
rho interpolation

Also, when I try to replace interpolationCellPoint with interpolationCell, it throws an error.

Thanks in advance

[vishwesh]

I couldn't resolve the interpolationCellPoint part but I revised them to interpolationCell in header files and the main C file and now it's working fine. Though I would grateful if someone could help me out with the interpolationCellPoint part.