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

solid-liquid two phase flow for two phase Euler Foam

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 30, 2017, 04:24
Question solid-liquid two phase flow for two phase Euler Foam
  #1
New Member
 
Kai Jiang
Join Date: Aug 2017
Posts: 19
Rep Power: 9
Jack_Kerouac is on a distinguished road
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: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 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
Jack_Kerouac is offline   Reply With Quote

Old   November 1, 2017, 01:45
Default error while running two phase euler foam
  #2
New Member
 
DEEPTI
Join Date: Sep 2017
Posts: 2
Rep Power: 0
deeptiof is on a distinguished road
while running the twophaseeulerfoam solver after some iterations of my solution i got these errors. please help me to resolve these issues.

#0 Foam::error::PrintStack(Foam::Ostream&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2 Uninterpreted:
#3 Foam:;DILUPreconditioner::calcReciprocalD(Foam::Fi eld<double>&, Foam::lduMatrix const&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#4 Foam:;DILUPreconditioner:;DILUPreconditioner(Foam: :lduMatrix::solver const&, Foam::dictionary const&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#5 Foam::lduMatrix::Preconditioner::addasymMatrixCons tructorToTable<Foam:;DILUPreconditioner>::New(Foam ::lduMatrix::solver const&, Foam::dictionary const&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#6 Foam::lduMatrix::Preconditioner::New(Foam::lduMatr ix::solver const&, Foam::dictionary const&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#7 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#8 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
#9 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam222/platforms/linuxGccDPOpt/bin/twoPhaseEulerFoam"
#10 Foam::fvMatrix<double>::solve() in "/opt/openfoam222/platforms/linuxGccDPOpt/bin/twoPhaseEulerFoam"
#11
in "/opt/openfoam222/platforms/linuxGccDPOpt/bin/twoPhaseEulerFoam"
#12 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#13
in "/opt/openfoam222/platforms/linuxGccDPOpt/bin/twoPhaseEulerFoam"
Floating point exception (core dumped)


Thank you
deeptiof is offline   Reply With Quote

Reply

Tags
drag forces, solid-liquid flows, twophaseeulerfoam


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
problem during mpi in server: expected Scalar, found on line 0 the word 'nan' muth OpenFOAM Running, Solving & CFD 3 August 27, 2018 05:18
Error - Solar absorber - Solar Thermal Radiation MichaelK CFX 12 September 1, 2016 06:15
[Commercial meshers] Using starToFoam clo OpenFOAM Meshing & Mesh Conversion 33 September 26, 2012 05:04
2 phase flow liquid flow solid particles sachin OpenFOAM Running, Solving & CFD 0 August 28, 2008 17:51
fluid flow fundas ram Main CFD Forum 5 June 17, 2000 22:31


All times are GMT -4. The time now is 03:39.