[remidemol]

Hi everyone,

I've run several times my case (a bubble column) using twoPhaseEulerFoam and mixtureKEpsilon turbulent model for both phases.
I want to change the turbulence model to the standard kEpsilon model. So I changed the turbulentProperties files. But it crashed at the first iteration.

Maybe the answer is in the forum but I didn't find it yet.

Does anyone have a clue to fix it ?

Edit : I have to had that I consider the free surface in my simulation, that's maybe the raison why it crashed.

Thanks in advance for your answer.

My version of OF is 3.0.1
Here is what I got in my terminal

Starting time loop

Courant Number mean: 0.000104375 max: 0.00166667
Max Ur Courant Number = 0.00166667
deltaT = 0.00119617
Time = 0.00119617

PIMPLE: iteration 1
MULES: Solving for alpha.air
MULES: Solving for alpha.air
alpha.air volume fraction = 0.0624943  Min(alpha.air) = 0  Max(alpha.air) = 1
Constructing momentum equations
smoothSolver:  Solving for e.air, Initial residual = 1, Final residual = 4.30053e-09, No Iterations 2
smoothSolver:  Solving for e.water, Initial residual = 0.998714, Final residual = 1.41482e-07, No Iterations 1000
min T.air 349.868
min T.water 350
GAMG:  Solving for p_rgh, Initial residual = 1, Final residual = 3.19227e-09, No Iterations 9
PIMPLE: iteration 2
MULES: Solving for alpha.air
MULES: Solving for alpha.air
alpha.air volume fraction = 0.0625092  Min(alpha.air) = 0  Max(alpha.air) = 1
Constructing momentum equations
smoothSolver:  Solving for e.air, Initial residual = 0.951929, Final residual = 8.62357e-10, No Iterations 2
smoothSolver:  Solving for e.water, Initial residual = 0.00241502, Final residual = 3.63505e-10, No Iterations 1
min T.air 350
min T.water 350
GAMG:  Solving for p_rgh, Initial residual = 0.000249446, Final residual = 3.19735e-09, No Iterations 5
PIMPLE: iteration 3
MULES: Solving for alpha.air
MULES: Solving for alpha.air
alpha.air volume fraction = 0.0625042  Min(alpha.air) = 0  Max(alpha.air) = 1
Constructing momentum equations
smoothSolver:  Solving for e.air, Initial residual = 0.998442, Final residual = 3.91661e-11, No Iterations 2
smoothSolver:  Solving for e.water, Initial residual = 0.00105975, Final residual = 2.61401e-10, No Iterations 1
min T.air 348.551
min T.water 350
GAMG:  Solving for p_rgh, Initial residual = 0.000489843, Final residual = 8.84302e-10, No Iterations 5
#0  "Foam::error::printStack(Foam::Ostream&)" in /Volumes/OpenFOAM/OpenFOAM-3.0.1/platforms/darwin64ClangDPInt32Opt/lib/libOpenFOAM.dylib
#1  "Foam::sigFpe::sigHandler(int)" in /Volumes/OpenFOAM/OpenFOAM-3.0.1/platforms/darwin64ClangDPInt32Opt/lib/libOpenFOAM.dylib
#2  "_sigtramp" in /usr/lib/system/libsystem_platform.dylib
#3  ? in /usr/lib/system/libsystem_platform.dylib
#4  "Foam::symGaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const" in /Volumes/OpenFOAM/OpenFOAM-3.0.1/platforms/darwin64ClangDPInt32Opt/lib/libOpenFOAM.dylib
#5  "Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const" in /Volumes/OpenFOAM/OpenFOAM-3.0.1/platforms/darwin64ClangDPInt32Opt/lib/libOpenFOAM.dylib
#6  "Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&)" in /Volumes/OpenFOAM/OpenFOAM-3.0.1/platforms/darwin64ClangDPInt32Opt/lib/libfiniteVolume.dylib
#7  "Foam::fvMatrix<double>::solve(Foam::dictionary const&)" in /Volumes/OpenFOAM/OpenFOAM-3.0.1/platforms/darwin64ClangDPInt32Opt/lib/libphaseCompressibleTurbulenceModels.dylib
#8  "Foam::fvMatrix<double>::solve()" in /Volumes/OpenFOAM/OpenFOAM-3.0.1/platforms/darwin64ClangDPInt32Opt/lib/libphaseCompressibleTurbulenceModels.dylib
#9  "Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<double> > const&)" in /Volumes/OpenFOAM/OpenFOAM-3.0.1/platforms/darwin64ClangDPInt32Opt/lib/libphaseCompressibleTurbulenceModels.dylib
#10  "Foam::RASModels::kEpsilon<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::PhaseCompressibleTurbulenceModel<Foam::phaseModel> > > >::correct()" in /Volumes/OpenFOAM/OpenFOAM-3.0.1/platforms/darwin64ClangDPInt32Opt/lib/libphaseCompressibleTurbulenceModels.dylib
#11  "Foam::twoPhaseSystem::correctTurbulence()" in /Volumes/OpenFOAM/OpenFOAM-3.0.1/platforms/darwin64ClangDPInt32Opt/lib/libcompressibleTwoPhaseSystem.dylib
#12  "main" in /Volumes/OpenFOAM/OpenFOAM-3.0.1/platforms/darwin64ClangDPInt32Opt/bin/twoPhaseEulerFoam
#13  "start" in /usr/lib/system/libdyld.dylib
#14  ? in /usr/lib/system/libdyld.dylib
Floating point exception: 8

[maalan]

Hi Rémi!

Did you find how to solve this problem?

Best,
Antonio

[remidemol]

Hi Maalan,

I didn't solve that problem with the twoPhaseEulerFoam solver but it works with multiphaseEulerFoam.

I suggest you to move on this solver in which you can define 2 phases or more than 2.

Best regards.

[maalan]

Great! I'm using OpenFOAM 2.3.x and I had to modify multiphaseEulerFoam to have RAS models also included.

Thanks Rémi!

Regards

[farzadpolytechnic]

hi remi and malaan

i have same problem with twoPhaseEulerFoam. so do you solve your problem by changing solver?
i really have limitation in my time
so do you think multiPhaseEulerFoam solve problem about modeling air bubble plume and using different turbulence model or i have to do something else?

thanks a lot
farzad

[atul_walia]

Hi, how do i define the proper values of drag and lift force in bubble column, is there any formula for that or i just have to choose a random value, please help me out,

[remidemol]

Hi Atul,

You can define drag, lift, turbulent dispersion, added-mass or wall lubrication models from file constant/phaseProperties

Few models have been already implemented, you can put a question mark after keyword "type" and then run twoPhaseEulerFoam (or other solver) to see the different models. You can also implement an other.