[Sugajen]

Hi again,

I am attaching my test case and the entire solver along with this post. There is something very crucial that I am going wrong that is causing the scalar solution to diverge. I am getting NaNs in every time step of the problem. But the pressure and velocity fields are looking good. Please do have a look and comment!

Thank you!
The cn term diverging in the code execution

Time = 0.074

smoothSolver:  Solving for Ux, Initial residual = 0.00219734, Final residual = 1.12576e-05, No Iterations 6
smoothSolver:  Solving for Uy, Initial residual = 0.01411, Final residual = 9.28678e-05, No Iterations 6
smoothSolver:  Solving for Uz, Initial residual = 0.0065224, Final residual = 2.29002e-05, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.00881422, Final residual = 4.36809e-05, No Iterations 6
time step continuity errors : sum local = 3.02676e-06, global = 9.44383e-07, cumulative = 9.34871e-06
smoothSolver:  Solving for cn, Initial residual = 0.286678, Final residual = 0.0192822, No Iterations 1000
ExecutionTime = 43.38 s  ClockTime = 44 s

Time = 0.075

smoothSolver:  Solving for Ux, Initial residual = 0.00212959, Final residual = 1.06745e-05, No Iterations 6
smoothSolver:  Solving for Uy, Initial residual = 0.0137242, Final residual = 8.92657e-05, No Iterations 6
smoothSolver:  Solving for Uz, Initial residual = 0.0061777, Final residual = 2.19587e-05, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.00850799, Final residual = 4.17185e-05, No Iterations 6
time step continuity errors : sum local = 2.88605e-06, global = 9.72972e-07, cumulative = 1.03217e-05
smoothSolver:  Solving for cn, Initial residual = 0.273087, Final residual = 0.0198734, No Iterations 1000
ExecutionTime = 45.44 s  ClockTime = 46 s

Time = 0.076

smoothSolver:  Solving for Ux, Initial residual = 0.0020727, Final residual = 1.01133e-05, No Iterations 6
smoothSolver:  Solving for Uy, Initial residual = 0.0133341, Final residual = 8.56583e-05, No Iterations 6
smoothSolver:  Solving for Uz, Initial residual = 0.00585649, Final residual = 2.0978e-05, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.00821889, Final residual = 4.02368e-05, No Iterations 6
time step continuity errors : sum local = 2.77443e-06, global = 9.80575e-07, cumulative = 1.13023e-05
smoothSolver:  Solving for cn, Initial residual = 0.256031, Final residual = 0.0226206, No Iterations 1000
ExecutionTime = 47.58 s  ClockTime = 49 s

Time = 0.077

smoothSolver:  Solving for Ux, Initial residual = 0.0020286, Final residual = 9.555e-06, No Iterations 6
smoothSolver:  Solving for Uy, Initial residual = 0.0129335, Final residual = 8.19352e-05, No Iterations 6
smoothSolver:  Solving for Uz, Initial residual = 0.00555411, Final residual = 1.99756e-05, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.00795094, Final residual = 3.91401e-05, No Iterations 6
time step continuity errors : sum local = 2.68237e-06, global = 9.69506e-07, cumulative = 1.22718e-05
smoothSolver:  Solving for cn, Initial residual = 0.241686, Final residual = 0.0271197, No Iterations 1000
ExecutionTime = 49.64 s  ClockTime = 51 s

Time = 0.078

smoothSolver:  Solving for Ux, Initial residual = 0.00199783, Final residual = 8.98516e-06, No Iterations 6
smoothSolver:  Solving for Uy, Initial residual = 0.0125206, Final residual = 7.81573e-05, No Iterations 6
smoothSolver:  Solving for Uz, Initial residual = 0.00526388, Final residual = 1.89463e-05, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.00769432, Final residual = 3.80934e-05, No Iterations 6
time step continuity errors : sum local = 2.5859e-06, global = 9.51297e-07, cumulative = 1.32231e-05
smoothSolver:  Solving for cn, Initial residual = 0.237521, Final residual = 0.0297912, No Iterations 1000
ExecutionTime = 51.7 s  ClockTime = 53 s

Time = 0.079

smoothSolver:  Solving for Ux, Initial residual = 0.00197628, Final residual = 8.42215e-06, No Iterations 6
smoothSolver:  Solving for Uy, Initial residual = 0.0121026, Final residual = 7.44988e-05, No Iterations 6
smoothSolver:  Solving for Uz, Initial residual = 0.00498768, Final residual = 1.79099e-05, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.00743783, Final residual = 3.67183e-05, No Iterations 6
time step continuity errors : sum local = 2.46585e-06, global = 9.34576e-07, cumulative = 1.41576e-05
smoothSolver:  Solving for cn, Initial residual = 0.235452, Final residual = 0.0325307, No Iterations 1000
ExecutionTime = 53.71 s  ClockTime = 55 s

Time = 0.08

smoothSolver:  Solving for Ux, Initial residual = 0.0019587, Final residual = 7.88607e-06, No Iterations 6
smoothSolver:  Solving for Uy, Initial residual = 0.011688, Final residual = 7.10823e-05, No Iterations 6
smoothSolver:  Solving for Uz, Initial residual = 0.00472905, Final residual = 1.68969e-05, No Iterations 6
GAMG:  Solving for p, Initial residual = 0.00717117, Final residual = 3.5029e-05, No Iterations 6
time step continuity errors : sum local = 2.32748e-06, global = 9.20884e-07, cumulative = 1.50785e-05
smoothSolver:  Solving for cn, Initial residual = 0.233661, Final residual = 0.0357598, No Iterations 1000
ExecutionTime = 55.9 s  ClockTime = 57 s