Phase field method