[Davyd]

Hello.
I am trying to do simulation of the gas-solid fluidized bed catalytic reactor. For this purpose I use tutorial reactingTwoPhaseEulerFoam--->laminar--->fluidisedBed. However, since it does not contain any files responsible for chemistry, I use chemistry files from combustion tutorial--->chemFoam--->h2. Of course I have modified them according to my case (reaction, species, conditions).
When I start simulation I get the ERROR described below, however it seems that I have defined the mole fractions for species. You can find the files responsible for chemistry attached to this thread.
Please, if someone knows how to overcome this issue, give me your advice. Thank you for your help!!!
Best,
Davyd.


/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v1912 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : _f3950763fe-20191219 OPENFOAM=1912
Arch : "LSB;label=32;scalar=64"
Exec : reactingTwoPhaseEulerFoam
Date : Jan 28 2021
Time : 15:44:49
Host : LAPTOP-LTH2ON9V
PID : 1922
I/O : uncollated
Case : /mnt/c/Users/David/Downloads/OpenFOAM/OpenFOAM-v1912/RUN/fluidisedBed_REACTING_comb_h2
nProcs : 1
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0


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


Reading g

Reading hRef
Creating phaseSystem

Selecting twoPhaseSystem interfaceCompositionPhaseChangeTwoPhaseSystem
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;
}

Selecting phaseModel for air: reactingPhaseModel
Selecting diameterModel for phase air: isothermal
Selecting thermodynamics package
{
type heRhoThermo;
mixture reactingMixture;
transport sutherland;
thermo janaf;
equationOfState perfectGas;
specie specie;
energy sensibleEnthalpy;
}

Selecting chemistryReader chemkinReader


--> FOAM FATAL ERROR:
Sum of mass fractions is zero for species 6(NO NH3 O2 CO2 N2 H2O)

From function void Foam::multiComponentMixture<ThermoType>::correctMa ssFractions() [with ThermoType = Foam::sutherlandTransport<Foam::species::thermo<Fo am::janafThermo<Foam:erfectGas<Foam::specie> >, Foam::sensibleEnthalpy> >]
in file lnInclude/multiComponentMixture.C at line 64.

FOAM exiting