A method for adaptive 3D mesh generation based on approximate Hessian recovery is proposed and analyzed. Several numerical results illustrate basic properties of the method.