[harsha_kulkarni]

Hello friends,

I'm doing simulation of two opposed jet microextractor using OpenFOAM. from one jet water is comming and fro other toluene. I'm using multiphaseInterFOAM and k-epsilon for turbulance and I'm getting floating point error . I'm attaching the geometry.
when I switch off the turbulance It woks.

inlet velocity are given as 1m/s
value of k-0.0015
value of epsilon -0.0454569

================================================== =============

Courant Number mean: 67.1235 max: 2.55783e+09
Interface Courant Number mean: 3.33317 max: 3.18604e+08
Time = 5e-05

MULES: Solving for alpha.water
water volume fraction, min, max = 7.64012e-06 -0.587206 1
MULES: Solving for alpha.toluene
toluene volume fraction, min, max = 7.63143e-06 -0.00422416 1
MULES: Solving for alpha.air
air volume fraction, min, max = -1.66329e+10 -1.31092e+34 1.63662e+34
Phase-sum volume fraction, min, max = -1.66329e+10 -1.31092e+34 1.63662e+34
MULES: Solving for alpha.water
water volume fraction, min, max = 8.08753e-06 -50.491 61.9137
MULES: Solving for alpha.toluene
toluene volume fraction, min, max = 8.07943e-06 -0.0364036 1
MULES: Solving for alpha.air
air volume fraction, min, max = 8.37159e+16 -2.13494e+41 1.86964e+41
Phase-sum volume fraction, min, max = 8.37159e+16 -2.13494e+41 1.86964e+41
MULES: Solving for alpha.water
water volume fraction, min, max = 8.53494e-06 -10980.1 17978.1
MULES: Solving for alpha.toluene
toluene volume fraction, min, max = 8.52743e-06 -0.318669 1
MULES: Solving for alpha.air
air volume fraction, min, max = 0 -1.55348e+49 1.02287e+49
Phase-sum volume fraction, min, max = 0 -1.55348e+49 1.02287e+49
MULES: Solving for alpha.water
water volume fraction, min, max = 8.98242e-06 -2.66078e+07 2.61686e+07
MULES: Solving for alpha.toluene
toluene volume fraction, min, max = 8.97543e-06 -3.42201 2.06948
MULES: Solving for alpha.air
air volume fraction, min, max = 5.02697e+33 -6.86083e+57 4.84214e+57
Phase-sum volume fraction, min, max = 5.02697e+33 -6.86083e+57 4.84214e+57
GAMG: Solving for p_rgh, Initial residual = 6.89188e-33, Final residual = 6.89188e-33, No Iterations 0
time step continuity errors : sum local = 6.60636e+07, global = -0.0168154, cumulative = -0.0255993
GAMG: Solving for p_rgh, Initial residual = 7.30448e-33, Final residual = 7.30448e-33, No Iterations 0
time step continuity errors : sum local = 7.00186e+07, global = 0.0165987, cumulative = -0.00900067
smoothSolver: Solving for epsilon, Initial residual = 0.245263, Final residual = 1.77353e-12, No Iterations 1
bounding epsilon, min: -3.01345e+22 max: 1.8657e+32 average: 2.29564e+27
smoothSolver: Solving for k, Initial residual = 0.999858, Final residual = 7.32211e-06, No Iterations 1000
bounding k, min: -504.816 max: 6.00298e+22 average: 8.58603e+20
ExecutionTime = 649.02 s ClockTime = 674 s

Courant Number mean: 5.34644e+07 max: 3.22865e+15
Interface Courant Number mean: 1.66442e+07 max: 3.22865e+15
Time = 6e-05

MULES: Solving for alpha.water
water volume fraction, min, max = 9106.17 -1.36203e+16 1.325e+16
MULES: Solving for alpha.toluene
toluene volume fraction, min, max = 9.42388e-06 -59600 167014
MULES: Solving for alpha.air
air volume fraction, min, max = -4.31814e+42 -1.21498e+67 1.30926e+67
Phase-sum volume fraction, min, max = -4.31814e+42 -1.21498e+67 1.30926e+67
MULES: Solving for alpha.water
water volume fraction, min, max = 5.0864e+08 -5.20075e+27 3.97036e+27
MULES: Solving for alpha.toluene
toluene volume fraction, min, max = 9.87405e-06 -2.04296e+12 2.24217e+12
MULES: Solving for alpha.air
air volume fraction, min, max = 4.63656e+51 -2.58722e+76 2.33644e+76
Phase-sum volume fraction, min, max = 4.63656e+51 -2.58722e+76 2.33644e+76
MULES: Solving for alpha.water
water volume fraction, min, max = 5.88807e+14 -8.99561e+35 6.87842e+35
MULES: Solving for alpha.toluene
toluene volume fraction, min, max = -0.00736009 -2.89893e+22 3.47649e+22
MULES: Solving for alpha.air
air volume fraction, min, max = -1.99139e+61 -4.72101e+85 5.41195e+85
Phase-sum volume fraction, min, max = -1.99139e+61 -4.72101e+85 5.41195e+85
MULES: Solving for alpha.water
water volume fraction, min, max = -1.80447e+26 -5.22978e+49 5.69986e+49
MULES: Solving for alpha.toluene
toluene volume fraction, min, max = 1.12489e+21 -2.58103e+45 2.60982e+45
MULES: Solving for alpha.air
air volume fraction, min, max = 0 -1.23016e+95 1.68637e+95
Phase-sum volume fraction, min, max = 0 -1.23016e+95 1.68637e+95
GAMG: Solving for p_rgh, Initial residual = 0.999984, Final residual = 8782.28, No Iterations 1000
time step continuity errors : sum local = 3.59764e+45, global = 8.32869e+31, cumulative = 8.32869e+31
GAMG: Solving for p_rgh, Initial residual = 1.77362e-66, Final residual = 1.77362e-66, No Iterations 0
time step continuity errors : sum local = 5.13781e+57, global = -2.39079e+41, cumulative = -2.39079e+41
[1] #0 [0] #0 [3] #0 [2] #0 Foam::error:rintStack(Foam::Ostream&)Foam::error :rintStack(Foam::Ostream&)Foam::error:rintStac k(Foam::Ostream&)Foam::error:rintStack(Foam::Ost ream&) at ??:?
[0] #1 Foam::sigFpe::sigHandler(int) at ??:?
[1] #1 Foam::sigFpe::sigHandler(int) at ??:?
[3] #1 Foam::sigFpe::sigHandler(int) at ??:?
[2] #1 Foam::sigFpe::sigHandler(int) at ??:?
[2] #2 Uninterpreted:
[2] #3 Foam::symGaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:?
[3] #2 Uninterpreted:
[3] #3 Foam::symGaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:?
[0] #2 Uninterpreted:
[0] #3 Foam::symGaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:?
[1] #2 Uninterpreted:
[1] #3 Foam::symGaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:?
[2] #4 Foam::symGaussSeidelSmoother::smooth(Foam::Field<d ouble>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
[0] #4 Foam::symGaussSeidelSmoother::smooth(Foam::Field<d ouble>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
[1] #4 Foam::symGaussSeidelSmoother::smooth(Foam::Field<d ouble>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
[3] #4 Foam::symGaussSeidelSmoother::smooth(Foam::Field<d ouble>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
[0] #5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[3] #5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[2] #5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[1] #5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[2] #6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
[3] #6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
[0] #6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
[1] #6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
[3] #7 at ??:?
[1] #7 at ??:?
[0] #7 at ??:?
[2] #7 Foam::fvMatrix<double>::solve(Foam::dictionary const&)Foam::fvMatrix<double>::solve(Foam::diction ary const&)Foam::fvMatrix<double>::solve(Foam::diction ary const&)Foam::fvMatrix<double>::solve(Foam::diction ary const&) at ??:?
[0] #8 Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<doubl e> > const&) at ??:?
[1] #8 Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<doubl e> > const&) at ??:?
[2] #8 Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<doubl e> > const&) at ??:?
[3] #8 Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<doubl e> > const&) at ??:?
[1] #9 Foam::incompressible::RASModels::kEpsilon::correct () at ??:?
[0] #9 Foam::incompressible::RASModels::kEpsilon::correct () at ??:?
[2] #9 Foam::incompressible::RASModels::kEpsilon::correct () at ??:?
[0] #10 at ??:?
[1] #10 at ??:?
[3] #9 Foam::incompressible::RASModels::kEpsilon::correct () at ??:?
[2] #10
at ??:?
[3] #10
================================================== ===============
can you tell me how to tackle this error

for analysis purpose I also want velocity data which is also not obtained in log file

[rcastilla]

Check where is k diverging. It can be a mesh problem...

[harsha_kulkarni]

but for singlePhase same model working for given k. when I'm including multiphase it is producing error.

[nero235]

Late to the party but for multiphase flow you encounter a problem that you can't specify an constant y+ for your case (if U would be constant). The dimensionless wall distance y+ depends on the kinematic viscosity nu, which for multiphase flow is different in each phase. Say you have 3 phases: air, water and oil. You will have 3 different nu and therefore is y+ different for each phase, even if you have the same velocity in each phase.

This problem can be solved by using a limited scheme (instead of corrected) for the

laplacianScheme, sndGradScheme

eg.:

Gauss linear limited 0.5

Regards, Sebastian