Run reactingTwoPhaseEulerFoam with IATE source option phaseChange , The error displays as follows
--> FOAM FATAL ERROR:
request for volScalarField iDmdt.Su.gasInLiquid from objectRegistry region0 failed
available objects of type volScalarField are
47
(
iDmdtSu.liquidInGas
T.liquid
alpha.gas
max((0.5*average(interpolate((alpha.gas+alpha.gas_0)))),residualAlpha)
continuityErrorSources.gas
thermo:alpha.gas
continuityErrorFlow.liquid
thermo:psi.gas
thermo:mu.liquid
iDmdtSu.gasInLiquid
k.liquid
thermo:mu.gas
nut.liquid
alpha.liquid
nut.gas
Tf.gasAndLiquid
Kd.gasAndLiquid
d.gas
e.gas
p_rgh
epsilon.gas
dpdt
air.gas
alphat.liquid
gh
thermo:rho.liquid
iDmdtSp.gasInLiquid
k.gas
alpha.gas_0
YSolvent.liquidInGas
continuityErrorFlow.gas
T.gas
continuityErrorSources.liquid
water.gas
kappai.gas
p
iDmdtSp.liquidInGas
p_0
water.liquid
thermo:psi.liquid
thermo:alpha.liquid
e.liquid
epsilon.liquid
air.liquid
thermo:rho.gas
Vm.gasAndLiquid
alphat.gas
)
From function const Type& Foam::objectRegistry::lookupObject(const Foam::word&, bool) const [with Type = Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>]
in file /OpenFOAM/OpenFOAM-v1912/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 463.
FOAM aborting
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const& Foam::objectRegistry::lookupObject<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> >(Foam::word const&, bool) const at ??:?
#3 Foam::diameterModels::IATEsources::phaseChange::R(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&) const at ??:?
#4 Foam::diameterModels::IATE::correct() at ??:?
#5 Foam::MovingPhaseModel<Foam::ThermoPhaseModel<Foam::phaseModel, Foam::rhoReactionThermo> >::correct() at ??:?
#6 Foam::phaseSystem::correct() at ??:?
#7 Foam::PhaseTransferPhaseSystem<Foam::TwoResistanceHeatTransferPhaseSystem<Foam::MomentumTransferPhaseSystem<Foam::twoPhaseSystem> > >::correct() at ??:?
#8 ? at ??:?
#9 __libc_start_main in /lib64/libc.so.6
#10 ? at ??:?
Aborted (core dumped)
The gas part of IATA setting for the case’s constant/phaseProperties
gas
{
diameterModel IATE;
IATECoeffs
{
dMax 1e-2;
dMin 1e-4;
residualAlpha 1e-6;
sources
(
phaseChange
{
pairName Su.gasInLiquid;
}
wakeEntrainmentCoalescence
{
Cwe 0.002;
}
randomCoalescence
{
Crc 0.04;
C 3;
alphaMax 0.75;
}
turbulentBreakUp
{
Cti 0.085;
WeCr 6;
}
);
}
residualAlpha 1e-6;
}
If the phaseChange source is removed, reactingTwoPhaseEulerFoam works.