CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Surge tank emptying using InterFAOM- Stops when surface reaches outlet (http://www.cfd-online.com/Forums/openfoam-solving/76470-surge-tank-emptying-using-interfaom-stops-when-surface-reaches-outlet.html)

 Ramnik May 26, 2010 10:42

Surge tank emptying using InterFAOM- Stops when surface reaches outlet

Hello,

I am a lil new to openFOAM. I am trying ti simulate the emptying process of a surge tank. Its like a tank with a atmospheric opening at top ant drain at bottom.
ves an err
Its initially filled with water. Now I am using RANS kE model for simulation. It runs fien till the surface touches the outlet then it gives an error.

file:///home/ramnik/OpenFOAM/run/SurgeTank/RASurgetank/post/final.png

MULES: Solving for alpha1
Liquid phase volume fraction = 0.670639 Min(alpha1) = 0 Max(alpha1) = 1.00568
MULES: Solving for alpha1
Liquid phase volume fraction = 0.670638 Min(alpha1) = 0 Max(alpha1) = 1.01137
MULES: Solving for alpha1
Liquid phase volume fraction = 0.670638 Min(alpha1) = 0 Max(alpha1) = 1.01706
MULES: Solving for alpha1
Liquid phase volume fraction = 0.670638 Min(alpha1) = 0 Max(alpha1) = 1.02275
DICPCG: Solving for p, Initial residual = 0.0992787, Final residual = 0.00352065, No Iterations 3
DICPCG: Solving for p, Initial residual = 0.0016073, Final residual = 7.41546e-05, No Iterations 21
DICPCG: Solving for p, Initial residual = 0.000173425, Final residual = 9.82588e-08, No Iterations 109
time step continuity errors : sum local = 0.282273, global = 0.00788512, cumulative = 0.00788587
DILUPBiCG: Solving for epsilon, Initial residual = 0.110028, Final residual = 2.53095e-10, No Iterations 6
bounding epsilon, min: -3.72537e+06 max: 1.27902e+20 average: 5.69326e+17
DILUPBiCG: Solving for k, Initial residual = 0.220076, Final residual = 3.50192e-10, No Iterations 3
bounding k, min: -2.36398e-06 max: 8.31733e+10 average: 1.65517e+08
ExecutionTime = 29.2 s ClockTime = 104 s

Courant Number mean: 13.1522 max: 23472.1
deltaT = 9.41571e-11
--> FOAM Warning :
From function Time::operator++()
in file db/Time/Time.C at line 682
Increased the timePrecision from 8 to 9 to distinguish between timeNames at time 5.85179
Time = 5.85178531

--> FOAM Warning :
From function Time::operator++()
in file db/Time/Time.C at line 682
Increased the timePrecision from 9 to 10 to distinguish between timeNames at time 5.85179
MULES: Solving for alpha1

I have specified my controlDict as

application interFoam;

startFrom startTime;

startTime 1;

stopAt endTime;

endTime 10;

deltaT 0.01;

writeInterval 0.05;

purgeWrite 0;

writeFormat ascii;

writePrecision 6;

writeCompression uncompressed;

timeFormat general;

timePrecision 6;

runTimeModifiable yes;

maxCo 0.5;

maxDeltaT 1;

I need help..