CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Particle-tracking coupling with other cases. (https://www.cfd-online.com/Forums/openfoam-solving/225350-particle-tracking-coupling-other-cases.html)

saidc. March 24, 2020 08:25
Particle-tracking coupling with other cases.
 
4 Attachment(s)
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 March 25, 2020 06:00
You can use icoUncoupledParcelFoam to solve particles with obtained flow field.

saidc. March 26, 2020 11:47
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 March 26, 2020 23:19
You can try different types of injections available.

saidc. March 27, 2020 11:19
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
Code:
/*--------------------------------*- 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.
Code:
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. March 28, 2020 09:19
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

Code:
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.


All times are GMT -4. The time now is 22:25.