# Surge tank emptying using interFoam- Stops when surface reaches outlet

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

 May 26, 2010, 10:46 Surge tank emptying using interFoam- Stops when surface reaches outlet #1 New Member   Ramnik Singh Join Date: Jun 2009 Location: Stuttgart, Germany Posts: 18 Rep Power: 13 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 outlet at bottom. 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. 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:perator++() 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:perator++() 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; writeControl adjustableRunTime; writeInterval 0.05; purgeWrite 0; writeFormat ascii; writePrecision 6; writeCompression uncompressed; timeFormat general; timePrecision 6; runTimeModifiable yes; adjustTimeStep on; maxCo 0.5; maxDeltaT 1; I need help.. Thanks in advance.. Ramnik

 March 14, 2012, 07:26 I need your help. #2 New Member   Seon Hye Join Date: Feb 2012 Location: South Korea Posts: 18 Rep Power: 10 Hi. I'm new to OpenFOAM and I met the same problem as you. Did you solve it? I use pimpleDyMFoam. Strangely, depending to the rotation speed, somtimes this error disappears in the same analysis model.

 Tags free surface, interfoam, outlet boundary, surge tank, vof