• 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 Soil and Groundwater Environment
• /
• v.13 no.1
• /
• pp.77-91
• /
• 2008

In this study, we developed a technique of applying DRASTIC, which is the most widely used tool for estimation of groundwater vulnerability to the aqueous phase contaminant infiltrated from the surface, and a groundwater flow model jointly to assess groundwater contamination potential. The developed technique is then applied to Buyeo-eup area in Buyeo-gun, Chungcheongnam-do, Korea. The input thematic data of a depth to water required in DRASTIC model is known to be the most sensitive to the output while only a few observations at a few time schedules are generally available. To overcome this practical shortcoming, both steady-state and transient groundwater level distributions are simulated using a finite difference numerical model, MODFLOW. In the application for the assessment of groundwater vulnerability, it is found that the vulnerability results from the numerical simulation of a groundwater level is much more practical compared to cokriging methods. Those advantages are, first, the results from the simulation enable a practitioner to see the temporally comprehensive vulnerabilities. The second merit of the technique is that the method considers wide variety of engaging data such as field-observed hydrogeologic parameters as well as geographic relief. The depth to water generated through geostatistical methods in the conventional method is unable to incorporate temporally variable data, that is, the seasonal variation of a recharge rate. As a result, we found that the vulnerability out of both the geostatistical method and the steady-state groundwater flow simulation are in similar patterns. By applying the transient simulation results to DRASTIC model, we also found that the vulnerability shows sharp seasonal variation due to the change of groundwater recharge. The change of the vulnerability is found to be most peculiar during summer with the highest recharge rate and winter with the lowest. Our research indicates that numerical modeling can be a useful tool for temporal as well as spatial interpolation of the depth to water when the number of the observed data is inadequate for the vulnerability assessments through the conventional techniques.

• Geophysics and Geophysical Exploration
• /
• v.13 no.1
• /
• pp.1-8
• /
• 2010

Seismic reflection surveying is one of the most widely used and effective techniques for coal seam structure delineation and risk mitigation for underground longwall mining. However, the ability of the method can be compromised by the presence of volcanic cover. This problem arises within parts of the Bowen and Sydney Basins of Australia and seismic surveying can be unsuccessful. As a consequence, such areas are less attractive for coal mining. Techniques to improve the success of seismic surveying over basalt flows are needed. In this paper, we use elastic wave-equation-based forward modelling techniques to investigate the effects and characteristics of seismic wave propagation under different settings involving changes in basalt properties, its thickness, lateral extent, relative position to the shot position and various forms of inhomogeneity. The modelling results suggests that: 1) basalts with high impedance contrasts and multiple flows generate strong multiples and weak reflectors; 2) thin basalts have less effect than thick basalts; 3) partial basalt cover has less effect than full basalt cover; 4) low frequency seismic waves (especially at large offsets) have better penetration through the basalt than high frequency waves; and 5) the deeper the coal seams are below basalts of limited extent, the less influence the basalts will have on the wave propagation. In addition to providing insights into the issues that arise when seismic surveying under basalts, these observations suggest that careful management of seismic noise and the acquisition of long-offset seismic data with low-frequency geophones have the potential to improve the seismic results.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.33 no.6
• /
• pp.2255-2265
• /
• 2013

This study presents the width function-based Clark model. To this end, rescaled width function with distinction between hillslope and channel velocity is used as time-area curve and then it is routed through linear storage within the framework of not finite difference scheme used in original Clark model but analytical expression of linear storage routing. There are three parameters focused in this study: storage coefficient, hillslope velocity and channel velocity. SCE-UA, one of the popular global optimization methods, is applied to estimate them. The shapes of resulting IUHs from this study are evaluated in terms of the three statistical moments of hydrologic response functions: mean, variance and the third moment about the center of IUH. The correlation coefficients to the three statistical moments simulated in this study against these of observed hydrographs were estimated at 0.995 for the mean, 0.993 for the variance and 0.983 for the third moment about the center of IUH. The shape of resulting IUHs from this study give rise to satisfactory simulation results in terms of the mean and variance. But the third moment about the center of IUH tend to be overestimated. Clark model proposed in this study is superior to the one only taking into account mean and variance of IUH with respect to skewness, peak discharge and peak time of runoff hydrograph. From this result it is confirmed that the method suggested in this study is useful tool to reflect the heterogeneity of drainage path and hydrodynamic parameters. The variation of statistical moments of IUH are mainly influenced by storage coefficient and in turn the effect of channel velocity is greater than the one of hillslope velocity. Therefore storage coefficient and channel velocity are the crucial factors in shaping the form of IUH and should be considered carefully to apply Clark model proposed in this study.

• The Journal of Korean Academy of Prosthodontics
• /
• v.50 no.1
• /
• pp.10-20
• /
• 2012

Purpose: The purpose of this study was to investigate the stress distribution in mandibular implant-supported overdentures and tooth-supported overdentures with telescopic crowns. Materials and methods: The assumption of this study was that there were 2, 3, 4 natural teeth and implants which are located in the second premolar and canine regions in various distributed conditions. The mandible, teeth (or implants and abutments), and connectors are modeled, and analyzed with the commercial software, ANSYS Version 10.1. Stress distribution was evaluated under 150 N vertical load bilaterally on 3 experimental conditions - between canine areas, canine and $2^{nd}$ premolars, 10 mm posterior to $2^{nd}$ premolars. Results: Overall, the case of the implant group showed more stress than the case of the teeth group in stress distribution to bone. In stress distribution to superstructures of tooth and implants, there was no significant difference between TH group and IM group and the highest stress appeared in TH-IV and IM-IV. The stress caused from bar was much higher than those of implant and tooth. TH group showed less stress than IM group in stress distribution to abutment teeth and implant. Conclusion: The results shows that it is crucial to make sure that distance between impact loading point and abutment tooth does not get too far apart, and if it does, it is at best to set abutment tooth on premolar tooth region. It will be necessary to conduct more experiments on effects on implants, natural teeth and bone, in order to apply these results to a clinical treatment.

• Journal of Radiation Protection and Research
• /
• v.41 no.1
• /
• pp.31-39
• /
• 2016

Background: It is necessary to consider the overall countermeasure for analysis of nuclear activities according to the increase of the nuclear facilities like nuclear power and reprocessing plants in the neighboring countries including China, Taiwan, North Korea, Japan and South Korea. South Korea and comprehensive nuclear-test-ban treaty organization (CTBTO) are now operating the monitoring instruments to detect radionuclides released into the air. It is important to estimate the origin of radionuclides measured using the detection technology as well as the monitoring analysis in aspects of investigation and security of the nuclear activities in neighboring countries. Materials and methods: A three-dimensional forward/backward trajectory model has been developed to estimate the origin of radionuclides for a covert nuclear activity. The developed trajectory model was composed of forward and backward modules to track the particle positions using finite difference method. Results and discussion: A three-dimensional trajectory model was validated using the measured data at Chernobyl accident. The calculated results showed a good agreement by using the high concentration measurements and the locations where was near a release point. The three-dimensional trajectory model had some uncertainty according to the release time, release height and time interval of the trajectory at each release points. An atmospheric dispersion model called long-range accident dose assessment system (LADAS), based on the fields of regards (FOR) technique, was applied to reduce the uncertainties of the trajectory model and to improve the detective technology for estimating the radioisotopes emission area. Conclusion: The detective technology developed in this study can evaluate in release area and origin for covert nuclear activities based on measured radioisotopes at monitoring stations, and it might play critical tool to improve the ability of the nuclear safety field.

• Economic and Environmental Geology
• /
• v.45 no.2
• /
• pp.71-78
• /
• 2012

Illite-smectite mixed layers (I-S) occurring authigenically in diagenetic and hydrothermal environments reacts toward more illite-rich phases as temperature and potassium ion concentration increase. For that reason, I-S is often used as geothermometry and/or geochronometry at the field of hydrocarbons or ore minerals exploration. Generally, I-S shows X-ray powder diffraction (XRD) patterns of ultra-thin lamellar structures, which consist of restricted numbers of sillicate layers (normally, 5 ~ 15 layers) stacked in parallel to a-b planes. This ultra-thinness is known to decrease I-S expandability (%S) rather than theoretically expected one (short-stacking effect). We attempt here to quantify the short stacking effect of I-S using the difference of two types of expandability: one type is a maximum expandability ($%S_{Max}$) of infinite stacks of fundamental particles (physically inseparable smallest units), and the other type is an expandability of finite particle stacks normally measured using X-ray powder diffraction (XRD) ($%S_{XRD}$). Eleven I-S samples from the Geumseongsan volcanic complex, Uiseong, Gyeongbuk, have been analyzed for measuring $%S_{XRD}$ and average coherent scattering thickness (CST) after size separation under 1 ${\mu}m$. Average fundamental particle thickness ($N_f$) and $%S_{Max}$ have been determined from $%S_{XRD}$ and CST using inter-parameter relationships of I-S layer structures. The discrepancy between $%S_{Max}$ and $%S_{XRD}$ (${\Delta}%S$) suggests that the maximum short-stacking effect happens approximately at 20 $%S_{XRD}$, of which point represents I-S layer structures consisting of ca. average 3-layered fundamental particles ($N_f{\approx}3$). As a result of inferring the $%S_{XRD}$ range of each Reichweite using the $%S_{XRD}$ vs. $N_f$ diagram of Kang et al. (2002), we can confirms that the fundamental particle thickness is a determinant factor for I-S Reichweite, and also that the short-stacking effect shifts the $%S_{XRD}$ range of each Reichweite toward smaller $%S_{XRD}$ values than those that can be theoretically prospected using junction probability.

• Journal of the Korea Concrete Institute
• /
• v.26 no.4
• /
• pp.441-448
• /
• 2014

From several researches, recently, it was found that using hollowed precast concrete (HPC) column made more compact concrete casting in joint region possible than using normal solid PC (Precast concrete) column. Therefore, the rigidity of joints can be improved like those of monolithic reinforced concrete (RC). After filling the hollow with grout concrete, however, it is expected that the HPC column behaviors like composite structure since PC element and grout concrete have different materials as well as there is a contact surface between two elements. These may affect the structural behavior and strength of the composite column. A compressive strength test was performed for the HPC column with parameter of hollow ratio for the case with and without grout in the hollow and the result is presented in this paper. The hollow ratios in the test are 35, 50 and 59% of whole section of column. Concentrated axial force was applied to top of the specimens supported as pin connection for both ends. In addition, finite element (FE) analysis was performed to simulate the failure behavior of HPC column for axial compression. As a result, it was found that the hollow ratio did not affect the initial stiffness of HPC filled with grout regardless of the strength difference of HPC and grout. However the strength was increased inversely corresponding to the hollow ratio. The structural capacity of HPC without grout closely related to the hollow size. Especially, the local collapse governs the overall failure when the thickness of HPC is too thin. Based on these effect, a suitable equation was suggested for calculation of the compressive strength of HPC column with or without grout. FE analysis considering the contact surface between HPC and grout produced a good result matched to the test result.

• Korean Journal of Agricultural Science
• /
• v.9 no.1
• /
• pp.357-370
• /
• 1982

Recent concern regarding price and availability of fossil fuels has spurred the interest in alternative sources for farm crop drying. Among the available options such as biomass energy, wind power, nuclear energy and solar energy etc., the increasing attention is being directed to the utilization of heat from solar energy especially for farm crop drying. Even though solar energy is dispersed over a large land area and only a relatively small amount of energy can be simply collected, the advantages of solar energy is that the energy is free, non-polluting. The study reported here was designed to help supply the informations for the development of simple and relatively inexpensive solar warehouse for farm crop drying and storage. Specifically, the objectives of this study were to determine the performance of the solar collector fabricated, to compare solar supplemented heat drying with natural air drying and to develop a simulation model of temperature in stored grain, which can be used to study the effects due to changes in ambient air temperature. For those above objectives, solar collector was fabricated from available materials. Corrugated steel galvanized sheet, painted flat black, was used as absorbers and clear 0.2mm polyethylene sheet was the cover material. The warehouse for rough rice drying and storage was constructed with concrete block, and the solar collector was used as the roof of warehouse instead of original roofing system of it. The results obtained in this study were as follows: 1. The thermal efficiency of the solar collector was average 26 percent and the overall heat transfer coefficient of the collector was approximately $25kJ/hr.m^2\;^{\circ}K$. 2. Solar heated air was sufficient to dry one cubic meter of rough rice from 23.5 to 15.0 percent in 7 days and natural air was able to dry the same amount of rough rice from 20.0 to 5 percent in l2 days. 3. Drying with solar heat reduced the required drying time to dry the same amount of rough rice into a half compared to natural air drying, but overdrying problems of the bottom layer were so severe that these problems should be thoroughly analyzed. 4. Simulation model of temperature in stored grain was developed and the results of predicted temperature agreed well with test results. 5. Based on those simulated temperature, changes in the grain-temperature were a large at the points of the wallside and the damage of the grain would be severe at the contact area of wall.

• Proceedings of the Korean Society for Rock Mechanics Conference
• /
• 2000.09a
• /
• pp.183-193
• /
• 2000

A series of studies such as geological logging data analysis, detailed geological survey, rock mass evaluation, in-situ and laboratory tests, rock strength and mechanical properties of the rock were concerned. The stability of the slope were carried out inorder to design the pit slope and individual benches using the stereographic projection analysis and numerical methods in Roto Pit of Pasir coal fetid. The bedding plane was one of the major discontinuities in the Roto Pit and the dip of which is about $60^{\circ}$in the northern part and $83^{\circ}$in the southern part. The dip of bedding becomes steeper from north to south. The plane and toppling failures are presented in many slopes. In laboratory test the average uniaxial compressive strength of mudstone was 9 MPa and that of weak sandstone was 10 MPa. In-situ test showed that the rocks of Roto north mining area are mostly weak enough to be classified in grade from R2(weak) to R3(medium strong weak) and the coal is classified in grades from R1(Very weak) to R2(Weak). The detailed stability analysis were carried out on 4 areas of Roto north(east, west, south and north), and 2 areas of Roto south(east and west). In this paper, the minimum factor of safety was set to 1.2 which is a general criterion for open pit mines. Using the stereographic projection analysis and the limit equilibrium method, slope angles were calculated as 30~$36^{\circ}$for a factor of safety greater than 1.2. Then these results were re-evaluated by numerical analysis using FLAC. The final slope angles were determined by rational described abode. A final slope of 34 degrees can guarantee the stability for the eastern part of the Roto north area, 33 degrees for the western part, 35 degrees for the northern part and 35 degrees for the southern part. For the Roto south area, 36 degrees was suggested for both sides of the pit. Once the pit slope is designed based on the stability analysis and the safety measures. the stability of 니ope should be checked periodically during the mining operations. Because the slope face will be exposed long time to the rain fall, a study such aspreventive measures against weathering and erosion is highly recommended to be implemented.