
[Sponsors] 
March 6, 2009, 04:05 
Hi Vishal,
can you give me

#41 
Senior Member
Tian
Join Date: Mar 2009
Location: Berlin, germany
Posts: 111
Rep Power: 10 
Hi Vishal,
can you give me a hint where did you find in the OpenFOAM code T.relax()? Thanks Bye Thomas
__________________
BIM HVACTool, The Green Building Simulation Tool for OpenFOAM, Energy Plus and Radiance. 

March 6, 2009, 08:28 
Hello Thomas
it seems that

#42 
Member
Leonardo Honfi Camilo
Join Date: Mar 2009
Location: Delft, Zuid Holland, The Netherlands
Posts: 48
Rep Power: 10 
Hello Thomas
it seems that I got here a little late, the link that you have posted has reached its limit, could you perhaps mail me the file? thanks leo 

March 6, 2009, 08:40 
Hi,
try this one:
http:/

#43 
Senior Member
Tian
Join Date: Mar 2009
Location: Berlin, germany
Posts: 111
Rep Power: 10 
__________________
BIM HVACTool, The Green Building Simulation Tool for OpenFOAM, Energy Plus and Radiance. 

March 6, 2009, 08:47 
thanks

#44 
Member
Leonardo Honfi Camilo
Join Date: Mar 2009
Location: Delft, Zuid Holland, The Netherlands
Posts: 48
Rep Power: 10 
thanks


March 6, 2009, 12:40 
Hi Tian,
I found such a rel

#45 
Senior Member
Vishal Nandigana
Join Date: Mar 2009
Location: Champaign, Illinois, U.S.A
Posts: 208
Rep Power: 11 
Hi Tian,
I found such a relaxation technique mentioned by Prof Hrv.. in this forum I would like to know precisely what each one exactly does in openFOAM... and which is more advantageous... http://www.cfdonline.com/cgibin/Op...ow.cgi?1/10111 Regards Vishal 

July 30, 2009, 06:28 

#46 
Senior Member
J. Cai
Join Date: Apr 2009
Posts: 180
Rep Power: 10 
Hi, dear Foamers, I also meet the same problem when use simpleFoam to do a simple case. I have try some suggestions in this thread, but fails. any other comments?
Thank you very much. have a nice day. Chiven Time = 8 DILUPBiCG: Solving for Ux, Initial residual = 0.0154965, Final residual = 3.56285e05, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.00436824, Final residual = 3.31732e05, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 0.00323331, Final residual = 3.00037e05, No Iterations 1 GAMG: Solving for p, Initial residual = 0.961782, Final residual = 0.00929941, No Iterations 3 time step continuity errors : sum local = 1.72947e+25, global = 2.12055e+19, cumulative = 2.12055e+19 DILUPBiCG: Solving for epsilon, Initial residual = 0.999551, Final residual = 0.0140238, No Iterations 1 bounding epsilon, min: 1.97709e+32 max: 1.13481e+50 average: 7.81416e+44 DILUPBiCG: Solving for k, Initial residual = 0.999999, Final residual = 0.0180508, No Iterations 1 bounding k, min: 1.56134e+22 max: 6.46051e+40 average: 1.40592e+36 ExecutionTime = 72.457 s ClockTime = 73 s Time = 9 DILUPBiCG: Solving for Ux, Initial residual = 0.258312, Final residual = 0.00394438, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.322752, Final residual = 0.00548298, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 0.134508, Final residual = 0.00229101, No Iterations 1 GAMG: Solving for p, Initial residual = 0.999988, Final residual = 0.00957924, No Iterations 10 time step continuity errors : sum local = 6.46819e+24, global = 6.00802e+18, cumulative = 2.72135e+19 DILUPBiCG: Solving for epsilon, Initial residual = 0.0935502, Final residual = 1.68919e17, No Iterations 1 bounding epsilon, min: 5.66454e+59 max: 1.62366e+65 average: 4.84566e+60 DILUPBiCG: Solving for k, Initial residual = 0.271442, Final residual = 0.00498093, No Iterations 1 bounding k, min: 4.51792e+38 max: 2.46246e+56 average: 5.53053e+51 ExecutionTime = 83.2949 s ClockTime = 84 s Time = 10 DILUPBiCG: Solving for Ux, Initial residual = 0.213457, Final residual = 0.000426911, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.194438, Final residual = 0.000360278, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 0.140209, Final residual = 0.000161624, No Iterations 1 GAMG: Solving for p, Initial residual = 0.00575245, Final residual = 5.53957e05, No Iterations 1 time step continuity errors : sum local = 6.46819e+49, global = 8.29418e+43, cumulative = 8.29418e+43 DILUPBiCG: Solving for epsilon, Initial residual = 2.16312e18, Final residual = 2.16312e18, No Iterations 0 bounding epsilon, min: 1.59155e21 max: 2.55066e+77 average: 2.55902e+72 DILUPBiCG: Solving for k, Initial residual = 0.00017126, Final residual = 6.14596e07, No Iterations 1 bounding k, min: 9.10313e+53 max: 3.20352e+90 average: 3.67758e+85 ExecutionTime = 90.6045 s ClockTime = 91 s Time = 11 DILUPBiCG: Solving for Ux, Initial residual = 0.0174846, Final residual = 0.000301742, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.00157757, Final residual = 2.24617e05, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 0.0368666, Final residual = 0.000505577, No Iterations 1 #0 _ZN4Foam5error10printStackERNS_7OstreamE0xb383f0 in "/home/g2/e090012/OpenFOAM/OpenFOAM1.5/lib/linuxIA64GccDPOpt/libOpenFOAM.so" #1 _ZN4Foam6sigFpe13sigFpeHandlerEi0xb01820 in "/home/g2/e090012/OpenFOAM/OpenFOAM1.5/lib/linuxIA64GccDPOpt/libOpenFOAM.so" #2 Uninterpreted: [0xa0000000000107e0] #3 _ZN4Foam17DICPreconditioner15calcReciprocalDERNS_5 FieldIdEERKNS_9lduMatrixE0xef5f6e in "/home/g2/e090012/OpenFOAM/OpenFOAM1.5/lib/linuxIA64GccDPOpt/libOpenFOAM.so" #4 _ZN4Foam11DICSmootherC1ERKNS_4wordERKNS_9lduMatrix ERKNS_10FieldFieldINS_5FieldEdEESB_RKNS_8UPtrListI KNS_17lduInterfaceFieldEEE0xefe2d0 in "/home/g2/e090012/OpenFOAM/OpenFOAM1.5/lib/linuxIA64GccDPOpt/libOpenFOAM.so" #5 _ZN4Foam9lduMatrix8smoother30addsymMatrixConstruct orToTableINS_11DICSmootherEE3NewERKNS_4wordERKS0_R KNS_10FieldFieldINS_5FieldEdEESE_RKNS_8UPtrListIKN S_17lduInterfaceFieldEEE0xeb2c00 in "/home/g2/e090012/OpenFOAM/OpenFOAM1.5/lib/linuxIA64GccDPOpt/libOpenFOAM.so" #6 _ZN4Foam9lduMatrix8smoother3NewERKNS_4wordERKS0_RK NS_10FieldFieldINS_5FieldEdEESB_RKNS_8UPtrListIKNS _17lduInterfaceFieldEEERNS_7IstreamE0xf1a020 in "/home/g2/e090012/OpenFOAM/OpenFOAM1.5/lib/linuxIA64GccDPOpt/libOpenFOAM.so" #7 _ZNK4Foam10GAMGSolver10initVcycleERNS_7PtrListINS_ 5FieldIdEEEES5_RNS1_INS_9lduMatrix8smootherEEE0xedc9e0 in "/home/g2/e090012/OpenFOAM/OpenFOAM1.5/lib/linuxIA64GccDPOpt/libOpenFOAM.so" #8 _ZNK4Foam10GAMGSolver5solveERNS_5FieldIdEERKS2_h0xed4f70 in "/home/g2/e090012/OpenFOAM/OpenFOAM1.5/lib/linuxIA64GccDPOpt/libOpenFOAM.so" #9 _ZN4Foam8fvMatrixIdE5solveERNS_7IstreamE0x1f46d20 in "/home/g2/e090012/OpenFOAM/OpenFOAM1.5/lib/linuxIA64GccDPOpt/libfiniteVolume.so" #10 main in "/home/g2/e090012/OpenFOAM/OpenFOAM1.5/applications/bin/linuxIA64GccDPOpt/simpleFoam" #11 __libc_start_main0x739090 in "/lib/tls/libc.so.6.1" #12 _start in "/home/g2/e090012/OpenFOAM/OpenFOAM1.5/applications/bin/linuxIA64GccDPOpt/simpleFoam" 

August 22, 2009, 10:27 

#47 
New Member
Krzysztof Przysowa
Join Date: Mar 2009
Location: Frimley, Surrey, United Kingdom
Posts: 13
Rep Power: 10 
Hi chiven,
there is obviously something causing time step continuity errors being sky high. I recall that I had a similar problem when, by mistake, in boundary definition velocity direction was from outlet to inlet. Of course it should be other way around. So, I would recommend to check your boundary conditions. Cheers 

August 23, 2009, 05:13 

#48 
Senior Member
J. Cai
Join Date: Apr 2009
Posts: 180
Rep Power: 10 
Hello, Krzysztof, thanks a lot for your reply. You are right, and now the problems are over.
Best regards, Chiven 

November 27, 2009, 08:16 

#49 
New Member
Michael Schwarz
Join Date: Nov 2009
Posts: 5
Rep Power: 9 
Hey Guys!!
I have the same problemes changing my solvers to GAMG like mayol had. FOAM FATAL IO ERROR : expected 'true' or 'false', found yes file: /opt/OpenFOAM/caelinux1.4.1/run/BLENDE/system/fvSolution:::scaleCorrection at line 41. From function operator>>(Istream&, bool&) in file primitives/bool/boolIO.C at line 72. FOAM exiting I just started working with OpenFoam and still be playing arround a bit. Here is my checkMesh log: Create polyMesh for time = constant Time = constant Mesh stats points: 27304 edges: 172646 faces: 283432 internal faces: 268924 cells: 138089 boundary patches: 3 point zones: 0 face zones: 0 cell zones: 0 Number of cells of each type: hexahedra: 0 prisms: 0 wedges: 0 pyramids: 0 tet wedges: 0 tetrahedra: 138089 polyhedra: 0 Checking topology... Boundary definition OK. Point usage OK. Upper triangular ordering OK. Topological cell zipup check OK. Face vertices OK. Faceface connectivity OK. Number of regions: 1 (OK). Checking patch topology for multiply connected surfaces ... Patch Faces Points Surface In 122 77 ok (not multiply connected) Out 172 102 ok (not multiply connected) Walls 14214 7137 ok (not multiply connected) Checking geometry... Domain bounding box: (0.6 0 0.6) (0.6 5 0.6) Boundary openness (4.08264e17 2.77916e17 3.03938e17) OK. Max cell openness = 1.48446e16 OK. Max aspect ratio = 4.30296 OK. Minumum face area = 1.32319e05. Maximum face area = 0.028895. Face area magnitudes OK. Min volume = 2.50715e08. Max volume = 0.00159632. Total volume = 5.60409. Cell volumes OK. Mesh nonorthogonality Max: 51.4236 average: 15.1122 Nonorthogonality check OK. Face pyramids OK. Max skewness = 0.791886 OK. Min/max edge length = 0.00478093 0.295993 OK. All angles in faces OK. All face flatness OK. Mesh OK. To me, it seems that the meshscale is fine. Any idea what causes the problems? Thanks for Help! 

November 28, 2009, 05:32 

#50  
Senior Member
Gijsbert Wierink
Join Date: Mar 2009
Posts: 383
Rep Power: 11 
Hi Michael,
Quote:
__________________
Regards, Gijs 

November 29, 2009, 13:28 

#51 
New Member
Michael Schwarz
Join Date: Nov 2009
Posts: 5
Rep Power: 9 
Hi Gijs,
sorry for the late reply. Have been away for the weekend... I tried to Put "true" there but the same errors occure. It is the same error for U,p,k and epsilon.... No idea what to do 

January 26, 2010, 07:46 
time step continuity errors...remain

#52 
New Member
George Pechlivanoglou
Join Date: Mar 2009
Posts: 4
Rep Power: 10 
Hello everyone!
I also get the time step continuity error and I am having a hard time solving the problem. My grid is built with gmsh and engrid (extrusion, boundary conditions, OpenFOAM export). It is unstructured, quasi 2D (i.e. sides are BL type = empty) and represents a multi element airfoil. the grid works perfectly with simpleFoam and komegaSST until 16° AoA. When I try computations at 25° AoA I get the time step continuity errors. I have tried all the proposals from this thread (e.g. reduction of relaxation factors, high initial epsilon, smaller time step) it still crashes. I am pretty sure that the problem has to do with the massive separation due to high AoA...but I have no idea what to do next. Any proposals ?? Thanx a lot for you help! cheers George 

January 28, 2010, 15:20 

#53 
Senior Member
Gijsbert Wierink
Join Date: Mar 2009
Posts: 383
Rep Power: 11 
Hi George,
I'm not an aerodynamicist or anything, but could it be related to the pressure correction method? SimpleFoam uses the SIMPLE algorithm (of course ), which can be quite sensitive to relaxations factors and time step. The PISO algorithm may be less sensitive and contains an extra correction for skewness (do you have a more skewed mesh for higher angle of attack?). Have you tried e.g. pisoFoam (or turbFoam) to test this or do you need to use simpleFoam? Well, don't know if this is of any help, but it just comes to mind.
__________________
Regards, Gijs 

January 28, 2010, 16:53 
SimpleFoam, Time step continuity errors

#54 
New Member
George Pechlivanoglou
Join Date: Mar 2009
Posts: 4
Rep Power: 10 
gwierink thanks for the tips
I am actually building my unstructured grids with Gmsh and then extrude them with EnGrid...meaning that the whole process is more or less standard for any airfoil I use. I therefore think that the grid itself is not really problematic. My checkMesh does not give me any errors...and the aspect ratios and expansion ratios are good. It might be the sensitivity of simpleFoam to relaxation and time step....but maybe I try to find some other possible solutions as well before I change the model. The problem is that I am running the computations via script which goes from one AoA to the other in order to get the airfoil polar at the end. I suspect that is I change solver between AoA then my polar will have sort of..random errors. Anyway, thanx again for the tips!! cheers George 

February 10, 2010, 22:57 

#55  
New Member
Jonathas Assunção de Castro
Join Date: Mar 2009
Posts: 11
Rep Power: 10 
Quote:
Could you send to me the following files, please? initConvergenceCheck.H convergenceCheck.H pEqn.H UEqn.H Tnahks 

April 11, 2010, 06:29 

#56 
New Member
Join Date: Feb 2010
Posts: 21
Rep Power: 9 
Hello All,
I am trying to simulate an impingingjet using the simpleFoam in OF1.6.x. I am trying to follow the threads and tried different things. But I am still getting the same error. I am posting my error and attaching my initial bc. Your support will be greatly appreciated. Time = 759 DILUPBiCG: Solving for Ux, Initial residual = 0.364807, Final residual = 0.0331038, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.401396, Final residual = 0.0185339, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 0.439448, Final residual = 0.0046591, No Iterations 1 GAMG: Solving for p, Initial residual = 0.569614, Final residual = 0.00428777, No Iterations 66 time step continuity errors : sum local = 3.14096e+42, global = 3.09629e+39, cumulative = 6.80133e+40 DILUPBiCG: Solving for epsilon, Initial residual = 1.08298e09, Final residual = 1.08298e09, No Iterations 0 DILUPBiCG: Solving for k, Initial residual = 0.000115294, Final residual = 2.84805e06, No Iterations 2 ExecutionTime = 47327.3 s ClockTime = 48069 s Time = 760 DILUPBiCG: Solving for Ux, Initial residual = 0.506818, Final residual = 0.00973123, No Iterations 2 DILUPBiCG: Solving for Uy, Initial residual = 0.603294, Final residual = 0.0122224, No Iterations 2 DILUPBiCG: Solving for Uz, Initial residual = 0.280386, Final residual = 0.00518198, No Iterations 2 GAMG: Solving for p, Initial residual = 0.0695246, Final residual = 0.000686704, No Iterations 30 time step continuity errors : sum local = 1.41992e+43, global = 4.71133e+41, cumulative = 5.39147e+41 DILUPBiCG: Solving for epsilon, Initial residual = 3.87187e11, Final residual = 3.87187e11, No Iterations 0 DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 0.0894909, No Iterations 1 bounding k, min: 7.27376e+51 max: 5.05817e+89 average: 5.17561e+84 ExecutionTime = 47386.5 s ClockTime = 48159 s Time = 761 DILUPBiCG: Solving for Ux, Initial residual = 0.417362, Final residual = 0.0145424, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.51114, Final residual = 0.0224704, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 0.376437, Final residual = 0.0288099, No Iterations 1 GAMG: Solving for p, Initial residual = 0.629327, Final residual = 0.00528483, No Iterations 8 time step continuity errors : sum local = 1.11456e+78, global = 4.87405e+75, cumulative = 4.87405e+75 DILUPBiCG: Solving for epsilon, Initial residual = 4.23728e09, Final residual = 4.23728e09, No Iterations 0 #0 Foam::error:rintStack(Foam::Ostream&) in "/home/jish/OpenFOAM/OpenFOAM1.6.x/lib/linux64GccDPOpt/libOpenFOAM.so" #1 Foam::sigFpe::sigFpeHandler(int) in "/home/jish/OpenFOAM/OpenFOAM1.6.x/lib/linux64GccDPOpt/libOpenFOAM.so" #2 ?? in "/lib/libc.so.6" #3 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/jish/OpenFOAM/OpenFOAM1.6.x/lib/linux64GccDPOpt/libOpenFOAM.so" #4 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/home/jish/OpenFOAM/OpenFOAM1.6.x/lib/linux64GccDPOpt/libfiniteVolume.so" #5 Foam::lduMatrix::solverPerformance Foam::solve<double>(Foam::tmp<Foam::fvMatrix<doubl e> > const&) in "/home/jish/OpenFOAM/OpenFOAM1.6.x/lib/linux64GccDPOpt/libincompressibleRASModels.so" #6 Foam::incompressible::RASModels::kEpsilon::correct () in "/home/jish/OpenFOAM/OpenFOAM1.6.x/lib/linux64GccDPOpt/libincompressibleRASModels.so" #7 main in "/home/jish/OpenFOAM/OpenFOAM1.6.x/applications/bin/linux64GccDPOpt/simpleFoam" #8 __libc_start_main in "/lib/libc.so.6" #9 _start at /build/buildd/eglibc2.10.1/csu/../sysdeps/x86_64/elf/start.S:116 Floating point exception thanks regards jish 

February 21, 2016, 06:44 

#57  
New Member
swaps
Join Date: Jul 2015
Posts: 1
Rep Power: 0 
Quote:


October 11, 2017, 11:06 

#58  
New Member
Andrea
Join Date: Sep 2017
Posts: 8
Rep Power: 2 
Quote:
I found that this solution helped keep the run stable. Thanks! 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Time step continuity  anja  OpenFOAM Running, Solving & CFD  29  April 10, 2017 02:18 
Time step size and max iterations per time step  pUl  FLUENT  31  August 21, 2015 04:46 
SELECTING TIME STEP SIZE, NUMBER OF TIME STEP  NITUL KALITA  FLUENT  2  November 22, 2012 09:28 
Speedup with GAMG for simplefoam forward Step  tutlhino  OpenFOAM Running, Solving & CFD  9  June 24, 2007 21:44 
Relation of computational time step with real time  Salman  Main CFD Forum  2  August 3, 2005 14:13 