Hexahedral Mesh Refinement Using an Error Sizing Function