CFD Online Discussion Forums - solid-liquid two phase flow for two phase Euler Foam

Hello!!

I'm a beginner to OpenFOAM and i'm working on twoPhaseEulerFoam and simulate solid-liquid two phase fluid. That is liquid (water or oil) with carrying solid particle (Diameter:0.0003m, the volume fraction:0.3) flows in a pipeline. It is worthing noting that i have deleted the lift force, virtual mass force and wall lubrication force from the solver code. In addition, i have deleted the temperature equation from the solver code in order to match the facts.

However the terminal always hinted the following error when i tried to run the case.

The error i met is the following:
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 4.1
Exec : twoPhaseEulerFoamHydrate
Date : Aug 30 2017
Time : 15:07:12
Host : "JK"
PID : 4190
Case : /home/jk/OpenFOAM/jk-4.1/run/paperCase/solid-liquid two phase flows
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
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.air
Selecting diameterModel for phase air: constant
Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 -0.33;
sigmak 1;
sigmaEps 1.3;
}

Selecting default blending method: none
Selecting dragModel for (particles in air): GidaspowErgunWenYu
Selecting swarmCorrection for (particles in air): none
Selecting swarmCorrection for (particles in air): none
Selecting swarmCorrection for (particles in air): none
Selecting heatTransferModel for (particles 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.00180506 max: 2.27071
Max Ur Courant Number = 0
Calculating field DDtU1 and DDtU2

Starting time loop

fieldAverage fieldAverage1:
Starting averaging at time 0

Courant Number mean: 0.00180506 max: 2.27071
Max Ur Courant Number = 0
Time = 0.001

PIMPLE: iteration 1
MULES: Solving for alpha.particles
MULES: Solving for alpha.particles
alpha.particles volume fraction = 0.3 Min(alpha.particles) = 0.3 Max(alpha.particles) = 0.3
Constructing momentum equations
GAMG: Solving for p_rgh, Initial residual = 0.991305, Final residual = 8.79876e-09, No Iterations 28
PIMPLE: iteration 2
MULES: Solving for alpha.particles
MULES: Solving for alpha.particles
smoothSolver: Solving for alpha.particles, Initial residual = 0.00150465, Final residual = 2.74931e-13, No Iterations 2
alpha.particles volume fraction = 0.315835 Min(alpha.particles) = -6.15607 Max(alpha.particles) = 1.72539
Constructing momentum equations
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? addr2line failed
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) addr2line failed
#4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<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&) addr2line failed
#5 Foam::dragModels::Ergun::CdRe() const addr2line failed
#6 Foam::dragModels::GidaspowErgunWenYu::CdRe() const addr2line failed
#7 Foam::dragModel::Ki() const addr2line failed
#8 Foam::dragModel::K() const addr2line failed
#9 Foam::BlendedInterfacialModel<Foam::dragModel>::K( ) const addr2line failed
#10 Foam::twoPhaseSystem::Kd() const addr2line failed
#11 ?
#12 __libc_start_main addr2line failed
#13 ?

************************************************** *******************************************
#3 and #5 indicates the denominator is 0 during the calculation of drag force. I have tried to modify the case, but still haven't found a solution.

Any help is greatly appreciated.
Best Regards
Kai Jiang