Abstract: In this paper we present a novel model for computing the oriented normal field on a point cloud. Differently from previous two-stage approaches, our method integrates the unoriented normal estimation and the consistent normal orientation into one variational framework. The normal field with consistent orientation is obtained by minimizing a combination of the Dirichlet energy and the coupled-orthogonality deviation, which controls the normals perpendicular to and continuously varying on the underlying shape. The variational model leads to solving an eigenvalue problem. If unoriented normal field is provided, the model can be modified for consistent normal orientation. We also present experiments which demonstrate that our estimates of oriented normal vectors are accurate for smooth point clouds, and robust in the presence of noise, and reliable for surfaces with sharp features, e.g., corners, ridges, close-by sheets and thin structures.