• 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.

• Journal of Wetlands Research
• /
• v.19 no.1
• /
• pp.132-143
• /
• 2017

It is essential to obtain periodic sediment discharge data in a river in order to minimize problems that may arise from the erosion, transport, and deposition of sediment. However, it is difficult to estimate sediment discharge by the sediment discharge measurement plan in Korea at present, and empirical fomulas or numerical models are used to replace them. This paper has applied the K-DRUM model, a grid-based rainfall-runoff-sediment model, to estimate sediment discharge and ensure the continuity of the data in the watershed. Discharge and sediment load in 17 watersheds were estimated and the applicability of the model was analyzed through comparisons with measured data. For quantitative evaluation, NSE, PBIAS and RSR items were used, and discharge results reflected the tendency of rainfall and showed high statistical value. In case of sediment discharge, the soil erosion process of the watershed is physically well reflected. When the calibration was performed using the measure data, the applicability seems to be excellent in estimating the continuous sediment discharge data in the real watershed.

• Computers and Concrete
• /
• v.28 no.5
• /
• pp.521-532
• /
• 2021

Present paper focuses on the modeling of size effect on the compressive strength of normal concrete with the application of Discrete Element Method (DEM). Test specimens with different size and shape were cast and uniaxial compressive strength test was performed on each sample. Five different concrete mixes were used, all belonging to a different normal strength concrete class (C20/25, C30/37, C35/45, C45/55, and C50/60). The numerical simulations were carried out by using the PFC 5 software, which applies rigid spheres and contacts between them to model the material. DEM modeling of size effect could be advantageous because the development of micro-cracks in the material can be observed and the failure mode can be visualized. The series of experiments were repeated with the model after calibration. The relationship of the parallel bond strength of the contacts and the laboratory compressive strength test was analyzed by aiming to determine a relation between the compressive strength and the bond strength of different sized models. An equation was derived based on Bazant's size effect law to estimate the parallel bond strength of differently sized specimens. The parameters of the equation were optimized based on measurement data using nonlinear least-squares method with SSE (sum of squared errors) objective function. The laboratory test results showed a good agreement with the literature data (compressive strength is decreasing with the increase of the size of the specimen regardless of the shape). The derived estimation models showed strong correlation with the measurement data. The results indicated that the size effect is stronger on concretes with lower strength class due to the higher level of inhomogeneity of the material. It was observed that size effect is more significant on cube specimens than on cylinder samples, which can be caused by the side ratios of the specimens and the size of the purely compressed zone. A limit value for the minimum size of DE model for cubes and cylinder was determined, above which the size effect on compressive strength can be neglected within the investigated size range. The relationship of model size (particle number) and computational time was analyzed and a method to decrease the computational time (number of iterations) of material genesis is proposed.

• Earthquakes and Structures
• /
• v.25 no.1
• /
• pp.27-41
• /
• 2023

This study focused on nonlinear effective stress site response analysis using two coupled constitutive models, that is, the DM model (Dafalias and Manzari 2004), which incorporated a simple plasticity sand model accounting for fabric change effects, and the PMDY03 model (Khosravifar et al. 2018), that is, a 3D model for earthquake-induced liquefaction triggering and postliquefaction response. A detailed parametric study was conducted to validate the effectiveness of nonlinear site response analysis and porewater pressure (PWP) generation through a true coupled formulation for assessing the initiation of liquefaction at ground level. The coupled models demonstrated accurate prediction of liquefaction triggering, which was in line with established empirical liquefaction triggering relations in published databases. Several limitations were identified in the evaluation of liquefaction using the cyclic stress method, despite its widespread implementation for calculating liquefaction triggering. Variations in shear stiffness, represented by changes in shear wave velocity (Vs1), exerted the most significant influence on site response. The study further indicated that substantial differences in response spectra between nonlinear total stress and nonlinear effective stress analyses primarily occurred when liquefaction was triggered or on the verge of being triggered, as shown by excess PWP ratios approaching unity. These differences diminished when liquefaction occurred towards the later stages of intense shaking. The soil response was predominantly influenced by the higher stiffness values present prior to liquefaction. A key contribution of this study was to validate the criteria used to assess the triggering of level-ground liquefaction using true coupled effective-stress constitutive models, while also confirming the reliability of numerical approximations including the PDMY03 and DM models. These models effectively captured the principal characteristics of liquefaction observed in field tests and laboratory experiments.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.37 no.1
• /
• pp.43-59
• /
• 2017

This study presents the prediction methodology of debris flow occurrence areas using the SINMAP model. Former studies used a single calibration region applying some of the soil test results to predict debris flow occurrence in SINMAP model, which couldn't subdivide the soil properties for the target areas. On the other hands, a multi-calibration region using a detailed soil map and soil strength parameters (c, ${\phi}$) for each soil series to make up for limitation of former studies is proposed. In this process, soils with soil erodibility factor (K) are classified into three types: 1) gravel and gravelly soil. 2) sand and sandy soil, and 3) silt and clay. In addition, T/R estimation method using mean elevation of target area instead of T/R method using actual occurrence time is suggested in this study. The suggested method is applied to Seobyeok-1 ri area, Bonghwa-gun where debris flow occurred. As a result of comparison between two T/R estimation method, both T/R estimations are almost equal. Therefore, the suggested methodologies in this study will contribute to set up the national-wide mitigation plan against debris flow occurrence.

• Journal of Soil and Groundwater Environment
• /
• v.20 no.1
• /
• pp.65-75
• /
• 2015

An artificial recharge test aimed at investigating transport characteristics of the injected water plume in a fractured rock aquifer was conducted. The test used an injection well for injecting tap water whose temperature and electrical conductivity were different from the groundwater. Temporal and depth-wise variation of temperature and electrical conductivity was monitored in both the injection well and a nearby observation well. A highly permeable fracture zone acting as the major pathway of groundwater flow was distinctively revealed in the monitoring data. A finite element subsurface flow and transport simulator (FEFLOW) was used to investigate sensitivity of the transport process to associated aquifer parameters. Simulated results showed that aperture thickness of the fracture and the hydraulic gradient of groundwater highly affected spatio-temporal variation of temperature and electrical conductivity of the injected water plume. The study suggests that artificial recharge of colder water in a fractured rock aquifer could create a thermal plume persistent over a long period of time depending on hydro-thermal properties of the aquifer as well as the amount of injected water.

• Advances in concrete construction
• /
• v.7 no.4
• /
• pp.241-247
• /
• 2019

In this paper the effect of bedding layer on the failure mechanism of rock in direct shear test has been investigated using particle flow code, PFC. For this purpose, firstly calibration of pfc2d was performed using Brazilian tensile strength. Secondly direct shear test consisting bedding layer was simulated numerically. Thickness of layers was 10 mm and rock bridge length was 10 mm, 40 mm and 60 mm. In each rock bridge length, bedding layer angles changes from $0^{\circ}$ to $90^{\circ}$ with increment of $15^{\circ}$. Totally 21 models were simulated and tested. The results show that two types of cracks develop within the model. Shear cracks and tensile cracks. Also failure pattern is affected by bridge length while shear strength is controlled by failure pattern. It's to be noted that bedding layer has not any effect on the failure pattern because the layer interface strength is too high.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.4
• /
• pp.875-886
• /
• 1994

This study is to develop an advanced storm runoff analysis program which takes geomorphological characteristics of watershed into consideration in determining model parameters. Basic concept of storm runoff modelling is based upon the kinematic wave theory. And numerical solution is obtained by the characteristic curve method. The storm runoff analysis program developed by this study is composed of multiple equivalent roughness sub-basins, each of which has two equivalent catchments on both side of a stream. Because it is based upon the stream-order of the Strahler system, the equivalent catchment-stream network reflects the stochastic geomorphological characteristics in the model parameter. Applicability and reliability of the storm runoff analysis program is evidenced by model calibration and verification process utilizing geographical and hydrological data of the Bocheong-river area which is a representative watershed of IHP projects in Korea. This study will hopefully contribute to hydrological calculation essentially required to understand water quality effect caused by regional development.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.46 no.4
• /
• pp.338-344
• /
• 2018

The lunar surface temperature is important as a environmental parameter for the thermal design of the lunar exploration vehicles such as orbital spacecraft, lander, and rovers. In this study, the temperature is numerically predicted through a simplified lumped system model for the energy conservation. The physical values required for the analysis of the energy equation are derived by considering the geometric shape, and the values presented in the previous research results. The areal specific heat, which is the most important thermo-physical property of the lumped system model, was extracted from the temperature measurements by the Diviner loaded on the LRO, and the value was predicted by calibration of the analytical model to the measurements. The predicted temperature distribution obtained through numerical integration has sufficient accuracy to be applied to the thermal design of the lunar exploration vehicles.

• The Journal of the Acoustical Society of Korea
• /
• v.23 no.2
• /
• pp.125-133
• /
• 2004

Most of the way shape estimation method using reference sources assume that the reference sources are in the farfield. That is, the reference sources are assumed to be far from the array. However, in applications of the array with reference sources, the reference sources are not far from the way, so that in practical ocean environments, the conventional method using farfield source model fail to estimate the positions of the hydrophones. In this paper, based on the nearfield source model, a subspace-based array shape estimation method was proposed. In the proposed method, nearfield reference source is modeled using the differential time delay at each hydrophone, and nearfield parameters are derived. Using these parameters, a subspace-based array shape estimation method that generalizes the existing farfield subspace fitting method which can work regardless of the range of the source is proposed. The Cramer-Rao lower bound for the proposed method is investigated. The results of the numerical experiments indicate that the proposed method performs well in estimating the shape of a perturbed way regardless of the ranges of the reference sources.