Strategies for the coupling of global and local crystal growth models

Abstract The modular coupling of existing numerical codes to model crystal growth processes will provide for maximum effectiveness, capability, and flexibility. However, significant challenges are posed to make these coupled models mathematically self-consistent and algorithmically robust. This paper presents sample results from a coupling of the CrysVUn code, used here to compute furnace-scale heat transfer, and Cats2D, used to calculate melt fluid dynamics and phase-change phenomena, to form a global model for a Bridgman crystal growth system. However, the strategy used to implement the CrysVUn–Cats2D coupling is unreliable and inefficient. The implementation of under-relaxation within a block Gauss–Seidel iteration is shown to be ineffective for improving the coupling performance in a model one-dimensional problem representative of a melt crystal growth model. Ideas to overcome current convergence limitations using approximations to a full Newton iteration method are discussed.

[1]  C. Martínez-Tomás,et al.  CdTe crystal growth process by the Bridgman method: numerical simulation , 2001 .

[2]  W. Shyy,et al.  Numerical simulation of CdTe vertical Bridgman growth , 1997 .

[3]  Gabriel Wittum,et al.  Numerical Treatment of Coupled Systems , 1995 .

[4]  D. E. Bornside,et al.  Toward an integrated analysis of czochralski growth , 1989 .

[5]  Arun Pandy,et al.  Fixed‐point convergence of modular, steady‐state heat transfer models coupling multiple scales and phenomena for melt–crystal growth , 2006 .

[6]  Marcel Crochet,et al.  Global modelling of heat transfer in crystal growth furnaces , 1990 .

[7]  Åke Björck,et al.  Numerical Methods , 1995, Handbook of Marine Craft Hydrodynamics and Motion Control.

[8]  K. Kakimoto,et al.  Partly three-dimensional global modeling of a silicon Czochralski furnace. I. Principles, formulation and implementation of the model , 2005 .

[9]  Charbel Farhat,et al.  Partitioned analysis of coupled mechanical systems , 2001 .

[10]  G. Müller,et al.  Process modeling of the industrial VGF growth process using the software package CrysVUN , 2000 .

[11]  Arun Pandy,et al.  Three-dimensional imperfections in a model vertical Bridgman growth system for cadmium zinc telluride , 2004 .

[12]  W. Shyy,et al.  Multi-zone simulation of the bridgman growth process of β-NiAl crystal , 1996 .

[13]  Miguel Cervera,et al.  On the computational efficiency and implementation of block-iterative algorithms for nonlinear coupled Problems , 1996 .

[14]  G. Müller,et al.  Crystal growth : from fundamentals to technology , 2004 .

[15]  Arun Pandy,et al.  Analysis of the growth of cadmium zinc telluride in an electrodynamic gradient freeze furnace via a self-consistent, multi-scale numerical model , 2005 .

[16]  Tony F. Chan An Approximate Newton Method for Coupled Nonlinear Systems , 1985 .

[17]  Georg Müller,et al.  Development of a new powerful computer code CrysVUN++ especially designed for fast simulation of bulk crystal growth processes , 1999 .

[18]  A. Mühlbauer,et al.  Numerical investigation of silicon melt flow in large diameter CZ-crystal growth under the influence of steady and dynamic magnetic fields , 2001 .

[19]  Yu.N. Makarov,et al.  Advances in the simulation of heat transfer and prediction of the melt-crystal interface shape in silicon CZ growth , 2004 .

[20]  H. Matthies,et al.  Partitioned but strongly coupled iteration schemes for nonlinear fluid-structure interaction , 2002 .

[21]  Hermann G. Matthies,et al.  Algorithms for strong coupling procedures , 2006 .

[22]  Georg Müller,et al.  Global simulation of heat transport, including melt convection in a Czochralski crystal growth process-combined finite element/finite volume approach , 1993 .