Problems with compressible twoPhaseEulerFoam

Sören Sander · October 26, 2012, 08:36

Hi here,

I imported my own mesh to the testcase of the compressible tPEF, adopted the BC and started a case. Anyways it always throws the same error:

Courant Number mean: 6.63685e-09 max: 0.00185723
Max Ur Courant Number = 0.00185723
deltaT = 1.2e-05
Time = 1.2e-05

DILUPBiCG: Solving for alpha1, Initial residual = 1, Final residual = 4.67011e-17, No Iterations 1
Dispersed phase volume fraction = 1 Min(alpha1) = 0 Max(alpha1) = 1
DILUPBiCG: Solving for alpha1, Initial residual = 0.0621701, Final residual = 2.18343e-18, No Iterations 1
Dispersed phase volume fraction = 0.999999 Min(alpha1) = 0 Max(alpha1) = 1
DILUPBiCG: Solving for epsilon, Initial residual = 0.0171299, Final residual = 9.50906e-06, No Iterations 87
DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 8.80289e-06, No Iterations 166
DILUPBiCG: Solving for T1, Initial residual = 1, Final residual = 2.2903e-14, No Iterations 2
DILUPBiCG: Solving for T2, Initial residual = 1, Final residual = 2.08691e-11, No Iterations 1
#0 Foam::error:

rintStack(Foam::Ostream&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/libc.so.6"
#3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:

perator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/compressibleTwoPhaseEulerFoam"
#5
in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/compressibleTwoPhaseEulerFoam"
#6 __libc_start_main in "/lib/libc.so.6"
#7
in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/compressibleTwoPhaseEulerFoam"
Floating point exception

Sören Sander · November 1, 2012, 03:58

I found that the problem was with the inlet boundary condition of U2. I tried to use a fixedValue of uniform (-1 0 0) instead of $internalField. If I change the BC to $interalField again and the solver runs without problems. However I want to specify an inlet velocity for phase 2. How can I do this?

October 26, 2012, 08:36	Problems with compressible twoPhaseEulerFoam	#1
Sören Sander Member Sören Sander Join Date: Sep 2012 Location: Germany Posts: 43 Rep Power: 13	Hi here, I imported my own mesh to the testcase of the compressible tPEF, adopted the BC and started a case. Anyways it always throws the same error: Courant Number mean: 6.63685e-09 max: 0.00185723 Max Ur Courant Number = 0.00185723 deltaT = 1.2e-05 Time = 1.2e-05 DILUPBiCG: Solving for alpha1, Initial residual = 1, Final residual = 4.67011e-17, No Iterations 1 Dispersed phase volume fraction = 1 Min(alpha1) = 0 Max(alpha1) = 1 DILUPBiCG: Solving for alpha1, Initial residual = 0.0621701, Final residual = 2.18343e-18, No Iterations 1 Dispersed phase volume fraction = 0.999999 Min(alpha1) = 0 Max(alpha1) = 1 DILUPBiCG: Solving for epsilon, Initial residual = 0.0171299, Final residual = 9.50906e-06, No Iterations 87 DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 8.80289e-06, No Iterations 166 DILUPBiCG: Solving for T1, Initial residual = 1, Final residual = 2.2903e-14, No Iterations 2 DILUPBiCG: Solving for T2, Initial residual = 1, Final residual = 2.08691e-11, No Iterations 1 #0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #2 in "/lib/libc.so.6" #3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" #4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:perator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/compressibleTwoPhaseEulerFoam" #5 in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/compressibleTwoPhaseEulerFoam" #6 __libc_start_main in "/lib/libc.so.6" #7 in "/opt/openfoam210/platforms/linux64GccDPOpt/bin/compressibleTwoPhaseEulerFoam" Floating point exception

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Blow of compressible solver while using K-epsilon model in openfoam	Amit Mathur	OpenFOAM	16	October 6, 2013 11:09
Compressible code problems	praveen	OpenFOAM Programming & Development	0	April 11, 2011 00:00
Fixed grid methods for compressible viscous flow	liujmljm	Main CFD Forum	1	November 7, 2010 17:54
compressible flow	maria teresa	FLUENT	1	September 7, 2007 16:58
urgent help needed with 2d compressible flow	James	FLUENT	2	June 20, 2007 04:22

November 1, 2012, 03:58		#2
Sören Sander Member Sören Sander Join Date: Sep 2012 Location: Germany Posts: 43 Rep Power: 13	I found that the problem was with the inlet boundary condition of U2. I tried to use a fixedValue of uniform (-1 0 0) instead of $internalField. If I change the BC to $interalField again and the solver runs without problems. However I want to specify an inlet velocity for phase 2. How can I do this?