[Cluap]

Hello everybody,

I'm currently simulating the injection of lagrangian particles in a water/air circuit : see the geometry file attached.

I use the coupling of the interfoam solver with lagrangian particles provided by swak4foam.


Here is my controlDict file :


libs (
"libOpenFOAM.so"
"libsimpleSwakFunctionObjects.so"
"libswakFunctionObjects.so"
"libgroovyBC.so"

"libsimpleFunctionObjects.so"
"libsimpleLagrangianFunctionObjects.so"
"libswakLagrangianCloudSourcesFunctionPlugin.s o"
);

application interFoam;

startFrom latestTime;

startTime 0;

stopAt endTime;

endTime 100;

deltaT 0.001;

writeControl adjustableRunTime;
//writeControl timeStep;

writeInterval 0.5;

purgeWrite 0;

writeFormat ascii;

writePrecision 6;

writeCompression uncompressed;

timeFormat general;

timePrecision 6;

runTimeModifiable yes;

adjustTimeStep yes;

maxCo 6;
maxAlphaCo 6;

maxDeltaT 1;


functions
{
rho
{
type expressionField;
expression "1000";
fieldName rho;
autowrite true;
}
mu
{
type expressionField;
expression "(1e-6)*rho";
fieldName mu;
autowrite true;
}
Particle
{
type evolveKinematicCloud;
cloudName ParticleCloud;
rhoName rho;
muName mu;
UName U;
}
}


And here my ParticleCloudProperties file :


solution
{
active on;
coupled false;
//coupled true;
transient true;
cellValueSourceCorrection on;

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

sourceTerms
{
schemes
{
rho explicit 1;
U explicit 1;
}
}

integrationSchemes
{
U Euler;
}
}

constantProperties
{
parcelTypeId 1;
rhoMin 1e-15;
minParticleMass 1e-15;
rho0 1.02;
epsilon0 10;
Pr 0.7;
pMin 0;
constantVolume 1;
youngsModulus 1e9;
poissonsRatio 0.35;
}

subModels
{
particleForces
{
gravity;
sphereDrag;
}

injectionModels
{
model1
{
type coneNozzleInjection;
SOI 0;
massTotal 58.3e-3;
parcelBasisType mass;
injectionMethod disc;
flowType constantVelocity;
//outerDiameter 0.0315;
outerDiameter 0.05;
innerDiameter 0;
duration 100;
position (-0.2 0 -0.05);
direction (1 0 0);
parcelsPerSecond 100000;
thetaInner constant 0.0;
thetaOuter constant 0.0;
sizeDistribution
{
type normal;
normalDistribution
{
minValue 1e-6;
maxValue 2e-6;
expectation 1.5e-6;
variance 0.5e-6;
}
}

flowRateProfile table
(
(0 1)
(1 1)
(2 1)
(3 1)
(100 1)
);
UMag 0;
}
}

dispersionModel none;

patchInteractionModel standardWallInteraction;

heatTransferModel none;
radiation off;
surfaceFilmModel none;
stochasticCollisionModel none;

standardWallInteractionCoeffs
{
type stick;
}
}


It seems to work well but then suddenly and randomly (it can occure after 5 or 10 or 20 seconds of simulation) the simulation crashes during the injection of the particles (before the crash both the injection and the solving of the particles is working very fine) with the message in the error.txt file attached.

Up to now, I haven'f found any solution to this problem, I have been working a lot on my mesh, on the injection position.... and I have no more idea, especially because this error appears randomly....

If anybody could help me on this issue it would be very appreciated !

Thanks in advance,

Paul