[saidc.]

Hi everyone,

Firstly, if this is an unnecessary question sorry for this. I did solve a case with using PISO algorithm similar to pitzDaily file. I did LES simulation now I want to take it to the next level with particle-tracking. You can see what I've done at the attachment. I just want to a few particle goes to the volume of my geometry with inlet (at the bottom) and flow to the outlet (at the top). I looked for a tutorial case like this but I couldn't find anything. Is there a way to couple my case with particle-tracking cases? I am open to any advice (or source).

By the way, I'm adding my files and some screenshots in case you want to know.

Best regards,
Said.

[ybapat]

You can use icoUncoupledParcelFoam to solve particles with obtained flow field.

[saidc.]

Hi ybapat,

I tried icoUncoupledParcelFoam to ran my own case and it worked just fine. As you see in the attachment I did enter the 11 particle's position with my hands. But I want to the particles to enter the system on their own, without entering their position. How can i do this? And I wonder, is there a way to see where these particles accumulate in the system?

Note: I tried with simpleFoam pitzDaily case, not LES.

Gif:
https://gfycat.com/filthyhardibadanmalimbe

Best regards,
Said.

[ybapat]

You can try different types of injections available.

[saidc.]

Hi again ybapat,

I ran the simulation by injection (patchInjection). But there is a problem. After a while, the cycle stops working before completion. Why is this happening? What should I do?


kinematicCloudProperties

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  6
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      kinematicCloudProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

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

    integrationSchemes
    {
        U               Euler;
    }
}

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

subModels
{
    particleForces
    {
        sphereDrag;
        gravity;
    }

    injectionModels
    {
        model1
        {
            type            patchInjection;
            parcelBasisType fixed;
            patchName       inlet;
            nParticle       1;
            SOI             0;
            U0              (20 0 0);
            parcelsPerSecond 5000;
            sizeDistribution
            {
                type        normal;
                normalDistribution
                {
                    expectation 650e-6;
                    variance 25e-6;
                    minValue 500e-6;
                    maxValue 800e-6;
                }
            }
            flowRateProfile constant 1;
            massTotal       0;
            duration 10;
        }
    }

    dispersionModel none;

    patchInteractionModel none;

    surfaceFilmModel none;

    stochasticCollisionModel none;

    collisionModel none;

    pairCollisionCoeffs
    {
        maxInteractionDistance  0.0001;

        writeReferredParticleCloud no;

        pairModel pairSpringSliderDashpot;

        pairSpringSliderDashpotCoeffs
        {
            useEquivalentSize   no;
            alpha               0.12;
            b                   1.5;
            mu                  0.52;
            cohesionEnergyDensity 0;
            collisionResolutionSteps 12;
        };

        wallModel    wallLocalSpringSliderDashpot;

        wallLocalSpringSliderDashpotCoeffs
        {
            useEquivalentSize no;
            collisionResolutionSteps 12;
            upperWall
            {
                youngsModulus   1e10;
                poissonsRatio   0.23;
                alpha           0.12;
                b               1.5;
                mu              0.43;
                cohesionEnergyDensity 0;
            }
            lowerWall
            {
                youngsModulus   1e10;
                poissonsRatio   0.23;
                alpha           0.12;
                b               1.5;
                mu              0.43;
                cohesionEnergyDensity 0;
            }
            frontAndBack
            {
                youngsModulus   1e10;
                poissonsRatio   0.23;
                alpha           0.12;
                b               1.5;
                mu              0.1;
                cohesionEnergyDensity 0;
            }
        };
    }
}


cloudFunctions
{}


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

It looks like this.

Time = 0.29

Evolving kinematicCloud

Solving 3-D cloud kinematicCloud

Cloud: kinematicCloud injector: model1
    Added 5 new parcels

Cloud: kinematicCloud
    Current number of parcels       = 1450
    Current mass in system          = 2.49819e-07
    Linear momentum                 = (1.93968e-07 4.25057e-08 5.0635e-12)
   |Linear momentum|                = 1.98571e-07
    Linear kinetic energy           = 1.09014e-07
    model1:
        number of parcels added     = 1450
        mass introduced             = 2.49819e-07
    Rotational kinetic energy       = 0

ExecutionTime = 260.48 s  ClockTime = 262 s

Time = 0.291

Evolving kinematicCloud

Solving 3-D cloud kinematicCloud

Cloud: kinematicCloud injector: model1
    Added 5 new parcels

Cloud: kinematicCloud
    Current number of parcels       = 1455
    Current mass in system          = 2.50729e-07
    Linear momentum                 = (1.9443e-07 4.27472e-08 1.03438e-12)
   |Linear momentum|                = 1.99073e-07
    Linear kinetic energy           = 1.09353e-07
    model1:
        number of parcels added     = 1455
        mass introduced             = 2.50729e-07
    Rotational kinetic energy       = 0

ExecutionTime = 260.49 s  ClockTime = 262 s

Time = 0.292

Evolving kinematicCloud

Solving 3-D cloud kinematicCloud

Cloud: kinematicCloud injector: model1
    Added 5 new parcels

Gifs:

https://gfycat.com/revolvingoldfashionedindigobunting
https://gfycat.com/poshfarhoatzin

[saidc.]

Hi everyone,

I think i did solve my problem. I just reduced the parcelsPerSecond and increased the duration.

Gif:

1) https://gfycat.com/energeticdecisivecurlew
2) https://gfycat.com/admirablemadaruanas

injectionModels
    {
        model1
        {
            type            patchInjection;
            parcelBasisType fixed;
            patchName       inlet;
            nParticle       1;
            SOI             0;
            U0              (10 0 0);
            parcelsPerSecond 1000;
            sizeDistribution
            {
                type        normal;
                normalDistribution
                {
                    expectation 650e-6;
                    variance 25e-6;
                    minValue 500e-6;
                    maxValue 800e-6;
                }
            }
            flowRateProfile constant 1;
            massTotal       0;
            duration 100;
        }
    }

Best regards,
Said.