Hi Fred, I don't think the the residual projection (p94,p95) is your problem! Please post your USERIC, USERBC, complete logfile. Stefan On Jan 19, 2010, at 8:55 AM, [email protected] wrote:

Hello,

when I run my simulation (periodic open-channel flow) I got the following error message at timestep 6: Step 6, t= 1.0043885E-01, DT= 2.1145021E-02, C= 0.203 1.7406E+01 3.5078E+00 Solving for fluid 6 6.1565E-10 6.1565E-10 1.0000E+00 0 alph1x 6 halpha 0 6 Hmholtz VELX: 1 4.3435E-10 6.1565E-10 1.0000E-08 6 2.6858E-09 NaN NaN 1 alph1x 6 halpha 1 NaN 6 5000 **ERROR**: Failed in HMHOLTZ: VELY NaN NaN 1.0000E-08 6 6.2743E+00 8.5949E-01 7.3000E+00 1 alph1x 6 halpha 1 2.0336E+00 6 Hmholtz VELZ: 76 8.4969E-09 8.5949E-01 1.0000E-08 reorthogo: 1 6 1 alpha: NaN 6 1 NaN NaN NaN alph12 6 U-Pres gmres: 120 NaN 1.0000E-05 NaN 1.5230E+01 6 DNORM, DIVEX NaN NaN 6 1.0044E-01 1.2954E+02 Fluid done

In the rea-file the parameter 94 and 95 are set to: P94 = 5.0 P95 = 5.0

If I set P94=0.0 the simulation is running without stopping at timestep 6. I tried to decrease P21 (DIVERGENCE) and P22 (HELMHOLTZ) similar to the turbChannel example, but without succes.

May you help me with this problem?

Thanks a lot.

Best regards, Fred

_______________________________________________ Nek5000-users mailing list [email protected] https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users