A system and method for a hybrid, variational, user-controlled, 3D mesh smoothing for orphaned shell meshes. The smoothing model is based on a variational combination of energy and equi-potential minimization theories. A variety of smoothing techniques for predicting a new location for the node-to-smooth are employed. Each node is moved according to a specific smoothing algorithm so as to keep element included angles, skew and distortion to a minimum. The variational smoother selection logic is based on nodal valency and element connectivity pattern of the node to smooth. Results show its consistency with both quadrilateral and quad-dominant meshes with a significant gain over conventional Laplacian schemes in terms of mesh quality, stability, user control and flexibility.