Second and fourth-order adaptive phase field modeling of fracture using PHT-splines in the framework of IGA.