[ssss]

Hi,

I've implemented a scalar Transport equation to pisoFoam, but I cannot obtain some values from it

Let's start with it:

In createFields.H I add the following:

   
Info<< "Reading field alphaT\n" << endl;
volScalarField alphaT
(
    IOobject
        (
            "alphaT",
            runTime.timeName(),
            mesh,
            IOobject::MUST_READ,
            IOobject::AUTO_WRITE
        ),
        mesh
);

dimensionedScalar prandtlTurbulent(laminarTransport.lookup("prandtlTurbulent"));

singlePhaseTransportModel laminarTransport(U, phi);

    autoPtr<incompressible::turbulenceModel> turbulence
    (
        incompressible::turbulenceModel::New(U, phi, laminarTransport)
    );

In my_pisoFoam.C

#include "turbulenceModel.H"
#include "createFields.H"
...
int main(){
...
        alphaT=turbulence->nut()/prandtlTurbulent;
        alphaT.correctBoundaryConditions();
...
}

In 0/alphaT

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

dimensions      [0 2 -1 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    inlet 
    {
        type            calculated; //search for it
        //value           uniform 0;

    }
    outlet
    {
        type            calculated; //search for it
        //value           uniform 0;
    }
    channelWall
    {
        type            zeroGradient;
    }
    symmetry
    {
        type            symmetryPlane;
    }

    frontAndBack
    {
        type            empty;
    }
}


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

In constant/transportProperties

prandtlTurbulent prandtlTurbulent [0 0 0 0 0 0 0] 0.85;

But whenever I open alphaT in every folderIteration I only see:

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

dimensions      [0 2 -1 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    channelWall
    {
        type            zeroGradient;
    }
    symmetry
    {
        type            symmetryPlane;
    }
    inlet
    {
        type            calculated;
        value           uniform 0;
    }
    outlet
    {
        type            calculated;
        value           uniform 0;
    }
    frontAndBack
    {
        type            empty;
    }
}


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

It only repeats the alphaT which is in 0/ . It doesn write anything.

Maybe it is because I choose this?

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      RASProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

RASModel        laminar;

turbulence      on;

printCoeffs     on;

RASModel = laminar.

Thank you very much.

[ssss]

I found the problem, if I choose the turbulence as laminar, he would assume that nut is 0. I will try to see how could I make a dns simulation with pisoFoam and then get nut somewhere

Thank you