3D Finite Volume Scheme for Czochralski Single Crystal Growth

Czochralski (Cz) technique, which is used for growing single crystals, has dominated the production of single crystals for electronic applications. The Cz growth process involves multiple phases, moving interface and three-dimensional behavior. Much has been done to study these phenomena by means of numerical methods as well as experimental observations. A three-dimensional curvilinear finite volume based algorithm has been developed to model the Cz process. A body-fitted transformation based approach is adopted in conjunction with a multizone adaptive grid generation (MAGG) technique to accurately handle the three-dimensional problems of phase-change in irregular geometries with free and moving surfaces. The multizone adaptive model is used to perform a three-dimensional simulation of the Cz growth of silicon single crystals.Since the phase change interface are irregular in shape and they move in response to the solution, accurate treatment of these interfaces is important from numerical accuracy point of view. The multizone adaptive grid generation (MAGG) is the appropriate scheme for this purpose. Another challenge encountered is the moving and periodic boundary conditions, which is essential to the numerical solution of the governing equations. Special treatments are implemented to impose the periodic boundary condition in a particular direction and to determine the internal boundary position and shape varying with the combination of ambient physicochemical transport process and interfacial dynamics. As indicated above that the applications and processes characterized by multi-phase, moving interfaces and irregular shape render the associated physical phenomena three-dimensional and unsteady. Therefore a generalized 3D model rather than a 2D simulation, in which the governing equations are solved in a general non-orthogonal coordinate system, is constructed to describe and capture the features of the growth process. All this has been implemented and validated by using it to model the low pressure Cz growth of silicon. Accuracy of this scheme is demonstrated by agreement of simulation data with available experimental data. Using the quasi-steady state approximation, it is shown that the flow and temperature fields in the melt under certain operating conditions become asymmetric and unsteady even in the absence of extrinsic sources of asymmetry. Asymmetry in the flow and temperature fields, caused by high shear initiated phenomena, affects the interface shape in the azimuthal direction thus results in the thermal stress distribution in the vicinity, which has serious implications from crystal quality point of view.

Perturbational finite volume scheme for the one-dimensional Navier-Stokes equations

Starting from the second-order finite volume scheme,though numerical value perturbation of the cell facial fluxes, the perturbational finite volume (PFV) scheme of the Navier-Stokes (NS) equations for compressible flow is developed in this paper. The central PFV scheme is used to compute the one-dimensional NS equations with shock wave.Numerical results show that the PFV scheme can obtain essentially non-oscillatory solution.

A Hybrid Scheme Based on Finite Element/Volume Methods for Two Immiscible Fluid Flows

We have successfully extended our implicit hybrid finite element/volume (FE/FV) solver to flows involving two immiscible fluids. The solver is based on the segregated pressure correction or projection method on staggered unstructured hybrid meshes. An intermediate velocity field is first obtained by solving the momentum equations with the matrix-free implicit cell-centered FV method. The pressure Poisson equation is solved by the node-based Galerkin FE method for an auxiliary variable. The auxiliary variable is used to update the velocity field and the pressure field. The pressure field is carefully updated by taking into account the velocity divergence field. This updating strategy can be rigorously proven to be able to eliminate the unphysical pressure boundary layer and is crucial for the correct temporal convergence rate. Our current staggered-mesh scheme is distinct from other conventional ones in that we store the velocity components at cell centers and the auxiliary variable at vertices. The fluid interface is captured by solving an advection equation for the volume fraction of one of the fluids. The same matrix-free FV method, as the one used for momentum equations, is used to solve the advection equation. We will focus on the interface sharpening strategy to minimize the smearing of the interface over time. We have developed and implemented a global mass conservation algorithm that enforces the conservation of the mass for each fluid.

High Order Multi-Moment Constrained Finite Volume Method. Part I: Basic Formulation

A new high-order finite volume method based on local reconstruction is presented in this paper. The method, so-called the multi-moment constrained finite volume (MCV) method, uses the point values defined within single cell at equally spaced points as the model variables (or unknowns). The time evolution equations used to update the unknowns are derived from a set of constraint conditions imposed on multi kinds of moments, i.e. the cell-averaged value and the point-wise value of the state variable and its derivatives. The finite volume constraint on the cell-average guarantees the numerical conservativeness of the method. Most constraint conditions are imposed on the cell boundaries, where the numerical flux and its derivatives are solved as general Riemann problems. A multi-moment constrained Lagrange interpolation reconstruction for the demanded order of accuracy is constructed over single cell and converts the evolution equations of the moments to those of the unknowns. The presented method provides a general framework to construct efficient schemes of high orders. The basic formulations for hyperbolic conservation laws in 1- and 2D structured grids are detailed with the numerical results of widely used benchmark tests. (C) 2009 Elsevier Inc. All rights reserved.

COLOR-HOLOGRAPHY USING THE ANGULAR SELECTIVITY OF VOLUME RECORDING MEDIA

A display hologram of an object can be recorded and reconstructed in three primary colors if the angular selectivity of volume recording media is exploited. Three holograms are recorded in the same medium, each at a different primary color. These three holograms are reconstructed by simultaneous illumination of the hologram with the original reference beams. By proper choice of the angles that the reference beams make to the hologram, it is possible to suppress strongly cross talk between the different reconstructions (e.g., the red object reconstruction in green light). The technique exhibits high resolution, high diffraction efficiency, and vivid colors. Through the addition of three holographically recorded volume gratings it is possible to reconstruct the hologram with a beam of white light. The saturation and brightness of each primary color in the reconstruction can be adjusted by selection of an appropriate thickness for the corresponding grating.

High-efficiency volume hologram recording with a pulsed signal beam

We present an efficient photorefractive volume hologram recording technique with a pulsed signal beam and continuous reference-beam illumination. The grating envelope can be simply controlled by manipulation of the duty cycle of the signal beam. Thus, for any grating coupling strength and different initial reference-signal intensity ratios, the diffraction efficiency can be maximized with this technique and can be greatly increased in comparison with that of the conventional recording technique. (C) 1998 Optical Society of America.

Optimizations for the uniformity of the non-volatile volume grating in doubly doped LiNbO3 crystals

The formation of the non-uniformity of the non-volatile volume grating in doubly doped LiNbO3 crystals is studied in detail. We find that the non-uniformity of the grating is mainly caused by strong ultraviolet light absorption, and the average saturation space-charge field is small and the diffraction efficiency is low as a result of the non-uniformity of the grating. In order to optimize the uniformity of the grating, we propose the recording scheme by using two sensitizing beams simultaneously from the two opposite sides of the crystals. Theoretical simulations and experimental verifications are performed. Results show that the well uniformed grating with high diffraction efficiency can be obtained by using this optimization scheme. (c) 2004 Elsevier B.V. All rights reserved.

Pulse shaping properties of volume holographic gratings in anisotropic media

Based on a modified coupled wave theory, the pulse shaping properties of volume holographic gratings (VHGs) in anisotropic media VHGs are studied systematically. Taking photorefractive LiNbO3 crystals as an example, the combined effect that the grating parameters, the dispersion and optical anisotropy of the crystal, the pulse width, and the polarization state of the input ultrashort pulsed beam (UPB) have on the pulse shaping properties are considered when the input UPB with arbitrary polarization state propagates through the VHG. Under the combined effect, the diffraction bandwidth, pulse profiles of the diffracted and transmitted pulsed beams, and the total diffraction efficiency are shown. The studies indicate that the properties of the shaping of the o and e components of the input UPB in the crystal are greatly different; this difference can be used for pulse shaping applications. (c) 2006 Optical Society of America.

Wavefront conversion by local volume holograms between cylindrical and plane waves with reconstruction at a different readout wavelength

Based on the two-dimensional coupled-wave theory, the wavefront conversion between cylindrical and plane waves by local volume holograms recorded at 632.8 nm and reconstructed at 800 nm is investigated. The proposed model can realize the 90 degrees holographic readout at a different readout wavelength. The analytical integral solutions for the amplitudes of the space harmonics of the field inside the transmission geometry are presented. The values of the off-Bragg parameter at the reconstructed process and the diffracted beam's amplitude distribution are analysed. In addition, the dependences of diffraction efficiency on the focal length of the recording cylindrical wave and on the geometrical dimensions of the grating are discussed. Furthermore, the focusing properties of this photorefractive holographic cylindrical lens are analysed.

Anisotropic diffraction of stratified volume holographic grating in photorefractive media

We have studied the anisotropic diffraction properties of the stratified volume holographic gratings recorded in photorefractive media using the anisotropic coupled wave theory. It is shown that the diffraction efficiency of such system exhibit the uniform periodic Bragg selectivity properties. In addition the dependence of the stratified volume holographic optical elements (SVHOEs) diffraction properties on the buffer-layer thickness, grating-layer thickness, number of modulation layers, and total thickness of system are discussed in detail. (c) 2005 Elsevier GrnbH. All rights reserved.

Diffraction properties of volume holographic gratings for an ultrashort pulsed beam with different polarization states readout

The diffraction properties of volume holographic gratings are studied when the gratings are illuminated by an ultrashort pulsed beam with different polarization states. The developed coupled wave theory of Kogelnik is used. Considering the dispersion effect of the grating media, solutions for the diffracted and transmitted intensities, diffraction efficiencies and the bandwidths of the gratings are given in transmission volume holographic gratings and reflection volume holographic gratings. The bandwidths of the gratings are reduced by the dispersion effect of the grating media. They also have different influences on the diffraction of an ultrashort pulsed beam with different polarization states. For different values of the ratio of the spectral bandwidth of the input pulse to that of the grating, the changes of the spectral and temporal distributions of the diffracted intensities, as well as the diffraction efficiencies of the gratings are shown.

Diffraction of an ultrashort pulsed beam with arbitrary polarization state from a volume holographic grating in LiNbO3 crystals

Based on a modified coupled wave theory of Kogelnik, we have studied the diffraction of an ultrashort pulsed beam with an arbitrary polarization state from a volume holographic grating in photorefractive LiNbO3 crystals. The results indicate that the diffracted intensity distributions in the spectral and temporal domains and the diffraction efficiency of the grating are both changed by the polarization state and spectral bandwidth of the input pulsed beam. A method is given of choosing the grating parameters and input conditions to obtain a large variation range of the spectral bandwidth of the diffracted pulsed beam with an appropriate diffraction efficiency. Our study presents a possibility of using a volume holographic grating recorded in anisotropic materials to shape a broadband ultrashort pulsed beam by modulating its polarization state.

Analysis of an ultrashort pulsed finite beam diffracted by volume gratings

We obtain analytical solutions of the coupled wave equations that describe the Bragg diffraction of ultrashort pulsed finite beams by a thick planar grating, using two-dimensional coupled wave theory. The diffraction properties for the case of an ultrashort pulsed finite beam with Gaussian profiles in both the time and spatial domains are investigated. The spectral bandwidth of the diffracted beam, the Bragg selectivity bandwidth and the diffraction efficiency of the volume grating are influenced by the geometry parameter and the input bandwidth. Therefore extra attention should be paid to designing optical elements based on volume gratings for use with ultrashort pulsed waves in applications of pulse shaping and processing.