[Canesin]

---------_EDIT_---------- SOLVED: was relaxing in outerloop

Hi all,

In an custom solver I need to add a value in the temperature..

I'm doing the equivalent of:

for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
    {
        tmp<fvScalarMatrix> TEqn
        (
            fvm::ddt(rho*cp, T)
          - fvm::laplacian(K, T)
        );
        
        TEqn().solve();
    }
    Info<< "TEMP: Min/max T: " << min(T).value() << " / " << max(T).value() << endl;
}

T += 6;

But when I visualize the results in paraview it looks like the temperature change has not "propagated", say. If I do T.write() after the temperature change I can see the difference, but in the next time step the solutions looks like the change has never happened.