This paper proposes a vertex-estimation-based, feature-preserving smoothing technique for meshes. A robust mesh smoothing operator called mean value coordinates flow is introduced to modify mean curvature flow and make it more stable. Also the paper proposes a three-pass vertex estimation based on bilateral filtering of local neighbors which is transferred from image processing settings and a Quasi-Laplacian operation, derived from the standard Laplacian operator, is performed to increase the smoothness order of the mesh rapidly whilst denoising meshes efficiently, preventing volume shrinkage as well as preserving sharp features of the mesh. Compared with previous algorithms, the result shows it is simple, efficient and robust.