Abstract

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.

The formulation is particularly useful for heterogeneous incompressible materials with discontinuous material properties, as it allows computing all the stabilization matrices at the element level. Details of the finite element formulation are given. The good behaviour of the new pressure Laplacian stabilization (PLS) technique is shown in simple but demonstrative examples of application. A very accurate solution was obtained in all cases in 2–3 iterations. Copyright © 2010 John Wiley & Sons, Ltd.

Abstract

We present a collection of stabilized finite element (FE) methods derived via first‐ and second‐order finite calculus (FIC) procedures. It is shown that several well known existing stabilized FE methods such as the penalty technique, the Galerkin Least Square (GLS) method, the Pressure Gradient Projection (PGP) method and the orthogonal sub‐scales (OSS) method are recovered from the general residual‐based FIC stabilized form. A new family of stabilized Pressure Laplacian Stabilization (PLS) FE methods with consistent nonlinear forms of the stabilization parameters are derived. The distinct feature of the family of PLS methods is that they are residual‐based, i.e. the stabilization terms depend on the discrete residuals of the momentum and/or the incompressibility equations. The advantages and disadvantages of the different stabilization techniques are discussed and several examples of application are presented.