[gustavoliveira]

Hello,


does anyone have run the dambreak case in the multiphaseeulerfoam solver but with a turbulence model? When I use a tubulence model (kepsilon or komegasst) occur this error:

...

GAMGPCG:  Solving for p_rgh, Initial residual = 4.06039e-05, Final residual = 2.92177e-08, No Iterations 4
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#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 ??:?
#4  Foam::symGaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
#5  Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam"
#8  Foam::fvMatrix<double>::solve() in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam"
#9  Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<double> > const&) at ??:?
#10  Foam::kOmegaSST<Foam::eddyViscosity<Foam::RASModel<Foam::PhaseCompressibleMomentumTransportModel<Foam::phaseModel> > >, Foam::PhaseCompressibleMomentumTransportModel<Foam::phaseModel> >::correct() at ??:?
#11  Foam::phaseSystem::correctTurbulence() at ??:?
#12  ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam"
#13  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#14  ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam"

It seems something in the BC of alphat.air, but I dont have sure and already have search a lot without found a solution.


Thanks!

[user_cyy]

Originally Posted by gustavoliveira

Hello,


does anyone have run the dambreak case in the multiphaseeulerfoam solver but with a turbulence model? When I use a tubulence model (kepsilon or komegasst) occur this error:

...

GAMGPCG:  Solving for p_rgh, Initial residual = 4.06039e-05, Final residual = 2.92177e-08, No Iterations 4
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#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 ??:?
#4  Foam::symGaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
#5  Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam"
#8  Foam::fvMatrix<double>::solve() in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam"
#9  Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<double> > const&) at ??:?
#10  Foam::kOmegaSST<Foam::eddyViscosity<Foam::RASModel<Foam::PhaseCompressibleMomentumTransportModel<Foam::phaseModel> > >, Foam::PhaseCompressibleMomentumTransportModel<Foam::phaseModel> >::correct() at ??:?
#11  Foam::phaseSystem::correctTurbulence() at ??:?
#12  ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam"
#13  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#14  ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam"

It seems something in the BC of alphat.air, but I dont have sure and already have search a lot without found a solution.


Thanks!

Have you solved your case ?

[gustavoliveira]

Originally Posted by user_cyy

Have you solved your case ?

Actually not, but using mixtureKEpsilon in both phases works. But I discover multiphaseEulerFoam is not apropriate to this kind type of problem.