Pedro Gomes

In places like `if (tkeNeeded) Energy += GetTke_Inf();` we would instead read tke from the turbulence solver.

@rois1995 since you are looking into DDES stuff can you take a look at this issue please?

One more https://www.cfd-online.com/Forums/su2/246710-laminar-cylinder-tutorial-results.html

Axisymmetry on a node-based code requires special treatment near the axis because some terms are 1/r and so if r is close to zero... If your geometry had the axis...

You could try turning off or simplifying the contributions to the Jacobian for that type of boundary, see CNSSolver.ccp around line 550. For example, leaving only the `Jacobian_i[nDim+1][nDim+1]` or none.

If you restart the heat transfer case (with h=0) from the 0 heat flux case, are the residuals the same or is there a jump?

There is a contribution to the Jacobian because the heat flux depends on the temperature at the node. Does it work better if you remove that contribution?

Thank you for investigating the issue, please open a pull request with the fix. You have not provided a case to reproduce the issue. So it's difficult for someone else...

This is a nice-to-have but there are 2600 uses of cout

A PR was just opened that seems to contain more progress on this topic.