
[Sponsors] 
Cannot run the code properly: very large time step continuity error 

LinkBack  Thread Tools  Search this Thread  Display Modes 
January 17, 2014, 04:40 
Cannot run the code properly: very large time step continuity error

#1 
New Member
Steven
Join Date: Dec 2013
Location: Perth
Posts: 15
Rep Power: 12 
Hi all,
I am still very new to OpenFOAM. Currently I'm using OpenFOAM 2.2.2 and I've finished doing all of the tutorials given in the user guide. I want to model a basic data center configuration. My supervisor recommended me to base my model on the "hotRoom" tutorial using buoyantBoussinesqSimpleFoam as the solver. Initially, I had been able to run the simulation smoothly using just natural convection. Now, I incorporate some cool air raised from the floor in order to introduce forced convection cooling to the room. I met no problem when generating the mesh, but I got problems when I run the code. I have tried to fix it a lot of times, but still end up with the same error messages as follow. Code:
Time = 6 DILUPBiCG: Solving for Ux, Initial residual = 0.192382, Final residual = 0.00139781, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.328512, Final residual = 0.00357734, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 0.281499, Final residual = 0.00392744, No Iterations 1 DILUPBiCG: Solving for T, Initial residual = 0.000289576, Final residual = 6.58994e06, No Iterations 1 DICPCG: Solving for p_rgh, Initial residual = 0.999852, Final residual = 27153.3, No Iterations 1001 time step continuity errors : sum local = 1.90649e+33, global = 4.30234e+27, cumulative = 4.30429e+27 DILUPBiCG: Solving for epsilon, Initial residual = 0.0588009, Final residual = 2.80475e06, No Iterations 1 DILUPBiCG: Solving for k, Initial residual = 0.086828, Final residual = 0.00649997, No Iterations 1 ExecutionTime = 62.06 s ClockTime = 64 s Time = 7 DILUPBiCG: Solving for Ux, Initial residual = 0.000934178, Final residual = 7.98045e06, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.00368795, Final residual = 5.94455e05, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 0.00220201, Final residual = 1.0352e05, No Iterations 1 DILUPBiCG: Solving for T, Initial residual = 0.318863, Final residual = 0.0149065, No Iterations 1 DICPCG: Solving for p_rgh, Initial residual = 1, Final residual = 0.0094375, No Iterations 19 time step continuity errors : sum local = 2.49073e+89, global = 3.75377e+88, cumulative = 3.75377e+88 #0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #2 Uninterpreted: #3 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #4 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libfiniteVolume.so" #5 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam222/platforms/linuxGccDPOpt/bin/buoyantBoussinesqSimpleFoam" #6 Foam::fvMatrix<double>::solve() in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libincompressibleRASModels.so" #7 Foam::incompressible::RASModels::kEpsilon::correct() in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libincompressibleRASModels.so" #8 in "/opt/openfoam222/platforms/linuxGccDPOpt/bin/buoyantBoussinesqSimpleFoam" #9 __libc_start_main in "/lib/i386linuxgnu/libc.so.6" #10 in "/opt/openfoam222/platforms/linuxGccDPOpt/bin/buoyantBoussinesqSimpleFoam" Floating point exception (core dumped) Could anybody please help me?? I have attached several files which I think are the main cause of the error. Kind regards, Steven epsilon.txt k.txt p_rgh.txt U.txt fvSolution.txt Last edited by wyldckat; February 22, 2014 at 05:14. Reason: Added [CODE][/CODE] 

January 17, 2014, 05:22 
apologize

#2 
New Member
Steven
Join Date: Dec 2013
Location: Perth
Posts: 15
Rep Power: 12 
I forgot to clear all of the comments I put inside the files I've attached.
My apologize for that. Just ignore them . Cheers 

January 17, 2014, 05:35 

#3 
Senior Member
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,377
Rep Power: 30 
Do your inlets and outlets allow for the same mass flow? I.e. are you putting in the same flux as you are taking out?
__________________
*On twitter @akidTwit *Spend as much time formulating your questions as you expect people to spend on their answer. 

January 17, 2014, 07:45 

#4 
New Member
Steven
Join Date: Dec 2013
Location: Perth
Posts: 15
Rep Power: 12 
Hi Anton,
Thanks for your reply. Yes, I have employed mass continuity equation for the air velocity at inlets and outlets. I set the inlets and outlets to a particular temperature, so that I can determine the density variation and hence the velocity for each of them (the area for inlets and outlets are equal in magnitude). However, I rounded up the value up to 4 significant figures. Does it affect the result and causing error? Cheers 

February 19, 2014, 09:07 

#5  
New Member
Himanshu Kapoor
Join Date: Sep 2013
Posts: 4
Rep Power: 13 
Quote:
Did you manage to get a solution to your problem of high time step errors, I am also facing the same errors. Code:
Time = 1 smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.0422091, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.030476, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.0306028, No Iterations 2 DICPCG: Solving for p, Initial residual = 1, Final residual = 0.134351, No Iterations 100 DICPCG: Solving for p, Initial residual = 0.691177, Final residual = 0.0125443, No Iterations 100 DICPCG: Solving for p, Initial residual = 0.055563, Final residual = 0.00557511, No Iterations 100 DICPCG: Solving for p, Initial residual = 0.0220327, Final residual = 0.00114917, No Iterations 100 DICPCG: Solving for p, Initial residual = 0.0159711, Final residual = 0.00917452, No Iterations 100 DICPCG: Solving for p, Initial residual = 0.0187368, Final residual = 0.00188332, No Iterations 100 time step continuity errors : sum local = 1.54575, global = 0.135485, cumulative = 0.135485 ExecutionTime = 46.8 s ClockTime = 47 s Time = 2 smoothSolver: Solving for Ux, Initial residual = 0.451668, Final residual = 0.028682, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.751694, Final residual = 0.0484893, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.714429, Final residual = 0.0431807, No Iterations 2 DICPCG: Solving for p, Initial residual = 0.22116, Final residual = 0.00557813, No Iterations 100 DICPCG: Solving for p, Initial residual = 0.249683, Final residual = 0.00707944, No Iterations 100 DICPCG: Solving for p, Initial residual = 0.0734626, Final residual = 0.00226722, No Iterations 100 DICPCG: Solving for p, Initial residual = 0.024608, Final residual = 0.00154313, No Iterations 100 DICPCG: Solving for p, Initial residual = 0.0104501, Final residual = 0.00039303, No Iterations 100 DICPCG: Solving for p, Initial residual = 0.00496492, Final residual = 0.000300146, No Iterations 100 time step continuity errors : sum local = 0.075618, global = 0.00192783, cumulative = 0.133557 ExecutionTime = 82.56 s ClockTime = 83 s Time = 3 smoothSolver: Solving for Ux, Initial residual = 0.193993, Final residual = 0.0125072, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.292282, Final residual = 0.0178138, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.248495, Final residual = 0.0144195, No Iterations 2 DICPCG: Solving for p, Initial residual = 0.257503, Final residual = 0.00678077, No Iterations 100 DICPCG: Solving for p, Initial residual = 0.364373, Final residual = 0.00410846, No Iterations 100 DICPCG: Solving for p, Initial residual = 0.072356, Final residual = 0.00315433, No Iterations 100 DICPCG: Solving for p, Initial residual = 0.0276443, Final residual = 0.00105369, No Iterations 100 DICPCG: Solving for p, Initial residual = 0.0124894, Final residual = 0.000895709, No Iterations 100 DICPCG: Solving for p, Initial residual = 0.00933875, Final residual = 0.000367333, No Iterations 100 time step continuity errors : sum local = 0.121702, global = 0.0140299, cumulative = 0.147587 ExecutionTime = 118.29 s ClockTime = 118 s Time = 4 smoothSolver: Solving for Ux, Initial residual = 0.12122, Final residual = 0.00895779, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 0.141129, Final residual = 0.00892482, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 0.119265, Final residual = 0.00762635, No Iterations 2 DICPCG: Solving for p, Initial residual = 0.256797, Final residual = 0.00561705, No Iterations 100 DICPCG: Solving for p, Initial residual = 0.135723, Final residual = 0.000590494, No Iterations 100 DICPCG: Solving for p, Initial residual = 0.0450871, Final residual = 0.00114046, No Iterations 100 DICPCG: Solving for p, Initial residual = 0.00622627, Final residual = 0.000147822, No Iterations 100 DICPCG: Solving for p, Initial residual = 0.00502869, Final residual = 6.77996e05, No Iterations 100 DICPCG: Solving for p, Initial residual = 0.00374468, Final residual = 3.62683e05, No Iterations 100 time step continuity errors : sum local = 0.210205, global = 0.0591626, cumulative = 0.20675 ExecutionTime = 153.99 s ClockTime = 154 s Time = 5 smoothSolver: Solving for Ux, Initial residual = 0.269691, Final residual = 0.00476426, No Iterations 4 smoothSolver: Solving for Uy, Initial residual = 0.254113, Final residual = 0.00455208, No Iterations 4 smoothSolver: Solving for Uz, Initial residual = 0.25261, Final residual = 0.00443524, No Iterations 4 DICPCG: Solving for p, Initial residual = 0.0143989, Final residual = 0.000155821, No Iterations 100 DICPCG: Solving for p, Initial residual = 0.000232032, Final residual = 6.96373e07, No Iterations 100 DICPCG: Solving for p, Initial residual = 0.00017752, Final residual = 1.75389e06, No Iterations 100 DICPCG: Solving for p, Initial residual = 7.55826e05, Final residual = 6.34178e07, No Iterations 100 DICPCG: Solving for p, Initial residual = 8.81612e05, Final residual = 8.39305e07, No Iterations 100 DICPCG: Solving for p, Initial residual = 7.70637e05, Final residual = 7.11628e07, No Iterations 100 time step continuity errors : sum local = 2.85601, global = 0.77367, cumulative = 0.980419 ExecutionTime = 190.23 s ClockTime = 190 s Time = 6 smoothSolver: Solving for Ux, Initial residual = 0.29682, Final residual = 0.00460685, No Iterations 4 smoothSolver: Solving for Uy, Initial residual = 0.293445, Final residual = 0.00429992, No Iterations 4 smoothSolver: Solving for Uz, Initial residual = 0.299224, Final residual = 0.00440848, No Iterations 4 DICPCG: Solving for p, Initial residual = 9.35095e05, Final residual = 1.09954e06, No Iterations 100 DICPCG: Solving for p, Initial residual = 4.91151e06, Final residual = 4.30099e08, No Iterations 100 DICPCG: Solving for p, Initial residual = 8.06364e06, Final residual = 1.08889e07, No Iterations 100 DICPCG: Solving for p, Initial residual = 4.85845e06, Final residual = 6.03338e08, No Iterations 100 DICPCG: Solving for p, Initial residual = 5.74193e06, Final residual = 7.29363e08, No Iterations 100 DICPCG: Solving for p, Initial residual = 5.16892e06, Final residual = 6.34448e08, No Iterations 100 time step continuity errors : sum local = 77.6375, global = 23.7415, cumulative = 24.7219 ExecutionTime = 226.48 s ClockTime = 227 s Time = 7 smoothSolver: Solving for Ux, Initial residual = 0.318025, Final residual = 0.00485944, No Iterations 4 smoothSolver: Solving for Uy, Initial residual = 0.319293, Final residual = 0.00491225, No Iterations 4 smoothSolver: Solving for Uz, Initial residual = 0.321538, Final residual = 0.0049102, No Iterations 4 DICPCG: Solving for p, Initial residual = 5.87032e06, Final residual = 6.7756e08, No Iterations 100 DICPCG: Solving for p, Initial residual = 2.36787e07, Final residual = 1.64937e09, No Iterations 100 DICPCG: Solving for p, Initial residual = 4.45469e07, Final residual = 4.81168e09, No Iterations 100 DICPCG: Solving for p, Initial residual = 2.28813e07, Final residual = 2.19573e09, No Iterations 100 DICPCG: Solving for p, Initial residual = 2.69278e07, Final residual = 2.80125e09, No Iterations 100 DICPCG: Solving for p, Initial residual = 2.37724e07, Final residual = 2.25178e09, No Iterations 100 time step continuity errors : sum local = 1284.9, global = 278.963, cumulative = 303.685 ExecutionTime = 262.78 s ClockTime = 263 s Time = 8 smoothSolver: Solving for Ux, Initial residual = 0.311316, Final residual = 0.00514459, No Iterations 4 smoothSolver: Solving for Uy, Initial residual = 0.315939, Final residual = 0.00511636, No Iterations 4 smoothSolver: Solving for Uz, Initial residual = 0.320675, Final residual = 0.00517196, No Iterations 4 DICPCG: Solving for p, Initial residual = 3.33222e07, Final residual = 6.01642e09, No Iterations 100 DICPCG: Solving for p, Initial residual = 1.04966e08, Final residual = 2.08452e10, No Iterations 100 DICPCG: Solving for p, Initial residual = 1.98828e08, Final residual = 5.56181e10, No Iterations 100 DICPCG: Solving for p, Initial residual = 1.1944e08, Final residual = 3.33406e10, No Iterations 100 DICPCG: Solving for p, Initial residual = 1.39477e08, Final residual = 3.59735e10, No Iterations 100 DICPCG: Solving for p, Initial residual = 1.24737e08, Final residual = 2.85011e10, No Iterations 100 time step continuity errors : sum local = 67282.6, global = 12498.2, cumulative = 12801.9 ExecutionTime = 298.53 s ClockTime = 299 s Time = 9 smoothSolver: Solving for Ux, Initial residual = 0.321623, Final residual = 0.00491377, No Iterations 4 smoothSolver: Solving for Uy, Initial residual = 0.324749, Final residual = 0.00493128, No Iterations 4 smoothSolver: Solving for Uz, Initial residual = 0.328151, Final residual = 0.0049728, No Iterations 4 DICPCG: Solving for p, Initial residual = 1.74395e08, Final residual = 1.91872e10, No Iterations 100 DICPCG: Solving for p, Initial residual = 6.3535e10, Final residual = 5.90779e12, No Iterations 100 DICPCG: Solving for p, Initial residual = 1.16615e09, Final residual = 1.47742e11, No Iterations 100 DICPCG: Solving for p, Initial residual = 5.85987e10, Final residual = 7.23812e12, No Iterations 100 DICPCG: Solving for p, Initial residual = 6.90344e10, Final residual = 1.02724e11, No Iterations 100 DICPCG: Solving for p, Initial residual = 6.25039e10, Final residual = 8.83031e12, No Iterations 100 time step continuity errors : sum local = 953861, global = 176059, cumulative = 188861 ExecutionTime = 335.09 s ClockTime = 336 s Time = 10 smoothSolver: Solving for Ux, Initial residual = 0.314423, Final residual = 0.00535151, No Iterations 4 smoothSolver: Solving for Uy, Initial residual = 0.319104, Final residual = 0.00540234, No Iterations 4 smoothSolver: Solving for Uz, Initial residual = 0.324228, Final residual = 0.00544633, No Iterations 4 DICPCG: Solving for p, Initial residual = 1.42342e08, Final residual = 2.40163e10, No Iterations 100 DICPCG: Solving for p, Initial residual = 3.85046e10, Final residual = 9.44379e12, No Iterations 100 DICPCG: Solving for p, Initial residual = 7.06354e10, Final residual = 2.04228e11, No Iterations 100 DICPCG: Solving for p, Initial residual = 4.03936e10, Final residual = 1.17704e11, No Iterations 100 DICPCG: Solving for p, Initial residual = 4.79926e10, Final residual = 1.53423e11, No Iterations 100 DICPCG: Solving for p, Initial residual = 4.42188e10, Final residual = 1.16098e11, No Iterations 100 time step continuity errors : sum local = 4.0144e+07, global = 5.41432e+06, cumulative = 5.60318e+06 ExecutionTime = 371.52 s ClockTime = 372 s Time = 11 smoothSolver: Solving for Ux, Initial residual = 0.317013, Final residual = 0.00480118, No Iterations 4 smoothSolver: Solving for Uy, Initial residual = 0.320907, Final residual = 0.00484383, No Iterations 4 smoothSolver: Solving for Uz, Initial residual = 0.324356, Final residual = 0.00487736, No Iterations 4 DICPCG: Solving for p, Initial residual = 3.42875e08, Final residual = 4.90087e10, No Iterations 100 DICPCG: Solving for p, Initial residual = 9.60001e10, Final residual = 1.46431e11, No Iterations 100 DICPCG: Solving for p, Initial residual = 2.01496e09, Final residual = 5.44019e11, No Iterations 100 DICPCG: Solving for p, Initial residual = 1.05359e09, Final residual = 1.75239e11, No Iterations 100 DICPCG: Solving for p, Initial residual = 1.24825e09, Final residual = 2.67017e11, No Iterations 100 DICPCG: Solving for p, Initial residual = 1.14462e09, Final residual = 2.11082e11, No Iterations 100 time step continuity errors : sum local = 6.98003e+08, global = 1.12935e+08, cumulative = 1.18538e+08 ExecutionTime = 407.69 s ClockTime = 408 s Last edited by wyldckat; February 22, 2014 at 05:14. Reason: Added [CODE][/CODE] 

February 21, 2014, 03:35 

#6 
New Member
Steven
Join Date: Dec 2013
Location: Perth
Posts: 15
Rep Power: 12 
Hi Himanshu,
Yes, finally I was able to solve the problem. If I'm not mistaken, I either put wrong boundary conditions or several boundary conditions did not suit well to the case I was simulating. However, Anton's suggestion about mass continuity also plays important role to make sure small continuity error. What are you simulating? Maybe just put more detail about your case since it's hard to guess what's going on just from your post. Cheers 

February 21, 2014, 20:49 

#7 
Member
Francis
Join Date: Jan 2014
Location: Toronto
Posts: 50
Rep Power: 12 
Recalculate your time step with ref to Courant number, speed and the cell size.
You should be fine. Chose the mesh size based on max U and minimum cell size or vise versa. Last edited by wyldckat; February 22, 2014 at 05:13. Reason: merged two posts, since they were posted 1 minute apart 

February 21, 2014, 21:09 

#8 
New Member
Steven
Join Date: Dec 2013
Location: Perth
Posts: 15
Rep Power: 12 
And also make sure your meshes are ok. Try to type "checkmesh" in command window and see the results. Last time I got problem because some of my meshes were not properly generated.


December 12, 2014, 09:16 

#9 
Member
Francis
Join Date: Jan 2014
Location: Toronto
Posts: 50
Rep Power: 12 
Also ensure uniform mesh near the movingWall. Grading of 1 in an interior domain created near the wall. Others can be graded.


December 14, 2014, 18:17 

#10 
Member
Francis
Join Date: Jan 2014
Location: Toronto
Posts: 50
Rep Power: 12 
Use laplaceFaceDecomposition in you dynamicMesh solver with motionU.


Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
AMI speed performance  danny123  OpenFOAM  21  October 24, 2020 04:13 
Transient simulation not converging  skabilan  OpenFOAM Running, Solving & CFD  14  December 16, 2019 23:12 
Mesquite  Adaptive mesh refinement / coarsening?  philippose  OpenFOAM Running, Solving & CFD  94  January 27, 2016 09:40 
mixerVesselAMI2D's mass is not balancing  sharonyue  OpenFOAM Running, Solving & CFD  6  June 10, 2013 09:34 
Compiling problems with hello worldC  fw407  OpenFOAM Installation  21  January 6, 2008 17:38 