Finite-element methods for analysis of the dynamics and control of Czochralski crystal growth

Numerical methods are presented for solution of the complex moving-boundary problem described by a thermal-capillary model for Czochralski crystal growth, which accounts for conduction through melt, crystal, and crucible and radiation between diffuse-gray body surfaces. Transients are included that are caused by energy transport, by changes in the shapes of the melt-crystal, melt-ambient phase boundaries and the moving crystal, and by the batchwise decrease of the melt volume in the crucible. Finite-element discretizations are used to approximate the moving boundaries and the energy equation in each phase. A two-level, implicit integration algorithm is presented for transient calculations. The temperature fields and moving boundaries are advanced in time by a trapezoid rule approximation with modified Newton's iterations to solve algebraic systems for effective ambient temperatures computed with diffuse-gray radiation. The implicit coupling between radiative exchange, interface shapes, and the temperature field is necessary for preserving the second-order accuracy of the integration method and is achieved by successive iterations between the radiation calculation and solution of the thermal capillary model. Analysis of a quasi-steady-state model (QSSM) demonstrates the inherent stability of the CZ process. Including either diffuse-gray radiation among crystal, melt, and crucible or a simple controller for maintaining constant radius can lead to oscillations in the crystal radius. The effects of these oscillations on batchwise crystal growth are addressed.

[1]  L. Petzold Differential/Algebraic Equations are not ODE's , 1982 .

[2]  H. Ettouney,et al.  A thermal-capillary mechanism for a growth rate limit in edge-defined film-fed growth of silicon sheets , 1986 .

[3]  G. C. Joyce,et al.  The weighing method of automatic Czochralski crystal growth: I. Basic theory , 1977 .

[4]  R. Sani,et al.  On the time-dependent solution of the incompressible Navier-Stokes equations in two and three dimensions , 1980 .

[5]  Lyle H. Ungar,et al.  Finite element methods for unsteady solidification problems arising in prediction of morphological structure , 1988, J. Sci. Comput..

[6]  A. B. Shapiro,et al.  FACET: a radiation view factor computer code for axisymmetric, 2D planar, and 3D geometries with shadowing , 1983 .

[7]  L. E. Scriven,et al.  Study of coating flow by the finite element method , 1981 .

[8]  C. W. Gear,et al.  Simultaneous Numerical Solution of Differential-Algebraic Equations , 1971 .

[9]  T. Surek,et al.  The growth of shaped crystals from the melt , 1980 .

[10]  Jeffrey J. Derby,et al.  On the quasi-steady-state assumption in modeling Czochralski crystal growth , 1988 .

[11]  Claes Johnson,et al.  Finite element methods for linear hyperbolic problems , 1984 .

[12]  Robert A. Brown,et al.  LU decomposition of matrices with augmented dense constraints , 1987 .

[13]  T. Surek,et al.  Theory of shape stability in crystal growth from the melt , 1976 .

[14]  William H. Raymond,et al.  Selective Damping in a Galerkin Method for Solving Wave Problems with Variable Grids , 1976 .

[15]  W. Zulehner,et al.  Czochralski growth of silicon , 1983 .

[16]  Hisham M. Ettouney,et al.  Finite-element methods for steady solidification problems , 1983 .

[17]  Daniel R. Lynch,et al.  Unified approach to simulation on deforming elements with application to phase change problems , 1982 .

[18]  G. W. Green,et al.  The meniscus in Czochralski growth , 1974 .

[19]  Jeffrey J. Derby,et al.  Radiative heat exchange in Czochralski crystal growth , 1987 .

[20]  Liang-Heng Chen,et al.  Bifurcation characteristics of nonlinear systems under conventional pid control , 1984 .

[21]  G. C. Joyce,et al.  The weighing method of automatic Czochralski crystal growth: II. Control equipment , 1977 .

[22]  Jeffrey J. Derby,et al.  Thermal-capillary analysis of Czochralski and liquid encapsulated Czochralski crystal growth: I. Simulation , 1986 .

[23]  Thomas J. R. Hughes,et al.  Finite element methods for first-order hyperbolic systems with particular emphasis on the compressible Euler equations , 1984 .

[24]  J. Derby,et al.  Thermal-capillary analysis of Czochralski and liquid encapsulated Czochralski crystal growth. II - Processing strategies , 1986 .

[25]  J. Derby,et al.  On the dynamics of Czochralski crystal growth , 1987 .

[26]  D. Hurle Control of diameter in Czochralski and related crystal growth techniques , 1977 .

[27]  G. Strang,et al.  An Analysis of the Finite Element Method , 1974 .

[28]  John M. Sullivan,et al.  Heat conservation in deforming element phase change simulation , 1985 .

[29]  B. Chalmers,et al.  The direction of growth of the surface of a crystal in contact with its melt , 1975 .

[30]  C. W. Gear,et al.  Numerical initial value problem~ in ordinary differential eqttations , 1971 .

[31]  Jeffrey J. Derby,et al.  Finite Element Analysis of a Thermal‐Capillary Model for Liquid Encapsulated Czochralski Growth , 1985 .