[kane]

I want to know the wall heat flux of Enclosure case.
Internal medium is gray and scattering gas.
I have to consider a scattering model.


first, I add this code.

scatterModel  constantScatter;

constantScatterCoeffs
{
   sigma        sigma   [ 0 -1 0 0 0 0 0 ] 0.5;
   C            C       [ 0 0 0 0 0 0 0 ] 0;
}

But I think it was not effective.
What should I do additionally??

++++oh, Is is not possible to consider scattering in fvdom model????

[dlahaye]

What do you mean by "not effective"?

[kane]

oh I'm sorry.
''NOT effective'' means that I added the code, but it did not work.
As a result, i could not deal with scattering in fvdom model.

[kane]

I'm really confused...

Is it not possible to deal with scattering in fvdom model??

[dlahaye]

I image that fvDOM implemented in OpenFoam does allow for scattering. Could you please be more descriptive than "not working". Do you obtain results that are incorrect or no results at all?

[kane]

thank you for commenting, dlahaye!!

I'm sorry for my abstract explanation.

First I made a solver that deal with only radiation heat transfer.
So I applied this solver in Enclosed case.
As a result, I gained wall heat fluxes such as qr,qin in Openfoam.

But, this time I want to get wall heat fluxes(qr,qin) using a radiation solver
added to scattering. I added below codes in Enclosed case- file(constant/radiationproperties) to deal with scattering.

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      radiationProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

radiation on;

radiationModel  fvDOM;

fvDOMCoeffs
{
    nPhi    10;         // azimuthal angles in PI/2 on X-Y.(from Y to X)
    nTheta  10;          // polar angles in PI (from Z to X-Y plane)
    convergence 1e-6;   // convergence criteria for radiation iteration
    maxIter 1000;         // maximum number of iterations
}

// Number of flow iterations per radiation iteration
solverFreq 1;

absorptionEmissionModel constantAbsorptionEmission;

constantAbsorptionEmissionCoeffs
{
   absorptivity    absorptivity    [ 0 -1 0 0 0 0 0 ]  1;
   emissivity      emissivity      [ 0 -1 0 0 0 0 0 ]  1;
   E               E               [ 1 -1 -3 0 0 0 0 ] 0;
}

scatterModel  constantScatter;

constantScatterCoeffs
{
   sigma        sigma   [ 0 -1 0 0 0 0 0 ] 0.5;
   C            C       [ 0 0 0 0 0 0 0 ] 0;
}

But it was the same results(qin,qr) as no scattering.
I thought that i could only add scatterModel codes in radiationproperties to
deal with scattering. What should i do additionally??

I attached my solver and case file. Please check it out.
'Enclosescatter' is case file that deal with scattering.
'radiationFoam' is solver file.

[kane]

Does anyone who have ever tried to run fvdom model case including scattering????

[wangchengan2003]

Dear kane,
If you know how to use incident radition 'G' in the solver? Because I want to calculate the temperaute field by 'T=pow(fvDOM->G()/(4*sigma),0.25);'

Best regards,

Chengan

[keitaro7_14]

If you analyze fvDOM.C, you can see that scatterModel is not considered (although scatterModel.H is loaded). However, in P1 model the variable scatter_ is used.

[sante_junior]

Originally Posted by keitaro7_14

If you analyze fvDOM.C, you can see that scatterModel is not considered (although scatterModel.H is loaded). However, in P1 model the variable scatter_ is used.

that's right