[alicharchi]

I am working on reactive flow and I do not understand the way reactingFOAM solves the Yi equations.

To be more specific, the question is how the source terms (reaction->R(Yi)) are discretized? I am quite sure that the Ode solver is used to calculate that term in reaction->correct() member function while the temporal term (fvm::ddt(rho, Yi)) is discretized by Euler ddt scheme.

I have noticed the R member function returns a pointer of type fvScalarMatrix.
Can anyone give some insight into the procedure of discretizing the source term?

Note that I am using OF 2.3.1.

Thanks

[kbeat]

step 1 The reaction-> correct() actually solves the stiff chemical ode to get R[i]
(order of integration is 10e-9 integrated upto ~ deltaT)

step 2 the species equation is solved based normal implicit scheme

Since the implementation allows interchangeable use of fvScalarMatric and volScalarField( difference being the boundary)