Numerical approximation for the phase-field transition system

The 1D phase-field transition system introduced by Caginalp to describe the moving boundary in melting problems is considered. It is discretized by finite differences and three algorithms are presented to solve the resulting nonlinear algebraic system: the Newton method, an improved Newton method with reduced system and a fractional step method. Numerical results are reported and a comparison between the algorithms is made.