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

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

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.19 no.2
• /
• pp.87-98
• /
• 2016

During the period between July 3 and 27 of 2009, water samples were collected from the Russian coast at a depth of 30m from 26 stations (including Ulleung and Japan basins) onboard the Russian survey vessel R/V Lavrentyev following 4 lines (D, R, E, and A). The samples were analyzed for nutrients and chlorophyll a contents. All parameters exhibited higher values in warm waters than in cold waters ($NH_4:1.8-fold$, $PO_4:1.8-fold$, $SiO_2:1.2-fold$, and chlorophyll-${\alpha}$:1.9-fold), except nitrates, which was 1.4-fold higher in cold waters than in warm waters. The horizontal distribution of ammonia, phosphate, and chlorophyll-${\alpha}$ was very similar to each other and showed the highest values in the waters near Russia, where a upwelling influence of cold current and bottom water prevails, while relatively low distribution was observed at the Ulleung Basin. On the other hand, nitrates showed the highest concentration at the Ulleung Basin, which is under the direct influence of the Tsushima warm water, and showed a gradual decrease northward. The N/P ratio showed the highest value in the Tsushima middle water, rather than in the North Korean Cold Water, the Tsushima Warm Water was the primary source of nitrate flow into the East Sea. However, the average concentration of phosphate in the warm waters was < $0.2{\mu}M$, thereby limiting phytoplankton growth, while a high concentration of phosphate in cold waters showed a direct correlation with chlorophyll-${\alpha}$. The results of principal component analysis for the identification of primary factors that influence the marine environment showed that principal component I was water temperature and principal component II was influenced chlorophyll-${\alpha}$ and nutrients. Therefore, Study area has greatest influenced by water temperature, and clearly distinct cold and warm water regions were observed in the East Sea.

• Journal of Bio-Environment Control
• /
• v.5 no.2
• /
• pp.215-235
• /
• 1996

The thermal analysis by mathematical model simulation makes it possible to reasonably predict heating and/or cooling requirements of certain greenhouses located under various geographical and climatic environment. It is another advantages of model simulation technique to be able to make it possible to select appropriate heating system, to set up energy utilization strategy, to schedule seasonal crop pattern, as well as to determine new greenhouse ranges. In this study, the control pattern for greenhouse microclimate is categorized as cooling and heating. Dynamic model was adopted to simulate heating requirements and/or energy conservation effectiveness such as energy saving by night-time thermal curtain, estimation of Heating Degree-Hours(HDH), long time prediction of greenhouse thermal behavior, etc. On the other hand, the cooling effects of ventilation, shading, and pad ||||&|||| fan system were partly analyzed by static model. By the experimental work with small size model greenhouse of 1.2m$\times$2.4m, it was found that cooling the greenhouse by spraying cold water directly on greenhouse cover surface or by recirculating cold water through heat exchangers would be effective in greenhouse summer cooling. The mathematical model developed for greenhouse model simulation is highly applicable because it can reflects various climatic factors like temperature, humidity, beam and diffuse solar radiation, wind velocity, etc. This model was closely verified by various weather data obtained through long period greenhouse experiment. Most of the materials relating with greenhouse heating or cooling components were obtained from model greenhouse simulated mathematically by using typical year(1987) data of Jinju Gyeongnam. But some of the materials relating with greenhouse cooling was obtained by performing model experiments which include analyzing cooling effect of water sprayed directly on greenhouse roof surface. The results are summarized as follows : 1. The heating requirements of model greenhouse were highly related with the minimum temperature set for given greenhouse. The setting temperature at night-time is much more influential on heating energy requirement than that at day-time. Therefore It is highly recommended that night- time setting temperature should be carefully determined and controlled. 2. The HDH data obtained by conventional method were estimated on the basis of considerably long term average weather temperature together with the standard base temperature(usually 18.3$^{\circ}C$). This kind of data can merely be used as a relative comparison criteria about heating load, but is not applicable in the calculation of greenhouse heating requirements because of the limited consideration of climatic factors and inappropriate base temperature. By comparing the HDM data with the results of simulation, it is found that the heating system design by HDH data will probably overshoot the actual heating requirement. 3. The energy saving effect of night-time thermal curtain as well as estimated heating requirement is found to be sensitively related with weather condition: Thermal curtain adopted for simulation showed high effectiveness in energy saving which amounts to more than 50% of annual heating requirement. 4. The ventilation performances doting warm seasons are mainly influenced by air exchange rate even though there are some variations depending on greenhouse structural difference, weather and cropping conditions. For air exchanges above 1 volume per minute, the reduction rate of temperature rise on both types of considered greenhouse becomes modest with the additional increase of ventilation capacity. Therefore the desirable ventilation capacity is assumed to be 1 air change per minute, which is the recommended ventilation rate in common greenhouse. 5. In glass covered greenhouse with full production, under clear weather of 50% RH, and continuous 1 air change per minute, the temperature drop in 50% shaded greenhouse and pad ＆ fan systemed greenhouse is 2.6$^{\circ}C$ and.6.1$^{\circ}C$ respectively. The temperature in control greenhouse under continuous air change at this time was 36.6$^{\circ}C$ which was 5.3$^{\circ}C$ above ambient temperature. As a result the greenhouse temperature can be maintained 3$^{\circ}C$ below ambient temperature. But when RH is 80%, it was impossible to drop greenhouse temperature below ambient temperature because possible temperature reduction by pad ||||&|||| fan system at this time is not more than 2.4$^{\circ}C$. 6. During 3 months of hot summer season if the greenhouse is assumed to be cooled only when greenhouse temperature rise above 27$^{\circ}C$, the relationship between RH of ambient air and greenhouse temperature drop($\Delta$T) was formulated as follows : $\Delta$T= -0.077RH＋7.7 7. Time dependent cooling effects performed by operation of each or combination of ventilation, 50% shading, pad ＆ fan of 80% efficiency, were continuously predicted for one typical summer day long. When the greenhouse was cooled only by 1 air change per minute, greenhouse air temperature was 5$^{\circ}C$ above outdoor temperature. Either method alone can not drop greenhouse air temperature below outdoor temperature even under the fully cropped situations. But when both systems were operated together, greenhouse air temperature can be controlled to about 2.0-2.3$^{\circ}C$ below ambient temperature. 8. When the cool water of 6.5-8.5$^{\circ}C$ was sprayed on greenhouse roof surface with the water flow rate of 1.3 liter/min per unit greenhouse floor area, greenhouse air temperature could be dropped down to 16.5-18.$0^{\circ}C$, whlch is about 1$0^{\circ}C$ below the ambient temperature of 26.5-28.$0^{\circ}C$ at that time. The most important thing in cooling greenhouse air effectively with water spray may be obtaining plenty of cool water source like ground water itself or cold water produced by heat-pump. Future work is focused on not only analyzing the feasibility of heat pump operation but also finding the relationships between greenhouse air temperature(T$_{g}$ ), spraying water temperature(T$_{w}$ ), water flow rate(Q), and ambient temperature(T$_{o}$).

• Environmental and Resource Economics Review
• /
• v.18 no.2
• /
• pp.279-309
• /
• 2009

Greenhouse gas emissions should be precisely forecast to reduce the emissions from industrial production processes. This study calculated the direct and indirect $CO_2$ emission intensities of 401 industries using the Input-Output tables 2003 and statistical data on the amount of energy use. This study had some limitations in drawing study findings because overseas data were used given the lack of domestic data. Other limiting factors included the oil distribution problems in the oil refinery sector, re-review of carbon neutral, and insufficient consideration of waste treatment. Nonetheless, this study is very meaningful since the direct and indirect $CO_2$ emission intensities of 401 industries were calculated. Specifically, this study considered from the zero-waste perspective the effects of waste, which attract interest worldwide since coke gas and gas from the steel industry are obtained as byproducts for the first time in Korea. According to the results of the analysis of $CO_2$ emission intensity per industry, typical industries whose indirect $CO_2$ emission intensity is high include crude steel making, Remicon, steel wire rods & track rail, cast iron, and iron reinforcing rods & bar steel. These industries produce products using the raw materials produced in the industrial sector whose $CO_2$ emission intensity is high. The representative industries whose direct $CO_2$ emission intensity is high include cement, pig iron, lime & plaster products, andcoal-based compounds. These industries extract raw ore from nature and refine them into raw materials that are useful in other industries. The findings in this study can be effectively used for the following case: estimation of target $CO_2$ emission reduction level reflecting each industrial sector's characteristics, calculation of potential emission reduction of each policy to reduce $CO_2$ emissions, identification of a firm's $CO_2$ emission level, and setting of the target level of emission reduction. Moreover, the findings in this study can be utilized widely in fields such as System of integrated Environmental and Economic Accounting(SEEA) and Material Flow Analysis(MFA) as the current topic of research in Korea.

• Journal of Wetlands Research
• /
• v.24 no.3
• /
• pp.204-212
• /
• 2022

The LID technique began to be applied in Korea after 2009, and LID facilities are installed and operated for rainwater management in business districts such as the Ministry of Environment, the Ministry of Land, Infrastructure and Transport, and LH Corporation, public institutions, commercial land, housing, parks, and schools. However, looking at domestic cases, the application cases and operation periods are insufficient compared to those outside the country, so appropriate design standards and measures for operation and maintenance are insufficient. In particular, LID facilities constructed using LID techniques need to maintain the environment inside LID facilities because hydrological and environmental effects are expressed by material circulation and energy flow. The LID facility is designed with the treatment capacity planned for the water circulation target, and the proper maintenance, vegetation, and soil conditions are periodically identified, and the efficiency is maintained as much as possible. In other words, the soil created in LID is a very important design element because LID facilities are expected to have effects such as water pollution reduction, flood reduction, water resource acquisition, and temperature reduction while increasing water storage and penetration capacity through water circulation construction. In order to maintain and manage the functions of LID facilities accurately, the current state of the facilities and the cycle of replacement and maintenance should be accurately known through various quantitative data such as soil contamination, snow removal effects, and vegetation criteria. This study was conducted to investigate the current status of LID facilities installed in Korea from 2009 to 2020, and analyze soil changes through the continuity and current status of LID facilities applied over the past 10 years after collecting soil samples from the soil layer. Through analysis of Saturn, organic matter, hardness, water contents, pH, electrical conductivity, and salt, some vegetation-type LID facilities more than 5 to 7 years after construction showed results corresponding to the lower grade of landscape design. Facilities below the lower level can be recognized as a point of time when maintenance is necessary in a state that may cause problems in soil permeability and vegetation growth. Accordingly, it was found that LID facilities should be managed through soil replacement and replacement.

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

• Economic and Environmental Geology
• /
• v.46 no.1
• /
• pp.51-61
• /
• 2013

Researches on the oil recovery enhancement using the nanotechnology has recently been studied in the United States. The previous researches has focused mainly on the flow characteristics of nanoparticles in porous media, and the stability of the nano-emulsion itself. However, the analysis did not deal with the size effects between a nano-emulsion and the pore size which has an important role when nano-emulsion flows in the porous media. In this research, nano-based emulsion was fabricated which is able to be applied for the enhanced oil recovery techniques and its characteristics was analyzed. In addition, in order to identify the characteristics of nano-emulsions flowing through the porous media, the size effect was analysed by filtering test. According to the results, when the emulsion was fabricated, SCA(Silane Coupling Agent) or PVA(Poly Vinyl Alcohol) are added to improve the stability of emulsion. As the ratio of the decane to water increased, the viscosity of emulsion and the droplet size also increased. For the filtering test at the atmospheric conditions, the droplet did not go through the filter; only the separated water from the emulsion was able to be filtered. This phenomenon occurred because the droplet was not able to overcome the capillary pressure. At the filtering test by suction pressure, most of the emulsion was filtered over the filter size of $60{\mu}m$. However, the ratio of filtration was rapidly degraded at less than $45{\mu}m$ filters. This is caused due to deformation and destruction of the droplet by strong shear stress when passing through the pore. The results from the study on the basic characteristic of nano-emulsion and filtering test will be expected to play as the important role for the fabrication of the stable nano-emulsion or the research on the recovery of residual oil in porous media.