Adaptive phase-field modeling of three-dimensional thermally-induced brittle fracture in COMSOL