“…The aforementioned codes resort either on finite differences/volumes method (e.g., GERMINAL [15] and TRANSURANUS [8]) or on the finite element method (e.g., CEDAR [13] and BISON [14,17]) to solve the equations governing the migration of porosity. The solution of the (pure) advection equation by the standard Galerkin Finite Element Method (G-FEM), resorting on a centered scheme for the discretization of gradient operator, is known to generate spurious spatial oscillations (see for example [20]).…”