[boffin5]

In a previous post regarding how to specify a temperature in a porous zone, TEMPEST answers his own question as follows. But as a relative newbie, I don't understand how to implement the TEqn.H file. Where does it go, and are there any programming changes that need to accompany it?

Originally Posted by Tempest

We can do this by using fvOptions and

TEqn.H

{
    fvScalarMatrix TEqn
    (
        fvm::ddt(T)
      + (1/porosity)*fvm::div(phi, T)
      - fvm::laplacian(nu/Pr, T)
    ==
    fvOptions(T)
    );

    TEqn.relax();

    fvOptions.constrain(TEqn);

    TEqn.solve();

    fvOptions.correct(T);

    rhok = 1.0 - beta*(T - TRef);
}

And in fvOptions we can add

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

s1
{
    type            scalarFixedValueConstraint;
    active          true;


    scalarFixedValueConstraintCoeffs
    {
        selectionMode   cellZone;
            cellZone         porousZone;
        fieldValues
        {
            T              350; \\the fixed temperature value
        }
    }
}


// ************************************************************************* //

[petros]

Hi,

All files ending in "<>Eqn.H" contain source code that defines a mathematical model that the solver will solve. All these files are included and compiled in the solver's application folder. For example, see buoyantBoussinesqSimpleFoam and the respective TEqn.H file, which you can use as a template for your own implementation.

Best,
Petros

[boffin5]

Thanks Petros!
From my quick look at buoyantBoussinesqSimpleFoam, it seems like it address liquids, not air. My need is to simulate airflow through a heated porous media.
Am I incorrect in my observation about the solver, or is there another more appropriate one?

[petros]

Hi,

Air is a fluid itself so it makes no difference.

I don't know what you want to simulate exactly, not to mention what assumptions your case involves, but you can have a brief description of every solver if you search for solvers dedicating to heat transfer. See here for example for heat transfer solvers in v2012.