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

constantEmissivity option in view factor radiation model

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 27, 2022, 13:30
Default constantEmissivity option in view factor radiation model
  #1
New Member
 
Join Date: Mar 2022
Posts: 7
Rep Power: 4
ngbidi is on a distinguished road
Hi Foamers,

Usually I set the option "constantEmissivity" in the view factor radiation model to true because I used the lookup model with constant emissivity in the past. Recently I implemented a submodel for boundaryProperties such that the emissivity changes over time so I started to look more into this option.
Apparently the emissivity still changes over time when setting constantEmissivity to true. This is because the function to update the emissivity is called regardless of the option:

In radiation::viewFactor::calculate()
Code:
 const tmp<scalarField> teb = boundaryRadiation.emissivity(patchID, bandI);
constantEmissivity is then used later to do some calculations:

Also in radiation::viewFactor::calculate()
Code:
            if (!constEmissivity_)
            {
                scalarSquareMatrix C(totalNCoarseFaces_, 0.0);

                for (label i=0; i<totalNCoarseFaces_; i++)
                {
                    for (label j=0; j<totalNCoarseFaces_; j++)
                    {
                        const scalar invEj = 1.0/E[j];
                        const scalar sigmaT4 =
                            physicoChemical::sigma.value()*T4[j];

                        if (i==j)
                        {
                            C(i, j) = invEj - (invEj - 1.0)*Fmatrix_()(i, j);
                            q[i] +=
                                (Fmatrix_()(i, j) - 1.0)*sigmaT4 + qrExt[j];
                        }
                        else
                        {
                            C(i, j) = (1.0 - invEj)*Fmatrix_()(i, j);
                            q[i] += Fmatrix_()(i, j)*sigmaT4;
                        }

                    }
                }

                Info<< "Solving view factor equations for band :"
                    << bandI << endl;

                // Negative coming into the fluid
                LUsolve(C, q);
            }
            else //Constant emissivity
            {
                // Initial iter calculates CLU and caches it
                if (iterCounter_ == 0)
                {
                    for (label i=0; i<totalNCoarseFaces_; i++)
                    {
                        for (label j=0; j<totalNCoarseFaces_; j++)
                        {
                            const scalar invEj = 1.0/E[j];
                            if (i==j)
                            {
                                CLU_()(i, j) =
                                    invEj-(invEj-1.0)*Fmatrix_()(i, j);
                            }
                            else
                            {
                                CLU_()(i, j) = (1.0 - invEj)*Fmatrix_()(i, j);
                            }
                        }
                    }

                    DebugInFunction << "\nDecomposing C matrix..." << endl;

                    LUDecompose(CLU_(), pivotIndices_);
                }

                for (label i=0; i<totalNCoarseFaces_; i++)
                {
                    for (label j=0; j<totalNCoarseFaces_; j++)
                    {
                        const scalar sigmaT4 =
                            constant::physicoChemical::sigma.value()*T4[j];

                        if (i==j)
                        {
                            q[i] +=
                            (Fmatrix_()(i, j) - 1.0)*sigmaT4  + qrExt[j];
                        }
                        else
                        {
                            q[i] += Fmatrix_()(i, j)*sigmaT4;
                        }
                    }
                }


                Info<< "Solving view factor equations for band : "
                    << bandI  << endl;


                LUBacksubstitute(CLU_(), pivotIndices_, q);
                iterCounter_ ++;
            }
I don't understand what is going on in this part. When testing both options, qr is twice as high when setting constantEmissivity to false.
Does anyone know more about this?

Thanks in advance!
ngbidi is offline   Reply With Quote

Reply

Tags
emissivity, radiation model, view factor model

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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
Error - Solar absorber - Solar Thermal Radiation MichaelK CFX 12 September 1, 2016 06:15
Ansys CFX problem: unexpected very high temperatures in premix laminar combustion faizan_habib7 CFX 4 February 1, 2016 18:00
Wrong flow in ratating domain problem Sanyo CFX 17 August 15, 2015 07:20
Waterwheel shaped turbine inside a pipe simulation problem mshahed91 CFX 3 January 10, 2015 12:19
Water subcooled boiling Attesz CFX 7 January 5, 2013 04:32


All times are GMT -4. The time now is 05:57.