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

Possible error in the BC of the wideBandDiffusive radiation model

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 27, 2017, 05:13
Default Possible error in the BC of the wideBandDiffusive radiation model
  #1
New Member
 
Adrien
Join Date: Apr 2013
Posts: 2
Rep Power: 0
Heyheymymy is on a distinguished road
Hi Foamers,

I would like to have feedback on the following problem I found using the fvDOM radiation model with the wideBandAbsorptionDiffusive model.

The BCs for solving the RTE are calculated at the wall through the wideBandDiffusiveRadiation BC.

This BC uses the radiation emitted by the wall. This is calculated using the emissivity of the wall (temissivity in the code below) and black body radiation at the temperature of the wall (Eb).

Code:
            refValue()[faceI] =
                (
                    Ir[faceI]*(1.0 - temissivity[faceI])
                  + temissivity[faceI]*Eb[faceI]
                )/pi;

            // Emmited heat flux from this ray direction
            Qem[faceI] = refValue()[faceI]*nAve[faceI]
In my understanding of the band model, Eb should be the fraction of the radiation emitted by the blackbody in the wavelength range of the considered band.

However it seems that the BC is instead using the full blackbody emissive power.
I related this to the fact that the blackBodyEmission.C performs a loop on the whole internal temperature but not on its the boundaries.

A simple test for that was to add a counter

Code:
    
else
    { 
   int j=0;
        forAll(T, i)
        { 
            j+=1;
            scalar T1 = fLambdaT(band[1]*T[i]);
            scalar T2 = fLambdaT(band[0]*T[i]);
            dimensionedScalar fLambdaDelta
            (
                "fLambdaDelta",
                dimless,
                T1 - T2
            );
            Eb()[i] = Eb()[i]*fLambdaDelta.value();
Info << "j=" << j << endl;
        }
        return Eb;
    }
It appeared that this doesn't loop in the temperature BC, thus Eb remains equal to the initialized value that is sigma T.

It seems to me that this is problematic for the BC, am I right?

Thank you for your help,

Adrien
Heyheymymy is offline   Reply With Quote

Reply

Tags
boundary condition error, radiation fvdom, wideband

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Overflow Error in Multiphase Modelling with Two Continuous Fluids ashtonJ CFX 6 August 11, 2014 14:32
How model radiation of both solid wall and fluid£¿ Harry Qiu FLUENT 2 February 4, 2013 00:04
problem in glass radiation model (DO model) wunan805 FLUENT 0 October 10, 2012 05:06
Help: Radiation model MANOJKUMAR FLUENT 4 November 24, 2005 03:45
definition difficulty-->DO radiation model Harry Qiu FLUENT 0 March 29, 2001 09:19


All times are GMT -4. The time now is 10:34.