[celestial]

Good day everybody

I see we have

const volScalarField& T = thermo->T();

in createFields.H of solver dbnsTurbFoam

(because thermo->T() is const )

which is really too bad as temperature can very well vary with time for a compressible fluid.

Now is that a bug or I'm missing something really basic here ?

Regards to all

[celestial]

Well, I decided to run the dbns tutorial on the naca airfoil and even though the initial T is 300 deg K everywhere on the mesh, it changes spatially with time in spite of it being a const reference. Somehow, the code designers found a way in C++ to change the memory location value. It is still mysterious why it should be a const reference when enthalpy h and pressure p are non constant references.