CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Pre-Processing

kinematicCloudProperties: keywords undefined

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 7, 2019, 09:14
Default kinematicCloudProperties: keywords undefined
  #1
New Member
 
Join Date: Sep 2019
Posts: 14
Rep Power: 6
johanning is on a distinguished road
Hello together

i have a problem executing an icoUncoupledKinematicParcel solver. I get the following error message.



FOAM FATAL IO ERROR:
[0] keyword Reflector is undefined in dictionary "/home/constant/kinematicCloudProperties/subModels/pairCollisionCoeffs/wallLocalSpringSliderDashpotCoeffs


My kinematicCloudProperties looks like:


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

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

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

integrationSchemes
{
U Euler;
}
}

constantProperties
{
rho0 1000;
youngsModulus 6e8; //E-Modul
poissonsRatio 0.35;

constantVolume false;
}

subModels
{
particleForces
{
sphereDrag;
gravity;
}

injectionModels
{
model1
{
type patchFlowRateInjection;
patchName Inlet;
duration 0.01;
concentration 0.25;
parcelConcentration 0.25;
parcelsPerSecond 20000;
meanParticleDiameter 0.001;
massTotal 0;
parcelBasisType fixed;
nParticle 70000;
SOI 0;
//positionsFile "kinematicCloudPositions"; //
// U0 ( 0 0 0 ); //
sizeDistribution
{
type normal;
// type fixedValue; //
normalDistribution
{
minValue 0.0001;
expectation 0.001;
variance 0.5;
maxValue 0.002;
}
// fixedValueDistribution //
// { //
// value 0.001; //
// }
}
}
}

stochasticCollisionModel none;

dispersionModel none;

patchInteractionModel standardWallInteraction;

heatTransferModel none;

surfaceFilmModel none;

collisionModel pairCollision;

radiation off;

pairCollisionCoeffs
{
// Maximum possible particle diameter expected at any time
maxInteractionDistance 0.002;

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;
"wall.*"
{
youngsModulus 1e10;
poissonsRatio 0.23;
alpha 0.12;
b 1.5;
mu 0.43;
cohesionEnergyDensity 0;
}
};
}

standardWallInteractionCoeffs
{
type rebound;
}
}







Can someone help me?





Thank you! Regards Johanning
johanning is offline   Reply With Quote

Old   May 14, 2020, 13:24
Default
  #2
Member
 
SimonStar's Avatar
 
Simon
Join Date: Sep 2019
Location: Germany
Posts: 51
Rep Power: 6
SimonStar is on a distinguished road
Hi,

good question - do you already have an solution? If yes I am very interessted. If no maybe you can check all the brackets? The error message
came in my case often when a bracket has been wrong. I know thats not the perfect answer but at least...

I would like to ask a question back to the community because I have a quite similar problem. I dont have an error message but it seems like there were the particles were not/wrong simulated. The log is the following - everywhere is zero:

Code:
Time = 0.0005

Evolving kinematicCloud

Solving 3-D cloud kinematicCloud
Cloud: kinematicCloud
    Current number of parcels       = 0
    Current mass in system          = 0
    Linear momentum                 = (0 0 0)
   |Linear momentum|                = 0
    Linear kinetic energy           = 0
    model1:
        number of parcels added     = 0
        mass introduced             = 0
    Rotational kinetic energy       = 0

ExecutionTime = 8.17 s  ClockTime = 9 s

End
For me that sounds not so good even if the simulation is running because his appears in every iteration

In the folder: [lastTimestep]/uniform/lagrangian/kinematicCloud are two files after the simulation:

cloudProperties:
Code:
processor0
 {
    particleCount   0;
 }
and kinematicCloudOutputProperties
Code:
injectionModel
{
    model1
    {
        massInjected    0;
        nInjections     245;
        parcelsAddedTotal 0;
        timeStep0       0;
    }
 }
So how can there be 245 Injections when the number of particles is zero?

Here is the kinematicCloudProperties only if you need it - but if I change it it doesn't change the result.
Code:
/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  7
     \\/     M anipulation  | icoUncoupledKinematicParcelFoam
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      kinematicCloudProperties;
}
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solution
{
    active          true;    //Activate or deactivate the particle cloud
    coupled         false;    //Enable or disable phase coupling
    transient       yes;    //Transient or steady-state solution
    //Enable/disable correction of momentum transferred to the Eulerian phase
    cellValueSourceCorrection off;
    maxCo           0.5;    //Maximum Courant (CFL) Number

    sourceTerms
    {
        schemes
        {
        //explicit or semiImplicit, relaxCoeff for each of the fields
        U semiImplicit 1;
        }
    }

    interpolationSchemes
    /*Available schemes are: cell, cellPatchConstrained, cellPoint
    cellPointFace, cellPointWallModified, pointMVC */
    {
        rho             cell;
        U               cellPoint;
        mu              cell;
    }

    integrationSchemes
    {
        U               Euler;
    }
}

// Properties of the particles (discrete phase)
constantProperties 
{
    parcelTypeId 1;

    rhoMin          1e-15;
    minParcelMass   1e-15;
    //particle density (kg/m3)
    rho0            964;
    //constant of proportionality relating the deformation to the 
    //applied stress for your particles (N/m2)
    youngsModulus   6e8;
    //ratio of traverse to axial strain for your particles (dimensionless)
    poissonsRatio   0.35;
    //MPPCIFoam needs only rho0 and alphaMax
    alphaMax        0.99;
}

// Specify applicable forces on the particles
subModels
{
    // Kräfte auf Partikel
    particleForces
    {
        sphereDrag;
        /*ErgunWenYuDrag //no clue
    {
            alphac alpha.air;
        }*/
        gravity;
    /*pressureGradient //no clue
    {
            U U;
        }*/

    /*paramagnetic {
            magneticSusceptibility -6.3e-9; //m^3/kg for graphite
            HdotGradH U;
    }*/

        /*virtualMass  {
            Cvm 0.5;
        }*/
        
        /*nonInertialFrame  {
            linearAccelerationName  linearAc;
            linearAcceleration      10;
            angularVelocityName     angVelo;
            angularVelocity         5;
            angularAccelerationName angAcc;
            angularAcceleration     5;
        }*/
    }

    injectionModels
    {
    model1
        {
    /*type - options are: patchFlowRateInjection; patchInjection; 
    cellZoneInjection; coneInjection, coneNozzleInjection, 
    fieldActivatedInjection, inflationInjection, injectionModel,
    kinematicLookupTableInjection, manualInjection, noInjection  */
            type            patchFlowRateInjection;

    //If type = manualInjection; or others?:
        //positionsFile    "kinematicCloudPositions";
        //massTotal        0;             //Total mass to be injected (kg)


    // If type = patchFlowRateInjection; or patchInjection;:
        //Name of patch
            patchName            inlet;
        //Total duration of particle injection
            duration        0.0005;

        // parcel concentration in my fluid?
            parcelConcentration    0.5;  
        concentration    0.1;  
        //Injected parcels/particles per second        
            parcelsPerSecond     100000;
            meanParticleDiameter 0.00005;
        flowRateProfile constant 1;

        //---------------------------------------------------------------------
        // options of how to represent your particles. besides 'fixed' by 'number' or 'mass'
            parcelBasisType     mass; 

        // number of particles to be injected or Number of particles per parcel?
            nParticle            10000;

            //total mass to be injected (kg)
            massTotal           0.5;

        // start of injection (time in s) 
            SOI                 0;     

        //Initial parcel/particle velocity (U0) 
        U0                  (2 0 0);      
              
        // Particle size variation - normal distribution or specified sizes     
            sizeDistribution
            {
        /*
                type        fixedValue; 
                fixedValueDistribution
                {
                    value   0.006;             
                    value   0.01;
                    value   0.001;
                    value   0.015;
                    value   0.004;
                    value   0.009;
                } */

        /*
        type        normal;
              normalDistribution    
        {
                  expectation   650e-6;
                    variance    25e-6;
                    minValue    500e-6;
                    maxValue    800e-6;
              } */

        
        type RosinRammler;
              RosinRammlerDistribution  
        {
                    minValue 200e-6;
                    maxValue 300e-6;
                    d        250e-6;
                    n        3;
              }
            }
        }
    }

    //Sub-models for particle behaviour
    //Turbulent dispersion models (Discrete Random Walk, Gradient Dispersion)
    // D I F F U S I O N  ?
    dispersionModel         none;

    //standardWallInteraction;
    patchInteractionModel     none;

    stochasticCollisionModel     none;

    localInteractionCoeffs
    {
        patches
        (
            wall       
        {
                type rebound;       //rebound, stick or escape
                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;
    }

    heatTransferModel         none;

    // Surface film model for dripping and film interaction 
    surfaceFilmModel         none; //absorb, bounce, splash

    collisionModel pairCollision;

    pairCollisionCoeffs
    {
        // Maximum possible particle diameter expected at any time -> sizeDistribution
        maxInteractionDistance  0.0001; //100 Mikrometer

        writeReferredParticleCloud no;

        pairModel pairSpringSliderDashpot;

        pairSpringSliderDashpotCoeffs
        {
            useEquivalentSize   no;
            alpha               0.12; //function of coefficient of restitution
            b                   1.5;  //Spring power: b = 1 (linear) or b = 3/2 (Hertzian theory)
            mu                  0.52; //Coefficient of friction
            cohesionEnergyDensity 0;
            collisionResolutionSteps 12;
        };

        wallModel    wallLocalSpringSliderDashpot;

        wallLocalSpringSliderDashpotCoeffs
        {
            useEquivalentSize no;
            collisionResolutionSteps 12;
            tank
            {
                youngsModulus   1e10;
                poissonsRatio   0.23;
                alpha           0.12;
                b               1.5;
                mu              0.43;
                cohesionEnergyDensity 0;
            }
            inlet
            {
                youngsModulus   1e10;
                poissonsRatio   0.23;
                alpha           0.12;
                b               1.5;
                mu              0.1;
                cohesionEnergyDensity 0;
            }
            outlet
            {
                youngsModulus   1e10;
                poissonsRatio   0.23;
                alpha           0.12;
                b               1.5;
                mu              0.1;
                cohesionEnergyDensity 0;
            }
            pipe
            {
                youngsModulus   1e10;
                poissonsRatio   0.23;
                alpha           0.12;
                b               1.5;
                mu              0.1;
                cohesionEnergyDensity 0;
            }
        };
    }

    stochasticCollisionModel none; //Particle collision: spring, slider+dash-pot

    radiation off;
}

// Use cloudFunctions to record particle tracks, calculate particle erosion ...
cloudFunctions
{
/*particleErosion   {
    functionObjectLibs ("libcloudFunctionObjects.so");
    enabled       true;
    outputControl outputTime;
    log           true;
    valueOutput   true;
    
    //yield stress for aluminium = 11000000 Pa or 11 MPa
    p             11000000;
    
    //Ratio of the depth of contact to the depth of cut (default value = 2 )
    psi 2;
    
    //Ratio of vertical to horizontal force components (2 for angular abrassive grains)
    K 2;
    
    patches 
    (
        moving-wall
    );
  }
*/
}
// Source: https://www.cfd-online.com/Forums/openfoam-pre-processing/174534-kinematiccloudproperties-details.html  and others

// ************************************************************************* //
I am like allways thankful for any hint.
Thanks
Simon

Last edited by SimonStar; May 19, 2020 at 04:33. Reason: Being more precise
SimonStar is offline   Reply With Quote

Old   May 15, 2020, 02:41
Default
  #3
Member
 
SimonStar's Avatar
 
Simon
Join Date: Sep 2019
Location: Germany
Posts: 51
Rep Power: 6
SimonStar is on a distinguished road
I have another thing - maybe a stupid question...

Could someone explain me what the difference between "parcelConcentration" and "concentration" is?

Thank you
Simon

Last edited by SimonStar; May 15, 2020 at 04:33.
SimonStar is offline   Reply With Quote

Old   May 20, 2021, 12:05
Default
  #4
Senior Member
 
Join Date: Jun 2020
Posts: 100
Rep Power: 6
atul1018 is on a distinguished road
Hey Simon

Quote:
Could someone explain me what the difference between "parcelConcentration" and "concentration" is?

Are you able to find the answer of this? I have looked in many threads but i am not able to understand it properly and also no tutorials are available for patchFlowRateInjection model.


Thanks
Atul
atul1018 is offline   Reply With Quote

Reply

Tags
icouncoupledkinematic, keyword, kinematiccloudproperties, no calculation

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
Second Derivative Zero - Boundary Condition fu-ki-pa OpenFOAM 11 March 27, 2021 04:28
LEMOS InflowGenerator r_gordon OpenFOAM Running, Solving & CFD 103 December 18, 2018 00:58
OpenFOAM 1.6-ext git installation on Ubuntu 11.10 x64 Attesz OpenFOAM Installation 45 January 13, 2012 12:38
OpenFOAM on MinGW crosscompiler hosted on Linux allenzhao OpenFOAM Installation 127 January 30, 2009 19:08
G95 + CGNS Bruno Main CFD Forum 1 January 30, 2007 00:34


All times are GMT -4. The time now is 14:51.