]>

Pseudo Code: (1) Formulate and initialize STLSFEM (2) Generate bubble fns using Jacobi polynomials (3) Introduce bubble fns into original set of nodal shape function using tensor product, element stiffness matrix size goes up from original $m$ to $p=m+b^n$. Where $n$ is number of dimensions. (4) While $(p\ge m)$$\{{/}^{\ast }$to get the original size of element stiffness matrix back${}^{\ast }/$ Apply Static Condensation() $p=p-1$; $\}$ (6) Set the time limit and convect the solution using linear solver (7) end

Algorithm 1