[QuentinB7]

Hello,

I am working on the simulation of combustion inside an oven. I use chtMultiRegionFoam but I have problems trying to integrate the combustion into my fluid region. I took the same files and conditions as the SmallPoolFire3D tutorial.

I find myself with this error:


--> FOAM FATAL ERROR:

request for volScalarField ph_rgh from objectRegistry Air failed
available objects of type volScalarField are

123
(
aLambda_0
I10
a
thermo:mu
thermosi
qem12
qin2
h
I1
fres_CO2
CO2
qr0
qin9
Qdot
G
ILambda_2_0
qr2
qin10
qr
dpdt
O2_0
ILambda_1_0
I7
ILambda_11_0
alphat
qin7
I3
qin15
CH4_0
qem9
qr7
qem13
wFuel
p
T
I8
qem6
qem1
ILambda_13_0
qr10
fres_N2
H2O
I14
ILambda_12_0
I12
I5
fres_H2O
qr3
qem10
ILambda_7_0
ILambda_4_0
ILambda_5_0
I13
fres_O2
H2O_0
fres_CH4
qin
qem5
I11
I15
ILambda_0_0
qem8
N2
nut
qin11
K
I4
qr14
ILambda_10_0
ILambda_3_0
K_0
CO2_0
ILambda_14_0
qr15
qr1
rho
rAU
k
qin6
I9
h_0
qem2
I6
O2
p_rgh
qem14
qin4
qin8
qr12
gh
qem11
ILambda_8_0
CH4
I2
qem0
qin1
qr9
ILambda_15_0
qr13
rho_0
qem
I0
qr5
qem4
qin13
qin14
qr11
qr6
ILambda_9_0
qr8
ILambda_6_0
bLambda_0
qin3
qem3
qem15
thermo:rho
qem7
qin12
qr4
epsilon
qin0
thermo:alpha
qin5
)


From function const Type& Foam:bjectRegistry::lookupObject(const Foam::word&) const [with Type = Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>]
in file /home/bulardq/OpenFOAM/OpenFOAM-6/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 193.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const& Foam:bjectRegistry::lookupObject<Foam::Geometric Field<double, Foam::fvPatchField, Foam::volMesh> >(Foam::word const&) const at ??:?
#3 Foam:rghTotalHydrostaticPressureFvPatchScalarFie ld::updateCoeffs() at ??:?
#4 Foam::fvMatrix<double>::fvMatrix(Foam::GeometricFi eld<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) at ??:?
#5 Foam::fv::EulerDdtScheme<double>::fvmDdt(Foam::Geo metricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#6 ? at ??:?
#7 ? at ??:?
#8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9 ? at ??:?
Abandon (core dumped)
bulardq@c-2014-074:~/Op


Has anyone ever come across this problem? Can I try to integrate the combustion into chtMultiRegionFoam without having to recreate a new application?

Thank you in advance for your answers.

[Yann]

Hello Quentin,

As far as I know, combustion isn't implemented in chtMultiRegionFoam.
Have a look at this thread : Can chtMultiRegionFoam do combustion, radiation and buoyancy?

About your error, the interesting part is here :

--> FOAM FATAL ERROR: 

    request for volScalarField ph_rgh from objectRegistry Air failed

The solver is looking for a variable named "ph_rgh" in the "Air" region and it seems it cannot find it.
Since you started from the smallPoolFire3D, the need for ph_rgh is probably related to the "prghTotalHydrostaticPressure" boundary condition for p_rgh.

First thing is to check your boundary conditions to see if you have the appropriate "ph_rgh" file in the right directory. (or modify your boundary conditions accordingly)


Yann

[QuentinB7]

Hello Yann,

Thank you for your quick reply.

My variable "ph_rgh" is in the good folder "0". At the level of my "p_rgh", "prghTotalHydrostaticPressure" is the boundary condition that I use.

There may be some changes to be made in terms of fvSchemes or fvSolution?


Quentin

[Yann]

It's not easy to help you without more information.
The ph_rgh file should be located in 0/Air/ph_rgh or if running in parallel : processorN/Air/ph_rgh.


Can you post a log file showing the error and if possible a minimal case with your setup and a run script to reproduce what you do?

[QuentinB7]

Here is the required information

https://filesender.renater.fr/?s=dow...b-f483edb60bf9

[clapointe]

Firefoam has a ph_rgh equation that's called as the silver spins up to computd a ph_rgh field that is then used by the accompanying boundary condition. If it is not computed (which it clearly is not here, hence the error), the boundary condition will return an error.

Caelan