|
[Sponsors] |
September 8, 2008, 06:25 |
hello friends,
|
#1 |
Senior Member
|
hello friends,
i am often getting this type of problem i dont know how to solve it please help me to rectify this error Time = 6.9 DILUPBiCG: Solving for Ux, Initial residual = 0.447215, Final residual = 0.00723408, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.432012, Final residual = 0.00989161, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 0.434152, Final residual = 0.00850153, No Iterations 1 DICPCG: Solving for p, Initial residual = 0.884161, Final residual = 157.723, No Iterations 1001 DICPCG: Solving for p, Initial residual = 0.00260549, Final residual = 8.82254e-07, No Iterations 153 DICPCG: Solving for p, Initial residual = 0.00121632, Final residual = 9.73498e-07, No Iterations 95 time step continuity errors : sum local = 5.82606e+86, global = 1.27825e+70, cumulative = 1.27825e+70 #0 Foam::error::printStack(Foam:stream&) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM. so" #1 Foam::sigFpe::sigFpeHandler(int) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM. so" #2 ?? in "/lib64/libc.so.6" #3 Foam::DILUPreconditioner::calcReciprocalD(Foam::Fi eld<double>&, Foam::lduMatrix const&) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM. so" #4 Foam::DILUPreconditioner::DILUPreconditioner(Foam: :lduMatrix::solver const&, Foam::Istream&) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM. so" #5 Foam::lduMatrix::preconditioner::addasymMatrixCons tructorToTable<foam::dilupreco nditioner>::New(Foam::lduMatrix::solver const&, Foam::Istream&) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM. so" #6 Foam::lduMatrix::preconditioner::New(Foam::lduMatr ix::solver const&, Foam::Istream&) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM. so" #7 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM. so" #8 Foam::fvMatrix<double>::solve(Foam::Istream&) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVol ume.so" #9 Foam::lduMatrix::solverPerformance Foam::solve<double>(Foam::tmp<foam::fvmatrix<doubl e> > const&) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libincompres sibleRASModels.so" #10 Foam::incompressible::RASModels::kEpsilon::correct () in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libincompres sibleRASModels.so" #11 main in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/applications/bin/linux64GccDPOpt /simpleFoam" #12 __libc_start_main in "/lib64/libc.so.6" #13 Foam::regIOobject::readIfModified() in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/applications/bin/linux64GccDPOpt /simpleFoam" Fließkommafehler junkers-selv_si:~/OpenFOAM1.5/selv_si-1.5/run/tutorials/simpleFoam/kabine_withp/ kabine_p_KEpsilon 5:05 > #0 Foam::error::printStack(Foam:stream&) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM. so" #1 Foam::sigFpe::sigFpeHandler(int) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM. so" #2 ?? in "/lib64/libc.so.6" #3 Foam::DILUPreconditioner::calcReciprocalD(Foam::Fi eld<double>&, Foam::lduMatrix const&) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM. so" #4 Foam::DILUPreconditioner::DILUPreconditioner(Foam: :lduMatrix::solver const&, Foam::Istream&) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM. so" #5 Foam::lduMatrix::preconditioner::addasymMatrixCons tructorToTable<foam::dilupreco nditioner>::New(Foam::lduMatrix::solver const&, Foam::Istream&) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM. so" #6 Foam::lduMatrix::preconditioner::New(Foam::lduMatr ix::solver const&, Foam::Istream&) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM. so" #7 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM. so" #8 Foam::fvMatrix<foam::vector<double> >::solve(Foam::Istream&) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/applications/bin/linux64GccDPOpt /simpleFoam" #9 Foam::lduMatrix::solverPerformance Foam::solve<foam::vector<double> >(Foam::tmp<foam::fvmatrix<foam::vector<double> > > const&) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/applications/bin/linux64GccDPOpt /simpleFoam" #10 main in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/applications/bin/linux64GccDPOpt /simpleFoam" #11 __libc_start_main in "/lib64/libc.so.6" #12 Foam::regIOobject::readIfModified() in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/applications/bin/linux64GccDPOpt /simpleFoam" [1] Fließkommafehler simpleFoam > log its simple foam delta T is 0.1 if u people need any more info please mail me i am awaiting for someone help with friendly, sivakumar |
|
September 8, 2008, 06:58 |
This means that your simulatio
|
#2 |
Member
Michael Wild
Join Date: Mar 2009
Location: Bern, Switzerland
Posts: 79
Rep Power: 17 |
This means that your simulation is diverging. If you look at the output just before OF crashes, you see:
time step continuity errors : sum local = 5.82606e+86, global = 1.27825e+70, cumulative = 1.27825e+70 And having read that you are using deltaT=0.1, honestly it's no wonder things blow up right into your face! Try something like 0.0001, which, depending on your case might still be too large. Unfortunately you don't show enough output to see a CFL number, but make sure that number does not grow too fast (i.e. grow from 0.5 to 10 in three time steps, or some such) and remains moderately low. |
|
September 8, 2008, 07:13 |
Hi, sivakumar,
First, I am
|
#3 |
Member
Ivan Lau
Join Date: Mar 2009
Location: Hong Kong
Posts: 56
Rep Power: 17 |
Hi, sivakumar,
First, I am 'relatively' new with OpenFOAM, but I have I think we simply need more information to help you. Since it is simpleFoam (steadyState problem), I think time step should not be the problem. In fact, many people including myself use delta T=1 instead. Did you checkMesh? change relaxation factor? numerical schemes? .... regards, iL |
|
September 8, 2008, 07:42 |
hi friends,
thank
|
#4 |
Senior Member
|
hi friends,
thank you for u r replay, i have given here some time step also please have look on it and then give me suggestion Time = 5.9 DILUPBiCG: Solving for Ux, Initial residual = 0.203474, Final residual = 0.00729043, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.208975, Final residual = 0.00742262, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 0.181622, Final residual = 0.00542371, No Iterations 1 DICPCG: Solving for p, Initial residual = 0.641968, Final residual = 9.61305e-07, No Iterations 361 DICPCG: Solving for p, Initial residual = 0.126658, Final residual = 9.84017e-07, No Iterations 307 DICPCG: Solving for p, Initial residual = 0.0280797, Final residual = 9.80951e-07, No Iterations 288 time step continuity errors : sum local = 3.55806e-08, global = -1.23831e-11, cumulative = 2.92999e-09 DILUPBiCG: Solving for epsilon, Initial residual = 0.00730311, Final residual = 6.51562e-05, No Iterations 1 bounding epsilon, min: -5.43626e-05 max: 159.868 average: 2.47117 DILUPBiCG: Solving for k, Initial residual = 0.0421587, Final residual = 0.00162023, No Iterations 1 bounding k, min: -1.72949e-06 max: 4.85201 average: 0.0675263 ExecutionTime = 36934.3 s ClockTime = 83185 s Time = 6 DILUPBiCG: Solving for Ux, Initial residual = 0.175165, Final residual = 0.0104977, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.174033, Final residual = 0.00756013, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 0.15981, Final residual = 0.0060471, No Iterations 1 DICPCG: Solving for p, Initial residual = 0.425168, Final residual = 9.66345e-07, No Iterations 344 DICPCG: Solving for p, Initial residual = 0.120603, Final residual = 9.60781e-07, No Iterations 302 DICPCG: Solving for p, Initial residual = 0.0293922, Final residual = 9.87513e-07, No Iterations 277 time step continuity errors : sum local = 2.74681e-08, global = -4.01837e-13, cumulative = 2.92959e-09 DILUPBiCG: Solving for epsilon, Initial residual = 0.0128378, Final residual = 0.000130153, No Iterations 1 bounding epsilon, min: -0.000754239 max: 607.271 average: 2.43784 DILUPBiCG: Solving for k, Initial residual = 0.0335742, Final residual = 0.00140461, No Iterations 1 bounding k, min: -3.99416e-06 max: 6.34245 average: 0.0667909 ExecutionTime = 37492 s ClockTime = 84583 s Time = 6.1 DILUPBiCG: Solving for Ux, Initial residual = 0.155243, Final residual = 0.00782558, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.173157, Final residual = 0.00700699, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 0.147001, Final residual = 0.00536645, No Iterations 1 DICPCG: Solving for p, Initial residual = 0.371839, Final residual = 9.31589e-07, No Iterations 354 DICPCG: Solving for p, Initial residual = 0.130531, Final residual = 9.31035e-07, No Iterations 319 DICPCG: Solving for p, Initial residual = 0.0423754, Final residual = 9.53662e-07, No Iterations 312 time step continuity errors : sum local = 1.93093e-08, global = 1.00267e-11, cumulative = 2.93962e-09 DILUPBiCG: Solving for epsilon, Initial residual = 0.00576475, Final residual = 0.000135716, No Iterations 1 bounding epsilon, min: -0.00264849 max: 619.858 average: 2.40535 DILUPBiCG: Solving for k, Initial residual = 0.0386964, Final residual = 0.00180239, No Iterations 1 bounding k, min: -2.20046e-06 max: 6.51581 average: 0.0660643 ExecutionTime = 38090.8 s ClockTime = 85806 s Time = 6.2 DILUPBiCG: Solving for Ux, Initial residual = 0.355442, Final residual = 0.017519, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.187081, Final residual = 0.00915445, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 0.579492, Final residual = 0.0206377, No Iterations 1 DICPCG: Solving for p, Initial residual = 0.598516, Final residual = 9.38731e-07, No Iterations 358 DICPCG: Solving for p, Initial residual = 0.999974, Final residual = 9.99365e-07, No Iterations 359 DICPCG: Solving for p, Initial residual = 0.249458, Final residual = 9.75732e-07, No Iterations 300 time step continuity errors : sum local = 0.0108284, global = 8.65464e-06, cumulative = 8.65758e-06 DILUPBiCG: Solving for epsilon, Initial residual = 0.129531, Final residual = 0.0111364, No Iterations 1 DILUPBiCG: Solving for k, Initial residual = 0.949008, Final residual = 0.0779747, No Iterations 1 ExecutionTime = 38700.6 s ClockTime = 87270 s Time = 6.3 DILUPBiCG: Solving for Ux, Initial residual = 0.403861, Final residual = 0.0258618, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.487313, Final residual = 0.0348382, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 0.421526, Final residual = 0.0214651, No Iterations 1 DICPCG: Solving for p, Initial residual = 0.948659, Final residual = 9.7634e-07, No Iterations 988 DICPCG: Solving for p, Initial residual = 0.999085, Final residual = 9.87094e-07, No Iterations 334 DICPCG: Solving for p, Initial residual = 0.231313, Final residual = 9.5807e-07, No Iterations 322 time step continuity errors : sum local = 1.75909, global = -0.370928, cumulative = -0.370919 DILUPBiCG: Solving for epsilon, Initial residual = 1, Final residual = 0.0051317, No Iterations 2 DILUPBiCG: Solving for k, Initial residual = 0.999994, Final residual = 0.00662627, No Iterations 1 ExecutionTime = 39651.6 s ClockTime = 89377 s Time = 6.4 DILUPBiCG: Solving for Ux, Initial residual = 0.980097, Final residual = 0.035499, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.973669, Final residual = 0.0330624, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 0.906473, Final residual = 0.0220345, No Iterations 1 DICPCG: Solving for p, Initial residual = 0.995252, Final residual = 9.79861e-07, No Iterations 837 DICPCG: Solving for p, Initial residual = 0.095584, Final residual = 9.85317e-07, No Iterations 562 DICPCG: Solving for p, Initial residual = 0.20073, Final residual = 9.70705e-07, No Iterations 644 time step continuity errors : sum local = 8.4294, global = -0.00271281, cumulative = -0.373632 DILUPBiCG: Solving for epsilon, Initial residual = 0.999265, Final residual = 0.0489265, No Iterations 1 DILUPBiCG: Solving for k, Initial residual = 0.998745, Final residual = 0.0434789, No Iterations 1 ExecutionTime = 40829.7 s ClockTime = 92024 s Time = 6.5 DILUPBiCG: Solving for Ux, Initial residual = 0.415975, Final residual = 0.0094284, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.375868, Final residual = 0.00964882, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 0.48942, Final residual = 0.00950915, No Iterations 1 DICPCG: Solving for p, Initial residual = 0.906716, Final residual = 7.84503e-05, No Iterations 1001 DICPCG: Solving for p, Initial residual = 0.321106, Final residual = 9.93537e-07, No Iterations 941 DICPCG: Solving for p, Initial residual = 0.284919, Final residual = 9.68267e-07, No Iterations 969 time step continuity errors : sum local = 395.228, global = -0.0742896, cumulative = -0.447922 DILUPBiCG: Solving for epsilon, Initial residual = 0.611765, Final residual = 0.00406887, No Iterations 2 bounding epsilon, min: -1.95427e+13 max: 5.50778e+20 average: 7.45698e+14 DILUPBiCG: Solving for k, Initial residual = 0.325837, Final residual = 0.0125421, No Iterations 1 ExecutionTime = 42536 s ClockTime = 104002 s Time = 6.6 DILUPBiCG: Solving for Ux, Initial residual = 0.369153, Final residual = 0.00714651, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.747845, Final residual = 0.0389383, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 0.425775, Final residual = 0.00888966, No Iterations 1 DICPCG: Solving for p, Initial residual = 0.985586, Final residual = 0.000145983, No Iterations 1001 DICPCG: Solving for p, Initial residual = 0.923968, Final residual = 0.000179442, No Iterations 1001 DICPCG: Solving for p, Initial residual = 0.979707, Final residual = 4.50676e-06, No Iterations 1001 time step continuity errors : sum local = 5.02621e+07, global = -21516.5, cumulative = -21516.9 DILUPBiCG: Solving for epsilon, Initial residual = 0.90413, Final residual = 0.0296732, No Iterations 1 bounding epsilon, min: -8.63445e+16 max: 4.20959e+29 average: 1.71743e+23 DILUPBiCG: Solving for k, Initial residual = 0.99993, Final residual = 4.65803e-05, No Iterations 1 bounding k, min: -1.20571e+10 max: 2.89022e+21 average: 7.184e+15 ExecutionTime = 44393.3 s ClockTime = 127169 s Time = 6.7 DILUPBiCG: Solving for Ux, Initial residual = 2.20736e-05, Final residual = 1.74502e-08, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.00112565, Final residual = 1.56897e-06, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 1.04156e-05, Final residual = 4.41427e-08, No Iterations 1 DICPCG: Solving for p, Initial residual = 0.967611, Final residual = 0.0060076, No Iterations 1001 DICPCG: Solving for p, Initial residual = 0.506632, Final residual = 0.00309892, No Iterations 1001 DICPCG: Solving for p, Initial residual = 0.0946733, Final residual = 0.00115806, No Iterations 1001 time step continuity errors : sum local = 3.95103e+28, global = 2.38804e+27, cumulative = 2.38804e+27 DILUPBiCG: Solving for epsilon, Initial residual = 1, Final residual = 0.0789002, No Iterations 14 DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 0.0371839, No Iterations 1 bounding k, min: -1.75461e+33 max: 2.09263e+45 average: 5.0326e+39 ExecutionTime = 46199 s ClockTime = 145249 s Time = 6.8 DILUPBiCG: Solving for Ux, Initial residual = 0.549612, Final residual = 0.0120315, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.500578, Final residual = 0.00733108, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 0.412698, Final residual = 0.00761554, No Iterations 1 DICPCG: Solving for p, Initial residual = 1, Final residual = 148.923, No Iterations 1001 DICPCG: Solving for p, Initial residual = 1.89001e-05, Final residual = 6.55384e-07, No Iterations 5 DICPCG: Solving for p, Initial residual = 5.99975e-05, Final residual = 8.21955e-07, No Iterations 5 time step continuity errors : sum local = 2.4247e+58, global = -5.25818e+50, cumulative = -5.25818e+50 DILUPBiCG: Solving for epsilon, Initial residual = 1, Final residual = 0.0351205, No Iterations 35 bounding epsilon, min: 1.70467e-20 max: 1.14436e+140 average: 2.63994e+134 DILUPBiCG: Solving for k, Initial residual = 1.41429e-09, Final residual = 1.41429e-09, No Iterations 0 ExecutionTime = 46858.8 s ClockTime = 147782 s Time = 6.9 DILUPBiCG: Solving for Ux, Initial residual = 0.447215, Final residual = 0.00723408, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.432012, Final residual = 0.00989161, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 0.434152, Final residual = 0.00850153, No Iterations 1 DICPCG: Solving for p, Initial residual = 0.884161, Final residual = 157.723, No Iterations 1001 DICPCG: Solving for p, Initial residual = 0.00260549, Final residual = 8.82254e-07, No Iterations 153 DICPCG: Solving for p, Initial residual = 0.00121632, Final residual = 9.73498e-07, No Iterations 95 time step continuity errors : sum local = 5.82606e+86, global = 1.27825e+70, cumulative = 1.27825e+70 #0 Foam::error::printStack(Foam:stream&) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM. so" #1 Foam::sigFpe::sigFpeHandler(int) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM. so" #2 ?? in "/lib64/libc.so.6" #3 Foam::DILUPreconditioner::calcReciprocalD(Foam::Fi eld<double>&, Foam::lduMatrix const&) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM. so" #4 Foam::DILUPreconditioner::DILUPreconditioner(Foam: :lduMatrix::solver const&, Foam::Istream&) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM. so" #5 Foam::lduMatrix::preconditioner::addasymMatrixCons tructorToTable<foam::dilupreco nditioner>::New(Foam::lduMatrix::solver const&, Foam::Istream&) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM. so" #6 Foam::lduMatrix::preconditioner::New(Foam::lduMatr ix::solver const&, Foam::Istream&) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM. so" #7 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM. so" #8 Foam::fvMatrix<double>::solve(Foam::Istream&) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVol ume.so" #9 Foam::lduMatrix::solverPerformance Foam::solve<double>(Foam::tmp<foam::fvmatrix<doubl e> > const&) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libincompres sibleRASModels.so" #10 Foam::incompressible::RASModels::kEpsilon::correct () in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libincompres sibleRASModels.so" #11 main in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/applications/bin/linux64GccDPOpt /simpleFoam" #12 __libc_start_main in "/lib64/libc.so.6" #13 Foam::regIOobject::readIfModified() in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/applications/bin/linux64GccDPOpt /simpleFoam" Fließkommafehler junkers-selv_si:~/OpenFOAM1.5/selv_si-1.5/run/tutorials/simpleFoam/kabine_withp/ kabine_p_KEpsilon 5:05 > #0 Foam::error::printStack(Foam:stream&) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM. so" #1 Foam::sigFpe::sigFpeHandler(int) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM. so" #2 ?? in "/lib64/libc.so.6" #3 Foam::DILUPreconditioner::calcReciprocalD(Foam::Fi eld<double>&, Foam::lduMatrix const&) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM. so" #4 Foam::DILUPreconditioner::DILUPreconditioner(Foam: :lduMatrix::solver const&, Foam::Istream&) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM. so" #5 Foam::lduMatrix::preconditioner::addasymMatrixCons tructorToTable<foam::dilupreco nditioner>::New(Foam::lduMatrix::solver const&, Foam::Istream&) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM. so" #6 Foam::lduMatrix::preconditioner::New(Foam::lduMatr ix::solver const&, Foam::Istream&) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM. so" #7 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM. so" #8 Foam::fvMatrix<foam::vector<double> >::solve(Foam::Istream&) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/applications/bin/linux64GccDPOpt /simpleFoam" #9 Foam::lduMatrix::solverPerformance Foam::solve<foam::vector<double> >(Foam::tmp<foam::fvmatrix<foam::vector<double> > > const&) in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/applications/bin/linux64GccDPOpt /simpleFoam" #10 main in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/applications/bin/linux64GccDPOpt /simpleFoam" #11 __libc_start_main in "/lib64/libc.so.6" #12 Foam::regIOobject::readIfModified() in "/home/junkers/selv_si/OpenFOAM1.5/OpenFOAM-1.5/applications/bin/linux64GccDPOpt /simpleFoam" [1] Fließkommafehler simpleFoam > log with friendly, sivakumar |
|
September 8, 2008, 08:07 |
hi ivan lau,
|
#5 |
Senior Member
|
hi ivan lau,
i didnt change any numerical scheme, i changed only turbulence model,k epsilon nonorthogonalcorretors 2 nothing else with friendly, sivakumar |
|
September 8, 2008, 16:40 |
hi, sivakumar,
I can tell t
|
#6 |
Member
Ivan Lau
Join Date: Mar 2009
Location: Hong Kong
Posts: 56
Rep Power: 17 |
hi, sivakumar,
I can tell that your number of iterations of p is no good. First, you can lower the relaxation factors under the file system/fvSolution. Especially p. The default number is 0.3. But, try to lower it to 0.2 or lower. It usually helps. Second, you can refine your mesh, and it helps me sometimes and check if the number of iterations of p can be lower. Run 'checkMesh' (post your result here). Thrid, you can change your delta T to 1. Is your mesh highly nonOrthgonal? Also, you may change your pRefValue is your problem open or close? If it is open, you can change your pRefValue 101325. I hope it helps. Regards, iL |
|
September 9, 2008, 04:32 |
hi,
i have checked my mesh
|
#7 |
Senior Member
|
hi,
i have checked my mesh here is the result junkers-selv_si:~/OpenFOAM1.5/selv_si-1.5/run/tutorials/simpleFoam/kabine_withp/ kabine_p_KOmega 9:12 > checkMesh /*---------------------------------------------------------------------------*\ | ========= | | | \ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \ / O peration | Version: 1.5 | | \ / A nd | Web: http://www.OpenFOAM.org | | \/ M anipulation | | \*---------------------------------------------------------------------------*/ Exec : checkMesh Date : Sep 09 2008 Time : 09:12:20 Host : junkers PID : 11585 Case : /home/junkers/selv_si/OpenFOAM1.5/selv_si-1.5/run/tutorials/simpleFoam/kabine_wi thp/kabine_p_KOmega nProcs : 1 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create polyMesh for time = constant Time = constant Mesh stats points: 896946 faces: 6308249 internal faces: 6030423 cells: 2823525 boundary patches: 17 point zones: 0 face zones: 0 cell zones: 0 Number of cells of each type: hexahedra: 0 prisms: 1044572 wedges: 0 pyramids: 0 tet wedges: 0 tetrahedra: 1778953 polyhedra: 0 Checking topology... Boundary definition OK. Point usage OK. Upper triangular ordering OK. Topological cell zip-up check OK. Face vertices OK. Face-face connectivity OK. Number of regions: 1 (OK). Checking patch topology for multiply connected surfaces ... Patch Faces Points Surface INLE1 1711 1616 ok (not multiply connected) INLE2 1496 1504 ok (not multiply connected) INLE3 1834 1670 ok (not multiply connected) INLE4 1731 1739 ok (not multiply connected) OUTL5 1214 1133 ok (not multiply connected) OUTL6 1214 1133 ok (not multiply connected) SYMP7 3615 2681 ok (not multiply connected) SYMP8 3868 2825 ok (not multiply connected) WALL9 1204 753 ok (not multiply connected) WALL10 36525 19097 ok (not multiply connected) WALL11 13890 7623 ok (not multiply connected) WALL12 9335 4768 ok (not multiply connected) WALL13 47426 24351 ok (not multiply connected) WALL14 169 135 ok (not multiply connected) WALL15 167 134 ok (not multiply connected) WALL16 1376 831 ok (not multiply connected) defaultFaces 151051 78077 ok (not multiply connected) Checking geometry... Domain bounding box: (-0.115964 -0.0297577 -1.49236) (3.1635 2.12222 -1.0045e-05) Boundary openness (-2.44902e-16 4.3874e-16 -1.36129e-16) OK. Max cell openness = 1.4646e-15 OK. Max aspect ratio = 82.173 OK. Minumum face area = 5.25232e-08. Maximum face area = 0.0451964. Face area magnitudes OK. Min volume = 5.54617e-11. Max volume = 0.00214196. Total volume = 7.40215. Cell volumes OK. Mesh non-orthogonality Max: 83.4504 average: 18.2635 *Number of severely non-orthogonal faces: 230. Non-orthogonality check OK. <<Writing 230 non-orthogonal faces to set nonOrthoFaces Face pyramids OK. ***Max skewness = 7.88912, 259 highly skew faces detected which may impair the quality of the results <<Writing 259 skew faces to set skewFaces All angles in faces OK. Face flatness (1 = flat, 0 = butterfly) : average = 0.998791 min = 0.902594 All face flatness OK. Failed 1 mesh checks. End with friendly, sivakumar |
|
September 9, 2008, 06:56 |
Hi, Sivakumar,
Well, you pr
|
#8 |
Member
Ivan Lau
Join Date: Mar 2009
Location: Hong Kong
Posts: 56
Rep Power: 17 |
Hi, Sivakumar,
Well, you probably should fix your mesh. How did you create your mesh? which program did you use? Cheers, iL |
|
September 9, 2008, 08:06 |
hi,
thank you, i am not cr
|
#9 |
Senior Member
|
hi,
thank you, i am not creating mesh, they are using star-CD to create mesh whats the problem there. can u explain the problem? with friendly, sivakumar |
|
September 9, 2008, 13:53 |
Hi, Sivakumar,
In your chec
|
#10 |
Member
Ivan Lau
Join Date: Mar 2009
Location: Hong Kong
Posts: 56
Rep Power: 17 |
Hi, Sivakumar,
In your checkMesh: *Number of severely non-orthogonal faces: 230. Non-orthogonality check OK. <<Writing 230 non-orthogonal faces to set nonOrthoFaces ***Max skewness = 7.88912, 259 highly skew faces detected which may impair the quality of the results <<Writing 259 skew faces to set skewFaces I think you probably have some very sharp non-orthogonal cells (non-parallel to axis). I know you can use paraFoam or some other program to view these problem cells and try to fix them. I personally didn't use this method, but you can find more details by searchng this forum. Or, you can simply find them by looking your mesh in paraFoam if they are obvious. I don't use Star-CD, so I can't recommend on it. Regards, iL |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Error : access : Unknown variable & Phase-domain | MKP | FLUENT | 0 | August 28, 2006 04:40 |
UNKNOWN ERROR | Robbie | FLUENT | 0 | July 15, 2005 20:27 |
Unknown Error | Robbie | FLUENT | 0 | July 12, 2005 18:05 |
another unknown Pre error | Louwrens | CFX | 3 | November 20, 2003 17:34 |
Unknown error in running star CD | Rahul | Siemens | 0 | July 5, 2003 11:03 |