[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,

[Severus]

Hello,
You can use fvc instead.
For example the continuity error is calculated in a similar fashion

 volScalarField contErr(fvc::div(phi));

See line 33: https://cpp.openfoam.org/v9/src_2fin...8H_source.html

Thanks

[Sereff]

Thank you Severus, and I guess if I use

fvc::div(phi,Yi)

then it would give me the convection of species Yi for all cells?

kind regards,

[Severus]

Hello,
Yes you should get whole field values (i.e. for all cells).

Thanks