[johanning]

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

[SimonStar]

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:

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:

processor0
 {
    particleCount   0;
 }

and kinematicCloudOutputProperties

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.

/*--------------------------------*- 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

[SimonStar]

I have another thing - maybe a stupid question...

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

Thank you
Simon

[atul1018]

Hey Simon

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