High accuracy solutions to energy gradient flows from material science models