[vs1]

Dear All,

Why createFields.H for simpleFoam does not read transport properties dict and nu but for icoFoam it does?

The code below is present in createFields.H for icoFoam but not in simpleFoam. Why is it so ?

Info<< "Reading transportProperties\n" << endl;

IOdictionary transportProperties
(
IOobject
(
"transportProperties",
runTime.constant(),
mesh,
IOobject::MUST_READ,
IOobject::NO_WRITE
)
);

dimensionedScalar nu
(
transportProperties.lookup("nu")
)


Best,
Vaibhav

[MSF]

Hi,

simpleFoam uses transportModels and icoFoam doesn't. That is why in simpleFoam there is something like

 singlePhaseTransportModel laminarTransport(U, phi);

instead of the lines of code you posted. More about transportModels can be found in http://www.tfd.chalmers.se/~hani/kur...nFoam%20v2.pdf .

Best,
Moritz

[Andrea1984]

Hi,

I understand it's an old thread but I figured it is better to ask here before opening a new one.

My understanding is that when we are dealing with simpleFoam and pimpleFoam we can't control directly from createFields.h what is read in from transportProperties since this file is now read when we instantiate our transport model, i.e.:

singlePhaseTransportModel laminarTransport(U, phi);

This is different from icoFoam where we read the transportProperties directly from createFields.H.

I am adding a passive scalar transport equation to pimpleFoam. The passive scalar equation is:

    alphat = turbulence->nut()/Sct;
    alphat.correctBoundaryConditions();

    volScalarField alphaEff("alphaEff", turbulence->nu()/Sc + alphat);

    fvScalarMatrix TEqn
    (
        fvm::ddt(T)
      + fvm::div(phi, T)
      - fvm::laplacian(alphaEff, T)
    );

In order to specify the laminar (Sc) and turbulent Schmidt (Sct) numbers for the passive scalar I have created two DimensionedScalar entries in createFields.H,i.e.:

dimensionedScalar Sc
(
   "Sc",
   dimless,
   0.71
);

dimensionedScalar Sct
(
   "Sct",
   dimless,
   0.85
);

Everything is working fine with this approach.

However it is a bit inconvenient that every time I want to change the values of Sc and Sct I have to do it manually in createFields.H and re-compile my solver. Is there any way to read these two parameters from transportProperties without having to modify the singlePhaseTransportModel class?

Many thanks,
Andrea

[Andrea1984]

I have actually managed to do achieve this by using:

// Schmidt number
dimensionedScalar Sc("Sc", dimless, laminarTransport);

// Turbulent Schmidt number
dimensionedScalar Prt("Sct", dimless, laminarTransport);

after

singlePhaseTransportModel laminarTransport(U, phi);

in createFields.H

[lizzy]

Hi,

In this case, the value of Sc or Sct are given in the transportProperties just like nu? Please correct me if I'm wrong.

Thanks,
Ran

Originally Posted by Andrea1984

I have actually managed to do achieve this by using:

// Schmidt number
dimensionedScalar Sc("Sc", dimless, laminarTransport);

// Turbulent Schmidt number
dimensionedScalar Prt("Sct", dimless, laminarTransport);

after

singlePhaseTransportModel laminarTransport(U, phi);

in createFields.H

[Andrea1984]

Yes, the only difference is that you don't have to specify a dimension set because both Sc and Sct are dimless.

Andrea