[miladrakhsha]

I am experiencing with the implementation of IBM in OpenFOAM. I am working with simpleIbFoam. It would be great if this method could be implemented in the future version of OpenFOAM ext by default as it is a very powerful method. However, I did not find the currrent implementation robust enough yet. My case terminates after a couple of iterations without any warning or any error message so I dont rally know where the problem is!!! My only gauss is that it has something to do with the "k" equation!
this is a small part of my log file:

Time = 6

DILUPBiCG:  Solving for Ux, Initial residual = 7.0369e-05, Final residual = 6.80687e-11, No Iterations 5
DILUPBiCG:  Solving for Uy, Initial residual = 0.000326453, Final residual = 2.10094e-10, No Iterations 5
DILUPBiCG:  Solving for Uz, Initial residual = 7.49995e-06, Final residual = 1.87974e-11, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.0588972, Final residual = 0.00623836, No Iterations 1000
time step continuity errors : sum local = 1.84801e-06, global = -2.04249e-09, cumulative = -4.51088e-09
DILUPBiCG:  Solving for omega, Initial residual = 0.0135086, Final residual = 3.48746e-05, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 0.0159819, Final residual = 5.767e-05, No Iterations 1
bounding k, min: -0.940922 max: 1.77288 average: 0.00776461
ExecutionTime = 198.19 s  ClockTime = 198 s

forces output:
    forces(pressure, viscous)((0.000193756 -2.15261e-05 0.000364652) (4.89027e-10 7.44523e-11 6.68369e-13))
    moment(pressure, viscous)((2.45693e-05 -3.77346e-05 -1.48994e-05) (-1.28583e-12 5.10029e-12 -2.34689e-11))

Field p min = -19.5442 max = 3.89681
Field U magnitude min = 0 max = 6.16393
forceCoeffs output:
    Cd = -0.013689
    Cl = 0.581642
    Cm = 0.059539

Time = 7

DILUPBiCG:  Solving for Ux, Initial residual = 6.97614e-05, Final residual = 9.57514e-11, No Iterations 5
DILUPBiCG:  Solving for Uy, Initial residual = 0.000280219, Final residual = 1.82e-10, No Iterations 5
DILUPBiCG:  Solving for Uz, Initial residual = 6.51344e-06, Final residual = 2.49325e-11, No Iterations 5
GAMG:  Solving for p, Initial residual = 0.0732228, Final residual = 0.00731821, No Iterations 871
time step continuity errors : sum local = 1.70086e-06, global = -2.23902e-09, cumulative = -6.7499e-09
DILUPBiCG:  Solving for omega, Initial residual = 0.00702349, Final residual = 3.47449e-05, No Iterations 1
DILUPBiCG:  Solving for k, Initial residual = 0.00908902, Final residual = 4.9913e-05, No Iterations 1
bounding k, min: -1.34002 max: 3.64086 average: 0.00851153

THAT IS IT, THIS IS THE END OF THE LOG FILE!!!!!!!!


I tried to use bounded schemes for k and omega in fvSchemes but it seems that it does not recognize bounded scheme !

I would really appreciate if somebody can help.
Thank you

----------------
Moderator note: This specific topic was later discussed here in more detail: http://www.cfd-online.com/Forums/ope...m-1-6-ext.html
The remaining thread is on the topic of porting to foam-extend 3.1