• Journal of the Institute of Electronics Engineers of Korea SD
• /
• v.41 no.6
• /
• pp.25-31
• /
• 2004

In this paper, we introduce kinetic Monte Carlo (kMC) methods for simulating diffusion process in nano-scale device fabrication. At first, we review kMC theory and backgrounds and give a simple point defect diffusion process modeling in thermal annealing after ion (electron) implantation into Si crystalline substrate to help understand kinetic Monte Carlo methods. kMC is a kind of Monte Carlo but can simulate time evolution of diffusion process through Poisson probabilistic process. In kMC diffusion process, instead of. solving differential reaction-diffusion equations via conventional finite difference or element methods, it is based on a series of chemical reaction (between atoms and/or defects) or diffusion events according to event rates of all possible events. Every event has its own event rate and time evolution of semiconductor diffusion process is directly simulated. Those event rates can be derived either directly from molecular dynamics (MD) or first-principles (ab-initio) calculations, or from experimental data.

• Geophysics and Geophysical Exploration
• /
• v.5 no.2
• /
• pp.84-92
• /
• 2002

In this study, the new modeling scheme has been developed for recently designed and tested electromagnetic survey, which adapts horizontal magnetic dipole with $1\;kHz\~1\;MHz$ frequency range as a source. The 2.5-D secondary field formulation in wavenumber domain was constructed using finite element method and verified through comparing results with layered-earth solutions calculated by integral equations. 2-D conductive- and resistive-block models were constructed for calculating electric field, magnetic field and impedance - the ratio of electric and magnetic fields which are orthogonal each other. This study showed that electric field and impedance are superior in identifying 2-D isolated-body model to magnetic field. In particular, impedance gives more stable results than electric field with similar spatial resolving power, because electric field is divided by magnetic field in impedance. Thus the impedance analysis which uses electric and magnetic fields together would give better result in imaging the shallow anomalies than conventional EM method.

• Journal of the Korean GEO-environmental Society
• /
• v.9 no.2
• /
• pp.39-46
• /
• 2008

Due to the growth in industrialization, potential hazards in subsurface environments are becoming increasingly significant. The extraction of the contaminant from the soil and movement of the water are restricted due to the low permeability and adsorption characteristics of the reclaimed soils. There are a number of approaches to in-situ remediation that are used in contaminated sites for removing contaminants. These include soil flushing, dual phase extraction, and soil vapor extraction. Among these techniques, soil flushing was the focus of the investigation in this paper. Incorporated technique with PVDs has been used for dewatering from fine-grained soils for the purpose of ground improvement by means of prefabricated vertical drain systems. The laboratory model tests were performed by using the flushing tracer solutions for silty soils and recorded the tracer concentration changes with the elapsed time and flow rates. The modeling was intended to predict the effectiveness and time dependence of the remediation process. Modeling has been performed on the extraction, considering tracer concentration and laboratory model test characteristics. The computer model used herein are SEEP/W and CTRAN/W, this 2-D finite element program allows for modeling to determine hydraulic head and pore water pressure distribution, efficiency of remediation for the subsurface environment. It is concluded that the coefficient of permeability of contaminated soil is related with vertical velocity and extracted flow rate. The vertical velocity and extracted flow rate have an effect on dispersivity and finally are played an important role in-situ soil remediation.

• Geophysics and Geophysical Exploration
• /
• v.24 no.2
• /
• pp.35-44
• /
• 2021

As the execution speed of Python is slower than those of other programming languages (e.g., C, C++, and FORTRAN), Python is not considered to be efficient for writing numerical geodynamic code that requires numerous iterations. Recently, many computational techniques, such as the Just-In-Time (JIT) compiler, have been developed to enhance the calculation speed of Python. Here, we developed two-dimensional (2D) numerical geodynamic code that was optimized for the JIT compiler, based on Python. Our code simulates mantle convection by combining the Particle-In-Cell (PIC) scheme and the finite element method (FEM), which are both commonly used in geodynamic modeling. We benchmarked well-known mantle convection problems to evaluate the reliability of our code, which confirmed that the root mean square velocity and Nusselt number obtained from our numerical modeling were consistent with those of the mantle convection problems. The matrix assembly and PIC processes in our code, when run with the JIT compiler, successfully achieved a speed-up 30× and 258× faster than without the JIT compiler, respectively. Our Python-based FEM-PIC code shows the high potential of Python for geodynamic modeling cases that require complex computations.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.28 no.7
• /
• pp.1259-1266
• /
• 2022

Recently, due to the specialization of structural design standards and evaluation methods, the classification rules are being integrated. A good example is the common international rules (CSR). However, detailed regulations are presented only for the cargo hold area where the longitudinal load is greatly applied, and no specific evaluation guidelines exist for the bow and stern structures. Structural design of the mentioned area is carried out depending on the design experience of the shipbuilder, and because no clear standard exists even in the classification, determining the root cause is difficult even if a structural damage problem occurs. In this study, an engineering-based solution was presented to identify the root cause of representative cases of buckling damage that occurs mainly in the stern. Buckling may occur at the panel wall owing to hull girder bending moment acting on the stern structure, and the plate thickness must be increased or vertical stiffeners must be added to increase the buckling rigidity. For structural strength verification based on finite element analysis modeling, reasonable solutions for load conditions, boundary conditions, modeling methods, and evaluation criteria were presented. This result is expected to be helpful in examining the structural strength of the stern part of similar carriers in the future.

• Computers and Concrete
• /
• v.31 no.4
• /
• pp.337-348
• /
• 2023

Shear wall is commonly used as a lateral force resisting system of concrete mid-rise and high-rise buildings, but it brings challenges in providing relatively large space throughout the building height. For this reason, the structure system where the upper structure with bearing, non-bearing and/or shear walls that sits on top of a transfer plate system supported by widely spaced columns at the lower stories is preferred in some regions, particularly in low to moderate seismic regions in Asia. A thick reinforced concrete (RC) plate has often been used as a transfer system, along with RC transfer girders; however, the RC plate becomes very thick for tall buildings. Applying the post-tensioning (PT) technique to RC plates can effectively reduce the thickness and reinforcement as an economical design method. Currently, a simplified model is used for numerical modeling of PT transfer plate, which does not consider the interaction of the plate and the upper structure. To observe the actual behavior of PT transfer plate under seismic loads, it is necessary to model whole parts of the structure and tendons to precisely include the interaction and the secondary effect of PT tendons in the results. This research evaluated the seismic behavior of shear wall-type residential buildings with PT transfer plates for the condition that PT tendons are included or excluded in the modeling. Three-dimensional finite element models were developed, which includes prestressing tendon elements, and response spectrum analyses were carried out to evaluate seismic forces. Two buildings with flat-shape and L-shape plans were considered, and design forces of shear walls and transfer columns for a system with and without PT tendons were compared. The results showed that, in some cases, excluding PT tendons from the model leads to an unrealistic estimation of the demands for shear walls sit on transfer plate and transfer columns due to excluding the secondary effect of PT tendons. Based on the results, generally, the secondary effect reduces shear force demand and axial-flexural demands of transfer columns but increases the shear force demand of shear walls. The results of this study suggested that, in addition to the effect of PT on the resistance of transfer plate, it is necessary to include PT tendons in the modeling to consider its effect on force demand.

• Journal of Korea Water Resources Association
• /
• v.41 no.11
• /
• pp.1123-1133
• /
• 2008

It needs to conceptualize watershed with triangular or rectangular elements and to analyze the changes in hydrological components of each element for distributed modeling of rainfall-runoff process. This study is the network analysis of watershed grid for flow routing occurred in each element when analyzing rainfall-runoff process by one-dimensional kinematic wave equation. Single flow direction from D8-method(deterministic eight-neighbors method) is used, and the information of flow direction and flow accumulation are used to determine the computation order of each element. The application theory of finite volume method is suggested for each flow direction pattern between elements, and it is applied it to calculate the flow of each grid. Network analysis method from this study is applied to GRM(Grid based Rainfall-runoff Model) which is physically based distributed rainfall-runoff model, and the results from simplified hypothetical watersheds are compared with $Vflo^{TM}$ to examine the reasonability of the method. It is applied to Jungrangcheon watershed in Han river for verification, and examination of the applicability to real site. The results from Jungrangcheon watershed show good agreement with measured hydrographs, and the application of the network analysis method to real site is proper.

• Journal of Korean Society of Steel Construction
• /
• v.25 no.3
• /
• pp.307-320
• /
• 2013

This paper proposes a new hybrid support structure by the driven piles which removes disadvantages of the existing type of support structure for offshore wind turbines. The hybrid type of support structure is combined with concrete cone and steel shaft, and is supported not only by gravity type foundations but also by driven piles. For three dimensional analysis of the huge and thick concrete structure, a solid-shell element that is capable of exact modeling and node interpolations of stresses is developed. By applying wave theory of stream function and solid-shell element in XSEA simulation software for fixed offshore wind turbines, a quasi-static analysis and natural frequency analysis of proposed support structure are performed with the environmental condition on Southwest Coast in Korea. In the result, lateral displacement is not exceed allowable displacement and a superiority of dynamic behavior of new hybrid support structure is validated by natural frequency analysis. Consequently, the hybrid support structure presented in this study has a structural stability enough to be applied on real-site condition in Korea. The optimized structures based on the preliminary design concept resulted in an efficient structure, which reasonably reduces fabrication costs.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.6 no.2
• /
• pp.459-470
• /
• 2014

During ship design, welding-induced distortions are roughly estimated as a function of the size of the component as well as the welding process and residual stresses are assumed to be locally in the range of the yield stress. Existing welding simulation methods are very complex and time-consuming and therefore not applicable to large structures like ships. Simplified methods for the estimation of welding effects were and still are subject of several research projects, but mostly concerning smaller structures. The main goal of this paper is the application of a multi-layer welding simulation to the block joint of a ship structure. When welding block joints, high constraints occur due to the ship structure which are assumed to result in accordingly high residual stresses. Constraints measured during construction were realized in a test plant for small-scale welding specimens in order to investigate their and other effects on the residual stresses. Associated welding simulations were successfully performed with fine-mesh finite element models. Further analyses showed that a courser mesh was also able to reproduce the welding-induced reaction forces and hence the residual stresses after some calibration. Based on the coarse modeling it was possible to perform the welding simulation at a block joint in order to investigate the influence of the resulting residual stresses on the behavior of the real structure, showing quite interesting stress distributions. Finally it is discussed whether smaller and idealized models of definite areas of the block joint can be used to achieve the same results offering possibilities to consider residual stresses in the design process.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.10 no.5 s.51
• /
• pp.49-61
• /
• 2006

The floor impact sound insulations are installed frequently for reducing the floor impact sound into the floor slab of the residential buildings in recent years. Therefore the analytical FE model considering the insulation is needed for the sound and vibration analysis of the floor and it is necessary to use a refined finite element model fer considering the large number of modes involved dynamic responses. So it is very difficult to use FE model because of the tiresome task for constructing the FE model, taking a lot of times for analysis and the impossibility of using the proportional damping. The efficient analysis and modeling method are proposed to the dynamic analysis for the floor with damping layer in this study. The floor slabs and finished layers are modeled individually and the spring elements that mean damping layers used to connect two parts. The dynamic analysis by the $Newmark-{\beta}$ method is performed to solve the non-proportional damping problem due to the damping coefficient of insulations.