• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 2006.04a
• /
• pp.87-91
• /
• 2006

본 연구의 대상 지역은 전라북도 부안군 변산면 격포리 일대이다. 이 지역의 해안 대수층 내에서 발생하고 있는 해수 침투 현상을 효과적으로 모사하고 지하수 양수정의 설치 및 운영의 영향을 정량적으로 평가하기 위하여 지형과 지층의 구성 등의 수리지질학적 특성을 종합적으로 고려한 지하수 유동 및 염분 이동 삼차원 수치 모델링을 수행하였다. 지하수 양수 전의 정상 상태 수치 모델링 결과에서는 지하수 유동 및 염분 이동의 측면에서 대상지역의 수리지질학적 특성이 잘 반영되어 나타났다. 또한 지하수 양수 시의 비정상 상태 수치 모델링 결과에서는 지형과 지층의 특성에 따라서 염분 이동의 양상이 특징적으로 나타났으며 양수정의 위치에 따른 염분 농도의 시간적 공간적 분포의 변화를 효과적으로 예측할 수 있었다. 이러한 결과들은 향후 연구지역 뿐만 아니라 국내의 다른 해안 지역에서 지하수 자원의 개발 시에 해수 침투의 양상을 파악하고 예측하는데 효과적으로 활용될 수 있을 것으로 기대된다.

• Proceedings of the Korean Radioactive Waste Society Conference
• /
• 2009.06a
• /
• pp.184-185
• /
• 2009

• Journal of Soil and Groundwater Environment
• /
• v.14 no.6
• /
• pp.1-18
• /
• 2009

A series of three-dimensional numerical simulations using a generalized multidimensional hydrodynamic dispersion numerical model is performed to simulate effectively and to evaluate quantitatively impacts of urbanization on density-dependent groundwater flow and salt transport in a coastal aquifer system, Suyeong-Gu, Busan, Korea. A series of steady-state numerical simulations of groundwater flow and salt transport before urbanization with material properties of geologic formations, which are established by numerical modeling calibrations considering all the urbanization factors, is performed first without considering all the urbanization factors. A series of transient-state numerical simulations of groundwater flow and salt transport after urbanization is then performed considering the urbanization factors individually and all together. Finally, the results of both numerical simulations are compared with each other and analyzed. The results of the numerical simulations show that density-dependent groundwater flow, salt transport, and seawater intrusion in the coastal aquifer system are intensively and extensively impacted by the urbanization factors. Especially, these urbanization factors result in the changes of the total groundwater volume and salt mass in the coastal aquifer system. However, such impacts of each urbanization factor are not spatially uniform but locally different.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.4
• /
• pp.469-476
• /
• 2007

A domain/boundary decomposition technique is applied to carry out efficient finite element analyses of 3-D contact problems. Appropriate penalty functions are selected for connecting an interface and contact interfaces with neighboring subdomains that satisfy continuity constraints. As a consequence, all the effective stiffness matrices have positive definiteness, and computational efficiency can be improved to a considerable degree. If necessary, any complex-shaped 3-D domain can be divided into several simple-shaped subdomains without considering the conformity of meshes along the interface. With a set of numerical examples, the basic characteristics of computational efficiency are investigated carefully.

• Journal of Soil and Groundwater Environment
• /
• v.13 no.5
• /
• pp.33-46
• /
• 2008

A series of three-dimensional numerical simulations using a generalized multidimensional hydrodynamic dispersion numerical model is performed to simulate effectively and to evaluate quantitatively impacts of fault existence on densitydependent groundwater flow and salt transport in coastal aquifer systems. A series of steady-state numerical simulations with calibration is performed first for an actual coastal aquifer system which contains a major fault. A series of steadystate numerical simulations is then performed for a corresponding coastal aquifer system which does not have such a major fault. Finally, the results of both numerical simulations are compared with each other and analyzed. The results of the numerical simulations show that the major fault produces hydrogeologically significant heterogeneity and true anisotropy in the actual coastal aquifer system, and density-dependent groundwater flow, salt transport, and seawater intrusion patterns in the coastal aquifer systems are intensively and extensively dependent upon the existence or absence of such a major fault. Especially, the major fault may act as a pathway for groundwater flow and salt transport along the direction parallel to its plane, while it may also behave as a barrier against groundwater flow and salt transport along the direction perpendicular to its plane.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.15 no.10 s.89
• /
• pp.961-968
• /
• 2004

In this paper, we present a set of numerical schemes to solve the Muller integral equation for the analysis of electromagnetic scattering from arbitrarily shaped three-dimensional dielectric bodies by applying the method of moments(Mon. The piecewise homogeneous dielectric structure is approximated by planar triangular patches. A set of the RWG(Rao, Wilton, Glisson) functions is used for expansion of the equivalent electric and magnetic current densities and a combination of the RWG function and its orthogonal component is used for testing. Numerical results for a dielectric sphere are presented and compared with solutions obtained using other formulations.

• Journal of Soil and Groundwater Environment
• /
• v.13 no.4
• /
• pp.8-21
• /
• 2008

A series of three-dimensional numerical simulations using a hydrodynamic dispersion numerical model is performed to analyze quantitatively impacts of layered heterogeneity of geologic media and groundwater pumping schemes on groundwater flow and salt transport in coastal aquifer systems. A two-layer heterogeneous coastal aquifer system composed of a lower sand layer (aquifer) and an upper clay layer (aquitard) and a corresponding single-layer homogeneous coastal aquifer system composed of an equivalent lumped material are simulated to evaluate impacts of layered heterogeneity on seawater intrusion. In addition, a continuous groundwater pumping scheme and two different periodical groundwater pumping schemes, which withdraw the same amount of groundwater during the total simulation time, are applied to the above two coastal aquifer systems to evaluate impacts of groundwater pumping schemes on seawater intrusion. The results of the numerical simulations show that the periodical groundwater pumping schemes have more significant adverse influences on groundwater flow and salt transport not only in the lower sand layer but also in the upper clay layer, and groundwater salinization becomes more intensified spatially and temporally as the pumping intensity is higher under the periodical groundwater pumping schemes. These imply that the continuous groundwater pumping scheme may be more suitable to minimize groundwater salinization due to seawater intrusion. The results of the numerical simulations also show that groundwater salinization in the upper clay layer occurs significantly different from that in the lower sand layer under the periodical groundwater pumping schemes. Such differences in groundwater salinization between the two adjacent layers may result from layered heterogeneity of the layered coastal aquifer system.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.27 no.4B
• /
• pp.379-387
• /
• 2002

In this paper, we present a procedure to obtain the transient scattering response from three-dimensional arbitrarily shaped and closed conducting bodies using time-domain magnetic field integral equation (TD-MFIE) with triangular patch functions. This approach results in accurate and comparably stable transient responses from conducting scatterers. Detailed mathematical steps are included, and several numerical results are presented and compared with results from a time-domain electric field integral equation (TD-EFIE) and the inverse courier transform solution of the frequency domain results.

• Journal of the Korean Geotechnical Society
• /
• v.30 no.3
• /
• pp.29-46
• /
• 2014

For circular rigid footings with a rough base on sand, combined vertical - horizontal loading capacity was studied by three-dimensional numerical modelling. A numerical model was implemented to simulate the swipe loading and the probe loading methods and an interpretation procedure was devised in order to eliminate the numerical error from the restricted mesh density. Using the Mohr-Coulomb plasticity model, the effect of friction angle was studied under the associated flow-rule condition. The swipe loading method, which is efficient in that the interaction diagram can be drawn with smaller number of analyses, was confirmed to give similar results with the probe loading method, which follows closely the load-paths applied to real structures. For circular footings with a rough base, the interaction diagram for combined vertical (V) - horizontal (H) loading and the inclination factor were barely affected by the friction angle. It was found that the inclination factors for strip and rectangular footings are applicable to circular footings. For high H/V ratios, the results by numerical modelling of this study were smaller than the results of previous studies. Discussions are made on the factors affecting the numerical results and the areas for further researches.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.35 no.9
• /
• pp.1041-1050
• /
• 2011

An elasto-plastic finite element method using the theory of strain gradient plasticity is proposed to evaluate the size dependency of structural plasticity that occurs when the configuration size decreases to micron scale. For this method, we suggest a low-order plane and three-dimensional displacement-based elements, eliminating the need for a high order, many degrees of freedom, a mixed element, or super elements, which have been considered necessary in previous researches. The proposed method can be performed in the framework of nonlinear incremental analysis in which plastic strains are calculated and averaged at nodes. These strains are then interpolated and differentiated for gradient calculation. We adopted a strain-gradient-hardening constitutive equation from the Taylor dislocation model, which requires the plastic strain gradient. The developed finite elements are tested numerically on the basis of typical size-effect problems such as micro-bending, micro-torsion, and micro-voids. With respect to the strain gradient plasticity, i.e., the size effects, the results obtained by using the proposed method, which are simple in their calculation, are in good agreement with the experimental results cited in previously published papers.