CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Simulation crashes when turbulence is on (reactingMultiphaseEulerFoam)

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 14, 2019, 08:38
Default Simulation crashes when turbulence is on (reactingMultiphaseEulerFoam)
  #1
New Member
 
Rémi Demol
Join Date: Mar 2016
Location: Nancy, France
Posts: 10
Rep Power: 5
remidemol is on a distinguished road
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.




Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
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
remidemol is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Long simulation time using LES turbulence model? tmy CONVERGE 2 September 7, 2018 09:53
Turbulence Model in axisymmetric nozzle simulation cric92 FLUENT 5 September 4, 2017 23:35
Simulation FPEs - turbulence for transient and steady-state? DaveR OpenFOAM Running, Solving & CFD 5 March 5, 2017 15:06
Turbulence postprocessing Mohsin FLUENT 2 October 3, 2016 14:18
Simulation seems to converge but crashes suddenly xxxx OpenFOAM 16 September 12, 2014 08:07


All times are GMT -4. The time now is 15:53.