[ali_atrian]

Hello
The following code is a part of the code written for a turbulent lift model (zanke lift model) which is used for lagrangian particles:

template<class CloudType>
Foam::vector Foam::ZankeLiftForce<CloudType>::n
(
    const label pID
) const
{
    dimensionedScalar Vsmall ("Vsmall",dimensionSet(0,-1,0,0,0),1E-08);
    volVectorField npV =fvc::grad(alphac_); 
    return npV[pID];    
}
    vector np=n(p.cell());
    value.Su() = fD *np*(1.0/2.5*pow(rnd*UrmsByUstar,2));

it compiles well but My solver can use this lift force only in series and not in parallel (with mpirun). In parallel run I receive these errors (even if all the particles are in a single processor)

[ubuntu:9557] *** An error occurred in MPI_Recv
[ubuntu:9557] *** on communicator MPI_COMM_WORLD
[ubuntu:9557] *** MPI_ERR_TRUNCATE: message truncated
[ubuntu:9557] *** MPI_ERRORS_ARE_FATAL: your MPI job will now abort
--------------------------------------------------------------------------
mpirun has exited due to process rank 1 with PID 9557 on
node ubuntu exiting improperly. There are two reasons this could occur:
1. this process did not call "init" before exiting, but others in
the job did. This can cause a job to hang indefinitely while it waits
for all processes to call "init". By rule, if one process calls "init",
then ALL processes must call "init" prior to termination.
2. this process called "init", but exited without calling "finalize".
By rule, all processes that call "init" MUST call "finalize" prior to
exiting or it will be considered an "abnormal termination"
This may have caused other processes in the application to be
terminated by signals sent by mpirun (as reported here).

I examined the problem and found that arises from the line:

volVectorField npV =fvc::grad(alphac_);

(although I can have :

volVectorField npV = alphac_* vector::one;

with no problem.)
So the main problem is because of fvc::grad. I don’t know the reason. Can any one tell me why alphac_ can not be used with fvc::grad !

[ali_atrian]

there is no answer to my problem yet ?!!
i investigated more :
the problem is due to this line in my code:

    vector np=n(p.cell());

this line should be executed as number as my particles (17000)
after passing the first particle this error appears:

[ubuntu:5427] *** An error occurred in MPI_Recv
[ubuntu:5427] *** on communicator MPI_COMM_WORLD
[ubuntu:5427] *** MPI_ERR_TRUNCATE: message truncated
[ubuntu:5427] *** MPI_ERRORS_ARE_FATAL: your MPI job will now abort

but my solver can not proceed any more for the second particle and i will receive the error:

mpirun has exited due to process rank 1 with PID 5427 on
node ubuntu exiting improperly. There are two reasons this could occur:

1. this process did not call "init" before exiting, but others in
the job did. This can cause a job to hang indefinitely while it waits
for all processes to call "init". By rule, if one process calls "init",
then ALL processes must call "init" prior to termination.

2. this process called "init", but exited without calling "finalize".
By rule, all processes that call "init" MUST call "finalize" prior to
exiting or it will be considered an "abnormal termination"

This may have caused other processes in the application to be
terminated by signals sent by mpirun (as reported here).

no one has no suggestion ?!!!
i can attach the whole code of my added lift force if neccessary

[Polli]

Hello,
i have exactly the same problem like you, do you sill found any solution for the mpi particle problem?