[infrared]

Hello,

I performed an openfoam case but I got an error message:

Create mesh for time = 0


PIMPLE: Operating solver in PISO mode

Reading thermophysical properties

Selecting thermodynamics package
{
type hePsiThermo;
mixture reactingMixture;
transport sutherland;
thermo janaf;
equationOfState perfectGas;
specie specie;
energy sensibleInternalEnergy;
}



--> FOAM FATAL IO ERROR: (openfoam-2312 patch=240220)
Attempt to return dictionary entry as a primitive

file: constant/thermophysicalProperties/chemistryReader at line 33.

From virtual Foam::ITstream& Foam::dictionaryEntry::stream() const
in file db/dictionary/dictionaryEntry/dictionaryEntry.C at line 72.

FOAM aborting

[stack trace]
=============
#1 Foam::error::simpleExit(int, bool) in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#2 Foam::dictionaryEntry::stream() const in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#3 ? in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libreactionThermophysicalModels.so
#4 ? in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libreactionThermophysicalModels.so
#5 ? in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libreactionThermophysicalModels.so
#6 ? in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libreactionThermophysicalModels.so
#7 Foam:siReactionThermo::addfvMeshConstructorToTab le<Foam::hePsiThermo<Foam:siReactionThermo, Foam::SpecieMixture<Foam::reactingMixture<Foam::su therlandTransport<Foam::species::thermo<Foam::jana fThermo<Foam:erfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > > > >::New(Foam::fvMesh const&, Foam::word const&) in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libreactionThermophysicalModels.so
#8 ? in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libreactionThermophysicalModels.so
#9 Foam:siReactionThermo::New(Foam::fvMesh const&, Foam::word const&) in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libreactionThermophysicalModels.so
#10 ? in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/bin/reactingFoam
#11 ? in /lib/x86_64-linux-gnu/libc.so.6
#12 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#13 ? in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/bin/reactingFoam
=============
Aborted

Here is the thermophysical Properties:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2312 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object thermophysicalProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

thermoType
{
type hePsiThermo;
mixture reactingMixture;
transport sutherland;
thermo janaf;
equationOfState perfectGas;
specie specie;
energy sensibleInternalEnergy;
}

CHEMKINFile "$FOAM_CASE/constant/chem.inp";
CHEMKINThermoFile "$FOAM_CASE/constant/therm.dat";

chemistryReader
{
type foamChemistryReader;
}

reactingMixture
{
specieCoeffs
{
O2
{
molWeight 32.0;
}
N2
{
molWeight 28.0;
}
H2
{
molWeight 2.016;
}
H2O
{
molWeight 18.01528;
}
CO2
{
molWeight 44.01;
}
CH4
{
molWeight 16.04;
}
// Add other species as needed
}

thermodynamicCoeffs
{
H2
{
Tlow 200.0;
Thigh 6000.0;
Tcommon 1000.0;
lowCpCoeffs (3.298124 0.000824 -0.000308 0.000101 -0.000018);
highCpCoeffs (2.991423 0.000007 -0.000001 0.000000 -0.000000);
molWeight 2.016;
}
O2
{
Tlow 200.0;
Thigh 6000.0;
Tcommon 1000.0;
lowCpCoeffs (3.782456 -0.000293 0.000636 -0.000276 0.000056);
highCpCoeffs (3.660960 -0.000070 0.000197 -0.000081 0.000016);
molWeight 32.0;
}
N2
{
Tlow 200.0;
Thigh 6000.0;
Tcommon 1000.0;
lowCpCoeffs (3.531005 -0.000123 0.000254 -0.000142 0.000032);
highCpCoeffs (2.952576 0.000016 -0.000008 0.000002 -0.000000);
molWeight 28.0;
}
H2O
{
Tlow 200.0;
Thigh 6000.0;
Tcommon 1000.0;
lowCpCoeffs (4.198640 -0.002036 0.000307 -0.000097 0.000016);
highCpCoeffs (3.033992 0.000030 -0.000016 0.000005 -0.000001);
molWeight 18.01528;
}
CO2
{
Tlow 200.0;
Thigh 6000.0;
Tcommon 1000.0;
lowCpCoeffs (2.356773 0.008984 -0.000007 0.000002 -0.000000);
highCpCoeffs (4.636511 0.000072 -0.000022 0.000005 -0.000001);
molWeight 44.01;
}
CH4
{
Tlow 200.0;
Thigh 6000.0;
Tcommon 1000.0;
lowCpCoeffs (5.149876 -0.013671 0.003131 -0.000658 0.000122);
highCpCoeffs (1.485849 0.002366 -0.000432 0.000083 -0.000013);
molWeight 16.04;
}
// Add other species as needed
}

transportCoeffs
{
H2
{
As 7.76644e-08;
Ts 20.33;
muCoeffs (2.21e-7 0.0 0.67);
kappaCoeffs (0.1805 0.0 1.1);
}
O2
{
As 3.45846e-08;
Ts 97;
muCoeffs (1.8e-7 0.0 0.67);
kappaCoeffs (0.025 0.0 1.1);
}
N2
{
As 6.55987e-08;
Ts 71.4;
muCoeffs (1.89e-7 0.0 0.67);
kappaCoeffs (0.026 0.0 1.1);
}
H2O
{
As 1.73478e-08;
Ts 809.1;
muCoeffs (2.0e-7 0.0 0.67);
kappaCoeffs (0.6065 0.0 1.1);
}
CO2
{
As 2.148e-08;
Ts 244;
muCoeffs (1.5e-7 0.0 0.67);
kappaCoeffs (0.0165 0.0 1.1);
}
CH4
{
As 2.593e-08;
Ts 148.6;
muCoeffs (1.27e-7 0.0 0.67);
kappaCoeffs (0.034 0.0 1.1);
}
// Add other species as needed
}
}


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

Do you have any ideas on how to fix this issue?

I really apreciate your help. Thank you.

[Yann]

Hello,

Are you sure of your thermophysicalProperties file syntax?
Where did you get this syntax? It looks pretty different from the tutorials.

Regards,
Yann

[infrared]

Thank you for your reply. Actually, I am a newbie using OpenFOAM for simulations. After reading the guidelines, I tried to establish the thermophysicalProperties syntax with AI as well. Is there anything wrong? Hopefully, I will get insights by asking in this forum.

Please give me some suggestions.

[Yann]

Don't use AI for this. Have a look at tutorials for the solver you want to use, and adjust the files to adapt it to the simulation you want to run.

OpenFOAM is shipped with tutorials, look there. If you look for stuff online, make sure you are looking at resources corresponding to the OpenFOAM version you are using. (other versions/branches might use different syntaxes)

Yann

[infrared]

Hello Yann,

Thank you for your suggestions. I decided to adopt gri case with the solver chemFoam from the tutorials. Then make necessary adjustment.

The directory structure is
0
|-- alphat
|-- CH4
|-- CO2
|-- H2O
|-- N2
|-- O2
|-- p
|-- rho
|-- T
|-- U
|-- Ydefault

constant
|-- polyMesh
|-- boundary
|-- faces
|-- neighbour
|-- owner
|-- points
|-- chemistryProperties
|-- combustionProperties
|-- initial condition
|-- reactions
|-- thermo.compressibleGas
|-- thermophysicalProperties
|-- turbulenceProperties

system
|-- blockMeshDict
|-- controlDict
|-- decomposeParDict
|-- fvSchemes
|-- fvSolution

chemkin
|-- chem.inp
|-- senk.inp
|-- senk.out
|-- therm.dat
|-- transportProperties

But it shows:
Creating base fields for time 0
Reading thermophysical properties

Selecting thermodynamics package
{
type hePsiThermo;
mixture reactingMixture;
transport sutherland;
thermo janaf;
energy sensibleEnthalpy;
equationOfState perfectGas;
specie specie;
}

Selecting chemistryReader foamChemistryReader
Selecting chemistry solver
{
solver ode;
method standard;
}

StandardChemistryModel: Number of species = 5 and reactions = 1
Selecting ODE solver seulex
Reading/calculating face flux field phi

[stack trace]
=============
#1 Foam::sigFpe::sigHandler(int) in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#2 ? in /lib/x86_64-linux-gnu/libc.so.6
#3 ? in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/bin/chemFoam
#4 ? in /lib/x86_64-linux-gnu/libc.so.6
#5 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#6 ? in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/bin/chemFoam
=============
Floating point exception

Do you any suggestions to counter this issue? Thank you for your help.

[infrared]

Hello Yann,

I would like to ask you about the combustion visualization. I adopted the 'counterFlowFlame2D' case from the tutorials and modified the domain and its conditions. In general, my case is a cylindrical domain with a single inlet placed at the top wall in an impinging setup. I successfully performed the new case, but the visualizations seem weird. Please see the attachments.





I use diesel (C12H26) as the inlet fuel, and the case is still ongoing since it has not reached the impingement yet, but I find the visualization strange. What do you think about it?

Thank you for your help.

Best regards,
IR

[Yann]

Hello IR,

It seems you forgot or failed to attach the screenshots!

Yann

[infrared]

Can you see this picture?

http://https://www.cfd-online.com/Fo...1-diesel-t.png

http://https://www.cfd-online.com/Fo...2-diesel-u.png

[Yann]

It seems I cannot