CRYSTAL
Silicon Crystallization Simulation

CRYSTAL is a FORTRAN90 program which simulates the process of creating a crystal of silicon.

CRYSTAL carries out a minimization, for which it calls TOMS611, a minimization algorithm from the ACM TOMS algorithm collection.

Languages:

CRYSTAL is available in a FORTRAN90 version.

Related Data and Programs:

CRYSTAL_PLOT, a FORTRAN90 program which creates graphic images of data from CRYSTAL or CRYSTAL_QED.

CRYSTAL_QED, a FORTRAN90 program which seeks to minimize a certain objective function during the growth of a silicon crystal under the Czochralski process.

TOMS611, a FORTRAN90 library which minimizes a function of several variables.

Reference:

1. Suhas Patankar,
   Numerical Heat Transfer and Fluid Flow,
   Taylor and Francis, 1980,
   ISBN: 0891165223,
   LC: QC320.P37
2. Fumio Shimura,
   Semiconductor Silicon Crystal Technology,
   Academic Press, 1989,
   LC: TK7871.85 S523.

Source Code:

• crystal.f90, the source code.
• crystal.sh, commands to compile, link and load the source code.
• crystal.txt, some development notes.
• crystal_dic.txt, a data dictionary.

Examples and Tests:

• crystal.txt, the output file.
• crystal_rs.txt, a restart file.

List of Routines:

• MAIN is the main program for CRYSTAL.
• DUMMY is a dummy subroutine needed as formal input to CRYFUN.
• CRYFUN evaluates the cost function for the minimizing software.
• ADAPT executes MAGG, the multizone adaptive grid generation.
• CUBIC constructs a cubic spline through data.
• DIFLOW computes the convection-diffusion coefficient.
• DOAREA computes the area of the control volumes.
• FINDP finds the boundary nodes for a Neumann BC.
• FLUX computes the value of the flux of various quantities.
• GAMSOR sets the coefficients for the transport problems.
• GRADNT calculates gradients at the primary nodes.
• INIDAT sets the initial values of certain data.
• INIGRD makes an initial assignment of the grid points XC, YC.
• INITL estimates quantities at control volume interfaces.
• MOVGRD calculates the new position of the interface and free surface.
• OUTPUT prints information about the current solution.
• PMOD extends pressure values from primary nodes to corner nodes.
• PRDAT prints out the initial values of certain data.
• RESID computes the linear equation residual at interior primary nodes.
• RSWRIT writes out restart information.
• SETCST computes a portion of the cost functional.
• SETGEO calculates various geometric quantities.
• SETUP calculates the coefficients of the equations, and solves them.
• SETX locates the primary nodes from the corner nodes.
• SOLVE1 sets or modifies the linear system to be handled by SOLVE2.
• SOLVE2 is the tridiagonal matrix solver.
• TIMESTAMP prints the current YMDHMS date as a time stamp.

You can go up one level to the FORTRAN90 source codes.