CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   sloshingTank2D3DoF: "Continuity error cannot be removed by adjusting the outflow." (http://www.cfd-online.com/Forums/openfoam-solving/64303-sloshingtank2d3dof-continuity-error-cannot-removed-adjusting-outflow.html)

juergenstarck May 6, 2009 15:50

sloshingTank2D3DoF: "Continuity error cannot be removed by adjusting the outflow."
 
I am asking for help to get the interDyMFoam/sloshingTank2D3DoF from OpenFOAM-1.5-dev working. Only the 0.05, 0.1 and 0.15 times are written, and then the solver aborts with the following end of the logfile:

[...]
Courant Number mean: 0.0345293 max: 0.21001 velocity magnitude: 10.9003
deltaT = 0.01
Time = 0.18

solidBodyMotionFunctions::SDA::transformation(): Time = 0.18 transformation: ((0 -0.0259902 -0.0127197) (0.998593 (0.0530292 0 0)))
Execution time for mesh.update() = 3.55076e-17 s
MULES: Solving for gamma
Liquid phase volume fraction = 0.337341 Min(gamma) = -1.10013e-52 Max(gamma) = 1
MULES: Solving for gamma
Liquid phase volume fraction = 0.337343 Min(gamma) = -1.08473e-52 Max(gamma) = 1
MULES: Solving for gamma
Liquid phase volume fraction = 0.337345 Min(gamma) = -1.06955e-52 Max(gamma) = 1
GAMG: Solving for pd, Initial residual = 0.025885, Final residual = 9.95525e-05, No Iterations 2
time step continuity errors : sum local = 1.10822e-05, global = -1.02059e-05, cumulative = -0.000310133
PCG: Solving for pd, Initial residual = 0.00105857, Final residual = 7.2366e-10, No Iterations 6
time step continuity errors : sum local = 1.03142e-05, global = -1.03142e-05, cumulative = -0.000320448
ExecutionTime = 0.99 s ClockTime = 1 s

Courant Number mean: 0.0339053 max: 0.206394 velocity magnitude: 10.7115
deltaT = 0.01
Time = 0.19

solidBodyMotionFunctions::SDA::transformation(): Time = 0.19 transformation: ((0 -0.0272624 -0.0141292) (0.998452 (0.0556206 0 0)))
Execution time for mesh.update() = 0 s
MULES: Solving for gamma
Liquid phase volume fraction = 0.337348 Min(gamma) = -1.05447e-52 Max(gamma) = 1
MULES: Solving for gamma
Liquid phase volume fraction = 0.33735 Min(gamma) = -1.0396e-52 Max(gamma) = 1
MULES: Solving for gamma
Liquid phase volume fraction = 0.337352 Min(gamma) = -1.02494e-52 Max(gamma) = 1


Continuity error cannot be removed by adjusting the outflow.
Please check the velocity boundary conditions and/or run potentialFoam to initialise the outflow.
Specified mass inflow : 170.15
Specified mass outflow : 170.15
Difference : 2.19495e-08
Adjustable mass outflow : 0


From function adjustPhi
(
surfaceScalarField& phi,
const volVectorField& U,
const volScalarField& p
)
in file cfdTools/general/adjustPhi/adjustPhi.C at line 128.

FOAM exiting

(end of logfile)

When I increase the closedDomainTol by several orders of magnitude from 1e-10 to 1e-2, I can run the simulation until t=1.6 before the same error occurs. How can the problem be fixed? The other sloshingTank 2D/3D and 3DoF/6DoF show the same problematic behaviour.

The platform is 32bit Linux. Maybe the tutorial had beed developed on 64bit? Would this make a difference in double precision accuracy?

regards,
Juergen

juergenstarck May 7, 2009 15:28

Hmm, if it is not so easy to get this problem fixed, I'd like to ask you to quickly run one of the interDyMFoam/sloshingTank* tutorials and tell me whether you get past the 0.15s output time. It only takes a few seconds of wall time until my simulation aborts.

Please state whether you run a 32/64 bit system and whether you have compiled with any nondefault configuration options.

regards,
Juergen

manu May 16, 2009 17:37

"Continuity error cannot be removed by adjusting the outflow."
 
Hi Juergen,

I experience exactly the same problem with the sloshingTank tutorials. I run OpenFOAM-1.5-dev on a 64bit-machine (Ubuntu). OpenFoam has been compiled with the default configuration options.

Have you been able to find out on your side what was going wrong?

Regards,
Manu

manu May 17, 2009 04:01

Hi,

It's possible to circumvent the problem by downloading the last release of OF-1.5-dev from svn repository (r1240).

Regards,
Manu


All times are GMT -4. The time now is 09:31.