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

Error using porousBafflePressure BC with twoPhaseEulerFoam

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 31, 2018, 09:47
Default Error using porousBafflePressure BC with twoPhaseEulerFoam
  #1
New Member
 
Matthias Bürgler
Join Date: Mar 2018
Posts: 3
Rep Power: 5
mbuergle is on a distinguished road
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

Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  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 ??:?
Attached Files
File Type: zip TestCase1.zip (25.2 KB, 8 views)
mbuergle is offline   Reply With Quote

Old   August 9, 2018, 02:18
Default
  #2
New Member
 
Ainal Hoque Gazi
Join Date: May 2018
Location: India
Posts: 27
Rep Power: 5
A H Gazi is on a distinguished road
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.
A H Gazi is offline   Reply With Quote

Old   August 14, 2018, 06:09
Default
  #3
New Member
 
Matthias Bürgler
Join Date: Mar 2018
Posts: 3
Rep Power: 5
mbuergle is on a distinguished road
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
mbuergle is offline   Reply With Quote

Old   August 14, 2018, 09:25
Default
  #4
New Member
 
Ainal Hoque Gazi
Join Date: May 2018
Location: India
Posts: 27
Rep Power: 5
A H Gazi is on a distinguished road
Quote:
Originally Posted by mbuergle View Post
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 .
A H Gazi is offline   Reply With Quote

Reply

Tags
porousbafflepressure, twophaseeulerfoam

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
twoPhaseEulerFoam: sudden crash hcen OpenFOAM Running, Solving & CFD 18 September 29, 2020 03:04
Oscillating a porousBafflePressure (cyclic) boundary mhasif02 OpenFOAM 0 April 10, 2018 06:45
twoPhaseEulerFoam fluidized bed wrong results FabianB OpenFOAM Running, Solving & CFD 0 June 14, 2017 03:48
Is twoPhaseEulerFoam applicable to 3D cases / delivering erroneous results? ThomasV OpenFOAM 0 November 11, 2013 08:10
twoPhaseEulerFoam freemankofi OpenFOAM 0 May 23, 2011 16:24


All times are GMT -4. The time now is 19:56.