• Journal of Korea Water Resources Association
• /
• v.42 no.3
• /
• pp.247-257
• /
• 2009

In this study, a distributed rainfall-runoff model is developed using a multi-directional flow allocation algorithm and the real-time runoff updating algorithm. The developed model consists of relatively simple governing equations of hydrologic processes in order to apply developed algorithms and to enhance the efficiency of computational time which is drawback of distributed model application. The variability of topographic characteristics and flow direction according to various spatial resolution were analyzed using DEM(Digital Elevation Model) data. As a preliminary process using fine resolution DEM data, a multi-directional flow allocation algorithm was developed to maintain detail flow information in distributed rainfall-runoff simulation which has strong advantage in computation efficiency and accuracy. Also, a real-time updating algorithm was developed to update current watershed condition. The developed model is able to hold the information of actual behavior of runoff process in low resolution simulation. Therefore it is expected the improvement of forecasting accuracy and computational efficiency.

• Journal of Korean Society for Geospatial Information Science
• /
• v.14 no.1 s.35
• /
• pp.49-55
• /
• 2006

Nowadays, flood hazard maps have been generated to minimize the damages from the flooding. To generate such flood hazard maps, LiDAR data can be used as data source with higher data accuracy. LiDAR data, however, requires relatively higher cost and longer processing time. In this background, this study proposed DEM generation using NGIS digital topographic maps. For that, breaklines were processed to count directions of water flows. In addition, the river profile data, unique data source to represent real topography of the river area, were integrated to the breaklines to generate DEM. City of Kuri in Kyunggi Province was selected for this study and 1:1,000 and 1:5,000 topographic maps were integrated to process breaklines and river profile data were also linked to generate DEM. The generated DEM showed relatively lower vertical accuracy from mixing 1:1,000 and 1:5,000 topographic maps since 1:1,000 topographic maps were not available for some portion of the area. However, the DEM generated demonstrated reasonable accuracy and resolution for flood map generation as well as higher cost saving effects. On the contrary, for more efficient utilization of NGIS topographic maps, periodic map updating needs to be made including technical consideration in building breaklines and applying interpolation methods.

• Computers and Concrete
• /
• v.5 no.4
• /
• pp.343-358
• /
• 2008

This article focuses on concrete structures submitted to impact loading and is aimed at predicting local damage in the vicinity of an impact zone as well as the global response of the structure. The Discrete Element Method (DEM) seems particularly well suited in this context for modeling fractures. An identification process of DEM material parameters from macroscopic data (Young's modulus, compressive and tensile strength, fracture energy, etc.) will first be presented for the purpose of enhancing reproducibility and reliability of the simulation results with DE samples of various sizes. The modeling of a large structure by means of DEM may lead to prohibitive computation times. A refined discretization becomes required in the vicinity of the impact, while the structure may be modeled using a coarse FE mesh further from the impact area, where the material behaves elastically. A coupled discrete-finite element approach is thus proposed: the impact zone is modeled by means of DE and elastic FE are used on the rest of the structure. The proposed approach is then applied to a rock impact on a concrete slab in order to validate the coupled method and compare computation times.

• Geomechanics and Engineering
• /
• v.35 no.5
• /
• pp.475-486
• /
• 2023

Ballast particles have an irregular shape and are discrete in nature. Due to the discrete nature of ballast, it exhibits complex mechanical behaviour under loading conditions. The discrete element method (DEM) can model the behaviour of discrete particles under a multitude of loading conditions. DEM is used in this paper to simulate a series of three-dimensional direct shear tests in order to investigate the shear behaviour of railway ballast and its interaction at the microscopic level. Particle flow code in three dimension (PFC3D) models the irregular shape of ballast particles as clump particles. To investigate the influence of particle size distribution (PSD), real PSD of Indian railway ballast specification IRS:GE:1:2004, China high-speed rail (HSR) and French rail specifications are generated. PFC3D built-in linear contact model is used to simulate the interaction of ballast particles under various normal stresses, shearing rate and shear box sizes. The results indicate how shear resistance and volumetric changes in ballast assembly are affected by normal stress, shearing rate, PSD and shear box size. In addition to macroscopic behaviour, DEM represents the microscopic behaviour of ballast particles in the form of particle displacement at different stages of the shearing process.

• Proceedings of the KSRS Conference
• /
• 2008.10a
• /
• pp.297-300
• /
• 2008

High-resolution mesoscale meteorological modeling requires more accurate and higher resolution digital elevation model (DEM) data. Shuttle Radar Topographic Mission (SRTM) has created 90 m DEM for entire globe and that is freely available for meteorological modeling and environmental applications. In this research, the effects of the topographic interpolation methods on high-resolution wind field simulation in the coastal regions were quantitatively analyzed using Weather Research and Forecasting (WRF) model with SRTM data. Sensitivity experiments with three different interpolation schemes (four-point bilinear, sixteen-point overlapping parabolic and nearest neighbor interpolation methods) were preformed using SRTM. In WRF modeling with sixteen-point overlapping parabolic interpolation, the coastal line and some small islands show more clearly than other cases. The maximum height of inland is around 140 meters higher, while the minimum of sea height is about 80 meter lower. As it concerns the results of each scheme it seems that the sixteen-point overlapping parabolic scheme indicates the well agreement with observed surface wind data. Consequently, topographic changes due to interpolation methods can lead to the significant influence on mesoscale wind field simulation of WRF modeling.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2005.03a
• /
• pp.751-758
• /
• 2005

In this study, the Discrete Element Code was applied to the analysis of falling rock blocks. The simulation was performed using the PFC2D computer code. Falling rock blocks should be applied as additional force to each others. The force affect the motion of falling rock blocks. This was used to find out the behavior of each blocks. This study revealed that the DEM can successfully capture the behavior of falling rock blocks.

• Proceedings of the Korea Society for Simulation Conference
• /
• 2001.05a
• /
• pp.125-125
• /
• 2001

실 지형정보를 이용한 3D 가상환경은 사용자에게 좀 더 현실에 가까운 교육환경을 제공한다. 3D 가상환경에 사용되는 지형정보는 mesh를 생성하기 위한 고도 data와 mapping을 위한 위성영상이나 항공사진 등이 사용된다. 고도 데이터는 DEM,DTED와 같은 데이터 포맷이 있는데 해상도에 따라 초단위 또는 M 단위로 다양하게 분류되어 있으며 위성영상이나 항공사진도 해상도에 따라 50M∼10Cm 까지 다양하여 사용목적에 맞는 데이터 선택이 필요하다. 고도데이터와 mapping 데이터를 이용하여 기본적인 3D 지형을 생성한 후에 안개나 비, 눈, 빛, 구름과 같은 기상환경을 시뮬레이션하거나 건물이나 이정표, 또는 텍스트 같은 사용자 정보를 Vector overlay 하여 좀 더 현실감 있는 3D 가상환경을 만들 수 있다. 최근에는 인터넷이 일반화되면서 네트웍을 통해 지명데이터를 전송하고 렌더링 하고자 하는 요구가 발생하고 있다. 그러나 3차원 가상환경을 위한 지형 데이터는 2D 데이터에 비해 크기가 크고 고사양의 하드웨어사양을 필요로 하여 네트웍을 통해 전송하고 랜더링 하기에는 여러 가지 제약이 따른다. 이러한 재약을 극복하기 위해 데이터를 한꺼번에 전송하지 않고 점진적으로 전송하고자 하는 연구가 많이 있어 왔으며 점진적 메쉬나 딜로니 규칙에 기반한 TIN 압축 점진적 시각화 기법, DEM 웨이블릿 변환을 적용한 저장, 전송 렌더링 하고자 하는 연구가 시도되어 왔다.

• Proceedings of the Korean Powder Metallurgy Institute Conference
• /
• 2006.09a
• /
• pp.187-188
• /
• 2006

In contrast with the Finite Element Method, the Discrete Element Method (DEM) takes explicitly into account the particulate nature of powders. DEM exhibits some drawbacks and many advantages. Simulations can be computationally expensive and they are only able to represent a volume element. However, these simulations have the great advantage of providing a wealth of information at the microstructural level. Here we demonstrate that the method is well suited for modelling, in coordination with FEM, the compaction of ceramic $UO_2$ particles that have been aggregated. Aggregates of individual ceramic crystallites that are strongly bonded together are represented by porous spheres.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.13 no.1
• /
• pp.136-146
• /
• 2021

This paper aims to numerically estimate the dynamic ice load on a conical structure. The Discrete Element Method (DEM) is employed to model the level ice as the assembly of numerous spherical particles. To mimic the realistic fracture mechanism of ice, the parallel bonding method is introduced. Cases with four different ice drifting velocities are considered in time domain. For validation, the statistics of time-varying ice forces and their frequencies obtained by numerical simulations are extensively compared against the physical model-test results. Ice properties are directly adopted from the targeted experimental test set up. The additional parameters for DEM simulations are systematically determined by a numerical three-point bending test. The findings reveal that the numerical simulation estimates the dynamic ice force in a reasonably acceptable range and its results agree well with experimental data.

• Korean Chemical Engineering Research
• /
• v.55 no.4
• /
• pp.456-466
• /
• 2017

This work investigated the particle morphology change to difference in milling media in a metal based composite fabrication process using a stirred ball mill with ball behavior of DEM simulation. A simulation of the three dimensional motion of grinding media in the stirred ball mill for the research of grinding mechanism to clarify the force, kinetic energy, and medium velocity of grinding media were calculated. In addition, the rotational speed of the stirred ball mill was changed to the experimental conditions for the composite fabrication, and change of the input energy was also calculated while changing the ball material, the flow velocity, and the friction coefficient under the same conditions. As the rotating speed of the stirred ball mill increased, the impact energy between the grinding media to media, media to wall, and media and the stirrer increased quantitatively. Also, we could clearly analyze the change of the particle morphology under the same experimental conditions, and it was found that the ball behavior greatly influences in the particle morphology changes.