• Journal of Korean Tunnelling and Underground Space Association
• /
• v.25 no.6
• /
• pp.403-421
• /
• 2023

As the number of underground structures has increased in recent decades, it has become crucial to predict geological hazards ahead of a tunnel face during tunnel construction. Consequently, this study developed a finite element (FE) numerical model to simulate electrical resistivity surveys in tunnel boring machine (TBM) operations for predicting mixed ground conditions in front of tunnel faces. The accuracy of the developed model was verified by comparing the numerical results not only with an analytical solution but also with experimental results. Using the developed model, a series of parametric studies were carried out to estimate the effect of geological conditions and sensor geometric configurations on electrical resistivity measurements. The results of these studies showed that both the interface slope and the difference in electrical resistivity between two different ground formations affect the patterns and variations in electrical resistivity observed during TBM excavation. Furthermore, it was revealed that selecting appropriate sensor spacing and optimizing the location of the electrode array were essential for enhancing the efficiency and accuracy of predictions related to mixed ground conditions. In conclusion, the developed model can serve as a powerful and reliable tool for predicting mixed ground conditions during TBM tunneling.

• Journal of the Korean Association of Geographic Information Studies
• /
• v.22 no.4
• /
• pp.39-58
• /
• 2019

This study performed the dam watershed storm runoff modeling using GPM(Global Precipitation Measurement) satellite rain and KIMSTORM2(KIneMatic wave STOrm Runoff Model 2) distributed model. For YongdamDam watershed(930㎢), three heavy rain events of 25th August 2014, 11th September 2017, and 26th June 2018 were selected and tested for 4 cases of spatial rainfalls such as (a) Kriging interpolated data using ground observed data at 7 stations, (b) original GPM data, (c) GPM corrected by CM(Conditional Merging), and GPM corrected by GDA(Geographical Differential Analysis). For the 4 kinds of data(Kriging, GPM, CM-GPM, and GDA-GPM), the KIMSTORM2 was calibrated respectively using the observed flood discharges at 3 water level gauge stations(Cheoncheon, Donghyang, and Yongdam) with parameters of initial soil moisture contents, stream Manning's roughness coefficient, and effective hydraulic conductivity. The total average Nash-Sutcliffe efficiency(NSE) for the 3 events and 3 stations was 0.94, 0.90, 0.94, and 0.94, determination coefficient(R²) was 0.96, 0.92, 0.97 and 0.96, the volume conservation index(VCI) was 1.03, 1.01, 1.03 and 1.02 for Kriging, GPM, CM-GPM, and GDA-GPM applications respectively. The CM-GPM and GDA-GPM showed better results than the original GPM application for peak runoff and runoff volume simulations, and they improved NSE, R², and VCI results.

• Tunnel and Underground Space
• /
• v.30 no.4
• /
• pp.359-381
• /
• 2020

Within the framework of DECOVALEX-2019 Task D, full-scale engineered barriers experiment (FEBEX) at Grimsel Test Site was numerically simulated to investigate an applicability of implemented Barcelona basic model (BBM) into TOUGH2-MP/FLAC3D simulator, which was developed for the prediction of the coupled thermo-hydro-mechanical behavior of bentonite buffer. And the calculated heater power, temperature, relative humidity, total stress, saturation, water content and dry density were compared with in situ data monitored in the various sections. In general, the calculated heater power and temperature provided a fairly good agreement with experimental observations, however, the difference between power of heater #1 and that of heater #2 could not captured in the numerical analysis. It is necessary to consider lamprophyre with low thermal conductivity around heater #1 and non-simplified installation progresses of bentonite blocks in the tunnel for better modeling results. The evolutions and distributions of relative humidity were well reproduced, but hydraulic model needs to be modified because the re-saturation process was relatively fast near the heaters. In case of stress evolutions due to the thermal and hydraulic expansions, the computed stress was in good agreement with the data. But, the stress is slightly higher than the measured in situ data at the early stage of the operation, because gap between rock mass and bentonite blocks have not been considered in the numerical simulations. The calculated distribution of saturation, water content, and dry density along the radial distance showed good agreement with the observations after the first and final dismantling. The calculated dry density near the center of the FEBEX tunnel and heaters were overestimated compared with the observations. As a result, the saturation and water content were underestimated with the measurements. Therefore, numerical model of permeability is needed to modify for the production of better numerical results. It will be possible to produce the better analysis results and more realistically predict the coupled THM behavior in the bentonite blocks by performing the additional studies and modifying the numerical model based on the results of this study.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.24 no.5
• /
• pp.628-636
• /
• 2018

Because large buoys are mainly made of steel, they are heavy and vulnerable to corrosion by sea water. This makes buoy installation and maintenance difficult. Moreover, vessel collision accidents with buoys and damage to vessels due to the material of buoys (e.g., steel) are reported every year. Recently, light buoys adopting eco-friendly and lightweight materials have come into the spotlight in order to solve the previously-mentioned problems. In Korea, a new lightweight buoy with a 7-Nautical Mile lantern adopting expanded polypropylene (EPP) and aluminum to create a buoyant body and tower structure, respectively, was developed in 2017. When these light buoys are operated in the ocean, the visibility and angle of light from the lantern installed on the light buoys changes, which may cause them to function improperly. Therefore, research on the performance of light buoys is needed since the weight distribution and motion characteristics of these new buoys differ from conventional models. In this study, stability estimation and motion analyses for newly-developed buoys under various environmental conditions considering a mooring line were carried out using ANSYS AQWA. Numerical simulations for the estimation of wind and current loads were performed using commercial CFD software, Siemens STAR-CCM+, to increase the accuracy of motion analysis. By comparing the estimated maximum significant motions of the light buoys, it was found that waves and currents were more influential in the motion of the buoys. And, the estimated motions of the buoys became larger as the sea state became worser, which might be the reason that the peak frequencies of the wave spectra got closer to those of the buoys.

• Journal of Korean Society of Forest Science
• /
• v.106 no.1
• /
• pp.77-86
• /
• 2017

This study was conducted to analyze the relationship among climatic factors and radial growth of Pinus densiflora and Larix kaempferi in South Korea. To determine the climate-growth relationship, cluster analysis was applied to group surveyed regions by the climatical similarity, and a dendroclimatological model was developed to predict radial growth for each climate group under the RCP 4.5 and RCP 8.5 scenarios for greenhouse gases. The cluster analysis showed four climatic clusters (Cluster 1~4) from 10 regions for P. densiflora and L. kaempferi. The dendroclimatological model was developed through climatic variables and standardized residual chronology for each climatic cluster of P. densiflora and L. kaempferi. Four climatic variables were used in the models for P. densiflora ($R^2$ values between 0.38 to 0.58). Two to five climatic variables were used in the models for L. kaempferi ($R^2$ values between 0.31 to 0.43). The growth simulations with two RCP climate-change scenarios(RCP 4.5 and RCP 8.5) were used for growth prediction. The radial growth of the Cluster 4 of P. densiflora, the mountainous region at high elevation, tend to increase, while those of cluster 2 and 3 of P. densiflora, the region of the hightest average temperature, tends to decrease. The radial growth of the Cluster 1 of L. kaempferi the region of the lowest minimum temperature, while that of Cluster 2, the region of the highest average temperature, tends to decrease. The radial growth of Cluster 3 of L. kaempferi, the region in the east coast and Cluster 4, the region at high elevation, tends to hold steady. The results of this study are expected to provide valuable information necessary for predicting changes in radial growth of Pinus densiflora and Larix kaempferi caused by climate change.

• Journal of Navigation and Port Research
• /
• v.42 no.6
• /
• pp.427-436
• /
• 2018

In this study, experimental and numerical analysis were performed on the survivability of a long pipe-type buoy structure in waves. The buoy structure is an articulated tower consisting of an upper structure, buoyancy module, and gravity anchor with long pipes forming the base frame. A series of experiment were performed in the ocean engineering basin of KRISO with the scaled model of 1/ 22 to evaluate the survivability of the buoy structure at West Sea in South Korea. Survival condition was considered as the wave of 50 year return period. Additional experiments were performed to investigate the effects of current and wave period. The factors considered for the evaluation of the buoy's survival were the pitch angle of the structure, anchor reaction force, and the number of submergence of the upper structure. Numerical simulations were carried out with the OrcaFlex, the commercial program for the mooring analysis, with the aim of performing mutual validation with the experimental results. Based on the evaluation, the behavior characteristics of the buoy structure were first examined according to the tidal conditions. The changes were investigated for the pitch angle and anchor reaction force at HAT and LAT conditions, and the results directly compared with those obtained from numerical simulation. Secondly, the response characteristics of the buoy structure were studied depending on the wave period and the presence of current velocity. Third, the number of submergence through video analysis was compared with the simulation results in relation to the submergence of the upper structure. Finally, the simulation results for structural responses which were not directly measured in the experiment were presented, and the structural safety discussed in the survival waves. Through a series of survivability evaluation studies, the behavior characteristics of the buoy structure were examined in survival waves. The vulnerability and utility of the buoy structure were investigated through the sensitivity studies of waves, current, and tides.

• Journal of Korea Water Resources Association
• /
• v.56 no.2
• /
• pp.75-89
• /
• 2023

As the Mekong River basin is a nationally shared river, it is difficult to collect precipitation data, and the quantitative and qualitative quality of the data sets differs from country to country, which may increase the uncertainty of hydrological analysis results. Recently, with the development of remote sensing technology, it has become easier to obtain grid-based precipitation products(GPPs), and various hydrological analysis studies have been conducted in unmeasured or large watersheds using GPPs. In this study, rainfall-runoff simulation in the Mekong River basin was conducted using the SWAT model, which is a quasi-distribution model with three satellite GPPs (TRMM, GSMaP, PERSIANN-CDR) and two GPPs (APHRODITE, GPCC). Four water level stations, Luang Prabang, Pakse, Stung Treng, and Kratie, which are major outlets of the main Mekong River, were selected, and the parameters of the SWAT model were calibrated using APHRODITE as an observation value for the period from 2001 to 2011 and runoff simulations were verified for the period form 2012 to 2013. In addition, using the ConvAE, a convolutional neural network model, spatio-temporal correction of original satellite precipitation products was performed, and rainfall-runoff performances were compared before and after correction of satellite precipitation products. The original satellite precipitation products and GPCC showed a quantitatively under- or over-estimated or spatially very different pattern compared to APHPRODITE, whereas, in the case of satellite precipitation prodcuts corrected using ConvAE, spatial correlation was dramatically improved. In the case of runoff simulation, the runoff simulation results using the satellite precipitation products corrected by ConvAE for all the outlets have significantly improved accuracy than the runoff results using original satellite precipitation products. Therefore, the bias correction technique using the ConvAE technique presented in this study can be applied in various hydrological analysis for large watersheds where rain guage network is not dense.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.37 no.1
• /
• pp.67-76
• /
• 2024

In this paper, we present a DIP-MLS testing method that combines digital image processing with a rigid body-based MLS differencing approach to measure mechanical variables and analyze the impact of target location and image resolution. This method assesses the displacement of the target attached to the sample through digital image processing and allocates this displacement to the node displacement of the MLS differencing method, which solely employs nodes to calculate mechanical variables such as stress and strain of the studied object. We propose an effective method to measure the displacement of the target's center of gravity using digital image processing. The calculation of mechanical variables through the MLS differencing method, incorporating image-based target displacement, facilitates easy computation of mechanical variables at arbitrary positions without constraints from meshes or grids. This is achieved by acquiring the accurate displacement history of the test specimen and utilizing the displacement of tracking points with low rigidity. The developed testing method was validated by comparing the measurement results of the sensor with those of the DIP-MLS testing method in a three-point bending test of a rubber beam. Additionally, numerical analysis results simulated only by the MLS differencing method were compared, confirming that the developed method accurately reproduces the actual test and shows good agreement with numerical analysis results before significant deformation. Furthermore, we analyzed the effects of boundary points by applying 46 tracking points, including corner points, to the DIP-MLS testing method. This was compared with using only the internal points of the target, determining the optimal image resolution for this testing method. Through this, we demonstrated that the developed method efficiently addresses the limitations of direct experiments or existing mesh-based simulations. It also suggests that digitalization of the experimental-simulation process is achievable to a considerable extent.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.26 no.3
• /
• pp.174-183
• /
• 2014

Seabed beneath and near coastal structures may undergo large excess pore water pressure composed of oscillatory and residual components in the case of long durations of high wave loading. This excess pore water pressure may reduce effective stress and, consequently, the seabed may liquefy. If liquefaction occurs in the seabed, the structure may sink, overturn, and eventually increase the failure potential. In this study, to evaluate the liquefaction potential on the seabed, numerical analysis was conducted using the expanded 2-dimensional numerical wave tank to account for an irregular wave field. In the condition of an irregular wave field, the dynamic wave pressure and water flow velocity acting on the seabed and the surface boundary of the composite breakwater structure were estimated. Simulation results were used as input data in a finite element computer program for elastoplastic seabed response. Simulations evaluated the time and spatial variations in excess pore water pressure, effective stress, and liquefaction potential in the seabed. Additionally, the deformation of the seabed and the displacement of the structure as a function of time were quantitatively evaluated. From the results of the analysis, the liquefaction potential at the seabed in front and rear of the composite breakwater was identified. Since the liquefied seabed particles have no resistance to force, scour potential could increase on the seabed. In addition, the strength decrease of the seabed due to the liquefaction can increase the structural motion and significantly influence the stability of the composite breakwater. Due to limitations of allowable paper length, the studied results were divided into two portions; (I) focusing on the dynamic response of structure, acceleration, deformation of seabed, and (II) focusing on the time variation in excess pore water pressure, liquefaction, effective stress path in the seabed. This paper corresponds to (II).

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.26 no.3
• /
• pp.160-173
• /
• 2014

Seabed beneath and near coastal structures may undergo large excess pore water pressure composed of oscillatory and residual components in the case of long durations of high wave loading. This excess pore water pressure may reduce effective stress and, consequently, the seabed may liquefy. If liquefaction occurs in the seabed, the structure may sink, overturn, and eventually increase the failure potential. In this study, to evaluate the liquefaction potential on the seabed, numerical analysis was conducted using the expanded 2-dimensional numerical wave tank to account for an irregular wave field. In the condition of an irregular wave field, the dynamic wave pressure and water flow velocity acting on the seabed and the surface boundary of the composite breakwater structure were estimated. Simulation results were used as input data in a finite element computer program for elastoplastic seabed response. Simulations evaluated the time and spatial variations in excess pore water pressure, effective stress, and liquefaction potential in the seabed. Additionally, the deformation of the seabed and the displacement of the structure as a function of time were quantitatively evaluated. From the results of the analysis, the liquefaction potential at the seabed in front and rear of the composite breakwater was identified. Since the liquefied seabed particles have no resistance to force, scour potential could increase on the seabed. In addition, the strength decrease of the seabed due to the liquefaction can increase the structural motion and significantly influence the stability of the composite breakwater. Due to limitations of allowable paper length, the studied results were divided into two portions; (I) focusing on the dynamic response of structure, acceleration, deformation of seabed, and (II) focusing on the time variation in excess pore water pressure, liquefaction, effective stress path in the seabed. This paper corresponds to (I).