[vidyadhar]

Hi,

Can any one let me know how to add particle positions file in a case in openfoam.

Actually i am working on a case using the solver icoLagrangianFoam in foam-extend 3.2.

The kinematiccloud properties file contains particle injection position as shown below:

ConeInjectionCoeffs
{
SOI 0.2; //0
parcelBasisType mass;
duration 0.6;//7.5
position (0 0.0024 0.0024);
direction (1 0 0);
parcelsPerSecond 10000; //500
volumeFlowRate constant 0.01;
Umag constant 5.0;
thetaInner constant 0.0;
thetaOuter constant 30.0; //90.0;
parcelPDF
{
pdfType RosinRammler;
RosinRammlerPDF
{
minValue 50.0e-06;
maxValue 100.0e-06;
d (75.0e-06);
n (0.5);
}
}

// Total mass to inject
massTotal massTotal [ 1 0 0 0 0] 2.0e-4;
}

When I am working in series, there is no problem and execution is ok.

But, when i am trying to run the solver in parallel, I am getting error message that "Cannot read particle positions file" as shown below:


mpirun -np 4 icoLagrangianFoam -parallel
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 3.2 |
| \\ / A nd | Web: http://www.foam-extend.org |
| \\/ M anipulation | For copyright notice see file Copyright |
\*---------------------------------------------------------------------------*/
Build : 3.2-334ba0562a2c
Exec : icoLagrangianFoam -parallel
Date : Sep 23 2016
Time : 13:04:17
Host : mech5
PID : 9539
CtrlDict : "/home/mech5/Vidyadhar/work on foamextend3.2/Re1000withparticles_parallelprocess/system/controlDict"
Case : /home/mech5/Vidyadhar/work on foamextend3.2/Re1000withparticles_parallelprocess
nProcs : 4
Slaves :
3
(
iitmech5.9540
iitmech5.9541
iitmech5.9542
)

Pstream initialized with:
nProcsSimpleSum : 16
commsType : blocking
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading transportProperties

Reading field p

Reading field U

Reading/calculating face flux field phi


Reading g
Constructing kinematicCloud
--> FOAM Warning :
From function Cloud<ParticleType>::initCloud(const bool checkClass)
in file /home/mech5/foam/foam-extend-3.2/src/lagrangian/basic/lnInclude/CloudTemplateIO.C at line 124
Cannot read particle positions file
"/home/mech5/Vidyadhar/work on foamextend3.2/Re1000withparticles_parallelprocess/processor0/0/lagrangian/kinematicCloud"
assuming the initial cloud contains 0 particles.
Selecting DispersionModel none
Selecting DragModel SphereDrag
Selecting InjectionModel ConeInjection
Constructing 3-D injection
Selecting pdfType RosinRammler
[mech5:9541] *** An error occurred in MPI_Recv
[mech5:9541] *** on communicator MPI_COMM_WORLD
[mech5:9541] *** MPI_ERR_TRUNCATE: message truncated
[mech5:9541] *** MPI_ERRORS_ARE_FATAL: your MPI job will now abort
--------------------------------------------------------------------------
mpirun has exited due to process rank 2 with PID 9541 on
node iitmech5 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).
--------------------------------------------------------------------------
[mech5:09538] 1 more process has sent help message help-mpi-errors.txt / mpi_errors_are_fatal
[mech5:09538] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages

I want to work in parallel but, i have the following queries:

1. How to add this particle positions file, which was not existing in the original case?
OR
2. How to rectify this error without adding the particle positions file.

Can anyone let me know this.

Thank you