Hierarchical hybrid grids: achieving TERAFLOP performance on large scale finite element simulations