• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.18 no.2_spc
• /
• pp.305-315
• /
• 2020

For geological disposal of radioactive wastes, a method was proposed to evaluate the radionuclide transport in the biosphere by calculating the elapsed time of nuclide migration. The radionuclides were supposed to be introduced from a natural barrier and reached a large surface water body following a groundwater flow in a shallow subsurface. The biosphere was defined as a shallow subsurface environment that included aquifers on a host rock. Using the proposed method, a calculation algorithm was established, and a computer code that implemented the algorithm was developed. The developed code was verified by comparing the simulation results of the simple cases with the results of the analytical solution and a public program, which has been widely used to evaluate the radiation dose using the radionuclide transport near the surface. A case study was constructed using the previous research for radionuclide transport from the hypothetical geological disposal repository. In the case study, the code calculated the mass discharge rate of radionuclide to a stream in the biosphere. Because the previous research only demonstrated the transport of radionuclides from the hypothetical repository to the host rock, the developed code in the present study could help identify the total transport of radionuclide along the complete pathway.

• The Journal of Engineering Geology
• /
• v.13 no.1
• /
• pp.51-65
• /
• 2003

A series of numerical simulations were carried out using STOMP( Subsurface Transport over Multiple Phase) simulator. The flow and distribution of LNAPL were analyzed in homogeneous fine and coarse sand. Vertical movement of LNAPL is faster in the coarser sand. But the total volume of LNAPL retained in the unsaturated zone is larger in the finer sand. A fine layer in the coarse sand domain is also simulated. The results showed that the retained LNAPL volume and shape are highly influenced by the Position of the fine layer. Flow and distributions of LNAPL were simulated when there were heterogeneous lenses in the sand domain. Water table fluctuation was also considered. In these cases, it was found that the heterogeneous lens was a barrier to LNAPL flow, and water table fluctuation stimulated the downward movement of retained LNAPL. The LNAPL flow and distribution observed in these numerical experiments show that in the subsurface environment, the behaviors of LNAPL highly depend on heterogeneities of unsaturated zone and the dynamic hydrogeologic condition such as water table fluctuation. These results can explain some of the complexity of LNAPL flow and distribution Patterns in LNAPL contaminated field sites.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.36 no.10
• /
• pp.979-988
• /
• 2012

The 3D structure of GDL for fuel cells was measured using high-resolution X-ray tomography in order to study material transport in the GDL. A computational algorithm has been developed to remove noise in the 3D image and construct 3D elements representing carbon fibers of GDL, which were used for both structural and fluid analyses. Changes in the pore structure of GDL under various compression levels were calculated, and the corresponding volume meshes were generated to evaluate the anisotropic permeability of gas within GDL as a function of compression. Furthermore, the transfer of liquid water and reactant gases was simulated by using the volume of fluid (VOF) and pore-network model (PNM) techniques. In addition, the simulation results of liquid water transport in GDL were validated by analogous experiments to visualize the diffusion of fluid in porous media. Through this research, a procedure for simulating the material transport in deformed GDL has been developed; this will help in optimizing the clamping force of fuel cell stacks as well as in determining the design parameters of GDL, such as thickness and porosity.

• Journal of Soil and Groundwater Environment
• /
• v.6 no.2
• /
• pp.69-81
• /
• 2001

The constitutive relation among capillary pressure, saturation and relative permeability should be predetermined in order to simulate immiscible water-gas flow in porous media. The relation between saturation and relative permeability becomes more important when the capillary force can be disregarded and viscous friction force governs the flow. In this study, a 2-dimensional finite difference numerical model was developed, in which the variation of viscosity with pressure and that of relative permeability with water saturation can be treated. Seven cases of parallel plate tests were performed in order to obtain the characteristic equation of relative permeability which would be used in. the developed numerical model. It was not possible, however, to match the curves of relative permeability from the plate tests with the existing emperical models. Consequently a logistic equation was proposed as a new emperical model. As this model was composed of the parameter involving aperture size, any aperture size of fracture can be applied to the model. For the purpose of verification, the characteristic equation of relative permeability was applied to the developed numerical model and the computed results were compared with those of plate test. As a result of application of numerical model, in order to check the field applicability, to single fracture surrounding an underground storage cavern, the simultaneous flow of water and propane gas was able to be simulated properly by the model.

• Journal of the Korean Society of Groundwater Environment
• /
• v.6 no.1
• /
• pp.23-32
• /
• 1999

A new numerical method using solution-adaptive grids (SAG) is developed to solve the Richards' equation (RE) for unsaturated flow in porous media. Using a grid generation technique, the SAG method automatically redistributes a fixed number of grid points during the flow process, so that more grid points are clustered in regions of large solution gradients. The method uses the coordinate transformation technique to employ a new transformed RE, which is solved with the standard finite difference method. The movement of grid points is incorporated into the transformed RE, and therefore all computation is performed on fixed grid points of the transformed domain without using any interpolation techniques. Thus, numerical difficulties arising from the movement of the wetting front during the infiltration process have been substantially overcome by the new method. Numerical experiments for an one-dimensional infiltration problem are presented to compare the SAG method to the modified Picard method using a fixed grid. Results show that accuracy of a SAG solution using 41 nodes is comparable with the solution of the fixed grid method using 201 nodes, while it requires only 50% of the CPU time. The global mass balance and the convergence of SAG solutions are strongly affected by the time step size (Δt) and the weighting parameter (${\gamma}$) used for generating solution-adaptive grids. Thus, the method requires automated readjustment of Δt and ${\gamma}$ to yield mass-conservative and convergent solutions, although it may increase computational costs. The method can be effective especially for simulating unsaturated flow and other transport problems involving the propagation of a sharp-front.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.23 no.1
• /
• pp.79-92
• /
• 2011

In this study, wave pressure of short-period gravity waves and tsunami acting on the upright section of the horizontal-slit type caisson placed on the impermeable or permeable seabed, which is a well-known permeable breakwater with a good wave controlling ability, are investigated via numerical simulations. Further, the permeable seabed was modeled as the porous media with porosity of 0.4. Using the numerical results, the effects of the seabed conditions on the wave pressure on the front wall and inside wall of the chamber have been studied. In the numerical simulations, short-period gravity waves and tsunami(solitary wave or bore) with the same amplitude to the gravity wave are considered. A numerical wave tank is used, which is able to consider a gas-liquid two-phase flow in the same calculation zone. Numerical results show that the wave pressure of the tsunami was 3~5 times higher than the short-period gravity waves acting on the front wall and it was 2~4 times higher than the short-period gravity waves acting on the inner wall.

• Economic and Environmental Geology
• /
• v.53 no.2
• /
• pp.213-220
• /
• 2020

When CO₂ foam is injected into the saline aquifer, the relative permeability of CO₂ decreases and its viscosity increases, thereby reducing mobility in porous media and ultimately improving CO₂ storge with enhanced sweep efficiency. In general, surfactants were used to fabricate CO₂ foam. Recently, nanoparticles have been used to form stable foam than surfactant. This paper introduces CO₂ storage technology using nanoparticle stabilized CO₂ foam. If the surface of the hydrophilic nanoparticles is partially modified into a CO₂-philic portion, the particles have an affinity for CO₂ and water, thus forming a stable CO₂ foam even in deep saline aquifers under high temperature and high salinity conditions, thereby it can be stored in the pores of the rock. In terms of economics, injection method using nanopaticle-stabilized CO₂ foam is more expensive than the conventional CO₂ injection, but it is estimated that it will have price competitiveness because the injection efficiency is improved. From an environmental point of view, it is possible to inject chemical substances such as surfactants and nanomaterials into aquifers or reservoirs for specific purposes such as pollutant removal and oil production. However, some studies have shown that nanoparticles and surfactants are toxic to aquatic animals, so environmentally proven substances should be used. Therefore, further research and development will be needed to study the production and injection of nanoparticle-stabilized CO₂ foam that are environmentally safe and economically reasonable.

• Journal of Conservation Science
• /
• v.33 no.6
• /
• pp.467-483
• /
• 2017

In this study, the degradation properties by temperature stress of $Araldite^{(R)}$ rapid-curing epoxy resin used for inorganic cultural heritages, was identified. The tensile and tensile shear strength of durability decreased for 12,624 hours at temperatures of $40{\sim}60^{\circ}C$. In terms of stability of external stress and temperature, the slow-curing epoxy was superior to the rapid-curing epoxy, and cultural heritage conservation plans should therefore consider the strength and stress properties of restoration materials. Color differences increased for 12,624 hours at temperatures of $40{\sim}60^{\circ}C$, and glossiness decreased. Both color and gloss stability were weak, which necessitates the improvement of optical properties. Thermal properties (weight loss, decomposition temperature, and glass transition temperature) of adhesives are linked to mechanical properties. Interfacial properties of the adherend and water vapor transmission rates of adhesives are linked to performance variation. For porous media (ceramics, brick, and stone), isothermal and isohumid environments are important. For outdoor artifacts on display in museums, changes in physical properties by exposure to varying environmental conditions need to be minimized. These results can be used as baseline data in the study of the degradation velocity and lifetime prediction of rapid-curing epoxy resin for the restoration of cultural heritages.

• The Journal of Engineering Geology
• /
• v.14 no.3 s.40
• /
• pp.301-311
• /
• 2004

This study is suggested a new dielectric tracer test method to understand geological structure of porous media and groundwater flow to use the dielectric constant which is one of electrical special quality of various geological materials. To measure their parameters, tracer material is made an ethanol mixing liquid(EML) having a same specific gravity of water. Also, soil materials are prepared a dielectric tracer test using the FDR system that could measure dielectric constant for saturated standard sand and river sand layers which have different initial porosity. To compare with their results, we discussed with the concentration variation of saline water having a saline concentration $3\%$ which is general tracer material by using the electro multi-meter system in the laboratory or field test. In two tracer experiment results, EML tracer test could confirm definitely EML concentration variation from each saturated soil layer as standard and river sands. However, tracer test of saline water $3\%$ concentration could not confirm permeating movement of water by degree of salinity change because these are settled at lower part column in a whole column area continuously. These causes are that specific gravity of saline water is heavier than water. That is, it could know that deposition of saline water is composed of lower part of soil column continuously independently of the direction of water into saturated soil material.

• Journal of Korean Society of Environmental Engineers
• /
• v.31 no.8
• /
• pp.655-662
• /
• 2009

In order to obtain knowledge on the design and operation of practical?scale Cylinder-Shaped Filters, pilot experiments were conducted to observe the effects of stepwise augmentation of production rate on water quality and clogging. A production rate increase from 0.8 L/min to 2.4 L/min did not appear to affect the removal efficiencies of BOD(Biochemical Oxygen Demand), SS(Suspended Solids), and turbidity, as the values were maintained around 80%, over 95%, and over 90% respectively;however, COD(Chemical Oxygen Demand) removal decreased from 44% to 29%. In addition, results indicated an inverse relationship between production rate and detention time required to remove unit contaminant concentration, the observed detention time in the filter to remove 1 mg-COD/L being 83 minutes for the production rate of 0.8 L/min and 45 minutes for the production rate of 2.4 L/min, suggesting that a relatively higher production rate is likely to be more advantageous in the purification of reservoir water when compared to a lower production rate. Clogging was observed to originate from the surface and advance to the center of the filter, and although clogging seemed to increase as the production rate increased, this did not cause any difficulties in normal functioning of the filter for more than 2 years of operation, suggesting that this filter system can be used effectively in the purification of reservoir water.