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

icoUncoupledKinematicParcelFoam problem with injection

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

Like Tree1Likes
  • 1 Post By Adamkat

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 14, 2017, 07:53
Cool icoUncoupledKinematicParcelFoam problem with injection
  #1
New Member
 
Adam Katz
Join Date: Feb 2017
Posts: 23
Rep Power: 9
Adamkat is on a distinguished road
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:
Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
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

Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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:

Code:
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
Adamkat is offline   Reply With Quote

Old   March 14, 2017, 13:26
Default
  #2
Member
 
Joshua
Join Date: Dec 2016
Location: St. Louis, Missouri
Posts: 91
Rep Power: 10
Joshua14 is on a distinguished road
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
Joshua14 is offline   Reply With Quote

Old   March 15, 2017, 05:17
Default
  #3
New Member
 
Adam Katz
Join Date: Feb 2017
Posts: 23
Rep Power: 9
Adamkat is on a distinguished road
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:
Code:
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:
Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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.
Attached Files
File Type: zip animation.ogv.zip (88.0 KB, 63 views)
sourav90 likes this.
Adamkat is offline   Reply With Quote

Reply

Tags
particle tracing, particle tracking

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
problem with Injection Mark_89 STAR-CCM+ 1 December 19, 2013 06:20
problem with DPM injection Hayder Mohammed Main CFD Forum 0 August 12, 2013 03:01
jou file command for the injection! asal FLUENT 0 June 4, 2013 05:15
Adiabatic and Rotating wall (Convection problem) ParodDav CFX 5 April 29, 2007 20:13
DPM injection problem akr FLUENT 5 November 6, 2006 11:34


All times are GMT -4. The time now is 00:17.