[jason]

Hello,


I was wondering if anyone has successfully specified a RAS model for all phases of a 3-phase flow?


I have oil/air/water in a pipe,



1) when I set all to kEpsilon I have an immediate crash at the first iteration



GAMG: Solving for p_rgh, Initial residual = 0.00296003, Final residual = 1.99966e-05, No Iterations 12
GAMG: Solving for p_rgh, Initial residual = 0.00122383, Final residual = 1.07626e-05, No Iterations 13
GAMGPCG: Solving for p_rgh, Initial residual = 1.92961e-05, Final residual = 8.79505e-08, No Iterations 10
#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<d ouble>&, 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::dict ionary const&) at ??:?
#7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam"
#8 Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<doubl e> > const&) at ??:?
#9 Foam::RASModels::kEpsilon<Foam::PhaseCompressibleM omentumTransportModel<Foam::phaseModel> >::correct() at ??:?
#10 Foam::phaseSystem::correctTurbulence() at ??:?
#11 ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam"
#12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#13 ? in "/opt/openfoam8/platforms/linux64GccDPInt32Opt/bin/multiphaseEulerFoam"
Floating point exception (core dumped)







2) when I set all to mixtureKEpsilon I have the following error



--> FOAM FATAL ERROR:
Call from a two-phase model in a multi-phase system.

From function const Foam::phaseModel& Foam::phaseSystem::otherPhase(const Foam::phaseModel&) const
in file ../phaseSystems/lnInclude/phaseSystemI.H at line 118.






3) When I set all to LaheyKEpsilon I have the following error


--> FOAM FATAL ERROR:
Call from a two-phase model in a multi-phase system.

From function const Foam::phaseModel& Foam::phaseSystem::otherPhase(const Foam::phaseModel&) const
in file ../phaseSystems/lnInclude/phaseSystemI.H at line 118.







Could advise on the error in (1)...


And have (2) and (3) not been implemented for 3-phase and more?




Jason

[Elia38]

Hi jason. Could you solve the problem?