CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Problems with compressible twoPhaseEulerFoam (https://www.cfd-online.com/Forums/openfoam-solving/108562-problems-compressible-twophaseeulerfoam.html)

 Sören Sander October 26, 2012 08:36

Problems with compressible twoPhaseEulerFoam

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::printStack(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::operator/<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 04: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?

 All times are GMT -4. The time now is 19:36.