A FAST NUMERICAL METHOD FOR A SIMPLIFIED PHASE FIELD MODEL

• Published in: Mathematical Methods In Scattering Theory And Biomedical Engineering pp. 208 - 215
• Main Authors: SFYRAKIS, CHR. A., DOUGALIS, V. A.
• Format: Book Chapter
• Language: English
• Published: WORLD SCIENTIFIC 01.08.2006
• Subjects: Applied Mathematics
• ISBN: 9812568603; 9812773193; 9789812773197; 9789814477598; 9814477591; 9789812568601
• DOI: 10.1142/9789812773197_0022

We consider the phase field model consisting of the system of p.d.e' s \begin{array}{l} q(\theta)\phi t = \nabla \cdot (A\nabla\phi) + f(\phi, u)\,,\\[2pt] u_t = \Delta u + [p(\phi)]_t\,, \end{array} where ϕ = ϕ(x,y,t) is the phase indicator function, θ = arctan $\frac{\phi_{y}}{\phi_{x}}$, u = u(x,y,t) is the temperature, q, p, and f are given scalar functions, and A is a 2 × 2 matrix of suitable functions of θ. This system describes the evolution of phase and temperature in a two phase medium, and is posed for t ≥ 0 on a rectangle in the x, y space with appropriate boundary and initial conditions. We solve the system numerically by a finite difference ADI method and show the results of relevant numerical experiments in the semi-anisotropic case, where A = I. In the case of an isotropic medium with q = 1, and A a diagonal matrix of functions of ϕ, we prove an error estimate of second order accuracy in space and time.