CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Particle-tracking coupling with other cases.

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree3Likes
  • 1 Post By ybapat
  • 1 Post By ybapat
  • 1 Post By saidc.

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 24, 2020, 09:25
Post Particle-tracking coupling with other cases.
  #1
Member
 
cal
Join Date: Feb 2020
Location: nowhere
Posts: 65
Rep Power: 6
saidc. is on a distinguished road
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.
Attached Images
File Type: jpg 5-U.jpg (42.6 KB, 17 views)
File Type: jpg 5-UGlyph.jpg (52.7 KB, 17 views)
File Type: jpg 5-UMean.jpg (36.7 KB, 15 views)
Attached Files
File Type: zip My_Files.zip (7.1 KB, 4 views)
saidc. is offline   Reply With Quote

Old   March 25, 2020, 07:00
Default
  #2
Senior Member
 
Yogesh Bapat
Join Date: Oct 2010
Posts: 102
Rep Power: 16
ybapat is on a distinguished road
You can use icoUncoupledParcelFoam to solve particles with obtained flow field.
saidc. likes this.
ybapat is offline   Reply With Quote

Old   March 26, 2020, 12:47
Default
  #3
Member
 
cal
Join Date: Feb 2020
Location: nowhere
Posts: 65
Rep Power: 6
saidc. is on a distinguished road
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.
saidc. is offline   Reply With Quote

Old   March 27, 2020, 00:19
Default
  #4
Senior Member
 
Yogesh Bapat
Join Date: Oct 2010
Posts: 102
Rep Power: 16
ybapat is on a distinguished road
You can try different types of injections available.
saidc. likes this.
ybapat is offline   Reply With Quote

Old   March 27, 2020, 12:19
Post
  #5
Member
 
cal
Join Date: Feb 2020
Location: nowhere
Posts: 65
Rep Power: 6
saidc. is on a distinguished road
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. is offline   Reply With Quote

Old   March 28, 2020, 10:19
Post
  #6
Member
 
cal
Join Date: Feb 2020
Location: nowhere
Posts: 65
Rep Power: 6
saidc. is on a distinguished road
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.
Farid likes this.
saidc. is offline   Reply With Quote

Reply

Tags
les, particle tracking, piso, pitzdaily

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Particle tracking file per time step in an Unsteady Particle Tracking Manu4CFD FLUENT 1 December 29, 2022 08:15
Particle size and Mesh Size; Particle tracking; Suman Sapkota CFX 11 August 12, 2018 20:39
dispersion model with lagragian particle tracking model for incompressible flows eelcovv OpenFOAM Running, Solving & CFD 54 April 10, 2018 10:36
Ubuntu 12.10 + openfoam2.2.0 ==> paraview error message peteryuan OpenFOAM Installation 6 August 18, 2013 19:00
Particle tracking prob, urgent. sakurabogoda CFX 1 March 11, 2013 22:11


All times are GMT -4. The time now is 21:28.