CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Programming & Development (https://www.cfd-online.com/Forums/openfoam-programming-development/)
-   -   Unknown patchTypefield "turbulentIntensityKineticEnergyInlet" (https://www.cfd-online.com/Forums/openfoam-programming-development/99273-unknown-patchtypefield-turbulentintensitykineticenergyinlet.html)

mabinty March 30, 2012 06:37

Unknown patchTypefield "turbulentIntensityKineticEnergyInlet"
 
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

Code:

const volSymmTensorField RS = RASModel->R();
with the error message:

HTML Code:

--> 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):

Code:

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


All times are GMT -4. The time now is 15:59.