CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

Finding calculation of T_ field

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 27, 2012, 09:02
Default Finding calculation of T_ field
  #1
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
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++;
}
And in the Debug-Function I am getting the face-values:

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;
                        }
                }
        }
}
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?


Tobi
Tobi is offline   Reply With Quote

Old   July 27, 2012, 09:38
Default
  #2
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi,

its me again. I find it out. There is no calculation of the scalar T field.
Well ... thanks for readying.

Tobi
Tobi is offline   Reply With Quote

Reply


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


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


All times are GMT -4. The time now is 03:09.