[Davyd]

Please could someone explain me the reason why the simulation is not running?

I am working with OpenFOAM-v1912 and using reactingTwoPhaseEulerFoam solver. I am applying reactingTwoPhaseEulerFoam tutorial (laminar fluidised bed) modified using chemistry from other reactingTwoPhaseEulerFoam tutorial - bubbleColumnEvaporatingReacting.

Here is the result I got:

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / 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
LAPTOP-LTH2ON9V:OpenFOAM/OpenFOAM-v1912/RUN/fluidisedBed_REACTING_bcolumn$


After that the simulation is stopped.

Who knows the reason or at least possible reason your suggestions are highly appreciated.

[geth03]

can you share your phaseProperties-file input?

[Davyd]

Sure I can:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v1912 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object phaseProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

type interfaceCompositionPhaseChangeTwoPhaseSystem;

phases (air particles);

air
{
type reactingPhaseModel;
diameterModel isothermal;
isothermalCoeffs
{
d0 1e-4;
p0 1e5;
}

Sct 1.4;
residualAlpha 0;
}

particles
{
type purePhaseModel;
diameterModel constant;
constantCoeffs
{
d 1e-4;
}
alphaMax 0.62;
residualAlpha 1e-6;
}


blending
{
default
{
type linear;
minFullyContinuousAlpha.air 0.4;
minPartlyContinuousAlpha.air 0;
minFullyContinuousAlpha.particles 0.4;
minPartlyContinuousAlpha.particles 0;
}

heatTransfer
{
type linear;
minFullyContinuousAlpha.air 0.4;
minPartlyContinuousAlpha.air 0;
minFullyContinuousAlpha.particles 0.4;
minPartlyContinuousAlpha.particles 0;
}

massTransfer
{
type linear;
minFullyContinuousAlpha.air 0.4;
minPartlyContinuousAlpha.air 0;
minFullyContinuousAlpha.particles 0.4;
minPartlyContinuousAlpha.particles 0;
}
}

surfaceTension
(
(air and particles)
{
type constant;
sigma 0;
}
);

aspectRatio
(
(particles in air)
{
type constant;
E0 1.0;
}
);

drag
(
(air in particles)
{
type GidaspowErgunWenYu;
residualRe 1e-3;
swarmCorrection
{
type none;
}
}
(particles in air)
{
type GidaspowErgunWenYu;
residualRe 1e-3;
swarmCorrection
{
type none;
}
}
);

virtualMass
(
(air in particles)
{
type constantCoefficient;
Cvm 0.5;
}
(particles in air)
{
type constantCoefficient;
Cvm 0.5;
}
);

interfaceComposition
(
(particles in air)
{
type saturated;
species ( NH3 );
Le 1.79;
saturationPressure
{
type ArdenBuck;
}
}
);

heatTransfer.air
(
(air in particles)
{
type spherical;
residualAlpha 1e-3;
}
(particles in air)
{
type RanzMarshall;
residualAlpha 1e-3;
}
);

heatTransfer.particles
(
(air in particles)
{
type RanzMarshall;
residualAlpha 1e-3;
}
(particles in air)
{
type spherical;
residualAlpha 1e-3;
}
);

massTransfer.air
(
(air in particles)
{
type spherical;
Le 1.79;
}
(particles in air)
{
type Frossling;
Le 1.79;
}
);

massTransfer.particles
(
(air in particles)
{
type Frossling;
Le 1.79;
}
(particles in air)
{
type spherical;
Le 1.79;
}
);

phaseTransfer
(
);

lift
(
);

wallLubrication
(
);

turbulentDispersion
(
);

// Minimum allowable pressure
pMin 10000;


// ************************************************** *********************** //

[geth03]

well, from your first post, i can see that the solver stops when the phaseproperties dict is been read.
i can't judge if you incorporated everything correctly though.

the simulation doesn't even start, so i won't think that its something related to BC's or IC's.

if something you entered is not in the options list, then openfoam should normally tell, which options are possible.

those crashes without any error output are hard to find the origin of the failure.

i suggest that you start with the example tutorial and then change the settings step by step.at every new step run the solver just to see if it stops working or if its continues. it will be easier to find the cause.

[Davyd]

Ok, thanks for your suggestion.