CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

reactingTwoPhaseEulerFoam modifying bubbleColumnEvaporatingReacting case

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

LinkBack Thread Tools Display Modes
Old   August 1, 2016, 05:19
Default reactingTwoPhaseEulerFoam modifying bubbleColumnEvaporatingReacting case
New Member
Simo Sepponen
Join Date: Jun 2016
Posts: 8
Rep Power: 3
sepponen is on a distinguished road
Hello everyone,

I am looking into reactingTwoPhaseEulerFoam and I am trying to find out how to implement some simple chemical reactions in to a case of modelling a circulating fluidised bed cold model. The case has been modeled succesfully with twoPhaseEulerFoam.

My starting point is to look at the tutorial bubbleColumnEvaporatingReacting (Openfoam 3.0.1) and replace the liquid phase with particle phase. I only changed the initial values and phaseproperties to represent the partclee phase. I did no modifications to the reactions in constant/chemistryProperties.gas, constant/reactions.gas and constant/thermo.gas. My idea was to first make the tutorial case work with particle phase instead of liquid phase and after that try to add a more meaningful reaction.

I managed to set up the case so that it starts calculation, but it stops as soon as it starts to calculate the volume fraction of CO.gas. Here is the log:

Create time

Create mesh for time = 0

PIMPLE: no residual control data found. Calculations will employ 5 corrector loops

Reading g

Reading hRef
Creating phaseSystem

Selecting twoPhaseSystem heatAndMomentumTransferTwoPhaseSystem
Selecting phaseModel for gas: reactingPhaseModel
Selecting diameterModel for phase gas: isothermal
Selecting combustion model PaSR<rhoChemistryCombustion>
Selecting chemistry type
    chemistrySolver EulerImplicit;
    chemistryThermo rho;

Selecting thermodynamics package
    type            heRhoThermo;
    mixture         reactingMixture;
    transport       sutherland;
    thermo          janaf;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleInternalEnergy;

Selecting chemistryReader foamChemistryReader
chemistryModel: Number of species = 5 and reactions = 1
    using integrated reaction rate
Calculating face flux field phi.gas
Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
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 phasePressure
    preAlphaExp     500;
    expMax          1000;
    alphaMax        0.62;
    g0              1000;

No MRF models present

No finite volume options present

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

Reading field p_rgh

Courant Number mean: 0.0007164125 max: 0.0007500000019

Starting time loop

Courant Number mean: 0.0007164125 max: 0.0007500000019
Max Ur Courant Number = 0.0007500000019
Time = 0.0001

PIMPLE: iteration 1
MULES: Solving for alpha.gas
MULES: Solving for alpha.gas
MULES: Solving for alpha.gas
alpha.gas volume fraction = 0.9579969532  Min(alpha1) = 0.5  Max(alpha1) = 1

CO.gas not found in table.  Valid entries: 0()

    From function HashTable<T, Key, Hash>::operator[](const Key&)
    in file /shared/apps/openfoam/openmpi-1.8.4-gcc/gcc-4.9.2/3.0.1/src/OpenFOAM/lnInclude/HashTableI.H at line 117.

FOAM exiting
Can someone explain the error? The case is attached in a zip-file. Thank you for your attention.

Attached Files
File Type: zip (91.1 KB, 13 views)
sepponen is offline   Reply With Quote

Old   June 2, 2017, 12:29
New Member
David Buentello
Join Date: Jan 2017
Posts: 7
Rep Power: 2
pyroWinter is on a distinguished road
Hello, I know this is an old post, but did you manage to get this working? I am currently having the same issue. I would greatly appreciate any feedback
pyroWinter is offline   Reply With Quote


Thread Tools
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
Is Playstation 3 cluster suitable for CFD work hsieh OpenFOAM 9 August 16, 2015 14:53
Superlinear speedup in OpenFOAM 13 msrinath80 OpenFOAM Running, Solving & CFD 18 March 3, 2015 06:36
CFD Post - How to check for case in an expression pilakin ANSYS 0 September 26, 2014 04:44
Transient case running with a super computer microfin FLUENT 0 March 31, 2009 11:20
Turbulent Flat Plate Validation Case Jonas Larsson Main CFD Forum 0 April 2, 2004 10:25

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