[mbuergle]

Dear OpenFoam Users,

I am trying to run a case with a porousBafflePressure BC with the twoPhaseEulerFoam solver (OpenFOAM v1712 on Ubuntu 16.04). However, when running the case I get the error message posted below. My search in this forum und Google have not lead to a solution so far. With interFoam the geometric setup with the porousBafflePressure BC worked.

I don't know if this information is of any help, but i got a very similar error message when i tried to implement the turbulentIntensityKineticEnergyInlet in Kombination mit turbulentMixingLengthDissipationRateInlet for k.air/k.water and epsilon.air/epsilon.water respectively.

Thank you in advance for your help!
Best regards,
Matthias

Steps to Reproduce:

1. blockMesh
2. topoSet
3. createBaffles -overwrite
4. setFields
5. twoPhaseEulerFoam | tee log

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1712                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : v1712
Arch   : "LSB;label=32;scalar=64"
Exec   : twoPhaseEulerFoam
Date   : May 31 2018
Time   : 13:34:38
Host   : "7a4463c6b203"
PID    : 118875
I/O    : uncollated
Case   : /home/mbuergle/openfoam/bottomoutlet/tPEF/TestCase1
nProcs : 1
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0


PIMPLE: max iterations = 5
    field U	: relTol 0, tolerance 1e-05
    field T	: relTol 0, tolerance 0.001
    field p_rgh	: relTol 0, tolerance 1e-05


Reading g

Reading hRef
Creating twoPhaseSystem

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

Calculating face flux field phi.air
Selecting diameterModel for phase air: constant
Selecting turbulence model type RAS
Selecting RAS turbulence model continuousGasKEpsilon
RAS
{
    RASModel        continuousGasKEpsilon;
    turbulence      on;
    printCoeffs     on;
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              0;
    sigmak          1;
    sigmaEps        1.3;
    alphaInversion  0.7;
}

Selecting thermodynamics package 
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectFluid;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Calculating face flux field phi.water
Selecting diameterModel for phase water: constant
Selecting turbulence model type RAS
Selecting RAS turbulence model continuousGasKEpsilon
RAS
{
    RASModel        continuousGasKEpsilon;
    turbulence      on;
    printCoeffs     on;
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              0;
    sigmak          1;
    sigmaEps        1.3;
    alphaInversion  0.7;
}

Selecting default blending method: linear
Selecting aspectRatioModel for (air in water): constant
Selecting aspectRatioModel for (water in air): constant
Selecting dragModel for (air and water): segregated
Selecting swarmCorrection for (air and water): none
Selecting dragModel for (air in water): SchillerNaumann
Selecting swarmCorrection for (air in water): none
Selecting dragModel for (water in air): SchillerNaumann
Selecting swarmCorrection for (water in air): none
Selecting virtualMassModel for (air in water): constantCoefficient
Selecting virtualMassModel for (water in air): constantCoefficient
Selecting heatTransferModel for (air in water): RanzMarshall
Selecting heatTransferModel for (water in air): RanzMarshall
Calculating field g.h

Reading field p_rgh

Creating field dpdt

Creating field kinetic energy K

No MRF models present

No finite volume options present

Courant Number mean: 0 max: 0
Max Ur Courant Number = 0
Calculating field DDtU1 and DDtU2


Starting time loop

surfaceFieldValue sFV_inletWater_phiAir:
    operation     = weightedSum
    total faces   = 500
    total area    = 0.05
    weight field  = alpha.air


surfaceFieldValue sFv_inletAir_phiAir:
    operation     = weightedSum
    total faces   = 200
    total area    = 0.02
    weight field  = alpha.air


surfaceFieldValue sFV_outlet_phiAir:
    operation     = weightedSum
    total faces   = 600
    total area    = 0.06
    weight field  = alpha.air


surfaceFieldValue sFV_inletWater_phiWater:
    operation     = weightedSum
    total faces   = 500
    total area    = 0.05
    weight field  = alpha.water


surfaceFieldValue sFV_inletAir_phiWater:
    operation     = weightedSum
    total faces   = 200
    total area    = 0.02
    weight field  = alpha.water


surfaceFieldValue sFV_outlet_phiWater:
    operation     = weightedSum
    total faces   = 600
    total area    = 0.06
    weight field  = alpha.water


Courant Number mean: 0 max: 0
Max Ur Courant Number = 0
deltaT = 0.00117647
Time = 0.00117647

PIMPLE: iteration 1
MULES: Solving for alpha.air
MULES: Solving for alpha.air
MULES: Solving for alpha.air
alpha.air volume fraction = 0.217166  Min(alpha.air) = 0  Max(alpha.air) = 1
MULES: Solving for alpha.air
MULES: Solving for alpha.air
MULES: Solving for alpha.air
alpha.air volume fraction = 0.217166  Min(alpha.air) = 0  Max(alpha.air) = 1
Constructing momentum equations
min T.air 293
min T.water 293


--> FOAM FATAL ERROR: 

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

2
(
turbulenceProperties.air
turbulenceProperties.water
)


    From function const Type& Foam::objectRegistry::lookupObject(const Foam::word&, bool) const [with Type = Foam::turbulenceModel]
    in file /home/buzz2/pawan/OpenFOAM/OpenFOAM-v1712/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 239.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::error::abort() at ??:?
#2  Foam::turbulenceModel const& Foam::objectRegistry::lookupObject<Foam::turbulenceModel>(Foam::word const&, bool) const at ??:?
#3  Foam::porousBafflePressureFvPatchField::updateCoeffs() at ??:?
#4  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::Boundary::updateCoeffs() at ??:?
#5  Foam::fvMatrix<double>::fvMatrix(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) at ??:?
#6  Foam::fv::EulerDdtScheme<double>::fvmDdt(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#7  Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::ddt<double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#8  ? at ??:?
#9  __libc_start_main in "/lib64/libc.so.6"
#10  ? at ??:?

[A H Gazi]

Hi..
i am facing kind of same problem.Did you solve your one?please let me know how??

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

Create mesh for time = 0


PIMPLE: no residual control data found. Calculations will employ 3 corrector loops


Reading g

Reading hRef
Creating twoPhaseSystem

Selecting thermodynamics package
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState rhoConst;
specie specie;
energy sensibleInternalEnergy;
}

Calculating face flux field phi.particles
Selecting diameterModel for phase particles: constant
Selecting turbulence model type RAS
Selecting RAS turbulence model kineticTheory
Selecting viscosityModel Gidaspow
Selecting conductivityModel Gidaspow
Selecting radialModel SinclairJackson
Selecting granularPressureModel Lun
Selecting frictionalStressModel JohnsonJackson
kineticTheoryCoeffs
{
equilibrium off;
e 0.8;
alphaMax 0.62;
alphaMinFriction 0.5;
residualAlpha 0.0001;
viscosityModel Gidaspow;
conductivityModel Gidaspow;
granularPressureModel Lun;
frictionalStressModel JohnsonJackson;
radialModel SinclairJackson;
JohnsonJacksonCoeffs
{
Fr 0.05;
eta 2;
p 5;
phi 28.5;
alphaDeltaMin 0.05;
}
}

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

Calculating face flux field phi.water
Selecting diameterModel for phase water: constant
Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
RAS
{
RASModel kEpsilon;
turbulence on;
printCoeffs on;
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 0;
sigmak 1;
sigmaEps 1.3;
}

Selecting default blending method: none
Selecting dragModel for (particles in water): GidaspowErgunWenYu
Selecting swarmCorrection for (particles in water): none
Selecting swarmCorrection for (particles in water): none
Selecting swarmCorrection for (particles in water): none
Selecting turbulentDispersionModel for (particles in water): constantCoefficient
Calculating field g.h

Reading field p_rgh

Creating field dpdt

Creating field kinetic energy K

No MRF models present

No finite volume options present

Courant Number mean: 0 max: 0
Max Ur Courant Number = 0
Calculating field DDtU1 and DDtU2


Starting time loop

fieldAverage fieldAverage1:
Starting averaging at time 0

Courant Number mean: 0 max: 0
Max Ur Courant Number = 0
Time = 0.002

PIMPLE: iteration 1
MULES: Solving for alpha.particles
MULES: Solving for alpha.particles
alpha.particles volume fraction = 1 Min(alpha.particles) = 1 Max(alpha.particles) = 1
Constructing momentum equations
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:perator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
#5 Foam::dragModels::Ergun::CdRe() const at ??:?
#6 Foam::dragModels::GidaspowErgunWenYu::CdRe() const at ??:?
#7 Foam::dragModel::Ki() const at ??:?
#8 Foam::dragModel::K() const at ??:?
#9 Foam::BlendedInterfacialModel<Foam::dragModel>::K( ) const at ??:?
#10 Foam::twoPhaseSystem::Kd() const at ??:?
#11 ? in "/opt/openfoam5/platforms/linux64GccDPInt32Opt/bin/twoPhaseEulerFoam"
#12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#13 ? in "/opt/openfoam5/platforms/linux64GccDPInt32Opt/bin/twoPhaseEulerFoam"
Floating point exception (core dumped)


thanks and regards.

[mbuergle]

Hi,
I did not manage to solve to problem with the porousBafflePressure BC but for me it worked to use the fanPressure BC instead. I don't know if this is of any help.
Best regards

[A H Gazi]

Originally Posted by mbuergle

Hi,
I did not manage to solve to problem with the porousBafflePressure BC but for me it worked to use the fanPressure BC instead. I don't know if this is of any help.
Best regards

Thanks for your reply.It may be funny but I have changed PBIC to GAMG in p_rgh and it is running well. Still in vain what is going on behind this .