
[Sponsors] 
March 23, 2016, 04:18 
Stuck in a Rut interDyMFoam!

#1 
New Member
Alpha Beta
Join Date: Mar 2016
Posts: 28
Rep Power: 2 
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_64linuxgnu/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_64linuxgnu/libc.so.6" #9 ? at ??:? Floating point exception (core dumped) 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.1119cac7e8750 Exec : interDyMFoam Date : Mar 23 2016 Time : 07:34:51 Host : "alvitodesktop" PID : 4904 Case : /media/alvito/50B8ED83B8ED67C4/sloshingTank3D nProcs : 1 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring runtime modified files using timeStampMaster allowSystemOperations : Allowing usersupplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 Selecting dynamicFvMesh solidBodyMotionFvMesh Selecting solidbody 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 Phase1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1 MULES: Solving for alpha.water Phase1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1 MULES: Solving for alpha.water Phase1 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.7707e18, cumulative = 3.7707e18 GAMGPCG: Solving for p_rgh, Initial residual = 0.320791, Final residual = 1.66912e09, No Iterations 14 time step continuity errors : sum local = 2.54153e08, global = 3.77621e18, cumulative = 5.5112e21 ExecutionTime = 8.53 s ClockTime = 9 s Courant Number mean: 9.53097 max: 1341 Interface Courant Number mean: 0 max: 0 deltaT = 1.86359e05 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 Phase1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1 MULES: Solving for alpha.water Phase1 volume fraction = 0.502155 Min(alpha.water) = 0 Max(alpha.water) = 1 MULES: Solving for alpha.water Phase1 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.26632e06, global = 1.82539e17, cumulative = 1.82594e17 GAMGPCG: Solving for p_rgh, Initial residual = 0.350105, Final residual = 1.49944e09, No Iterations 9 time step continuity errors : sum local = 5.5288e13, global = 1.8254e17, cumulative = 3.65134e17 ExecutionTime = 12.49 s ClockTime = 13 s Courant Number mean: 0.00334541 max: 0.395015 Interface Courant Number mean: 2.97421e05 max: 0.0853419 deltaT = 2.09918e05 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 Phase1 volume fraction = 0.502155 Min(alpha.water) = 1.33512e22 Max(alpha.water) = 1 MULES: Solving for alpha.water Phase1 volume fraction = 0.502155 Min(alpha.water) = 1.3318e22 Max(alpha.water) = 1 MULES: Solving for alpha.water Phase1 volume fraction = 0.502155 Min(alpha.water) = 1.32848e22 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.83318e07, global = 6.91323e23, cumulative = 1.08717e16 GAMGPCG: Solving for p_rgh, Initial residual = 0.958908, Final residual = 8.78e10, No Iterations 13 time step continuity errors : sum local = 1.21724e13, global = 2.32278e22, cumulative = 1.08717e16 ExecutionTime = 5135.04 s ClockTime = 5170 s Courant Number mean: 9.05004e05 max: 0.614533 Interface Courant Number mean: 7.31745e08 max: 0.00101864 deltaT = 2.48734e105 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 Phase1 volume fraction = 0.502155 Min(alpha.water) = 1.32518e22 Max(alpha.water) = 1 MULES: Solving for alpha.water Phase1 volume fraction = 0.502155 Min(alpha.water) = 1.32187e22 Max(alpha.water) = 1 MULES: Solving for alpha.water Phase1 volume fraction = 0.502155 Min(alpha.water) = 1.31856e22 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.47764e07, global = 4.32756e22, cumulative = 1.08717e16 GAMGPCG: Solving for p_rgh, Initial residual = 0.956138, Final residual = 8.04641e10, No Iterations 14 time step continuity errors : sum local = 1.08039e13, global = 1.24837e23, cumulative = 1.08717e16 ExecutionTime = 5140.02 s ClockTime = 5175 s Courant Number mean: 9.14019e05 max: 0.623251 Interface Courant Number mean: 7.45456e08 max: 0.00101048 deltaT = 1.99546e105 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 Phase1 volume fraction = 0.502155 Min(alpha.water) = 1.66064e22 Max(alpha.water) = 1 MULES: Solving for alpha.water Phase1 volume fraction = 0.502155 Min(alpha.water) = 1.31193e22 Max(alpha.water) = 1 MULES: Solving for alpha.water Phase1 volume fraction = 0.502155 Min(alpha.water) = 1.30861e22 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.8281e07, global = 2.26007e23, cumulative = 1.08717e16 GAMGPCG: Solving for p_rgh, Initial residual = 0.958879, Final residual = 8.74106e10, No Iterations 13 time step continuity errors : sum local = 1.21032e13, global = 2.05206e22, cumulative = 1.08717e16 ExecutionTime = 5144.76 s ClockTime = 5180 s Courant Number mean: 9.03909e05 max: 0.614533 Interface Courant Number mean: 7.30126e08 max: 0.00101815 deltaT = 1.62356e105 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 Phase1 volume fraction = 0.502155 Min(alpha.water) = 1.92573e22 Max(alpha.water) = 1 MULES: Solving for alpha.water Phase1 volume fraction = 0.502155 Min(alpha.water) = 1.302e22 Max(alpha.water) = 1 MULES: Solving for alpha.water Phase1 volume fraction = 0.502155 Min(alpha.water) = 1.2987e22 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.47244e07, global = 7.20103e22, cumulative = 1.08716e16 GAMGPCG: Solving for p_rgh, Initial residual = 0.956105, Final residual = 8.81804e10, No Iterations 14 time step continuity errors : sum local = 1.18219e13, global = 3.23608e23, cumulative = 1.08716e16 ExecutionTime = 5149.86 s ClockTime = 5185 s Courant Number mean: 9.12917e05 max: 0.62325 Interface Courant Number mean: 7.43825e08 max: 0.00101 deltaT = 1.30249e105 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 Phase1 volume fraction = 0.502155 Min(alpha.water) = 1.29538e22 Max(alpha.water) = 1 MULES: Solving for alpha.water Phase1 volume fraction = 0.502155 Min(alpha.water) = 1.29207e22 Max(alpha.water) = 1 MULES: Solving for alpha.water Phase1 volume fraction = 0.502155 Min(alpha.water) = 1.69442e22 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.82306e07, global = 5.58656e23, cumulative = 1.08716e16 GAMGPCG: Solving for p_rgh, Initial residual = 0.958851, Final residual = 8.69625e10, No Iterations 13 time step continuity errors : sum local = 1.2026e13, global = 3.44824e22, cumulative = 1.08717e16 ExecutionTime = 5154.7 s ClockTime = 5190 s Courant Number mean: 9.02822e05 max: 0.614534 Interface Courant Number mean: 7.2852e08 max: 0.00101766 deltaT = 1.05974e105 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 Phase1 volume fraction = 0.502155 Min(alpha.water) = 1.28545e22 Max(alpha.water) = 1 MULES: Solving for alpha.water Phase1 volume fraction = 0.502155 Min(alpha.water) = 1.94065e22 Max(alpha.water) = 1 MULES: Solving for alpha.water Phase1 volume fraction = 0.502155 Min(alpha.water) = 1.49905e22 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.46729e07, global = 3.3235e22, cumulative = 1.08716e16 Awaiting your reply, sincerely XoitX Last edited by wyldckat; March 26, 2016 at 06:50. Reason: Changed [QUOTE][/QUOTE] to [CODE][/CODE] 

March 23, 2016, 05:31 

#2 
Member
Join Date: Aug 2014
Location: Berlin
Posts: 58
Rep Power: 3 
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. 

March 23, 2016, 07:46 

#3  
New Member
Alpha Beta
Join Date: Mar 2016
Posts: 28
Rep Power: 2 
Hello Phoenix
My boundary conditions in polymesh file are: Quote:
the u and p_rgh files are also similar: U Quote:
Quote:
If so is there any error? Thank you for your patience 

March 23, 2016, 09:40 

#4 
Member
Jan
Join Date: Jul 2009
Location: Hamburg
Posts: 55
Rep Power: 12 
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 E105. How do you define your motion?
Best regards, Jan 

March 23, 2016, 11:38 

#5  
New Member
Alpha Beta
Join Date: Mar 2016
Posts: 28
Rep Power: 2 
Greetings JNSN
Here is the dynamicmeshdict: Quote:
Quote:
Awaiting your reply Regards XoitX 

March 24, 2016, 03:34 

#6 
Member
Join Date: Sep 2014
Location: Germany
Posts: 55
Rep Power: 3 
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 nonorthogonality. Have you checked if your case runs without dynamic mesh (staticFvMesh)? 

March 24, 2016, 04:30 

#7  
New Member
Alpha Beta
Join Date: Mar 2016
Posts: 28
Rep Power: 2 
Dear TobM
Here is the checkMesh output: Quote:
However, I guess I need a powerful computer for this? Awaiting your reply Sincerely XoitX 

March 24, 2016, 04:58 

#8 
New Member

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.72856E05(=0.05/1341). Regards, kkk 

March 24, 2016, 08:13 

#9 
New Member
Alpha Beta
Join Date: Mar 2016
Posts: 28
Rep Power: 2 

March 24, 2016, 08:23 

#10 
Member
Join Date: Aug 2014
Location: Berlin
Posts: 58
Rep Power: 3 
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. 

March 24, 2016, 08:56 

#11  
New Member
Alpha Beta
Join Date: Mar 2016
Posts: 28
Rep Power: 2 
Greetings Phoenix
Just a heads up: First few lines of interFoam: It's still running though Quote:
https://drive.google.com/file/d/0B2_...ew?usp=sharing 

March 24, 2016, 11:35 

#12 
New Member

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 1e5 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 

March 24, 2016, 14:33 

#13  
New Member
Alpha Beta
Join Date: Mar 2016
Posts: 28
Rep Power: 2 
Dear kkk
Quote:
Quote:
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 

March 25, 2016, 04:34 

#14 
New Member

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 

March 25, 2016, 08:09 

#15 
New Member
Alpha Beta
Join Date: Mar 2016
Posts: 28
Rep Power: 2 
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 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
interDyMFoam with VOF + 6DOF instable  pbalz  OpenFOAM Running, Solving & CFD  9  May 15, 2016 23:36 
modifying interDyMFoam for floatingObject  Elisabeth_ofoam  OpenFOAM Programming & Development  7  June 11, 2014 08: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 13:21 