|
[Sponsors] |
July 27, 2012, 09:02 |
Finding calculation of T_ field
|
#1 |
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51 |
Hi all,
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: Code:
template<class MixtureType> void Foam::hPdfThermo<MixtureType>::correct() { // force the saving of the old-time values psi_.oldTime(); if (counter == propertyUpdate) { Info << "Updating look-up table extractions..." << endl; update(); counter = 0; } if (counter_mass_fractions == massFractionsUpdate) { Info << "Updating mass fraction extractions..." << endl; updateMassFractions(); counter_mass_fractions = 0; } Debug(); // added by myself calculate(); counter++; counter_mass_fractions++; } Code:
template<class MixtureType> void Foam::hPdfThermo<MixtureType>::Debug() { Info << "Debug" << endl; forAll(csi_.boundaryField(), patchi) { if(isA<wallHeatTransferFvPatchScalarField>(T_.boundaryField()[patchi])) { const fvPatchScalarField & pcsi = csi_.boundaryField()[patchi]; const fvPatchScalarField & pT = T_.boundaryField()[patchi]; forAll(pcsi, facei) { Info << "T = " << pT[facei] << endl; } } } } 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? Tobi |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
problems after decomposing for running | alessio.nz | OpenFOAM | 7 | March 5, 2021 04:49 |
Starting field averaging using libFunctionObject after certain time | eelcovv | OpenFOAM Programming & Development | 25 | December 7, 2015 22:28 |
how to impose fixed value at a point (or region) during calculation ? | Cyp | OpenFOAM Programming & Development | 25 | February 29, 2012 02:42 |
Phase field: calculation of the |grad f| div(grad f/|grad f|) term | Alucard | OpenFOAM | 0 | November 2, 2011 09:44 |
Need help evaluating field by interpolation with a "lookup table" | JohnRHawk | OpenFOAM Programming & Development | 2 | June 28, 2011 13:27 |