We present a stabilized numerical formulation for incompressible continua based on a higher‐order Finite Calculus (FIC) approach and the finite element method. The focus of the paper is on the derivation of a stabilized form for the mass balance (incompressibility) equation. The simpler form of the momentum equations neglecting the non‐linear convective terms, which is typical for incompressible solids, Stokes flows and Lagrangian flows is used for the sake of clarity. The discretized stabilized mass balance equation adds to the standard divergence of velocity term a pressure Laplacian and an additional boundary term. The boundary term is relevant for the accuracy of the numerical solution, especially for free surface flow problems. The Laplacian and boundary stabilization terms are multiplied by non‐linear parameters that have an extremely simple expression in terms of element sizes, the pressure and the discrete residuals of the incompressibility equation and the momentum equations, thus ensuring the consistency of the method. The stabilized formulation allows solving the incompressible problem iteratively using an equal‐order interpolation for the velocities (or displacements) and the pressure, which are the only unknowns. The use of additional pressure gradient projection variables, typical of many stabilized methods, is unnecessary.

We present a stabilized numerical formulation for incompressible continua based on a higher‐order Finite Calculus (FIC) approach and the finite element [...]