[Adamkat]

Hi all,
This is my first LPT and I can't seem to figure it out. It is quite frustrating, help will be much appreciated.

So far I ran simpleFoam on my model, then I run icoUncoupledKinematicParcelFoam but it doesn't time-loop (and doesn't create lagrangian particles). The output I get is:

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

Create mesh for time = 402


Reading g
Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type laminar
Constructing kinematicCloud kinematicCloud
Constructing particle forces
    Selecting particle force sphereDrag
Constructing cloud functions
    none
Constructing particle injection models
Creating injector: model1
Selecting injection model patchInjection
    Constructing 3-D injection
    Choosing nParticle to be a fixed value, massTotal variable now does not determine anything.
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
Selecting collision model none
Reading non-inertial frame fields

Starting time loop

End

I think it might be either the controlDict or the injectionModel so here they are:

controlDict

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

application     icoUncoupledKinematicParcelFoam;

startFrom       latestTime;

startTime       0;

stopAt          endTime;

endTime         1;

deltaT          1e-3;

writeControl    timeStep;

writeInterval   1;

purgeWrite      0;

writeFormat     binary;

writePrecision  8;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;

functions
{
    #includeFunc residuals
}
// ************************************************************************* //

And the injectionModel:

subModels
{
    particleForces
    {
        sphereDrag;
        
    }

    injectionModels
    {
    model1
        {
            type            patchInjection;
            parcelBasisType fixed;
            patchName       Inlet;
        nParticle       1;
            parcelsPerSecond 300000;
            massTotal  0;

            U0              (0 4.55 0);
            flowRateProfile constant 1;
        SOI 0;
            duration        1; // NOTE: set to 1 for steady state
        
        sizeDistribution
            {
                type        fixedValue;
                fixedValueDistribution
                {
                    value   15.5e-6;
                }
            }
       
            
    }
    }


    dispersionModel none;

    patchInteractionModel none;

    surfaceFilmModel none;

    stochasticCollisionModel none;

    collisionModel none;

    
}

Just to be super clear, I have a channel with water flowing in it and I wish to inject particles at the inlet at a constant rate (and size distribution).

Thanks a lot

[Joshua14]

Adamkat,

It does appear to be both files. Your controlDict says that the simulation ends at time = 1. When creating your mesh (and taking your most recent time) the simulation pulls time = 402 (from your log). As 402 > 1, the code reads this as your simulation is over.

Change your end time from 1 to 403, and your SOI for your injector from 0 to 402. This should have your Lagrangian solver run from 402 to 403.

Hope this helps,
Joshua

[Adamkat]

It worked!
Thanks Joshua!

I had to do some modifications to make it really work like adding wall rebound (see code below).

It is still missing one part, rotation. The particles must have rotation but for some reason they do not rotate.

What am I doing wrong?

I see this when I try to plot angular momentum (all zero) . It can also be seen in the runtime output:

Time = 402.00998

Evolving kinematicCloud

Solving 3-D cloud kinematicCloud
Cloud: kinematicCloud
    Current number of parcels       = 28
    Current mass in system          = 5.2629444e-11
    Linear momentum                 = (-8.1501672e-11 -1.4583216e-11 -8.5672026e-13)
   |Linear momentum|                = 8.2800524e-11
    Linear kinetic energy           = 2.8511479e-10
    model1:
        number of parcels added     = 297
        mass introduced             = 5.5824803e-10
    Parcel fate (number, mass)      : patch Walls
      - escape                      = 0, 0
      - stick                       = 0, 0
    Parcel fate (number, mass)      : patch Inlet
      - escape                      = 0, 0
      - stick                       = 0, 0
    Parcel fate (number, mass)      : patch Outlet
      - escape                      = 269, 5.0561818e-10
      - stick                       = 0, 0
    Rotational kinetic energy       = 0

ExecutionTime = 2.65 s  ClockTime = 3 s

I'm attaching the code from kinematicCloudProperties below:

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

solution
{
    active          true;
    coupled         false;
    transient       yes;
    cellValueSourceCorrection off;
    maxCo           0.3;

    sourceTerms
    {
        schemes
        {
        }
    }

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

    integrationSchemes
    {
        U               Euler;
    }
}


constantProperties
{
    rho0            964;
    youngsModulus   6e8;
    poissonsRatio   0.35;
}


subModels
{
    particleForces
    {
        sphereDrag;
    SaffmanMeiLiftForce
    {
     U U;
    }
        
    }

    injectionModels
    {
    model1
        {
            type            patchInjection;
            parcelBasisType fixed;
            patchName       Inlet;
        nParticle       1;
            parcelsPerSecond 30000;
            massTotal  0;

            U0              (0 4.55 0);
            flowRateProfile constant 1;
        SOI 402;
            duration        1; // NOTE: set to 1 for steady state
        
        sizeDistribution
            {
                type        fixedValue;
                fixedValueDistribution
                {
                    value   15.5e-6;
                }
            }
       
            
    }
    }


    dispersionModel none;

    patchInteractionModel none;

    surfaceFilmModel none;

    stochasticCollisionModel none;

    collisionModel none;
    
    patchInteractionModel localInteraction;

    localInteractionCoeffs
    {
        patches
        (
            Walls
            {
                type rebound;
                e    0.97;
                mu   0.09;
            }
            Inlet
        {
        type rebound;
                e    0.97;
                mu   0.09;
        }
        Outlet
        {
        type escape;
        }
        );
    }

    StandardWallInteractionCoeffs
    {
        type rebound;
        e    0.97;
        mu   0.09;
    }
}


cloudFunctions
{}


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

I also attached as a bonus a little animation of what I am talking about.