[mabinty]

dear all,

i created a utility for post processing where i manipulate variables of the RAS turbulence model (in OF-1.7.1). when i apply it to a case using the "turbulentIntensityKineticEnergyInlet" BC for k, the utility crashes at the line

const volSymmTensorField RS = RASModel->R();

with the error message:

--> FOAM FATAL ERROR: 
Unknown patchTypefield type turbulentIntensityKineticEnergyInlet

Valid patchField types are :

30
(
advective
calculated
compressible::kqRWallFunction
cyclic
directMapped
directionMixed
empty
fixedGradient
fixedInternalValue
fixedMeanValue
fixedValue
freestream
inletOutlet
mixed
oscillatingFixedValue
outletInlet
outletMappedUniformInlet
partialSlip
processor
sliced
slip
symmetryPlane
timeVaryingMappedFixedValue
timeVaryingUniformFixedValue
timeVaryingUniformInletOutlet
turbulentInlet
uniformFixedValue
waveTransmissive
wedge
zeroGradient
)


    From function fvPatchField<Type>::New(const word&, const fvPatch&, const DimensionedField<Type, volMesh>&)
    in file /home/opencfd/OpenFOAM/OpenFOAM-1.7.x/src/finiteVolume/lnInclude/newFvPatchField.C at line 53.

FOAM exiting

i had a look at the definition of R (src/turbulenceModels/compressible/RAS/kEpsilon/kEpsilon.C):

tmp<volSymmTensorField> kEpsilon::R() const
{
    return tmp<volSymmTensorField>
    (
        new volSymmTensorField
        (
            IOobject
            (
                "R",
                runTime_.timeName(),
                mesh_,
                IOobject::NO_READ,
                IOobject::NO_WRITE
            ),
            ((2.0/3.0)*I)*k_ - (mut_/rho_)*dev(twoSymm(fvc::grad(U_))),
            k_.boundaryField().types()
        )
    );
}

i assume (cause i don t really know what the last line in the IOobject definition means), that through the line "k_.boundaryField().types()" R checks k BCs and cannot find "turbulentIntensityKineticEnergyInlet". i cannot really understand this behaviour as "turbulentIntensityKineticEnergyInlet" is part of the finiteVolume lib, which i included in the option file of my utility and is also present in the option file of the turbulence model lib.

does anybody know what i m doing wrong here?

thanks in adcance!

beste regards,
aram