Turbulence Model

Submitted by rqwang on Thu, 04/14/2011 - 20:35

Hi cfdemers,

I tried to change the turbulence model from RAS to LES in ErgunTestMPI.

I got error message like below:

// run_parallel_cfdemSolverPiso_ErgunTestMPI_CFDDEM_shared //

/home/ruo/OpenFOAM/ruo-1.6/run/PFliggghts/CFD

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 1.6-53b7f692aa41
Exec : cfdemSolverPiso_shared -parallel
Date : Apr 14 2011
Time : 14:26:56
Host : ruo-ubuntu
PID : 7759
Case : /home/ruo/OpenFOAM/ruo-1.6/run/PFliggghts/CFD
nProcs : 2
Slaves :
1
(
ruo-ubuntu.7760
)

Pstream initialized with:
floatTransfer : 0
nProcsSimpleSum : 0
commsType : nonBlocking
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

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

Create mesh for time = 0

Reading field p

Reading physical velocity field U
Note: only if voidfraction at boundary is 1, U is superficial velocity!!!

Reading momentum exchange field Ksl

Reading voidfraction field voidfraction = (Vgas/Vparticle)

Creating dummy density field rho = 1

Reading particle velocity field Us

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type LESModel
Selecting LES turbulence model dynMixedSmagorinsky
scaleSimilarityCoeffs
{
filter simple;
ce 1.048;
}
dynSmagorinskyCoeffs
{
filter simple;
ce 1.048;
}

dynMixedSmagorinskyCoeffs
{
filter simple;
ce 1.048;
}

Reading g
[1]
[1]
[1]
request for turbulenceModel RASProperties from objectRegistry region0 failed
available objects of type turbulenceModel are
(
LESProperties
)
#0 Foam::error::printStack(Foam::Ostream&) in "/home/ruo/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/ruo/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Foam::Ostream& Foam::operator<< (Foam::Ostream&, Foam::errorManip) in "/home/ruo/OpenFOAM/ruo-1.6/applications/bin/linuxGccDPOpt/cfdemSolverPiso_shared"
#3 Foam::incompressible::turbulenceModel const& Foam::objectRegistry::lookupObject(Foam::word const&) const in "/home/ruo/OpenFOAM/ruo-1.6/lib/linuxGccDPOpt/liblagrangianCFDEM_shared.so"
#4 Foam::cfdemCloud::cfdemCloud(Foam::fvMesh const&) in "/home/ruo/OpenFOAM/ruo-1.6/lib/linuxGccDPOpt/liblagrangianCFDEM_shared.so"
#5
in "/home/ruo/OpenFOAM/ruo-1.6/applications/bin/linuxGccDPOpt/cfdemSolverPiso_shared"
#6 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#7
in "/home/ruo/OpenFOAM/ruo-1.6/applications/bin/linuxGccDPOpt/cfdemSolverPiso_shared"
[1]
[1]
[1] From function objectRegistry::lookupObject(const word&) const
[1] in file /home/ruo/OpenFOAM/OpenFOAM-1.6/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 140.
[1]
FOAM parallel run aborting
[1]
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 1 in communicator MPI_COMM_WORLD
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------
[0]
[0]
[0]
request for turbulenceModel RASProperties from objectRegistry region0 failed
available objects of type turbulenceModel are

1
(
LESProperties
)
#0 Foam::error::printStack(Foam::Ostream&) in "/home/ruo/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/ruo/OpenFOAM/OpenFOAM-1.6/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Foam::Ostream& Foam::operator<< (Foam::Ostream&, Foam::errorManip) in "/home/ruo/OpenFOAM/ruo-1.6/applications/bin/linuxGccDPOpt/cfdemSolverPiso_shared"
#3 Foam::incompressible::turbulenceModel const& Foam::objectRegistry::lookupObject(Foam::word const&) const in "/home/ruo/OpenFOAM/ruo-1.6/lib/linuxGccDPOpt/liblagrangianCFDEM_shared.so"
#4 Foam::cfdemCloud::cfdemCloud(Foam::fvMesh const&) in "/home/ruo/OpenFOAM/ruo-1.6/lib/linuxGccDPOpt/liblagrangianCFDEM_shared.so"
#5
in "/home/ruo/OpenFOAM/ruo-1.6/applications/bin/linuxGccDPOpt/cfdemSolverPiso_shared"
#6 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#7
in "/home/ruo/OpenFOAM/ruo-1.6/applications/bin/linuxGccDPOpt/cfdemSolverPiso_shared"
[0]
[0]
[0] From function objectRegistry::lookupObject(const word&) const
[0] in file /home/ruo/OpenFOAM/OpenFOAM-1.6/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 140.
[0]
FOAM parallel run aborting
[0]
--------------------------------------------------------------------------
mpirun has exited due to process rank 1 with PID 7760 on
node ruo-ubuntu exiting without calling "finalize". This may
have caused other processes in the application to be
terminated by signals sent by mpirun (as reported here).
--------------------------------------------------------------------------

Could anyone help me?

Thanks.

rq

cgoniva's picture

cgoniva | Fri, 04/15/2011 - 08:15

Hi Rq,

Solver and testcase were neither designed nor tested for LES modelling of the flow.

Try to run the CFD part of the testcase with a pure CFD solver (without coupling) e.g. pisoF...
check whether all dictioanry entries are provided
do this in serial to get clearer output messages
try again with the coupled solver

please keep us posted about your progress.

Cheers,
Chris

rqwang | Fri, 04/15/2011 - 16:02

Hi Chris,

Thank you for your reply.

After setting up the default divScheme to Gaussian linear, the CFD part can run without error.

However, the same error message occurred if I run the coupled solver again.

I suspect in cfdemCloud.C the solve still tries to read the "RASProperties" (See the line of 75).

That looks what the error message indicates. Any idea?

Thanks.

Rq

rqwang | Mon, 04/18/2011 - 20:50

Hi,

To update, I changed "RASProperties" to "LESProperties" in cfdemCloud.C and re-compiled.

Also the "fix couple" command is changed to "fix cfd all couple/cfd/force every mpi".

LES seems to work.

Thanks.

Rq

cgoniva's picture

cgoniva | Tue, 04/19/2011 - 09:20

Hi Rq,

that is excellent news. Thank you for your update!
I will try to figure out if I can include that optional turbulence model to the next update!

Cheers,
Christoph

cgoniva's picture

cgoniva | Tue, 04/19/2011 - 13:27

Hi Rq,

I did some changes to the source files - now switching between RAS and LES should be straight forward without recompiling.

You need to set the switches in CFD/constant/turbulenceProperties as well as in CFD/constant/couplingProperties.

Please give it a try and keep us posted about the results.

A real test case for that feature would be nice!

Cheers,
Chris

rqwang | Fri, 04/22/2011 - 21:55

Hi Chris,

I tested the LES model in the updated version.

It works! Thank you for your quick update!

rq