Finding calculation of T_ field
I want to implement the wallHeatFlux BC for the T-Field in my solver. For now I can use the BC but the calculation of the T faces of the boundary is not working. I get still the values of the initial Temperature (its like i have a fixedValue BC).
How do I test the values of the T_ faces?
After calculate of hEqn.H there is the thermo.correct() code.
I had a look at the correct function and there I found that:
// force the saving of the old-time values
if (counter == propertyUpdate)
Info << "Updating look-up table extractions..." << endl;
counter = 0;
if (counter_mass_fractions == massFractionsUpdate)
Info << "Updating mass fraction extractions..." << endl;
counter_mass_fractions = 0;
Debug(); // added by myself
And in the Debug-Function I am getting the face-values:
Info << "Debug" << endl;
const fvPatchScalarField & pcsi = csi_.boundaryField()[patchi];
const fvPatchScalarField & pT = T_.boundaryField()[patchi];
Info << "T = " << pT[facei] << endl;
Now my question is, how can I figure out why the BC wallHeatTransfer is not working, or is it working and I made some errors?
Normally I should get different values of the face-values (or not?) ...
And if I get the enthalpy equation. Where is the solver solving the temperature eqn or calculate it from the enthalpy?