• Journal of the Korean Geotechnical Society
• /
• v.30 no.11
• /
• pp.25-37
• /
• 2014

Clay or sand does not exist alone but various sizes of soil are mixed in the field. In this study, the effect of water content on large deformation of such mixed soils is studied by using soil column tests and a particle method. A soil column with 7 cm in diameter and 13 cm in height, which was made out of kaolinite with sand content of 0, 10, 25, or 50%, was tested for large deformation. Its deformation was monitored with time. While increasing its water content from 40, 60, to 80%, a total of 12 types of soil column tests were carried out. The particle method simulated their deformation with time. A maximum plastic shear modulus, which was a function of undrained shear strength and plasticity index for soils with different water contents, was associated with soil viscosity to simulate large deformation of soil column. When a sand content of soil column was constant, the deformation of soil column increased with increasing water content. When a water content of soil column was constant, large deformation occurred with increasing the sand content. The maximum deformation, which was 22 cm in diameter, was observed in the case of water content of 80% and sand content of 50%. The particle method was able to relatively well simulate such large deformation and stress change of soils.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.48-54
• /
• 2004

The effects of injector spacing s and injector diameter d on mixing are numerically investigated in supersonic combustor with perpendicular injection behind a backward-facing step. Simulations are reported for airstream Mach number of 2.4. Parameters are changed on following 4 cases to investigate the effects of injector configuration on mixing efficiency $\eta_m$. In the case of varying d or s, dynamic pressure ratio $Rq(=(pu^2)_j/(pu^2)_a)$ is also varied to keep bulk equivalence ratio $\Phi({\oe})Rq.d^2/s)$ constant. (l) Injector spacing s is varied at constant $\Phi$=0.5, 1, 2 for injector diameter d=6mm. In the case of $\Phi$=1, $\eta_m$ has its maximum value at s=24mm. The reason is that increase of $\eta_m$. , by widening spacing at Rq=constant competes with decrease of $\eta_m$ by increasing Rq at s=constant. When spacing is narrow, the flow field of vicinity of injector becomes two-dimensional because adjacent jets interferes each other. By widening spacing, air is easily entrained by three-dimensional effect. This mechanism also appears in the case of $\Phi$=0.5, 2 for d=6mm, and $\eta_m$. reaches its maximum value at s=24mm for $\Phi$=0.5 and at s=42mm for $\Phi$=2. (2) In the case of injector diameter d varied at $\Phi$=1 for s=30mm, $\eta_m$. has its maximum value at d=3mm. The reason is that decrease of $\eta_m$ by increasing injector diameter competes with increase of $\eta_m$ by decreasing Rq at d=constant.(3) In the case of s varied at $\Phi$=0.5, 1,2 for d=3mm, the injector spacing at which mixing efficiency has its maximum value is s= 18mm for $\Phi$=0.5, s=24mm for $\Phi$=1, s=24mm for $\Phi$=2. Therefore it is found that d=3mm and s=24mm can be optimum configuration over a range of $\Phi$=0.5~2.(4) The effect of h on the optimum spacing is investigated. s is varied for d=6mm at step height h=4, 6, 8mm. The simulation results do not show significant change on the step height.

• Journal of Soil and Groundwater Environment
• /
• v.11 no.5
• /
• pp.1-8
• /
• 2006

The objective of this study is to assess impacts of sorption heterogeneity on the transport of leachate leaked from unlined landfill sites and is accomplished by examining the results from a series of Monte-Carlo simulations. For random distribution coefficient ($K_{d}$) fields with four different levels of heterogeneity ranging from homogeneous to highly heterogeneous, the transport of leachate was investigated by linking a saturated flow model with a contaminant transport model. Impacts of a chemical heterogeneity were evaluated using point statistics values such as mean, standard deviation, and coefficient of variation of the concentration obtained at monitoring wells from 100 Monte-Carlo trials. Inspection of point statistics shows that the distribution of distribution coefficient in the landfill site proves to be an important parameter in controlling leachate concentrations. In comparison to homogeneous sorption, heterogeneous $K_{d^-}$ fields produce the variability in the leachate concentration for different realizations. The variability increases significantly as the variance in the $K_{d^-}$ field and the travel time between source and monitoring well increase. These outcomes indicate that use of a constant homogeneous $K_{d}$ value for predicting the transport of leachate can result in significant error, especially when variability in $K_{d}$ is high.

• Journal of the Korean Institute of Electrical and Electronic Material Engineers
• /
• v.36 no.4
• /
• pp.413-417
• /
• 2023

Gallium Oxide (Ga₂O₃) is preferred as a material for next generation power semiconductors. The Ga₂O₃ should solve the disadvantages of low thermal resistance characteristics and difficulty in forming an inversion layer through p-type ion implantation. However, Ga₂O₃ is difficult to inject p-type ions, so it is being studied in a heterojunction structure using p-type oxides, such as NiO, SnO, and Cu₂O. Research the lateral-type FET structure of NiO/Ga₂O₃ heterojunction under the Gate contact using the Sentaurus TCAD simulation. At this time, the VG-ID and VD-ID curves were identified by the thickness of the Epi-region (channel) and the doping concentration of NiO of 1×10¹⁷ to 1×10¹⁹ cm-3. The increase in Epi region thickness has a lower threshold voltage from -4.4 V to -9.3 V at I_D = 1×10^-8 mA/mm, as current does not flow only when the depletion of the PN junction extends to the Epi/Sub interface. As an increase of NiO doping concentration, increases the depletion area in Ga₂O₃ region and a high electric field distribution on PN junction, and thus the breakdown voltage increases from 512 V to 636 V at I_D =1×10^-3 A/mm.

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

The utilization of underground spaces in relation to tunnels and energy/waste storage is on the rise. To ensure the stability of underground spaces, it is crucial to reinforce rock fractures and discontinuities. Discontinuities, such as joints, can weaken the strength of the rock and lead to groundwater inflow into underground spaces. In order to enhance the strength and stability of the area around these discontinuities, rock grouting techniques are employed. However, during rock grouting, it is impossible to visually confirm whether the grouting material is being smoothly injected as intended. Without proper injection, the expected increases in strength, durability, and degree of consolidation may not be achieved. Therefore, it is necessary to predict in advance whether the grouting material is being injected as designed. In this study, we aimed to assess the injection performance based on injection variables such as the water/cement mixture ratio, injection pressure, and injection flow using UDEC (Universal Distinct Element Code) numerical program. Additionally, numerical results were validated by the lab experiment. The results of this study are expected to help optimize variables such as injection material properties, injection time, and pump pressure in the grouting design in the field.

• Korean Journal of Agricultural and Forest Meteorology
• /
• v.11 no.4
• /
• pp.192-205
• /
• 2009

The unexpected wind over the Mt. Hwawang on 9 February 2009 was deadly when many spectators were watching a traditional event to burn dried grasses and the fire went out of control due to the wind. We analyzed the fatal wind based on wind flow simulations over a digitized complex terrain of the mountain with a localized heating area using a three dimensional computational fluid dynamics model, CFD_NIMR_SNU (Computational Fluid Dynamics_National Institute of Meteorological Research_Seoul National University). Three levels of fire intensity were simulated: no fire, $300^{\circ}C$ and $600^{\circ}C$ of surface temperature at the site on fire. The surface heat accelerated vertical wind speed by as much as $0.7\;m\;s^{-1}$ (for $300^{\circ}C$) and $1.1\;m\;s^{-1}$ (for $600^{\circ}C$) at the center of the fire. Turbulent kinetic energy was increased by the heat itself and by the increased mechanical force, which in turn was generated by the thermal convection. The heating together with the complex terrain and strong boundary wind induced the unexpected high wind conditions with turbulence at the mountain. The CFD_NIMR_SNU model provided valuable analysis data to understand the consequences of the fatal mountain fire. It is suggested that the place of fire was calm at the time of the fire setting due to the elevated terrain of the windward side. The suppression of wind was easily reversed when there was fire, which caused updraft of hot air by the fire and the strong boundary wind. The strong boundary wind in conjunction with the fire event caused the strong turbulence, resulting in many fire casualties. The model can be utilized in turbulence forecasting over a small area due to surface fire in conjunction with a mesoscale weather model to help fire prevention at the field.

• Journal of Korean Society on Water Environment
• /
• v.23 no.4
• /
• pp.467-473
• /
• 2007

Although the dominant land use at the Imha-dam watershed is forest areas, soil erosion has been increasing because of intensive agricultural activities performed at the fields located along the stream for easy-access to water supply and relatively favorable topography. In addition, steep topography at the Imha-dam watershed is also contributing increased soil erosion and sediment loads. At the Imha-dam watershed, outflow has increased sharply by the typhoons Rusa and Maemi in 2002, 2003 respectively. In this study, the Soil and Water Assessment Tool (SWAT) model was evaluated for simulation of flow and sediment behaviors with long-term temporal and spatial conditions. The precipitation data from eight precipitation observatories, located at Ilwol, Subi and etc., were used. There was no significant difference in monthly rainfall for 8 locations. However, there was slight differences in rainfall amounts and patterns in 2003 and 2004. The topographical map at 1:5000 scale from the National Geographic Information Institute was used to define watershed boundaries, the detailed soil map at 1:25,000 scale from the National Institute of Highland Agriculture and the land cover data from the Korea Institute of Water and Environment were used to simulate the hydrologic response and soil erosion and sediment behaviors. To evaluate hydrologic component of the SWAT model, calibration was performed for the period from Jan. 2002 to Dec. 2003, and validation for Jan. 2004 to Apr. 2005. The $R^2$ value and El value were 0.93 and 0.90 respectively for calibration period, and the $R^2$ value and El value for validation were 0.73 and 0.68 respectively. The $R^2$ value and El value of sediment yield data with the calibrated parameters was 0.89 and 0.84 respectively. The comparisons with the measured data showed that the SWAT model is applicable to simulate hydrology and sediment behaviors at Imha dam watershed. With proper representation of the Best Management Practices (BM Ps) in the SWAT model, the SWAT can be used for pre-evaluation of the cost-effective and sustainable soil erosion BMPs to solve sediment issues at the Imha-dam watershed. In Korea, the Universal Soil Loss Equation (USLE) has been used to estimate the soil loss for over 30 years. However, there are limitations in the field scale mdel, USLE when applied for watershed. Also, the soil loss changes temporarily and spatially, for example, the Imha-dam watershed. Thus, the SW AT model, capable of simulating hydrologic and soil erosion/sediment behaviors temporarily and spatially at watershed scale, should be used to solve the muddy water issues at the Imha-dam watershed to establish more effective muddy water reduction countermeasure.

• Journal of Environmental Impact Assessment
• /
• v.30 no.1
• /
• pp.35-48
• /
• 2021

In this study, the assessment of field applicability of in-situ remediation of nitrate-contaminated groundwater located in Yesan-gun was performed. Zero-valent iron/bio composite media injected PRB (Permeable Reactive Barrier) and monitoring well were installed in the contaminated groundwater site and monitored main remediation indicators during the PRB operation. Nitrate, nitrite, ammonia, Fe ion, TOC, and turbidity were analyzed and the diversity and population of microorganism in the PRB installed site were investigated for the verification of effect of injected PRB. In the study site where is an agricultural area, a river flows from west to east that forms a river boundary and the southern area has an impermeable sector. It was found that nitrate flows into the river, which is similar as groundwater flow. Simulation result for the fate of nitrate in groundwater showed steady state of nitrate arrived after 3~5 years passed. However, it is just to consider current conditions with no additional input of contaminant source, if additional input of contaminant source occurs contamination dispersion and time for steady state are expected to be increased. The monitoring results showed that Fe ion, TOC and turbidity in groundwater were not clearly changed in concentration after PRB installation, which indicates adaptability of the injected PRB for remediation of groundwater with no additional harmful effect to water quality. The concentration of nitrate maintained less than 5mg/L until 42 days after PRB installation and recovered its initial concentration after 84 days passed and showed termination of reactivity of injected zero-valent iron/bio composite media for removal nitrate. Nitrite and ammonia ions found after installation of PRB indicates reductive removal of nitrate. And the outstanding increase of microorganism diversity and population of Betaproteobacteria Class which includes denitrification microorganism explains biologically reductive removal of nitrate in injected PRB.

• Journal of Soil and Groundwater Environment
• /
• v.10 no.6
• /
• pp.20-31
• /
• 2005

Three dimensional Monte-Carlo simulations of non-ergodic transport of a lion-reactive solute plume by steady-state groundwater flow under a uniform mean velocity in isotropic heterogeneous aquifers were conducted. The log-normally distributed hydraulic conductivity, K(x), is modeled as a random field. Significant efforts are made to reduce tile simulation uncertainties. Ensemble averages of the second spatial moments of the plume and plume centroid variances were simulated with 1600 Monte Carlo runs for three variances of log K, ${\sigma}_Y^2=0.09,\;0.23$, and 0.46, and three dimensionless lengths of line plume sources normal to the mean velocity. The simulated second spatial moment and the plume centroid variance in longitudinal direction fit well to the first order theoretical results while the simulated transverse moments are generally larger than the first order results. The first order theoretical results significantly underestimated the simulated dimensionless transverse moments for the aquifers of large ${\sigma}_Y^2$ and large dimensionless time. The ergodic condition for the second spatial moments is far from reaching in all cases simulated, and transport In transverse directions may reach ergodic condition much slower than that in longitudinal direction. The evolution of the contaminant transported in a heterogeneous aquifer is not affected by the shape of the initial plume but affected mainly by the degree of the heterogeneity and the size of the initial plume.

• Tunnel and Underground Space
• /
• v.28 no.5
• /
• pp.400-425
• /
• 2018

This study presents the research results and current status of the DECOVALEX-2019 project Task B. Task B named 'Fault slip modelling' is aiming at developing a numerical method to simulate the coupled hydro-mechanical behavior of fault, including slip or reactivation, induced by water injection. The first research step of Task B is a benchmark simulation which is designed for the modelling teams to familiarize themselves with the problem and to set up their own codes to reproduce the hydro-mechanical coupling between the fault hydraulic transmissivity and the mechanically-induced displacement. We reproduced the coupled hydro-mechanical process of fault slip using TOUGH-FLAC simulator. The fluid flow along a fault was modelled with solid elements and governed by Darcy's law with the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by creating interface elements between two separating rock blocks in FLAC3D. A methodology to formulate the hydro-mechanical coupling relations of two different hydraulic aperture models and link the solid element of TOUGH2 and the interface element of FLAC3D was suggested. In addition, we developed a coupling module to update the changes in geometric features (mesh) and hydrological properties of fault caused by water injection at every calculation step for TOUGH-FLAC simulator. Then, the transient responses of the fault, including elastic deformation, reactivation, progressive evolutions of pathway, pressure distribution and water injection rate, to stepwise pressurization were examined during the simulations. The results of the simulations suggest that the developed model can provide a reasonable prediction of the hydro-mechanical behavior related to fault reactivation. The numerical model will be enhanced by continuing collaboration and interaction with other research teams of DECOLVAEX-2019 Task B and validated using the field data from fault activation experiments in a further study.