[sivakumar]

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

[mwild]

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.

[ivanwhlau]

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

[sivakumar]

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

[sivakumar]

hi ivan lau,

i didnt change any numerical scheme,
i changed only turbulence model,k epsilon
nonorthogonalcorretors 2 nothing else


with friendly,
sivakumar

[ivanwhlau]

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

[sivakumar]

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

[ivanwhlau]

Hi, Sivakumar,

Well, you probably should fix your mesh. How did you create your mesh? which program did you use?

Cheers,
iL

[sivakumar]

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

[ivanwhlau]

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