[Jack_Kerouac]

Hello!
I'm a beginner to OpenFOAM and i'm working on twoPhaseEulerFoam and simulate solid-liquid two phase flow. That is the liquid (water or oil) with carring the solid particles (diameter:0.0003m, the volume fraction:0.3 ) flows in a pipeline. It is wothing noting that i have deleted the lift force, virtual mass force, wall lubrication and the energy 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 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:rintStack(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: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&) 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 main source of the error comes from the denominator is zero during the calculation of drag forces. I have tried to modify the case, but still haven't found the solution.

Thank you very much.
Kai Jiang

[tomf]

Hi, from your log:

alpha.particles volume fraction = 0.315835 Min(alpha.particles) = -6.15607 Max(alpha.particles) = 1.72539

This indicates that the volume fraction of particles becomes < 0 and > 1 in some areas of the mesh, which is probably why the next equation fails. This may be due to bad mesh, bad initial conditions or unstable numerical schemes.

Regards,
Tom

[Jack_Kerouac]

Thanks so much for response. Does the wrong initial condition of phase fraction cause the problem?

[tomf]

It could be, it could also be the time step since you have a Courant number larger dan 1. Maybe start with a smaller time step?

I cannot help if you do not share the details of the case (like your 0 folder, fvSchemes/fvSolution/controlDict, checkMesh result)

Regards,
Tom

[Jack_Kerouac]

Thank you again. I'm uploading my case.

Regards,
Kai Jiang

[Jack_Kerouac]

Hi,

the case always can't upload successfully. I will send an email to you.

Regards
Kai Jiang