Numerical Tests of a Phase Field Model with Second Order Accuracy

Numerical computations are performed for a recently derived phase field model for the interface between two phases. The rigorous results indicate that solutions to this new phase field model should converge more rapidly than traditional ones to solutions of the corresponding sharp interface (free boundary) formulation for sufficiently small values of the approximation parameter e representing the thickness of the interfacial region. In particular, the distance between the sharp interface of the limiting model and the zero level set of the phase function in the phase field model is of order e 2 rather than e. Numerical computations within a three-dimensional spherically symmetric setting compare the computed solutions of this new model with the known exact solutions for the limiting free boundary problem and confirm the second order accuracy predictions of the theory for sufficiently small e. The sets of parameters include those of succinonitrile used in dendritic experiments.

[1]  W. Carter,et al.  Vector-valued phase field model for crystallization and grain boundary formation , 1998 .

[2]  Fernando Reitich,et al.  Local existence and uniqueness of solutions of the Stefan Problem with surface tension and kinetic undercooling , 1992 .

[3]  Gunduz Caginalp,et al.  Phase field equations in the singular limit of sharp interface problems , 1992 .

[4]  Department of Physics,et al.  Universal dynamics of phase-field models for dendritic growth , 1999 .

[5]  G. Caginalp,et al.  RAPIDLY CONVERGING PHASE FIELD MODELS VIA SECOND ORDER ASYMPTOTICS , 2005 .

[6]  G. W. Young,et al.  Comparison of Asymptotic Solutions of a Phase-Field Model to a Sharp-Interface Model , 2001, SIAM J. Appl. Math..

[7]  ON THE CLASSICAL SOLUTION OF THE MULTIDIMENSIONAL STEFAN PROBLEM FOR QUASILINEAR PARABOLIC EQUATIONS , 1981 .

[8]  Gunduz Caginalp,et al.  The role of microscopic anisotropy in the macroscopic behavior of a phase boundary , 1986 .

[9]  Halil Mete Soner,et al.  Convergence of the phase-field equations to the mullins-sekerka problem with kinetic undercooling , 1995 .

[10]  Ricardo H. Nochetto,et al.  Convergence Past Singularities for a Fully Discrete Approximation of Curvature-Drive Interfaces , 1997 .

[11]  M. Chipot.  Linear Parabolic Equations , 2000 .

[12]  A. Karma,et al.  Quantitative phase-field modeling of dendritic growth in two and three dimensions , 1996 .

[13]  Robert Almgren,et al.  Second-Order Phase Field Asymptotics for Unequal Conductivities , 1999, SIAM J. Appl. Math..

[14]  G. Caginalp,et al.  Efficient computation of a sharp interface by spreading via phase field methods , 1989 .

[15]  L. I. Rubinshteĭn.  The Stefan Problem , 1971 .

[16]  G. Caginalp,et al.  A rapidly converging phase field model , 2006 .

[17]  M. Glicksman,et al.  Dendritic Growth tip velocities and radii of curvature in microgravity , 1999 .

[18]  Gunduz Caginalp,et al.  Convergence of the phase field model to its sharp interface limits , 1998, European Journal of Applied Mathematics.

[19]  Peter W. Bates,et al.  Convergence of the Cahn-Hilliard equation to the Hele-Shaw model , 1994 .

[20]  Jean-Baptiste Lully,et al.  The collected works , 1996 .

[21]  G. Caginalp,et al.  Computation of sharp phase boundaries by spreading: the planar and spherically symmetric cases , 1991 .

[22]  G. Caginalp.  An analysis of a phase field model of a free boundary , 1986 .

[23]  R. Kobayashi.  Modeling and numerical simulations of dendritic crystal growth , 1993 .

[24]  Martin E. Glicksman,et al.  Dendritic growth-A test of theory , 1976 .

[25]  William Mendenhall,et al.  Introduction to Probability and Statistics , 1968 .

[26]  A. A. Wheeler,et al.  Thermodynamically-consistent phase-field models for solidification , 1992 .

[27]  A. M. Meirmanov,et al.  The Stefan Problem , 1992 .

[28]  Harald Garcke,et al.  Second order phase field asymptotics for multi-component systems , 2006 .

[29]  Initial conditions implied by t exp 1/2 solidification of a sphere with capillarity and interfacial kinetics , 1988 .

[30]  J. Cahn,et al.  A microscopic theory for antiphase boundary motion and its application to antiphase domain coasening , 1979 .

[31]  Gunduz Caginalp,et al.  Computations of Dendrites in 3-D and Comparison with Microgravity Experiments , 2003 .

[32]  R. Sekerka,et al.  Stability of a Planar Interface During Solidification of a Dilute Binary Alloy , 1964 .