| sepponen |
August 1, 2016 05:19 |
reactingTwoPhaseEulerFoam modifying bubbleColumnEvaporatingReacting case
1 Attachment(s)
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:
Code:
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 |
