• Journal of Soil and Groundwater Environment
• /
• v.22 no.6
• /
• pp.66-73
• /
• 2017

Empirical erosion models like Universal Soil Loss Equation (USLE) models have been widely used to make spatially distributed soil erosion vulnerability maps. Even if the models detect vulnerable sites relatively well utilizing big data related to climate, geography, geology, land use, etc within study domains, they do not adequately describe the physical process of soil erosion on the ground surface caused by rainfall or overland flow. In other words, such models are still powerful tools to distinguish the erosion-prone areas at large scale, but physics-based models are necessary to better analyze soil erosion and deposition as well as the eroded particle transport. In this study a physics-based soil erosion modeling system was developed to produce both runoff and sediment yield time series at watershed scale and reflect them in the erosion and deposition maps. The developed modeling system consists of 3 sub-systems: rainfall pre-processor, geography pre-processor, and main modeling processor. For modeling system validation, we applied the system for various erosion cases, in particular, rainfall-runoff-sediment yield simulation and estimation of probable maximum sediment (PMS) correlated with probable maximum rainfall (PMP). The system provided acceptable performances of both applications.

• Ecology and Resilient Infrastructure
• /
• v.7 no.4
• /
• pp.336-344
• /
• 2020

The prediction and evaluation of driftwood accumulation around river-crossing structures are essential because driftwood accumulation increases during flood disasters. In this study, the driftwood accumulation and behavior around bridge piers were evaluated via a numerical model that could be employed to analyze three-dimensional turbulent flow and driftwood motion. The moving particle semi-implicit-based model for driftwood motion was sensitive to the number of spheres. The numerical results showed that the approach velocity and the ratio of driftwood length to pier width were the key factors influencing driftwood accumulation, whereas the driftwood density had only a minor influence. Overall, it is expected that this study will contribute to the development of improved risk evaluation indexes for assessing driftwood accumulation around river-crossing structures.

• Advances in concrete construction
• /
• v.13 no.1
• /
• pp.83-99
• /
• 2022

The two, NBD and SCB tests using gypsum circular discs each containing a single notch have been experimentally accomplished in a rock mechanics laboratory. These specimens have also been numerically modelled by a two-dimensional particle flow which is based on Discrete Element Method (DEM). Each testing specimen had a thickness of 5 cm with 10 cm in diameter. The specimens' lengths varied as 2, 3, and 4 cm; and the specimens' notch angles varied as 0°, 45° and 90°. Similar semi-circular gypsum specimens were also prepared each contained one edge notch with angles 0° or 45°. The uniaxial testing machine was used to perform the experimental tests for both NBD and SCB gypsum specimens. At the same time, the numerical simulation of these tests were performed by PFC2D. The experimental results showed that the failure mechanism of rocks is mainly affected by the orientations of joints with respect to the loading directions. The failure mechanism and fracturing patterns of the gypsum specimens are directly related to the final failure loading. It has been shown that the number of induced tensile cracks showing the specimens' tensile behavior, and increases by decreasing the length and angle of joints. It should be noted that the fracture toughness of rocks' specimens obtained by NBD tests was higher than that of the SCB tests. The fracture toughness of rocks usually increases with the increasing of joints' angles but increasing the joints' lengths do not change the fracture toughness. The numerical solutions and the experimental results for both NDB and SCB tests give nearly similar fracture patterns during the loading process.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2015.05a
• /
• pp.217-217
• /
• 2015

The objective of this study was mainly to evaluate the water resources potential of Lake Tana Basin (LTB) by using Soil and Water Assessment Tool (SWAT). From SWAT simulation of LTB, about 5236 km2 area of LTB is gauged watershed and the remaining 9878 km2 area is ungauged watershed. For calibration of model parameters, four gauged stations were considered namely: Gilgel Abay, Gummera, Rib, and Megech. The SWAT-CUP built-in techniques, particle swarm optimization (PSO) and generalized likelihood uncertainty estimation (GLUE) method was used for calibration of model parameters and PSO method were selected for the study based on its performance results in four gauging stations. However the level of sensitivity of flow parameters differ from catchment to catchment, the curve number (CN2) has been found the most sensitive parameters in all gauged catchments. To facilitate the transfer of data from gauged catchments to ungauged catchments, clustering of hydrologic response units (HRUs) were done based on physical similarity measured between gauged and ungauged catchment attributes. From SWAT land use/ soil use/slope reclassification of LTB, a total of 142 HRUs were identified and these HRUs are clustered in to 39 similar hydrologic groups. In order to transfer the optimized model parameters from gauged to ungauged catchments based on these clustered hydrologic groups, this study evaluates three parameter transfer schemes: parameters transfer based on homogeneous regions (PT-I), parameter transfer based on global averaging (PT-II), and parameter transfer by considering Gilgel Abay catchment as a representative catchment (PT-III) since its model performance values are better than the other three gauged catchments. The performance of these parameter transfer approach was evaluated based on values of Nash-Sutcliffe efficiency (NSE) and coefficient of determination (R2). The computed NSE values was found to be 0.71, 0.58, and 0.31 for PT-I, PT-II and PT-III respectively and the computed R2 values was found to be 0.93, 0.82, and 0.95 for PT-I, PT-II, and PT-III respectively. Based on the performance evaluation criteria, PT-I were selected for modelling ungauged catchments by transferring optimized model parameters from gauged catchment. From the model result, yearly average stream flow for all homogeneous regions was found 29.54 m3/s, 112.92 m3/s, and 130.10 m3/s for time period (1989 - 2005) for region-I, region-II, and region-III respectively.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.10 no.4
• /
• pp.295-303
• /
• 2012

Based on the data observed and analyzed on a groundwater flow system in the KURT (KAERI Underground Research Tunnel) site, the transport of radionuclides, which were assumed to be released at the supposed position, was calculated on the time-domain. A groundwater pathway from the release position to the surface was identified by simulating the groundwater flow model with the hydrogeological characteristics measured from the field tests in the KURT site. The elapsed time when the radionuclides moved through the pathway is evaluated using TDRW (Time Domain Random Walk) method for simulating the transport on the time-domain. Some retention mechanisms, such as radioactive decay, equilibrium sorption, and matrix diffusion, as well as the advection-dispersion were selected as the factors to influence on the elapsed time. From the simulation results, the effects of the sorption and matrix diffusion, determined by the properties of the radionuclides and underground media, on the transport of the radionuclides were analyzed and a decay chain of the radionuclides was also examined. The radionuclide ratio of the mass discharge into the surface environment to the mass released from the supposed repository did not exceed $10^{-3}$, and it decreased when the matrix diffusion were considered. The method used in this study could be used in preparing the data on radionuclide transport for a safety assessment of a geological disposal facility because the method could evaluate the travel time of the radionuclides considering the transport retention mechanism.

• Journal of the Korean Institute of Intelligent Systems
• /
• v.22 no.6
• /
• pp.722-728
• /
• 2012

Urban sewer systems has a limitation of capacity of rainwater storage and problem of occurrence of untreated sewage, so adopting a storage facility for sewer flooding prevention and urban non-point pollution reduction has a big attention. The Korea Ministry of Environment has recently introduced a new concept of "multi-functional storage facility", which is crucial not only in preventive stormwater management but also in dealing with combined sewer overflow and sanitary sewer discharge, and also has been promoting its adoption. However, reserving a space for a single large-scale storage facility might be difficult especially in urban areas. Thus, decentralized construction of small- and midium-sized storage facilities and its operation have been introduced as an alternative way. In this paper, we propose a model predictive control scheme for an optimized operation of distributed storage facilities and sewer networks. To this aim, we first describe the mathematical model of each component of networks system which enables us to analyze its detailed dynamic behavior. Second, overflow locations and volumes will be predicted based on the developed network model with data on the external inflow occurred at specific locations of the network. MPC scheme based on the introduced particle swarm optimization technique then produces the optimized the gate setting for sewer network flow control, which minimizes sewer flooding and maximizes the potential storage capacity. Finally, the operational efficacy of the proposed control scheme is demonstrated by simulation study with virtual rainstorm event.

• Transactions of the Korean hydrogen and new energy society
• /
• v.22 no.3
• /
• pp.386-401
• /
• 2011

Mathematical models for various steps in coal gasification reactions were developed and applied to investigate the effects of operation parameters on dynamic behavior of gasification process. Chemical reactions considered in these models were pyrolysis, volatile combustion, water shift reaction, steam-methane reformation, and char gasification. Kinetics of heterogeneous reactions between char and gaseous agents was based on Random pore model. Momentum balance and Stokes' law were used to estimate the residence time of solid particles (char) in an up-flow reactor. The effects of operation parameters on syngas composition, reaction temperature, carbon conversion were verified. Parameters considered here for this purpose were $O_2$-to-coal mass ratio, pressure of reactor, composition of coal, diameter of char particle. On the basis of these parametric studies some quantitative parameter-response relationships were established from both dynamic and steady-state point of view. Without depending on steady state approximation, the present model can describe both transient and long-time limit behavior of the gasification system and accordingly serve as a proto-type dynamic simulator of coal gasification process. Incorporation of heat transfer through heterogenous boundaries, slag formation and steam generation is under progress and additional refinement of mathematical models to reflect the actual design of commercial gasifiers will be made in the near futureK.

• International Journal of Highway Engineering
• /
• v.8 no.1 s.27
• /
• pp.139-152
• /
• 2006

Many experimental and numerical approaches have been developed to evaluate paving materials and to predict pavement response and distress. Micromechanical simulation modeling is a technology that can reduce the number of physical tests required in material formulation and design and that can provide more details, e.g., the internal stress and strain state, and energy evolution and dissipation in simulated specimens with realistic microstructural features. A clustered distinct element modeling (DEM) approach was implemented In the two-dimensional particle flow software package (PFC-2D) to study the complex behavior observed in asphalt mixture fracturing. The relationship between continuous and discontinuous material properties was defined based on the potential energy approach. The theoretical relationship was validated with the uniform axial compression and cantilever beam model using two-dimensional plane strain and plane stress models. A bilinear cohesive displacement-softening model was implemented as an intrinsic interface and applied for both homogeneous and heterogeneous fracture modeling in order to simulate behavior in the fracture process zone and to simulate crack propagation. A disk-shaped compact tension test (DC(T)) with heterogeneous microstructure was simulated and compared with the experimental fracture test results to study Mode I fracture. The realistic arbitrary crack propagation including crack deflection, microcracking, crack face sliding, crack branching, and crack tip blunting could be represented in the fracture models. This micromechanical modeling approach represents the early developmental stages towards a 'virtual asphalt laboratory,' where simulations of laboratory tests and eventually field response and distress predictions can be made to enhance our understanding of pavement distress mechanisms, such its thermal fracture, reflective cracking, and fatigue crack growth.

• Geomechanics and Engineering
• /
• v.11 no.3
• /
• pp.345-359
• /
• 2016

In order to investigate the influence of the interfacial angel on failure characteristics and mechanism of combined coal-rock mass, 35 uniaxial/biaxial compressive simulation tests with 5 different interfacial angels of combined coal-rock samples were conducted by PFC2D software. The following conclusions are drawn: (1) The compressive strength and cohesion decrease with the increase of interfacial angle, which is defined as the angle between structure plane and the exterior normal of maximum principal plane, while the changes of elastic modulus and internal friction angle are not obvious; (2) The impact energy index $K_E$ decreases with the increase of interfacial angle, and the slip failure of the interface can be predicted based on whether the number of acoustic emission (AE) hits has multiple peaks or not; (3) There are four typical failure patterns for combined coal-rock samples including I (V-shaped shear failure of coal), II (single-fracture shear failure of coal), III (shear failure of rock and coal), and IV (slip rupture of interface); and (4) A positive correlation between interfacial angle and interface effect is shown obviously, and the interfacial angle can be divided into weak-influencing scope ($0-15^{\circ}$), moderate-influencing scope ($15-45^{\circ}$), and strong-influencing scope (> $45^{\circ}$), respectively. However, the confining pressure has a certain constraint effect on the interface effect.

• Wind and Structures
• /
• v.23 no.2
• /
• pp.127-142
• /
• 2016

Aerodynamic effects, such as drag force and flow-induced vibration (FIV), on civil engineering structures can be minimized by optimally modifying the structure shape. This work investigates the turbulent wake of a square prism with its faces modified into a sinusoidal wave along the spanwise direction using three-dimensional large eddy simulation (LES) and particle image velocimetry (PIV) techniques at Reynolds number $Re_{Dm}$ = 16,500-22,000, based on the nominal width ($D_m$) of the prism and free-stream velocity ($U_{\infty}$). Two arrangements are considered: (i) the top and bottom faces of the prism are shaped into the sinusoidal waves (termed as WSP-A), and (ii) the front and rear faces are modified into the sinusoidal waves (WSP-B). The sinusoidal waves have a wavelength of $6D_m$ and an amplitude of $0.15D_m$. It has been found that the wavy faces lead to more three-dimensional free shear layers in the near wake than the flat faces (smooth square prism). As a result, the roll-up of shear layers is postponed. Furthermore, the near-wake vortical structures exhibit dominant periodic variations along the spanwise direction; the minimum (i.e., saddle) and maximum (i.e., node) cross-sections of the modified prisms have narrow and wide wakes, respectively. The wake recirculation bubble of the modified prism is wider and longer, compared with its smooth counterpart, thus resulting in a significant drag reduction and fluctuating lift suppression (up to 8.7% and 78.2%, respectively, for the case of WSP-A). Multiple dominant frequencies of vortex shedding, which are distinct from that of the smooth prism, are detected in the near wake of the wavy prisms. The present study may shed light on the understanding of the underlying physical mechanisms of FIV control, in terms of passive modification of the bluff-body shape.