CFD Online Discussion Forums - Stuck in a Rut- interDyMFoam!

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)

- - Stuck in a Rut- interDyMFoam! (https://www.cfd-online.com/Forums/openfoam-solving/168486-stuck-rut-interdymfoam.html)

Stuck in a Rut- interDyMFoam!

Sloshing in a cylindrical tank(with baffles):
(linear oscillatory motion)
I was able to do everything up to the point of solving.
I ran the solver interDyMFoam in both parallel and in serial and it always stops at the same point.
I have searched a lot to solve this problem but I can't find the solution.
I think it has something to do with the fvSchmes file, I'm not sure.
The problem is I'm not good at the schemes part or coding.. so if you could kindly help me out it would be really great.
I've been stuck on this for 2 days now and I need to get out of this rut.
The error message I get after the solver stops is this

Code:

m > serial.log

#0  Foam::error::printStack(Foam::Ostream&) at ??:?

#1  Foam::sigFpe::sigHandler(int) at ??:?

#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"

#3  double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?

#4  Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?

#5  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?

#6  Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?

#7  ? at ??:?

#8  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"

#9  ? at ??:?

Floating point exception (core dumped)

Turbulence model set to laminar
used a .unv file (ideasunvtofoam) mesh seems to be fine.
checkMesh ran ok

The log file:

Code:

/*---------------------------------------------------------------------------*\

| =========                 |                                                 |

| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |

|  \\    /   O peration     | Version:  3.0.1                                 |

|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |

|    \\/     M anipulation  |                                                 |

\*---------------------------------------------------------------------------*/

Build  : 3.0.1-119cac7e8750

Exec   : interDyMFoam

Date   : Mar 23 2016

Time   : 07:34:51

Host   : "alvito-desktop"

PID    : 4904

Case   : /media/alvito/50B8ED83B8ED67C4/sloshingTank3D

nProcs : 1

sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

fileModificationChecking : Monitoring run-time modified files using timeStampMaster

allowSystemOperations : Allowing user-supplied system call operations



// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

Create time



Create mesh for time = 0



Selecting dynamicFvMesh solidBodyMotionFvMesh

Selecting solid-body motion function oscillatingLinearMotion

Applying solid body motion to entire mesh



PIMPLE: Operating solver in PISO mode



Reading field p_rgh



Reading field U



Reading/calculating face flux field phi



Reading transportProperties



Selecting incompressible transport model Newtonian

Selecting incompressible transport model Newtonian

Selecting turbulence model type laminar



Reading g



Reading hRef

Calculating field g.h



No MRF models present



No finite volume options present



GAMGPCG:  Solving for pcorr, Initial residual = 0, Final residual = 0, No Iterations 0

time step continuity errors : sum local = 0, global = 0, cumulative = 0

Reading/calculating face velocity Uf



Courant Number mean: 0 max: 0



Starting time loop



Courant Number mean: 0 max: 0

Interface Courant Number mean: 0 max: 0

deltaT = 0.05

Time = 0.05



PIMPLE: iteration 1

solidBodyMotionFunctions::oscillatingLinearMotion::transformation(): Time = 0.05 transformation: ((0 0 0.249974) (1 (0 0 0)))

Execution time for mesh.update() = 0.1 s

MULES: Solving for alpha.water

Phase-1 volume fraction = 0.502155  Min(alpha.water) = 0  Max(alpha.water) = 1

MULES: Solving for alpha.water

Phase-1 volume fraction = 0.502155  Min(alpha.water) = 0  Max(alpha.water) = 1

MULES: Solving for alpha.water

Phase-1 volume fraction = 0.502155  Min(alpha.water) = 0  Max(alpha.water) = 1

GAMG:  Solving for p_rgh, Initial residual = 1, Final residual = 0.0054145, No Iterations 5

time step continuity errors : sum local = 0.0363415, global = 3.7707e-18, cumulative = 3.7707e-18

GAMGPCG:  Solving for p_rgh, Initial residual = 0.320791, Final residual = 1.66912e-09, No Iterations 14

time step continuity errors : sum local = 2.54153e-08, global = -3.77621e-18, cumulative = -5.5112e-21

ExecutionTime = 8.53 s  ClockTime = 9 s



Courant Number mean: 9.53097 max: 1341

Interface Courant Number mean: 0 max: 0

deltaT = 1.86359e-05

Time = 0.0500186



PIMPLE: iteration 1

solidBodyMotionFunctions::oscillatingLinearMotion::transformation(): Time = 0.0500186 transformation: ((0 0 0.250067) (1 (0 0 0)))

Execution time for mesh.update() = 0.1 s

MULES: Solving for alpha.water

Phase-1 volume fraction = 0.502155  Min(alpha.water) = 0  Max(alpha.water) = 1

MULES: Solving for alpha.water

Phase-1 volume fraction = 0.502155  Min(alpha.water) = 0  Max(alpha.water) = 1

MULES: Solving for alpha.water

Phase-1 volume fraction = 0.502155  Min(alpha.water) = 0  Max(alpha.water) = 1

GAMG:  Solving for p_rgh, Initial residual = 0.991987, Final residual = 0.00864878, No Iterations 3

time step continuity errors : sum local = 1.26632e-06, global = -1.82539e-17, cumulative = -1.82594e-17

GAMGPCG:  Solving for p_rgh, Initial residual = 0.350105, Final residual = 1.49944e-09, No Iterations 9

time step continuity errors : sum local = 5.5288e-13, global = -1.8254e-17, cumulative = -3.65134e-17

ExecutionTime = 12.49 s  ClockTime = 13 s



Courant Number mean: 0.00334541 max: 0.395015

Interface Courant Number mean: 2.97421e-05 max: 0.0853419

deltaT = 2.09918e-05

Time = 0.0500396



...

...

...

.

.

.





PIMPLE: iteration 1

solidBodyMotionFunctions::oscillatingLinearMotion::transformation(): Time = 0.0507003 transformation: ((0 0 0.253475) (1 (0 0 0)))

Execution time for mesh.update() = 0.1 s

MULES: Solving for alpha.water

Phase-1 volume fraction = 0.502155  Min(alpha.water) = -1.33512e-22  Max(alpha.water) = 1

MULES: Solving for alpha.water

Phase-1 volume fraction = 0.502155  Min(alpha.water) = -1.3318e-22  Max(alpha.water) = 1

MULES: Solving for alpha.water

Phase-1 volume fraction = 0.502155  Min(alpha.water) = -1.32848e-22  Max(alpha.water) = 1

GAMG:  Solving for p_rgh, Initial residual = 0.956983, Final residual = 0.00452823, No Iterations 4

time step continuity errors : sum local = 4.83318e-07, global = -6.91323e-23, cumulative = 1.08717e-16

GAMGPCG:  Solving for p_rgh, Initial residual = 0.958908, Final residual = 8.78e-10, No Iterations 13

time step continuity errors : sum local = 1.21724e-13, global = -2.32278e-22, cumulative = 1.08717e-16

ExecutionTime = 5135.04 s  ClockTime = 5170 s



Courant Number mean: 9.05004e-05 max: 0.614533

Interface Courant Number mean: 7.31745e-08 max: 0.00101864

deltaT = 2.48734e-105

Time = 0.0507003



PIMPLE: iteration 1

solidBodyMotionFunctions::oscillatingLinearMotion::transformation(): Time = 0.0507003 transformation: ((0 0 0.253475) (1 (0 0 0)))

Execution time for mesh.update() = 0.1 s

MULES: Solving for alpha.water

Phase-1 volume fraction = 0.502155  Min(alpha.water) = -1.32518e-22  Max(alpha.water) = 1

MULES: Solving for alpha.water

Phase-1 volume fraction = 0.502155  Min(alpha.water) = -1.32187e-22  Max(alpha.water) = 1

MULES: Solving for alpha.water

Phase-1 volume fraction = 0.502155  Min(alpha.water) = -1.31856e-22  Max(alpha.water) = 1

GAMG:  Solving for p_rgh, Initial residual = 0.960139, Final residual = 0.00508782, No Iterations 4

time step continuity errors : sum local = 5.47764e-07, global = 4.32756e-22, cumulative = 1.08717e-16

GAMGPCG:  Solving for p_rgh, Initial residual = 0.956138, Final residual = 8.04641e-10, No Iterations 14

time step continuity errors : sum local = 1.08039e-13, global = 1.24837e-23, cumulative = 1.08717e-16

ExecutionTime = 5140.02 s  ClockTime = 5175 s



Courant Number mean: 9.14019e-05 max: 0.623251

Interface Courant Number mean: 7.45456e-08 max: 0.00101048

deltaT = 1.99546e-105

Time = 0.0507003



PIMPLE: iteration 1

solidBodyMotionFunctions::oscillatingLinearMotion::transformation(): Time = 0.0507003 transformation: ((0 0 0.253475) (1 (0 0 0)))

Execution time for mesh.update() = 0.1 s

MULES: Solving for alpha.water

Phase-1 volume fraction = 0.502155  Min(alpha.water) = -1.66064e-22  Max(alpha.water) = 1

MULES: Solving for alpha.water

Phase-1 volume fraction = 0.502155  Min(alpha.water) = -1.31193e-22  Max(alpha.water) = 1

MULES: Solving for alpha.water

Phase-1 volume fraction = 0.502155  Min(alpha.water) = -1.30861e-22  Max(alpha.water) = 1

GAMG:  Solving for p_rgh, Initial residual = 0.956953, Final residual = 0.00452787, No Iterations 4

time step continuity errors : sum local = 4.8281e-07, global = 2.26007e-23, cumulative = 1.08717e-16

GAMGPCG:  Solving for p_rgh, Initial residual = 0.958879, Final residual = 8.74106e-10, No Iterations 13

time step continuity errors : sum local = 1.21032e-13, global = -2.05206e-22, cumulative = 1.08717e-16

ExecutionTime = 5144.76 s  ClockTime = 5180 s



Courant Number mean: 9.03909e-05 max: 0.614533

Interface Courant Number mean: 7.30126e-08 max: 0.00101815

deltaT = 1.62356e-105

Time = 0.0507003



PIMPLE: iteration 1

solidBodyMotionFunctions::oscillatingLinearMotion::transformation(): Time = 0.0507003 transformation: ((0 0 0.253475) (1 (0 0 0)))

Execution time for mesh.update() = 0.1 s

MULES: Solving for alpha.water

Phase-1 volume fraction = 0.502155  Min(alpha.water) = -1.92573e-22  Max(alpha.water) = 1

MULES: Solving for alpha.water

Phase-1 volume fraction = 0.502155  Min(alpha.water) = -1.302e-22  Max(alpha.water) = 1

MULES: Solving for alpha.water

Phase-1 volume fraction = 0.502155  Min(alpha.water) = -1.2987e-22  Max(alpha.water) = 1

GAMG:  Solving for p_rgh, Initial residual = 0.960114, Final residual = 0.00508802, No Iterations 4

time step continuity errors : sum local = 5.47244e-07, global = -7.20103e-22, cumulative = 1.08716e-16

GAMGPCG:  Solving for p_rgh, Initial residual = 0.956105, Final residual = 8.81804e-10, No Iterations 14

time step continuity errors : sum local = 1.18219e-13, global = 3.23608e-23, cumulative = 1.08716e-16

ExecutionTime = 5149.86 s  ClockTime = 5185 s



Courant Number mean: 9.12917e-05 max: 0.62325

Interface Courant Number mean: 7.43825e-08 max: 0.00101

deltaT = 1.30249e-105

Time = 0.0507003



PIMPLE: iteration 1

solidBodyMotionFunctions::oscillatingLinearMotion::transformation(): Time = 0.0507003 transformation: ((0 0 0.253475) (1 (0 0 0)))

Execution time for mesh.update() = 0.1 s

MULES: Solving for alpha.water

Phase-1 volume fraction = 0.502155  Min(alpha.water) = -1.29538e-22  Max(alpha.water) = 1

MULES: Solving for alpha.water

Phase-1 volume fraction = 0.502155  Min(alpha.water) = -1.29207e-22  Max(alpha.water) = 1

MULES: Solving for alpha.water

Phase-1 volume fraction = 0.502155  Min(alpha.water) = -1.69442e-22  Max(alpha.water) = 1

GAMG:  Solving for p_rgh, Initial residual = 0.956922, Final residual = 0.00452752, No Iterations 4

time step continuity errors : sum local = 4.82306e-07, global = 5.58656e-23, cumulative = 1.08716e-16

GAMGPCG:  Solving for p_rgh, Initial residual = 0.958851, Final residual = 8.69625e-10, No Iterations 13

time step continuity errors : sum local = 1.2026e-13, global = 3.44824e-22, cumulative = 1.08717e-16

ExecutionTime = 5154.7 s  ClockTime = 5190 s



Courant Number mean: 9.02822e-05 max: 0.614534

Interface Courant Number mean: 7.2852e-08 max: 0.00101766

deltaT = 1.05974e-105

Time = 0.0507003



PIMPLE: iteration 1

solidBodyMotionFunctions::oscillatingLinearMotion::transformation(): Time = 0.0507003 transformation: ((0 0 0.253475) (1 (0 0 0)))

Execution time for mesh.update() = 0.1 s

MULES: Solving for alpha.water

Phase-1 volume fraction = 0.502155  Min(alpha.water) = -1.28545e-22  Max(alpha.water) = 1

MULES: Solving for alpha.water

Phase-1 volume fraction = 0.502155  Min(alpha.water) = -1.94065e-22  Max(alpha.water) = 1

MULES: Solving for alpha.water

Phase-1 volume fraction = 0.502155  Min(alpha.water) = -1.49905e-22  Max(alpha.water) = 1

GAMG:  Solving for p_rgh, Initial residual = 0.960088, Final residual = 0.00508821, No Iterations 4

time step continuity errors : sum local = 5.46729e-07, global = -3.3235e-22, cumulative = 1.08716e-16

I read somewhere that it is caused by dividing by 0.. but I can't find out where.. nor do I know any sort of coding to move forward.
Awaiting your reply,

sincerely
XoitX

Hi,

i usually get this error message when there is something wrong with boundary conditions. As an example:

When i set the volume fraction to zero at outlet for the water phase but initalize the water field around that boundary to be water with setFieldsDict this "dividing by zero" error occured. So check again your BC's and maybe even post them. Could be the error.

Good luck.

Hello Phoenix
My boundary conditions in polymesh file are:-

Quote:

2
(
body
{
type patch;
nFaces 11384;
startFace 234726;
}
walls
{
type wall;
nFaces 1140;
startFace 246110;
}
)

It's the same config for sloshing tack 3d in multiphase/interdymfoam tutorials

the u and p_rgh files are also similar:

U

Quote:

dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{
body
{
type movingWallVelocity;
value uniform (0 0 0);
}

walls
{
type movingWallVelocity;
value uniform (0 0 0);
}
}

p_rgh

Quote:

dimensions [1 -1 -2 0 0 0 0];

internalField uniform 0;

boundaryField
{
body
{
type fixedFluxPressure;
}

walls
{
type fixedFluxPressure;
}
}

Is this what you were referring to by boundary conditions?
If so is there any error?
Thank you for your patience:)

In first time step you already have an extreme max Courant number. Therefore your timestep drops. Before the simulation exits, you have a time step of E-105. How do you define your motion?

Best regards,
Jan

Greetings JNSN
Here is the dynamicmeshdict:

Quote:

dynamicFvMesh solidBodyMotionFvMesh;

solidBodyMotionFvMeshCoeffs
{
solidBodyMotionFunction oscillatingLinearMotion;
oscillatingLinearMotionCoeffs
{
amplitude (0 0 10); //moves 10 meters left and right in z dir'n
omega 0.5; // rads/sec
}
}

And here is the controldict file:

Quote:

application interDyMFoam;

startFrom latestTime;

startTime 0;

stopAt endTime;

endTime 15;

deltaT 0.05;

writeControl adjustableRunTime;

writeInterval 0.05;

purgeWrite 0;

writeFormat ascii;

writePrecision 6;

writeCompression compressed;

timeFormat general;

timePrecision 6;

runTimeModifiable yes;

adjustTimeStep yes;

maxCo 0.5;
maxAlphaCo 0.5;
maxDeltaT 1;

/*functions
{
probes
{
type probes;
functionObjectLibs ("libsampling.so");
outputControl timeStep;
outputInterval 1;
probeLocations
(
(0 9.95 19.77)
(0 -9.95 19.77)
);
fixedLocations false;
fields
(
p
);
}
}*/

I can see the Courant no. getting small.. but I thought only if courant number exceeds 1 it would be a problem. How do I solve the problem in this case?
Awaiting your reply
Regards XoitX

Your initial timestep deltaT is way too big. Reduce it to something like 0.0001 or even 0.00001. When the Courant number allows it your timestep will get bigger by itself.
Decrease your maxAlphaCo to at least 0.2.
Can you post the output of checkMesh? The interFoam family is extremely picky about the mesh quality, especially non-orthogonality.
Have you checked if your case runs without dynamic mesh (staticFvMesh)?

Dear TobM
Here is the checkMesh output:

Quote:

No I haven't tried static mesh as there is no change in alpha at all times except when being moved(dynamic).. So there really isn't anything to calculate. Or am I wrong? :O
However, I guess I need a powerful computer for this?
Awaiting your reply
Sincerely
XoitX

Hi

It seems that initial deltaT = 0.05 is the value of being set in controlDict.
How about trying to change the value of deltaT in controlDict.
Max courant number became 1341 when deltaT was 0.05. So i think it's better you set deltaT around 3.72856E-05(=0.05/1341).

Regards,
kkk

Hello kkk

Quote:

Originally Posted by kkk (Post 591426)

So i think it's better you set deltaT around 3.72856E-05(=0.05/1341).

!!!
Would't that take forever to calculate?!! :OOoooOO
Regards
XoitX

Hi,

can you please try to solve with interFoam solver? You have a pure tet mesh, not sure this can work with interDyMFoam.

And can you set in boundary conditions:

2
(
body
{
type wall;
nFaces 11384;
startFace 234726;
}
walls
{
type wall;
nFaces 1140;
startFace 246110;
}
)

Regards.

Edit: And try with a lower inital time step size. Sth like 0.0001 sec or so.

Greetings Phoenix
Just a heads up:
First few lines of interFoam:
It's still running though

Quote:

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 3.0.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 3.0.1-119cac7e8750
Exec : interFoam -parallel
Date : Mar 24 2016
Time : 15:40:58
Host : "alvito-desktop"
PID : 10387
Case : /media/alvito/50B8ED83B8ED67C4/Bryce_Project/sloshingTank3D
nProcs : 2
Slaves : 1("alvito-desktop.10388")
Pstream initialized with:
floatTransfer : 0
nProcsSimpleSum : 0
commsType : nonBlocking
polling iterations : 0
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

PIMPLE: Operating solver in PISO mode

Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Reading transportProperties

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type laminar

Reading g

Reading hRef
Calculating field g.h

No MRF models present

No finite volume options present

GAMGPCG: Solving for pcorr, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0, global = 0, cumulative = 0
Courant Number mean: 0 max: 0

Starting time loop

Courant Number mean: 0 max: 0
Interface Courant Number mean: 0 max: 0
deltaT = 0.000119904
Time = 0.000119904

PIMPLE: iteration 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1
GAMG: Solving for p_rgh, Initial residual = 1, Final residual = 0.00511966, No Iterations 2
time step continuity errors : sum local = 1.29239e-07, global = 2.11544e-23, cumulative = 2.11544e-23
GAMG: Solving for p_rgh, Initial residual = 0.109575, Final residual = 0.000787705, No Iterations 2
time step continuity errors : sum local = 3.9982e-08, global = 3.31833e-24, cumulative = 2.44727e-23
GAMG: Solving for p_rgh, Initial residual = 0.0210864, Final residual = 0.000195089, No Iterations 2
time step continuity errors : sum local = 9.93983e-09, global = -4.40716e-24, cumulative = 2.00655e-23
GAMGPCG: Solving for p_rgh, Initial residual = 0.00481213, Final residual = 6.21575e-10, No Iterations 8
time step continuity errors : sum local = 3.18062e-14, global = 3.6494e-24, cumulative = 2.37149e-23
ExecutionTime = 4.09 s ClockTime = 5 s

Courant Number mean: 8.81196e-06 max: 0.00637932
Interface Courant Number mean: 0 max: 0
deltaT = 0.000143747
Time = 0.000263651

PIMPLE: iteration 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1
GAMG: Solving for p_rgh, Initial residual = 0.0100374, Final residual = 2.79838e-05, No Iterations 3
time step continuity errors : sum local = 2.03401e-09, global = 1.1647e-23, cumulative = 3.5362e-23
GAMG: Solving for p_rgh, Initial residual = 0.00350694, Final residual = 2.15058e-05, No Iterations 2
time step continuity errors : sum local = 1.56032e-09, global = -1.00173e-23, cumulative = 2.53447e-23
GAMG: Solving for p_rgh, Initial residual = 0.00097833, Final residual = 7.74116e-06, No Iterations 2
time step continuity errors : sum local = 5.62299e-10, global = 1.65071e-23, cumulative = 4.18517e-23
GAMGPCG: Solving for p_rgh, Initial residual = 0.000413048, Final residual = 1.22852e-09, No Iterations 6
time step continuity errors : sum local = 8.91737e-14, global = 1.58678e-23, cumulative = 5.77195e-23
ExecutionTime = 7.41 s ClockTime = 8 s

Courant Number mean: 2.28006e-05 max: 0.0168241
Interface Courant Number mean: 0 max: 0
deltaT = 0.000172098
Time = 0.000435749

PIMPLE: iteration 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1
GAMG: Solving for p_rgh, Initial residual = 0.00972878, Final residual = 6.18735e-05, No Iterations 3
time step continuity errors : sum local = 6.33826e-09, global = 1.46233e-23, cumulative = 7.23428e-23
GAMG: Solving for p_rgh, Initial residual = 0.00422462, Final residual = 3.51462e-05, No Iterations 2
time step continuity errors : sum local = 3.60219e-09, global = 2.92093e-23, cumulative = 1.01552e-22
GAMG: Solving for p_rgh, Initial residual = 0.000880323, Final residual = 4.78069e-06, No Iterations 3
time step continuity errors : sum local = 4.95246e-10, global = -3.27012e-23, cumulative = 6.8851e-23
GAMGPCG: Solving for p_rgh, Initial residual = 0.000351567, Final residual = 4.12172e-10, No Iterations 7
time step continuity errors : sum local = 4.21981e-14, global = -8.50198e-25, cumulative = 6.80008e-23
ExecutionTime = 10.9 s ClockTime = 12 s

Courant Number mean: 4.36634e-05 max: 0.0329283
Interface Courant Number mean: 0 max: 0
deltaT = 0.000205661
Time = 0.00064141

PIMPLE: iteration 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1
GAMG: Solving for p_rgh, Initial residual = 0.00888878, Final residual = 6.77533e-05, No Iterations 3
time step continuity errors : sum local = 9.75243e-09, global = -1.47627e-23, cumulative = 5.32381e-23
GAMG: Solving for p_rgh, Initial residual = 0.00436175, Final residual = 3.95747e-05, No Iterations 2
time step continuity errors : sum local = 5.70282e-09, global = 8.34848e-23, cumulative = 1.36723e-22
GAMG: Solving for p_rgh, Initial residual = 0.00101417, Final residual = 6.6542e-06, No Iterations 3
time step continuity errors : sum local = 9.61811e-10, global = 3.19043e-24, cumulative = 1.39913e-22
GAMGPCG: Solving for p_rgh, Initial residual = 0.000491046, Final residual = 9.82267e-10, No Iterations 8
time step continuity errors : sum local = 1.46041e-13, global = -2.33797e-23, cumulative = 1.16534e-22
ExecutionTime = 14.28 s ClockTime = 15 s

Courant Number mean: 7.3881e-05 max: 0.0565653
Interface Courant Number mean: 0 max: 0
deltaT = 0.000245565
Time = 0.000886975

PIMPLE: iteration 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1
GAMG: Solving for p_rgh, Initial residual = 0.00713191, Final residual = 5.78411e-05, No Iterations 3
time step continuity errors : sum local = 1.1727e-08, global = -2.21931e-23, cumulative = 9.43404e-23
GAMG: Solving for p_rgh, Initial residual = 0.00385601, Final residual = 1.44462e-05, No Iterations 3
time step continuity errors : sum local = 2.9303e-09, global = 1.21584e-23, cumulative = 1.06499e-22
GAMG: Solving for p_rgh, Initial residual = 0.00124761, Final residual = 8.49988e-06, No Iterations 3
time step continuity errors : sum local = 1.72191e-09, global = -2.61221e-23, cumulative = 8.03768e-23
GAMGPCG: Solving for p_rgh, Initial residual = 0.00080284, Final residual = 1.3976e-09, No Iterations 8
time step continuity errors : sum local = 2.88e-13, global = -1.76416e-23, cumulative = 6.27352e-23
ExecutionTime = 17.52 s ClockTime = 19 s

Courant Number mean: 0.000117394 max: 0.0898264
Interface Courant Number mean: 0 max: 0
deltaT = 0.00029409
Time = 0.00118106

PIMPLE: iteration 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1
GAMG: Solving for p_rgh, Initial residual = 0.00586822, Final residual = 4.45309e-05, No Iterations 3
time step continuity errors : sum local = 1.2795e-08, global = 6.07239e-24, cumulative = 6.88076e-23
GAMG: Solving for p_rgh, Initial residual = 0.00346099, Final residual = 1.85211e-05, No Iterations 3
time step continuity errors : sum local = 5.32098e-09, global = 4.15688e-23, cumulative = 1.10376e-22
GAMG: Solving for p_rgh, Initial residual = 0.00182272, Final residual = 1.48574e-05, No Iterations 3
time step continuity errors : sum local = 4.27616e-09, global = 1.16329e-22, cumulative = 2.26706e-22
GAMGPCG: Solving for p_rgh, Initial residual = 0.00165449, Final residual = 6.69773e-10, No Iterations 8
time step continuity errors : sum local = 1.96504e-13, global = -3.93319e-23, cumulative = 1.87374e-22
ExecutionTime = 20.9 s ClockTime = 22 s

Courant Number mean: 0.000180825 max: 0.135475
Interface Courant Number mean: 5.96811e-07 max: 0.006106
deltaT = 0.000351215
Time = 0.00153228

PIMPLE: iteration 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1
GAMG: Solving for p_rgh, Initial residual = 0.00578503, Final residual = 4.33386e-05, No Iterations 3
time step continuity errors : sum local = 1.7517e-08, global = -9.87361e-23, cumulative = 8.86375e-23
GAMG: Solving for p_rgh, Initial residual = 0.00407819, Final residual = 2.87623e-05, No Iterations 3
time step continuity errors : sum local = 1.16215e-08, global = -1.52364e-22, cumulative = -6.37262e-23
GAMG: Solving for p_rgh, Initial residual = 0.00341846, Final residual = 3.02377e-05, No Iterations 3
time step continuity errors : sum local = 1.22393e-08, global = -1.6968e-22, cumulative = -2.33406e-22
GAMGPCG: Solving for p_rgh, Initial residual = 0.00368881, Final residual = 3.91778e-10, No Iterations 9
time step continuity errors : sum local = 1.59619e-13, global = -2.00196e-22, cumulative = -4.33602e-22
ExecutionTime = 24.54 s ClockTime = 26 s

Courant Number mean: 0.000271888 max: 0.195903
Interface Courant Number mean: 2.57549e-06 max: 0.00756495
deltaT = 0.000417825
Time = 0.00195011

PIMPLE: iteration 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1
GAMG: Solving for p_rgh, Initial residual = 0.00773477, Final residual = 6.39528e-05, No Iterations 3
time step continuity errors : sum local = 3.60454e-08, global = 1.46799e-22, cumulative = -2.86803e-22
GAMG: Solving for p_rgh, Initial residual = 0.00682457, Final residual = 5.72334e-05, No Iterations 3
time step continuity errors : sum local = 3.21907e-08, global = 5.24412e-23, cumulative = -2.34362e-22
GAMG: Solving for p_rgh, Initial residual = 0.00729836, Final residual = 6.74478e-05, No Iterations 3
time step continuity errors : sum local = 3.80747e-08, global = 7.10508e-23, cumulative = -1.63311e-22
GAMGPCG: Solving for p_rgh, Initial residual = 0.00847307, Final residual = 8.52216e-10, No Iterations 9
time step continuity errors : sum local = 4.84005e-13, global = 2.05198e-22, cumulative = 4.18877e-23
ExecutionTime = 28.14 s ClockTime = 29 s

Courant Number mean: 0.000401187 max: 0.286757
Interface Courant Number mean: 7.78858e-06 max: 0.00944518
deltaT = 0.000490305
Time = 0.00244041

PIMPLE: iteration 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1
GAMG: Solving for p_rgh, Initial residual = 0.013688, Final residual = 0.000126205, No Iterations 3
time step continuity errors : sum local = 9.63735e-08, global = 1.45868e-22, cumulative = 1.87756e-22
GAMG: Solving for p_rgh, Initial residual = 0.0140624, Final residual = 0.000127639, No Iterations 3
time step continuity errors : sum local = 9.73861e-08, global = 1.61451e-22, cumulative = 3.49207e-22
GAMG: Solving for p_rgh, Initial residual = 0.0163994, Final residual = 0.000154724, No Iterations 3
time step continuity errors : sum local = 1.1858e-07, global = -6.13579e-22, cumulative = -2.64372e-22
GAMGPCG: Solving for p_rgh, Initial residual = 0.0196176, Final residual = 1.66686e-09, No Iterations 9
time step continuity errors : sum local = 1.30425e-12, global = -2.5582e-22, cumulative = -5.20192e-22
ExecutionTime = 31.57 s ClockTime = 33 s

The case folder in case you're interested:
https://drive.google.com/file/d/0B2_...ew?usp=sharing

Hi

I tried your case with interDyMFoam and found that the velociry became so high locally.
I think you have to modify the mesh.
To consider moving wall velocity, I also think around 1e-5 deltaT is valid.
I wonder if there is no way except increasing parallel number.
Anyway I am very interested in your case.
I will try to do somehow.

Regards,
kkk

Dear kkk

Quote:

I wonder if there is no way except increasing parallel number.

This is why I almost gave up. I can't shift to 2D as it wouldn't account for the curvature in the 3rd dimension. :(

Quote:

Anyway I am very interested in your case.
I will try to do somehow.

This is only the beginning, but a crucial first step. :)
I will be optimizing the baffle plates to reduce the sloshing effect.
There is a long journey ahead, but first I need to setup the case. :O
I'm trying a lot of ways to solve this!!
Anyway thanks for the support :D
Regards
XoitX

Dear xoitx

I tried to modify mesh manually.
Then, calculation have done with hex mesh model.
There was no so high velocity locally.

I tried to do snappyHexMesh too. And now I am trying test running.
The progress of calculation is not bad.
It will take around 20 ~ 30 hours to run untill 15 sec with 60 cores.

I attache a case folder. You can create hex mesh model with it.
Please execute ./snappyHexAndOtherOperations then, it will do all things except running.
Sorry I couldn't upload stl model because of the file limitation.
You have to prepare it yourself or if you inform me some cloud url, then I will send you.

I hope this will help you.
Just give it try.

Regards,
kkk

Dear kkk
I have checked the files you gave me. Thank you so much for your time.
I had the .stl file from earlier which I'm using now. I tried to use snappyhexmesh before but I gave up as I had no control over it.I'm currently running the script you sent to the .stl file.
The pics you sent on top are modified versions?Did you remove the curved wall?
I can see why though.. It must have reduced a lot of computational time.. I guess.:P
I sent a private message to you.. there are some things I want to discuss about the case. :)
Thanks
Regards
XoitX