|
[Sponsors] |
March 23, 2016, 04:18 |
Stuck in a Rut- interDyMFoam!
|
#1 |
New Member
Alpha Beta
Join Date: Mar 2016
Posts: 28
Rep Power: 10 |
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) 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 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 |
Senior Member
Join Date: Aug 2014
Location: Germany
Posts: 292
Rep Power: 14 |
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: 10 |
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 |
Senior Member
Jan
Join Date: Jul 2009
Location: Hamburg
Posts: 143
Rep Power: 20 |
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 |
|
March 23, 2016, 11:38 |
|
#5 | ||
New Member
Alpha Beta
Join Date: Mar 2016
Posts: 28
Rep Power: 10 |
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: 88
Rep Power: 12 |
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)? |
|
March 24, 2016, 04:30 |
|
#7 | |
New Member
Alpha Beta
Join Date: Mar 2016
Posts: 28
Rep Power: 10 |
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.72856E-05(=0.05/1341). Regards, kkk |
|
March 24, 2016, 08:13 |
|
#9 |
New Member
Alpha Beta
Join Date: Mar 2016
Posts: 28
Rep Power: 10 |
||
March 24, 2016, 08:23 |
|
#10 |
Senior Member
Join Date: Aug 2014
Location: Germany
Posts: 292
Rep Power: 14 |
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: 10 |
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 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 |
|
March 24, 2016, 14:33 |
|
#13 | ||
New Member
Alpha Beta
Join Date: Mar 2016
Posts: 28
Rep Power: 10 |
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: 10 |
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 | Search this Thread |
Display Modes | |
|
|
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 |