Accessing reaction information from PaSR.C
 

Hello!

I want to modify the standard PaSR to calculate reaction rates based solely on turbulent mixing, without using Arrhenius data (Eddy-dissipation model). The current PaSR+ODEChemistryModel has no "direct dependency" on turbulent mixing time scale, in that the turbulent time scale appears only in the reactive volume scaling factor, kappa. The actual reaction rates are calculated based on Arrhenius data in ODEChemistryModel.

I figured out that probably the easiest way to proceed would be to modify the PaSR files, by copy-pasting the original PaSR folder to and by giving it a different name.

I would like to have an access to the reaction data (list of reactions, stoichiometric coefficients, etc.), as is done in the standard ODEChemistryModel.

What is the easiest way to access this reaction data?

I thought I could simply use a pointer to the ODEChemistryModel (this->pChemistry_ in PaSR). ODEChemistryModel.H has the following line in the list of public member functions:
Therefore, I thought, I could get the list of reactions simply with the following line:
However, as a result the compiler prints the following errors:
Any ideas how to access the reaction information from PaSR?

Thanks,

Tatu

I know it has been a while but did you figure this out and do you remember how?

Which version of OF are you using?
And you want to access `reactions` in the solver or in the library?

I am using OpenFOAM v5.0 and I would like to access reactions at the solver level. I am trying to modify the reactingFoam solver. Any insights will be helpful?

I think it is difficult to access reactions in solver.
Because: chemistryPtr_ in ChemistryCombustion class is a Protected data, means that it can be only used in itself or its derived class, but not its object.
I plotted a UML for this issue as in my bolg: https://openfoam.top/chemistry/. Sorry for the language is Chinese and the version is OF-dev, but the UML figure may be a little help for you.
UML

Yes, I agree. The UML figure is quite helpful. Thank you.

I would like to elaborate a bit on my problem and would really like your input on that.

1. Currently,the reaction rate is evaluated in "omega" function of the chemistryModel.C file. This reaction rate uses the concentration of the reactants and products as well as the reaction rate (k) calculated from the reaction rate models.

2. Instead, I need to calculate the reaction rate using turbulence parameters. The issue is that turbulence parameters cannot be accessed by the chemistryPtr_.

3. I have tried passing turbulence parameters as function arguments but they somehow end up giving an error due to some templates that the chemistryModel is using.

Any help on how to call the eddy dissipation energy, and kinematic viscosity fields in the chemistryModel will be really helpful.