CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   nut wall functions for incompressible RAS (http://www.cfd-online.com/Forums/openfoam-solving/101591-nut-wall-functions-incompressible-ras.html)

stilljourney May 8, 2012 10:55

nut wall functions for incompressible RAS
 
Hi guys! I'm new to openFOAM, and it's my first thread here. :)

I'm trying to run an existed openfoam case (it is written in 1.6.0) using openFOAM 2.1.0. I'm using LaunderSharmaKE turbulence model. The 0/nut file has some entries as:

fixedWalls
{
type nutWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0;
}

I guess this is where the problem generates. When I run simpleFoam solver, it gives out the error message:

--> FOAM FATAL IO ERROR:
Unknown patchField type nutWallFunction for patch type wall

Valid patchField types are :

65
(
advective
atmBoundaryLayerInletEpsilon
buoyantPressure
calculated
codedFixedValue
codedMixed
cyclic
cyclicAMI
cyclicSlip
directionMixed
empty
epsilonWallFunction
fan
fanPressure
fixedFluxPressure
fixedGradient
fixedInternalValue
fixedPressureCompressibleDensity
fixedValue
freestream
freestreamPressure
inletOutlet
inletOutletTotalTemperature
kappatJayatillekeWallFunction
kqRWallFunction
mapped
mappedField
mappedFixedInternalValue
mappedFixedPushedInternalValue
mixed
nonuniformTransformCyclic
nutLowReWallFunction
nutTabulatedWallFunction
nutURoughWallFunction
nutUSpaldingWallFunction
nutUWallFunction
nutkRoughWallFunction
nutkWallFunction
omegaWallFunction
oscillatingFixedValue
outletInlet
outletMappedUniformInlet
partialSlip
processor
processorCyclic
rotatingTotalPressure
sliced
slip
symmetryPlane
syringePressure
timeVaryingMappedFixedValue
totalPressure
totalTemperature
turbulentHeatFluxTemperature
turbulentInlet
turbulentIntensityKineticEnergyInlet
turbulentMixingLengthDissipationRateInlet
turbulentMixingLengthFrequencyInlet
uniformDensityHydrostaticPressure
uniformFixedValue
uniformTotalPressure
waveSurfacePressure
waveTransmissive
wedge
zeroGradient
)

I checked $FOAM_SRC/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFuntions, there is a folder named nutWallFunctions, inside which list many wall functions.

So my question is which wall function I should use to replace the "nutWallFunction" entry. Would anyone give me some hints?

Thanks a lot!

uli July 24, 2012 19:04

hi all

I am facing exactly the same problem. I am trying to run a case, which works in OF1.7.x, in OF2.1.1 and get.

Code:

--> FOAM FATAL IO ERROR:
Unknown patchField type nutWallFunction for patch type wall

Valid patchField types are :

69
(
advective
atmBoundaryLayerInletEpsilon
buoyantPressure
calculated
codedFixedValue
codedMixed
cyclic
cyclicAMI
cyclicSlip
directionMixed
empty
epsilonWallFunction
fan
fanPressure
fixedFluxPressure
fixedGradient
fixedInternalValue
fixedPressureCompressibleDensity
fixedValue
freestream
freestreamPressure
inletOutlet
inletOutletTotalTemperature
kappatJayatillekeWallFunction
kqRWallFunction
mapped
mappedField
mappedFixedInternalValue
mappedFixedPushedInternalValue
mixed
multiphaseFixedFluxPressure
nonuniformTransformCyclic
nuSgsUSpaldingWallFunction
nutLowReWallFunction
nutTabulatedWallFunction
nutURoughWallFunction
nutUSpaldingWallFunction
nutUWallFunction
nutkAtmRoughWallFunction
nutkRoughWallFunction
nutkWallFunction
omegaWallFunction
oscillatingFixedValue
outletInlet
outletMappedUniformInlet
partialSlip
phaseHydrostaticPressure
processor
processorCyclic
rotatingTotalPressure
sliced
slip
symmetryPlane
syringePressure
timeVaryingMappedFixedValue
totalPressure
totalTemperature
turbulentHeatFluxTemperature
turbulentInlet
turbulentIntensityKineticEnergyInlet
turbulentMixingLengthDissipationRateInlet
turbulentMixingLengthFrequencyInlet
uniformDensityHydrostaticPressure
uniformFixedValue
uniformTotalPressure
waveSurfacePressure
waveTransmissive
wedge
zeroGradient
)


file: /home/myname/OpenFOAM/myname-2.1.1/run/sa/mycase/0/nut::boundaryField::circular from line 26 to line 30.

    From function fvPatchField<Type>::New(const fvPatch&, const DimensionedField<Type, volMesh>&, const dictionary&)
    in file /home/myname/OpenFOAM/OpenFOAM-2.1.1/src/finiteVolume/lnInclude/fvPatchFieldNew.C at line 135.

FOAM exiting

The case is about vortex shedding in the wake of a circular cylinder at Re=90000

Any suggestions regarding the wallfunction that should be used since "nutWallFunction" is no more available?

mturcios777 July 24, 2012 19:09

Try nutkWallFunction or nutUWallFunction; nutkWallFunction replicates OF 1.5.x behaviour. If you look at the code you can see that the different between them is the way the yplus is calculated. The doxygen documentation is your friend here:

http://foam.sourceforge.net/docs/cpp/a08745.html

uli July 24, 2012 20:09

Thanks for the quick response.

So the initial condition for nut is given in the same way:

Code:

    circular
    {
        type            nutkWallFunction;
        Cmu            0.09;
        kappa          0.41;
        E              9.8;
        value          uniform 0;
    }

?

mturcios777 July 25, 2012 12:03

Yes. When in doubt, find a working case from the tutorials that has something similar to what you want, then copy that.


All times are GMT -4. The time now is 20:41.