We set up simulation principles,and build simulation method and codes with laminar kinetic energy turbulent model,in which natural,bypass and separation translation are considered.Time derivatives of three-dimensional Navier-Stokes equations are preconditioned with Weiss-Smith preconditioning matrix.Effect of turbulent eddy viscosity is considered.Time split,spatial split and boundary conditions are specially treated to diagonal dominant matrix.Boundary layer and airfoil flow simulations are available.It is shown that the method is efficient and robust for basic aerodynamic flows.It provides reference to more complex flows.