Symmetrization of matrix technique in a finite element context

A finite element Galerkin-based formulation of the mass conservation and momentum equations can require, if convective type terms are retained in the coefficient matrix, a non-symmetric solver. The resulting increase in core storage for efficient utilization of CPU time can be considerable. The current paper advocates a simple symmetrization of matrix technique, at element level which results in a considerable reduction in core requirement. The increase in CPU time required when solving linear systems of equations is considerable. However, for nonlinear systems the penalty can be negligible.