|
[Sponsors] |
![]() |
![]() |
#1 |
New Member
Join Date: Sep 2019
Posts: 14
Rep Power: 6 ![]() |
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 |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Member
Simon
Join Date: Sep 2019
Location: Germany
Posts: 51
Rep Power: 6 ![]() |
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 In the folder: [lastTimestep]/uniform/lagrangian/kinematicCloud are two files after the simulation: cloudProperties: Code:
processor0 { particleCount 0; } Code:
injectionModel { model1 { massInjected 0; nInjections 245; parcelsAddedTotal 0; timeStep0 0; } } 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 // ************************************************************************* // Thanks Simon Last edited by SimonStar; May 19, 2020 at 04:33. Reason: Being more precise |
|
![]() |
![]() |
![]() |
![]() |
#3 |
Member
Simon
Join Date: Sep 2019
Location: Germany
Posts: 51
Rep Power: 6 ![]() |
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. |
|
![]() |
![]() |
![]() |
![]() |
#4 | |
Senior Member
Join Date: Jun 2020
Posts: 100
Rep Power: 6 ![]() |
Hey Simon
Quote:
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 |
||
![]() |
![]() |
![]() |
Tags |
icouncoupledkinematic, keyword, kinematiccloudproperties, no calculation |
|
|
![]() |
||||
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 |