• Journal of the Korean Crystal Growth and Crystal Technology
• /
• v.10 no.3
• /
• pp.177-184
• /
• 2000

The effect of the buyancy and thermocapillarity for differnent aspect ratio of flow field on melt motion and mass transfer has been numerically investigated in magnetic Czochralski crystal growth of silicon. During the process of crystal growth, the melt depth of crucible reduces so the aspect ratio of flow field also reduces. Therefore the shape of magnetic field of the flow field changes and the flow pattern also changes significantly. Together with the melt flow which forms the Marangoni convection (or thermocapillary flow) that comes from the inside the flow field, a flow circulation is observed near the corner close both to the crucible wall and the free surface. Due to this circulation, buoyancy effect has been turned out to be local rather than global. As the aspect ratio decreases, the radial component of the magnetic field prevails compared with the axial component in the flow field. Under the influence of this magnetic field, the melt flow and the temperature distribution in a meridional plane tend to depend on the radial position. As the aspect ratio decreases, the temperature gradient near the edge of the crystal decreases yielding smaller thermocapillarity, and the oxygen concentration near the crystal and the oxygen incorporation rate also decrease.

• Proceedings of the Korea Association of Crystal Growth Conference
• /
• 1996.06a
• /
• pp.179-200
• /
• 1996

The intrinsic instabilities of fluid flow occurred in the melt of the Czochralski crystal growth system Czochralski method, asymmetric flow patterns and temperature profiles in the melt have been studied by many researchers. The idea that the non-symmetric structure of the growing equipment is responsible for the asymmetric profiles is usually accepted at the first time. However further researches revealed that some intrinsic instabilities not related to the non-symmetric equipment structure in the melt could also appear. Ristorcelli had pointed out that there are many possible causes of instabilities in the melt. The instabilities appears because of the coupling effects of fluid flow and temperature profiles in the melt. Among the instabilities, the B nard type instabilities with no or low crucible rotation rates are analyzed by the visualizing experiments using X-ray radiography and the 3-D numerical simulation in this study. The velocity profiles in the Silicon melt at different crucible rotation rates were measured using X-ray radiography method using tungsten tracers in the melt. The results showed that there exits two types of fluid flow mode. One is axisymmetric flow, the other is asymmetric flow. In the axisymmetric flow, the trajectory of the tracers show torus pattern. However, more exact measurement of the axisymmetrc case shows that this flow field has small non-axisymmetric components of the velocity. When fluid flow is asymmetric, the tracers show random motion from the fixed view point. On the other hand, when the observer rotates to the same velocity of the crucible, the trajectory of the tracer show a rotating motion, the center of the motion is not same the center of the melt. The temperature of a point in the melt were measured using thermocouples with different rotating rates. Measured temperatures oscillated. Such kind of oscillations are also measured by the other researchers. The behavior of temperature oscillations were quite different between at low rotations and at high rotations. Above experimental results means that the fluid flow and temperature profiles in the melt is not symmetric, and then the mode of the asymmetric is changed when rotation rates are changed. To compare with these experimental results, the fluid flow and temperature profiles at no rotation and 8 rpm of crucible rotation rates on the same size of crucible is calculated using a 3-dimensional numerical simulation. A finite different method is adopted for this simulation. 50×30×30 grids are used. The numerical simulation also showed that the velocity and flow profiles are changed when rotation rates change. Futhermore, the flow patterns and temperature profiles of both cases are not axisymmetric even though axisymmetric boundary conditions are used. Several cells appear at no rotation. The cells are formed by the unstable vertical temperature profiles (upper region is colder than lower part) beneath the free surface of the melt. When the temperature profile is combined with density difference (Rayleigh-B nard instability) or surface tension difference (Marangoni-B nard instability) on temperature, cell structures are naturally formed. Both sources of instabilities are coupled to the cell structures in the melt of the Czochralski process. With high rotation rates, the shape of the fluid field is changed to another type of asymmetric profile. Because of the velocity profile, isothermal lines on the plane vertical to the centerline change to elliptic. When the velocity profiles are plotted at the rotating view point, two vortices appear at the both sides of centerline. These vortices seem to be the main reason of the tracer behavior shown in the asymmetric velocity experiment. This profile is quite similar to the profiles created by the baroclinic instability on the rotating annulus. The temperature profiles obtained from the numerical calculations and Fourier transforms of it are quite similar to the results of the experiment. bove esults intend that at least two types of intrinsic instabilities can occur in the melt of Czochralski growing systems. Because the instabilities cause temperature fluctuations in the melt and near the crystal-melt interface, some defects may be generated by them. When the crucible size becomes large, the intensity of the instabilities should increase. Therefore, to produce large single crystals with good quality, the behavior of the intrinsic instabilities in the melt as well as the effects of the instabilities on the defects in the ingot should be studied. As one of the cause of the defects in the large diameter Silicon single crystal grown by the

• Journal of Biomedical Engineering Research
• /
• v.28 no.3
• /
• pp.397-409
• /
• 2007

Owing to needs of biomechanical comprehension and analysis to obtain various medical treatment designs which are related with the spine in order to cure and diagnose LBP patients, the FE modeling and nonlinear analysis of lumbosacrum including a partial ilium and iliolumbar ligaments, were carried out. First, we investigated whether the geometrical configuration of vertebrae displayed by DICOM slice files is regular and normal condition. After constructing spinal vertebrae including a partial ilium, a sacrum and five lumbars (from L1 to L5)with anatomical shape reconstructed using softwares such as image modeler and CAD modeler, we added iliolumbar ligaments, lumbar ligaments, discs and facet joints, etc.. And also, we assigned material property and discretized the model using proper finite element types, thus it was completely modeled through the above procedure. For the verification of each segment, average sagittal ROM, average coronal ROM and average transversal ROM under various loading conditions(${\pm}10Nm$), average vertical displacement under compression(400N), ALL(Anterior Longitudinal Ligament) and PLL(Posterior Longitudinal Ligament) force at L12 level, strains of seven ligaments on sagittal plane at L45 level and maximal strain of disc fibers according to various loading conditions at L45 level, etc., they were compared with experimental results. For the verification of multilevel-lumbosacrum spine including partial ilium and iliolumbar ligaments, the cases with and without iliolumbar ligaments were compared with ROM of experiment. The results were obtained from analysis of the verified FE model as follows: I) Iliolumbar ligaments played a stabilizing role as mainly posterior iliolumbar ligaments under flexion and as both posterior and anterior iliolumbar ligaments of one side under lateral bending. 2) The iliolumbar ligaments decreased total ROM of 1-8% in total model according to various motion conditions, which changed facet contact forces of L5S level by approximately 0.8-1.4 times and disc forces of L5S level by approximately 0.8-1.5 times more than casewithout ilioligaments, under various loading conditions. 3) The force of lower discs such as L45 and L5S was bigger than upper discs under flexion, left and right bending and left and right twisting, except extension. 4) It was predicted that strains of posterior ligaments among iliolumbar ligaments would produce the maximum 16% under flexion and the maximum 10% under twisting. 5) It's expected that this present model applies to the development and design of artificial disc, since it was comparatively in agreement with the experimental datum.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.7 no.4
• /
• pp.192-198
• /
• 2004

A numerical method is developed to solve a two-dimensional diffraction problem for a body located in a sediment pocket where a heavier muddy water is trapped. In the present study, the wave exciting forces acting on a submerged body in the water-sediment interface by an incident wave is investigate. It is assumed that the heavier mud is trapped locally in a sediment pocket. A mathematical formulation is made in the scope of the potential theory. The fluid is assumed to be inviscid, incompressible and its motion irrotational. The boundary conditions on the unknown free surface and interface are linearized. As a method of solution, the localized finite-element method is adopted. In the method, the computation domain is reduced by utilizing the complete set of analytic solutions known in the infinite subdomain to be truncated by introduction of an appropriate juncture conditions. The main advantage of this method is that any complex geometry of the boundaries can be easily accommodated. Computations are carried out for mono-chromatic plane progressive surface waves normally incident on the domain. Numerical results are compared with those obtained by Lassiter based on Schwingers variational method. Good Agreements are obtained in general. Another numerical computations are made for the cases with and without a body in the sediment pocket.

• Korean Journal of Digital Imaging in Medicine
• /
• v.16 no.1
• /
• pp.13-19
• /
• 2014

Purpose: FLAIR image is beneficial for the diagnosis of various bran diseases including ischemic CVS, brain tumors and infections. However the border between the legion of brain metastasis and surrounding edema may not be clear. Therefore, this study aims to investigate the practical benefits of delayed imaging by comparing the image from a patient with brain metastasis before a contrast enhancement and the image 10 minutes after a contrast enhancement. Materials and methods: Of the 92 people who underwent MRI brain metastases in suspected patients 13 people in three patients there is no video to target the 37 people confirmed cases, and motion artifacts brain metastases in our hospital June-December 2013, 18 people measurement position except for the three incorrect patient (male: 11 people, female: 7 people, average age: 60 years) in the target, test equipment, 3.0T MR System (ACHIEVA Release, Philips, I was 8ChannelSENSE Head Coil use Best, and the Netherlands). TR 11000 ms, TE 125 ms, TI2800 ms, Slice Thickness 5 mm, gap 5 mm, is a Slice number 21, the parameters of the 3D FFE, T2 FLAIR variable that was used to test, TR 8.1 ms, TE 3.7 ms, Slice number 240 I set to. The experiment was conducted by acquiring the FLAIR prior to contrast enhancement (heretofore referred to as Pre FLAIR), and acquiring the 3D FFE CE five minutes after the contrast enhancement, and recomposing the images in an axial plane of S/T 3mm, G 0mm (heretofore referred to as MPR TRA CE). Using the FLAIR 10 minutes after the contrast enhancement (heretofore referred to as Post FLAIR) and Pi-View, a retrospective study was conducted. Using MRIcro on the image of a patient confirmed for his diagnosis, the images before and after the contrast media, as well as the CNR and SNR of the MPR TRA CE images of the lesion and the site absent of lesion were compared and analyzed using a one-way analysis of variance. Results: CNR for Pre FLAIR and Post FLAIR were 34.35 and 60.13, respectively, with MPR TRA CE at 23.77 showing no significant difference (p<0.050). Post-experiment analysis shows a difference between Pre FLAIR and Post FLAIR in terms of CNR (p<0.050), but no difference in CNR between Post FLAIR and MPR TRA CE (p>0.050), indicating that the contrast media had an effect only on Pre FLAIR and Post FLAIR. The SNR for the normal site Pre FLAIR was 106.43, and for the lesion site 140.79. Post FLAIR for the normal site was 107.79, and for the lesion site 167.91. MPR TRA CE for the normal site was 140.23 and for the lesion site 183.19, showing significant difference (p<0.050), and post-experiment analysis shows that there was a difference in SNR only on the lesion sites for Pre FLAIR and Post FLAIR (p<0.050). There was no difference in SNR between the normal site and lesion site for Post FLAIR and MPR TRA CE, indicating no effect from the contrast media (p>0.050). Conclusions: This experiment shows that Post FLAIR has a higher contrast than Pre FLAIR, and a higher SNR for lesions, It was not not statistically significant and MPR TRA CE but CNR came out high. Inspection of post-contrast which is used in a high magnetic field is frequently used images of 3D T1 but, since the signal of the contrast medium and the blood flow is included, this method can be diagnostic accuracy is reduced, it is believed that when used in combination with Post FLAIR, and that can provide video information added to the diagnosis of brain metastases.