In the adaptive mesh generation, the space mesh should be adequate to the surface mesh. When the analytical surface representation is not known, additional surface information may be extracted from triangular surface meshes. We describe a new surface reconstruction method which uses approximate Hessian of a piecewise linear function representing the discrete surface. Efficiency of the proposed method is illustrated with two CFD applications.