[vcattoir]

Hello Foamers,

I am currently trying to simulate a compressible supersonic flow in a nozzle using LES and with rhoCentralFoam as a solver. I am using oneEqEddy for the LES and have been trying to use wall functions but not successfully for the moment. I initialize my fields with result of RANS simulation on the same meshing (i refined it more and used mapFields).

Whatever I do, i end up with this error :

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.x                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.3.x-6b18d98df2bc
Exec   : rhoCentralFoam
Date   : Jan 09 2017
Time   : 13:07:10
Host   : "lemaitre2.cism.ucl.ac.be"
PID    : 20285
Case   : /home/users/v/c/vcattoir/Nozzle
nProcs : 1
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0.003

Reading thermophysical properties

Selecting thermodynamics package 
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       sutherland;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Reading field U

Creating turbulence model

Selecting turbulence model type LESModel
Selecting LES turbulence model oneEqEddy
Selecting LES delta type smooth
Selecting LES delta type cubeRootVol


--> FOAM FATAL ERROR: 

    request for volScalarField mut from objectRegistry region0 failed
    available objects of type volScalarField are

14
(
thermo:mu
thermo:psi
rho
k
delta
thermo:psi_0
p
T
rhoE
e
alphaSgs
geometricDelta
muSgs
thermo:alpha
)


    From function objectRegistry::lookupObject<Type>(const word&) const
    in file /home/nl/OpenFOAM/OpenFOAM-2.3.x/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 198.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) in "/usr/local/centFOAM/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::error::abort() in "/usr/local/centFOAM/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const& Foam::objectRegistry::lookupObject<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> >(Foam::word const&) const in "/usr/local/centFOAM/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#3  Foam::compressible::alphatWallFunctionFvPatchScalarField::updateCoeffs() in "/usr/local/centFOAM/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libcompressibleRASModels.so"
#4  Foam::fvPatchField<double>::evaluate(Foam::UPstream::commsTypes) in "/usr/local/centFOAM/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/bin/rhoCentralFoam"
#5  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::evaluate() in "/usr/local/centFOAM/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/bin/rhoCentralFoam"
#6  Foam::compressible::LESModels::oneEqEddy::updateSubGridScaleFields() in "/usr/local/centFOAM/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libcompressibleLESModels.so"
#7  Foam::compressible::LESModels::oneEqEddy::oneEqEddy(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::fluidThermo const&, Foam::word const&, Foam::word const&) in "/usr/local/centFOAM/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libcompressibleLESModels.so"
#8  Foam::compressible::LESModel::adddictionaryConstructorToTable<Foam::compressible::LESModels::oneEqEddy>::New(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::fluidThermo const&, Foam::word const&) in "/usr/local/centFOAM/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libcompressibleLESModels.so"
#9  Foam::compressible::LESModel::New(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::fluidThermo const&, Foam::word const&) in "/usr/local/centFOAM/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libcompressibleLESModels.so"
#10  Foam::compressible::turbulenceModel::addturbulenceModelConstructorToTable<Foam::compressible::LESModel>::NewturbulenceModel(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::fluidThermo const&, Foam::word const&) in "/usr/local/centFOAM/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libcompressibleLESModels.so"
#11  Foam::compressible::turbulenceModel::New(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::fluidThermo const&, Foam::word const&) in "/usr/local/centFOAM/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/lib/libcompressibleTurbulenceModel.so"
#12  
 in "/usr/local/centFOAM/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/bin/rhoCentralFoam"
#13  __libc_start_main in "/lib64/libc.so.6"
#14  
 in "/usr/local/centFOAM/OpenFOAM/OpenFOAM-2.3.x/platforms/linux64GccDPOpt/bin/rhoCentralFoam"
Abandon (core dumped)

Which i understand is incorrect use of wall function because it looks for a mut when only muSgs is available but I am currently trying without any wall functions and all 0 files on walls are set at calculated. I want to begin without any wall function and add them to be sure it still works but could not get the most simple case to at least begin.

I may be missing a thing or two but cannot find any reliable documentation on wall functions on LES and how to choose them.

find $FOAM_SRC/turbulenceModels -name wallFunctions

only outputs:

/usr/local/centFOAM/OpenFOAM/OpenFOAM-2.3.x/src/turbulenceModels/incompressible/RAS/backwardsCompatibility/wallFunctions
/usr/local/centFOAM/OpenFOAM/OpenFOAM-2.3.x/src/turbulenceModels/incompressible/RAS/derivedFvPatchFields/wallFunctions
/usr/local/centFOAM/OpenFOAM/OpenFOAM-2.3.x/src/turbulenceModels/compressible/RAS/backwardsCompatibility/wallFunctions
/usr/local/centFOAM/OpenFOAM/OpenFOAM-2.3.x/src/turbulenceModels/compressible/RAS/derivedFvPatchFields/wallFunctions

Nothing for LES, only RANS

Could any of you help me, it would be greatly appreciated ?

Thank you

Vcattoir

[Hillie]

Hi Vcattoir,

Made any progress already?
The wallfunctions for LES are the same ones you will use for RANS. So you just have to use the eddy viscosity wall model that is available for compressible flows.

Just out of interest, what are you trying to to with your nozzle simulation?

[vcattoir]

Thank you for your answer. I am still blocked at the same error and have no clue how to find where and why it looks for mut.

I am doing this LES simulation in preparation for the simulation I will have to do for my master these which is a 3D LES simulation of a supersonic ejector. I do that because LES with OpenFOAM has never been tried in my departement and we wanted to try an easy and well known case before the real one.

Vcattoir

[Hillie]

hey,

Cool, I am working on nozzle flows as well for my PhD. I am quite a similar position, since I am the first one in our group to get my head around some LES simulations.

Regarding your error, somewhere you are trying to request your mut field.
If you are still having issue, I am happy to have a quick look over it if you want to.

[vcattoir]

Hi,
I finally found my mistake after nearly 2 months !
I still have a problem using wall functions. I do not understand what you mean by the same one that in RAS, i used mut in RANS, not muSgs and so on. For example, I am now trying to use muSgsUSpaldingWallFunction which is absent from the Valid patchField types OpenFOAM gives me.
I have not succeeded to find an example on the internet and do not understand why I would not work.

PS: I am using OpenFOAM 2.3

[Hillie]

Hi,

I am not to sure. I started working with OpenFoam after V3.0 where all the wall functions have been integrated into a nut field. So I am not to sure why it doesn't work, sorry.

[vcattoir]

I have done a first version that works but would like to try yours using OF3. What do you put in your nut, i thought than in compressible, mut was concerned, not nut. I tried with simple "zeroGradient" but have no clue what it should be.

Thanks for the reply

vcattoir

[Hillie]

If you use mut or nut is not that important I guess, since it is just a normalization with density.
If you want to run without a wall model the eddy viscosity should be zero at the wall, i.e just a fixed value BC.

For the wall functions, till now I have been running with nutUWallFunction.
I haven't run a case yet with a high enough wall resolution to asses if this generates a good boundary layer profile though.