[Sereff]

Dear Foamers,

I am working on a solver for a chemical reaction.

fvScalarMatrix YiEqn
(
    fvm::ddt(Yi)
  + fvm::div(phi, Yi)
  - fvm::laplacian(Deff, Yi)
 ==
    fvOptions(Yi)
  + reaction->R(i)
);

And for one of the reactions, in order to calculate the source term R(i), I would need the information of the convection of the species of each cell with its surrounding cells and boundaries. I am aware that there's

fvm::div(phi, Yi)

but I need the information as a volScalarField or GeometricField<scalar>. And the adorementioned term is operating on a fvScalarMatrix, which contains a lot more information and I am not sure how to extract the information i want. Could someone help me on this issue?

kind regards,