Reply To: Reactive ADE for gas flow
for Neumann boundary conditions that are not zero you could use the setInterpolatedConvectionBoundary or the local version setLocalInterpolatedConvectionBoundary that are made for Navier Stokes equations. You have to check if this also fits for your ADE application.
If you want to use a source-termed ADE use a field SOURCE and then calculate in a Post Processor (you can also include it in your coupling Post Processor) the source term with the computeRho method and set for each voxel your source term. The SOURCE field is then used by the SourcedAdvectionDiffusionBGKdynamics.