CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   DamBreak case multiphaseEulerFoam with RAS (https://www.cfd-online.com/Forums/openfoam-solving/237098-dambreak-case-multiphaseeulerfoam-ras.html)

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!

user_cyy December 3, 2021 04:28

Have you solved your case ?
 
Quote:

Originally Posted by gustavoliveira (Post 807211)
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!

Have you solved your case ?

gustavoliveira December 3, 2021 19:14

Quote:

Originally Posted by user_cyy (Post 817815)
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.


All times are GMT -4. The time now is 10:57.