CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

LES of a nozzle with rhoCentralFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   January 9, 2017, 08:28
Default LES of a nozzle with rhoCentralFoam
  #1
New Member
 
-
Join Date: Oct 2016
Posts: 3
Rep Power: 2
vcattoir is on a distinguished road
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 :

Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  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.
Code:
find $FOAM_SRC/turbulenceModels -name wallFunctions
only outputs:
Code:
/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
vcattoir is offline   Reply With Quote

Old   February 8, 2017, 18:54
Default
  #2
New Member
 
Hilbert
Join Date: Aug 2015
Location: Australia
Posts: 29
Rep Power: 3
Hillie is on a distinguished road
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?
Hillie is offline   Reply With Quote

Old   February 10, 2017, 07:53
Default
  #3
New Member
 
-
Join Date: Oct 2016
Posts: 3
Rep Power: 2
vcattoir is on a distinguished road
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
vcattoir is offline   Reply With Quote

Old   February 12, 2017, 18:37
Default
  #4
New Member
 
Hilbert
Join Date: Aug 2015
Location: Australia
Posts: 29
Rep Power: 3
Hillie is on a distinguished road
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.
Hillie is offline   Reply With Quote

Old   February 15, 2017, 12:37
Default
  #5
New Member
 
-
Join Date: Oct 2016
Posts: 3
Rep Power: 2
vcattoir is on a distinguished road
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
vcattoir is offline   Reply With Quote

Old   February 15, 2017, 18:18
Default
  #6
New Member
 
Hilbert
Join Date: Aug 2015
Location: Australia
Posts: 29
Rep Power: 3
Hillie is on a distinguished road
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.
Hillie is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
LES problems for nozzle geometry Bogdan90 OpenFOAM Running, Solving & CFD 0 December 7, 2015 07:41
calculate a nozzle flow using rhoCentralFoam hg2lf OpenFOAM 0 October 25, 2012 21:26
Always crash when solve a C-D nozzle flow field using rhoCentralFoam hawklion OpenFOAM Running, Solving & CFD 0 March 9, 2011 07:13
Always crash when solve a C-D nozzle flow field using rhoCentralFoam hawklion OpenFOAM 3 March 8, 2011 20:03
compressible flow in a counterflow nozzle d.vamsidhar FLUENT 0 November 24, 2005 02:45


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