A primal formulation for the Helmholtz decomposition

In 1999, Jean-Paul Caltagirone and Jerome Breil have developed in their paper [Caltagirone, J. Breil, Sur une methode de projection vectorielle pour la resolution des equations de Navier-Stokes, C.R. Acad. Sci. Paris 327(Serie II b) (1999) 1179-1184] a new method to compute a divergence-free velocity. They have used the grad(div) operator to extract the solenoidal part of a given vector field. In this contribution we explain how this method can be considered as a real Helmholtz decomposition and we present a stable approximation in the framework of spectral methods. Numerical results are presented to illustrate the efficiency of this approach. (c) 2007 Elsevier Inc. All rights reserved.