• Journal of the Korean Geotechnical Society
• /
• v.39 no.7
• /
• pp.57-67
• /
• 2023

Slope failures during rainfall have been observed in mountainous areas of South Korea as a result of the presence of solar power facilities. The seepage behavior and pore pressure distribution differ from typical slopes due to the presence of impermeable solar panels, and the load imposed by the solar power structures also affects the slope behavior. This study aims to develop a method for evaluating the stability of slopes with solar power facilities and to analyze vulnerable points by considering the maximum slope displacement. To assess the slope stability and predict behavior while considering rainfall seepage, a combined seepage analysis and finite difference method numerical analysis were employed. For the selected site, various variables were assumed, including parameters related to the Soil Water Characteristic Curve, strength parameters that satisfy the Mohr-Coulomb failure criterion, soil properties, and topographic factors such as slope angle and bedrock depth. The factors with the most significant influence on the factor of safety (FOS) were identified. The presence of solar power facilities was found to affect the seepage distribution and FOS, resulting in a decreasing trend due to rainfall seepage. The maximum displacement points were concentrated near the upper (crest) and lower (toe) sections of the slope.

• Journal of Dental Rehabilitation and Applied Science
• /
• v.17 no.4
• /
• pp.283-305
• /
• 2001

The purpose of this study was to analyze the stress distribution of condylar regions and edentulous mandible with implant-supported cantilever prostheses on the certain conditions, such as amount of load, location of load, direction of load, fixation or non-fixation on the condylar regions. Three dimensional finite element analysis was used for this study. FEM model was created by using commercial software, ANSYS(Swanson, Inc., U.S.A.). Fixed model which was fixed on the condylar regions was modeled with 74323 elements and 15387 nodes and spring model which was sprung on the condylar regions was modeled with 75020 elements and 15887 nodes. Six Br${\aa}$nemark implants with 3.75 mm diameter and 13 mm length were incorporated in the models. The placement was 4.4 mm from the midline for the first implant; the other two in each quardrant were 6.5 mm apart. The stress distribution on each model through the designed mandible was evaluated under 500N vertical load, 250N horizontal load linguobuccally, buccal 20 degree 250N oblique load and buccal 45 degree 250N oblique load. The load points were at 0 mm, 10 mm, 20 mm along the cantilever prostheses from the center of the distal fixture. The results were as follows; 1. The stress distribution of condylar regions between two models showed conspicuous differences. Fixed model showed conspicuous stress concentration on the condylar regions than spring model under vertical load only. On the other hand, spring model showed conspicuous stress concentration on the condylar regions than fixed model under 250N horizontal load linguobuccally, buccal 20 degree 250N oblique load and buccal 45 degree 250N oblique load. 2. Fixed model showed stress concentration on the posterior and mesial side of working and balancing condylar necks but spring model showed stress concentration on the posterior and mesial side of working condylar neck and the posterior and lateral side of balancing condylar neck under vertical load. 3. Fixed model showed stress concentration on the posterior and lateral side of working condylar neck and the anterior and mesial side of balancing condylar neck but spring model showed stress concentration on the anterior sides of working and balancing condylar necks under horizontal load linguobuccally. 4. Fixed model showed stress concentration on the posterior side of working condylar neck and the posterior and lateral side of balancing condylar neck but spring model showed stress concentration on the anterior side of working condylar neck and the anterior and lateral side of balancing condylar neck under buccal 20 degree oblique load. 5. Fixed model showed stress concentration on the anterior and lateral side of working condylar neck and the posterior and mesial side of balancing condylar neck but spring model showed stress concentration on the anterior side of working condylar neck and the anterior and lateral side of balancing condylar neck under buccal 45 degree oblique load.. 6. The stress distribution of bone around implants between two models revealed difference slightly. In general, magnitude of Von Mises stress was the greatest at the bone around the most distal implant and the progressive decrease more and more mesially. Under vertical load, the stress values were similar between implant neck and superstructure vertically, besides the greatest on the distal side horizontally. 7. Under horizontal load linguobuccally, buccal 20 degree oblique load and buccal 45 degree oblique load, the stress values were the greatest on the implant neck vertically, and great on the labial and lingual sides horizontally. After all, it was considered that spring model was an indispensable condition for the comprehension of the stress distributions of condylar regions.

• The Journal of Korean Academy of Prosthodontics
• /
• v.46 no.5
• /
• pp.479-489
• /
• 2008

Statement of the problem: Under anatomical limitations on maxillary posterior region, a poor crown-to root ratio acting on dental implants can result in undesirable stress in surrounding bone, which in turn can cause bone defects and eventual failure of implants. Purpose: The purpose is to compare stress distribution due to different crown-root ratio and effect of splinting between natural teeth and implants in maxillary molar area under different loads. Material and methods: Analysis of stress arising supporting bone of the natural teeth and the implant was made with 3-dimensional finite element method. The model simulated naturel teeth was made with 2nd premolar and 1st molar in the maxillary molar region (Model T). The model simulated implants placed on same positions with two parallel implants of Straumann Dental Implant cemented abutment (Model I). Each model was designed in different crown-root ratio (0.7:1, 1:1, 1.25:1) and set cement type gold crown to make it non-splinted or splinted. After that, 300 N force was loaded to each model in five ways (Load 1: middle of occlusal table, Load 2: middle of buccal cusp, Load 3: middle of lingual cusp, Load 4: horizontal load to buccal cusp of anterior abutment only, Load 5: horizontal load to middle of buccal cusp of each abutment), and stress distribution was analyzed. Results and conclusion: On all occasions, stress was concentrated at the cervical region of the implant. Under load 1, 2 and 3, stress was not increased even when crown-root ratio increases, but under load 4 and 5, when crown-root ratio increases, stress also increased. There was difference in stress values between natural teeth and implants when crown-root ratio gradually increases; In case of natural teeth, splinting decreased stress under vertical and horizontal loads. In case of implants, splinting decreased stress under vertical loads 1,2 and 3, but increased maximal stress under loads 2 and 3. Under horizontal loads, splinting decreased stress, however the effect of splinting decreased under load 5 than load 4. Furthermore, the stress was increased, when crown-root ratio is 1.25:1. Clinical implications: This limited finite element study suggests that the stress on supporting bone may be increased under non-axial loads and poor crown-root ratio. Under poor crown-root ratio, excessive stress was generated at the cervical region of the implant, and decreased splinting effect for stress distribution, which can be related to clinical failure.

• Journal of Dental Rehabilitation and Applied Science
• /
• v.25 no.4
• /
• pp.375-390
• /
• 2009

Statement of problem: Implant supported overdenture is accepted widely as a way to restore edentulous ridge providing better retention and support of dentures. Various types of attachment for overdenture have been developed. Purpose: The purpose of this study was to investigate the influence of attachment type in implant overdentures on the biomechanical stress distribution in the surrounding bone, prosthesis and interface between implant and bone. Material and methods: Finite element analysis method was used. Average CT image of mandibular body(Digital $Korea^{(R)}$, KISTI, Korea) was used to produce a mandibular model. Overdentures were placed instead of mandibular teeth and 2mm of mucosa was inserted between the overdenture and mandible. Two implants($USII^{(R)}$, Osstem, Korea) were placed at both cuspid area and 4 types of overdenture were fabricated ; ball and socket, Locator, magnet and bar type. Load was applied on the from second premolar to second molar tooth area. 6 times of finite element analyses were performed according to the direction of the force $90^{\circ}$, $45^{\circ}$, $0^{\circ}$ and unilateral or bilateral force applied. The stress at interface between implants and bone, and prosthesis and the bone around implants ware compared using von Mises stress. The results were explained with color coded graphs based on the equivalent stress to distinguish the force distribution pattern and the site of maximum stress concentration. Results: Unilateral loading showed that connection area between implant fixture and bar generated maximum stress in bar type overdentures. Bar type produced 100 Mpa which means the most among 4 types of attachments. Bilateral loading, however, showed that bar type was more stable than other implants(magnet, ball and socket). 26 Mpa of bar type was about a half of other types on overdenture under $90^{\circ}$ bilateral loading. Conclusions: In any directions of stress, bar type was proved to be the most vulnerable type in both implants and overdentures. Interface stress did not show any significant difference in stress distribution pattern.

• 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 the Korean Society of Hazard Mitigation
• /
• v.3 no.3 s.10
• /
• pp.151-163
• /
• 2003

In this study, the algorithm of groundwater flow process was established for koreanized groundwater program development dealing with the geographic and geologic conditions of the aquifer have dynamic behaviour in groundwater flow system. All the input data settings of the 3-DFM model which is developed in this study are organized in Korean, and the model contains help function for each input data. Thus, it is designed to get detailed information about each input parameter when the mouse pointer is placed on the corresponding input parameter. This model also is designed to easily specify the geologic boundary condition for each stratum or initial head data in the work sheet. In addition, this model is designed to display boxes for input parameter writing for each analysis condition so that the setting for each parameter is not so complicated as existing MODFLOW is when steady and unsteady flow analysis are performed as well as the analysis for the characteristics of each stratum. Descriptions for input data are displayed on the right side of the window while the analysis results are displayed on the left side as well as the TXT file for this results is available to see. The model developed in this study is a numerical model using finite differential method, and the applicability of the model was examined by comparing and analyzing observed and simulated groundwater heads computed by the application of real recharge amount and the estimation of parameters. The 3-DFM model is applied in this study to Sehwa-ri, and Songdang-ri area, Jeju, Korea for analysis of groundwater flow system according to pumping, and obtained the results that the observed and computed groundwater head were almost in accordance with each other showing the range of 0.03 - 0.07 error percent. It is analyzed that the groundwater flow distributed evenly from Nopen-orum and Munseogi-orum to Wolang-bong, Yongnuni-orum, and Songja-bong through the computation of equipotentials and velocity vector using the analysis result of simulation which was performed before the pumping started in the study area. These analysis results show the accordance with MODFLOW's.

• Geophysics and Geophysical Exploration
• /
• v.11 no.4
• /
• pp.318-325
• /
• 2008

Assigning an exact boundary condition is of great importance in three-dimensional (3D) magnetotelluric (MT) modeling, in which no source is considered in a computing domain. This paper presents a 3D MT modeling algorithm utilizing a Dirichlet condition for a 2D host. To compute boundary values for a model with a 2D host, we need to conduct additional 2D MT modeling. The 2D modeling consists of transverse magnetic and electric modes, which are determined from the relationship between the polarization of plane wave and the strike direction of the 2D structure. Since the 3D MT modeling algorithm solves Maxwell's equations for electric fields using the finite difference method with a staggered grid that defines electric fields along cell edges, electric fields are calculated at the same place in the 2D modeling. The algorithm developed in this study can produce reliable MT responses for a 3D model with a 2D host.

• The Journal of Engineering Geology
• /
• v.24 no.2
• /
• pp.261-271
• /
• 2014

An almost permanent anchor (friction type) is resistant to ground deformation due to the friction between the soil and grout at a fixed length from the anchor body. The purpose of this study is to calculate the force of bearing resistance for a bearing anchor in enlarged boreholes. We conducted analytical and numerical analyses, along with laboratory testing, to find the quantities of bearing resistance prior to grouting in EBA (Enlarged Bearing Anchor) construction. The force of bearing resistance from the analytical method was defined as a function of general borehole diameter, expanded borehole diameter, and soil unconfined compressive strength. We also employed the Flac 3D finite difference numerical modeling code to analyze the bearing resistance of the soil conditions. We then created a laboratory experimental model to measure bearing resistance and carried out a pull-out test. The results of these three analyses are presented here, and a regression analysis was performed between bearing resistance and uniaxial compression strength. The laboratory results yield the strongest bearing resistance, with reinforcement 28.5 times greater than the uniaxial compression strength; the analytical and numerical analyses yielded values of 13.3 and 9.9, respectively. This results means that bearing resistance of laboratory test appears to be affected by skin friction resistance. To improve the reliability of these results, a comparison field study is needed to verify which results (analytical, numerical, or laboratory) best represent field observations.

• Journal of the Korean Geotechnical Society
• /
• v.30 no.1
• /
• pp.49-63
• /
• 2014

Bearing capacity factor $N_{\gamma}$ and shape factor were studied for rigid strip and circular footings with a rough base on sand by numerical modelling considering the effect of dilation angle. The numerical model was developed with an explicit finite difference code. Loading procedures and interpretation methods were devised in order to shorten the running time while eliminating the exaggeration of the reaction caused by the explicit scheme. Using the Mohr-Coulomb plasticity model with associated (${\psi}={\phi}$) and nonassociated (${\psi}$ < ${\phi}$) flow-rules, the bearing capacity factor $N_{\gamma}$ was evaluated for various combinations of internal friction angles and dilation angles. Bearing capacity factor decreased as the dilation angle was reduced from the associated condition. An equation applicable to typical sands was proposed to evaluate the relative bearing capacity for the nonassociated condition compared to the associated condition on which most bearing capacity factor equations are based. The shape factor for the circular footing varied substantially when the plane-strain effect was taken into account for the strip footing. The numerical results of this study showed closer trends with the previous experimental results when the internal friction angle was increased for the strip footing. Discussions are made on the reason that previous equations for the shape factor give different results and recommendations are made for the appropriate design shape factor.

• Proceedings of the Korean Vacuum Society Conference
• /
• 2012.08a
• /
• pp.184-184
• /
• 2012

Recently, InGaN/GaN multi-quantum well grown on GaN pyramid structures have attracted much attention due to their hybrid characteristics of quantum well, quantum wire, and quantum dot. This gives us broad band emission which will be useful for phosphor-free white light emitting diode. On the other hand, by using quantum dot emission on top of the pyramid, site selective single photon source could be realized. However, these structures still have several limitations for the single photon source. For instance, the quantum efficiency of quantum dot emission should be improved further. As detection systems have limited numerical aperture, collection efficiency is also important issue. It has been known that micro-cavities can be utilized to modify the radiative decay rate and to control the radiation pattern of quantum dot. Researchers have also been interested in nano-cavities using localized surface plasmon. Although the plasmonic cavities have small quality factor due to high loss of metal, it could have small mode volume because plasmonic wavelength is much smaller than the wavelength in the dielectric cavities. In this work, we used localized surface plasmon to improve efficiency of InGaN qunatum dot as a single photon emitter. We could easily get the localized surface plasmon mode after deposit the metal thin film because lnGaN/GaN multi quantum well has the pyramidal geometry. With numerical simulation (i.e., Finite Difference Time Domain method), we observed highly enhanced decay rate and modified radiation pattern. To confirm these localized surface plasmon effect experimentally, we deposited metal thin films on InGaN/GaN pyramid structures using e-beam deposition. Then, photoluminescence and time-resolved photoluminescence were carried out to measure the improvement of radiative decay rate (Purcell factor). By carrying out cathodoluminescence (CL) experiments, spatial-resolved CL images could also be obtained. As we mentioned before, collection efficiency is also important issue to make an efficient single photon emitter. To confirm the radiation pattern of quantum dot, Fourier optics system was used to capture the angular property of emission. We believe that highly focused localized surface plasmon around site-selective InGaN quantum dot could be a feasible single photon emitter.