[Astrodan]

Dear everyone,

I'm currently trying to run simulations with reacting species. Species are transported using a simple passive scalar transport equation, reaction terms are introduced as source terms into the equation (own solver based on icoFoam).

The thing I haven't figured out is how to declare the outlet boundary condition for the species e.g. in a pipe flow. Assuming my reaction is finished halfway the pipe, zeroGradient would work. But if the reaction still takes place at the outlet, I would need something like a reactionSpeedGradient BC, which sets the gradient based on the local reaction rate (which does exist as a volScalarField in the solver).

Is there something similar, or would I have to write my own BC? And if so, what base BC could you suggest me to use?

Thanks,
-Timm