|May 20, 2016, 05:04||
How are chemical reactions stored after parsing?
Join Date: May 2016
Posts: 2Rep Power: 0
I'm currently trying to implement an enhancement for the modifiedReactingFoam-Solver to allow it to properly calculate surface reactions defined by the user.
I want the solver to allow for surface reactions to be defined and "attached" to any boundary within the domain and I want them to be proper chemical reactions that can for example create surface adsorbed species that stick to the boundaries or to generate reaction products which go right back into the adjecent fluid domain.
In order to do this, I need to solve some equations to get the reaction rates right and that's where I'm currently running into trouble.
Does somebody, by any chance, know how the chemical reactions and species are stored after parsing? I'm less interested in the storage of arrhenius parameters for now. I rather need to know how and where a chemical reaction and it's stoichometric coefficients are stored, so the solver later knows which of the species actually react within a given reaction.
Let's say I have five species:
A, B, C, D, E
and two reactions:
1) A + B = C
2) A + 2*D = B + E
The best case scenario for me would be if, upon parsing, the reactions from the reactions dictionary would generate a structure looking like this:
1 1 0 0 0 1 0 0 2 0
0 0 1 0 0 0 1 0 0 1
The key feature I am asking for here is some structure that stores the chemical reactions in a way that has each reaction contain stoichometric coefficients for all species whether they are actually affected by the reaction or not. If a specie doesn't occur in a certain reaction it simply gets stoichometric cofficient 0 as would be the case for specie D in the first reaction in my example above.The first reaction rewritten with what I have in mind would look something like this:
1*A + 1*B + 0*C + 0*D + 0*E = 0*A + 0*B + 1*C + 0*D + 0*E
So if anybody knows whether reactions are already stored like that and where I can find them, that would be really helpful. Otherwise I will implement some code that arranges them in a scheme similar to that myself, but I don't want to implement something that's already there just because I can't find it
Many thanks in advance. Even if somebody can tell me that this is currently not the case, but can tell me how chemical reactions are stored instead, that would be much appreciated.
Thanks to anyone reading the whole post,
|May 20, 2016, 11:27||
Join Date: Apr 2009
Location: Winnipeg, Canada
Posts: 436Rep Power: 14
Try looking into the ODE reaction solvers. They need to figure that out. I'm thinking chemistrySolver, chemistryModel. Look for where they implement derivatives and jacobian. Just off the top of my head.
Follow me on twitter @DavidGaden
|May 20, 2016, 12:05||
Join Date: May 2016
Posts: 2Rep Power: 0
I did examine the chemistryModel and EulerImplicit, which is used to build and solve the differential equation from the chemical reactions.
The algorithm in place iterates through all the reactions and sorts the coefficients "by hand" on each iteration to build the jacobian however.
From what I saw in the code and found through testing (Info << is my friend) it seems that the reactions are stored in a list and each element of this list has a component "rhs" and "lhs" which in itself is a list (containing the species that are present in this reactions' rhs/lhs), which contains yet another list for each specie appearing on the left/right hand side, which then contains the specie-ID and stoichometric coefficient.
So while I was able to rebuild each chemical reaction from those elements, it was not in a format I'd prefer. The problem I have with this is that I first need to determine whether a reaction affects a certain specie or not, if I want to calculate anything like net consumption rates for individual species.
I decided to build my own function for now, which creates the matrices I described above before any calculations are done at all, so I have to deal with this sorting process only once. So for now I have a solution and can continue my work and if I don't stumble upon something already implemented in OF, I'll just leave it in
|Thread||Thread Starter||Forum||Replies||Last Post|
|multi phase flow with chemical reactions in one phase||Habibfateh||OpenFOAM Programming & Development||0||February 10, 2014 08:33|
|Chemical reactions in openFOAM||nakul||OpenFOAM||2||April 22, 2013 11:51|
|inconsistency in the chemical reactions in Reacting Foam?||lfgmarc||OpenFOAM Running, Solving & CFD||1||October 10, 2011 11:11|
|How to specify chemical reactions in certain areas||Esti||CD-adapco||0||September 24, 2004 06:03|
|Finite rate chemical reactions||Marc Segovia||CD-adapco||3||November 28, 2001 05:21|