Continuum equations governing thermal non-equilibrium modeling of steady natural convection inside wavy enclosures with the effect of thermal radiation are developed. These equations account for such effects as the inter-phase heat transfer coefficient effect, the thermal radiation effect, the modified conductivity ratio effect and the Rayleigh number effect. Finite difference method is employed to solve these equations and comparisons between previous published works are presented. Numerical results for the flow and heat transfer for the fluid and solid phases are obtained for various combinations of the physical parameters. Graphical and tabular results illustrating interesting features of the physics of the problem are presented and discussed.