An efficient maximum bound principle preserving p-adaptive operator-splitting method for three-dimensional phase field shape transformation model. (15th August 2022)