reactionsSensitivityAnalysis does not work in OpenFOAM-v1806
 

Dear foamers,



I installed the version 1806 of openfoam by following step by step the following link:

https://www.openfoam.com/download/install-source.php


I did that as there is a new post processing tool in this version which is called reactionsSensitivityAnalysis:
(please have a look at https://www.openfoam.com/releases/op...processing.php)


And there is a tutorial case in /OpenFOAM/OpenFOAM-v1806/tutorials/combustion/reactingFoam/laminar/counterFlowFlame2D/ which includes this functionObjetc in the controlDict:



functions
{
sensitivityAnalysis
{
functionObjectLibs ("libutilityFunctionObjects.so");
type psiReactionsSensitivityAnalysis;
outputControl outputTime;
enabled true;
}
}


However when I run the case, I have the following error :





From function static Foam::autoPtr<Foam::functionObject> Foam::functionObject::New(const Foam::word&, const Foam::Time&, const Foam::dictionary&)
105 in file db/functionObjects/functionObject/functionObject.C at line 108.
106 From function bool Foam::functionObjectList::read()
107 in file db/functionObjects/functionObjectList/functionObjectList.C at line 832
108
109 --> while loading function object 'sensitivityAnalysis'
110
111
112 Starting time loop
113
114 --> FOAM Warning :
115 Unknown function type psiReactionsSensitivityAnalysis
116
117 Valid function types :
118
119 22
120 (
121 abort
122 coded
123 ensightWrite
124 patchProbes
125 probes
126 psiReactionThermoMoleFractions
127 psiSpecieReactionRates
128 removeRegisteredObject
129 residuals
130 rhoReactionThermoMoleFractions
131 rhoSpecieReactionRates
132 runTimeControl
133 setTimeStep
134 sets
135 surfMeshes
136 surfaces
137 systemCall
138 thermoCoupleProbes
139 timeActivatedFileUpdate
140 vtkWrite
141 writeDictionary
142 writeObjects
143 )









As I can see in /home/maryam/OpenFOAM/OpenFOAM-v1806/src/functionObjects/field/Make/file, the library name is libfieldFunctionObjects and not libutilityFunctionObjects. I change that and I have :




--> FOAM Warning :
72 Function object only applicable to single cell cases
73
74 From function Foam::functionObjects::reactionsSensitivityAnalysi s<chemistryType>::reactionsSensitivityAnalysis(con st Foam::word&, const Foam::Time&, const Foam::dictionary&) [with chemistryType = Foam: :BasicChemistryModel<Foam::psiReactionThermo>]
75 in file reactionSensitivityAnalysis/reactionsSensitivityAnalysis.C at line 204.
76 From function bool Foam::functionObjectList::read()
77 in file db/functionObjects/functionObjectList/functionObjectList.C at line 832




Anyone has experience with this reaction sensitivity analysis tool ? Is it working ? It is so strange that the tutorial case does not work.



Anyone can help me with that ? Am I doing a mistake or is there a bug ?



Thanks in advance

No one can help me with this post process tool ?

I tried the tutorial case of openFoam, the gri in solver/combustion/chemFoam. If you put :


functions
{
sensitivityAnalysis
{
libs ("libfieldFunctionObjects.so");
type psiReactionsSensitivityAnalysis;
writeControl writeTime;
}
}


in controlDict, it works fine during the run time. But if you remove the postProcessing folder generated during the run time and you try to generate the same folder after the simulation is finished by running chemFoam -postProcess -func psiReactionsSensitivityAnalysis ,
you have :





--> FOAM Warning :
From function static bool Foam::functionObjectList::readFunctionObject(const Foam::string&, Foam::dictionary&, Foam::HashSet<Foam::wordRe>&, const Foam::word&)
in file db/functionObjects/functionObjectList/functionObjectList.C at line 276
Cannot find functionObject file psiReactionsSensitivityAnalysis
Time = 0


--> FOAM FATAL ERROR:
Cannot find file "points" in directory "polyMesh" in times "0" down to constant

From function virtual Foam::IOobject Foam::fileOperation::findInstance(const Foam::IOobject&, Foam::scalar, const Foam::word&) const
in file global/fileOperations/fileOperation/fileOperation.C at line 876.

FOAM exiting



Is it normal ? Isn't it a bug ?

I have also met this bug on the post-processing psiSpecieReactionRates, which is used to calculate reaction rate of each species in each reaction pathway. When testing the reactingFoam tutorial, counterFlowFlame2D, It works well in run-time postprocessing. While I remove the generated postprocessing file, and use the command：reactingFoam -postProcess, it reports:
Create time

Create mesh for time = 0


PIMPLE: Operating solver in PISO mode

--> FOAM Warning :
Function object only applicable to single cell cases

From function Foam::functionObjects::reactionsSensitivityAnalysi s<chemistryType>::reactionsSensitivityAnalysis(con st Foam::word&, const Foam::Time&, const Foam::dictionary&) [with chemistryType = Foam::psiChemistryModel]
in file reactionSensitivityAnalysis/reactionsSensitivityAnalysis.C at line 204.
From function bool Foam::functionObjectList::read()
in file db/functionObjects/functionObjectList/functionObjectList.C at line 840


--> FOAM Warning :

request for psi chemistryProperties from objectRegistry region0 failed
available objects of type psi are
0()

From function const Type& Foam::objectRegistry::lookupObject(const Foam::word&, bool) const [with Type = Foam::psiChemistryModel]
in file /gpfs01/home/xingjiangkuan/OpenFOAM/OpenFOAM-v1712/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 239.
From function bool Foam::functionObjectList::read()




I am also confused with this problem. Anyone has some suggestions on solving this issue. We can discuss and solve this bug together.