[eretief]

Hi

I am trying to develop a custom boundary condition in OpenFOAM. The boundary condition value is calculated from a polynomial expression written in terms of the heat flux at the boundary. I am using the reactingFoam solver.

Can anyone perhaps tell me how exactly I could access the necessary variables (e.g. alphaEff) at the boundary to calculate the heat flux at the boundary?

Thanks

[chriss85]

Likely by using the thermo class, using something like this:

const fluidThermo& thermo_(mesh.lookupObject<fluidThermo>("thermophysicalProperties"))