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

Lagrangian particles injection with interFoam and swak4foam

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 12, 2018, 11:37
Default Lagrangian particles injection with interFoam and swak4foam
  #1
New Member
 
Join Date: Sep 2014
Posts: 11
Rep Power: 8
Cluap is on a distinguished road
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
Attached Images
File Type: png geometry.png (18.3 KB, 50 views)
Attached Files
File Type: txt error.txt (6.3 KB, 26 views)
Cluap is offline   Reply With Quote

Reply

Tags
injection, interfoam, lagrangian, particles, swak4foam

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
Multiphase Solid Particle Tracking alexlupo OpenFOAM Running, Solving & CFD 112 January 30, 2021 22:29
InterDyMFoam+simpleFunctionObject Elham OpenFOAM Running, Solving & CFD 5 July 10, 2017 11:59
[swak4Foam] Validation of swak4Foam for lagrangian particle tracking nero235 OpenFOAM Community Contributions 13 March 28, 2017 05:45
[swak4Foam] and setFields for interFoam Saideep OpenFOAM Community Contributions 4 September 4, 2015 19:20
[swak4Foam] Interface contour - interFoam PrzemekPL OpenFOAM Community Contributions 4 January 30, 2013 14:09


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