[Davyd]

Please, could someone help me with explanation of this error:

Selecting chemistryReader foamChemistryReader
elements not defined in "/mnt/c/Users/David/downloads/OpenFOAM/OpenFOAM-v1912/RUN/fluidisedBed_REACTING/constant/reactions.air".

I have described the reactions.air file based on the bubbleColumnEvaporatingReacting tutorial which is inside RAS directory of a reactingTwoPhaseEulerFoam tutorials. The description is as follows:

species
(
NO
NH3
O2
N2
H2O
AIR
);

reactions
{
waterGasShift
{
type reversibleArrheniusReaction;

reaction "NO^4 + NH3^4 + O2^1 = N2^4 + H2O^6";

A 3.85e6;
beta 0;
Ta 0.6764e4;
}
}


I am trying to run the simulation using reactingTwoPhaseEulerFoam for a fluidised bed reactor. Although it names "reacting" There is no chemistry inside, so I am trying to add chemistry files from bubbleColumnEvaporatingReacting and change them according to my case.


Thank you in advance for any help.

[HPE]

Hi,

- Is there any chance for you to attach the `entire` log file, please, including the error?

[Davyd]

Thanks for you reply.
Now I have passed through this problem and I have another problem the reason of which is not clear for me as well:

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / 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 : May 20 2020
Time : 19:52:24
Host : LAPTOP-LTH2ON9V
PID : 2578
I/O : uncollated
Case : /mnt/c/Users/David/Downloads/OpenFOAM/OpenFOAM-v1912/RUN/fluidisedBed_REACTING_bcolumn
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 air: reactingPhaseModel
Selecting diameterModel for phase air: isothermal
Selecting thermodynamics package
{
type heRhoThermo;
mixture reactingMixture;
transport sutherland;
thermo janaf;
equationOfState perfectGas;
specie specie;
energy sensibleInternalEnergy;
}

Selecting chemistryReader foamChemistryReader
elements not defined in "/mnt/c/Users/David/Downloads/OpenFOAM/OpenFOAM-v1912/RUN/fluidisedBed_REACTING_bcolumn/constant/reactions.air"
Calculating face flux field phi.air
Selecting turbulence model type laminar
Selecting laminar stress model Stokes
Selecting combustion model PaSR
Selecting chemistry solver
{
solver EulerImplicit;
method standard;
}

StandardChemistryModel: Number of species = 7 and reactions = 1
using integrated reaction rate
Selecting phaseModel for particles: purePhaseModel
Selecting diameterModel for phase particles: constant
Selecting thermodynamics package
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo eConst;
equationOfState perfectFluid;
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

Selecting default blending method: linear
Selecting heatTransfer blending method: linear
Selecting massTransfer blending method: linear
Selecting surfaceTensionModel for (air and particles): constant
Selecting aspectRatioModel for (particles in air): constant
Selecting dragModel for (particles in air): GidaspowErgunWenYu
Selecting swarmCorrection for (particles in air): none
Selecting swarmCorrection for (particles in air): none
Selecting swarmCorrection for (particles in air): none
Selecting virtualMassModel for (particles in air): constantCoefficient
Selecting virtualMassModel for (air in particles): constantCoefficient
Selecting heatTransferModel for (particles in air): RanzMarshall
Selecting heatTransferModel for (air in particles): spherical
Selecting heatTransferModel for (particles in air): spherical
Selecting heatTransferModel for (air in particles): RanzMarshall
davyd@LAPTOP-LTH2ON9V:/mnt/c/Users/David/Downloads/OpenFOAM/OpenFOAM-v1912/RUN/fluidisedBed_REACTING_bcolumn$

After that simulation is stopped.
Maybe you can help with this issue.

[Davyd]

When I adjust the phaseProperties file in another way I got such king of error:

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / 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 : May 21 2020
Time : 09:37:11
Host : LAPTOP-LTH2ON9V
PID : 2616
I/O : uncollated
Case : /mnt/c/Users/David/Downloads/OpenFOAM/OpenFOAM-v1912/RUN/fluidisedBed_REACTING_bcolumn_Copy_2020-05-20
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 multiComponentMixture;
transport const;
thermo eConst;
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 const;
thermo eConst;
equationOfState rhoConst;
specie specie;
energy sensibleInternalEnergy;
}

Selecting chemistryReader foamChemistryReader
elements not defined in "/mnt/c/Users/David/Downloads/OpenFOAM/OpenFOAM-v1912/RUN/fluidisedBed_REACTING_bcolumn_Copy_2020-05-20/constant/reactions.air"
Calculating face flux field phi.air
Selecting turbulence model type laminar
Selecting laminar stress model Stokes
Selecting combustion model PaSR
Selecting chemistry solver
{
solver EulerImplicit;
method standard;
}

StandardChemistryModel: Number of species = 7 and reactions = 1
using integrated reaction rate
No MRF models present

Selecting default blending method: none
Selecting surfaceTensionModel for (particles and air): constant
Selecting aspectRatioModel for (particles in air): constant
Selecting aspectRatioModel for (air in particles): constant
Selecting dragModel for (particles in air): GidaspowErgunWenYu
Selecting swarmCorrection for (particles in air): none
Selecting swarmCorrection for (particles in air): none
Selecting swarmCorrection for (particles in air): none
Selecting virtualMassModel for (particles in air): constantCoefficient
Selecting heatTransferModel for (particles in air): spherical
Selecting heatTransferModel for (particles in air): RanzMarshall
Selecting interfaceCompositionModel for (particles in air): saturated<heRhoThermo<multiComponentMixture<const< eConst<rhoConst<specie>>,sensibleInternalEnergy>>> ,heRhoThermo<reactingMixture<const<eConst<rhoConst <specie>>,sensibleInternalEnergy>>>>
Selecting saturationModel: ArdenBuck
Selecting massTransferModel for (particles in air): spherical
Selecting massTransferModel for (particles in air): Frossling
Calculating field g.h

Reading field p_rgh

Courant Number mean: 0.00287522 max: 0.003295

Starting time loop

fieldAverage fieldAverage1:
Restarting averaging for fields:
U.particles: starting averaging at time 0
U.air: starting averaging at time 0
alpha.particles: starting averaging at time 0
p: starting averaging at time 0

Courant Number mean: 0.00287522 max: 0.003295
Max Ur Courant Number = 0.003295
Time = 0.0002

PIMPLE: iteration 1
MULES: Solving for alpha.particles
MULES: Solving for alpha.particles
alpha.particles volume fraction = 0.129996 Min(alpha1) = 0 Max(alpha1) = 0.52
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigSegv::sigHandler(int) at ??:?
#2 ? in /lib/x86_64-linux-gnu/libc.so.6
#3 Foam::InterfaceCompositionModel<Foam::heRhoThermo< Foam::rhoReactionThermo, Foam::SpecieMixture<Foam::multiComponentMixture<Fo am::constTransport<Foam::species::thermo<Foam::eCo nstThermo<Foam::rhoConst<Foam::specie> >, Foam::sensibleInternalEnergy> > > > >, Foam::heRhoThermo<Foam::rhoReactionThermo, Foam::SpecieMixture<Foam::reactingMixture<Foam::co nstTransport<Foam::species::thermo<Foam::eConstThe rmo<Foam::rhoConst<Foam::specie> >, Foam::sensibleInternalEnergy> > > > > >:(Foam::word const&) const at ??:?
#4 Foam::InterfaceCompositionPhaseChangePhaseSystem<F oam::PhaseTransferPhaseSystem<Foam::TwoResistanceH eatTransferPhaseSystem<Foam::MomentumTransferPhase System<Foam::twoPhaseSystem> > > >::massTransfer() const at ??:?
#5 ? at ??:?
#6 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#7 ? at ??:?
Segmentation fault (core dumped)
davyd@LAPTOP-LTH2ON9V:/mnt/c/Users/David/Downloads/OpenFOAM/OpenFOAM-v1912/RUN/fluidisedBed_REACTING_bcolumn_Copy_2020-05-20$

[ERCI]

HI，Davyd
although i havnt see you constant/reactions.air flie

I think in this flie you havnt put element in first

like this ——





elements
(
O
C
H
N
);


species
(
O2
H2O
CH4
CO2
N2
);

reactions
{
methaneReaction
{
type irreversibleArrheniusReaction;
reaction "CH4 + 2O2 = CO2 + 2H2O";
A 5.2e16;
beta 0;
Ta 14906;
}
}
you can see how to define element in Openfoam web
mybe i am wrong ,

hahaha