Peclet number for multiphase systems
 

Dear Foamers,

I am simulating a gas-liquid system using reactingMultiphaseEulerFoam and I would like to have the Peclet number for each phase or only one phase (gas).

I included the following function object in the controlDict:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
peclet
{
type PecletNo;
libs ("libfieldFunctionObjects.so");
writeControl adjustableRunTime;
writeInterval 1e-6;
log false;
}

// ************************************************** *********************** //


But I got the following error:

"--> FOAM FATAL ERROR:

request for turbulenceModel turbulenceProperties from objectRegistry region0 failed
available objects of type turbulenceModel are

2
(
turbulenceProperties.liquid
turbulenceProperties.gas
)


From function const Type& Foam::objectRegistry::lookupObject(const Foam::word&) const [with Type = Foam::turbulenceModel]
in file /home/ubuntu/OpenFOAM/OpenFOAM-6/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 193.

FOAM aborting

#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::turbulenceModel const& Foam::objectRegistry::lookupObject<Foam::turbulenc eModel>(Foam::word const&) const at ??:?
#3 Foam::functionObjects::PecletNo::calc() at ??:?
#4 Foam::functionObjects::fieldExpression::execute() at ??:?
#5 Foam::functionObjectList::execute() at ??:?
#6 ? in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/reactingMultiphaseEulerFoam"
#7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8 ? in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/reactingMultiphaseEulerFoam"
Aborted (core dumped)
"

I was wondering if you could help me solve this problem.

Kind Regards,
Sut

You request for a field which doesnt exist, can you share the other dicts as well?

Dear HPE,
Thanks for your reply.
Could you please tell me which dictionary do you mean?

I guess the problem is with being multiphase. I tried to use pisoFoam solver and changed *.gas to *, e.g. U.gas ---> U and this worked. So, I guess the problem is that or 1) the PecletNo function object is unable to work with multiphase or 2) I should insert some additional lines in the peclet dict to ask it to choose only one phase, e.g. gas or liquid.

Anyway, please tell me which dictionaries do you want to see.

Kind Regards,
Sut