[Marshak]

In fireFoam solver, the radiation source term is included in hs equation ..file YhsEqn.H as follows:

fvScalarMatrix hsEqn
(
fvm::ddt(rho, hs)
+ mvConvection->fvmDiv(phi, hs)
- fvm::laplacian(turbulence->alphaEff(), hs)
==
DpDt
+ dQ
+ radiation->Shs(thermo)
+ parcels.Sh(hs)
+ surfaceFilm.Sh()
);

I want to write out the radiation source term 'radiation->Shs(thermo)'. Can someone tell what changes should I make in the fireFoam solver to write the radiation source term.

[jherb]

You have to add a new field (probably in createFields.H) which is written every timestep. Then copy radiotion-Shs(thermo) to this field before the timestep is finished.

[zqlhzx]

Hi Marshak,
I have the some problem with you!I also want to write out the radiation source term 'radiation->Shs(thermo)'.Have you solved it?If it is solved,could tell how to do it?Thanks in advance!

[Marshak]

'jhreb' I tried with createfields.H but it is not working...

[jherb]

Here is what I did in a comparable case:
Add a new field to createFields.H

    volScalarField Shs(
        IOobject
        (
            "Shs",
            runTime.timeName(),
            mesh,
            IOobject::NO_READ,
            IOobject::AUTO_WRITE // this tells OpenFOAM to store the field after each timestep
        ),
        dimensionedScalar("zero", dimensionSet(0, 0, 1, 0, 0), 0.0) // here you have to set the correct dimensions of your field
    );

Then in YhsEqn.H, either after hsEqn.solve() or thermo.correct() add the following:

if (pimple.finalIter())
        {
            Shs = radiation->Shs(thermo);
        }

The whole code is not tested. But it should be a starting point.

[dzi]

Hello,
fyi I put a related radiation question to another thread
http://www.cfd-online.com/Forums/ope...blem-5.html#99
thanks dirk

[Marshak]

Shs() in radiation::radiationModel is defined as fvScalarMatrix. How can a fvScalarmatrix be written as a volScalarField?

[jherb]

You are right. I my case what I really used was something like

rHS1 = fvc::reconstruct
            (
                (
                  - ghf*fvc::snGrad(rho)
                ) * mesh.magSf()
            );

But here there is a quite simple solution (OpenFOAM version 2.2.2):

    if (pimple.finalIter())
    {
        Info<< radiation->Sh(thermo);
    }

I guess in your version its ->Shs(thermo). This will output the coefficients of the matrix. If there are only diagonal elements (as it looks like), then you can add ->D() to the command.

I didn't test it, but it is probably possible to assign these diagonal elements to an volume scalar field, which could be written to disk and then loaded into Paraview.

Originally Posted by Marshak

Shs() in radiation::radiationModel is defined as fvScalarMatrix. How can a fvScalarmatrix be written as a volScalarField?

[davidbarreiro]

Hi everyone!

I've got a similar problem with a modification of the kinematicSingleLayer model for reactingParcelFilFoam. I added a pressure term to model the surface tension and it's working quite well. The problem is that I want to visualize this term and I'm not being able to save it. I've tried to modify the original NO_WRITE for AUTO_WRITE but it doesn't work. Any idea??

tmp<volScalarField> myKinematicSingleLayer::pu()
{
    return tmp<volScalarField>
    (
        new volScalarField 
        (
            IOobject
            (
                typeName + ":pu",
                time_.timeName(),
                regionMesh(),
                IOobject::NO_READ,
                IOobject::AUTO_WRITE
            ),
            pPrimary_                  // pressure (mapped from primary region)
          - pSp_                           // accumulated particle impingement
          - fvc::laplacian(sigma_, delta_) // surface tension
          - sigma_/hDisj_*(cos(theta_/180*3.141593)-1)*exp(-delta_/hDisj_) //Disjoining pressure
        )
    );
}

Thank you in advance,

David

[Madeleine C]

Originally Posted by davidbarreiro

Hi everyone!

I've got a similar problem with a modification of the kinematicSingleLayer model for reactingParcelFilFoam. I added a pressure term to model the surface tension and it's working quite well. The problem is that I want to visualize this term and I'm not being able to save it. I've tried to modify the original NO_WRITE for AUTO_WRITE but it doesn't work. Any idea??

tmp<volScalarField> myKinematicSingleLayer::pu()
{
    return tmp<volScalarField>
    (
        new volScalarField 
        (
            IOobject
            (
                typeName + ":pu",
                time_.timeName(),
                regionMesh(),
                IOobject::NO_READ,
                IOobject::AUTO_WRITE
            ),
            pPrimary_                  // pressure (mapped from primary region)
          - pSp_                           // accumulated particle impingement
          - fvc::laplacian(sigma_, delta_) // surface tension
          - sigma_/hDisj_*(cos(theta_/180*3.141593)-1)*exp(-delta_/hDisj_) //Disjoining pressure
        )
    );
}

Thank you in advance,

David

Dear David,
Did you find a solution to this?
I am attempting to write out a source term and got stuck.