[juergenstarck]

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]

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]

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]

Hi,

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

Regards,
Manu