Combined heat and mass transfer process by natural convection from a wavy vertical surface immersed in a fluid-saturated semi-infinite porous medium due to Soret and Dufour effects for Forchheimer extended non-Darcy model has been analyzed. A similarity transformation followed by a wavy to flat surface transformation is applied to the governing coupled non-linear partial differential equations, and they are reduced to boundary layer equations. The obtained boundary layer equations are solved by finite difference scheme based on the Keller-Box approach in conjunction with block-tridiagonal solver. Detailed simulations are carried out for a wide range of parameters like Groshof number (Gr*), Lewis number (Le), Buoyancy ratio (B), Wavy wall amplitude (a), Soret number (S (r) ), and Dufour number (D (f) ). Comparison tables local and average Nusselt (Nu) number, local and average Sherwood (Sh) number plots are presented.