NUMERICAL MODELING OF FLOW OVER A RIGID WAVY SURFACE BY LES

A flow field over a rigid wavy surface is numerically generated. The g eneral concept of the large eddy simulation that captures the large sc ale flow structure is utilized. A coordinate transformation scheme is developed to transform any wavy surface to a flat surface. The numeric al computations are carried out in the omega - psi domain. The subgrid scales are resolved by utilizing a Smagorinsky subgrid model. The fou rth-order central difference and the third-order upwinding schemes are used for the diffusive and convective terms, respectively. Adams-Bash forth and Dufort-Frankel schemes are used to solve the equations of th e fluid motion. The computational scheme ensured a non-penetration and a no-slip boundary condition. The velocity profile, and the separatio n and reattachment locations are found to agree well with an experimen tal investigation.