The 3-D modeling of ocean finite element method(OFEM) using $k-{varepsilon}$ turbulent model and tetrahedron grids has been used to investigate the internal wave generation during the expansion of the deep water from the open sea to the shelf with a simple shape, which can be widely used in the fields of submarine development, ocean environment and meteorology, etc. In this paper, the detailed configuration of internal wave with its length and height and also the distribution of salinity and turbulent kinematic energy, etc. were derived. It is hoped that this OFEM method can be successfully applied to the numerical calculation of internal wave for and the oceanographic problems (tidal flows around underwater hill, plateau, Georges Bank, etc.) and ocean engineering problems(flow past artificial sea reefs) in future.