Abstract

The Navier–Stokes equations written in Laplace form are often the starting point of many numerical methods for the simulation of viscous flows. Imposing the natural boundary conditions of the Laplace form or neglecting the viscous contributions on free surfaces are traditionally considered reasonable and harmless assumptions. With these boundary conditions any formulation derived from integral methods (like finite elements or finite volumes) recovers the pure Laplacian aspect of the strong form of the equations. This approach has also the advantage of being convenient in terms of computational effort and, as a consequence, it is used extensively. However, we have recently discovered that these resulting Laplacian formulations violate a basic axiom of continuum mechanics: the principle of objectivity. In the present article we give an accurate account about these topics. We also show that unexpected differences may sometimes arise between Laplace discretizations and divergence discretizations.