A new method for direct numerical simulation of incompressible Navier-Stokes equations is studied in the paper. The compact finite difference and the non-linear terms upwind compact finite difference schemes on non-uniform meshes in x and y directions are developed respectively. With the Fourier spectral expansion in the spanwise direction, three-dimensional N-S equation are converted to a system of two-dimensional explicit-implicit The treatment of equations. The third-order mixed scheme is employed the three-dimensional for time integration. non-reflecting outflow boundary conditions is presented, which is important for the numerical simulations of the problem of transition in boundary layers, jets, and mixing layer. The numerical results indicate that high accuracy, stabilization and efficiency are achieved by the proposed numerical method. In addition, a theory model for the coherent structure in a laminar boundary layer is also proposed, based on which the numerical method is implemented to the non-linear evolution of coherent structure. It is found that the numerical results of the distribution of Reynolds stress, the formation of high shear layer, and the event of ejection and sweeping, match well with the observed characteristics of the coherent structures in a turbulence boundary layer.
The numerical solution of the stable basic flow on a 3-D boundary layer is obtained by using local ejection, local suction, and combination of local ejection and suction to simulate the local rough wall. The evolution of 3-D disturbance T-S wave is studied in the spatial processes, and the effects of form and distribution structure of local roughness on the growth rate of the 3-D disturbance wave and the flow stability are discussed. Numerical results show that the growth of the disturbance wave and the form of vortices are accelerated by the 3-D local roughness. The modification of basic flow owing to the evolvement of the finite amplitude disturbance wave and the existence of spanwise velocity induced by the 3-D local roughness affects the stability of boundary layer. Propagation direction and phase of the disturbance wave shift obviously for the 3-D local roughness of the wall. The flow stability characteristics change if the form of the 2-D local roughness varies.
