CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Stuck in a Rut- interDyMFoam!

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

Like Tree1Likes
  • 1 Post By kkk

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 23, 2016, 04:18
Default Stuck in a Rut- interDyMFoam!
  #1
New Member
 
Alpha Beta
Join Date: Mar 2016
Posts: 28
Rep Power: 10
xoitx is on a distinguished road
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

Last edited by wyldckat; March 26, 2016 at 06:50. Reason: Changed [QUOTE][/QUOTE] to [CODE][/CODE]
xoitx is offline   Reply With Quote

Old   March 23, 2016, 05:31
Default
  #2
Senior Member
 
Join Date: Aug 2014
Location: Germany
Posts: 292
Rep Power: 14
BlnPhoenix is on a distinguished road
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.
BlnPhoenix is offline   Reply With Quote

Old   March 23, 2016, 07:46
Default
  #3
New Member
 
Alpha Beta
Join Date: Mar 2016
Posts: 28
Rep Power: 10
xoitx is on a distinguished road
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
xoitx is offline   Reply With Quote

Old   March 23, 2016, 09:40
Default
  #4
Senior Member
 
JNSN's Avatar
 
Jan
Join Date: Jul 2009
Location: Hamburg
Posts: 143
Rep Power: 20
JNSN is on a distinguished road
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
JNSN is offline   Reply With Quote

Old   March 23, 2016, 11:38
Default
  #5
New Member
 
Alpha Beta
Join Date: Mar 2016
Posts: 28
Rep Power: 10
xoitx is on a distinguished road
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
xoitx is offline   Reply With Quote

Old   March 24, 2016, 03:34
Default
  #6
Member
 
Join Date: Sep 2014
Location: Germany
Posts: 88
Rep Power: 12
TobM is on a distinguished road
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)?
TobM is offline   Reply With Quote

Old   March 24, 2016, 04:30
Default
  #7
New Member
 
Alpha Beta
Join Date: Mar 2016
Posts: 28
Rep Power: 10
xoitx is on a distinguished road
Dear TobM
Here is the checkMesh output:
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 : checkMesh
Date : Mar 24 2016
Time : 11:24:15
Host : "lex"
PID : 14782
Case : /home/bryce/Scrap/BRYCE/170316/with_two_patches/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 polyMesh for time = 0

Time = 0

Mesh stats
points: 693636
faces: 7659568
internal faces: 7466056
cells: 3781406
faces per cell: 4
boundary patches: 2
point zones: 0
face zones: 0
cell zones: 0

Overall number of cells of each type:
hexahedra: 0
prisms: 0
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 3781406
polyhedra: 0

Checking topology...
Boundary definition OK.
Cell to face addressing OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces...
Patch Faces Points Surface topology
body 187031 93738 ok (non-closed singly connected)
walls 6481 3447 ok (non-closed singly connected)

Checking geometry...
Overall domain bounding box (-0.9 -0.9 -0.362048) (0.9 0.9 6.06205)
Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
Mesh has 3 solution (non-empty) directions (1 1 1)
Boundary openness (1.92732e-17 -1.79585e-15 -3.48972e-16) OK.
Max cell openness = 2.84739e-16 OK.
Max aspect ratio = 7.63881 OK.
Minimum face area = 1.5321e-05. Maximum face area = 0.00360165. Face area magnitudes OK.
Min volume = 2.90864e-08. Max volume = 5.85732e-05. Total volume = 15.7076. Cell volumes OK.
Mesh non-orthogonality Max: 61.1335 average: 14.476
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 0.872569 OK.
Coupled point location match (average 0) OK.

Mesh OK.

End
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
xoitx is offline   Reply With Quote

Old   March 24, 2016, 04:58
Default
  #8
kkk
New Member
 
kenyoh
Join Date: Apr 2012
Location: Tokyo Japan
Posts: 21
Blog Entries: 1
Rep Power: 14
kkk is on a distinguished road
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
kkk is offline   Reply With Quote

Old   March 24, 2016, 08:13
Default
  #9
New Member
 
Alpha Beta
Join Date: Mar 2016
Posts: 28
Rep Power: 10
xoitx is on a distinguished road
Hello kkk

Quote:
Originally Posted by kkk View Post
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
xoitx is offline   Reply With Quote

Old   March 24, 2016, 08:23
Default
  #10
Senior Member
 
Join Date: Aug 2014
Location: Germany
Posts: 292
Rep Power: 14
BlnPhoenix is on a distinguished road
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.
BlnPhoenix is offline   Reply With Quote

Old   March 24, 2016, 08:56
Default
  #11
New Member
 
Alpha Beta
Join Date: Mar 2016
Posts: 28
Rep Power: 10
xoitx is on a distinguished road
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
xoitx is offline   Reply With Quote

Old   March 24, 2016, 11:35
Default
  #12
kkk
New Member
 
kenyoh
Join Date: Apr 2012
Location: Tokyo Japan
Posts: 21
Blog Entries: 1
Rep Power: 14
kkk is on a distinguished road
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
Attached Images
File Type: png aaa.png (115.6 KB, 58 views)
kkk is offline   Reply With Quote

Old   March 24, 2016, 14:33
Default
  #13
New Member
 
Alpha Beta
Join Date: Mar 2016
Posts: 28
Rep Power: 10
xoitx is on a distinguished road
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
Regards
XoitX
xoitx is offline   Reply With Quote

Old   March 25, 2016, 04:34
Default
  #14
kkk
New Member
 
kenyoh
Join Date: Apr 2012
Location: Tokyo Japan
Posts: 21
Blog Entries: 1
Rep Power: 14
kkk is on a distinguished road
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
Attached Images
File Type: png remesh.png (121.2 KB, 52 views)
File Type: gif wallPressure.gif (179.0 KB, 48 views)
Attached Files
File Type: zip sloshingTank3D_kkk.zip (17.0 KB, 15 views)
xoitx likes this.
kkk is offline   Reply With Quote

Old   March 25, 2016, 08:09
Default
  #15
New Member
 
Alpha Beta
Join Date: Mar 2016
Posts: 28
Rep Power: 10
xoitx is on a distinguished road
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
xoitx is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
interDyMFoam with VOF + 6DOF instable pbalz OpenFOAM Running, Solving & CFD 11 October 9, 2020 06:19
modifying interDyMFoam for floatingObject Elisabeth_ofoam OpenFOAM Programming & Development 7 June 11, 2014 09:42
[ANSYS Meshing] "I shouldnt be disconnecting stuck triangles" sshah ANSYS Meshing & Geometry 0 March 15, 2012 14:33
interDyMFoam with GGI in 3d stawrogin OpenFOAM Running, Solving & CFD 2 January 5, 2011 03:17
error using interDyMFoam with kOmegaSST to simulate sloshing anmartin OpenFOAM Running, Solving & CFD 0 July 20, 2010 14:21


All times are GMT -4. The time now is 01:03.