|
[Sponsors] |
TransportProperties and chemkinReader problem |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
New Member
Rohmawati Iis
Join Date: Apr 2014
Posts: 20
Rep Power: 12 ![]() |
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: ![]() ![]() ![]() #8 ? in /usr/lib/openfoam/openfoam2312/platforms/linux64GccDPInt32Opt/lib/libreactionThermophysicalModels.so #9 Foam: ![]() #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. |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Senior Member
Yann
Join Date: Apr 2012
Location: France
Posts: 1,169
Rep Power: 27 ![]() |
Hello,
Are you sure of your thermophysicalProperties file syntax? Where did you get this syntax? It looks pretty different from the tutorials. Regards, Yann |
|
![]() |
![]() |
![]() |
![]() |
#3 |
New Member
Rohmawati Iis
Join Date: Apr 2014
Posts: 20
Rep Power: 12 ![]() |
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. |
|
![]() |
![]() |
![]() |
![]() |
#4 |
Senior Member
Yann
Join Date: Apr 2012
Location: France
Posts: 1,169
Rep Power: 27 ![]() |
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 |
|
![]() |
![]() |
![]() |
![]() |
#5 |
New Member
Rohmawati Iis
Join Date: Apr 2014
Posts: 20
Rep Power: 12 ![]() |
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. |
|
![]() |
![]() |
![]() |
![]() |
#6 |
New Member
Rohmawati Iis
Join Date: Apr 2014
Posts: 20
Rep Power: 12 ![]() |
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 |
|
![]() |
![]() |
![]() |
![]() |
#7 |
Senior Member
Yann
Join Date: Apr 2012
Location: France
Posts: 1,169
Rep Power: 27 ![]() |
Hello IR,
It seems you forgot or failed to attach the screenshots! Yann |
|
![]() |
![]() |
![]() |
![]() |
#8 |
New Member
Rohmawati Iis
Join Date: Apr 2014
Posts: 20
Rep Power: 12 ![]() |
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 |
|
![]() |
![]() |
![]() |
![]() |
#9 |
Senior Member
Yann
Join Date: Apr 2012
Location: France
Posts: 1,169
Rep Power: 27 ![]() |
It seems I cannot
|
|
![]() |
![]() |
![]() |
|
|