Adding Radiation to a turbulent solver
 

Hello all.
I hope you are absolutely fine.

I have a solver that solves (air flow)turbulent mixed-convection in a vertical parallel plates channel. I want to add radiation to my solver. How can I do it?! what are the best ways to add radiation to my solver? my solver is steady, incompressible and with the use of bousinessq approximation for density.

I use openFoam 2.1.1. I heard that in new versions, radiation is included in all the buoyant solvers, is it right?

My other question is that is it right to use radiation along with bousinessq approximation for density?

Thank you and best regards.
Masoud