3d multiPhaseEulerFoam error

Hamzeh_Mirab · January 12, 2015, 09:16

I wanna run my case but an error appear in first time step:

Courant Number mean: 0.00205988 max: 0.00490449
Time = 0.01

MULES: Solving for alpha.oil
oil volume fraction, min, max = 0 0 1
MULES: Solving for alpha.water
water volume fraction, min, max = 1 0 1
Phase-sum volume fraction, min, max = 1 0 1.5
MULES: Solving for alpha.oil
oil volume fraction, min, max = 0 0 1
MULES: Solving for alpha.water
water volume fraction, min, max = 1 0 1
Phase-sum volume fraction, min, max = 1 0 1.5
#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::multiphaseSystem::nu() const at ??:?
#6 Foam::incompressible::LESModels::laminar::nuSgs() const at ??:?
#7 Foam::incompressible::LESModel::nut() const at ??:?
#8
at ??:?
#9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10
at ??:?
Floating point exception (core dumped)

whats the problem?

nimasam · January 12, 2015, 10:36

turn off turbulence and run your case again

, it seems it relates to your turbulence setup

Hamzeh_Mirab · January 13, 2015, 11:34

thank you,
my lesmodel is laminar. how can I turn off turbulence?

kwardle · January 20, 2015, 15:21

Looks to me like something might be wrong with your initialization of the alpha fields. Why does the phase sum show 1.5? Are you sure everything is consistent? It is not exactly failing due to turbulence, but rather when calculating the mixture viscosity nu() for the turbulenceModel--which it will use in either case for laminar or turbulent. Check your initialization.

Hamzeh_Mirab · January 21, 2015, 02:31

Thank you for your response ,
It fixed by setFields and as you truly mentioned void fraction was not consistent.However, my new problem is that my case is 3d and I wanna view 3d output but dispersed phase dont enter in inner domain and only dispersed in boundaries.You can view my case in link below :
http://www.cfd-online.com/Forums/ope...eulerfoam.html

January 12, 2015, 09:16	3d multiPhaseEulerFoam error	#1
Hamzeh_Mirab New Member Hamzeh Join Date: Oct 2014 Posts: 11 Rep Power: 11	I wanna run my case but an error appear in first time step: Courant Number mean: 0.00205988 max: 0.00490449 Time = 0.01 MULES: Solving for alpha.oil oil volume fraction, min, max = 0 0 1 MULES: Solving for alpha.water water volume fraction, min, max = 1 0 1 Phase-sum volume fraction, min, max = 1 0 1.5 MULES: Solving for alpha.oil oil volume fraction, min, max = 0 0 1 MULES: Solving for alpha.water water volume fraction, min, max = 1 0 1 Phase-sum volume fraction, min, max = 1 0 1.5 #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::multiphaseSystem::nu() const at ??:? #6 Foam::incompressible::LESModels::laminar::nuSgs() const at ??:? #7 Foam::incompressible::LESModel::nut() const at ??:? #8 at ??:? #9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #10 at ??:? Floating point exception (core dumped) whats the problem?

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
[OpenFOAM.org] compile error in dynamicMesh and thermophysicalModels libraries	NickG	OpenFOAM Installation	3	December 30, 2019 00:21
[blockMesh] blockMesh with double grading.	spwater	OpenFOAM Meshing & Mesh Conversion	92	January 12, 2019 09:00
[swak4Foam] GroovyBC the dynamic cousin of funkySetFields that lives on the suburb of the mesh	gschaider	OpenFOAM Community Contributions	300	October 29, 2014 18:00
OpenFOAM without MPI	kokizzu	OpenFOAM Installation	4	May 26, 2014 09:17
DecomposePar links against liblamso0 with OpenMPI	jens_klostermann	OpenFOAM Bugs	11	June 28, 2007 17:51

January 12, 2015, 10:36		#2
nimasam Senior Member Nima Samkhaniani Join Date: Sep 2009 Location: Tehran, Iran Posts: 1,266 Blog Entries: 1 Rep Power: 24	turn off turbulence and run your case again , it seems it relates to your turbulence setup __________________ My Personal Website (http://nimasamkhaniani.ir/) Telegram channel (https://t.me/cfd_foam)

January 13, 2015, 11:34		#3
Hamzeh_Mirab New Member Hamzeh Join Date: Oct 2014 Posts: 11 Rep Power: 11	thank you, my lesmodel is laminar. how can I turn off turbulence?

January 20, 2015, 15:21		#4
kwardle Senior Member Kent Wardle Join Date: Mar 2009 Location: Illinois, USA Posts: 219 Rep Power: 21	Looks to me like something might be wrong with your initialization of the alpha fields. Why does the phase sum show 1.5? Are you sure everything is consistent? It is not exactly failing due to turbulence, but rather when calculating the mixture viscosity nu() for the turbulenceModel--which it will use in either case for laminar or turbulent. Check your initialization.

January 21, 2015, 02:31		#5
Hamzeh_Mirab New Member Hamzeh Join Date: Oct 2014 Posts: 11 Rep Power: 11	Thank you for your response , It fixed by setFields and as you truly mentioned void fraction was not consistent.However, my new problem is that my case is 3d and I wanna view 3d output but dispersed phase dont enter in inner domain and only dispersed in boundaries.You can view my case in link below : http://www.cfd-online.com/Forums/ope...eulerfoam.html