[remidemol]

Hello everyone,


I try to model a fluidized bed with reactingMultiphaseEulerFoam (3 phases) and I can run this simulation properly without turbulence model.


Unfortunately when I add a RAS model (kEpsilon) the simulation crashes.
It seems working well during the first time step but the simulations stops before begining the next time step. And no error il the log file



I have read this problem could reveal a missing property but everything seems ok. I try to modify the BC, the discretization schemes...



Does anyone have a idea on what is going on ? I put the log file below.



Thank you.

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


PIMPLE: No convergence criteria found


PIMPLE: No corrector convergence criteria found
        Calclations will do 3 corrections


Reading g

Reading hRef
Creating phaseSystem

Selecting multiphaseSystem basicMultiphaseSystem
Selecting phaseModel for particles: purePhaseModel
Selecting diameterModel for phase particles: constant
Selecting thermodynamics package 
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState rhoConst;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Calculating face flux field phi.particles
Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
bounding k.particles, min: 0 max: 0 average: 0
bounding epsilon.particles, min: 0 max: 0 average: 0
RAS
{
    RASModel        kEpsilon;
    turbulence      on;
    printCoeffs     on;
    kEpsilon
    {
    }
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              0;
    sigmak          1;
    sigmaEps        1.3;
}

Selecting phaseModel for air: purePhaseModel
Selecting diameterModel for phase air: constant
Selecting thermodynamics package 
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Calculating face flux field phi.air
Selecting turbulence model type laminar
Selecting laminar stress model Stokes
Selecting phaseModel for tracer: purePhaseModel
Selecting diameterModel for phase tracer: constant
Selecting thermodynamics package 
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState rhoConst;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Calculating face flux field phi.tracer
Selecting turbulence model type laminar
Selecting laminar stress model Stokes
No MRF models present

Selecting default blending method: none
Selecting surfaceTensionModel for (tracer and air): constant
Selecting surfaceTensionModel for (particles and air): constant
Selecting dragModel for (tracer in air): SyamlalOBrien
Selecting swarmCorrection for (tracer in air): none
Selecting dragModel for (particles in air): SyamlalOBrien
Selecting swarmCorrection for (particles in air): none
Selecting virtualMassModel for (tracer in air): constantCoefficient
Selecting virtualMassModel for (particles in air): constantCoefficient
Selecting heatTransferModel for (tracer in air): RanzMarshall
Selecting heatTransferModel for (particles in air): RanzMarshall
Calculating field g.h

Reading field p_rgh

Courant Number mean: 1e-05 max: 1e-05

Starting time loop

Courant Number mean: 1e-05 max: 1e-05
deltaT = 5.99988e-06
Time = 5.99988e-06

PIMPLE: Iteration 1
MULES: Solving for alpha.particles
MULES: Solving for alpha.air
MULES: Solving for alpha.tracer
particles fraction, min, max = 0.16 0 0.6
air fraction, min, max = 0.84 0.4 1
tracer fraction, min, max = 0 0 0
Phase-sum volume fraction, min, max = 1 1 1
MULES: Solving for alpha.particles
MULES: Solving for alpha.air
MULES: Solving for alpha.tracer
particles fraction, min, max = 0.16 0 0.6
air fraction, min, max = 0.84 0.4 1
tracer fraction, min, max = 0 0 0
Phase-sum volume fraction, min, max = 1 1 1
MULES: Solving for alpha.particles
MULES: Solving for alpha.air
MULES: Solving for alpha.tracer
particles fraction, min, max = 0.16 0 0.6
air fraction, min, max = 0.84 0.4 1
tracer fraction, min, max = 0 0 0
Phase-sum volume fraction, min, max = 1 1 1
Constructing momentum equations
smoothSolver:  Solving for e.particles, Initial residual = 0.918774, Final residual = 1.19097e-05, No Iterations 1000
smoothSolver:  Solving for e.air, Initial residual = 1, Final residual = 2.04943e-09, No Iterations 1
smoothSolver:  Solving for e.tracer, Initial residual = 0.00137557, Final residual = 0, No Iterations 1
particles min/max T 600 - 600
air min/max T 300 - 600
tracer min/max T 600 - 600
GAMG:  Solving for p_rgh, Initial residual = 1.84941e-05, Final residual = 2.84545e-15, No Iterations 1
PIMPLE: Iteration 2
MULES: Solving for alpha.particles
MULES: Solving for alpha.air
MULES: Solving for alpha.tracer
particles fraction, min, max = 0.16 0 0.6
air fraction, min, max = 0.84 0.4 1
tracer fraction, min, max = 0 0 0
Phase-sum volume fraction, min, max = 1 1 1
MULES: Solving for alpha.particles
MULES: Solving for alpha.air
MULES: Solving for alpha.tracer
particles fraction, min, max = 0.16 0 0.6
air fraction, min, max = 0.84 0.4 1
tracer fraction, min, max = 0 0 0
Phase-sum volume fraction, min, max = 1 1 1
MULES: Solving for alpha.particles
MULES: Solving for alpha.air
MULES: Solving for alpha.tracer
particles fraction, min, max = 0.16 0 0.6
air fraction, min, max = 0.84 0.4 1
tracer fraction, min, max = 0 0 0
Phase-sum volume fraction, min, max = 1 1 1
Constructing momentum equations
smoothSolver:  Solving for e.particles, Initial residual = 0.482435, Final residual = 1.20091e-05, No Iterations 1000
smoothSolver:  Solving for e.air, Initial residual = 0.472375, Final residual = 1.37178e-10, No Iterations 1
smoothSolver:  Solving for e.tracer, Initial residual = 1, Final residual = 0, No Iterations 1
particles min/max T 600 - 600
air min/max T 300 - 600
tracer min/max T 599.993 - 600
GAMG:  Solving for p_rgh, Initial residual = 8.75304e-06, Final residual = 2.56878e-15, No Iterations 1
PIMPLE: Iteration 3
MULES: Solving for alpha.particles
MULES: Solving for alpha.air
MULES: Solving for alpha.tracer
particles fraction, min, max = 0.16 0 0.6
air fraction, min, max = 0.84 0.4 1
tracer fraction, min, max = 0 0 0
Phase-sum volume fraction, min, max = 1 1 1
MULES: Solving for alpha.particles
MULES: Solving for alpha.air
MULES: Solving for alpha.tracer
particles fraction, min, max = 0.16 0 0.6
air fraction, min, max = 0.84 0.4 1
tracer fraction, min, max = 0 0 0
Phase-sum volume fraction, min, max = 1 1 1
MULES: Solving for alpha.particles
MULES: Solving for alpha.air
MULES: Solving for alpha.tracer
particles fraction, min, max = 0.16 0 0.6
air fraction, min, max = 0.84 0.4 1
tracer fraction, min, max = 0 0 0
Phase-sum volume fraction, min, max = 1 1 1
Constructing momentum equations
smoothSolver:  Solving for e.particles, Initial residual = 0.0803322, Final residual = 2.31035e-06, No Iterations 1000
smoothSolver:  Solving for e.air, Initial residual = 0.00124982, Final residual = 3.93254e-11, No Iterations 1
smoothSolver:  Solving for e.tracer, Initial residual = 0.585601, Final residual = 0, No Iterations 1
particles min/max T 600 - 600
air min/max T 300 - 600
tracer min/max T 599.945 - 600
GAMG:  Solving for p_rgh, Initial residual = 7.00756e-08, Final residual = 1.91392e-15, No Iterations 1

[louisgag]

Hello Rémi,
did you try turning on turbulence only after a few laminar timesteps?

[chengineer]

I have exactly the same issue. Have you found a solution to this problem?

[geth03]

in my case the solver crashes when the conti phase void fraction approaches 0.
i only use RANS for the conti phase though.