• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.26 no.3
• /
• pp.174-183
• /
• 2014

Seabed beneath and near coastal structures may undergo large excess pore water pressure composed of oscillatory and residual components in the case of long durations of high wave loading. This excess pore water pressure may reduce effective stress and, consequently, the seabed may liquefy. If liquefaction occurs in the seabed, the structure may sink, overturn, and eventually increase the failure potential. In this study, to evaluate the liquefaction potential on the seabed, numerical analysis was conducted using the expanded 2-dimensional numerical wave tank to account for an irregular wave field. In the condition of an irregular wave field, the dynamic wave pressure and water flow velocity acting on the seabed and the surface boundary of the composite breakwater structure were estimated. Simulation results were used as input data in a finite element computer program for elastoplastic seabed response. Simulations evaluated the time and spatial variations in excess pore water pressure, effective stress, and liquefaction potential in the seabed. Additionally, the deformation of the seabed and the displacement of the structure as a function of time were quantitatively evaluated. From the results of the analysis, the liquefaction potential at the seabed in front and rear of the composite breakwater was identified. Since the liquefied seabed particles have no resistance to force, scour potential could increase on the seabed. In addition, the strength decrease of the seabed due to the liquefaction can increase the structural motion and significantly influence the stability of the composite breakwater. Due to limitations of allowable paper length, the studied results were divided into two portions; (I) focusing on the dynamic response of structure, acceleration, deformation of seabed, and (II) focusing on the time variation in excess pore water pressure, liquefaction, effective stress path in the seabed. This paper corresponds to (II).

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.26 no.3
• /
• pp.160-173
• /
• 2014

Seabed beneath and near coastal structures may undergo large excess pore water pressure composed of oscillatory and residual components in the case of long durations of high wave loading. This excess pore water pressure may reduce effective stress and, consequently, the seabed may liquefy. If liquefaction occurs in the seabed, the structure may sink, overturn, and eventually increase the failure potential. In this study, to evaluate the liquefaction potential on the seabed, numerical analysis was conducted using the expanded 2-dimensional numerical wave tank to account for an irregular wave field. In the condition of an irregular wave field, the dynamic wave pressure and water flow velocity acting on the seabed and the surface boundary of the composite breakwater structure were estimated. Simulation results were used as input data in a finite element computer program for elastoplastic seabed response. Simulations evaluated the time and spatial variations in excess pore water pressure, effective stress, and liquefaction potential in the seabed. Additionally, the deformation of the seabed and the displacement of the structure as a function of time were quantitatively evaluated. From the results of the analysis, the liquefaction potential at the seabed in front and rear of the composite breakwater was identified. Since the liquefied seabed particles have no resistance to force, scour potential could increase on the seabed. In addition, the strength decrease of the seabed due to the liquefaction can increase the structural motion and significantly influence the stability of the composite breakwater. Due to limitations of allowable paper length, the studied results were divided into two portions; (I) focusing on the dynamic response of structure, acceleration, deformation of seabed, and (II) focusing on the time variation in excess pore water pressure, liquefaction, effective stress path in the seabed. This paper corresponds to (I).

• Journal of the Korean GEO-environmental Society
• /
• v.19 no.4
• /
• pp.5-16
• /
• 2018

In the current study, the engineering behaviour of prebored and precast steel pipe piles was examined from a series of full-scale field measurements by conducting static pile load tests, dynamic pile load tests (EOID and restrike tests) and Class-A and Class-C1 type numerical analysis. The study includes the pile load - settlement relations, allowable pile capacity and shear stress transfer mechanism. Compared to the allowable pile capacity obtained from the static pile load tests, the dynamic pile load tests and the numerical simulation showed surprisingly large variations. Overall among these the restrike tests displayed the best results, however the reliability of the predictions from the numerical analysis was lower than those estimated from the dynamic pile load tests. The allowable pile capacity obtained from the EOID tests and the restrike tests indicated 20.0%-181.0% (avg: 69.3%) and 48.2%-181.1% (avg: 92.1%) of the corresponding measured values from the static pile loading tests, respectively. Furthermore, the computed results from the Class-A type analysis showed the largest scatters (37.1%-210.5%, avg: 121.2%). In the EOID tests, a majority of the external load were carried by the end bearing pile capacity, however, similar skin friction and end bearing capacity in magnitude were mobilised in the restrike tests. The measured end bearing pile capacity from the restrike tests were smaller than was measured from the EOID tests. The present study has revealed that if the impact energy is not sufficient in a restrike test, the end bearing pile capacity most likely will be underestimated. The shear stresses computed from the numerical analysis deviated substantially from the measured pile force distributions. It can be concluded that the engineering behaviour of the pile is heavily affected if a slime layer exists near the pile tip, and that the smaller the stiffness of the slime and the thicker the slime, the greater the settlement of the pile.

• 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

• 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.

• 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.

• 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.

• The Journal of the Petrological Society of Korea
• /
• v.6 no.2
• /
• pp.133-149
• /
• 1997

In the metapelites around the Ogcheon granite, the metamorphic grade increases from the biotite zone through the andalusite zone to the sillimanite zone towards the intrusion contact. In the metabasites around the Boeun granite, the metamorphic grade increases from transitional zone between the greenchist and amphibolite facies through the amphibolite facies to the upper amphibolite facies towards the intrusion contact. In the Doiri area locating near the intrusion contact of the Boeun granite, sillimanite- and andalusite-bearing metapelites are found with in 500 m away from the contact. The evidence described above indicates that the Ogcheon and Boeun granites caused low-P/T type contact metamorphism to the country rocks. The P-T condition of contact metamorphism due to the intrusion of the Ogcheon granite is $540{\pm}40^{circ}C, 2.8{\pm}0.9$ kb. The temperature condition of contact metamorphism due to the intrusion of the Boeun granite is $698{\pm}28^{\circ}C$. The wide compositional range of amphibole and plagioclase in the metabasites around the Boeun granite is due to the immisibility gab of amphibole and plagioclase and unstable relict composition resulted from an incomplete metamorphic reaction. The compositional range of stable amphibole and plagioclase decreases as a metamorphic grade increases due to a close of immiscibility gab. The thermal effect of contact metamorphism due to the intrusion of the Ogcheon and Boeun granites, are calculated using the CONTACT2 program based on a two dimensional finite difference method. In order to estimate the thermal effect of an introduced pluton, a circle with 10 km diameter and a triangle with 20 km side are used for the intrusion geometries of the Ogcheon granite and the Boeun granite, respectively. The results from the field and modeling studies suggest that the intrusion temperatures of the Ogcheon granite close to $800^{\circ}C$ and the intrusion temperature of the Boeun granite is higher than $1000^{\circ}C$. However, the intrusion temperatures can be lower than the suggested temperature, if the geothermal gradient prior to the intrusion of the Ogcheon and Boeun granites was higher than the normal continental grothermal gradient.

• Journal of KIISE:Computer Systems and Theory
• /
• v.30 no.5_6
• /
• pp.324-329
• /
• 2003

The public-key cryptosystems such as Diffie-Hellman Key Distribution and Elliptical Curve Cryptosystems are built on the basis of the operations defined in GF(2$^{m}$ ):addition, subtraction, multiplication and multiplicative inversion. It is important that these operations should be computed at high speed in order to implement these cryptosystems efficiently. Among those operations, as being the most time-consuming, multiplicative inversion has become the object of lots of investigation Formant's theorem says $\beta$$^{-1}$ =$\beta$$^{2}$sup m/-2/, where $\beta$$^{-1}$ is the multiplicative inverse of $\beta$$\in$GF(2$^{m}$ ). Therefore, to compute the multiplicative inverse of arbitrary elements of GF(2$^{m}$ ), it is most important to reduce the number of times of multiplication by decomposing 2$^{m}$ -2 efficiently. Among many algorithms relevant to the subject, the algorithm proposed by Itoh and Tsujii[2] has reduced the required number of times of multiplication to O(log m) by using normal basis. Furthermore, a few papers have presented algorithms improving the Itoh and Tsujii's. However they have some demerits such as complicated decomposition processes[3,5]. In this paper, in the case of 2$^{m}$ -2, which is mainly used in practical applications, an efficient algorithm is proposed for computing the multiplicative inverse at high speed by using both the factorization formula x$^3$-y$^3$=(x-y)(x$^2$＋xy＋y$^2$) and normal basis. The number of times of multiplication of the algorithm is smaller than that of the algorithm proposed by Itoh and Tsujii. Also the algorithm decomposes 2$^{m}$ -2 more simply than other proposed algorithms.