[particlecount_0]

Hi,


I am having a problem with the icoUncoupledKinematicParcelFoam solver. I am using a simple 3d Geometry of a square with a side length of 0.2m. Inside the square is a small obstacle. The mesh has about 160000 points. The left side of the square is my inlet and the right side my outlet. For generating the U-field needed for the icoUncoupledKinematicParcelFoam solver I used the simpleFoam solver. I have used different velocity boundary conditions ranging from 0.02m/s to 0.0002m/s. Up until now everything works fine.


When I am running the icoUncoupledKinematicParcelFoam solver I get now error messages but the particlecount stays 0, even though I should have injected particles. All generated cloud_properties files read:

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "0.5/uniform/lagrangian/kinematicCloud";
    object      cloudProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

processor0
{
    particleCount   0;
}

console output:

sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0


Reading g
Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type laminar
Selecting laminar stress model Stokes
Constructing kinematicCloud kinematicCloud
Constructing particle forces
    none
Constructing cloud functions
    none
Constructing particle injection models
Creating injector: model1
Selecting injection model patchInjection
    Constructing 3-D injection
Selecting distribution model fixedValue
Selecting dispersion model none
Selecting patch interaction model none
Selecting stochastic collision model none
Selecting surface film model none
Selecting U integration scheme Euler
Reading non-inertial frame fields
Courant Number mean: 0.0317746 max: 0.213496

Starting time loop

Time = 0.5

Evolving kinematicCloud
ExecutionTime = 1.43 s  ClockTime = 2 s

Time = 1

Evolving kinematicCloud
ExecutionTime = 1.84 s  ClockTime = 2 s

Time = 1.5

Evolving kinematicCloud
ExecutionTime = 2.23 s  ClockTime = 2 s

Time = 2

Evolving kinematicCloud
ExecutionTime = 2.64 s  ClockTime = 3 s

Time = 2.5

Evolving kinematicCloud
ExecutionTime = 3.05 s  ClockTime = 3 s

Time = 3

Evolving kinematicCloud
ExecutionTime = 3.45 s  ClockTime = 4 s

Time = 3.5

Evolving kinematicCloud
ExecutionTime = 3.86 s  ClockTime = 4 s

Time = 4

Evolving kinematicCloud
ExecutionTime = 4.26 s  ClockTime = 4 s

Time = 4.5

Evolving kinematicCloud
ExecutionTime = 4.67 s  ClockTime = 5 s

Time = 5

Evolving kinematicCloud

Solving 3-D cloud kinematicCloud

Cloud: kinematicCloud injector: model1
    Added 10 new parcels

Cloud: kinematicCloud
    Current number of parcels       = 10
    Current mass in system          = 0.008
    Linear momentum                 = (0 0 0)
   |Linear momentum|                = 0
    Linear kinetic energy           = 0
    model1:
        number of parcels added     = 10
        mass introduced             = 0.008

 ExecutionTime = 5.23 s  ClockTime = 5 s

kinematicCloudProperties:

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      kinematicCloudProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solution
{
    active          true;
    coupled         false;
    transient       no;
    calcFrequency   10;   // calc pro cycle
    cellValueSourceCorrection off;
    maxCo           0.3;
    maxTrackTime 5;

    interpolationSchemes
    {
        rho             cell;
        U               cellPoint;
        mu              cell;
    }

    integrationSchemes
    {
        U               Euler;
    }
}

constantProperties
{
    rho0            1000;
    youngsModulus   1e8;
    poissonsRatio   0.35;
}

subModels
{
    particleForces
    {
        //sphereDrag;
        //gravity;
    }

    injectionModels
    {
        /*model1
        {
            type            manualInjection;
            massFlowRate        0.004583;
            massTotal       0;
            parcelBasisType fixed;
            nParticle       1;
            SOI             0;
            //positionsFile   "kinematicCloudPositions";
            U0              (0 0 0);
            sizeDistribution
            {
                type        fixedValue;
                fixedValueDistribution
                {
                    value   0.006;
                }
            }
        }*/
        
        model1
    {
        type patchInjection;
        parcelBasisType mass;
        patchName inlet;            
        U0 (0 0 0);
        //nParticle 10;    
        parcelsPerSecond 10;
        sizeDistribution
                {
                    type        fixedValue;
                    fixedValueDistribution
                    {
                        value   0.000001;
                    }
                }
                flowRateProfile constant 1;
                //massTotal 0.3;  //Choosing nParticle to be a fixed value, massTotal variable now             does not determine anything.
                SOI 0;
                duration 1;
                massFlowRate        0.8e-02;//0.0001;//0.004583;
        }
    }

    dispersionModel none;
    
    patchInteractionModel none;

    //patchInteractionModel standardWallInteraction;

    //standardWallInteractionCoeffs
    //{
      //      type        rebound; // stick, escape
     //       e           1;       // optional - elasticity coeff
    //        mu          0;       // optional - restitution coeff
    //};
    
    surfaceFilmModel none;

    stochasticCollisionModel none;

    collisionModel pairCollision;

controlDict:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  5                                     |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application     icoUncoupledKinematicParcelFoam;

startFrom       startTime;

startTime       0;

stopAt          endTime;

endTime         10;

deltaT          0.5;

writeControl    timeStep;

writeInterval   1;

purgeWrite      0;

writeFormat     ascii;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;

//functions
//{
//    #includeFunc streamlines
//}

// ************************************************************************* //

I would be very grateful if someone could point me in the right direction!

[HPE]

"console output" shows that 10 parcels were injected into the domain. therefore, there are definitely parcels in your domain unless they convected out of the domain. "particleCount" quantity seems to be taken from a processor file - maybe in the first processor data group, there is no parcel. maybe parcels reside in other processors. if the case is small, I would run it in serial with interim time-steps- just to observe it.