On hp-adaptive solution of thin shells of revolution

We present an a posteriori hp-adaptive algorithm for shells. In an hp-solver there are two ways to increase accuracy: one can either refine the mesh (h-step) or increase the degree of the polynomial (p-step). In the adaptive setting it is also necessary to have the capability to reverse already made decisions. Our p-approach is influenced by work of Houston and Süli on the estimation of Sobolev regularity and analycity. Due to the possibility of numerical locking we initialize our algorithm only after initial probing of the energy distribution of the solution.