[KeiJun]

Dear Sir,

I am trying to calculate a system including chemical reaction using reactingFoam solver, but I don't understand how to use this solver.
Could you teach me the following questions?

1). Is it possible to describe a reaction rate equation in the constant/reactions file by myself, not choosing the "type"?
And moreover, I want to describe parameters of the rate equation too.
(Parameters are the function of pressure.)

2). Are "reaction->R(Yi) in YEqn.H calculated from the reaction rate equation? If so, can I use the result calculated by my described equation (above 1).) as R(Yi)?

3). Can I also use the result calculated by my described equation (above 1).) as "reaction->Sh()" in EEqn.H?

4). Calculation I want to do is not combustion reaction, but is the "constant/combustionProperties" file needed for the calculation?

5). Where is alpha (=k/cp) calculated? I think α must be calculated to calculated "αeff".


Thanks in advance.

KeiJun

[KeiJun]

Dear Sir,

I keep trying to customize reactingFoam solver, and what I don't understand how to describe remains.

I want to describe the reaction rate equation to UEqn.H, not refering the "type" on constant/reactions file.

The reaction rate equation is
∂qi/∂t = ki(qi* - qi),
where qi is the reaction amount of the specie i per catalyst (the unit of qi is mol/kg), and qi* is the equilibrium reaction amount per catalyst.

I can compile myReactingFoam.C, in which I add the above rate equation.

However, when I run the tutorials, the error of reaction type has come.

Then I want to make the reaction rate file mensioned above in src/thermophysicalModels/specie/reaction/reactionRate(?), not to emerge the error.
Alternatively, I want not to refer Reaction.C during running a tutorials, not to refer the type on constant/reactions.

Could you teach me how to do?

I am having a very hard time now.

I am really looking forward to your reply.


Sincerely,

KeiJun