gustavoliveira |
June 30, 2021 14:05 |
DamBreak case multiphaseEulerFoam with RAS
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:
Code:
...
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!
|