CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Pre-Processing

using chemkin

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 15, 2014, 14:30
Default using chemkin
  #1
New Member
 
Josh Dawson
Join Date: Oct 2014
Location: Justin, TX
Posts: 19
Rep Power: 11
JMDag2004 is on a distinguished road
Hello everyone! I'm attempting to use chemkin but I'm getting the following error:

jmdag1981@jmdag1981-Gazelle-Professional:~/OpenFOAM/jmdag1981-2.3.0/run/EAPDR/EAPDR$ blockMesh
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.3.0-f5222ca19ce6
Exec : blockMesh
Date : Dec 15 2014
Time : 13:23:09
Host : "jmdag1981-Gazelle-Professional"
PID : 5050
Case : /home/jmdag1981/OpenFOAM/jmdag1981-2.3.0/run/EAPDR/EAPDR
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Creating block mesh from
"/home/jmdag1981/OpenFOAM/jmdag1981-2.3.0/run/EAPDR/EAPDR/constant/polyMesh/blockMeshDict"
Creating curved edges
Creating topology blocks
Creating topology patches

Creating block mesh topology

Check topology

Basic statistics
Number of internal faces : 10
Number of boundary faces : 40
Number of defined boundary faces : 40
Number of undefined boundary faces : 0
Checking patch -> block consistency

Creating block offsets
Creating merge list .

Creating polyMesh from blockMesh
Creating patches
Creating cells
Creating points with scale 0.46
Block 0 cell size :
i : 0.00248768 .. 0.000497535
j : 0.000685672 .. 0.00205702
k : 0.0046 .. 0.0046

Block 1 cell size :
i : 0.00129373 .. 0.00129373
j : 0.000685672 .. 0.00205702
k : 0.0046 .. 0.0046

Block 2 cell size :
i : 0.00126691 .. 0.00190036
j : 0.000533293 .. 0.00159988
k : 0.0046 .. 0.0046

Block 3 cell size :
i : 0.00108267 .. 0.00216534
j : 0.000736823 .. 0.000736823
k : 0.0046 .. 0.0046

Block 4 cell size :
i : 0.00108267 .. 0.00216534
j : 0.00100004 .. 0.00100004
k : 0.0046 .. 0.0046

Block 5 cell size :
i : 0.00043858 .. 0.00131574
j : 0.00100004 .. 0.00100004
k : 0.0046 .. 0.0046

Block 6 cell size :
i : 0.000254255 .. 0.000762764
j : 0.000533293 .. 0.00159988
k : 0.0046 .. 0.0046

Block 7 cell size :
i : 0.000254255 .. 0.000762764
j : 0.00050048 .. 0.00050048
k : 0.0046 .. 0.0046

Block 8 cell size :
i : 0.000254255 .. 0.000762764
j : 0.000736823 .. 0.000736823
k : 0.0046 .. 0.0046

Block 9 cell size :
i : 0.000254255 .. 0.000762764
j : 0.00100004 .. 0.00100004
k : 0.0046 .. 0.0046


Adding cell zones
0 zoneZero
1 zoneOne
2 zoneTwo
3 zoneThree
4 zoneFour
5 zoneFive
6 zoneNine
7 zoneEight
8 zoneSeven
9 zoneSix

Writing cell zones as cellSets

Writing polyMesh
----------------
Mesh Information
----------------
boundingBox: (0 0 -0.0023) (0.500002 0.089999 0.0023)
nPoints: 123102
nCells: 60634
nFaces: 243452
nInternalFaces: 120352
----------------
Patches
----------------
patch 0 (start: 120352 size: 72) name: inletLeftBoundaryZoneZero
patch 1 (start: 120424 size: 208) name: fixedWallsTopBoundaryZoneZero_ZoneTwo
patch 2 (start: 120632 size: 20) name: symmetryPlaneBottomBoundaryZoneZero
patch 3 (start: 120652 size: 60) name: wedgeFaceBottomBoundaryZoneOne
patch 4 (start: 120712 size: 128) name: detonationTubeCasingBottomBoundaryZoneTwo
patch 5 (start: 120840 size: 2) name: backFaceDetonationTubeCasingLeftBoundaryZoneEight
patch 6 (start: 120842 size: 50) name: fixedWallInitiatorTubeTopBoundaryZoneFive
patch 7 (start: 120892 size: 5) name: fixedWallInitiatorTubeLeftBoundaryZoneThree
patch 8 (start: 120897 size: 50) name: symmetryPlaneInitiatorTubeBottomBoundaryZoneFive
patch 9 (start: 120947 size: 128) name: fixedWallDetonationTubeTopBoundaryZoneThree
patch 10 (start: 121075 size: 19) name: fixedWallDetonationTubeLeftBoundaryZoneThree
patch 11 (start: 121094 size: 128) name: symmetryPlaneDetonationTubeBottomBoundaryZoneFour
patch 12 (start: 121222 size: 98) name: outletRightBoundaryZoneSix_ZoneNine
patch 13 (start: 121320 size: 432) name: fixedWallsTopBoundaryZoneNine
patch 14 (start: 121752 size: 432) name: symmetryPlaneBottomBoundaryZoneSix
patch 15 (start: 122184 size: 121268) name: frontAndBack

End

jmdag1981@jmdag1981-Gazelle-Professional:~/OpenFOAM/jmdag1981-2.3.0/run/EAPDR/EAPDR$ setFields
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.3.0-f5222ca19ce6
Exec : setFields
Date : Dec 15 2014
Time : 13:23:16
Host : "jmdag1981-Gazelle-Professional"
PID : 5052
Case : /home/jmdag1981/OpenFOAM/jmdag1981-2.3.0/run/EAPDR/EAPDR
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading setFieldsDict

Setting field region values
Adding all cells of cellZone zoneFive ...
Found matching zone zoneFive with 250 cells.
Setting internal values of volScalarField H2
Setting internal values of volScalarField O2
Setting internal values of volScalarField N2
Setting internal values of volVectorField U
Adding all cells of cellZone zoneFour ...
Found matching zone zoneFour with 640 cells.
Setting internal values of volScalarField H2
Setting internal values of volScalarField O2
Setting internal values of volScalarField N2
Setting internal values of volVectorField U
Adding all cells of cellZone zoneThree ...
Found matching zone zoneThree with 2432 cells.
Setting internal values of volScalarField H2
Setting internal values of volScalarField O2
Setting internal values of volScalarField N2
Setting internal values of volVectorField U
Adding all cells of cellZone zoneZero ...
Found matching zone zoneZero with 1440 cells.
Setting internal values of volScalarField O2
Setting internal values of volScalarField N2
Adding all cells of cellZone zoneOne ...
Found matching zone zoneOne with 4320 cells.
Setting internal values of volScalarField O2
Setting internal values of volScalarField N2
Adding all cells of cellZone zoneTwo ...
Found matching zone zoneTwo with 9216 cells.
Setting internal values of volScalarField O2
Setting internal values of volScalarField N2
Adding all cells of cellZone zoneNine ...
Found matching zone zoneNine with 31104 cells.
Setting internal values of volScalarField O2
Setting internal values of volScalarField N2
Adding all cells of cellZone zoneEight ...
Found matching zone zoneEight with 864 cells.
Setting internal values of volScalarField O2
Setting internal values of volScalarField N2
Adding all cells of cellZone zoneSeven ...
Found matching zone zoneSeven with 8208 cells.
Setting internal values of volScalarField O2
Setting internal values of volScalarField N2
Adding all cells of cellZone zoneSix ...
Found matching zone zoneSix with 2160 cells.
Setting internal values of volScalarField O2
Setting internal values of volScalarField N2

End

jmdag1981@jmdag1981-Gazelle-Professional:~/OpenFOAM/jmdag1981-2.3.0/run/EAPDR/EAPDR$ reactingFoam
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.3.0-f5222ca19ce6
Exec : reactingFoam
Date : Dec 15 2014
Time : 13:23:31
Host : "jmdag1981-Gazelle-Professional"
PID : 5053
Case : /home/jmdag1981/OpenFOAM/jmdag1981-2.3.0/run/EAPDR/EAPDR
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


Reading g
Creating reaction model

Selecting combustion model laminar<psiChemistryCombustion>
Selecting chemistry type
{
chemistrySolver ode;
chemistryThermo psi;
}

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

Selecting chemistryReader chemkinReader


--> FOAM FATAL ERROR:
Unknown chemistryReader type chemkinReader

Valid chemistryReader types are:
1(foamChemistryReader)

From function chemistryReader::New(const dictionary&, speciesTable&)
in file lnInclude/chemistryReader.C at line 55.

FOAM exiting

Not sure what I'm doing incorrectly. I've got a chem.inp and therm.dat file respectively in a chemkin folder. I've run through tutorials and searched on this forum but can't find anything regarding this error. Most likely it's an easy fix, I wouldn't be shocked if I am overlooking something very simple. Certainly possible I am being a complete bonehead here! Any help would be greatly appreciated!
JMDag2004 is offline   Reply With Quote

Old   December 15, 2014, 16:38
Default
  #2
New Member
 
Josh Dawson
Join Date: Oct 2014
Location: Justin, TX
Posts: 19
Rep Power: 11
JMDag2004 is on a distinguished road
Can anyone help me? I'm not sure if I am posting these questions correctly, maybe there is something else I am doing incorrectly...but I am not getting any response. Even if someone doesn't have the time to give a complete answser, at least a nudge in the right direction maybe? Thanks in advance!
JMDag2004 is offline   Reply With Quote

Old   March 8, 2016, 22:38
Default
  #3
New Member
 
Thien Xuan Dinh
Join Date: Jul 2014
Location: Japan
Posts: 8
Rep Power: 11
dinh is on a distinguished road
Hi
It is to late, but please check energy option
dinh is offline   Reply With Quote

Reply

Tags
chemistryreader, chemkinreader


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
how to find distribution field of reaction rates when CHEMKIN is enabled? EdgarasM FLUENT 1 February 18, 2019 09:41
creating PDF Table from imported chemkin files idomeneo FLUENT 1 November 2, 2016 05:23
Chemkin - Detonation shrikant96 Main CFD Forum 2 July 19, 2016 04:11
Chemkin Library (CKLIB) juzer_700 OpenFOAM Programming & Development 0 February 20, 2013 00:41
Question on CHEMKIN INTERFACE Rui Yang Phoenics 7 January 30, 2002 04:04


All times are GMT -4. The time now is 09:22.