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

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?