[sepponen]

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
phasePressureCoeffs
{
    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


--> FOAM FATAL ERROR:
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.

Simo

[pyroWinter]

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